Recently solar flares exhibiting a double J-shaped ribbons in the lower solar atmosphere have been paid increasing attention in the context of extending the two-dimensional standard flare model to three dimensions, as motivated by the spatial correlation between photospheric current channels and flare ribbons. Here we study the electric currents through the photospheric area swept by flare ribbons (termed synthesized ribbon area or SRA), with a sample of 71 two-ribbon flares, of which 36 are J-shaped. Electric currents flowing through one ribbon are highly correlated with those through the other, therefore belonging to the same current system. The non-neutrality factor of this current system is independent of the flare magnitude, implying that both direct and return currents participate in flares. J-shaped flares are distinct from non-J-shaped flares in the following aspects: 1) Electric current densities within J-shaped SRA are significantly smaller than those within non-J-shaped SRA, but J-shaped SRA and its associated magnetic flux are also significantly larger. 2) Electric currents through SRA are positively correlated with the flare magnitude, but J-shaped flares show stronger correlation than non-J-shaped flares. 3) The majority (75\%) of J-shaped flares are eruptive, while the majority (86\%) of non-J-shaped flares are confined; accordingly, hosting active regions of J-shaped flares are more likely to be sigmoidal than non-J-shaped flares. Thus, J-shaped flares constitute a distinct subset of two-ribbon flares, probably the representative of eruptive ones. Further, we found that combining SRA and its associated magnetic flux has the potential to differentiate eruptive from confined flares.

Large-scale solar eruptions are believed to have a magnetic flux rope as the core structure. However, it remains elusive as to how the flux rope builds up and what triggers its eruption. Recent observations found that a prominence erupted following multiple episodes of "flux feeding". During each episode, a chromospheric fibril rose and merged with the prominence lying above. In this letter, we carried out 2.5-dimensional magnetohydrodynamic (MHD) numerical simulations to investigate whether the flux-feeding mechanism can explain such an eruption. The simulations demonstrate that the discrete emergence of small flux ropes can initiate eruptions by feeding axial flux into the preexistent flux rope until its total axial flux reaches a critical value. The onset of the eruption is dominated by an ideal MHD process. Our simulation results corroborate that the flux feeding is a viable mechanism to cause the eruption of solar magnetic flux ropes.

Star-forming galaxies are huge reservoirs of cosmic rays (CRs) and these CRs convert a significant fraction of their energy into $\gamma$-rays by colliding with the interstellar medium (ISM). Several nearby star-forming galaxies have been detected in GeV-TeV $\gamma$-rays. It is also found that the $\gamma$-ray luminosities in 0.1-100 GeV correlate well with indicators of star formation rates of the galaxies, such as the total infrared (IR) luminosity. In this paper, we report a systematic search for possible $\gamma$-ray emission from galaxies in the IRAS Revised Bright Galaxies Sample, using 11.4 years of $\gamma$-ray data taken by the Fermi Large Area Telescope (LAT). Two new galaxies, M33 and Arp 299, are detected significantly. The two galaxies are consistent with the empirical correlation between the $\gamma$-ray luminosity and total infrared luminosity, suggesting that their $\gamma$-ray emissions should mainly originate from CRs interacting with ISM. Nevertheless, there is a tentative evidence that the flux of the $\gamma$-ray emission from Arp~299 is variable. If the variability is true, part of the emission from Arp 299 should originate from the obscured AGN in this interacting galaxy system. In addition, we find that the $\gamma$-ray excess from M33 is located at the northeast region of the galaxy, where a supergiant H II region, NGC604, resides. This indicates that some bright star-forming regions in spiral galaxies could play a dominant role in the galaxy in producing $\gamma$-ray emission.

Solar flares and coronal mass ejections (CMEs) are closely coupled through magnetic reconnection. CMEs are usually accelerated impulsively within the low solar corona, synchronized with the impulsive flare energy release. We investigate the dynamic evolution of a fast CME and its associated X2.8 flare occurring on 2013 May 13. The CME experiences two distinct phases of enhanced acceleration, an impulsive one with a peak value of ~5 km s$^{-2}$ followed by an extended phase with accelerations up to 0.7 km s$^{-2}$. The two-phase CME dynamics is associated with a two-episode flare energy release. While the first episode is consistent with the "standard" eruption of a magnetic flux rope, the second episode of flare energy release is initiated by the reconnection of a large-scale loop in the aftermath of the eruption and produces stronger nonthermal emission up to $\gamma$-rays. In addition, this long-duration flare reveals clear signs of ongoing magnetic reconnection during the decay phase, evidenced by extended HXR bursts with energies up to 100--300 keV and intermittent downflows of reconnected loops for >4 hours. The observations reveal that the two-step flare reconnection substantially contributes to the two-phase CME acceleration, and the impulsive CME acceleration precedes the most intense flare energy release. The implications of this non-standard flare/CME observation are discussed.

We present the results of deep Chandra and XMM-Newton observations of a complex merging galaxy cluster Abell 2256 (A2256) that hosts a spectacular radio relic. The temperature and metallicity maps show clear evidence of a merger between the western subcluster and the primary cluster. We detect five X-ray surface brightness edges. Three of them near the cluster center are cold fronts (CFs): CF1 is associated with the infalling subcluster; CF2 is located in the east of the primary cluster and CF3 is to the west of the primary cluster core. The other two edges at cluster outskirts are shock fronts (SFs): SF1 near the radio relic in the NW has Mach numbers derived from the temperature and the density jumps respectively of $M_T=1.62\pm0.12$ and $M_\rho=1.23\pm0.06$; SF2 in the SE has $M_T=1.54\pm0.05$ and $M_\rho=1.16\pm0.13$. In the region of the radio relic there is no evidence for the correlation between X-ray and radio substructures, from which we estimate an upper limit for the inverse Compton emission, and therefore set a lower limit on the magnetic field ($\sim$ 450 kpc from primary cluster center) of $B>1.0\ \mu$G for a single power-law electron spectrum or $B>0.4\ \mu$G for a broken power-law electron spectrum. We propose a merger scenario including a primary cluster, a subcluster, and a group. Our merger scenario accounts for the X-ray edges, diffuse radio features, and galaxy kinematics, as well as projection effects.

White-light images from Heliospheric Imager-1 (HI1) onboard the Solar Terrestrial Relations Observatory (STEREO) provide 2-dimensional (2D) global views of solar wind transients traveling in the inner heliosphere from two perspectives. How to retrieve the hidden three-dimensional (3D) features of the transients from these 2D images is intriguing but challenging. In our previous work (Li et al., 2018), a 'correlation-aided' method is developed to recognize the solar wind transients propagating along the Sun-Earth line based on simultaneous HI1 images from two STEREO spacecraft. Here the method is extended from the Sun-Earth line to the whole 3D space to reconstruct the solar wind transients in the common field of view of STEREO HI1 cameras. We demonstrate the capability of the method by showing the 3D shapes and propagation directions of a coronal mass ejection (CME) and three small-scale blobs during 3-4 April 2010. Comparing with some forward modeling methods, we found our method reliable in terms of the position, angular width and propagation direction. Based on our 3D reconstruction result, an angular distorted, nearly North-South oriented CME on 3 April 2010 is revealed, manifesting the complexity of a CME's 3D structure.

Measurements of cosmic-ray (CR) positron fraction by PAMELA and other experiments have found an excess above 10 GeV relative to the standard predictions for secondary production in the interstellar medium (ISM). Although the excess has been mostly suggested to arise from some primary sources of positrons (such as pulsars and or annihilating dark matter particles), the almost constant flux ratio of $e^{+}/ \bar{p}$ argues for an alternative possibility that the excess positrons and antiprotons up to the highest energies are secondary products generated in hadronic interactions. Recently, Yang \& Aharonian (2019) revisit this possibility by assuming the presence of an additional population of CR nuclei sources. Here we examine this secondary product scenario using the \texttt{DRAGON} code, where the radiative loss of positrons is taken into account consistently. We confirm that the CR proton spectrum and the antiproton data can be explained by assuming the presence of an additional population of CR sources. However, the corresponding positron spectrum deviates from the measured data significantly above 100 GeV due to the strong radiative cooling. This suggests that, although hadronic interactions can explain the antiproton data, the corresponding secondary positron flux is still not enough to account for the AMS data. Hence contribution from some primary positron sources, such as pulsars or dark matter, is non-negligible.

Deep observation of the High Energy Stereoscopic System (HESS) on the most extended pulsar wind nebula HESS J1825-137 reveals an enhanced energy-dependent morphology, providing useful information on the particle transport mechanism in the nebula. We find that the energy-dependent morphology is consistent with a diffusion-dominated transport of electrons/positrons. It provides an alternative possible interpretation for the unusually large spatial extent (i.e., $\gtrsim 100$pc) of the nebula, which could then be attributed to the diffusion of escaping electrons/positrons from a compact plerion. The influence of various model parameters on the energy-dependent extent of the nebula is studied in the diffusion-dominated scenario. We also show that the energy-dependent morphology of the nebula may also be used to study the spin-down history of the pulsar.

The concept of the Solar Ring mission was gradually formed from L5/L4 mission concept, and the proposal of its pre-phase study was funded by the National Natural Science Foundation of China in November 2018 and then by the Strategic Priority Program of Chinese Academy of Sciences in space sciences in May 2019. Solar Ring mission will be the first attempt to routinely monitor and study the Sun and inner heliosphere from a full 360-degree perspective in the ecliptic plane. The current preliminary design of the Solar Ring mission is to deploy six spacecraft, grouped in three pairs, on a sub-AU orbit around the Sun. The two spacecraft in each group are separated by about 30 degrees and every two groups by about 120 degrees. This configuration with necessary science payloads will allow us to establish three unprecedented capabilities: (1) determine the photospheric vector magnetic field with unambiguity, (2) provide 360-degree maps of the Sun and the inner heliosphere routinely, and (3) resolve the solar wind structures at multiple scales and multiple longitudes. With these capabilities, the Solar Ring mission aims to address the origin of solar cycle, the origin of solar eruptions, the origin of solar wind structures and the origin of severe space weather events. The successful accomplishment of the mission will advance our understanding of the star and the space environment that hold our life and enhance our capability of expanding the next new territory of human.

Very high energy photons from cosmological gamma-ray bursts (GRBs) are expected to interact with extragalactic background light (EBL) and produce electron-positron pairs when they propagate through intergalactic medium (IGM). These relativistic pairs will then up-scatter cosmic microwave background (CMB) photons and emit secondary GeV emission. Meanwhile, the motion of these pairs are deflected by intergalactic magnetic field (IGMF), so the secondary GeV photons arrive later than the primary emission. It has been suggested that the properties of the secondary GeV emission can be used to constrain IGMF. Recently, TeV gamma-ray emission has been detected, for the first time, from a GRB (GRB 190114C) by the MAGIC telescope and its steep ${\rm \gamma-ray}$ spectrum shows a clear evidence of absorption by EBL. We then constrain the IGMF with the GeV flux limit obtained from the $Fermi$-LAT observations. We find a limit of $>10^{-19.5}$ G for the coherence length of $\lambda \leq 1$ Mpc. Although this limit is weaker than that obtained by using blazars, it represents the first limit from ${\rm \gamma-ray}$ observations of GRBs, which provides an independent constraint on IGMF. We also find that, for transient ${\rm \gamma-ray}$ sources, one can choose a favorable time window to search for the echo emission at a particular energy.

The interaction between a supernova ejecta and the circum-stellar medium drives a strong shock wave which accelerates particles (i.e., electrons and protons). The radio and X-ray emission observed after the supernova explosion constitutes the evidence of the electron acceleration. The accelerated protons are expected to produce GeV-TeV gamma-ray emission via $pp$ collisions, but the flux is usually low since only a small fraction of the supernova kinetic energy is converted into the shock energy at the very early time. The low gamma-ray flux of the nearest supernova explosion, SN 1987A, agrees with this picture. Here we report a serendipitous discovery of a fading GeV gamma-ray source in spatial coincidence with the second nearest supernova--SN 2004dj from our gamma-ray survey of nearby star-forming galaxies with Fermi-LAT. The total gamma-ray energy released by SN 2004dj is about $6\times10^{47}{\rm erg}$. We interpret this gamma-ray emission arising from the supernova ejecta interacting with a surrounding high-density shell, which decelerates the ejecta and converts ~1% of the ejecta's kinetic energy to relativistic protons. In addition, our gamma-ray survey of nearby star-forming galaxies discovers GeV emissions from two star-forming galaxies, i.e., Arp 299 and M33, for the first time.

Primordial non-Gaussianity (PNG) is a key probe of the origins of primordial fluctuations in the early universe. It has been shown that multi-tracer measurements of large-scale structure can produce high-precision measurements of PNG. Future line intensity mapping surveys are well-suited to these measurements owing to their ability to rapidly survey large volumes and access the large scales at which PNG becomes important. In this paper, we explore for the first time how multi-tracer PNG measurements with intensity mapping surveys depend on the sub-galactic scale physics which drives line emission. We consider an example cross-correlation between CO maps, and find that the choice of astrophysical models has a substantial impact on $f_{\mathrm{NL}}$ measurements at sufficiently high signal-to-noise. This in effect creates a coupling between horizon-scale PNG measurements and the molecular cloud-scale interstellar medium. We discuss how these effects depend on noise level and survey design. We further find that the cross-correlation shot noise, an effect nearly unique to intensity mapping measurements, plays an important role in multi-tracer analyses and cannot be neglected.

