Free Access
Issue
A&A
Volume 633, January 2020
Article Number A83
Number of page(s) 17
Section The Sun
DOI https://doi.org/10.1051/0004-6361/201936306
Published online 15 January 2020

© ESO 2020

1. Introduction

The blocking of cosmic rays by the Sun, as first proposed in Clark (1957), is an interesting subject in the investigation of cosmic-ray properties. In particular, it has been predicted that interactions of cosmic rays with the Sun would result in gamma- and neutrino-production that would eventually be detectable (see Seckel et al. 1991). Even cosmic-ray electrons were predicted to produce a halo component in gamma-rays. Today, gamma-ray signatures can indeed be measured: the halo component has been verified at an intensity level as expected by the prediction in Seckel et al. (1991; also see Orlando & Strong 2008; Abdo et al. 2011). The disk component, however, was found to be an order of magnitude higher than expected (Ng et al. 2016), a discrepancy that still could not be resolved in the literature. The neutrino signal that accompanies the disk component has been recalculated over the past few years (Argüelles et al. 2017; Ng et al. 2017; Edsjö et al. 2017) as it is now an important background for dark-matter searches with neutrino telescopes from the direction of the Sun. Nisa et al. (2019) also provides a general discussion of the possibilities of using the Sun as a laboratory for astroparticle physics.

The proper modeling of the cosmic-ray component itself only became important in recent years when ground-based detectors started detecting the shadow of the Sun in cosmic rays produced by their interactions with the solar surface and corona. The Tibet AS-Gamma experiment was the first to demonstrate that the cosmic-ray Sun shadow varies with time, in correlation with the variation of the sunspot number (Amenomori et al. 2013). As the latter is, in turn, correlated with magnetic field strength and structure, these measurements were the first to prove an effect of cosmic-ray deflection in the inner heliospheric field on the Sun shadow. Most recently, a similar study at higher energies ⟨ECR ⟩ ∼40 TeV was performed with the IceCube detector (Aartsen et al. 2019) and a general temporal variation of the cosmic-ray Sun shadow was verified. A correlation with the sunspot number, however, could not be validated yet because the solar cycle during the time of the measurements (2010/2011–2014/2015) was not variable enough to allow for a statistically significant statement (Bos 2017). The HAWC collaboration (Enriquez-Rivera & Lara 2016; Nisa et al. 2017) also observed a deficit of cosmic rays from the direction of the Sun for the years 2013 and 2014 at energies between 1 TeV and 142 TeV. A study of the variability of the Sun shadow could not yet, however, be presented. In the coming years, different experiments are expected to be able to probe the variation of the cosmic-ray Sun shadow with the solar magnetic field across a broad energy range and with a better time resolution.

It is the aim of this paper to investigate which effects are present and how the cosmic-ray shadow varies with time in the ∼5 TeV to ∼316 TeV energy range based on numerical calculations of particle propagation. Our studies are specifically aimed at a quantification of the energy behavior of the Sun shadow. While the Tibet, IceCube, and HAWC measurements present results at different energies, they are not comparable at this point for two central reasons: (1) the Sun shadow is presented on the event level, and therefore, the individual measurements for the shadow cannot be compared directly; (2) measurements are not necessarily performed during the same year and solar cycle. Thus, this study is aimed at carrying out the first theoretical investigation of the energy behavior in order to make a prediction on whether and how observatories could study cosmic-ray propagation in even greater detail. Given the high level of statistics achieved with Tibet, IceCube, and HAWC, this should certainly be possible in the future.

The propagation of cosmic rays can generally be modeled by using either an ensemble- or a single-particle approach. The many-particle approach is typically applied in diffusive regimes in which the gyro-radius rg of the particles is small compared to the system size L, such that a high number of interactions with the turbulent magnetic field for particle trajectories in the considered volume is expected (see, e.g., Schlickeiser 2002). For Galactic cosmic-ray propagation, with a magnetic field strength of ∼μG and assuming a particle energy of 10 TeV, the resulting gyroradius of rg ≈ 10−4 pc is, for example, much smaller than the diameter of the Milky Way of roughly 30 kpc. As a consequence, particles will complete a high number of gyrations and a diffusive description is appropriate. Over recent decades, sophisticated simulation schemes have been developed for Galactic cosmic-ray propagation. In particular, these efforts have resulted in the propagation tool GALPROP (see, e.g., Strong & Moskalenko 1998; Strong et al. 2010; Vladimirov et al. 2011) and newer propagation software, such as DRAGON (Evoli et al. 2008), which is capable of addressing questions of anisotropy, along with PICARD (Kissmann 2014), and CRPropa 3.1 (Merten et al. 2017). For the problem studied in this paper, due to the resulting small number of gyrations within the solar corona, a diffusive description of the propagation of > 10 TeV cosmic rays around the Sun would not be appropriate: the diffusive regime has not been reached in that case and a different approach needs to be applied. In this context, the ballistic (single-particle) approach, in which the equation of motion is solved explicitly for each particle, is suitable. This approach works well for cases in which the gyro-radii are comparable to or larger than the system size, allowing for the computation of the actual trajectories of TeV cosmic rays. A well-structured and comprehensive state-of-the-art tool originally designed for extragalactic propagation is CRPropa (Alves Batista et al. 2013). As we applied different scales in our simulations, we developed a dedicated code for the simulations. In the solar environment, the gyro-radius of rg ≈ 0.5 R (E = 10 TeV, B = 1 G) is on the same order of magnitude as the extent of the studied region of approximately 1 R to 10 R. For the purposes of efficiency, a back-tracing scheme has been set up in which anti-particles are propagated toward the Sun.

This paper is organized as follows: in Sect. 2, some general properties and the modeling of the solar magnetic field are discussed. Afterwards, we describe the details of the numerical propagation, the weighting scheme, and the final analysis quantifying the strength of the Sun shadow (Sect. 3). The results obtained are then presented in Sect. 4, followed by a summary (Sect. 5) and an outlook (Sect. 6).

2. The solar magnetic field

In this section, we summarize some general properties of the solar magnetic field (see Sect. 2.1), followed by the details of the implementation used in this work (see Sect. 2.2).

2.1. General properties

The Sun goes through a cycle of approximately 11 years during which the solar magnetic field varies in strength and structure. At the beginning of each cycle, the magnetic field strength is minimal and the magnetic field near the solar surface can approximately be described by a dipole field which, with increasing radial distance, gets deformed into a radial field by the solar wind. After five to six years, the magnetic field strength reaches its maximum and the dipolar structure disappears. During this period of time, the structure is highly irregular and also changes on the timescale of solar rotation (e.g., Priest 1982). After a total of approximately 11 years, the structure returns to a dipole structure with reversed polarity. Measurements show that the maxima and minima of solar magnetic activity correlate with the maxima and minima of the sunspot number (see, e.g., Hoyt & Schatten 1998; Solanki et al. 2006). Figure 1 shows the sunspot number from 2007 through 2017 over the analyzed time frames. Since our aim is to qualitatively compare our results to the IceCube measurements reported by Aartsen et al. (2019), we define a season as the period consisting of the months November through February of each Antarctic summer. The magnetic field strength at the solar surface is measured using observatories on Earth or in space. Such measurements are called magnetograms and are produced for each Carrington rotation (CR).

