Detection of Solar Neutrons and Solar Neutron Decay Protons

1. Introduction

In association with solar flares, ions are accelerated to high energies and arrive on the Earth. These high-energy ions are called solar cosmic rays, districting from galactic cosmic rays. However, it still needs to be clarified what energies the galactic cosmic rays are accelerated, where they are accelerated, and their acceleration mechanism. The acceleration mechanism for the solar cosmic rays produced by solar flares remains as a puzzle to be understood. Solar flares are used to be classified into two types: the impulsive flare and the gradual flare. For more details on this classification, refer to the paper written by Don Reames in this special issue [1,2].

Solar cosmic rays produced by the gradual flares have been attributed to the gradual acceleration of ions by their reciprocal motion across the shock inside the Coronal Mass Ejection (CME) that is observed in association with the event. However, the hypothesis of acceleration by shock waves fails to explain the observed short-term acceleration of energetic ions in impulsive flares. To illustrate this phenomenon, Holman, Litvinenko, Somov proposed the electric field acceleration model in the 1995 [3,4,5], which produces high-energy ions in a shorter duration. However, there is still no clear observational evidence supporting this model.

Fisk has proposed the wave-riding theory in order to explain the low-energy abundance of ³He and Fe ions [6,7]. However, this model is classified into a 2nd-order Fermi acceleration mechanism [8], which means it can’t quickly accelerate particles to high energies. Therefore, it is crucial to capture as many solar cosmic rays associated with the impulsive flare as possible. To do this, we must detect high-energy solar cosmic rays and simultaneously must observe the X-ray, gamma-ray and ultraviolet images on the solar surface. This is what we refer to as the solar version of the multi-messenger astronomy.

Observation of high-energy solar particle beams requires a great deal of ingenuity due to the presence of a magnetic field between the Sun and Earth, which is created by the solar dynamo activity and has a strength of approximately 5 nT near the Earth. The magnetic field causes the effect of the trajectory of charged particles from the Sun to the Earth. They are transported in the Parker field and transferred to the Earth not straightforwardly. Therefore, this leads to an orbit dependence for charged particles on the arrival time of the particles, longer than the light that travels straightforwardly. As a result, these data cannot be compared with the optical phenomena such as gamma-rays, ultraviolet rays and X-rays on the solar surface. To avoid the time lag problem, we used to compare the data with optical phenomena on the solar surface in case high-energy neutral particles such as gamma-rays and neutrons are detected. Therefore, the detection of these neutral particle rays is expected to clarify the acceleration processes of protons and electrons in the impulsive flare.

Neutrons have mass, which means they take longer time to travel through solar-terrestrial space than light. However, this issue can be resolved by measuring the neutron’s energy. Figure 1 displays the relationship between the energy of a neutron and the flight time traveling through 1 astronomical unit (1 au) of space. For example, a neutron with an energy of 100 MeV will lag behind light by approximately 660 seconds, whereas a neutron with an energy of 400 MeV will have a delay time of 200 seconds. At an energy level of 1000 MeV, the time difference between the neutron and light is only one minute. When the detector threshold is set at 30 MeV, which will be discussed in the next chapter, the delay time is 1,460 seconds (approximately 24 minutes), and at 40 MeV, it turns out as 1,280 seconds (approximately 21 minutes). If solar neutrons are produced on the solar surface in a $\delta \left(t\right)$-functional time distribution, they should be detected by our instruments within about 20 minutes.

However, during their journey over 1 au, neutrons decay into protons, electrons, and antineutrinos. Therefore, we must take this effect into account when determining the flux. For example, at an energy level of 70 MeV, 70 % of neutrons decay during their 1 au flight, and the probability of survival is around 30 %. On the other hand, at an energy level of 1 GeV, 92.5 % of neutrons do not decay and reach the Earth.

Unlike high-energy gamma-rays, neutrons are only produced by collisions of ions with solar atmosphere. Therefore, information on pure ions can be obtained from them, but the above correction is associated. Gamma-rays are also generated by electron bremsstrahlung and also by the decay of neutral mesons (${\pi }^0$) produced by the nuclear collisions of neutrons with the solar atmosphere. Hence, it will lead to an unclear problem of whether gamma-rays are produced by ions or electrons.

Figure 2 provides an example of the first detection of solar neutrons [9]. The Gamma-Ray Spectrometer onboard the Solar Maximum Mission satellite detected the signal, and the neutron monitor located at Jungfraujoch high altitude confirmed the increase of neutrons after one minute later. The peak at 11:43:30 UT was caused by a gamma-ray flash, while the signal arrived one minute later and is estimated to be 1 GeV with the energy of neutrons.

In Figure 3, the event of May 24, 1990, is presented. The Soviet satellite GRANAT received a strong gamma-ray flash due to this event, which had a gamma-ray peak of 10.6-109.5 MeV at 20:48:10 UT. At the same time, a neutron monitor located in North America recorded an arrival of massive neutrons with a rise time at 20:49 UT, with a time difference of approximately one minute [10,11]. Moreover, an accelerated charged ion beam approached the Earth’s vicinity 18 minutes later at 21:07 UT.