The orientation, chirality, and dynamics of solar eruptive filaments is a key to understanding the magnetic field of coronal mass ejections (CMEs) and therefore to predicting the geoeffectiveness of CMEs arriving at Earth. However, confusion and contention remain over the relationship between the filament chirality, magnetic helicity, and sense of rotation during eruption. To resolve the ambiguity in observations, in this paper, we used stereoscopic observations to determine the rotation direction of filament apex and the method proposed by Chen et al. (2014) to determine the filament chirality. Our sample of 12 eruptive active-region filaments establishes a strong one-to-one relationship, i.e., during the eruption, sinistral/dextral filaments (located in the southern/northern hemisphere) rotate clockwise/counterclockwise when viewed from above, and corroborates a weak hemispheric preference, i.e., a filament and related sigmoid both exhibit a forward (reverse) S shape in the southern (northern) hemisphere, which suggests that the sigmoidal filament is associated with a low-lying magnetic flux rope with its axis dipped in the middle. As a result of rotation, the projected S shape of a filament is anticipated to be reversed during eruption.

A second peak in extreme ultraviolet sometimes appears during the gradual phase of solar flares, which is known as EUV late phase (ELP). Stereotypically ELP is associated with two separated sets of flaring loops with distinct sizes, and it has been debated whether ELP is caused by additional heating or extended plasma cooling in the longer loop system. Here we carry out a survey of 55 M-and-above GOES-class flares with ELP during 2010--2014. Based on the flare-ribbon morphology, these flares are categorized as circular-ribbon (19 events), two-ribbon (23 events), and complex-ribbon (13 events) flares. Among them, 22 events (40%) are associated with coronal mass ejections, while the rest are confined. An extreme ELP, with the late-phase peak exceeding the main-phase peak, is found in 48% of two-ribbon flares, 37% of circular-ribbon flares, and 31% of complex-ribbon flares, suggesting that additional heating is more likely present during ELP in two-ribbon than in circular-ribbon flares. Overall, cooling may be the dominant factor causing the delay of the ELP peak relative to the main-phase peak, because the loop system responsible for the ELP emission is generally larger than, and well separated from, that responsible for the main-phase emission. All but one of the circular-ribbon flares can be well explained by a composite "dome-plate" quasi-separatrix layer (QSL). Only half of these show a magnetic null point, with its fan and spine embedded in the dome and plate, respectively. The dome-plate QSL, therefore, is a general and robust structure characterizing circular-ribbon flares.

Prominence bubbles are cavities rising into quiescent prominences from below. The bubble-prominence interface is often the active location for the formation of plumes, which flow turbulently into quiescent prominences. Not only the origin of prominence bubbles is poorly understood, but most of their physical characteristics are still largely unknown. Here, we investigate the dynamical properties of a bubble, which is observed since its early emergence beneath the spine of a quiescent prominence on 20 October 2017 in the H$\alpha$ line-center and in $\pm$0.4 angstrom line-wing wavelengths by the 1-m New Vacuum Solar Telescope. We report the prominence bubble to be exhibiting a disparate morphology in the H$\alpha$ line-center compared to its line-wings' images, indicating a complex pattern of mass motion along the line-of-sight. Combining Doppler maps with flow maps in the plane of sky derived from a Nonlinear Affine Velocity Estimator, we obtained a comprehensive picture of mass motions revealing a counter-clockwise rotation inside the bubble; with blue-shifted material flowing upward and red-shifted material flowing downward. This sequence of mass motions is interpreted to be either outlining a kinked flux rope configuration of the prominence bubble or providing observational evidence of the internal kink instability in the prominence plasma.

Observations by HAWC and Milagro have detected spatially extended TeV sources surrounding middle-aged ($t\sim100-400 \,{\rm kyr}$) pulsars like Geminga and PSR B0656+14, which have been named "TeV Halos", representing very extended TeV pulsar wind nebulae (PWNe) powered by relatively old pulsars. A few more HAWC-detected sources have been suggested to be TeV halo candidates. In this paper, we search for possible GeV counterparts of three TeV halo candidates with Fermi Large Area Telescopes. We detect a new spatially extended GeV source in the vicinity of the TeV halo candidate 2HWC J1912+099, which is also detected by HESS (HESS J1912+101). We find that the size of the GeV source is significantly larger than that of the TeV emission measured by HESS, and a spatial template characteristic of a PWN can fit the GeV data. We suggest that the GeV source is an extended PWN powered by the central middle-aged pulsar PSR J1913+1011. This discovery favors the TeV halo scenario for the TeV source 2HWC J1912+099 (HESS J1912+101), although the possible shell-like morphology measured by HESS challenges this interpretation. Alternatively, the TeV emission could be dominated by a supernova remnant (SNR) via the hadronic process. Future multi-wavelength observations of the source and more precise measurements of the spatial profile of the TeV emission will be useful to distinguish between the two scenarios.

In this letter we propose that the X-ray and the TeV observations in the vicinity of Geminga can be understood in the framework of anisotropic diffusion of injected electrons/positrons. This interpretation only requires the turbulence in the vicinity of Geminga to be sub-Alfv\'enic with the local mean magnetic field direction approximately aligned with our line of sight towards Geminga, without invoking extreme conditions for the environment, such as an extremely small diffusion coefficient and a weak magnetic field of strength $<1\mu$G as suggested in previous literature.

We describe the design and assembly of the LUX-ZEPLIN experiment, a direct detection search for cosmic WIMP dark matter particles. The centerpiece of the experiment is a large liquid xenon time projection chamber sensitive to low energy nuclear recoils. Rejection of backgrounds is enhanced by a Xe skin veto detector and by a liquid scintillator Outer Detector loaded with gadolinium for efficient neutron capture and tagging. LZ is located in the Davis Cavern at the 4850' level of the Sanford Underground Research Facility in Lead, South Dakota, USA. We describe the major subsystems of the experiment and its key design features and requirements.

Recently, very high-energy photons above 100 GeV were reported to be detected from GRB~190114C and GRB~180720B at, respectively, 100-1000 s and 10 hours after the burst. We model the available broad-band data of both GRBs with the synchrotron plus synchrotron self-Compton (SSC) emission of the afterglow shocks. We find that the sub-TeV emission of GRB~180720B can be interpreted as the SSC emission from afterglow shocks expanding in a constant density circum-burst medium. The SSC emission of GRB~190114C dominates over the synchrotron component from GeV energies at $\sim100$ s, which can explain the possible hard spectrum of the GeV emission at this time. The extrapolated flux of this SSC component to sub-TeV energies can explain the high-significance detection of GRB~190114C by the MAGIC telescope. The parameter values (such as the circum-burst density and shock microphysical parameters) in the modeling are not unusual for both GRBs, implying that the detection of sub-TeV photons from these two bursts should be attributed to their large burst energies and low redshifts.

Identifying the counterparts of submillimetre (submm) galaxies (SMGs) in multiwavelength images is a critical step towards building accurate models of the evolution of strongly star-forming galaxies in the early Universe. However, obtaining a statistically significant sample of robust associations is very challenging due to the poor angular resolution of single-dish submm facilities. Recently, a large sample of single-dish-detected SMGs in the UKIDSS UDS field, a subset of the SCUBA-2 Cosmology Legacy Survey (S2CLS), was followed up with the Atacama Large Millimeter/submillimeter Array (ALMA), which has provided the resolution necessary for identification in optical and near-infrared images. We use this ALMA sample to develop a training set suitable for machine-learning (ML) algorithms to determine how to identify SMG counterparts in multiwavelength images, using a combination of magnitudes and other derived features. We test several ML algorithms and find that a deep neural network performs the best, accurately identifying 85 per cent of the ALMA-detected optical SMG counterparts in our cross-validation tests. When we carefully tune traditional colour-cut methods, we find that the improvement in using machine learning is modest (about 5 per cent), but importantly it comes at little additional computational cost. We apply our trained neural network to the GOODS-North field, which also has single-dish submm observations from the S2CLS and deep multiwavelength data but little high-resolution interferometric submm imaging, and we find that we are able to classify SMG counterparts for 36/67 of the single-dish submm sources. We discuss future improvements to our ML approach, including combining ML with spectral energy distribution-fitting techniques and using longer wavelength data as additional features.

A high-energy muon neutrino event, IceCube-170922A, was recently discovered in both spatial and temporal coincidence with a gamma-ray flare of the blazar TXS 0506+056. It has been shown, with standard one-zone models, that neutrinos can be produced in the blazar jet via hadronic interactions, but with a flux which is mostly limited by the X-ray data. In this work, we explore the neutrino production from TXS 0506+056 by invoking two physically distinct emission zones in the jet, separated by the broad line region (BLR). Using the Doppler-boosted radiation of the BLR as the target photon field, the inner zone accounts for the neutrino and gamma-ray emission via $p\gamma$ interactions and inverse Compton scattering respectively, while the outer zone produces the optical and X-ray emission via synchrotron and synchrotron self-Compton processes. The different conditions of the two zones allow us to suppress the X-ray emission from the electromagnetic cascade, and set a much higher upper limit on the muon neutrino flux (i.e., $\sim 10^{-11}\rm erg~cm^{-2}s^{-1}$) than in one-zone models. We compare, in detail, our scenario with one-zone models discussed in the literature, and argue that differentiating between such scenarios will become possible with next generation neutrino telescopes, such as IceCube-Gen2.

The Gamma-Ray Integrated Detectors (GRID) is a space mission concept dedicated to monitoring the transient gamma-ray sky in the energy range from 10 keV to 2 MeV using scintillation detectors onboard CubeSats in low Earth orbits. The primary targets of GRID are the gamma-ray bursts (GRBs) in the local universe. The scientific goal of GRID is, in synergy with ground-based gravitational wave (GW) detectors such as LIGO and VIRGO, to accumulate a sample of GRBs associated with the merger of two compact stars and study jets and related physics of those objects. It also involves observing and studying other gamma-ray transients such as long GRBs, soft gamma-ray repeaters, terrestrial gamma-ray flashes, and solar flares. With multiple CubeSats in various orbits, GRID is unaffected by the Earth occultation and serves as a full-time and all-sky monitor. Assuming a horizon of 200 Mpc for ground-based GW detectors, we expect to see a few associated GW-GRB events per year. With about 10 CubeSats in operation, GRID is capable of localizing a faint GRB like 170817A with a 90% error radius of about 10 degrees, through triangulation and flux modulation. GRID is proposed and developed by students, with considerable contribution from undergraduate students, and will remain operated as a student project in the future. The current GRID collaboration involves more than 20 institutes and keeps growing. On August 29th, the first GRID detector onboard a CubeSat was launched into a Sun-synchronous orbit and is currently under test.

It has long been suggested that nearby pulsars within $\sim 1 \,{\rm kpc}$ are the leading candidate of the 10-500 GeV cosmic-ray positron excess measured by PAMELA and other experiments. The recent measurement of surface brightness profile of TeV nebulae surrounding Geminga and PSR~B0656+14 by the High-Altitude Water Cherenkov Observatory (HAWC) suggests inefficient diffusion of particles from the sources, giving rise to a debate on the pulsar interpretation of the cosmic-ray positron excess. Here we argue that GeV observations provide more direct constraints on the positron density in the TeV nebulae in the energy range of 10-500 GeV and hence on the origin of the observed positron excess. Motivated by this, we search for GeV emission from the TeV nebulae with the \textsl{Fermi} Large Area Telescope (LAT). No spatially-extended GeV emission is detected from these two TeV nebulae in the framework of two-zone diffusion spatial templates, suggesting a relatively low density of GeV electrons/positrons in the TeV nebulae. A joint modelling of the data from HAWC and \textsl{Fermi}-LAT disfavors Geminga and PSR~B0656+14 as the dominant source of the positron excess at $\sim 50-500$ GeV for the usual Kolmogorov-type diffusion, while for an energy-independent diffusion, a dominant part of the positron excess contributed by them cannot be ruled out by the current data.

The Large High Altitude Air Shower Observatory (LHAASO) project is a new generation multi-component instrument, to be built at 4410 meters of altitude in the Sichuan province of China, with the aim to study with unprecedented sensitivity the spec trum, the composition and the anisotropy of cosmic rays in the energy range between 10$^{12}$ and 10$^{18}$ eV, as well as to act simultaneously as a wide aperture (one stereoradiant), continuously-operated gamma ray telescope in the energy range between 10$^{11}$ and $10^{15}$ eV. The experiment will be able of continuously surveying the TeV sky for steady and transient sources from 100 GeV to 1 PeV, t hus opening for the first time the 100-1000 TeV range to the direct observations of the high energy cosmic ray sources. In addition, the different observables (electronic, muonic and Cherenkov/fluorescence components) that will be measured in LHAASO will allow to investigate origin, acceleration and propagation of the radiation through a measurement of energy spec trum, elemental composition and anisotropy with unprecedented resolution. The remarkable sensitivity of LHAASO in cosmic rays physics and gamma astronomy would play a key-role in the comprehensive general program to explore the High Energy Universe. LHAASO will allow important studies of fundamental physics (such as indirect dark matter search, Lorentz invariance violation, quantum gravity) and solar and heliospheric physics. In this document we introduce the concept of LHAASO and the main science goals, providing an overview of the project.