thumbnail Fig. 1.

Sunspot numbers taken from SILSO World Data Center (2007–2017) and averaged over the four months of each season are shown for a ten-year time period. The seasons 2007-08 and 2008-09 report a low sunspot number, while in seasons 2011-12, 2012-13, 2013-14, and 2014-15 a higher number is detected, before decreasing again afterwards. This approximately describes the solar activity cycle of 11 years (see, e.g., Solanki et al. 2006).

2.2. Implementation

The Global Oscillation Network Group (GONG, Harvey et al. 1996) provides integral synoptic magnetograms, which we used in order to extrapolate to a full three-dimensional magnetic field using different models. Here we used the Potential Field Source Surface (PFSS) model (Schatten et al. 1969; Altschuler & Newkirk 1969). The magnetic field was computed using the FDIPS software (Toth et al. 2016), which allows the user to define parameters, such as the grid resolution at which the field is to be computed or the so-called source surface radius RSS at which the magnetic field is assumed to be purely radial. More details can be found in such works as Wiengarten et al. (2014). This procedure results in field components on a spherical grid equidistant in the radius r, longitude φ, and the sine of the co-latitude ϑ. Here, we chose RSS = 2.5 R. A typical result for the coronal magnetic field between the solar surface and source surface can be seen in Fig. 2. The magnetic field data used for Fig. 2 has already been employed in Wiengarten et al. (2014), where a different visualization of the magnetic field can also be found (page 8, Fig. 6). As Carrington Rotation 2060 was comparably quiet in terms of solar activity, most field lines appear in an ordered, large-scale structure. It is only close to the surface, where the magnetic field strength becomes significantly larger in some places, that more unordered, small-scale structures are visible. For times of higher solar activity, as, for example, between 2011 and 2015, these unordered, small-scale structures become more relevant regarding cosmic-ray propagation through the solar magnetic field. For r >  RSS the magnetic field is implemented as purely radial and with a magnitude decreasing proportional to the inverse square of the distance. This implementation omits angular components of the magnetic field as present in the Parker spiral model Parker (1958). However, we restricted our study to a region close to the Sun in which these angular components do not contribute significantly to the deflection of the tracked particles (compare Sect. 3.1). We also refrained from modeling the turbulent component of the heliospheric magnetic field since it is not relevant for the > 10 TeV particles studied in this work; it would only become important in the case of a diffusive behavior and as argued in the introduction, we are investigating scales that are comparable to the gyroradii of the particles. Thus, propagation is assumed to be ballistic. An alternative approach for extrapolating from magnetograms to a three-dimensional solar magnetic field is the Current Sheet Source Surface (CSSS) model Zhao & Hoeksema (1995), which additionally uses a so-called cusp surface located at a cusp radius Rcp. Beyond this radius, all field lines are assumed to be open. In both the PFSS and the CSSS model, field lines are assumed to be purely radial at r >  RSS.

thumbnail Fig. 2.

Magnetic field strength on the solar surface and selected magnetic field lines for Carrington rotations 2077 (December 2008 with low solar activity, left) and 2158 (December 2014 close to the solar activity maximum, right), based on GONG magnetogram data and the PFSS model. Field lines are shown between the solar surface and the source surface at 2.5 solar radii, with blue indicating outgoing (Br >  0), brown indicating incoming (Br <  0) polarity. Magnetic field strength on the solar surface is shown using a rainbow color scale. It can clearly be seen that the dipolar structure which typically prevails during solar minimum turns into a more irregular field structure towards solar maximum.

3. Simulation of the Sun shadow

The basic setup of our simulation framework includes the solar magnetic field and anti-particles that are traced back in the vicinity of the Sun. The back-tracing ensures that only those particles that propagate in the direction of Earth after they have traversed the magnetic field are considered (those particles, in turn, are relevant for the cosmic-ray Sun shadow observed at Earth).

The Sun shadow is determined by defining a quadratic, planar injection window with a side length of four solar radii, located at a distance of five solar radii from the Sun’s center and arranged such that it is parallel to the yz-plane of the Heliocentric Earth equatorial (HEEQ) coordinate system, meaning that it is parallel to the Sun’s rotation axis and perpendicular to the intersection of the solar equator and the solar central meridian as seen from Earth (compare Hapgood 1992). For the sake of simplicity, we assumed the solar rotational axis to be perpendicular to the ecliptic, thus neglecting the Sun’s obliquity of 7.25° throughout this paper.

This injection window is uniformly filled with an equidistant, rectangular grid of N = 1002 = 104 starting points of anti-particles. The initial momentum vector of a given anti-particle i points along a straight line from a virtual point on the x-axis at x = 1 AU to the anti-particle’s starting point xi(0), where 1 ≤ i ≤ N. In other words, the point of divergence of the initial directions of all anti-particles is at the position of the Earth. An example of the shadow picture of an anti-proton (|Z| = 1) with an energy E = 40 TeV for December 2014 (Carrington rotation CR2158) is shown in Fig. 3.

thumbnail Fig. 3.

Schematic view of the simulation setup (distances not to scale). Starting points of anti-particles as indicated by the blue bubbles are arranged in the injection window. Timing information for the trajectories of the anti-particles is given in color code. The starting points of those anti-particles hitting the Sun are removed from the plot. In this way, the cosmic-ray deficit caused by the Sun is simulated using back-tracing.

3.1. Computation of particle trajectories

In general, the relativistic equation of motion for a particle of rest mass m under the influence of an external force F is given by

F = d d t ( γ ( v ) m v ) . $$ \begin{aligned} \boldsymbol{F} = \frac{\mathrm{d}}{\mathrm{d} t} \left( \gamma (v) \, m \, \boldsymbol{v} \right). \end{aligned} $$(1)

In this equation, v is the velocity and γ is the Lorentz factor. The motional electric field at radius r is of order rΩB, with the solar angular rotational speed Ω. Comparing this to the (v × B) term due to the Lorentz force shows a ratio of rΩ/c. Since this ratio is on the order of 10−5, we assume that the particles passing through the solar magnetic field are only affected by the Lorentz force. And since the latter cannot change v in magnitude because it is perpendicular to it, the Lorentz factor is constant and thus unaffected by the time derivative. With that, the equation of motion becomes

d v ( t ) d t = q γ m ( v ( t ) × B ) with d x ( t ) d t = v ( t ) , $$ \begin{aligned} \frac{\mathrm{d} \boldsymbol{v}(t)}{\mathrm{d} t} = \frac{q}{\gamma m} \left(\boldsymbol{v}(t)\times \boldsymbol{B}\right) \quad \mathrm{with} \quad \frac{\mathrm{d} \boldsymbol{x}(t)}{\mathrm{d} t} = \boldsymbol{v}(t), \end{aligned} $$(2)