Figure 4-a and b show examples of observations received in Mt. Norikura [12,13]. The observation dates are June 4 1991 (Figure 4a) and June 6 1991 (Figure 4b), respectively. The figures show, from top to bottom, the time distribution (1-minute values) of neutrons above 40 MeV observed by the 1 m² Solar Neutron Telescope, the 36 m² Muon Telescope, the Neutron Monitor [14,15], the CGRO satellite, and the 80 GHz radio waves observed by the Nobeyama Solar Radio Telescope. It shows that high-energy electrons were accelerated for a short time on the solar surface and that protons and ions were simultaneously accelerated to high energies. The second one, Figure 4b, shows that the high-energy electrons and ions were again accelerated in the X10 class flare two days later. The fifth graph from the top is the results observed by the neutron monitor located at Haleakala, Hawaii [13]. At the Mt. Norikura neutron telescope, a significant signal increase is observed in the highest energy channel (> 390 MeV). The Nobeyama data shows that the acceleration continued for about five minutes. It would be impressive to know that the difference between the rise time of the Haleakala data and the peak time of the Nobeyama radio signal is one minute.

On September 7th, 2005, two examples were received in the American continent: one at Mt. Chacaltaya in Bolivia (Figure 5a, on the left side) and the other at Mt. Sierra Negra in Mexico (Figure 5b, on the right side). In these figures (Figure 5a and Figure 5b), the start time of ion acceleration is illustrated, which was assumed to be at 17:36:40 UT, based on the X-ray data above >50 keV from the GEOTAIL satellite. The time is shown by the green dashed line in the figures. Then, the energy spectrum was obtained [16]. Although the instantaneous acceleration model cannot explain a part of the tail, a weak signal with an intensity of about 1/10 is recorded for about 10 minutes on the X-ray production time profile. Likely, ions were also produced during this time profile and generated weakly.

On November 7th, 2004, a solar neutron signal was detected over the American continent [17]. However, this event differs from the event on September 7th, 2005. The excess observed in Bolivia may be produced by the solar neutrons as the event of 2005. On the other hand, the event observed at Mt. Sierra Negra is significantly different from the 2005 event. We present the data observed by Mt. Sierra Negra in Figure 6. The data showed an increase in the signal that lasted for 78 minutes. However, if they were produced by solar neutrons, signals are expected to last less than 24 minutes. Therefore, different interpretations are necessary for the two data observed in Bolivia and Mexico. In this paper, we have a working hypothesis that the data observed at Mt. Sierra Negra are due to Solar Neutron Decay Protons. Hereafter, we may call them as SNDP.

This paper is organized as follows. Chapter 2 presents the detection methods and sensitivities of the instruments for solar neutron detection installed in Bolivia, Mexico and Finland. Chapter 3 presents the converting process of the signals observed at Mt. Chacaltaya in Bolivia, into differential spectra, and Chapter 4 describes the flux observed by the detector located at Mt. Sierra Negra. In chapter 5, we explain the general situation of the heliosphere at 16:00 UT on 7 November 2004. Chapter 6 presents the results of the anti-proton orbit analysis launched from 20 km above the ground. These results of the analysis show that the solar neutron decay proton model is the most likely probable model to explain the excesses observed at Mt. Sierra Negra. Chapter 7 presents the results of analysing the signals received by the GOES satellite and the neutron monitor located at Oulu in Finland. After discussions in chapter 8, we finally summarize these stories in chapter 9 and present conclusions.

2. Solar Neutron Detectors

In the early days of cosmic-ray research, passive methods had been used to detect energetic particles, such as GLEs, emitted by solar flares. For this purpose, a stable and long-term operating instrument is traditionally used. Initially, the ionization chamber was selected as the instrument [18], but it was later replaced by the neutron monitor for continuous observation of cosmic rays [19,20,21]. In 1970, for the long-term stable observation of cosmic ray intensity, a detector consisting of a 36 m² plastic scintillator was constructed at Mt. Norikura (2,778 m). The detector is possible to measure the arrival direction of the muon with the use of the two-layer structure. Then, another detector to measure the energy deposit of charged particles in the scintillator was carried out for the first time in 1990. In 1973, a carpet array for air shower observations with a container of liquid scintillators was started at Baksan laboratory in the USSR.

In Figure 7, we present a photo of the neutron monitor that consists of the BF₃ gas-filled proportional counter tube in the center. The BF₃ counter is surrounded by the lead and further covered with paraffin on the outside. When incoming neutrons collide with lead nuclei, several neutrons are emitted. The process is known as the nuclear evaporation process. The average energy of these emitted neutrons is around 20 MeV, which is still high for the detection of the signals by the BF₃ counter. Therefore, they are further collided by the hydrogen and the carbon nuclei in the outer paraffin. Then, the thermal neutrons are guided into the proportional counter, where $\alpha$-rays are produced that are detected as a signal by the BF₃ counter.

To determine the detection efficiency of the neutron monitor, we conducted an experiment using the neutron beam from the accelerator at the Centre for Nuclear Physics, Osaka University. We have found the detection efficiency as 0.20 for $E_n$=100 MeV and 0.42 for $E_n$=350 MeV. More details about these beam experiments are published in the paper [22].

When neutrons enter into the plastic scintillator, they interact with protons inside the material through the charge exchange process, a typical feature of nuclear force. This interaction makes neutrons eject forward as high-energy protons. The incident neutron transfers energy into the proton and is captured in the carbon nucleus. If an incident neutron collides with a hydrogen nucleus in the plastic scintillator, a proton is ejected forward. The energy of the incident neutron can be determined by measuring the energy of the forward-ejected proton with emission angle.