Recently, the High Altitude Water Cherenkov (HAWC) collaboration reported the discovery of the TeV halo around the Geminga pulsar. The TeV emission is believed to originate from inverse Compton scattering of pulsar-injected electrons/positrons off cosmic microwave background photons. In the mean time, these electrons should inevitably radiate X-ray photons via the synchrotron radiation, providing a useful constraint on the magnetic field in the TeV halo. In this work, we analyse the data of XMM-Newton and Chandra, and obtain an upper limit for the diffuse X-ray flux in a region of $600''$ around the Geminga pulsar, which is at a level of $\lesssim 10^{-14}\rm erg\,cm^{-2}s^{-1}$. Through a numerical modelling on both the X-ray and the TeV observations assuming isotropic diffusion of injected electrons/positrons, we find the magnetic field inside the TeV halo is required to be $<1\mu$G, which is significantly weaker than the typical magnetic field in the interstellar medium. The weak magnetic field together with the small diffusion coefficient inferred from HAWC's observation implies that the Bohm limit of particle diffusion may probably have been achieved in the TeV halo. We also discuss alternative possibilities for the weak X-ray emission, such as the hadronic origin of the TeV emission or a specific magnetic field topology, in which a weak magnetic field and a very small diffusion coefficient might be avoided.

Similar to the Sun, other stars shed mass and magnetic flux via ubiquitous quasi-steady wind and episodic stellar coronal mass ejections (CMEs). We investigate the mass loss rate via solar wind and CMEs as a function of solar magnetic variability represented in terms of sunspot number and solar X-ray background luminosity. We estimate the contribution of CMEs to the total solar wind mass flux in the ecliptic and beyond, and its variation over different phases of the solar activity cycles. The study exploits the number of sunspots observed, coronagraphic observations of CMEs near the Sun by SOHO/LASCO, in situ observations of the solar wind at 1 AU by WIND, and GOES X-ray flux during solar cycle 23 and 24. We note that the X-ray background luminosity, occurrence rate of CMEs and ICMEs, solar wind mass flux, and associated mass loss rates from the Sun do not decrease as strongly as the sunspot number from the maximum of solar cycle 23 to the next maximum. Our study confirms a true physical increase in CME activity relative to the sunspot number in cycle 24. We show that the CME occurrence rate and associated mass loss rate can be better predicted by X-ray background luminosity than the sunspot number. The solar wind mass loss rate which is an order of magnitude more than the CME mass loss rate shows no obvious dependency on cyclic variation in sunspot number and solar X-ray background luminosity. These results have implications to the study of solar-type stars.

Multi-wavelength observations show that Abell 1367 (A1367) is a dynamically young cluster, with at least two subclusters merging along the SE-NW direction. With the wide-field XMM-Newton mosaic of A1367, we discover a previously unknown merger shock at the NW edge of the cluster. We estimate the shock Mach number from the density and temperature jumps as $M_{\rho}=1.21\pm0.08$ and $M_T=1.60\pm0.07$, respectively. This shock region also corresponds to a radio relic discovered with the VLA and GBT, which could be produced by the shock re-acceleration of pre-existing seed relativistic electrons. We suggest that some of the seed relativistic electrons originate from late-type, star-forming galaxies in this region.

Although many high-energy neutrinos detected by the IceCube telescope are believed to have anextraterrestrial origin, their astrophysical sources remain a mystery. Recently, an unprecedenteddiscovery of a high-energy muon neutrino event coincident with a multiwavelength flare from ablazar, TXS 0506+056, shed some light on the origin of the neutrinos. It is usually believed that ablazar is produced by a relativistic jet launched from an accreting supermassive black hole (SMBH).Here we show that the high-energy neutrino event can be interpreted by the inelastic hadronuclearinteractions between the accelerated cosmic-ray protons in the relativistic jet and the dense gasclouds in the vicinity of the SMBH. Such a scenario only requires a moderate proton power in thejet, which could be much smaller than that required in the conventional hadronic model whichinstead calls upon the photomeson process. Meanwhile, the flux of the multiwavelength flare fromthe optical to gamma-ray band can be well explained by invoking a second radiation zone in thejet at a larger distance to the SMBH. In our model, the neutrino emission lasts a shorter time thanthe multiwavelength flare so the neutrino event is not necessarily correlated with the flare but it is probably accompanied by a spectrum hardening above a few GeV.

The Sun's atmosphere is frequently disrupted by coronal mass ejections (CMEs), coupled with flares and energetic particles. In the standard picture, the coupling is explained by magnetic reconnection at a vertical current sheet connecting the flare loops and the CME, with the latter embedding a helical magnetic structure known as flux rope. As it jumps upward due to instabilities or loss of equilibrium, the flux rope stretches the overlying coronal loops so that oppositely directed field is brought together underneath, creating the current sheet. However, both the origin of flux ropes and their nascent paths toward eruption remain elusive. Here we present an observation of how a stellar-sized CME bubble evolves continuously from plasmoids, mini flux ropes that are barely resolved, within half an hour. The eruption initiates when plasmoids springing from a vertical current sheet merge into a leading plasmoid occupying the upper tip of the current sheet. Rising at increasing speed to stretch the overlying loops, this leading plasmoid then expands impulsively into the CME bubble, in tandem with hard X-ray bursts. This observation illuminates for the first time a complete CME evolutionary path that has the capacity to accommodate a wide variety of plasma phenomena by bridging the gap between micro-scale dynamics and macro-scale activities.

It is often envisaged that dense filament material lies in the dips of magnetic field lines belonging to either a sheared arcade or a magnetic flux rope. But it is also debated which configuration correctly depicts filaments' magnetic structure, due to our incapacity to measure the coronal magnetic field. In this paper, we address this issue by employing mass motions in an active-region filament to diagnose its magnetic structure. The disturbance in the filament was driven by a surge initiated at the filament's eastern end in the NOAA active region 12685, which was observed by the 1-m New Vacuum Solar Telescope (NVST) in the H$\alpha$ line center and line wing ($\pm0.4$~{\AA}). Filament material predominately exhibits two kinds of motions, namely, rotation about the spine and longitudinal oscillation along the spine. The former is evidenced by antisymmetric Doppler shifts about the spine; the latter features a dynamic barb with mass extending away from the H$\alpha$ spine until the transversal edge of the EUV filament channel. The longitudinal oscillation in the eastern section of the filament is distinct from that in the west, implying that the underlying field lines have different lengths and curvature radii. The composite motions of filament material suggest a double-decker host structure with mixed signs of helicity, comprising a flux rope atop a sheared-arcade system.

It is well accepted that a magnetic flux rope (MFR) is a critical component of many coronal mass ejections (CMEs), yet how it evolves toward eruption remains unclear. Here we investigate the continuous evolution of a pre-existing MFR, which is rooted in strong photospheric magnetic fields and electric currents. The evolution of the MFR is observed by the Solar Terrestrial Relations Observatory (STEREO) and the Solar Dynamics Observatory (SDO) from multiple viewpoints. From STEREO's perspective, the MFR starts to rise slowly above the limb five hours before it erupts as a halo CME on 2012 June 14. In SDO observations, conjugate dimmings develop on the disk, simultaneously with the gradual expansion of the MFR, suggesting that the dimmings map the MFR's feet. The evolution comprises a two-stage gradual expansion followed by another stage of rapid acceleration/eruption. Quantitative measurements indicate that magnetic twist of the MFR increases from 1.0 +/- 0.5 to 2.0 +/- 0.5 turns during the five-hour expansion, and further increases to about 4.0 turns per AU when detected as a magnetic cloud at 1 AU two day later. In addition, each stage is preceded by flare(s), implying reconnection is actively involved in the evolution and eruption of the MFR. The implications of these measurements on the CME initiation mechanisms are discussed.

We study the requirement on the jet power in the conventional p-$\gamma$ models (photopion production and Bethe-Heitler pair production) for TeV BL Lac objects. We select a sample of TeV BL Lac objects whose SEDs are difficult to be explained by the one-zone leptonic model. Based on the relation between the p-$\gamma$ interaction efficiency and the opacity of $\gamma\gamma$ absorption, we find that the detection of TeV emission poses upper limits on the p-$\gamma$ interaction efficiencies in these sources and hence minimum jet powers can be derived accordingly. We find that the obtained minimum jet powers exceed the Eddington luminosity of the supermassive black holes. Implications for the accretion mode of the supermassive black hole in these BL Lac objects and the origin of their TeV emissions are discussed.

Various studies have implied the existence of a gaseous halo around the Galaxy extending out to 100 kpc. Galactic cosmic rays (CRs) that propagate to the halo, either by diffusion or by convection with the possibly existing large-scale Galactic wind, can interact with the gas therein and produce gamma-rays via proton-proton collision. We calculate the cosmic ray distribution in the halo and the gamma-ray flux, and explore the dependence of the result on model parameters such as diffusion coefficient, CR luminosity, CR spectral index. We find that the current measurement of isotropic gamma-ray background at $\lesssim$TeV with Fermi Large Area Telescope already approaches a level that can provide interesting constraints on the properties of Galactic cosmic ray (e.g., with CR luminosity $L_{CR}\leq 10^{41}$erg/s). We also discuss the possibilities of the Fermi bubble and IceCube neutrinos originating from the proton-proton collision between cosmic rays and gas in the halo, as well as the implication of our results for the baryon budget of the hot circumgalactic medium of our Galaxy. Given that the isotropic gamma-ray background is likely to be dominated by unresolved extragalactic sources, future telescopes may extract more individual sources from the IGRB, and hence put even more stringent restriction on the relevant quantities (such as Galactic cosmic ray luminosity and baryon budget in the halo) in the presence of a turbulent halo that we consider.

Solar eruptions occurring at different places within a relatively short time interval are considered to be sympathetic. However, it is difficult to determine whether there exists a cause and effect between them. Here we study a failed and a successful filament eruption following an X1.8-class flare on 2014 December 20, in which slipping-like magnetic reconnections serve as a key causal link among the eruptions. Reconnection signatures and effects are: at both sides of the filament experiencing the failed eruption, serpentine ribbons extend along chromospheric network to move away from the filament, while a hot loop apparently grows above it; at the filament undergoing the successful eruption, overlying cold loops contract, while coronal dimming appears at both sides even before the filament eruption. These effects are understood by reconnections continually transforming magnetic fluxes overlying one filament to the other, which adjusts how the magnetic field decays with increasing height above the filaments in opposite trends, therefore either strengthening or weakening the magnetic confinement of each filament.

Large-scale propagating fronts are frequently observed during solar eruptions, yet it is open whether they are waves or not, partly because the propagation is modulated by coronal structures, whose magnetic field we still cannot measure. However, when a front impacts coronal structures, an opportunity arises for us to look into the magnetic properties of both interacting parties in the low-$\beta$ corona. Here we studied large-scale EUV fronts accompanying three coronal mass ejections (CMEs), each originating from a kinking rope-like structure in the NOAA active region (AR) 12371. These eruptions were homologous and the surrounding coronal structures remained stationary. Hence we treated the events as one observed from three different viewing angles, and found that the primary front directly associated with the CME consistently transmits through 1) a polar coronal hole, 2) the ends of a crescent-shaped equatorial coronal hole, leaving a stationary front outlining its AR-facing boundary, and 3) two quiescent filaments, producing slow and diffuse secondary fronts. The primary front also propagates along an arcade of coronal loops and slows down due to foreshortening at the far side, where local plasma heating is indicated by an enhancement in 211 {\AA} (Fe XIV) but a dimming in 193 {\AA} (Fe XII) and 171 {\AA} (Fe IX). The strength of coronal magnetic field is therefore estimated to be $\sim\,$2 G in the polar coronal hole and $\sim\,$4 G in the coronal arcade neighboring the active region. These observations substantiate the wave nature of the primary front and shed new light on slow fronts.

Three-dimensional magnetic topology is crucial to understanding the explosive release of magnetic energy in the corona during solar flares. Much attention has been given to the pre-flare magnetic topology to identify candidate sites of magnetic reconnection, yet it is unclear how the magnetic reconnection and its attendant topological changes shape the eruptive structure and how the topology evolves during the eruption. Here we employed a realistic, data-constrained magnetohydrodynamic simulation to study the evolving magnetic topology for an X9.3 eruptive flare that occurred on 2017 September 6. The simulation successfully reproduces the eruptive features and processes in unprecedented detail. The numerical results reveal that the pre-flare corona contains multiple twisted flux systems with different connections, and during the eruption, these twisted fluxes form a coherent flux rope through tether-cutting-like magnetic reconnection below the rope. Topological analysis shows that the rising flux rope is wrapped by a quasi-separatrix layer, which intersects itself below the rope, forming a topological structure known as hyperbolic flux tube, where a current sheet develops, triggering the reconnection. By mapping footpoints of the newly-reconnected field lines, we are able to reproduce both the spatial location and, for the first time, the temporal separation of the observed flare ribbons, as well as the dynamic boundary of the flux rope's feet. Futhermore, the temporal profile of the total reconnection flux is comparable to the soft X-ray light curve. Such a sophisticated characterization of the evolving magnetic topology provides important insight into the eventual understanding and forecast of solar eruptions.