with x and q denoting the location and the charge of the particle. The equation of motion (2) is solved numerically using the so called Boris Push method (Boris 1970). Also see Birdsall & Langdon (1985) for more information. The main advantage of the Boris Push method as compared to other standard methods like the Runge-Kutta-based Cash-Karp method is the conservation of energy for the case that there are no electric fields (which is shown in Qin et al. 2013). Also its relatively short computing time compared to Runge-Kutta algorithms with comparable accuracy is beneficial. Given the initial conditions x(t0) and v(t0), v(t0 + Δt) and x(t0 + Δt) can be determined as follows:

l = q B Δ t 2 γ m , $$ \begin{aligned}&\boldsymbol{l} = \frac{q\boldsymbol{B}\Delta t}{2\gamma m}, \end{aligned} $$(3)

s = 2 l 1 + | l | 2 , $$ \begin{aligned}&\boldsymbol{s} = \frac{2\boldsymbol{l}}{1+|\boldsymbol{l}|^2}, \end{aligned} $$(4)

v = v ( t 0 ) + v ( t 0 ) × l , $$ \begin{aligned}&\boldsymbol{v}^{\prime } = \boldsymbol{v}_{\perp }(t_{0}) + \boldsymbol{v}_{\perp }(t_{0}) \times \boldsymbol{l}, \end{aligned} $$(5)

v ( t 0 + Δ t ) = v ( t 0 ) + v × s , $$ \begin{aligned}&\boldsymbol{v}(t_{0} + \Delta t) = \boldsymbol{v}(t_{0}) + \boldsymbol{v}^{\prime } \times \boldsymbol{s}, \end{aligned} $$(6)

x ( t 0 + Δ t ) = x ( t 0 ) + v ( t 0 + Δ t ) Δ t $$ \begin{aligned}&\boldsymbol{x}(t_{0} + \Delta t) = \boldsymbol{x}(t_{0}) + \boldsymbol{v}(t_{0} + \Delta t) \; \Delta t \end{aligned} $$(7)

(and analogously for all following time steps), where v denotes the component of v perpendicular to B.

3.2. Simulation tests

Various tests of increasing complexity have been conducted prior to the actual runs presented in the ensuing sections in order to confirm the proper operation and reliability of our code. From those tests, we present the noteworthy example of a dipolar magnetic field. This exercise is useful because its results will help us to understand and interpret the ones of the actual Solar magnetic field, for which a dipole is often used as a crude approximation.

Our test setup consists of an ensemble of mono-energetic protons (charge +e, mass mp) being deflected in the magnetic field

B = ( p · r r 3 ) = 3 r ( p · r ) r 5 p r 3 $$ \begin{aligned} \boldsymbol{B} = -\nabla \left(\frac{\boldsymbol{p} \cdot \boldsymbol{r}}{r^3} \right) = \frac{3 \, \boldsymbol{r} \, (\boldsymbol{p} \cdot \boldsymbol{r})}{r^5} - \frac{\boldsymbol{p}}{r^3} \end{aligned} $$(8)

of a dipole with magnetic moment p = pez and strength p= B dip R 3 $ p = B_{\rm dip} \, R_{\odot}^3 $, such that the magnetic field strength at the equator ( x 2 + y 2 = R 2 $ x^2+y^2=R_{\odot}^2 $, z = 0) is exactly equal to Bdip = 1 G. Normalizing all lengths to R and speeds to light speed c, we obtain

d v ¯ d t ¯ = η v ¯ × B ¯ with η = Z e B 0 c R γ m p c 2 20.88 Z B 0 [ G ] E kin [ TeV ] $$ \begin{aligned} \frac{\mathrm{d} \boldsymbol{\bar{v}}}{\mathrm{d} {\bar{t}}} = \eta \, \boldsymbol{\bar{v}} \times \boldsymbol{\bar{B}} \quad \mathrm{with} \quad \eta = \frac{Z \, e \, B_0 \, c \, R_{\odot }}{\gamma \, m_{\rm p} \, c^2} \approx 20.88\,\frac{Z \, B_0 \, [\mathrm{G} ]}{E_{\rm kin} \, \mathrm{[TeV]} } \end{aligned} $$(9)

(in which dimensionless quantities are denoted by a bar and normalization constants with a lower index “0”) as the dimensionless form of Eq. (2). We may, thus, easily identify η, which is proportional to the normalization field strength divided by the particle’s (unchanging) rigidity, as the only external parameter of relevance in this case. For protons in a static field, it is sufficient to vary only the particle energy between different simulations.

Figure 4 shows shadow maps obtained in this way for two specific energies. Somewhat surprisingly, we see that as opposed to what might have been expected in view of previous observations showing the Sun’s shadow diminishing with increased magnetic activity, the shadow’s total area may, in fact, exceed that of the solar disk considerably at certain energies. Figure 5 summarizes the relative shadow size s as a function of the particle energy employed in each simulation. Here the shadow size is determined as the total shadow area as seen in green color in Fig. 4, divided by the area of the solar disk when projected onto the injection window. We see that as expected, this ratio approaches unity for very large energies (or very small field strengths) because trajectories are, then, almost straight lines tracing the geometrical, circular shape of the projected Sun. Interestingly, we may also confirm that the shadow size does in fact increase toward smaller energies, reaching a maximum of about 185 percent of the geometrical disk area before declining again. But it is only at comparatively low energies that this value decreases below unity as the shape of the shadow starts to break apart to ultimately dissolve beyond recognition as it fragments into smaller and smaller “islands”.

thumbnail Fig. 4.

Shadow maps for mono-energetic protons in the magnetic field of a dipole, using the example of Ekin = 3 TeV (left) and Ekin = 12 TeV (right), and 100 × 100 particles in each case. Green squares represent the starting points of particle trajectories intersecting the solar surface. Shadow-to-disk ratio amounts to 0.787 for the left case and 1.516 on the right. To guide the eye, blue circles indicate projected radii of 1 R (solid) and RSS = 2.5 R (dashed). We note that for the right plot, the lateral half-length d of the quadratic injection window is decreased from d = 1.5° to d = 0.8° for improved resolution.

thumbnail Fig. 5.

Total shadow area s as function of particle energy for dipolar field (8), normalized to the area of the geometrical disk. Two arrows indicate respective energy values of 3 and 12 TeV, for which shadow images are shown in Fig. 4. The respective side length 2 d of the injection window was chosen as 3.0° (1.6°, 1.0°) for E ∈ [1, 10] ([10, 20], [20, 100]) TeV to optimize spatial resolution, while at the same time keeping the field of view large enough to avoid losing any “shadow particles.”

This simple example may serve to illustrate that at least for large-scale magnetic fields, the shadow may not necessarily become less pronounced with increasing energy but may, in fact, get amplified considerably compared to the case of a non-magnetized shadow-casting obstacle. We come back to this phenomenon in Sect. 4.3.

3.3. Models of the primary cosmic-ray spectrum and composition

For the energy spectrum and composition of the simulations, we used the model from Gaisser (2012) with a mixed extragalactic component (hereinafter, HGm model) and from Gaisser & Honda (2002; hereinafter, GH model).