Photographs of current equipment are shown in Figure 8 and Figure 9.Figure 8 presents a photograph of the solar neutron detector installed on Mt. Chacaltaya (5,250 m) in Bolivia, where four vessels with an area of 1 m² stacked 40 cm thick plastic scintillator were prepared inside (Figure 8a). Arrival of neutron signals is recognized by the photomultiplier mounted on the top of the vessels; muons above 80 MeV pass through these scintillators without stopping. This character was used for the energy calibration of the detector. Figure 8b shows a picture of the anti-counter, which was prepared to surround the four vessels of the plastic scintillator. Photomultiplier tubes on both sides are used to determine the signals whether the incident charged particles have passed through these anti-counters or not.

Figure 9 represents a solar-neutron telescope that has been installed on Mt. Sierra Negra, Mexico (4,600 m) [23]. The principle of operation is similar to that of the detector in Bolivia. The thickness of the plastic scintillator is 30 cm (Figure 9a). Muons with energies above 60 MeV can pass through the layer of plastic scintillator. The outer part of the plastic scintillator container is covered by an anti-counter gondola made of proportional counters. Figure 9b) illustrates the installation of the proportional counter for the anti-counter.

In contrast to the apparatus at Mt. Chacaltaya, four additional layers of proportional counters are placed under the plastic scintillator box with an area of 4 m² to discriminate the incident direction. The topmost part is covered with a 5 mm thick lead plate designed to convert gamma-rays into electrons and distinguish them from neutrons. This is because the intensity of background cosmic rays is almost equal to that of muons and electrons at an altitude of 4,600 m. This reduces the background noise from gamma-rays.

The neutron detection efficiency of the instrument was also measured with the use of a neutron beam from the accelerator at the Centre for Nuclear Physics, Osaka University. The values are 0.11 at $E_n$=100 MeV and 0.21 at $E_n$=400 MeV. Details are available in [24]. The detection efficiency of neutrons in the higher energy range than 400 MeV has also been determined by Monte Carlo simulation [25].

The photograph of the 12NM64-type neutron monitor located in Mt. Norikura Cosmic Ray Observatory is shown in Figure 7. The monitor successfully received solar neutrons in association with solar flares on 4 June 1991 and 6 June 1991, and it is still in operation. Figure 7b) shows the part of the power supply to the BF₃ counter. This photograph shows the signal extraction section.

3. Flux of Solar Neutrons observed at Mt. Chacaltaya

For the understanding the event observed at Mt. Sierra Negra, It is important to understand the solar neutron event observed at Mt. Chacaltaya. This chapter focuses on the determination process of the observed solar neutron flux at the Chacaltaya observatory.

There were important data on the hard X-rays obtained by the detectors on board the satellite for the events that occurred on May 24, 1990, September 7, 2005, June 4, 1991, and June 6, 1991. However, in the present case of the November 7, 2004 event, we will use X-ray data obtained from the GOES satellite to estimate the start time of particle acceleration. The green plot in Figure 10 represents the derivative of the GOES X-ray flux. Assuming that the ions were rapidly accelerated at the first peak, which occurred at 15:46:00 UT when the derivative shows its first maximum, and taking into account the fact that 1 GeV protons lag behind light by about 1 minute, the 20-minute excess seen in Figure 11 ➁ would be consistent with the hypothesis that the solar neutrons created this excess. The statistical significance of the increment is 5.0$\sigma$.

According to the data represented by the green line in Figure 1, it can be inferred that neutrons with energies exceeding 40 MeV reach the Earth within 1,250 seconds, which is approximately 20 minutes. The observed time of the increase is consistent with the anticipated 20-minute growth in the 40 MeV threshold energy channel. This assumption is also corroborated by the rise time in the channel for neutrons above 120 MeV and 160 MeV, as illustrated in Figure 12. The increase of these channels occurs approximately 10-12 minutes from 15:47:00 UT and is statistically significant at 3.7 $\sigma$ and 3.2 $\sigma$, respectively. It is observed that neutrons at $E_n$=120 MeV are delayed by 10 minutes when compared to light. The increased times of these channels are in the time distribution expected from the assumption that neutrons were created on the solar surface at 15:46:00 UT. On the other hand, the assumption that ions were accelerated to high energies at either the 15:53 UT or 15:57 UT peak in Figure 10 does not explain the time distribution of neutrons observed at Mt. Chacaltaya. Therefore, this hypothesis shall not be adopted.

Assuming that neutrons are emitted simultaneously from the solar surface, it’s possible to draw an energy spectrum from the flight time difference in arrival time between the Sun and Earth. The result is shown in Figure 13. However, Figure 13 presents the number of neutrons ($\Delta I$) incident in the energy interval ($\Delta E$), divided by the energy width ($\Delta E$), which hasn’t yet been corrected for the energy dependence of the neutron detection efficiency. Two further corrections are necessary: one for the energy dependence of the detection efficiency of the detector itself and another factor for the attenuation of neutrons in the atmosphere. The energy dependence of the detection efficiency was measured using the accelerator beam at the Centre for Nuclear Physics, Osaka University, and the value was applied [24]. For the energy dependence of neutron attenuation in the atmosphere, the value calculated by Shibata was used[26], which agrees well with that calculated using GEANT4 [23]. These two corrections were made to Figure 13, and the resulting differential spectrum is shown in Figure 14.