The photospheric response to solar flares, also known as coronal back reaction, is often observed as sudden flare-induced changes in vector magnetic field and sunspot motions. However, it remains obscure whether evolving flare ribbons, the flare signature closest to the photosphere, are accompanied by changes in vector magnetic field therein. Here we explore the relationship between the dynamics of flare ribbons in the chromosphere and variations of magnetic fields in the underlying photosphere, using high-resolution off-band H-alpha images and near-infrared vector magnetograms of the M6.5 flare on 2015 June 22 observed with the 1.6 m Goode Solar Telescope. We find that changes of photospheric fields occur at the arrival of the flare ribbon front, thus propagating analogously to flare ribbons. In general, the horizontal field increases and the field lines become more inclined to the surface. When ribbons sweep through regions that undergo a rotational motion, the fields transiently turn more vertical with decreased horizontal field and inclination angle, and then restore and/or become more horizontal than before the ribbon arrival. The ribbon propagation decelerates near the sunspot rotation center, where the vertical field becomes permanently enhanced. Similar magnetic field changes are discernible in magnetograms from the Helioseismic and Magnetic Imager (HMI), and an inward collapse of coronal magnetic fields is inferred from the time sequence of non-linear force-free field models extrapolated from HMI magnetograms. We conclude that photospheric fields respond nearly instantaneously to magnetic reconnection in the corona.

The Heliospheric Imagers (HI) on board the two spacecraft of the Solar Terrestrial Relations Observatory (STEREO) provided white-light images of transients in the solar wind from dual perspectives from 2007 to 2014. In this paper, we develop a new method to identify and locate the transients automatically from simultaneous images from the two inner telescopes, known as HI-1, based on a correlation analysis. Correlation coefficient (cc) maps along the Sun-Earth line are constructed for the period from 1 Jan 2010 to 28 Feb 2011. From the maps, transients propagating along the Sun-Earth line are identified, and a 27-day periodic pattern is revealed, especially for small-scale transients. Such a periodicity in the transient pattern is consistent with the rotation of the Sun's global magnetic structure and the periodic crossing of the streamer structures and slow solar wind across the Sun-Earth line, and this substantiates the reliability of our method and the high degree of association between the small-scale transients of the slow solar wind and the coronal streamers. Besides, it is suggested by the cc map that small-scale transients along the Sun-Earth line are more frequent than large-scale transients by a factor of at least 2, and that they quickly diffused into background solar wind within about 40 Rs in terms of the signal-to-noise ratio of white-light emissions. The method provides a new tool to reconstruct inhomogeneous structures in the heliosphere from multiple perspectives.

Vela X is a nearby pulsar wind nebula (PWN) powered by a $\sim 10^4$ year old pulsar. Modeling of the spectral energy distribution of the Vela X PWN has shown that accelerated electrons have largely escaped from the confinement, which is likely due to the disruption of the initially confined PWN by the SNR reverse shock. The escaped electrons propagate to the earth and contribute to the measured local cosmic-ray (CR) electron spectrum. We find that the escaped CR electrons from Vela X would hugely exceed the measured flux by HESS at $\sim 10$ TeV if the standard diffusion coefficient for the interstellar medium is used. We propose that the diffusion may be highly inefficient around Vela X and find that a spatially-dependent diffusion can lead to CR flux consistent with the HESS measurement. Using a two-zone model for the diffusion around Vela X, we find that the diffusion coefficient in the inner region of a few tens of parsecs should be $<10^{28}{\rm cm^2 s^{-1}}$ for $\sim10$ TeV CR electrons, which is about two orders of magnitude lower than the standard value for ISM. Such inefficient diffusion around PWN resembles the case of the Geminga and Monogem PWNe, suggesting that inefficient diffusion may be common in the vicinity of PWNe spanning a wide range of ages.

Recently, a high-energy neutrino event IceCube-170922A in the spatial and temporal coincidence with the flaring gamma-ray blazar TXS 0506+056 was reported. A neutrino outburst between September 2014 and March 2015 was discovered in the same direction by a further investigation of $9.5$ years of IceCube data, while the blazar is in a quiescent state during the outburst with a gamma-ray flux only about one-fifth of the neutrino flux. In this letter, we propose the neutrino outburst originates from the interaction between a relativistic jet and a dense gas cloud which may be formed via the tidally disrupted envelope of a red giant being blown by the impact of the jet. Gamma-ray photons and electron/positron pairs that are produced correspondingly will induce electromagnetic cascades. Comptonization of the cascade emission inside the cloud forms an X-ray photon field with Wien distribution. GeV flux is suppressed due to the absorption by the Comptonized photon field and, as a result, a hard spectrum above 10 GeV is formed. The gamma-ray spectrum predicted in our model is consistent with the Fermi-LAT data of TXS 0506+056.

Many galaxy clusters have giant halos of non-thermal radio emission, indicating the presence of relativistic electrons in the clusters. Relativistic protons may also be accelerated by merger and/or accretion shocks in galaxy clusters. These cosmic-ray (CR) electrons and/or protons are expected to produce gamma-rays through inverse-Compton scatterings or inelastic $pp$ collisions respectively. Despite of intense efforts in searching for high-energy gamma-ray emission from galaxy clusters, conclusive evidence is still missing so far. Here we report the discovery of $\ge 200$ MeV gamma-ray emission from the Coma cluster direction with an unbinned likelihood analysis of the 9 years of {\it Fermi}-LAT Pass 8 data. The gamma-ray emission shows a spatial morphology roughly coincident with the giant radio halo, with an apparent excess at the southwest of the cluster. Using the test statistic analysis, we further find tentative evidence that the gamma-ray emission at the Coma center is spatially extended. The extended component has an integral energy flux of $\sim 2\times 10^{-12}{\rm \ erg\ cm^{-2}\ s^{-1}}$ in the energy range of 0.2 - 300 GeV and the spectrum is soft with a photon index of $\simeq-2.7$. Interpreting the gamma-ray emission as arising from CR proton interaction, we find that the volume-averaged value of the CR to thermal pressure ratio in the Coma cluster is about $\sim 2\%$. Our results show that galaxy clusters are likely a new type of GeV gamma-ray sources, and they are probably also giant reservoirs of CR protons.

Centaurus A (Cen A) is the nearest radio galaxy discovered as a very-high-energy (VHE; 100 GeV-100 TeV) $\gamma$-ray source by the High Energy Stereoscopic System (H.E.S.S.). It is a faint VHE $\gamma$-ray emitter, though its VHE flux exceeds both the extrapolation from early Fermi-LAT observations as well as expectations from a (misaligned) single-zone synchrotron-self Compton (SSC) description. The latter satisfactorily reproduces the emission from Cen A at lower energies up to a few GeV. New observations with H.E.S.S., comparable in exposure time to those previously reported, were performed and eight years of Fermi-LAT data were accumulated to clarify the spectral characteristics of the $\gamma$-ray emission from the core of Cen A. The results allow us for the first time to achieve the goal of constructing a representative, contemporaneous $\gamma$-ray core spectrum of Cen A over almost five orders of magnitude in energy. Advanced analysis methods, including the template fitting method, allow detection in the VHE range of the core with a statistical significance of 12$\sigma$ on the basis of 213 hours of total exposure time. The spectrum in the energy range of 250 GeV-6 TeV is compatible with a power-law function with a photon index $\Gamma=2.52\pm0.13_{\mathrm{stat}}\pm0.20_{\mathrm{sys}}$. An updated Fermi-LAT analysis provides evidence for spectral hardening by $\Delta\Gamma\simeq0.4\pm0.1$ at $\gamma$-ray energies above $2.8^{+1.0}_{-0.6}$ GeV at a level of $4.0\sigma$. The fact that the spectrum hardens at GeV energies and extends into the VHE regime disfavour a single-zone SSC interpretation for the overall spectral energy distribution (SED) of the core and is suggestive of a new $\gamma$-ray emitting component connecting the high-energy emission above the break energy to the one observed at VHE energies.

In a strong magnetic field, a high-energy photon can be absorbed and then produce an electron-positron pair. The produced electron/positron will in turn radiate a high-energy photon via synchrotron radiation, which then initiates a cascade. We built a one-dimensional Monte-Carlo code to study the development of the cascade especially after it reaches the saturated status, when almost all the energy of the primary particles transfers to the photons. The photon spectrum in this status has a cut-off due to the absorption by magnetic fields, which is much sharper than the exponential one. Below the cut-off, the spectral energy distribution (SED) manifest itself as a broken power-law with a spectral index of $0.5$ and $0.125$, respectively, below and above the broken energy. The SED can be fitted by a simple analytical function, which is solely determined by the product of the cascade scale $R$ and the magnetic field perpendicular to the motion of the particle B_{\perp}, with an accuracy better than 96\%. The similarity of the spectrum to that from the cascade in an isotropic black-body photon field is also studied.

Spectral lines are among the most powerful signatures for dark matter (DM) annihilation searches in very-high-energy $\gamma$-rays. The central region of the Milky Way halo is one of the most promising targets given its large amount of DM and proximity to Earth. We report on a search for a monoenergetic spectral line from self-annihilations of DM particles in the energy range from 300 GeV to 70 TeV using a two-dimensional maximum likelihood method taking advantage of both the spectral and spatial features of signal versus background. The analysis makes use of Galactic Center (GC) observations accumulated over ten years (2004 - 2014) with the H.E.S.S. array of ground-based Cherenkov telescopes. No significant $\gamma$-ray excess above the background is found. We derive upper limits on the annihilation cross section $\langle\sigma v\rangle$ for monoenergetic DM lines at the level of $\sim4\times10^{-28}$ cm$^{3}$s$^{-1}$ at 1 TeV, assuming an Einasto DM profile for the Milky Way halo. For a DM mass of 1 TeV, they improve over the previous ones by a factor of six. The present constraints are the strongest obtained so far for DM particles in the mass range 300 GeV - 70 TeV. Ground-based $\gamma$-ray observations have reached sufficient sensitivity to explore relevant velocity-averaged cross sections for DM annihilation into two $\gamma$-ray photons at the level expected from the thermal relic density for TeV DM particles.

A search for new supernova remnants (SNRs) has been conducted using TeV gamma-ray data from the H.E.S.S. Galactic plane survey. As an identification criterion, shell morphologies that are characteristic for known resolved TeV SNRs have been used. Three new SNR candidates were identified in the H.E.S.S. data set with this method. Extensive multiwavelength searches for counterparts were conducted. A radio SNR candidate has been identified to be a counterpart to HESS J1534-571. The TeV source is therefore classified as a SNR. For the other two sources, HESS J1614-518 and HESS J1912+101, no identifying counterparts have been found, thus they remain SNR candidates for the time being. TeV-emitting SNRs are key objects in the context of identifying the accelerators of Galactic cosmic rays. The TeV emission of the relativistic particles in the new sources is examined in view of possible leptonic and hadronic emission scenarios, taking the current multiwavelength knowledge into account.

Aims. We study gamma-ray emission from the shell-type supernova remnant (SNR) RX J0852.0$-$4622 to better characterize its spectral properties and its distribution over the SNR. Methods. The analysis of an extended High Energy Spectroscopic System (H.E.S.S.) data set at very high energies (E > 100 GeV) permits detailed studies, as well as spatially resolved spectroscopy, of the morphology and spectrum of the whole RX J0852.0$-$4622 region. The H.E.S.S. data are combined with archival data from other wavebands and interpreted in the framework of leptonic and hadronic models. The joint Fermi-LAT-H.E.S.S. spectrum allows the direct determination of the spectral characteristics of the parent particle population in leptonic and hadronic scenarios using only GeV-TeV data. Results. An updated analysis of the H.E.S.S. data shows that the spectrum of the entire SNR connects smoothly to the high-energy spectrum measured by Fermi-LAT. The increased data set makes it possible to demonstrate that the H.E.S.S. spectrum deviates significantly from a power law and is well described by both a curved power law and a power law with an exponential cutoff at an energy of $E_\mathrm{cut} = (6.7 \pm 1.2_\mathrm{stat} \pm 1.2_\mathrm{syst})$ TeV. The joint Fermi-LAT-H.E.S.S. spectrum allows the unambiguous identification of the spectral shape as a power law with an exponential cutoff. No significant evidence is found for a variation of the spectral parameters across the SNR, suggesting similar conditions of particle acceleration across the remnant. A simple modeling using one particle population to model the SNR emission demonstrates that both leptonic and hadronic emission scenarios remain plausible. It is also shown that at least a part of the shell emission is likely due to the presence of a pulsar wind nebula around PSR J0855$-$4644.

Various observations are revealing the widespread occurrence of fast and powerful winds in active galactic nuclei (AGNs) that are distinct from relativistic jets, likely launched from accretion disks and interacting strongly with the gas of their host galaxies. During the interaction, strong shocks are expected to form that can accelerate non-thermal particles to high energies. Such winds have been suggested to be responsible for a large fraction of the observed extragalactic gamma-ray background (EGB) in the GeV-TeV range and the diffuse neutrino background in the PeV range, via the decay of neutral and charged pions generated in inelastic $pp$ collisions between protons accelerated by the forward shock and the ambient gas. However, previous studies did not properly account for processes such as adiabatic losses that may reduce the gamma-ray and neutrino fluxes significantly. We evaluate the production of gamma-rays and neutrinos by AGN-driven winds in some detail by modeling their hydrodynamic and thermal evolution, including the effects of their two-temperature structure. We find that they can only account for less than $\sim 30$% of the EGB flux, as otherwise the model would violate the independent upper limit derived from the diffuse isotropic gamma-ray background. If the neutrino spectral index is steep with $\Gamma\gtrsim 2.2$, a severe tension with the isotropic gamma-ray background would arise as long as the winds contribute more than $20$% of the IceCube neutrino flux in the $10-100$TeV range. Nevertheless, at energies $\gtrsim100$~TeV, we find that the IceCube neutrino flux may still be accountable by AGN-driven winds if the spectral index is as small as $\Gamma\sim2.0-2.1$. The detectability of gamma-ray point sources also provides important constraints on such scenarios.