3.4. Passing probability and weighting scheme

In order to mimic the rotation of the solar magnetic field within a Carrington rotation, which has a synodic rotation period of ∼27.3 days, a total of 36 calculations were performed, revolving the injection plane (instead of the solar magnetic field) in steps of 10° around the rotation axis of the Sun. To quantify how likely it is for a particle with energy E and atomic number Z coming from direction i to pass the Sun without hitting its surface, we then defined the passing probability p of a specific Carrington rotation CR as the number npass of angular steps for which the particle passes the Sun, divided by the total number ntot = 36 of angular steps:

p i , Z , E , CR = n pass n tot · $$ \begin{aligned} p_{i, Z, E, \mathrm{CR} } = \frac{n_{\mathrm{pass} }}{n_{\mathrm{tot} }}\cdot \end{aligned} $$(10)

This was done for the nuclei of hydrogen, helium, nitrogen, aluminum, and iron with particle energies of 101.0, 101.6, and 102.2 TeV (corresponding to 10, ∼40 and ∼158 TeV). Particles with higher energies were omitted because of their lower particle flux. We ran our computations for the Carrington rotations approximately over the time window between November and February of each season since it is only between these months that the IceCube detector observes the cosmic-ray flux from a window around the Sun, considering the seasons 2007-08 through 2016-17.

The result for each set of atomic number and energy was weighted according to its abundance based on the HGm model (see Sect. 3.3). Assuming that the calculated Sun shadow does not vary considerably within an energy range of ±0.3 in log(E), we obtained the weighting factors gZ, E by integrating the energy spectra in an integration interval from 10x − 0.3 TeV to 10x + 0.3 TeV for all mentioned energies and atomic numbers. Using these weighting factors, the weighted arithmetic mean of all sets of energy and atomic number was determined for each Carrington rotation according to

p ¯ i , CR = Z , E p i , Z , E , CR · g Z , E Z , E g Z , E · $$ \begin{aligned} \bar{p}_{i, \mathrm{CR} } = \frac{\sum \limits _{Z, E}p_{i, Z, E, \mathrm{CR} }\cdot g_{Z, E}}{\sum \limits _{Z, E}g_{Z, E}}\cdot \end{aligned} $$(11)

In the last step, the mean of p ¯ i , CR $ \bar{p}_{i, \mathrm{CR}} $ for all Carrington rotations and the corresponding standard deviation were determined.

3.5. A quantitative measure for the shadow depth

In order to quantify the calculated shadowing effect of the Sun for a particular Carrington rotation, we compared the computed shadow size with that of an unmagnetized object with the Sun’s size and shape. Since each particle would hit the solar surface with an average probability of 1 p ¯ i , CR $ 1-\bar{p}_{i, \mathrm{CR}} $, we expect a total of N i p ¯ i , CR $ N-\sum\limits_{i} \bar{p}_{i, \mathrm{CR}} $ particles to contribute to the shadow. With RS ≡ arctan[R/(1 AU)] ≈ 0.27° as the angular radius of the Sun as seen from Earth, the angular area (solid angle) covered by the solar disk amounts to π R S 2 $ \pi R_{\mathrm{S}}^2 $. Since the injection area of size ω ≡ (4RS)2 contains N evenly distributed starting positions, (π R S 2 /ω)N=(π/16)N $ (\pi R_{\mathrm{S}}^2/\omega) N = (\pi/16) N $ shadowed events are expected from an unmagnetized object with the Sun’s size and shape. The normalized net shadowing effect for one Carrington rotation can therefore be quantified using the ratio

s CR = N i p ¯ i , CR ( π R S 2 / ω ) N = 16 π ( 1 1 N i p ¯ i , CR ) · $$ \begin{aligned} s_{\mathrm{CR} } = \frac{N-\sum \limits _{i} \bar{p}_{i, \mathrm{CR} }}{(\pi R_{\mathrm{S} }^2/\omega ) N} = \frac{16}{\pi } \bigg ( 1- \frac{1}{N} \sum \limits _{i} \bar{p}_{i, \mathrm{CR} } \bigg )\cdot \end{aligned} $$(12)

In particular, a ratio of sCR = 1 would mean that the shadow’s size is equivalent to a geometrical shadowing by an unmagnetized Sun. As the last step, the mean of those ratios

s = 1 n CR CR s CR $$ \begin{aligned} s = \frac{1}{n_{\mathrm{CR} }} \sum _{\mathrm{CR} } s_{\mathrm{CR} } \end{aligned} $$(13)

for the Carrington rotations considered was calculated. The uncertainty of s was calculated as the standard deviation of the ratios of the different Carrington rotations.

4. Results

The passing probabilities (compare Sect. 3.4) are visualized on a two-dimensional grid in longitude and latitude. The passing probability for each bin is represented by a color, with darker colors representing smaller passing probabilities and lighter colors representing larger ones. A passing probability of 1.0 in a specific bin means that every particle from that bin passed the Sun unimpeded, while a passing probability of 0.0 means that every particle hit the Sun. Two examples for resulting shadow plots can be seen in Fig. 6. For the purposes of readability, shadow plots related to the following sections are presented in Appendix A. For numerical studies of the shadow, the shadow size ratio s as described in Sect. 3.5 was used.

thumbnail Fig. 6.

Examples of visualizations of the calculated Sun shadow, where 2008-09 (top; year of least solar activity) and 2014-15 (bottom; year of most solar activity) seasons are shown, using the energy spectrum and composition according to the GH model.

4.1. Sun shadow as a function of mass number

In Fig. 7, the shadow for December 2015 can be seen at E = 40 TeV and for all different nuclei that were simulated. As expected, the deflection of particles increases with increasing atomic number (decreasing rigidity). In Fig. 8 the ratio s can be seen for pure proton, pure iron, and a mixed composition according to the HGm model for the years from 2007 through 2017. Some obvious features can be seen from the figures: First of all, as expected, the hydrogen shadow is much more focused and less fragmented as compared to the heavier components. Fragmentation and spread is increasing with increasing charge Z. The full shadow plots, displayed in the appendix in Figs. A.1 and A.2, also reveal that the measure for the shadow, s, is globally smaller for iron than for proton. In addition, the differences between the seasons are less pronounced for iron than for protons. The picture in the real world is expected to be proton-dominated up to PeV-energies, as the cosmic-ray spectrum itself is dominated by hydrogen (Wiebel-Sooth et al. 1998). A detailed investigation of the “true” Sun shadow is presented in the following subsections.

thumbnail Fig. 7.

Calculated shadow at E = 40 TeV for the different nuclei simulated.

thumbnail Fig. 8.

Calculated Sun shadow size ratio from 2007 through 2017 for pure proton, pure iron, and a mixed composition according to the HGm model.

4.2. Sun shadow as a function of energy