The observed fluxes at Mt. Chacaltaya atmospheric top are $4\times {10}^5/\left({\mathrm{m}}^2\mathrm{MeV}\right)$ and $12/\left({\mathrm{m}}^2\mathrm{MeV}\right)$ for $E_n$=100 MeV and $E_n$=1,000 MeV, respectively. When approximated by a power function, it can be approximated by $E^{-3.1}\mathrm{d}E$ in the range of $E_n=100-1000$ MeV. The integrated spectrum is shown in Figure 15: the flux at Mt. Chacaltaya atmospheric top is estimated as $4\times {10}^5/{\mathrm{m}}^2$ above 100 MeV and $2.8\times {10}^3/{\mathrm{m}}^2$ above 1,000 MeV. The top of the atmosphere has a flux of 400,000 high-energy neutrons per m². The next chapter will attempt to interpret the phenomenon using these values.

4. Excess observed at Mt. Sierra Negra

This chapter will discuss the results of observation at Mt. Sierra Negra, one of the main subjects of this paper. In the results of Mt. Sierra Negra (Figure 6), we notice two typical features: (a) Firstly, the excess of the events started almost at the same time with nearly the same intensity. The threshold energy of the detector at Mt. Chacaltaya is set at over 40 MeV, while the detector at Mt. Sierra Negra is set at above 30 MeV. The actual intensities at Mt. Sierra Negra and Mt. Chacaltaya were 2,300/5 min. and 2,550/5 min. respectively. However, there was a difference between the atmospheric pressure between the two sites with approximately $217\mathrm{g}/{\mathrm{cm}}^2$, as shown in Figure 16. According to the calculations on the attenuation of solar neutrons in the atmosphere[26], the fluxes at Mt. Sierra Negra are expected to be attenuated by 1/10 and 1/7 for $T_n$ = 100 MeV and 200 MeV neutrons, respectively. However, the observed fluxes of both sites were observed with almost the same intensity. This observation result contradicts the prediction of the simulation.