Solar flare precursors depict constrained rate of energy release contrasting the imminent rapid energy release which calls for different regime of plasma processes to be at play. Due to subtle emission during the precursor phase, its diagnostics remain delusive, revealing either the non-thermal electrons (NTEs) or the thermal conduction to be the driver. In this regard, we investigate the chromospheric response during various phases of a B6.4 flare on August 20, 2005. Spatio-temporal investigation of flare ribbon enhancement during the precursor phase, carried out using spectra-images recorded in several wavelength positions on the H-alpha line profile, revealed its delayed response (180 seconds) compared to the X-ray emission, as well as sequential increment in the width of the line-profile which are indicative of a slow heating process. However, energy contained in the H-alpha emission during the precursor phase reach as high as 80% of that estimated during the main phase. Additionally, the plasma hydrodynamics during the precursor phase, as resulted from the application of a single-loop one-dimensional model, revealed the presence of power-law extension in the model generated X-ray spectra, with flux lower than the RHESSI background. Therefore, our multi-wavelength diagnostics and hydrodynamical modeling of the precursor emission indicates the role of a two-stage process. Firstly, reconnection triggered NTEs, although too small in flux to overcome the observational constraints, thermalize in the upper chromosphere. This leads to the generation of a slow conduction front which causes plasma heating during the precursor phase.

We study five sequential solar flares (SOL2015-08-07) occurring in Active Region 12396 observed with the Goode Solar Telescope (GST) at the BBSO, complemented by IRIS and SDO observations. The main flaring region is an arch filament system (AFS) consisting of multiple bundles of dark filament threads enclosed by semi-circular flare ribbons. We study the magnetic configuration and evolution of the active region by constructing coronal magnetic field models based on SDO/HMI magnetograms using two independent methods, i.e., the nonlinear force-free field (NLFFF) extrapolation and the flux rope insertion method. The models consist of multiple flux ropes with mixed signs of helicity, i.e., positive (negative) in the northern (southern) region, which is consistent with the GST observations of multiple filament bundles. The footprints of quasi-separatrix layers (QSLs) derived from the extrapolated NLFFF compare favorably with the observed flare ribbons. An interesting double-ribbon fine structure located at the east border of the AFS is consistent with the fine structure of the QSL's footprint. Moreover, magnetic field lines traced along the semi-circular footprint of a dome-like QSL surrounding the AFS are connected to the regions of significant helicity and Poynting flux injection. The maps of magnetic twist show that positive twist became dominant as time progressed, which is consistent with the injection of positive helicity before the flares. We hence conclude that these circular shaped flares are caused by 3D magnetic reconnection at the QSLs associated with the AFS possessing mixed signs of helicity.

The disintegration of solar filaments via mass drainage is a frequently observed phenomenon during a variety of filament activities. It is generally considered that the draining of dense filament material is directed by both gravity and magnetic field, yet the detailed process remains elusive. Here we report on a partial filament eruption during which filament material drains downward to the surface not only along the filament's legs, but to a remote flare ribbon through a fan-out curtain-like structure. It is found that the magnetic configuration is characterized by two conjoining dome-like quasi-sepratrix layers (QSLs). The filament is located underneath one QSL dome, whose footprint apparently bounds the major flare ribbons resulting from the filament eruption, whereas the remote flare ribbon matches well with the other QSL dome's far-side footprint. We suggest that the interaction of the filament with the overlying QSLs results in the splitting and disintegration of the filament.

The high-energy (> 100MeV) emission observed by Fermi-LAT during the prompt phase of some luminous gamma-ray bursts (GRBs) could arise from the cascade induced by interactions between accelerated protons and the radiation field of GRBs. The photomeson process, which is usually suggested to operate in such a hadronic explanation, requires a rather high proton energy (> 10^17eV) for an efficient interaction. However, whether GRBs can accelerate protons to such a high energy is far from guaranteed, although they have been suggested as the candidate source for ultrahigh-energy cosmic rays. In this work, we revisit the hadronic model for the prompt high-energy emission of GRBs with a smaller maximum proton energy than the usually adopted value estimated from the Bohm condition. In this case, the Bethe-Heitler pair production process becomes comparably important or even dominates over the photomeson process. We show that with a relatively low maximum proton energy with a Lorentz factor of 10^5 in the comoving frame, the cascade emission can still reproduce various types of high-energy spectrum of GRBs. For most GRBs without high-energy emission detected, the maximum proton energy could be even lower and relax the constraints on the parameters of GRB jet resulting from the fact of non-detection of GRB neutrinos by IceCube.

The solar corona is frequently disrupted by coronal mass ejections (CMEs), whose core structure is believed to be a flux rope made of helical magnetic field. This has become a "standard" picture although it remains elusive how the flux rope forms and evolves toward eruption. While 1/3 of the ejecta passing through spacecrafts demonstrate a flux-rope structure, the rest have complex magnetic fields. Are they originating from a coherent flux rope, too? Here we investigate the source region of a complex ejecta, focusing on a flare precursor with definitive signatures of magnetic reconnection, i.e., nonthermal electrons, flaring plasma, and bi-directional outflowing blobs. Aided by nonlinear force-free field modeling, we conclude that the reconnection occurs within a system of multiple braided flux ropes with different degree of coherency. The observation signifies the importance of internal structure and dynamics in understanding CMEs and in predicting their impacts on Earth.

Magnetic flux rope (MFR) is the core structure of the greatest eruptions, i.e., the coronal mass ejections (CMEs), on the Sun, and magnetic clouds are post-eruption MFRs in interplanetary space. There is a strong debate about whether or not a MFR exists prior to a CME and how the MFR forms/grows through magnetic reconnection during the eruption. Here we report a rare event, in which a magnetic cloud was observed sequentially by four spacecraft near Mercury, Venus, Earth and Mars, respectively. With the aids of a uniform-twist flux rope model and a newly developed method that can recover a shock-compressed structure, we find that the axial magnetic flux and helicity of the magnetic cloud decreased when it propagated outward but the twist increased. Our analysis suggests that the `pancaking' effect and `erosion' effect may jointly cause such variations. The significance of the `pancaking' effect is difficult to be estimated, but the signature of the erosion can be found as the imbalance of the azimuthal flux of the cloud. The latter implies that the magnetic cloud was eroded significantly leaving its inner core exposed to the solar wind at far distance. The increase of the twist together with the presence of the erosion effect suggests that the post-eruption MFR may have a high-twist core enveloped by a less-twisted outer shell. These results pose a great challenge to the current understanding on the solar eruptions as well as the formation and instability of MFRs.

PKS 0625-354 (z=0.055) was observed with the four H.E.S.S. telescopes in 2012 during 5.5 hours. The source was detected above an energy threshold of 200 GeV at a significance level of 6.1$\sigma$. No significant variability is found in these observations. The source is well described with a power-law spectrum with photon index $\Gamma =2.84 \pm 0.50_{stat} \pm 0.10_{syst}$ and normalization (at $E_0$=1.0 TeV) $N_0(E_0)=(0.58 \pm 0.22_{stat} \pm 0.12_{syst})\times10^{-12}$ TeV$^{-1}$cm$^{-2}$s$^{-1}$. Multi-wavelength data collected with Fermi-LAT, Swift-XRT, Swift-UVOT, ATOM and WISE are also analysed. Significant variability is observed only in the Fermi-LAT $\gamma$-ray and Swift-XRT X-ray energy bands. Having a good multi-wavelength coverage from radio to very high energy, we performed a broadband modelling from two types of emission scenarios. The results from a one zone lepto-hadronic, and a multi-zone leptonic models are compared and discussed. On the grounds of energetics, our analysis favours a leptonic multi-zone model. Models associated to the X-ray variability constraint supports previous results suggesting a BL Lac nature of PKS 0625-354, with, however, a large-scale jet structure typical of a radio galaxy.

The X-ray emission from the jets in Active Galactic Nuclei (AGN) carries important information on the distributions of relativistic electrons and magnetic fields on large scales. We reanalyze archival Chandra observations on the jet of M87 from 2000 to 2016 with a total exposure of 1460 kiloseconds to explore the X-ray emission characteristics along the jet. We investigate the variability behaviours of the nucleus and the inner jet component HST-1, and confirm indications for day-scale X-ray variability in the nucleus contemporaneous to the 2010 high TeV gamma-ray state. HST-1 shows a general decline in X-ray flux over the last few years consistent with its synchrotron interpretation. We extract the X-ray spectra for the nucleus and all knots in the jet, showing that they are compatible with a single power-law within the X-ray band. There are indications of the resultant X-ray photon index to exhibit a trend, with slight but significant index variations ranging from $\simeq 2.2$ (e.g. in knot D) to $\simeq 2.4-2.6$ (in the outer knots F, A, and B). When viewed in a multi-wavelength context, a more complex situation is arising. Fitting the radio to X-ray spectral energy distributions (SEDs) assuming a synchrotron origin, we show that a broken power-law electron spectrum with break energy $E_b$ around $1~(300\mu G/B)^{1/2}$ TeV allows a satisfactorily description of the multi-band SEDs for most of the knots. However, in the case of knots B, C and D we find indications that an additional high energy component is needed to adequately reproduce the broadband SEDs. We discuss the implications and suggest that a stratified jet model may account for the differences.

Shell-type supernova remnants (SNRs) are considered prime candidates for the acceleration of Galactic cosmic rays (CRs) up to the knee of the CR spectrum at $\mathrm{E} \approx \mathrm{3}\times \mathrm{10}^\mathrm{15}$ eV. Our Milky Way galaxy hosts more than 350 SNRs discovered at radio wavelengths and at high energies, of which 220 fall into the H.E.S.S. Galactic Plane Survey (HGPS) region. Of those, only 50 SNRs are coincident with a H.E.S.S source and in 8 cases the very high-energy (VHE) emission is firmly identified as an SNR. The H.E.S.S. GPS provides us with a legacy for SNR population study in VHE $\gamma$-rays and we use this rich data set to extract VHE flux upper limits from all undetected SNRs. Overall, the derived flux upper limits are not in contradiction with the canonical CR paradigm. Assuming this paradigm holds true, we can constrain typical ambient density values around shell-type SNRs to $n\leq 7~\textrm{cm}^\textrm{-3}$ and electron-to-proton energy fractions above 10~TeV to $\epsilon_\textrm{ep} \leq 5\times 10^{-3}$. Furthermore, comparisons of VHE with radio luminosities in non-interacting SNRs reveal a behaviour that is in agreement with the theory of magnetic field amplification at shell-type SNRs.

Context. Recently, the high-energy (HE, 0.1-100 GeV) $\gamma$-ray emission from the object LMC P3 in the Large Magellanic Cloud (LMC) has been discovered to be modulated with a 10.3-day period, making it the first extra-galactic $\gamma$-ray binary. Aims. This work aims at the detection of very-high-energy (VHE, >100 GeV) $\gamma$-ray emission and the search for modulation of the VHE signal with the orbital period of the binary system. Methods. LMC P3 has been observed with the High Energy Stereoscopic System (H.E.S.S.); the acceptance-corrected exposure time is 100 h. The data set has been folded with the known orbital period of the system in order to test for variability of the emission. Energy spectra are obtained for the orbit-averaged data set, and for the orbital phase bin around the VHE maximum. Results. VHE $\gamma$-ray emission is detected with a statistical significance of 6.4 $\sigma$. The data clearly show variability which is phase-locked to the orbital period of the system. Periodicity cannot be deduced from the H.E.S.S. data set alone. The orbit-averaged luminosity in the $1-10$ TeV energy range is $(1.4 \pm 0.2) \times 10^{35}$ erg/s. A luminosity of $(5 \pm 1) \times 10^{35}$ erg/s is reached during 20% of the orbit. HE and VHE $\gamma$-ray emissions are anti-correlated. LMC P3 is the most luminous $\gamma$-ray binary known so far.

In this paper, we present the detailed analysis of recurrent homologous jets originating from an emerging negative magnetic flux at the edge of an Active Region. The observed jets show multi-thermal features. Their evolution shows high consistence with the characteristic parameters of the emerging flux, suggesting that with more free magnetic energy, the eruptions tend to be more violent, frequent and blowout-like. The average temperature, average electron number density and axial speed are found to be similar for different jets, indicating that they should have been formed by plasmas from similar origins. Statistical analysis of the jets and their footpoint region conditions reveals a strong positive relationship between the footpoint-region total 131 {\AA} intensity enhancement and jets' length/width. Stronger linearly positive relationships also exist between the total intensity enhancement/thermal energy of the footpoint regions and jets' mass/kinetic/thermal energy, with higher cross-correlation coefficients. All the above results, together, confirm the direct relationship between the magnetic reconnection and the jets, and validate the important role of magnetic reconnection in transporting large amount of free magnetic energy into jets. It is also suggested that there should be more free energy released during the magnetic reconnection of blowout than of standard jet events.