In Fig. 9, the shadow depth s, calculated as described above, is displayed for the three different energies. Each data point contains contributions by all five elements that were simulated with a composition according to the HGm model. As expected, the points for E = 10 TeV overall show the smallest s, indicating a weaker shadow, while the points for E = 158 TeV overall show the largest s, indicating a stronger shadow. In addition, the separation between the energies becomes stronger for seasons with high solar activity (season 2011-12 through season 2015-16). This is also expected as the effect of increasing magnetic activity and, in turn, a stronger magnetic field impacts low-energy particles the most. In Figs. A.3A.5, the full shadow plots can be seen for the different energies. In agreement with the numerical results, the shadow becomes more disk-like for higher energies, yielding s values close to unity, while it shows a distinct temporal variation at lower energies, resulting in a smearing of the shadow in seasons with high solar activity. In Figs. 10 and 11, the calculated s is plotted for all 15 different combinations of energy and mass number for the seasons 2008-09 and 2014-15, respectively. These seasons were chosen as respective examples for periods with notably low and high solar activity (compare with Fig. 1). For season 2014-15, the behavior of s as a function of energy and mass number is as expected: the shadow is stronger for high-energy particles with low atomic numbers and weaker for low-energy particles with high atomic numbers. For season 2008-09, however, the picture is different. While for iron nuclei, s increases with increasing energy, there is a peak structure for aluminum and nitrogen nuclei, and s even decreases for helium and hydrogen nuclei. In order to better understand this effect, we reduce the particle properties to only one parameter in the following section.

thumbnail Fig. 9.

Calculated Sun shadow in the years from 2007 through 2017 for the different energies with a mixed composition according to the HGm model.

thumbnail Fig. 10.

Calculated Sun shadow in season 2008-09 for the different energies and mass numbers. Solar activity was comparatively low during this season (compare with Fig. 1).

thumbnail Fig. 11.

Calculated Sun shadow in season 2014-15 for the different energies and mass numbers. Solar activity was comparatively high during this season (compare with Fig. 1).

4.3. Sun shadow as a function of rigidity

In this section, we use rigidity R ∼ E/Z in order to reduce the dependency of the shadow on energy and atomic number to one parameter. In Figs. 12 and 13, the calculated s is plotted as a function of rigidity for the respective seasons 2008-09 and 2014-15. The expected equivalence of particles with the same rigidity is confirmed for both seasons. As an example, the middle blue square in Fig. 13, representing ∼101.6 TeV nitrogen nuclei, shows an s comparable to that of ∼102.2 TeV iron nuclei, the rightmost orange pentagon. With rigidities of ∼5700 GV and ∼6100 GV respectively, these particles are expected to show a similar behavior with respect to their gyro-radii. For season 2008-09, however, s does not increase monotonically, as it does in season 2014-15, but, rather, it peaks at a certain rigidity between 103.5 and 104.0 GV. A possible explanation for this is that s seems to be significantly larger than 1.0 for years with low solar activity (compare Sect. 4.1), that is, in years when the magnetic field is approximately described by a simple dipole. It seems possible that particles above a certain rigidity no longer form the shadow typical for a dipole field, but rather “converge” towards the shadow of an unmagnetized disk, which, by definition, would result in values of s ≈ 1. Particularly in view of the non-monotonous dependence of shadow size on rigidity that was found for the dipolar test case (see Sect. 3.2) and illustrated in Fig. 5, the transition between these two regions – s for a dipole field versus s for an unmagnetized disk – can possibly explain the peak visible in Fig. 10 as well.

thumbnail Fig. 12.

Calculated Sun shadow in season 2008-09 as a function of rigidity.

thumbnail Fig. 13.

Calculated Sun shadow in season 2014-15 as a function of rigidity.

Thus, our simulations indicate that an energy-dependent study of the Sun shadow by cosmic-ray observatories would be highly beneficial. From these studies, we expect the observatories to see the dipole-like energy dependence, with a monotonic rise up to a peak energy. In the peak, we expect the shadow effect to be up to 40% stronger than the geometrical-disk shadow. At energies above this peak, we expect the shadow to converge to the size of the geometrical one. In years of high solar activity, we predict a different behavior, which is a purely monotonic increase of the shadow depth, starting at low energies with a shadow much shallower than the geometric one. For all energies, the shadow is expected to be shallower than the geometrical one. Toward high energies, it converges to the geometrical shadow depth again.

4.4. Sun shadow for different cosmic-ray flux models

In Fig. 14, the calculated s is plotted for the HGm and GH models for the seasons 2007-08 through 2016-17. While there is a systematic offset in the total magnitude of s, its temporal variation is very similar for both models, suggesting that the shape of the temporal variation does not depend significantly on the primary cosmic-ray model used. In Figs. A.6 and A.7 (see appendix), the calculated passing probabilities and full shadow plots for the HGm and GH models for seasons 2007-08 through 2016-17 can be seen.

thumbnail Fig. 14.

Comparison of s for calculated Sun shadow from 2007 through 2017 among the HGm and GH models.

4.5. Sun shadow compared to solar activity

Figure 15 shows the temporal behavior of s (assuming the HGm model) and the sunspot number averaged over the time periods considered in this work (data taken from SILSO World Data Center 2007–2017). We note that the time range of the Carrington rotations considered for a specific year and the time range in which the sunspot numbers are measured are slightly different since the duration of a Carrington rotation is always shorter than the duration of a calendar month. Figure 16 shows the correlation of s with the sunspot number, together with a linear fit of the form

s = a + b · SSN . $$ \begin{aligned} s = a + b \cdot \mathrm{SSN} . \end{aligned} $$(14)

thumbnail Fig. 15.

Temporal behavior of shadow ratio s (assuming the HGm model) and sunspot number from 2007 through 2017. Sunspot numbers are taken from SILSO World Data Center (2007–2017).

thumbnail Fig. 16.

Correlation of ratio s and sunspot number from 2007 through 2017. Sunspot numbers are taken from SILSO World Data Center (2007–2017).

While a is reasonable close to 1.0, which is expected as per the discussion in Sects. 4.2 and 4.3, there is currently no straightforward interpretation of b. For the Spearman’s rank correlation coefficient ρ of s and the sunspot number, we find ρSpearman = −0.86.

5. Summary

In this paper, we investigate cosmic-ray transport across the Sun’s corona in the energy range between ∼5 TeV and ∼316 TeV and for hydrogen, helium, nitrogen, aluminum, and iron nuclei. For that purpose, particles with three different energies were studied, namely 101.0, 101.6, and 102.2 TeV (corresponding to 10, ∼40, and ∼158 TeV). We chose these supporting points based on their being equidistant in log ECR. The energy range was selected to reflect the range that most detectors measuring the cosmic-ray Sun shadow (such as IceCube, Tibet, ARGO-YBJ, HAWC, and others) are sensitive to. The data sample used by IceCube has a median energy of ⟨ECR ⟩ ∼40 TeV, with 68% of the events between 11 and 200 TeV. The mass numbers of the nuclei are chosen to reflect the most abundant mass groups, with nitrogen representing the CNO and aluminum representing the MgAlSi element group. By choosing these mass numbers and energies, we are able to compute a generic prediction for the expected time variability of the cosmic-ray Sun shadow that can be compared to measurements by the experiments mentioned above. We associate the above-mentioned time variability with the time-changing solar magnetic field and its influence on the shape and depth of the Sun shadow of cosmic rays studied by ground-based cosmic-ray experiments such as IceCube, Tibet, ARGO-YBJ, and HAWC. In the future, these observatories will be able to further improve their spatial and temporal resolution, possibly contributing to investigations of the solar cycle via indirect measurements of magnetic field properties.