(b) The second question arises from the duration of the event. The threshold energy of the scintillator at Mt. Sierra Negra is set at over 30 MeV (S1 channel). This means when neutrons are produced impulsively on the solar surface with any energy spectrum, they will reach the detector within 1,500 s (25 min) (as shown by the orange line in Figure 1. However, the excess of the signal continues to 55 minutes. The observed results presented in (a) and (b) may reject a hypothesis that solar neutrons produced the incremental component for 55 minutes at Mt. Sierra Negra. These facts give rise to a serious question of what caused the excess at Mt. Sierra Negra.

Therefore, we propose here that the excess at Mt. Sierra Negra was produced by protons resulting from the decay of solar neutrons. However, the expected energy of these protons is only a few GeV, which raises two issues: (1) can such low-energy protons penetrate the magnetosphere? (2) Can they reach the instruments located at high mountains without being attenuated by the Earth’s atmosphere?

To investigate the first issue, anti-protons were launched from 20 km above Mt. Sierra Negra to see if they could reach the leading edge of the magnetosphere (magnetopause). The second issue was investigated through the Monte Carlo simulation using GEANT4. These results are reported here. The Solar Neutron Decay Proton (SNDP) hypothesis can solve the duration problem of (b), because the neutron signal, which has been changed to a charged proton, is trapped in the magnetic field of the heliosphere (10 nT at present case) and transported to the tip of the Earth’s magnetosphere, which may extend the reception time.

The cut-off rigidity ($R_{\mathrm{c}}$) of Mexico City was predicted as 8.61 GV [27]. However, previous studies by Smart, Shea, and Flückiger have shown that charged particles with energies below the rigidity ($R_{\mathrm{c}}$) can also intrude [28]. Recent observations by PAMELA and AMS satellites indicate that the proton energy spectrum in the 1 GeV-11 GeV range reflects the cosmic ray flux in the energy region above rigidity. Nonetheless, the low-energy component is also observed to penetrate the region below rigidity[29,30,31,32].

Antiprotons were launched 20 km above the observation point, changing emission angle by one degree in both the zenith and azimuthal angles of $0^{\circ }\sim {90}^{\circ }$ and $0^{\circ }\sim {360}^{\circ }$, respectively. A total of 32,761 ($181\times 181=32,761$) antiprotons were launched with different energies. We then studied whether these antiprotons could pass through the magnetosphere and reach the magnetopause. The results of this study are presented in Figure 17. Taking account of an observational condition that protons incident at 20 km above the sky from a horizontal direction will not reach the observation point, antiprotons incident from 20-40 degrees from the zenith are marked as white circles (○). The black circles (●) represent the antiprotons entering above the Earth’s atmosphere within 20 degrees from the zenith. The symbols (●) represent nearly vertical incidence (${\theta }_{\mathrm{NS}}<{20}^{\circ }$ and ${\theta }_{\mathrm{EW}}<{20}^{\circ }$). The open triangles present the proton flux from the heliocentric side (dayside) (∆). The calculation results in Figure 17 show that protons above 5 GeV can penetrate the magnetosphere. It was also indicated that the flux at $E_p$=7 GeV was reduced to about 1/10 of that at $E_p$=20 GeV. These results are consistent with the observations of PAMELA [31,32] and AMS [30]. These attenuation rates in the magnetosphere were applied to the incident SNDP.

To determine whether protons entering the atmosphere above the observation point could reach the observation point, we had to conduct an investigation. We found that incident protons were energy-fragmented and absorbed by nuclear interactions in the atmosphere. According to Figure 18, the signal of low-energy protons can reach the altitude of Mt. Sierra Negra. If the energy of the incident proton is $E_p$=7 GeV, the energy of the parent incident particle is propagated into neutrons and gamma-rays ($E_n$, $E_{\gamma }>$30 MeV). The ratio is 0.7:0.3, but they are expected to be detected as a single signal. In other words, this suggests that a proton of $E_p$=7 GeV incident perpendicularly above the atmosphere can be detected at Mt. Sierra Negra, at an altitude of 4,600 m, with a detection efficiency of almost 100 %. However, the incidence at larger angles is more strongly attenuated, as shown in Figure 19. The attenuation rate was determined through GEANT4 calculations.

To summarize, the Earth’s magnetic field causes the signal to weaken by a factor of 1/100 during the passage in the magnetosphere. However, in the atmospheric passage, protons with $E_p$=7 GeV have almost no reduction effect. Suppose the ions were accelerated to 10 GeV or higher during the 2004.11.7 flare. In that case, when they struck the solar surface, neutrons were produced with energies of a few GeV, and they were emitted into solar-terrestrial space. The protons that resulted from their decay were then received at Mt Sierra Negra. The next chapter will investigate if the flux can be explained by the SNDP assumption without contradictions.

5. Solar-Terrestrial Environment around 2004.11.7

As discussed in the previous chapter, the signal detected at Mt. Sierra Negra might have resulted from the decay of solar neutrons. In this chapter, we investigate whether the observations made at Mt. Chacaltaya and Mt. Sierra Negra can be reconciled with this hypothesis. To do so, we need to know the condition of the plasma in the heliosphere around 16:00 UT on November 7, 2004.

The GOES satellite data obtained on November 7 2004, is shown in Figure 20. The neutron monitor data, displayed in the bottom graph, reveals that the intensity of the neutron monitor decreased in sync with the sudden magnetic field change at around 18:30 UT. This phenomenon is known as Forbush Decrease (FD), which is informed by the arrival of Coronal Mass Ejection (CME) at the Earth. It’s important to note here that the X2.0 flare at 15:42 UT did not create this FD, but it was caused by the M9.3 solar flare that occurred on the solar surface a day earlier at 00:34 UT on November 6. The FD associated with the CME generated by the 15:42 UT flare reached the Earth’s vicinity at 19:20 UT on November 9, and there is observational evidence to support this feature. The radio stars are used by the solar wind detector located at the ground in order to know the plasma motion of the heliosphere. Figure 21 illustrates the observation data for (22:00 UT on November 7, 2004 - 07:00 UT on November 8, 2004) and (22:00 UT on November 8, 2004 - 07:00 UT on November 9, 2004) detected by the Nagoya University solar wind detector, respectively. The latter data in Figure 21(b) clearly shows the high-speed (1,000 km/s) flow of the CME rope generated by the X2.0 flare.

Figure 22 depicts the Sun-Earth space around 16:00 UT. ACE satellite measured the internal magnetic field of CME1, which showed about 150 nT when it passed the L1 point [33]. The width of the rope of CME1, which can be determined by the transit time, helps to determine the maximum radius of the gyroscopic motion of the charged particles. Based on this, the maximum momentum of the charged particles trapped in the CME loop can be calculated at approximately P = 60 GeV/c. This information is crucial for the following discussion.

If the charged particles have smaller momentum than 60 GeV/c, they will be reflected by “the wall of the magnetic loop,” or some ions will be trapped inside the magnetic loop [34]. In the case of the first SNDP event in 1981, the neutrons from the Sun decayed to protons in a “free space”, i.e.,in the case of the first SNDP found, there was a neutron horizon, but no CME rope. Therefore, the heliosphere was more straightforward [35,36,37,38].

The speed of the magnetic rope of CME1 is calculated. It was created at 00:34 UT on November 6th, and they passed the Earth at 19:00 UT on November 7th. The time of flight is $1.5\times {10}^5$ seconds. Dividing the distance from the sun to the Earth (1 au = $1.5\times {10}^{11}$ meters) by the time of flight gives the propagation speed of the CME, which is precisely as 1,000 km/s. This speed is classified as the high-speed solar wind. This implies that at 15:42 UT, CME1 had not yet reached the vicinity of the Earth, and its tip was still about 0.066 au or $1\times {10}^{10}$ meters away from the Earth’s sun-facing side. Figure 22 illustrates this situation.

6. Fluxes between Solar Neutrons and Solar Neutron Decay Protons

This chapter determines at first (a) the expected flux of SNDPs flying to the Earth’s atmosphere. Then (b) determine the flux of SNDP received at Mt. Sierra Negra. After (c) we search for the conditions under which charged particles can pass through the magnetosphere. Then (d) In order to know actual number of entrance of SNDPs in the magnetosphere, we will estimate possible entrance area on the surface of the heliopause and (e) we will examine the matching condition between incoming protons and out-going anti protons at the heliopause. The key parameter may be the momentum vector of both particles. Then (f) we will estimate actual area of the entrance window that satisfies these conditions (1) $P_x$ > 0, (2) $\theta$ < 40^∘, and (3) $P_x$ > $P_y$.

(a)Let us determine the number of neutrons that decay in space with an energy of a few GeV. The probability that a solar neutron with $T_n$=6 GeV will decay while travelling $1\times {10}^{10}$ m is represented by the blue dotted line in Figure 23, and the value is estimated as 0.00510. Additionally, Figure 23 shows the expected values of the integral spectra of solar neutrons with $E_n$=1-10 GeV. The predicted values for the slopes of -2.75, -3.75 and -4.75 are shown, assuming the integral spectrum is expressed by a power function. As demonstrated in Figure 17, SNDP below 5 GeV cannot pass through the magnetosphere, and they do not reach the Earth. The expected solar neutron decay proton flux is estimated at approximately 0.11 protons per unit area (/m²) at the tip of the magnetosphere. This low flux of solar neutrons is comming from the reason that only 0.5% of neutrons at $T_n$ = 6∼7 GeV decay, because of the short neutron flight distance, such as 0.066 au.

(b) Next, we challenge a subject why such a weak signal is yet to be detected. Before we delve into this matter, let’s determine the flux of events observed at Mt. Sierra Negra. The average 5-minute reading of the channel (S1) with threshold energy above 30 MeV is 834,000 events. Meanwhile, an excess was observed in 11 bins (55 min), as depicted in Figure 6. The total excess component was 23,900 events. Since the total instrument area at Mt. Sierra Negra is 4 m², the flux per unit area (m²) was approximately 6,000 events/m². In comparison with the anticipated value of the extended Mt. Chacaltaya, the integral spectrum is estimated as 0.1 /m² around 7 GeV, and the observed value at Mt. Sierra Negra is approximately 60,000 times larger than the expected value. To understand this discrepancy with present hypothesis, analyzing the orbits of antiprotons launched 20 km above the observation point on the Earth may hold the key.

(c) To understand this problem, we launched antiprotons at various angles 20 km above Mt. Sierra Negra. We aimed to study the behavior of the antiprotons and determine where they would reach the magnetopause. The results of our analysis indicate that when the energy of the antiprotons exceeds 6 GeV, a significant number of them can reach the magnetopause. This is illustrated in Figure 24. For instance, when the energy of the antiprotons is 7 GeV, 0.36 % of them can reach the magnetopause. The coordinates of the antiprotons are distributed over a wide area of the magnetosphere tip, which is about 8 $R_{\mathrm{E}}$. This wide distribution may contribute to the boosting effect of the signal. It means protons at the end of the magnetosphere from the outside may be connected to the “several penetration routes” to the Earth, if they reach this region. Our antiproton calculations indicate multiple points where this is the case. It would be adequate to call it “funnel effect” or “Horn effect”. This effect may enhance very few fluxes of the high-energy SNDPs.

Figure 25 shows the distribution of arrival points for the antiproton orbit with $E_p$=6 GeV that meet the following conditions. Firstly, the antiproton must be oriented towards the Sun, meaning the momentum’s X component must be positive ($P_x>0$). Here, we are discussing the GSE coordinate. Secondly, the launch angle must be within 40 degrees from 20 km above the Earth’s surface in order to ensure that the incident signal of protons should not be absorbed by the Earth’s atmosphere before they reach the instruments. Thirdly, we imposed another condition, $P_x>P_y$, which ensures that the SNDP is trapped and transported by the solar wind. On the day of the observation, the 2D map of the solar wind shows that it was crossing the orbit of the Earth at about 50^∘ (Figure 26). Since the seed of solar wind has v=400 km/s by Parker’s equation $tan\theta =r\omega /v$[39], we will get $\theta ={50}^{\circ }$.

(d) Let’s consider that the Solar Neutron Decay Protons (SNDP) arrive uniformly at the front edge of the Earth’s magnetosphere (magnetopause), which is approximately 8 times the radius of the Earth (8 $R_{\mathrm{E}}$). Suppose we estimate the surface area of the magnetosphere’s leading edge when observed from the sun. In that case, we can find the total surface area as $S=\pi {\left(8R_{\mathrm{E}}\right)}^2$, where $R_{\mathrm{E}}$ is the Earth’s radius since solar neutrons don’t arrive on the magnetosphere’s night side (magneto-tail) and are therefore excluded from the calculation. Then, the calculated total area of the surface is estimated as $8.2\times {10}^{15}{\mathrm{m}}^2$. Based on the previous analyses, the expected value of SNDP is ∼0.1 neutrons/m². Multiplying this value by the surface area, we can estimate that the total flux of SNDP arriving at the tip of the Earth’s magnetosphere is $8.2\times {10}^{14}$ neutrons (for $\gamma$= -2.75). Although this number is significant, not all SNDPs can reach the Earth’s surface, and only a fraction is expected to enter the atmosphere. However, the initial values are large enough to boost SNDP (boosting factor).

(e) We estimate the likelihood of an antiproton penetrating the magnetosphere. First, we consider the probability of an antiproton emitted within 40 degrees from the zenith at 20 km altitude and 6 GeV whether they arrive at the magnetopause or not. From the simulation, this probability is estimated as only $1.4\times {10}^{-3}$. The next factor is to consider if the escape direction of antiprotons from the magnetosphere matches with the direction of the solar wind arrival, which must align with the angle of travel of the SNDP. The angle $\theta$ from the Earth’s orbital plane (Y-axis) must be at least ($3\pi /4\sim \pi$) radians to satisfy this restriction on the momentum space. Figure 25 suggests only one point may meet this condition. Furthermore, the injection conditions depend on the coupling between the Parker field and the magnetic field of the magnetosphere. Assuming one state-of-the-art example, a damping factor 1/32761 (=$3.1\times {10}^{-5}$) could be adopted for the incident particle.

(f) The dotted line in Figure 25 represents the area that satisfies conditions (1), (2), and (3). The size of this region is approximately $2{\left(R_{\mathrm{E}}\right)}^2$, equivalent to $8\times {10}^{13}{\mathrm{m}}^2$. Assuming a power of $\gamma$= -2.75, the expected average number of 0.11 SNDPs per unit area results in a flux of arrival protons of $9.2\times {10}^{12}$. Since there is one possible candidate of antiproton orbit among these SNDPs that matches the arrival direction of the Parker field in the momentum space, the results of calculation for -2.75, -3.75 and -4.75 may provide the following results:

$(1.4\times {10}^{-3})\times 0.075\times (3.1\times {10}^{-5})\times (9.2\times {10}^{12})=3.0\times {10}^4$ for $\gamma =-2.75$;

$(1.4\times {10}^{-3})\times 0.075\times (3.1\times {10}^{-5})\times (1.6\times {10}^{12})=6400$ for $\gamma =-3.75$;

$(1.4\times {10}^{-3})\times 0.075\times (3.1\times {10}^{-5})\times (2.7\times {10}^{11})=890$ for $\gamma =-4.75$.

To get the final result, you need to multiply 0.075, the factor coming from the detection efficiency when the entrance angle is tilted from the zenith (∼0.5) and taking into account the detection efficiency of the detector itself (∼0.15). If the spectrum is soft above 1 GeV, the observed value of 6,000/m² at Mt. Sierra Negra can be approached.

7. Comparison with Oulu data and GOES data

In the previous chapter, we presented the idea that solar neutron decay protons could be responsible for the observed increase in the signal at Mt. Sierra Negra recorded immediately after the X2.0 flare. In this chapter, we will delve into the 5.2 $\sigma$ increase signal detected simultaneously on the neutron monitor at Oulu, Finland, which is located near the Arctic Circle, as well as the 34-minute increase observed on the GOES 12 satellite.

The neutron monitor 9NM64 is located in Oulu, Finland. The geographical coordinate of Oulu is 65^∘N 25.5^∘E, with an elevation of 15 m above sea level and a Rigidity of 0.78 GV [27]. Due to its proximity to the Arctic Circle, low-energy protons can enter the monitor without being affected by the geomagnetic field. Figure 27 shows the observed 1-minute values, which indicate an increase in counts starting at 15:59 UT, reaching a peak at 16:00 UT, and ending at 16:04 UT. During these 4 minutes, there was an increase of 1,030 counts, 6.5 $\sigma$ higher beyond the background fluctuation of 157 counts/minute. When converting this excess to per unit area for comparison, the increase was estimated as ($28.5\pm 2$) pcs/m² min.

Let us consider an SNDP with $E_n$=3 GeV. According to Figure 23, its decay probability is 0.01. The predicted flux from Figure 23 is 280 pcs/m², which means an expected value of 2.8 /m² is given, taking account of the decay probability. This value is ten times higher than the observed value. However, the flux of protons is attenuated when they pass through the atmosphere. Figure 28 shows that the attenuation rate is about 1/50. Considering this rate, the intensity differs from the observed value by a factor of about 500. Nonetheless, the amplification effect discussed in the previous Mt Sierra Negra event is at work.

Figure 29 shows the distribution of 400 MeV antiprotons reaching at 15 $R_{\mathrm{E}}$ when launched from above Oulu. We can recognize an incident window on the night side (the blue circle). The area of the window at 8 $R_{\mathrm{E}}$ is estimated to be $0.024R_{\mathrm{E}}^2=1\times {10}^{12}{\mathrm{m}}^2$.

In other words, the 1 m² measuring instrument on the ground is spread over $1\times {10}^{12}{\mathrm{m}}^2$ at 8 $R_{\mathrm{E}}$. This may act as a “boosting factor”. On the other hand, due to the limitation of the coordinate and momentum space to satisfy the penetrating condition, the “attenuation factor” arises as to be $5\times {10}^{-10}$ in total. Therefore, the product can be estimated to be about ∼ 500 times $(1\times {10}^{12}\times 5\times {10}^{-10}=500)$. For a detailed process, please refer to another paper [17].

In this section, we will discuss the excess observed in GOES 10, 11, and 12, as shown in Figure 30 (a) and (b). When we compare the data between GOES 12 and GOES 10 and 11, we notice a clear difference near 16:00 UT, the time of the flare. On the other hand, the increase at 18:30 UT is due to the arrival of CME1 near the Earth, the same as the increase at ➂ in Figure 11. GOES 10 was located at 75^∘E, while GOES 11 was situated at 104^∘E. This position was exactly local time 11 am (LT) and 9 am (LT) respectively. GOES 10 was looking almost directly at the Sun, while GOES 11 viewed the Sun from an angle of 45^∘. The flux above 10 MeV was recorded by GOES 11 as 1.2 /${\mathrm{cm}}^2·\mathrm{s}·\mathrm{str}$. GOES 12, on the other hand, was situated at 135^∘E and the local time was 7:00 am (LT), just morning. The flux of GOES 11 was converted for comparison with the data obtained from the Alpine instrument. The intensity was estimated to be $2.3\times {10}^7/{\mathrm{m}}^2$ as an integrated value. Comparing this value with the integrated value in Figure 15, it is consistent with the flux above 20 MeV when the line of integration curve in Figure 15 is extended from 100 MeV to lower energies. However, Figure 15 does not include the probability of neutrons decaying during their flight between the Sun and Earth. It is only the intensity above the Earth.

The GOES 11 was received the signal between 4 and 14 minutes after 15:46:00 UT at the flare start time at different channels (> 10 MeV). On the other hand, the 4-9 MeV channel of GOES 12 received the incleased signal for 34 minutes. Assuming that solar neutrons were produced instantaneously, the energies of each neutron are estimated to be in the range of 70-300 MeV (for GOES 11) and 18-300 MeV (for GOES 12), respectively; the survival probability of a 70 MeV neutron is 0.3, while that of a 20 MeV neutron is 0.1. Considering these facts, it is highly likely that the GOES 11 and 12 satellites observed solar neutrons and SNDPs. If it was the signal of SNDP, it should have lasted for 55 minutes, the same as Sierra Negra observation time, so we assume that the signal recorded in the detectors of GOES 12 satellite was induced by a charged particle after a solar neutron flew in and was converted to a proton in the satellite’s instrument. In the following table (Table 1), we summarize the results.

8. Discussions

We have been discussing the solar flare on November 7th, 2004. There are reported very few examples of SNDP observations. The first discovery of SNDP was made on the flare of June 3rd, 1981 [35,36,37,38], followed by events on May 24th, 1990[10,11], and October 19th, 1989 [40]. Present November 7th, 2004 event is the fourth case. In the case of the October 19th, 1989 event, solar neutron decay electrons were also observed and analyzed in detail by one of the present authors [41].

Two of these cases were detected by satellite observations. On the other hand, the October 19, 1989 event was recorded by only the neutron monitors located at the surface level. The authors pointed out that the increase of the neutron monitors located at Deep River and Goose Bay was induced by the arrival of solar neutron decay protons because these detectors are located at the asymptotic directions of the anti-proton trajectory at the magnetopause. On the other hand, if the abrupt excess after two minutes of the X13 big flare at 13:00 UT was produced by the direct incidence of solar neutrons, any positive signals may be recorded in the neutron monitor located at Tsumeb, Namibia (17.35^∘E, 19.12^∘S). Still, no excess was recorded in the detector. Therefore, the authors concluded that the abrupt increase of the signals recorded in the neutron monitors must be produced by the arrival of solar neutron decay protons. The November 7 2004, event may be the second example recorded by the detectors on the ground.

The November 7th, 2004 event provided significant new information that differs from previous events. By comparing the SOHO satellite images taken before and after the 195 nm wavelength at 15:36 UT and 15:48 UT, we can see the formation of an arch connecting the two active regions. This is shown by the white arches in Figure 31. The formation of this arch is significant because it is believed that the ions accelerated to high energies above this loop, hit the solar surface, and produced solar neutrons and gamma-rays.

Figure 32(a) presents the solar surface around 16:06 UT and 16:10 UT observed by the RHESSI sectrometer in the band of 25-50 keV. The brightest point (282, 101) on the solar coordinate coincides with the white arrow shown in Fugure 31 (b). It would be interesting to know that the hottest point moved to (322,74) on the time 16:25-16:29 UT. Figure 32(b) was observed by the channel of energy band of 100-300 keV.

Figure 33 displays the energy spectrum detected by the KONUS satellite [42]. The KONUS satellite was primarily designed to detect gamma-ray bursts, but intense X-rays from the Sun also trigger the sensor. This is how it was received. You may notice line gamma-rays of 2.223 MeV produced by the annihilation process of neutrons with hydrogen targets forming the deuterium element. Also, you may notice another lines of gamma-rays in the region of 4.438 MeV and 6.129 MeV. They may be produced by the de-excitation process of carbon and oxygen nuclei in the solar atmosphere.

9. Summary and Conclusion

On November 7th, 2004, a giant flare of magnitude X2.0 occurred in the solar atmosphere, as a result, ions were accelerated to high energies greater than 10 GeV. The accelerated primary protons collided with material on the solar surface, producing line gamma-rays and neutrons. These particles were detected by instruments aboard the GOES and KONUS satellites as neutrons and line gamma-rays, respectively.

When neutrons with a few GeV travel 0.06 au in the solar-terrestrial space, they decay and produce secondary protons. These secondary protons arrive at the magnetopause, and some of them enter the magnetosphere and arrive over the Earth’s surface. Even though the flux of the secondary protons is low, the Funnel or Horn effect, which collects them over a wide area of the magnetosphere, results in significant received signals. This information implies that ions on the solar surface get accelerated to over 10 GeV, associated with an impulsive flare.