We analyze the high-resolution observations from the 1.6m telescope at Big Bear Solar Observatory that cover an active region filament. Counter-streaming motions are clearly observed in the filament. The northern end of the counter-streaming motions extends to a light bridge, forming a spectacular circulation pattern around a sunspot, with clockwise motion in the blue wing and counterclockwise motion in the red wing as observed in Halpha off-bands. The apparent speed of the flow is around 10 to 60 km/s in the filament, decreasing to 5 to 20 km/s in the light bridge. The most intriguing results are the magnetic structure and the counter-streaming motions in the light bridge. Similar to those in the filament, magnetic fields show a dominant transverse component in the light bridge. However, the filament is located between opposite magnetic polarities, while the light bridge is between strong fields of the same polarity. We analyze the power of oscillations with the image sequences of constructed Dopplergrams, and find that the filament's counter-streaming motion is due to physical mass motion along fibrils, while the light bridge's counter-streaming motion is due to oscillation in the direction along the line-of-sight. The oscillation power peaks around 4 minutes. However, the section of the light bridge next to the filament also contains a component of the extension of the filament in combination with the oscillation, indicating that some strands of the filament are extended to and rooted in that part of the light bridge.

We search for high-energy gamma-ray emission from the binary neutron star merger GW170817 with the H.E.S.S. Imaging Air Cherenkov Telescopes. The observations presented here have been obtained starting only 5.3h after GW170817. The H.E.S.S. target selection identified regions of high probability to find a counterpart of the gravitational wave event. The first of these regions contained the counterpart SSS17a that has been identified in the optical range several hours after our observations. We can therefore present the first data obtained by a ground-based pointing instrument on this object. A subsequent monitoring campaign with the H.E.S.S. telescopes extended over several days, covering timescales from 0.22 to 5.2 days and energy ranges between $270\,\mathrm{GeV}$ to $8.55\,\mathrm{TeV}$. No significant gamma-ray emission has been found. The derived upper limits on the very-high-energy gamma-ray flux for the first time constrain non-thermal, high-energy emission following the merger of a confirmed binary neutron star system.

Coronal magnetic flux ropes are closely related to large-scale solar activities. Using a 2.5-dimensional time-dependent ideal magnetohydrodynamic (MHD) model in Cartesian coordinates, we carry out numerical simulations to investigate the evolution of a magnetic system consisting of a flux rope embedded in a fully-closed quadrupolar magnetic field with different photospheric flux distributions. It is found that, when the photospheric flux is not concentrated too much towards the polarity inversion line (PIL) and the constraint exerted by the background field is not too weak, the equilibrium states of the system are divided into two branches: the rope sticks to the photosphere for the lower branch and levitates in the corona for the upper branch. These two branches are connected by an upward catastrophe (from the lower branch to the upper) and a downward catastrophe (from the upper branch to the lower). Our simulations reveal that there exist both upward and downward catastrophes in quadrupolar fields, which may be either force-free or non-force-free. The existence and the properties of these two catastrophes are influenced by the photospheric flux distribution, and a downward catastrophe is always paired with an upward catastrophe. Comparing the decay indices in catastrophic and non-catastrophic cases, we infer that torus unstable might be a necessary but not sufficient condition for a catastrophic system.

The rotational motion of solar jets is believed to be a signature of the untwisting process resulting from magnetic reconnection, which takes place between twisted closed magnetic loops (i.e., magnetic flux ropes) and open magnetic field lines. The identification of the pre-existing flux rope, and the relationship between the twist contained in the rope and the number of turns the jet experiences, are then vital in understanding the jet-triggering mechanism. In this paper, we will perform a detailed analysis of imaging, spectral and magnetic field observations of four homologous jets, among which the fourth one releases a twist angle of 2.6$\pi$. Non-linear force free field extrapolation of the photospheric vector magnetic field before the jet eruption presents a magnetic configuration with a null point between twisted and open fields - a configuration highly in favor of the eruption of solar jets. The fact that the jet rotates in the opposite sense of handness to the twist contained in the pre-eruption photospheric magnetic field, confirms the unwinding of the twist by the jet's rotational motion. Temporal relationship between jets' occurrence and the total negative flux at their source region, together with the enhanced magnetic submergence term of the photospheric Poynting flux, shows that these jets are highly associated with local magnetic flux cancellation.

The H.E.S.S. collaboration has discovered a new very high energy (VHE, E $>$ 0.1 TeV) $\gamma$-ray source, HESS J1741-302, located in the Galactic plane. Despite several attempts to constrain its nature, no plausible counterpart has been found so far at X-ray and MeV/GeV $\gamma$-ray energies, and the source remains unidentified. An analysis of 145-hour of observations of HESS J1741-302 at VHEs has revealed a steady and relatively weak TeV source ($\sim$1$\%$ of the Crab Nebula flux), with a spectral index of $\Gamma$ = 2.3 $\pm$ 0.2$_{\text{stat}}$ $\pm$ 0.2$_{\text{sys}}$, extending to energies up to 10 TeV without any clear signature of a cut-off. In a hadronic scenario, such a spectrum implies an object with particle acceleration up to energies of several hundred TeV. Contrary to most H.E.S.S. unidentified sources, the angular size of HESS J1741-302 is compatible with the H.E.S.S. point spread function at VHEs, with an extension constrained to be below 0.068$^{\circ}$ at a 99$\%$ confidence level. The $\gamma$-ray emission detected by H.E.S.S. can be explained both within a hadronic scenario, due to collisions of protons with energies of hundreds of TeV with dense molecular clouds, and in a leptonic scenario, as a relic pulsar wind nebula, possibly powered by the middle-aged (20 kyr) pulsar PSR B1737-30. A binary scenario, related to the compact radio source 1LC 358.266+0.038 found to be spatially coincident with the best fit position of HESS J1741-302, is also envisaged.

Index of H.E.S.S. conference proceedings to the 35th ICRC, Busan, Korea

The origin of the prompt high-energy ($>100$MeV) emission of Gamma-ray Bursts (GRBs), detected by the Large Area Telescope (LAT) on board the Fermi Gamma-ray Space Telescope, is still under debate, for which both the external shock origin and internal dissipation origin have been suggested. In the internal dissipation scenario, the high energy emission is expected to exhibit significant temporal variability, tracking the keV/MeV fast variable behavior. Here, we report a detailed analysis on the Fermi data of GRB~170214A, which is sufficiently bright in the high energy to enable a quantitative analysis of the correlation between high-energy emission and keV/MeV emission with high statistics. Our result shows a clear temporal correlation between high-energy and keV/MeV emission in the whole prompt emission phase as well as in two decomposed short time intervals. Such correlation behavior is also found in some other bright LAT GRBs, i.e., GRB 080916C, 090902B and 090926A. For these GRBs as well as GRB 090510, we also find the rapid temporal variability in the high-energy emission. We thus conclude that the prompt high-energy emission in these bright LAT GRBs should be due to internal origin.

Very high-energy gamma-rays (VHE, E>100 GeV) propagating over cosmological distances can interact with the low-energy photons of the extragalactic background light (EBL) and produce electron-positron pairs. The transparency of the universe to VHE gamma-rays is then directly related to the spectral energy distribution (SED) of the EBL. The observation of features in the VHE energy spectra of extragalactic sources allows the EBL to be measured, which otherwise is very difficult to determine. An EBL-model independent measurement of the EBL SED with the H.E.S.S. array of Cherenkov telescopes is presented. It is obtained by extracting the EBL absorption signal from the reanalysis of high-quality spectra of blazars. From H.E.S.S. data alone the EBL signature is detected at a significance of 9.5 sigma, and the intensity of the EBL obtained in different spectral bands is presented together with the associated gamma-ray horizon.

We report observations of an eruptive X2.8 flare on 2013 May 13, which shows two distinct episodes of energy release in the impulsive phase. The first episode is characterized by the eruption of a magnetic flux rope, similar to the energy-release process in most standard eruptive flares. While the second episode, which is stronger than the first normal one and shows enhanced high-energy X-ray and even $\gamma$-ray emissions, is closely associated with magnetic reconnection of a large-scale loop in the aftermath of the eruption. The reconnection inflow of the loop leg is observed in the Solar Dynamics Observatory (SDO)/Atmospheric Imaging Assembly (AIA) 304 \r{A} passband and accelerates towards the reconnection region to a speed as high as $\sim$130 km/s. Simultaneously the corresponding outflow jets are observed in the AIA hot passbands with a speed of $\sim$740 km/s and mean temperature of $\sim$14 MK. RHESSI observations show a strong burst of hard X-ray (HXR) and $\gamma$-ray emissions with hard electron spectra of $\delta\approx3$, exhibiting a soft-hard-harder behavior. A distinct altitude decrease of the HXR loop-top source coincides with the inward swing of the loop leg observed in the AIA 304 \r{A} passband, which is suggested to be related to the coronal implosion. This fast inflow of magnetic flux contained in the loop leg greatly enhances the reconnection rate and results in very efficient particle acceleration in the second-step reconnection, which also helps to achieve a second higher temperature peak up to T $\approx$ 30 MK.

The large jet kinetic power and non-thermal processes occurring in the microquasar SS 433 make this source a good candidate for a very high-energy (VHE) gamma-ray emitter. Gamma-ray fluxes have been predicted for both the central binary and the interaction regions between jets and surrounding nebula. Also, non-thermal emission at lower energies has been previously reported. We explore the capability of SS 433 to emit VHE gamma rays during periods in which the expected flux attenuation due to periodic eclipses and precession of the circumstellar disk periodically covering the central binary system is expected to be at its minimum. The eastern and western SS433/W50 interaction regions are also examined. We aim to constrain some theoretical models previously developed for this system. We made use of dedicated observations from MAGIC and H.E.S.S. from 2006 to 2011 which were combined for the first time and accounted for a total effective observation time of 16.5 h. Gamma-ray attenuation does not affect the jet/medium interaction regions. The analysis of a larger data set amounting to 40-80 h, depending on the region, was employed. No evidence of VHE gamma-ray emission was found. Upper limits were computed for the combined data set. We place constraints on the particle acceleration fraction at the inner jet regions and on the physics of the jet/medium interactions. Our findings suggest that the fraction of the jet kinetic power transferred to relativistic protons must be relatively small to explain the lack of TeV and neutrino emission from the central system. At the SS433/W50 interface, the presence of magnetic fields greater 10$\mu$G is derived assuming a synchrotron origin for the observed X-ray emission. This also implies the presence of high-energy electrons with energies up to 50 TeV, preventing an efficient production of gamma-ray fluxes in these interaction regions.

In this paper, we identified the magnetic source locations of 142 quasi- homologous (QH) coronal mass ejections (CMEs), of which 121 are from solar cycle (SC) 23 and 21 from SC 24. Among those CMEs, 63% originated from the same source location as their predecessor (defined as S-type), while 37% origi- nated from a different location within the same active region as their predecessor (defined as D-type). Their distinctly different waiting time distribution, peaking around 7.5 and 1.5 hours for S- and D-type CMEs, suggests that they might involve different physical mechanisms with different characteristic time scales. Through detailed analysis based on non-linear force free (NLFF) coronal mag- netic field modeling of two exemplary cases, we propose that the S-type QH CMES might involve a recurring energy release process from the same source location (by magnetic free energy replenishment), whereas the D-type QH CMEs can happen when a flux tube system disturbed by a nearby CME.

Context. Supra-arcade downflows (SADs) have been frequently observed during the gradual phase of solar flares near the limb. In coronal emission lines sensitive to flaring plasmas, they appear as tadpole-like dark voids against the diffuse fan-shaped "haze" above, flowing toward the well-defined flare arcade. Aims. We aim to investigate the evolution of SADs' thermal properties, and to shed light on the formation mechanism and physical processes of SADs. Methods. We carefully studied several selected SADs from two flare events and calculated their differential emission measures (DEMs) as well as DEM-weighted temperatures using data obtained by the Atmospheric Imaging Assembly (AIA) onboard the Solar Dynamic Observatory. Results. Our analysis shows that SADs are associated with a substantial decrease in DEM above 4 MK, which is 1-3 orders of magnitude smaller than the surrounding haze as well as the region before or after the passage of SADs, but comparable to the quiet corona. There is no evidence for the presence of the SAD-associated hot plasma (> 20 MK) in the AIA data, and this decrease in DEM does not cause any significant change in the DEM distribution as well as the DEM-weighted temperature, which supports this idea that SADs are density depletion. This depression in DEM rapidly recovers in the wake of the SADs studied, generally within a few minutes, suggesting that they are discrete features. In addition, we found that SADs in one event are spatio-temporally associated with the successive formation of post-flare loops along the flare arcade.

We studied the background field for 60 two-ribbon flares of M-and-above classes during 2011--2015. These flares are categorized into two groups, i.e., \emph{eruptive} and \emph{confined} flares, based on whether a flare is associated with a coronal mass ejection or not. The background field of source active regions is approximated by a potential field extrapolated from the $B_z$ component of vector magnetograms provided by the Helioseismic and Magnetic Imager. We calculated the decay index $n$ of the background field above the flaring polarity inversion line, and defined a critical height $h_\mathrm{crit}$ corresponding to the theoretical threshold ($n_\mathrm{crit}=1.5$) of the torus instability. We found that $h_\mathrm{crit}$ is approximately half of the distance between the centroids of opposite polarities in active regions, and that the distribution of $h_\mathrm{crit}$ is bimodal: it is significantly higher for confined flares than for eruptive ones. The decay index increases monotonously with increasing height for 86\% (84\%) of the eruptive (confined) flares but displays a saddle-like profile for the rest 14\% (16\%), which are found exclusively in active regions of multipolar field configuration. Moreover, $n$ at the saddle bottom is significantly smaller in confined flares than that in eruptive ones. These results highlight the critical role of background field in regulating the eruptive behavior of two-ribbon flares.