For the purpose mentioned above, we performed a numerical integration of particle trajectories in the solar magnetic field of the past ten years. Magnetogram data were obtained from GONG (Harvey et al. 1996) and interpreted in the PFSS model using the FDIPS code (Toth et al. 2016). For each year from 2007 through 2017, the respective four months, January, February, November, and December, were simulated. We chose to simulate these months as they correspond to the time period per year in which the Sun is visible from the South Pole. This way, our results are easily applicable to that of the IceCube detector, corresponding to the relevant energy range and observation times. It should be noted, however, that no detector effects, such as a limited acceptance or a limited angular resolution, have been taken into account.

The results of our simulations are provided in a two-dimensional representation. The variation of the cosmic-ray Sun shadow is verified in this simulation, and as expected, the effect is strongest at the lowest simulated energy of 10 TeV. The variation of the shadow furthermore shows that in times of low solar activity, the shadow is relatively compact near the Sun’s center. In times of high solar activity, however, the shadow tends to be washed out. It needs to be investigated by IceCube and other cosmic-ray experiments whether it is possible to detect such features of systematic deviation from a dipole structure despite the limitations associated with angular resolution. Such a measurement would make it possible to indirectly study the structure of the solar magnetic field and, thus, to contribute to its understanding. In this sense, our findings of the different energy behavior of the shadow for low activity (dipole-like) and high activity (monotonous) are very interesting as they would pose another possibility for the detection of a dipole-like shadow behavior. We can make the clear prediction that in times of low solar activity, the energy-behavior should resemble that of a pure dipole, with a shadow that becomes stronger with energy up to a maximum value, in which it exceeds the depth of the geometrical shadow. For the years studied here, the shadow is expected to become 40% deeper than the geometrical shadow. Above this peak energy, the shadow is expected to converge to the geometrical size. This behavior stands in contrast to years of high activity. Here we find a monotonous increase of the shadow depth with energy, converging toward the geometrical size, never exceeding the geometrical shadow. Future studies will show if this is an effect of the unordered field of the Sun, which is much more significant in years of high solar activity as compared to the low activity years of the Sun.

We can further show that while the shadow is still significantly distorted at 40 TeV, the deviation from an undisturbed shadow becomes extremely small at 158 TeV (compare with Sect. 4.2). Thus, we expect effects of the inner heliospheric field on cosmic-ray propagation to become negligible with regard to the variation of the Sun shadow at these high energies. This result is in agreement with recent HAWC measurements (Nisa et al. 2017), which show a symmetric, deep shadow at energies around 100 TeV. Future investigations of several years of data with HAWC will help to further quantify these findings.

The absolute normalization of the shadow depth in our simulation is not directly comparable to the IceCube data as we did not take into account detection effects like a limited acceptance or a limited angular resolution. In comparing the relative change in the deficit between the different years, however, our simulations shows a similar temporal variation of the shadow depth as the data taken with IceCube (Aartsen et al. 2019).

6. Outlook

We consider this paper as a first step toward achieving a detailed understanding of the effects of the solar magnetic field variation on the Sun’s cosmic-ray shadow as measured by different experiments on Earth. In the future, there will be opportunities for further work, which either need improved data from the experiments or require modeling approaches that are more computationally involved.

Magnetic field modeling. One aspect concerns modeling the inner heliospheric magnetic field, which is done on the basis of available data. The data are collected by spacecraft, such as Helios (Schwenn & Marsch 2012), which do not simultaneously provide data from the far side. We used the PFSS ansatz (Schatten et al. 1969; Altschuler & Newkirk 1969) in order to model the global coronal magnetic field, but it is desirable to compare these results to other models. The Tibet group, for instance, favors the CSSS model (Amenomori et al. 2013) because in using this model, they found better agreement with the measured data. A similar comparison between these two models and IceCube data would be an asset for a better understanding of the solar coronal magnetic field. In addition to the existing studies by Tibet, a comparison of different coronal magnetic field models on the analysis level using IceCube data would be very useful for evaluating the validity of different models. However, as these models both assume a purely radial magnetic field at a defined boundary sphere, no fundamental differences between these two are expected. Ultimately, it will be desirable to probe more complex models like force-free field approximations (see, e.g., Tadesse et al. 2014) or even full magnetohydrodynamic (MHD) simulations of the solar magnetic field.

Quantitative comparison with data. In this paper, we compare our results regarding the time variability of the Sun shadow with current IceCube data, which is available for five years (see Aartsen et al. 2019). In the future, IceCube will extend its data analysis to more years of data for the same solar cycle. Also, HAWC will be able to provide information in a similar energy range. As a result, a comprehensive comparison of data and simulation will be possible: first, by deducing the effects on the shadow caused by the magnetic field from simulations only; and, second, by comparing data and simulation on the detector level. Future arrays like IceCube-Gen2 (IceCube-Gen2 Collaboration 2014) will improve such studies even further, possibly enabling an investigation of the properties of the magnetic field through the combination of simulation and detection of the cosmic-ray shadow with better angular and time resolution.

Integrating the particle propagation into a detector environment is necessary in order to test which magnetic field configuration fits the data best. This needs to be done for each instrument individually and it can only be done through the collaborations themselves. Studies as those we present here, without including detector effects, are still desirable as this is the best way to find effects that arise from the magnetic field. The influence of coronal mass ejections (CMEs) on the cosmic-ray Sun shadow seen by different experiments is also something that can be studied, as shown in Amenomori et al. (2018).

Toward a complete model of transport and interaction. Furthermore, adding particle interactions in the solar corona and in the Earth’s atmosphere to the simulations will be another step that would increase our understanding of the phenomenon of cosmic-ray shadowing.

Acknowledgments

We would like to thank Tobias Wiengarten for assistance with the magnetic field data, as well as Malina Reytemeier, Paul Evenson and the entire IceCube collaboration, as well as members of the Ruhr Astroparticle Plasma Physics (RAPP) Center for helpful discussions on the subject. We are further grateful for the support from the BMBF (05A14PC1) and from MERCUR (St-2014-0040). This work utilizes data from the National Solar Observatory Integrated Synoptic Program, which is operated by the Association of Universities for Research in Astronomy, under a cooperative agreement with the National Science Foundation and with additional financial support from the National Oceanic and Atmospheric Administration, the National Aeronautics and Space Administration, and the United States Air Force. The GONG network of instruments is hosted by the Big Bear Solar Observatory, High Altitude Observatory, Learmonth Solar Observatory, Udaipur Solar Observatory, Instituto de Astrofísica de Canarias, and Cerro Tololo Interamerican Observatory.

