Open Access
Issue
A&A
Volume 694, February 2025
Article Number A58
Number of page(s) 7
Section The Sun and the Heliosphere
DOI https://doi.org/10.1051/0004-6361/202453144
Published online 31 January 2025

© The Authors 2025

Licence Creative CommonsOpen Access article, published by EDP Sciences, under the terms of the Creative Commons Attribution License (https://creativecommons.org/licenses/by/4.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

This article is published in open access under the Subscribe to Open model. Subscribe to A&A to support open access publication.

1. Introduction

Solar flares, the most powerful release of energy in the entire Solar System, emit radiation across the entire electromagnetic spectrum, from radio, to visible, up to hard X-rays (HXRs) and γ-rays (Fletcher et al. 2011). Through the release of “free” magnetic energy, high-energy electrons, protons, and ions are accelerated in the low corona, and can either travel along magnetic field lines toward the solar surface or escape into interplanetary space. The downward flare-accelerated particles, upon interaction with the much denser chromosphere and photosphere, lose most of their energy by collision (i.e., heating) or by radiation. In the HXR and γ-ray ranges, part of this radiation is generated by semi-relativistic and relativistic (≳300 − 400 keV) electrons through continuum bremsstrahlung emission. Additionally, protons and ions produce spectral lines in the γ-ray range through nuclear reactions and particle decay processes. For instance, upon collision with ambient ions, different nuclear de-excitation lines are produced from ∼0.8 MeV to ∼10 MeV, such as 12C at 4.438 MeV and 16O at 6.129 MeV. These broad de-excitation lines merge into a quasi-continuum dominating the electron bremsstrahlung continuum (Vilmer et al. 2011). Because the de-excitation after collision has a timescale of 10−12 s (Vilmer et al. 2011), these lines are known as “prompt” lines, as there is no a measurable delay and they are thus emitted almost instantaneously. The neutron-capture line, emitted at 2.223 MeV upon the capture of a neutron by ambient hydrogen (forming deuterium), has a delay between the generation of the neutron and its capture (e.g., Murphy et al. 2003; Vilmer et al. 2011). The reason for this delay is the thermalization of the neutrons, because the scattering cross section of a flare-produced neutron is larger than that of a neutron-capture (Hua & Lingenfelter 1987). Therefore, a flare-produced neutron first undergoes a number of collisions (i.e., thermalization) before being captured and consequently emitting the 2.223 MeV line. The timescale of the generation of the neutron-capture line is about 100 s (Murphy et al. 2007). However, shorter timescales have been found, such as in Lysenko et al. (2019), which found a delay of about 10 s. The capture probability is highest for neutrons produced in the photosphere, where they are least able to escape or decay before they are thermalized and captured (Hua & Lingenfelter 1987; Rieger 1989).

The history of HXR imaging observations of flares is marked by a successful heritage of space-based imaging telescopes. Some examples, among others, include the Hard X-ray Telescope (HXT; Kosugi et al. 1992) aboard the Yohkoh mission, the Reuven Ramaty High-Energy Solar Spectroscopic Imager (RHESSI; Lin et al. 2002), and the currently operational Spectrometer/Telescope for Imaging X-rays (STIX; Krucker et al. 2020), as well as the Hard X-ray Imager (HXI; Zhang et al. 2019) aboard the Advanced Space-based Solar Observatory (ASO-S Gan et al. 2019). Imaging γ-rays in flares, instead, is notoriously much more challenging. These challenges arise from the combination of technological difficulties in detecting and imaging γ-rays, which require significantly more mass in the space instrument than HXRs, and counting statistics issues due to the steep decrease in the spectrum from HXRs to γ-rays. This decrease amounts to approximately two to three orders of magnitude in flux for every tenfold increase in energy. Consequently, only flares that are highly efficient in accelerating high-energy protons and ions can produce detectable γ-rays, typically corresponding to the most energetic events (e.g., Shih et al. 2009). Currently, no operational solar γ-ray imaging telescope exists for routine flare observations. Nevertheless, efforts in developing future instruments are being undertaken in order to address this gap in our observational capabilities (e.g., Ryan et al. 2023).

Despite the challenges surrounding the detection of γ-rays, protons and ions constitute key components in flares. For instance, based on the analysis of a few dozen γ-ray producing flares, it was found that the energy contained in > 1 MeV ions lies in the range of 1029 − 33 ergs (e.g., Ramaty et al. 1995, 1997; Miller et al. 1997; Murphy et al. 1997; Vilmer et al. 2011). This shows that, beyond the electrons probed in the HXR range, diagnosing ions in the γ-ray range in flares is of the utmost importance, as the energy content of flare-accelerated electrons and ions are comparable (Emslie et al. 2012).

The first-ever and, to date, the only γ-ray images of solar flares were obtained by RHESSI. Even though RHESSI observed countless energetic flares, imaging of γ-ray emission was only possible for about half a dozen events (Lin et al. 2003; Hurford et al. 2006; Krucker et al. 2008; Vilmer et al. 2011). Nevertheless, many important results have been obtained and we refer to Vilmer et al. (2011) for a comprehensive review. The first-ever γ-ray image of a flare, of the X4 GOES class, occurred on July 23, 2002, and is reported in Lin et al. (2003). In this paper, the authors showed that the centroid of the ion-produced 2.223 MeV neutron-capture line emission is located approximately 20 ± 6 arcsec away from the electron bremsstrahlung source locations, suggesting that the acceleration and/or propagation of the ions is different from that of the electrons. This result was further corroborated by Hurford et al. (2006), who analyzed the source location of the 2.223 MeV neutron-capture γ-ray line in three large RHESSI flares of 2003 (October 28, October 29, and November 2), and reported a systematic displacement from the corresponding 0.2–0.3 MeV electron-bremsstrahlung emission footpoints of 10 to 25 arcsec. According to the standard flare model, this spatial separation remains puzzling, as both electrons and protons are assumed to be accelerated in the same region. Therefore, the question of why they propagate along different magnetic loops (and consequently precipitate at different locations) is still unresolved and could not be investigated with new observations. In this paper, we revisit the X25 GOES class flare that occurred on October 28, 20031, by incorporating additional information obtained from visible light observations.

Visible light observations have become crucial in studying the precipitation sites of flare-accelerated electrons, as over the past decade significant progress has been made in this context. Visible light enhancements associated with flares, also referred to as white-light (WL) emission, are observed as either continuum enhancements (e.g., Hudson et al. 1992, 2006; Martínez Oliveros et al. 2011, 2012; Krucker et al. 2015; Kleint et al. 2016; Kuhar et al. 2016; Jurčák et al. 2018) and/or strong emission from photospheric spectral lines (e.g., Abramenko & Baranovsky 2004; Yurchyshyn et al. 2017)2. By analyzing a flare close to the limb, Martínez Oliveros et al. (2012) found that the WL and HXR centroids closely match, which suggests that the observed WL emission mechanism is directly linked to the energy deposition by flare-accelerated electrons. A similar result was reported by Krucker et al. (2015), who analyzed three limb events. In addition, Kuhar et al. (2016) conducted a statistical analysis revealing that HXRs and WL emissions are correlated not only spatially but also in terms of peak time, intensity, and energy deposition, which highlights the crucial role of flare-accelerated electrons in WL formation. Therefore, WL observations in flares provide valuable insights into electron precipitation sites. To the best of the authors’ knowledge, no observational studies have yet investigated if there is any relationship between the source location of γ-rays (i.e., the proton precipitation sites) and WL flare ribbons.

In this paper, we revisit the extraordinary X25 GOES class flare that occurred on October 28, 2003, which produced the most intense γ-ray signal ever imaged in solar flares. This event has been extensively investigated in numerous studies (e.g., Bianda et al. 2005; Hurford et al. 2006; Kiener et al. 2006; Su et al. 2006; Zharkova & Zharkov 2007; Krucker et al. 2008; Trottet et al. 2008; Kurt et al. 2024). Our objective is to analyze the γ-ray source locations of the 2.223 MeV neutron-capture line and contextualize it with WL emissions. This approach allows us to examine, for the first time, the correlation between WL emission and the proton precipitation sites. In Sect. 2, the data analysis and an overview of the event are presented. In Sect. 3, we report results and discussions. Finally, our conclusions are drawn in Sect. 4.

2. Selected event and data analysis

2.1. Event overview

To investigate the relationship between the γ-ray source locations of the 2.223 MeV neutron-capture line and WL emission, we focused on the X25 GOES class flare of October 28, 2003. This event, originating from AR 10486 near the disk center, produced the most intense γ-ray signal ever imaged in a solar flare. This flare has been observed by RHESSI. Although the spacecraft passage through the South Atlantic Anomaly prevented the rising and peak of the impulsive phase from being recorded, the exceptional magnitude of the flare allowed for significant HXR and γ-ray flux observations during the decaying impulsive phase. RHESSI fully captured this latter phase, as is illustrated in Fig. 1.

thumbnail Fig. 1.

Time histories of the October 28, 2003, flare from RHESSI (top) and INTEGRAL/SPI (bottom). The gray area highlights the RHESSI coverage of the flare, which started right after the end of the peak of the impulsive phase at 11:06:20.

The spectrometer instrument (SPI; Vedrenne et al. 2003) aboard the International Gamma Ray Astrophysics Laboratory (INTEGRAL; Winkler et al. 2003) provided comprehensive time profile observations for the entire flare event. In this study, we use SPI time profiles as a reference for the impulsive phase. The Global Oscillation Network Group telescope (GONG; Harvey et al. 1996) in Cerro Tololo, Chile, also captured the impulsive phase. GONG observations are crucial for this study, as they provide information on the WL emission during both impulsive and decay phases, and thus reveal the electron precipitation sites. Fig. 2 presents an overview of the GONG continuum intensity observations. The pre-flare subtracted images clearly show the evolution of flare ribbons.

thumbnail Fig. 2.

Overview of the GONG intensity observations of the October 28, 2003, flare. The leftmost panel shows the continuum intensity image of the region where the flare happened. From the second to the fifth panel, we show the time evolution (around the peak of the impulsive phase at 1 minute time resolution) of the pre-flare subtracted continuum intensity images. From the pre-flare subtracted images, we can clearly see the flare ribbons. The percent level of the maximum increase in WL emission due to the flare is 35% of the pre-flare intensity.

2.2. Data analysis

RHESSI data were acquired using the RHESSI software within IDL SolarSoftware (SSWIDL). Image reconstruction was performed using the CLEAN algorithm (Högbom 1974) with natural weighting. For HXR images (200–300 keV), detectors 5 through 9 were used, while γ-ray images (2218–2228 keV) employed only detectors 6 and 9, which are associated with the thickest grids. Both HXR and γ-ray image reconstruction was done by considering counts from the rear detector segment only.

The INTEGRAL/SPI data were taken from Kiener et al. (2006). The GONG continuum intensity observations have been downloaded from the NSO GONG data archive3. These data product are available at one minute resolution. To obtain pre-flare subtracted continuum intensity images, we subtracted the observations taken one minute prior to flare onset from all subsequent observations. The Transition Region and Coronal Explorer (TRACE; Handy et al. 1999) 195 Å extreme ultraviolet (EUV) observations were downloaded from the virtual solar observatory portal4.

3. Results and dicussions

3.1. White-light and extreme ultraviolet observations

The aim of this subsection is to compare EUV and WL images at two distinct times, as is illustrated in Fig. 3. The upper panel displays EUV and WL images near the peak of the impulsive phase (11:04:05 UT). Notably, the bright EUV ribbons align well with the GONG WL contours. This combined EUV and WL emission clearly indicates the electron precipitation site, where energy is deposited into the lower atmosphere, chromosphere, and/or photosphere. To aid visualization, semicircles perpendicular to the solar surface have been drawn connecting the ribbons, highlighting the newly created flare loops.

thumbnail Fig. 3.

EUV and WL images of the flare at two different time instances. In both panels, the TRACE 195 Å image is plotted as the background, with red contours overlaid to indicate the flare ribbons derived from pre-flare subtracted GONG intensity images (see Fig. 2). The two different times correspond to a time near the peak of the impulsive phase (top) and a time during the decay phase (bottom). To represent the flare loops connected to the ribbons, we draw as a proxy semicircles perpendicular to the solar surface, where yellow semicircles denote flare loops formed during the impulsive phase, while the green ones represent flare loops created during the decay phase. The selection of the footpoints for the proxy loops was driven to have several loops along the ribbons connecting the WL sources without much twist in them. The selection is somewhat arbitrary and intended as a visual guide.

A similar pattern is observed at the later time (11:09:02 UT). The energy release and particle acceleration in this flare is long-lasting, which creates successively larger loops, as is illustrated in the bottom panel of Fig. 3, as well as new ribbons at different locations. These ribbons are associated with newly reconnected field lines (depicted as green loops), which align well with both EUV and WL emissions. Furthermore, the loops formed near the impulsive phase peak (shown in orange) now exhibit coronal emission in EUV, which have cooled down to become visible in the TRACE 194 Å passband.

The EUV and WL observations of this event showcase a nice example that aligns with the standard flare reconnection model. In this model, magnetic field lines initially interact with field lines at a specific altitude, followed by subsequent reconnection at higher altitudes (e.g., Qiu et al. 2009). In addition, Fig. 3 illustrates how GONG WL observations effectively agree with the electron precipitation site, corroborating what is reported in previous studies (e.g., Krucker et al. 2015; Kuhar et al. 2016).

3.2. Hard X-ray and γ-ray observations

In light of the recent advances in understanding WL generation and its correlation with HXRs, the comparison between WL and γ-ray emissions is central to this paper. By revisiting the γ-ray imaging of the October 28, 2003, event, we gained important insights into two distinct concepts: the challenge of comparing images at different energies reconstructed over extended time periods, and the varying contributions of different flare phases to the generation of the neutron-capture line. Along with the results obtained by revisiting the γ-ray imaging of this event, these two concepts are discussed in the following.

3.2.1. The extended time integration of previous studies

As a first step, we checked that the current version of the RHESSI software produces consistent images for RHESSI rear-segment images, as was reported in Hurford et al. (2006). The RHESSI imaging software had improved over the years, with the main change being the introduction of visibilities (Piana et al. 2007). Due to these improvements, the reconstructed images are not exactly the same as reported in Hurford et al. (2006), but the centroid location of the γ-ray sources are the same within the errors bars obtained from forward fitting. With the current software, by repeating the same analysis as was done in Hurford et al. (2006), the resulting center positions for the two γ-ray sources reported are: (xE, yE) = (−141″, −388″) for the eastern source and (xW, yW) = (−75″, −357″) for the western source. These positions align closely with the centroid positions (one sigma uncertainty of ±5″) reported by Hurford et al. (2006): ( x E H , y E H ) = ( 145 , 380 ) $ (x_{\mathrm{E}}^{\mathrm{H}},\mathit{y}_{\mathrm{E}}^{\mathrm{H}})=(-145\prime\prime,-380\prime\prime) $ for the eastern source and ( x W H , y W H ) = ( 74 , 354 ) $ (x_{\mathrm{W}}^{\mathrm{H}},\mathit{y}_{\mathrm{W}}^{\mathrm{H}})=(-74\prime\prime,-354\prime\prime) $ for the western source. The differences in coordinate values are up to 8″, within the error range.

The main difference of our analysis compared to previously published RHESSI γ-ray imaging results is in the selected time intervals and the new consideration of WL images. To enhance counting statistics as much as possible, the papers published in the past reconstructed images averaged over the entire duration of the flare. Comparing reconstructed RHESSI images of such long-time integrations is not as straightforward as it seems. The limited dynamic range and resolution of these reconstructions can be misleading if the relative time evolution in HXRs and γ-rays is different. To illustrate this, we assume a simple source geometry of a single, common source for electron and proton precipitation on a flare ribbon that moves from A to B during the integration time. Integrating over a long duration, the derived RHESSI source location is the flux-weighted average along the line between A to B. If the γ-rays decay slower than the HXRs, the gamma-ray source centroid is closer to B than for the HXR centroid. Hence, it appears that electrons and protons are precipitating on different locations, but this could be solely because the time profiles decay differently. Hence, the previously published results that electron and protons precipitate at different locations only holds if the electron and proton precipitation rates have the same time evolution. As both HXRs and γ-rays are used for imaging reconstruction, and their profiles are obviously different (refer to Fig. 1), diagnostics of source locations is, particularly for long integration times, ambiguous.

3.2.2. Different contributions to the neutron-capture line

In addition to the challenge of comparing sources at different energies reconstructed over extended time periods, there is the intrinsic complication that the neutron-capture line at 2.223 MeV is not emitted instantaneously. The neutron-capture line emission is not instantaneous due to the thermalization process preceding capture, resulting in a time evolution with a relatively long decay. Consequently, careful evaluation of the contributions from various flare phases is essential during the reconstruction of γ-ray neutron-capture line images. Therefore, in our analysis, we propose a new approach that uses three main improvements:

  • Much shorter integration times are used;

  • Selected time intervals derived from observations of prompt γ-ray lines in combination with simulations of neutron-capture process (Murphy et al. 2003);

  • Comparison of the γ-ray images to WL observations by GONG at the time of the peak of the prompt γ-ray emission.

Murphy et al. (2003) showed simulations of the actual delay of the neutron-capture line at 2.223 MeV produced by a δ-function injections of protons. The obtained time profiles depended on several parameters such as the spectral index of the accelerated protons, the level of pitch angle scattering, the relative 3He abundances, and the flare location relative to the observer. The top panel of Fig. 4 shows a typical profile taken from Murphy et al. (2003). To qualitatively understand how the 2.223 MeV profile is produced, we used the observed prompt line profile from INTEGRAL/SPI as a measure of when protons enter the chromosphere or photosphere and convolved it with the typical profile from Murphy et al. (2003). The resulting profile is shown in the second panel of Fig. 4 in black compared with the observed 2.223 MeV signal from INTEGRAL/SPI. The two profiles show a similar behavior in time, confirming that we have a qualitative agreement. By varying the parameters of the simulation, a potentially better fit could be found. For example, selecting a simulation from Murphy et al. (2003) with a sharper, early peak will likely give a better agreement around peak time where the currently used profile smears out individual peaks. However, such an analysis is not the focus of our paper. Here, we have only used the derived profiles to check how well the 100 s approximation of the generation of the 2.223 MeV signal holds for this flare. The black curve is separated into two components: the blue curve is the part of the 2.223 MeV signal that is produced by protons that precipitate during the part of the decay time when RHESSI was observing. The red curve is the signal from protons that precipitate at times when RHESSI was not observing. Hence, 100 s after the start of the RHESSI observations, the 2.223 MeV signal contains about 20% of the signal from the peak and 80% from the decay. This shows that the 100 s delay is a good approximation, considering the limited dynamic range of RHESSI imaging in the γ-ray range. Nevertheless, it is also clear that this approximation does not necessarily hold for all flares. For a flare with a with a weaker decay signal relative to the flare peak, the 2.223 MeV signal after 100 s would still be dominated by the flare peak.

thumbnail Fig. 4.

Relative timing of proton precipitation and creation of the 2.223 MeV neutron-capture line emission. (Top panel) Time profile of the 2.223 MeV line emission, assuming a δ-function injection of protons at time zero (for details, see Murphy et al. 2003). (Second panel) Time profiles of the INTEGRAL/SPI prompt lines, which is the same as the black curve in Fig. 1. Red indicates periods when RHESSI was not observing (peak of the impulsive phase), while blue highlights RHESSI observation times (decaying impulsive phase). (Third panel) Time profiles of the INTEGRAL/SPI 2.223 MeV neutron-capture line (black histogram). Red indicates the contribution to the 2.223 MeV neutron-capture line of the times prior to RHESSI observations, and hence the peak of the impulsive phase, while blue indicates the contribution from the decaying impulsive phase. (Bottom panel) Relative contribution to the 2.223 MeV neutron-capture line from the peak of the impulsive phase as a function of time, calculated from the ratio of the black curve in the middle panel divided by the red one. The vertical lines give the edges of the time intervals used for the images presented in Fig. 5.

3.2.3. Hard X-ray and γ-ray imaging results

Our previously introduced qualitative model shows that for the first 100 s after the start of the RHESSI observations, the 2.223 MeV signal mainly derives from protons that are accelerated at the peak of the flare (bottom panel of Fig. 4). Hence, by comparing the 2.223 MeV image reconstructed from this first part with the WL sources at the peak time, we can directly compare the precipitation site of flare-accelerated protons relative to the flare ribbons.

The left panel in Fig. 5 shows the comparison, revealing that the 2.223 MeV image produced mainly by protons from the peak time of the flare nicely agrees in location with WL. The limited dynamic range does not allow us to make a very detailed comparison; nevertheless, the strongest 2.223 MeV source is clearly associated with the easternmost part of the WL flare ribbon. The HXR images for this time interval are of course missing, but the general close association of WL and HXR flare ribbons sources reported in the literature strongly suggests that electron and proton precipitation happens on flare ribbons. Therefore, there is no evidence from RHESSI imaging that electrons and protons are precipitating at different locations.

thumbnail Fig. 5.

Multiwavelength imaging overview of the October 23, 2003, event at two different times. In both panels, the TRACE 195 Å image is plotted as a green background, with gray contours representing flare ribbons derived from pre-flare subtracted GONG intensity images. Red contours indicate RHESSI HXR images, while blue contours show RHESSI γ-ray images. The left panel depicts a time near the impulsive phase peak, and the right panel illustrates a time during the decaying impulsive phase. The time integration intervals for the RHESSI γ-ray images have been adjusted based on the contribution of the different flare phases to the neutron-capture line (see Sect. 3 for details).

The right panel of Fig. 5 is similar to the original figure from Hurford et al. (2006) comparing a time-shifted γ-ray image with a HXR image, but for a significant shorter integration time of 150 s compared to the 21.67 minutes (1300 s) considered in Hurford et al. (2006). Similarly to during the peak time, we can clearly see agreement between γ-ray and WL images. For the shorter integration time, the ribbon motion is slow compared to the angular resolution of the HXR and γ-ray image shown in Fig. 5. From TRACE observations, the ribbon motion within the time interval is between 5 and 7 arcsec. Hence, the comparison between HXR and γ-ray centroids is not significantly affected by the different time evolution at these two energy ranges. Comparing the centroids at the two energy ranges give the same location, within errors. A forward fit to the visibilities gives a difference of up to 1.5 sigma (∼8 arcsec) for both eastern and western sources. Hence, within errors, HXR and γ-ray centroids are the same and associated with the WL ribbons.

3.3. Comparison with other γ-ray flares

The event considered in this paper is the sole occurrence in the entire history of solar observations for which it is possible to reconstruct γ-ray images at distinct time intervals. Consequently, this analysis cannot be replicated for the other events documented in the literature where γ-ray imaging was possible. For these latter cases, only extended integration times are possible, and the published results, including the separation of the γ-ray and HXR centroid positions, are influenced by the different time evolutions.

For all the other flares (i.e., the July 23, 2002, event in Lin et al. (2003), the October 29, 2003, and November 2, 2003, events in Hurford et al. (2006), and the January 20, 2005, event in Krucker et al. 2008), GONG provides WL observations. We examined the WL ribbon evolution and found that, for all events, the flare ribbons intersect the centroid position of all γ-ray sources at some point during their temporal evolution, within one-sigma uncertainty of the reconstruction. This finding aligns with our previously described interpretation.

4. Conclusions

In this study, we have revisited the exceptional X25 GOES class flare of October 28, 2003, observed by RHESSI, which produced the most intense γ-ray signal ever imaged of a solar flare. Despite this flare containing more total 2.223 MeV counts than all other RHESSI flares combined, only the decay phase just after the flare peak time was observed. By optimizing the integration time for image reconstruction using simulations of the neutron capture process by Murphy et al. (2003), we could establish two time ranges for γ-ray imaging: one representing the flare peak time and a second during the start of the decay phase. We found a clear spatial correlation between the γ-ray sources and the HXR and WL images. We present two key conclusions: (1) the electron and proton precipitation sites coincide, at least within the RHESSI resolution, and (2) protons therefore contribute to WL formation, a factor often overlooked in previous studies. Both electrons and protons are found to be correlated with the WL and EUV flare ribbons. The source motion seen in WL and EUV between the flare peak time and the start of the decay is also seen in the γ-ray images, corroborating that protons as well as electrons deposit their energy in the flare ribbons.

The previously reported spatial separation between electron and proton precipitation sites is a result of a difference in the relative time evolution of the γ-ray and HXR signals in combination with the motion of the sources as the flare ribbons separate. Hence, the conclusions in Lin et al. (2003) and Hurford et al. (2006) that electrons and protons precipitate on different field lines are an overstatement. While RHESSI shows that protons precipitate into the flare ribbons, the available RHESSI observations do not have the diagnostic capabilities to investigate details of the precipitation sites along each of the flare ribbons. By comparing the WL ribbons with the γ-ray sources of the other four events documented in the literature where γ-ray imaging was possible, we found that these findings are not exclusive to the October 28, 2003, event but applicable for these other events too.

The agreement between γ-ray and WL sources highlights the potential importance of protons in the WL emission mechanism. Previous observations reported a clear correlation between WL flare ribbons and the precipitation sites of electrons. Furthermore, HXR and WL flux are correlated as well with correlation coefficients up to ∼0.7 (Kuhar et al. 2016). Nevertheless, there is a large scatter in the correlation, with individual values differing by up to a factor of three from the fitted value. Part of the scatter could be attributed to protons that contribute to the WL formation process.

In the γ-ray spectrum of solar flares, the neutron-capture at 2.223 MeV is the most prominent line. Therefore, it is the most attractive feature for image reconstruction, since solar flare γ-rays are orders of magnitude less intense than HXRs. Consequently, distinguishing between the contributions of the peak impulsive phase and the decay phase in this reconstruction process is of the utmost importance. This separation prevents averaging visibilities that can differ significantly. As there is currently no solar-dedicated γ-ray imager in space, future space missions should carefully consider this aspect in the design phase of the instrument.

This lack of dedicated solar γ-ray imaging capabilities in space underscores the need for new mission concepts to make progress in understanding the importance of energetic protons in the flare energy release process. The balloon payload GRIPS (Duncan et al. 2013) is the only planned imager, currently scheduled for its second flight in 2025 or 2026. Nevertheless, a dedicated space-born instrument needs to be envisioned. The Large Imaging Spectrometer for Solar Accelerated Nuclei (LISSAN; Ryan et al. 2023), as it was originally conceived – to be part of a suite of instruments for an ESA M-class mission concept (SPARK; Reid et al. 2023) – is such a concept.


1

The October 28, 2003, event was previously categorized as of the X17 GOES class. However, considering the recently established correction factor of previous GOES missions (Hudson et al. 2024), the class of this event is estimated to be X25.

2

In this paper, we use the term WL in order to refer to generic photospheric emission enhancement, without differentiating between line or continuum enhancement. Using this term, we want to emphasize the impact of flare-accelerated electrons or ions to the lower atmosphere.

Acknowledgments

The authors acknowledge Albert Shih, Brian Dennis, and Gordon Hurford for the constructive discussions and critical comments, which enhanced the quality of this paper. We also would like to thank Marina Battaglia for commenting on the final draft of the paper. RHESSI was a NASA small explorer mission lead and operated by UC Berkeley in close collaboration with GSFC. GONG data were acquired by instruments operated by NISP/NSO/AURA/NSF with contribution from NOAA. AFB is supported by the Swiss National Science Foundation Grant 200020_213147.

References

  1. Abramenko, V. I., & Baranovsky, E. A. 2004, Sol. Phys., 220, 81 [NASA ADS] [CrossRef] [Google Scholar]
  2. Bianda, M., Benz, A. O., Stenflo, J. O., Küveler, G., & Ramelli, R. 2005, A&A, 434, 1183 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  3. Duncan, N., Shih, A., Hurford, G., et al. 2013, in Solar Physics and Space Weather Instrumentation V, eds. S. Fineschi, & J. Fennelly, SPIE Conf. Ser., 8862, 88620W [NASA ADS] [CrossRef] [Google Scholar]
  4. Emslie, A. G., Dennis, B. R., Shih, A. Y., et al. 2012, ApJ, 759, 71 [Google Scholar]
  5. Fletcher, L., Dennis, B. R., Hudson, H. S., et al. 2011, Space Sci. Rev., 159, 19 [Google Scholar]
  6. Gan, W.-Q., Zhu, C., Deng, Y.-Y., et al. 2019, Res. Astron. Astrophys., 19, 156 [Google Scholar]
  7. Handy, B. N., Acton, L. W., Kankelborg, C. C., et al. 1999, Sol. Phys., 187, 229 [Google Scholar]
  8. Harvey, J. W., Hill, F., Hubbard, R. P., et al. 1996, Science, 272, 1284 [Google Scholar]
  9. Högbom, J. A. 1974, A&AS, 15, 417 [Google Scholar]
  10. Hua, X. M., & Lingenfelter, R. E. 1987, Sol. Phys., 107, 351 [NASA ADS] [CrossRef] [Google Scholar]
  11. Hudson, H. S., Acton, L. W., Hirayama, T., & Uchida, Y. 1992, PASJ, 44, L77 [NASA ADS] [Google Scholar]
  12. Hudson, H. S., Wolfson, C. J., & Metcalf, T. R. 2006, Sol. Phys., 234, 79 [NASA ADS] [CrossRef] [Google Scholar]
  13. Hudson, H., Cliver, E., White, S., et al. 2024, Sol. Phys., 299, 39 [NASA ADS] [CrossRef] [Google Scholar]
  14. Hurford, G. J., Krucker, S., Lin, R. P., et al. 2006, ApJ, 644, L93 [Google Scholar]
  15. Jurčák, J., Kašparová, J., Švanda, M., & Kleint, L. 2018, A&A, 620, A183 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  16. Kiener, J., Gros, M., Tatischeff, V., & Weidenspointner, G. 2006, A&A, 445, 725 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  17. Kleint, L., Heinzel, P., Judge, P., & Krucker, S. 2016, ApJ, 816, 88 [CrossRef] [Google Scholar]
  18. Kosugi, T., Sakao, T., Masuda, S., et al. 1992, PASJ, 44, L45 [NASA ADS] [Google Scholar]
  19. Krucker, S., Hurford, G. J., MacKinnon, A. L., Shih, A. Y., & Lin, R. P. 2008, ApJ, 678, L63 [Google Scholar]
  20. Krucker, S., Saint-Hilaire, P., Hudson, H. S., et al. 2015, ApJ, 802, 19 [CrossRef] [Google Scholar]
  21. Krucker, S., Hurford, G. J., Grimm, O., et al. 2020, A&A, 642, A15 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  22. Kuhar, M., Krucker, S., Martínez Oliveros, J. C., et al. 2016, ApJ, 816, 6 [Google Scholar]
  23. Kurt, V. G., Veronig, A. M., Fleishman, G. D., et al. 2024, A&A, 686, A195 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  24. Lin, R. P., Dennis, B. R., Hurford, G. J., et al. 2002, Sol. Phys., 210, 3 [Google Scholar]
  25. Lin, R. P., Krucker, S., Hurford, G. J., et al. 2003, ApJ, 595, L69 [NASA ADS] [CrossRef] [Google Scholar]
  26. Lysenko, A. L., Anfinogentov, S. A., Svinkin, D. S., Frederiks, D. D., & Fleishman, G. D. 2019, ApJ, 877, 145 [NASA ADS] [CrossRef] [Google Scholar]
  27. Martínez Oliveros, J. C., Couvidat, S., Schou, J., et al. 2011, Sol. Phys., 269, 269 [CrossRef] [Google Scholar]
  28. Martínez Oliveros, J.-C., Hudson, H. S., Hurford, G. J., et al. 2012, ApJ, 753, L26 [CrossRef] [Google Scholar]
  29. Miller, J. A., Cargill, P. J., Emslie, A. G., et al. 1997, J. Geophys. Res., 102, 14631 [NASA ADS] [CrossRef] [Google Scholar]
  30. Murphy, R. J., Share, G. H., Grove, J. E., et al. 1997, ApJ, 490, 883 [NASA ADS] [CrossRef] [Google Scholar]
  31. Murphy, R. J., Share, G. H., Hua, X. M., et al. 2003, ApJ, 595, L93 [Google Scholar]
  32. Murphy, R. J., Kozlovsky, B., Share, G. H., Hua, X. M., & Lingenfelter, R. E. 2007, ApJS, 168, 167 [CrossRef] [Google Scholar]
  33. Piana, M., Massone, A. M., Hurford, G. J., et al. 2007, ApJ, 665, 846 [CrossRef] [Google Scholar]
  34. Qiu, J., Gary, D. E., & Fleishman, G. D. 2009, Sol. Phys., 255, 107 [NASA ADS] [CrossRef] [Google Scholar]
  35. Ramaty, R., Mandzhavidze, N., Kozlovsky, B., & Murphy, R. J. 1995, ApJ, 455, L193 [NASA ADS] [Google Scholar]
  36. Ramaty, R., Mandzhavidze, N., Barat, C., & Trottet, G. 1997, ApJ, 479, 458 [NASA ADS] [CrossRef] [Google Scholar]
  37. Reid, H. A. S., Musset, S., Ryan, D. F., et al. 2023, Aerospace, 10, 1034 [NASA ADS] [CrossRef] [Google Scholar]
  38. Rieger, E. 1989, Sol. Phys., 121, 323 [NASA ADS] [CrossRef] [Google Scholar]
  39. Ryan, D. F., Musset, S., Reid, H. A. S., et al. 2023, Aerospace, 10, 985 [NASA ADS] [CrossRef] [Google Scholar]
  40. Shih, A. Y., Lin, R. P., & Smith, D. M. 2009, ApJ, 698, L152 [Google Scholar]
  41. Su, Y. N., Golub, L., Van Ballegooijen, A. A., & Gros, M. 2006, Sol. Phys., 236, 325 [NASA ADS] [CrossRef] [Google Scholar]
  42. Trottet, G., Krucker, S., Lüthi, T., & Magun, A. 2008, ApJ, 678, 509 [NASA ADS] [CrossRef] [Google Scholar]
  43. Vedrenne, G., Roques, J. P., Schönfelder, V., et al. 2003, A&A, 411, L63 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  44. Vilmer, N., MacKinnon, A. L., & Hurford, G. J. 2011, Space Sci. Rev., 159, 167 [NASA ADS] [CrossRef] [Google Scholar]
  45. Winkler, C., Courvoisier, T. J. L., Di Cocco, G., et al. 2003, A&A, 411, L1 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  46. Yurchyshyn, V., Kumar, P., Abramenko, V., et al. 2017, ApJ, 838, 32 [NASA ADS] [CrossRef] [Google Scholar]
  47. Zhang, Z., Chen, D.-Y., Wu, J., et al. 2019, Res. Astron. Astrophys., 19, 160 [CrossRef] [Google Scholar]
  48. Zharkova, V. V., & Zharkov, S. I. 2007, ApJ, 664, 573 [CrossRef] [Google Scholar]

All Figures

thumbnail Fig. 1.

Time histories of the October 28, 2003, flare from RHESSI (top) and INTEGRAL/SPI (bottom). The gray area highlights the RHESSI coverage of the flare, which started right after the end of the peak of the impulsive phase at 11:06:20.

In the text
thumbnail Fig. 2.

Overview of the GONG intensity observations of the October 28, 2003, flare. The leftmost panel shows the continuum intensity image of the region where the flare happened. From the second to the fifth panel, we show the time evolution (around the peak of the impulsive phase at 1 minute time resolution) of the pre-flare subtracted continuum intensity images. From the pre-flare subtracted images, we can clearly see the flare ribbons. The percent level of the maximum increase in WL emission due to the flare is 35% of the pre-flare intensity.

In the text
thumbnail Fig. 3.

EUV and WL images of the flare at two different time instances. In both panels, the TRACE 195 Å image is plotted as the background, with red contours overlaid to indicate the flare ribbons derived from pre-flare subtracted GONG intensity images (see Fig. 2). The two different times correspond to a time near the peak of the impulsive phase (top) and a time during the decay phase (bottom). To represent the flare loops connected to the ribbons, we draw as a proxy semicircles perpendicular to the solar surface, where yellow semicircles denote flare loops formed during the impulsive phase, while the green ones represent flare loops created during the decay phase. The selection of the footpoints for the proxy loops was driven to have several loops along the ribbons connecting the WL sources without much twist in them. The selection is somewhat arbitrary and intended as a visual guide.

In the text
thumbnail Fig. 4.

Relative timing of proton precipitation and creation of the 2.223 MeV neutron-capture line emission. (Top panel) Time profile of the 2.223 MeV line emission, assuming a δ-function injection of protons at time zero (for details, see Murphy et al. 2003). (Second panel) Time profiles of the INTEGRAL/SPI prompt lines, which is the same as the black curve in Fig. 1. Red indicates periods when RHESSI was not observing (peak of the impulsive phase), while blue highlights RHESSI observation times (decaying impulsive phase). (Third panel) Time profiles of the INTEGRAL/SPI 2.223 MeV neutron-capture line (black histogram). Red indicates the contribution to the 2.223 MeV neutron-capture line of the times prior to RHESSI observations, and hence the peak of the impulsive phase, while blue indicates the contribution from the decaying impulsive phase. (Bottom panel) Relative contribution to the 2.223 MeV neutron-capture line from the peak of the impulsive phase as a function of time, calculated from the ratio of the black curve in the middle panel divided by the red one. The vertical lines give the edges of the time intervals used for the images presented in Fig. 5.

In the text
thumbnail Fig. 5.

Multiwavelength imaging overview of the October 23, 2003, event at two different times. In both panels, the TRACE 195 Å image is plotted as a green background, with gray contours representing flare ribbons derived from pre-flare subtracted GONG intensity images. Red contours indicate RHESSI HXR images, while blue contours show RHESSI γ-ray images. The left panel depicts a time near the impulsive phase peak, and the right panel illustrates a time during the decaying impulsive phase. The time integration intervals for the RHESSI γ-ray images have been adjusted based on the contribution of the different flare phases to the neutron-capture line (see Sect. 3 for details).

In the text

Current usage metrics show cumulative count of Article Views (full-text article views including HTML views, PDF and ePub downloads, according to the available data) and Abstracts Views on Vision4Press platform.

Data correspond to usage on the plateform after 2015. The current usage metrics is available 48-96 hours after online publication and is updated daily on week days.

Initial download of the metrics may take a while.