The diffuse very high-energy (VHE, > 100 GeV) gamma-ray emission observed in the central 200 pc of the Milky Way by H.E.S.S. was found to follow the dense matter distribution in the Central Molecular Zone (CMZ) up to a longitudinal distance of about 130 pc to the Galactic Centre (GC), where the flux rapidly decreases. Recent sensitive H.E.S.S. analysis revealed that the cosmic-ray (CR) density profile drops with the distance to the centre, making data compatible with a steady cosmic PeVatron at the GC. In this paper, we extend this analysis to obtain for the first time a detailed characterisation of the correlation with matter and to search for additional features and individual gamma-ray sources in the inner 200 pc. Taking advantage of 250 hours of H.E.S.S. data and improved analysis techniques we perform a detailed morphology study of the diffuse VHE emission observed from the GC ridge and reconstruct its total spectrum. We show that the emission correlated with dense matter covers the full CMZ and that its flux is about half the total diffuse emission flux. We also detect some emission at higher latitude likely produced by hadronic collisions of CRs in less dense regions of the GC interstellar medium. We detect an additional emission component centred on the GC and extending over about 15 pc that is consistent with the existence of a strong CR density gradient and confirms the presence of a CR accelerator at the very centre of our Galaxy. We show that the spectrum of the full ridge diffuse emission is compatible with the one previously derived from the central regions, suggesting that a single population of particles fills the entire CMZ. Finally, we report the discovery of a VHE gamma-ray source near the GC radio arc and argue that it is produced by the pulsar wind nebula candidate G0.13-0.11.

We report the intriguing large-scale dynamic phenomena associated with the M6.5 flare~(SOL2015-06-22T18:23) in NOAA active region 12371, observed by RHESSI, Fermi, and the Atmospheric Image Assembly (AIA) and Magnetic Imager (HMI) on the Solar Dynamic Observatory (SDO). The most interesting feature of this event is a third ribbon (R3) arising in the decay phase, propagating along a dimming channel (seen in EUV passbands) towards a neighboring sunspot. The propagation of R3 occurs in the presence of hard X-ray footpoint emission, and is broadly visible at temperatures from 0.6~MK to over 10~MK through the Differential Emission Measure (DEM) analysis. The coronal loops then undergo an apparent slipping motion following the same path of R3, after a $\sim$ 80~min delay. To understand the underlying physics, we investigate the magnetic configuration and the thermal structure of the flaring region. Our results are in favor of a slipping-type reconnection followed by the thermodynamic evolution of coronal loops. In comparison with those previously reported slipping reconnection events, this one proceeds across a particularly long distance ($\sim$60 Mm) over a long period of time ($\sim$50 min), and shows two clearly distinguished phases: the propagation of the footpoint brightening driven by nonthermal particle injection and the apparent slippage of loops governed by plasma heating and subsequent cooling.

The origin of the extended X-ray emission in the large-scale jets of active galactic nuclei (AGNs) poses challenges to conventional models of acceleration and emission. Although the electron synchrotron radiation is considered the most feasible radiation mechanism, the formation of the continuous large-scale X-ray structure remains an open issue. As astrophysical jets are expected to exhibit some turbulence and shearing motion, we here investigate the potential of shearing flows to facilitate an extended acceleration of particles and evaluate its impact on the resultant particle distribution. Our treatment incorporates systematic shear and stochastic second-order Fermi effects. We show that for typical parameters applicable to large-scale AGN jets, stochastic second-order Fermi acceleration, which always accompanies shear particle acceleration, can play an important role in facilitating the whole process of particle energization. We study the time-dependent evolution of the resultant particle distribution in the presence of second-order Fermi acceleration, shear acceleration, and synchrotron losses using a simple Fokker--Planck approach and provide illustrations for the possible emergence of a complex (multicomponent) particle energy distribution with different spectral branches. We present examples for typical parameters applicable to large-scale AGN jets, indicating the relevance of the underlying processes for understanding the extended X-ray emission and the origin of ultrahigh-energy cosmic rays.

It has been suggested that high-energy gamma-ray emission ($>100{\rm MeV}$) of nearby star-forming galaxies may be produced predominantly by cosmic rays colliding with the interstellar medium through neutral pion decay. The pion-decay mechanism predicts a unique spectral signature in the gamma-ray spectrum, characterized by a fast rising spectrum and a spectral break below a few hundreds of MeV. We here report the evidence of a spectral break around 500 MeV in the disk emission of Large Magellanic Cloud (LMC), which is found in the analysis of the gamma-ray data extending down to 60 MeV observed by {\it Fermi}-Large Area Telescope. The break is well consistent with the pion-decay model for the gamma-ray emission, although leptonic models, such as the electron bremsstrahlung emission, cannot be ruled out completely.

Runaway stars form bow shocks by ploughing through the interstellar medium at supersonic speeds and are promising sources of non-thermal emission of photons. One of these objects has been found to emit non-thermal radiation in the radio band. This triggered the development of theoretical models predicting non-thermal photons from radio up to very-high-energy (VHE, E $\geq 0.1$TeV) gamma rays. Subsequently, one bow shock was also detected in X-ray observations. However, the data did not allow discrimination between a hot thermal and a non-thermal origin. Further observations of different candidates at X-ray energies showed no evidence for emission at the position of the bow shocks either. A systematic search in the \textit{Fermi}-LAT energy regime resulted in flux upper limits for 27 candidates listed in the E-BOSS catalogue. Here we perform the first systematic search for VHE gamma-ray emission from bow shocks of runaway stars. Using all available archival H.E.S.S. data we search for very-high-energy gamma-ray emission at the positions of bow shock candidates listed in the second E-BOSS catalogue release. Out of the 73 bow shock candidates in this catalogue, 32 have been observed with H.E.S.S. None of the observed 32 bow shock candidates in this population study show significant emission in the H.E.S.S. energy range. Therefore, flux upper limits are calculated in five energy bins and the fraction of the kinetic wind power that is converted into VHE gamma rays is constrained. Emission from stellar bow shocks is not detected in the energy range between 0.14 and 18 TeV. The resulting upper limits constrain the level of VHE gamma-ray emission from these objects down to 0.1-1 % of the kinetic wind energy.

Solar flares are generally believed to be powered by free magnetic energy stored in the corona, but the build up of coronal energy alone may be insufficient for the imminent flare occurrence. The flare onset mechanism is a critical but less understood problem, insights into which could be gained from small-scale energy releases known as precursors, which are observed as small pre-flare brightenings in various wavelengths, and also from certain small-scale magnetic configurations such as the opposite polarity fluxes, where magnetic orientation of small bipoles is opposite to that of the ambient main polarities. However, high-resolution observations of flare precursors together with the associated photospheric magnetic field dynamics are lacking. Here we study precursors of a flare using unprecedented spatiotemporal resolution of the 1.6 m New Solar Telescope, complemented by novel microwave data. Two episodes of precursor brightenings are initiated at a small-scale magnetic channel (a form of opposite polarity fluxes) with multiple polarity inversions and enhanced magnetic fluxes and currents, lying near the footpoints of sheared magnetic loops. The low-atmospheric origin of these precursor emissions is corroborated by microwave spectra. We propose that the emerging magnetic channel field interacts with the sheared arcades to cause precursor brightenings at the main flare core region. These high-resolution results provide evidence of low-atmospheric small-scale energy release and possible relationship to the onset of the main flare.

The nine-year H.E.S.S. Galactic Plane Survey (HGPS) yielded the most uniform observation scan of the inner Milky Way in the TeV gamma-ray band to date. The sky maps and source catalogue of the HGPS allow for a systematic study of the population of TeV pulsar wind nebulae found throughout the last decade. To investigate the nature and evolution of pulsar wind nebulae, for the first time we also present several upper limits for regions around pulsars without a detected TeV wind nebula. Our data exhibit a correlation of TeV surface brightness with pulsar spin-down power $\dot{E}$. This seems to be caused both by an increase of extension with decreasing $\dot{E}$, and hence with time, compatible with a power law $R_\mathrm{PWN}(\dot{E}) \sim \dot{E}^{-0.65 \pm 0.20}$, and by a mild decrease of TeV gamma-ray luminosity with decreasing $\dot{E}$, compatible with $L_{1-10\,\mathrm{TeV}} \sim \dot{E}^{0.59 \pm 0.21}$. We also find that the offsets of pulsars with respect to the wind nebula centres with ages around 10 kyr are frequently larger than can be plausibly explained by pulsar proper motion and could be due to an asymmetric environment. In the present data, it seems that a large pulsar offset is correlated with a high apparent TeV efficiency $L_{1-10\,\mathrm{TeV}}/\dot{E}$. In addition to 14 HGPS sources considered as firmly identified pulsar wind nebulae and 5 additional pulsar wind nebulae taken from literature, we find 10 HGPS sources that form likely TeV pulsar wind nebula candidates. Using a model that subsumes the present common understanding of the very high-energy radiative evolution of pulsar wind nebulae, we find that the trends and variations of the TeV observables and limits can be reproduced to a good level, drawing a consistent picture of present-day TeV data and theory.

Muon pairs can be produced in the annihilation of ultra-high energy (UHE, $E \gtrsim 10^{18} \,\mathrm{eV}$) photons with low energy cosmic background radiation in the intergalactic space, giving birth to neutrinos. Although the branching ratio of muon pair production is low, products of other channels, which are mainly electron/positron pairs, will probably transfer most of their energies into the new generated UHE photon in the subsequent interaction with the cosmic background radiation via Compton scattering in deep Klein-Nishina regime. The regeneration of these new UHE photons then provides a second chance to produce the muon pairs, enhancing the neutrino flux. We investigate the neutrino production in the propagation of UHE photons in the intergalactic space at different redshifts, considering various competing processes such as pair production, double pair production for UHE photons, and triplet production and synchrotron radiation for UHE electrons. Regarding the least energetic outgoing particles as energy loss, we obtain the effective penetration length of the leading particle, as well as energy loss rate including the neutrino emission rate in the cascade process. Finally, we find that an extra component of UHE neutrinos will arise from the propagation of UHE cosmic rays due to the generated UHE photons and electron/positrons. However, the flux of this component is quite small, with a flux of at most $10\%$ of that of the conventional cosmogenic neutrino at a few EeV, in the absence of a strong intergalactic magnetic field and a strong cosmic radio background. The precise contribution of extra component depends on several factors, e.g., cosmic radio background, intergalactic magnetic field, the spectrum of proton, which will be discussed in this work.

Supernova remnants exhibit shock fronts (shells) that can accelerate charged particles up to very high energies. In the past decade, measurements of a handful of shell-type supernova remnants in very-high-energy gamma rays have provided unique insights into the acceleration process. Among those objects, RX$~$J1713.7-3946 (also known as G347.3-0.5) has the largest surface brightness, allowing us in the past to perform the most comprehensive study of morphology and spatially resolved spectra of any such very-high-energy gamma-ray source. Here we present extensive new H.E.S.S. measurements of RX$~$J1713.7-3946, almost doubling the observation time compared to our previous publication. Combined with new improved analysis tools, the previous sensitivity is more than doubled. The H.E.S.S. angular resolution of $0.048^\circ$ ($0.036^\circ$ above 2 TeV) is unprecedented in gamma-ray astronomy and probes physical scales of 0.8 (0.6) parsec at the remnant's location. The new H.E.S.S. image of RX$~$J1713.7-3946 allows us to reveal clear morphological differences between X-rays and gamma rays. In particular, for the outer edge of the brightest shell region, we find the first ever indication for particles in the process of leaving the acceleration shock region. By studying the broadband energy spectrum, we furthermore extract properties of the parent particle populations, providing new input to the discussion of the leptonic or hadronic nature of the gamma-ray emission mechanism.

Since only the magnetic conditions at the photosphere can be routinely observed in current observations, it is of great significance to find out the influences of photospheric magnetic conditions on solar eruptive activities. Previous studies about catastrophe indicated that the magnetic system consisting of a flux rope in a partially open bipolar field is subject to catastrophe, but not if the bipolar field is completely closed under the same specified photospheric conditions. In order to investigate the influence of the photospheric magnetic conditions on the catastrophic behavior of this system, we expand upon the 2.5 dimensional ideal magnetohydrodynamic (MHD) model in Cartesian coordinates to simulate the evolution of the equilibrium states of the system under different photospheric flux distributions. Our simulation results reveal that a catastrophe occurs only when the photospheric flux is not concentrated too much toward the polarity inversion line and the source regions of the bipolar field are not too weak; otherwise no catastrophe occurs. As a result, under certain photospheric conditions, a catastrophe could take place in a completely closed configuration whereas it ceases to exist in a partially open configuration. This indicates that whether the background field is completely closed or partially open is not the only necessary condition for the existence of catastrophe, and that the photospheric conditions also play a crucial role in the catastrophic behavior of the flux rope system.