References

  1. Aartsen, M. G., Ackermann, M., Adams, J., et al. 2019, ApJ, 872, 133 [NASA ADS] [CrossRef] [Google Scholar]
  2. Abdo, A. A., Ackermann, M., Ajello, M., et al. 2011, ApJ, 734, 116 [NASA ADS] [CrossRef] [Google Scholar]
  3. Altschuler, M. D., & Newkirk, G. 1969, Sol. Phys., 9, 131 [Google Scholar]
  4. Alves Batista, R., Erdmann, M., Evoli, C., et al. 2013, ArXiv e-prints [arXiv:1307.2643] [Google Scholar]
  5. Amenomori, M., Bi, X. J., Chen, D., et al. 2013, Phys. Rev. Lett., 111, 011101 [NASA ADS] [CrossRef] [Google Scholar]
  6. Amenomori, M., Bi, X. J., Chen, D., et al. 2018, ApJ, 860, 13 [NASA ADS] [CrossRef] [Google Scholar]
  7. Argüelles, C. A., de Wasseige, G., Fedynitch, A., & Jones, B. J. P. 2017, JCAP, 7, 024 [NASA ADS] [CrossRef] [Google Scholar]
  8. Birdsall, C. K., & Langdon, A. B. 1985, Plasma Physics via Computer Simulation (McGraw-Hill Book Company), 13 [Google Scholar]
  9. Boris, J. 1970, Proceedings of the 4th Conference on Numerical Simulation of Plasmas (NRL), 3 [Google Scholar]
  10. Bos, F. 2017, PhD Thesis, Ruhr-Universität Bochum, Germany [Google Scholar]
  11. Clark, G. W. 1957, Phys. Rev., 108, 450 [NASA ADS] [CrossRef] [Google Scholar]
  12. Edsjö, J., Elevant, J., Enberg, R., & Niblaeus, C. 2017, JCAP, 6, 033 [NASA ADS] [CrossRef] [Google Scholar]
  13. Enriquez-Rivera, O., & Lara, A. 2016, Proceedings, 34th International Cosmic Ray Conference (ICRC2015), The Hague, The Netherlands, 099 [Google Scholar]
  14. Evoli, C., Gaggero, D., Grasso, D., & Maccione, L. 2008, JCAP, 2008, 018 [NASA ADS] [CrossRef] [Google Scholar]
  15. Gaisser, T. K. 2012, Astropart. Phys., 35, 801 [NASA ADS] [CrossRef] [Google Scholar]
  16. Gaisser, T. K., & Honda, M. 2002, Annu. Rev. Nucl. Part. Sci., 52, 153 [NASA ADS] [CrossRef] [Google Scholar]
  17. Hapgood, M. 1992, Planet. Space Sci., 40, 711 [Google Scholar]
  18. Harvey, J. W., Hill, F., Hubbard, R. P., et al. 1996, Science, 272, 1284 [NASA ADS] [CrossRef] [PubMed] [Google Scholar]
  19. Hoyt, D. V., & Schatten, K. H. 1998, Sol. Phys., 179, 189 [NASA ADS] [CrossRef] [Google Scholar]
  20. IceCube-Gen2 Collaboration (Aartsen, M. G., et al.) 2014, ArXiv e-prints [arXiv:1412.5106] [Google Scholar]
  21. Kissmann, R. 2014, Astropart. Phys., 55, 37 [NASA ADS] [CrossRef] [Google Scholar]
  22. Merten, L., Becker Tjus, J., Fichtner, H., Eichmann, B., & Sigl, G. 2017, JCAP, 6, 046 [Google Scholar]
  23. Ng, K. C. Y., Beacom, J. F., Peter, A. H. G., & Rott, C. 2016, Phys. Rev. D, 94, 023004 [NASA ADS] [CrossRef] [Google Scholar]
  24. Ng, K. C. Y., Beacom, J. F., Peter, A. H. G., & Rott, C. 2017, Phys. Rev. D, 96, 103006 [NASA ADS] [CrossRef] [Google Scholar]
  25. Nisa, M. U., Hampel, Z., & HAWC Collaboration 2017, APS Meeting Abstracts, Y3.006 [Google Scholar]
  26. Nisa, M. U., Beacom, J. F., BenZvi, S. Y., et al. 2019, BAAS, 51, 194 [Google Scholar]
  27. Orlando, E., & Strong, A. W. 2008, A&A, 480, 847 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  28. Parker, E. N. 1958, ApJ, 128, 664 [Google Scholar]
  29. Priest, E. 1982, Solar Magnetohydrodynamics (Netherlands: Springer) [CrossRef] [Google Scholar]
  30. Qin, H., Zhang, S., Xiao, J., et al. 2013, Phys. Plasmas, 20, 084503 [NASA ADS] [CrossRef] [Google Scholar]
  31. Schatten, K. H., Wilcox, J. M., & Ness, N. F. 1969, Sol. Phys., 6, 442 [Google Scholar]
  32. Schlickeiser, R. 2002, Cosmic Ray Astrophysics (Heidelberg: Springer) [CrossRef] [Google Scholar]
  33. Schwenn, R., & Marsch, E. 2012, in Physics of the Inner Heliosphere I: Large-Scale Phenomena (Berlin, Heidelberg: Springer), Phys. Chem. Space [Google Scholar]
  34. Seckel, D., Stanev, T., & Gaisser, T. K. 1991, ApJ, 382, 652 [NASA ADS] [CrossRef] [Google Scholar]
  35. SILSO World Data Center 2007–2017, International Sunspot Number Monthly Bulletin and Online Catalogue [Google Scholar]
  36. Solanki, S. K., Inhester, B., & Schüssler, M. 2006, Rep. Prog. Phys., 69, 563 [NASA ADS] [CrossRef] [Google Scholar]
  37. Strong, A. W., & Moskalenko, I. V. 1998, ApJ, 509, 212 [NASA ADS] [CrossRef] [Google Scholar]
  38. Strong, A. W., Porter, T. A., Digel, S. W., et al. 2010, ApJ, 722, L58 [NASA ADS] [CrossRef] [Google Scholar]
  39. Tadesse, T., Wiegelmann, T., MacNeice, P. J., et al. 2014, Sol. Phys., 289, 831 [NASA ADS] [CrossRef] [Google Scholar]
  40. Toth, G., van der Holst, B., & Huang, Z. 2016, Astrophysics Source Code Library [record ascl:1606.011] [Google Scholar]
  41. Vladimirov, A. E., Digel, S. W., Jóhannesson, G., et al. 2011, Comput. Phys. Commun., 182, 1156 [NASA ADS] [CrossRef] [Google Scholar]
  42. Wiebel-Sooth, B., Biermann, P. L., & Meyer, H. 1998, A&A, 330, 389 [NASA ADS] [Google Scholar]
  43. Wiengarten, T., Kleimann, J., Fichtner, H., et al. 2014, ApJ, 788, 80 [Google Scholar]
  44. Zhao, X., & Hoeksema, J. T. 1995, J. Geophys. Res.: Space Phys., 100, 19 [NASA ADS] [CrossRef] [Google Scholar]

Appendix A: Shadow plots

A.1. Sun shadow as a function of mass number

In Figs. A.1 and A.2, the Sun shadow can be seen for proton and iron nuclei only.

thumbnail Fig. A.1.

Calculated Sun shadow from 2007 through 2017 for protons only.

thumbnail Fig. A.2.

Calculated Sun shadow from 2007 through 2017 for iron nuclei only. We call attention to the narrow color scale.

A.2. Sun shadow as a function of energy

In Figs. A.3A.5, the Sun shadow can be seen at 10, 40 and 158 TeV, respectively.

thumbnail Fig. A.3.

Calculated Sun shadow from 2007 through 2017 for 10 TeV. Each plot contains contributions by all five elements that were simulated with a composition according to the HGm model.

thumbnail Fig. A.4.

Calculated Sun shadow from 2007 through 2017 for 40 TeV. Each plot contains contributions by all five elements that were simulated with a composition according to the HGm model.

thumbnail Fig. A.5.

Calculated Sun shadow from 2007 through 2017 for 158 TeV. Each plot contains contributions by all five elements that were simulated with a composition according to the HGm model.

A.3. Sun shadow for different cosmic-ray flux models

In Figs. A.6 and A.7, the Sun shadow can be seen for the HGm and GH model, respectively.

thumbnail Fig. A.6.

Calculated Sun shadow from 2007 through 2017 using the energy spectrum and composition according to the HGm model.

thumbnail Fig. A.7.

Calculated Sun shadow from 2007 through 2017 using the energy spectrum and composition according to the GH model.

All Figures

thumbnail Fig. 1.

Sunspot numbers taken from SILSO World Data Center (2007–2017) and averaged over the four months of each season are shown for a ten-year time period. The seasons 2007-08 and 2008-09 report a low sunspot number, while in seasons 2011-12, 2012-13, 2013-14, and 2014-15 a higher number is detected, before decreasing again afterwards. This approximately describes the solar activity cycle of 11 years (see, e.g., Solanki et al. 2006).

In the text
thumbnail Fig. 2.

Magnetic field strength on the solar surface and selected magnetic field lines for Carrington rotations 2077 (December 2008 with low solar activity, left) and 2158 (December 2014 close to the solar activity maximum, right), based on GONG magnetogram data and the PFSS model. Field lines are shown between the solar surface and the source surface at 2.5 solar radii, with blue indicating outgoing (Br >  0), brown indicating incoming (Br <  0) polarity. Magnetic field strength on the solar surface is shown using a rainbow color scale. It can clearly be seen that the dipolar structure which typically prevails during solar minimum turns into a more irregular field structure towards solar maximum.

In the text
thumbnail Fig. 3.

Schematic view of the simulation setup (distances not to scale). Starting points of anti-particles as indicated by the blue bubbles are arranged in the injection window. Timing information for the trajectories of the anti-particles is given in color code. The starting points of those anti-particles hitting the Sun are removed from the plot. In this way, the cosmic-ray deficit caused by the Sun is simulated using back-tracing.

In the text
thumbnail Fig. 4.

Shadow maps for mono-energetic protons in the magnetic field of a dipole, using the example of Ekin = 3 TeV (left) and Ekin = 12 TeV (right), and 100 × 100 particles in each case. Green squares represent the starting points of particle trajectories intersecting the solar surface. Shadow-to-disk ratio amounts to 0.787 for the left case and 1.516 on the right. To guide the eye, blue circles indicate projected radii of 1 R (solid) and RSS = 2.5 R (dashed). We note that for the right plot, the lateral half-length d of the quadratic injection window is decreased from d = 1.5° to d = 0.8° for improved resolution.

In the text
thumbnail Fig. 5.

Total shadow area s as function of particle energy for dipolar field (8), normalized to the area of the geometrical disk. Two arrows indicate respective energy values of 3 and 12 TeV, for which shadow images are shown in Fig. 4. The respective side length 2 d of the injection window was chosen as 3.0° (1.6°, 1.0°) for E ∈ [1, 10] ([10, 20], [20, 100]) TeV to optimize spatial resolution, while at the same time keeping the field of view large enough to avoid losing any “shadow particles.”

In the text
thumbnail Fig. 6.

Examples of visualizations of the calculated Sun shadow, where 2008-09 (top; year of least solar activity) and 2014-15 (bottom; year of most solar activity) seasons are shown, using the energy spectrum and composition according to the GH model.

In the text
thumbnail Fig. 7.

Calculated shadow at E = 40 TeV for the different nuclei simulated.

In the text
thumbnail Fig. 8.

Calculated Sun shadow size ratio from 2007 through 2017 for pure proton, pure iron, and a mixed composition according to the HGm model.

In the text
thumbnail Fig. 9.

Calculated Sun shadow in the years from 2007 through 2017 for the different energies with a mixed composition according to the HGm model.

In the text
thumbnail Fig. 10.

Calculated Sun shadow in season 2008-09 for the different energies and mass numbers. Solar activity was comparatively low during this season (compare with Fig. 1).

In the text
thumbnail Fig. 11.

Calculated Sun shadow in season 2014-15 for the different energies and mass numbers. Solar activity was comparatively high during this season (compare with Fig. 1).

In the text
thumbnail Fig. 12.

Calculated Sun shadow in season 2008-09 as a function of rigidity.

In the text
thumbnail Fig. 13.

Calculated Sun shadow in season 2014-15 as a function of rigidity.

In the text
thumbnail Fig. 14.

Comparison of s for calculated Sun shadow from 2007 through 2017 among the HGm and GH models.

In the text
thumbnail Fig. 15.

Temporal behavior of shadow ratio s (assuming the HGm model) and sunspot number from 2007 through 2017. Sunspot numbers are taken from SILSO World Data Center (2007–2017).

In the text
thumbnail Fig. 16.

Correlation of ratio s and sunspot number from 2007 through 2017. Sunspot numbers are taken from SILSO World Data Center (2007–2017).

In the text
thumbnail Fig. A.1.

Calculated Sun shadow from 2007 through 2017 for protons only.

In the text
thumbnail Fig. A.2.

Calculated Sun shadow from 2007 through 2017 for iron nuclei only. We call attention to the narrow color scale.

In the text
thumbnail Fig. A.3.

Calculated Sun shadow from 2007 through 2017 for 10 TeV. Each plot contains contributions by all five elements that were simulated with a composition according to the HGm model.

In the text
thumbnail Fig. A.4.

Calculated Sun shadow from 2007 through 2017 for 40 TeV. Each plot contains contributions by all five elements that were simulated with a composition according to the HGm model.

In the text
thumbnail Fig. A.5.

Calculated Sun shadow from 2007 through 2017 for 158 TeV. Each plot contains contributions by all five elements that were simulated with a composition according to the HGm model.

In the text
thumbnail Fig. A.6.

Calculated Sun shadow from 2007 through 2017 using the energy spectrum and composition according to the HGm model.

In the text
thumbnail Fig. A.7.

Calculated Sun shadow from 2007 through 2017 using the energy spectrum and composition according to the GH model.

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.