Recent HESS observations of the ~200 pc scale diffuse gamma-ray emission from the central molecular zone (CMZ) suggest the presence of a PeV cosmic-ray accelerator (PeVatron) located in the inner 10 pc region of the Galactic Center. Interestingly, the gamma-ray spectrum of the point-like source (HESS J1745-290) in the Galactic Center shows a cutoff at ~10 TeV, implying a cutoff around 100 TeV in the cosmic-ray proton spectrum. Here we propose that the gamma-ray emission from the inner and the outer regions may be explained self-consistently by run-away protons from a single, yet fading accelerator. In this model, gamma rays from the CMZ region are produced by protons injected in the past, while gamma rays from the inner region are produced by protons injected more recently. We suggest that the blast wave formed in a tidal disruption event (TDE) caused by the supermassive black hole (Sgr A*) could serve as such a fading accelerator. With typical parameters of the TDE blast wave, gamma-ray spectra of both the CMZ region and HESS J1745-290 can be reproduced simultaneously. Meanwhile, we find that the cosmic-ray energy density profile in the CMZ region may also be reproduced in the fading accelerator model when appropriate combinations of the particle injection history and the diffusion coefficient of cosmic rays are adopted.

The addition of a 28 m Cherenkov telescope (CT5) to the H.E.S.S. array extended the experiment's sensitivity to lower energies. The lowest energy threshold is obtained using monoscopic analysis of data taken with CT5, providing access to gamma-ray energies below 100 GeV. Such an extension of the instrument's energy range is particularly beneficial for studies of Active Galactic Nuclei with soft spectra, as expected for those at a redshift > 0.5. The high-frequency peaked BL Lac objects PKS 2155-304 (z = 0.116) and PG 1553+113 (0.43 < z < 0.58) are among the brightest objects in the gamma-ray sky, both showing clear signatures of gamma-ray absorption at E > 100 GeV interpreted as being due to interactions with the extragalactic background light (EBL). Multiple observational campaigns of PKS 2155-304 and PG 1553+113 were conducted during 2013 and 2014 using the full H.E.S.S. II instrument. A monoscopic analysis of the data taken with the new CT5 telescope was developed along with an investigation into the systematic uncertainties on the spectral parameters. The energy spectra of PKS 2155-304 and PG 1553+113 were reconstructed down to energies of 80 GeV for PKS 2155-304, which transits near zenith, and 110 GeV for the more northern PG 1553+113. The measured spectra, well fitted in both cases by a log-parabola spectral model (with a 5.0 sigma statistical preference for non-zero curvature for PKS 2155-304 and 4.5 sigma for PG 1553+113), were found consistent with spectra derived from contemporaneous Fermi-LAT data, indicating a sharp break in the observed spectra of both sources at E ~ 100 GeV. When corrected for EBL absorption, the intrinsic H.E.S.S. II mono and Fermi-LAT spectrum of PKS 2155-304 was found to show significant curvature. For PG 1553+113, however, no significant detection of curvature in the intrinsic spectrum could be found within statistical and systematic uncertainties.

Aims: Following the detection of the fast radio burst FRB150418 by the SUPERB project at the Parkes radio telescope, we aim to search for very-high energy gamma-ray afterglow emission. Methods: Follow-up observations in the very-high energy gamma-ray domain were obtained with the H.E.S.S. imaging atmospheric Cherenkov telescope system within 14.5 hours of the radio burst. Results: The obtained 1.4 hours of gamma-ray observations are presented and discussed. At the 99 % C.L. we obtained an integral upper limit on the gamma-ray flux of (E>350 GeV) < 1.33 x 10^-8 m^-2s^-1. Differential flux upper limits as function of the photon energy were derived and used to constrain the intrinsic high-energy afterglow emission of FRB 150418. Conclusions: No hints for high-energy afterglow emission of FRB 150418 were found. Taking absorption on the extragalactic background light into account and assuming a distance of z = 0.492 based on radio and optical counterpart studies and consistent with the FRB dispersion, we constrain the gamma-ray luminosity at 1 TeV to L < 5.1 x 10^47 erg/s at 99% C.L.

We report on the tornado-like evolution of a quiescent prominence on 2014 November 1. The eastern section of the prominence first rose slowly transforming into an arch-shaped structure as high as ~150 Mm above the limb; the arch then writhed moderately in a left-handed sense, while the originally dark prominence material became in emission in the Fe IX 171~{\AA} passband, and a braided structure appeared at the eastern edge of the warped arch. The unraveling of the braided structure was associated with a transient brightening in EUV and apparently contributed to the formation of a curtain-like structure (CLS). The CLS consisted of myriads of thread-like loops rotating counterclockwise about the vertical if viewed from above. Heated prominence material was observed to slide along these loops and land outside the filament channel. The tornado was eventually disintegrated and the remaining material flew along a left-handed helical path of approximately a full turn, as corroborated through stereoscopic reconstruction, into the cavity of the stable, western section of the prominence. We suggest that the tornado-like evolution of the prominence was governed by the helical kink instability, and that the CLS formed through magnetic reconnections between the prominence field and the overlying coronal field.

Studying the temporal variability of BL Lac objects at the highest energies provides unique insights into the extreme physical processes occurring in relativistic jets and in the vicinity of super-massive black holes. To this end, the long-term variability of the BL Lac object PKS 2155-304 is analyzed in the high (HE, 100 MeV < E < 300 GeV) and very high energy (VHE, E > 200 GeV) gamma-ray domain. Over the course of ~9 yr of H.E.S.S observations the VHE light curve in the quiescent state is consistent with a log-normal behavior. The VHE variability in this state is well described by flicker noise (power-spectral-density index {\ss}_VHE = 1.10 +0.10 -0.13) on time scales larger than one day. An analysis of 5.5 yr of HE Fermi LAT data gives consistent results ({\ss}_HE = 1.20 +0.21 -0.23, on time scales larger than 10 days) compatible with the VHE findings. The HE and VHE power spectral densities show a scale invariance across the probed time ranges. A direct linear correlation between the VHE and HE fluxes could neither be excluded nor firmly established. These long-term-variability properties are discussed and compared to the red noise behavior ({\ss} ~ 2) seen on shorter time scales during VHE-flaring states. The difference in power spectral noise behavior at VHE energies during quiescent and flaring states provides evidence that these states are influenced by different physical processes, while the compatibility of the HE and VHE long-term results is suggestive of a common physical link as it might be introduced by an underlying jet-disk connection.

A search for dark matter line-like signals was performed in the vicinity of the Galactic Centre by the H.E.S.S. experiment on observational data taken in 2014. An unbinned likelihood analysis was developed to improve the sensitivity to line-like signals. The upgraded analysis along with newer data extend the energy coverage of the previous measurement down to 100 GeV. The 18 h of data collected with the H.E.S.S. array allow one to rule out at 95% CL the presence of a 130 GeV line (at $l = -1.5^{\circ}, b = 0^{\circ}$ and for a dark matter profile centered at this location) previously reported in Fermi-LAT data. This new analysis overlaps significantly in energy with previous Fermi-LAT and H.E.S.S. results. No significant excess associated with dark matter annihilations was found in the energy range 100 GeV to 2 TeV and upper limits on the gamma-ray flux and the velocity weighted annihilation cross-section are derived adopting an Einasto dark matter halo profile. Expected limits for present and future large statistics H.E.S.S. observations are also given.

Typical solar flares display two quasi-parallel, bright ribbons on the chromosphere. In between is the polarity inversion line (PIL) separating concentrated magnetic fluxes of opposite polarity in active regions (ARs). Intriguingly a series of flares exhibiting X-shaped ribbons occurred at the similar location on the outskirts of NOAA AR 11967, where magnetic fluxes were scattered, yet three of them were alarmingly energetic. The X shape, whose center coincided with hard X-ray emission, was similar in UV/EUV, which cannot be accommodated in the standard flare model. Mapping out magnetic connectivities in potential fields, we found that the X morphology was dictated by the intersection of two quasi-separatrix layers, i.e., a hyperbolic flux tube (HFT), within which a separator connecting a double null was embedded. This topology was not purely local but regulated by fluxes and flows over the whole AR. The nonlinear force-free field model suggested the formation of a current layer at the HFT, where the current dissipation can be mapped to the X-shaped ribbons via field-aligned heat conduction. These results highlight the critical role of HFTs in 3D magnetic reconnection and have important implications for astrophysical and laboratory plasmas.

During their propagation through intergalactic space, ultrahigh energy cosmic rays (UHECRs) interact with the background radiation fields. These interactions give rise to energetic electron/positron pairs and photons which in turn feed electromagnetic cascades, contributing to the isotropic gamma-ray background (IGRB). The gamma-ray flux level generated in this way highly depends upon the UHECR propagation distance, as well as the evolution of their sources with redshift. Recently, the Fermi-LAT collaboration reported that the majority of the total extragalactic gamma-ray flux originates from extragalactic point sources. This posits a stringent upper limit on the IGRB generated via UHECR propagation, and subsequently constrains their abundance in the distant Universe. Focusing on the contribution of UHECR at energies below the ankle within a narrow energy band ($(1-4)\times 10^{18}$eV), we calculate the diffuse gamma-ray flux generated through UHECR propagation, normalizing the total cosmic ray energy budget in this band to that measured. We find that in order to not over-produce the new IGRB limit, a local "fog" of UHECR produced by nearby sources may exist, with a possible non-negligible contribution from our Galaxy. Following the assumption that a given fraction of the observed IGRB at 820GeV originates from UHECR, we obtain a constraint on the maximum distance for the majority of their sources. With other unresolved source populations still contaminating the new IGRB limit, and UHECR above the ankle invariably contributing also to this background, the results presented here are rather conservative.

The supernova remnant (SNR) W49B originated from a core-collapse supernova that occurred between one and four thousand years ago, and subsequently evolved into a mixed-morphology remnant, which is interacting with molecular clouds (MC). $\gamma$-ray observations of SNR/MC associations are a powerful tool to constrain the origin of Galactic cosmic-rays, as they can probe the acceleration of hadrons through their interaction with the surrounding medium and subsequent emission of non-thermal photons. The detection of a $\gamma$-ray source coincident with W49B at very high energies (VHE; E > 100 GeV) with the H.E.S.S. Cherenkov telescopes is reported together with a study of the source with 5 years of Fermi-LAT high energy $\gamma$-ray (0.06 - 300 GeV) data. The smoothly-connected combined source spectrum, measured from 60 MeV to multi-TeV energies, shows two significant spectral breaks at $304\pm20$ MeV and $8.4_{-2.5}^{+2.2}$ GeV, the latter being constrained by the joint fit from the two instruments. The detected spectral features are similar to those observed in several other SNR/MC associations and are found to be indicative of $\gamma$-ray emission produced through neutral-pion decay.

Magnetic flux ropes (MFRs) are one kind of fundamental structures in the solar physics, and involved in various eruption phenomena. Twist, characterizing how the magnetic field lines wind around a main axis, is an intrinsic property of MFRs, closely related to the magnetic free energy and stableness. So far it is unclear how much amount of twist is carried by MFRs in the solar atmosphere and in heliosphere and what role the twist played in the eruptions of MFRs. Contrasting to the solar MFRs, there are lots of in-situ measurements of magnetic clouds (MCs), the large-scale MFRs in interplanetary space, providing some important information of the twist of MFRs. Thus, starting from MCs, we investigate the twist of interplanetary MFRs with the aid of a velocity-modified uniform-twist force-free flux rope model. It is found that most of MCs can be roughly fitted by the model and nearly half of them can be fitted fairly well though the derived twist is probably over-estimated by a factor of 2.5. By applying the model to 115 MCs observed at 1 AU, we find that (1) the twist angles of interplanetary MFRs generally follow a trend of about $0.6\frac{l}{R}$ radians, where $\frac{l}{R}$ is the aspect ratio of a MFR, with a cutoff at about $12\pi$ radians AU$^{-1}$, (2) most of them are significantly larger than $2.5\pi$ radians but well bounded by $2\frac{l}{R}$ radians, (3) strongly twisted magnetic field lines probably limit the expansion and size of MFRs, and (4) the magnetic field lines in the legs wind more tightly than those in the leading part of MFRs. These results not only advance our understanding of the properties and behavior of interplanetary MFRs, but also shed light on the formation and eruption of MFRs in the solar atmosphere. A discussion about the twist and stableness of solar MFRs are therefore given.

The inner region of the Milky Way halo harbors a large amount of dark matter (DM). Given its proximity, it is one of the most promising targets to look for DM. We report on a search for the annihilations of DM particles using $\gamma$-ray observations towards the inner 300 parsecs of the Milky Way, with the H.E.S.S. array of ground-based Cherenkov telescopes. The analysis is based on a 2D maximum likelihood method using Galactic center (GC) data accumulated by H.E.S.S. over the last 10 years (2004-2014), and does not show any significant $\gamma$-ray signal above background. Assuming Einasto and Navarro-Frenk-White DM density profiles at the GC, we derive upper limits on the annihilation cross section $\langle \sigma v\rangle$. These constraints are the strongest obtained so far in the TeV DM mass range and improve upon previous limits by a factor 5. For the Einasto profile, the constraints reach $\langle \sigma v\rangle$ values of $\rm 6\times10^{-26} cm^3s^{-1}$ in the $W^+W^-$ channel for a DM particle mass of 1.5 TeV, and $\rm 2\times10^{-26} cm^3s^{-1}$ in the $\tau^+\tau^-$ channel for 1 TeV mass. For the first time, ground-based $\gamma$-ray observations have reached sufficient sensitivity to probe $\langle \sigma v\rangle$ values expected from the thermal relic density for TeV DM particles.