Simulation of the isotropic ultra-high energy photons flux in the solar magnetic field

Both the lack of observation of ultra-high energy (UHE) photons and the limitations of the state-of-the-art methodology being applied for their identification motivate studies on alternative approaches to the relevant simulations and the related observational strategies. One of such new approaches is proposed in this report and it concerns new observables allowing indirect identification of UHE photons through cosmic ray phenomena composed of many spatially correlated extensive air showers or primary cosmic rays observed at one time. The study is based on simulations of interactions of UHE photons with the magnetic field of the Sun using the PRESHOWER program with some essential modifications. One of the expected results of such interactions is a generation of cosmic ray ensembles (CRE) in the form of very thin and very elongated cascades of secondary photons of energies spanning the whole cosmic ray energy spectrum. Upon entering the Earth's atmosphere, these cascades or their parts may generate uniquely characteristic walls of spatially correlated extensive air showers, and the effect is expected also in cases when primary UHE photons are not directed towards the Earth. Particle distributions in these multi-primary UHE photon footprints are expected to have thicknesses of the order of meters and elongations reaching even hundreds of millions kilometers, making them potentially observable with a global, multi-experiment approach, including re-exploring of the historical data, with the expected event rate exceeding the capabilities of even very large cosmic ray observatories. The methods described in this report allow for simulating potentially observable quantities related to CRE induced by UHE photons: densities, energy spectra and geographical orientations of particles at the top of the Earth's atmosphere.

Abstract: Both the lack of observation of ultra-high energy (UHE) photons and the limitations of the stateof-the-art methodology being applied for their identification motivate studies on alternative approaches to the relevant simulations and the related observational strategies. One of such new approaches is proposed in this report and it concerns new observables allowing indirect identification of UHE photons through cosmic ray phenomena composed of many spatially correlated extensive air showers or primary cosmic rays observed at one time. The study is based on simulations of interactions of UHE photons with the magnetic field of the Sun using the PRESHOWER program with some essential modifications. One of the expected results of such interactions is a generation of cosmic ray ensembles (CRE) in the form of very thin and very elongated cascades of secondary photons of energies spanning the whole cosmic ray energy spectrum. Upon entering the Earth's atmosphere, these cascades or their parts may generate uniquely characteristic walls of spatially correlated extensive air showers, and the effect is expected also in cases when primary UHE photons are not directed towards the Earth. Particle distributions in these multi-primary UHE photon footprints are expected to have thicknesses of the order of meters and elongations reaching even hundreds of millions kilometers, making them potentially observable with a global, multi-experiment approach, including re-exploring of the historical data, with the expected event rate exceeding the capabilities of even very large cosmic ray observatories. The methods described in this report allow for simulating potentially observable quantities related to CRE induced by UHE photons: densities, energy spectra and geographical orientations of particles at the top of the Earth's atmosphere. A measurement of at least one of these quantities would be equivalent to a confirmation of the existence of UHE photons which would give an insight to fundamental physics processes at unprecedentedly high energies, far beyond the reach of man-made accelerators. On the other hand, a lack of such an observation would allow for further constraining of these fundamental processes with the physically new upper limits on UHE photon fluxes. The novel advantage of such non-observation results would lay in the purely electrodynamical character of the underlying simulations which are fully independent on extrapolations of hadronic interaction models by many orders of magnitude. Such extrapolations are necessary in the UHE photon identification methods based on the analyses of properties of individual extensive air showers presently used to determine the UHE photon upper limits.

Introduction
The nature of the cosmic rays of the highest energies is one of the unsolved problems in physics. Their sources and mechanisms of acceleration remain unknown for us, as they are the subject of ongoing research (see e.g. Ref. [1] and the references therein). Therefore, the detailed analysis of physical phenomena at the highest energy regime may improve our understanding of the universe and enable us to probe modern physical theories. As far as we understand the interactions of subatomic particles described by the Standard Model, there should be the cut-off of the cosmic rays energy spectrum, due to the interactions of cosmic radiation components with cosmic microwave background radiation. At cosmic rays energies about 5-6×10 19 eV the collisions of protons in the cosmic ray flux with cosmic microwave background photons may cause the production of the short-lived particle ∆ + , which decays into a nucleon with the energy lower than the initial proton energy. This phenomenon is called the Greisen-Zatsepin-Kuzmin (GZK) limit (see [2] and [3]). In turn, heavier nuclei of cosmic radiation at energies about 10 20 eV may undergo photodisintegration and in this way also lower their energy. Therefore, the maximum distance from the source for the possible observation of the highest energy cosmic rays should be limited. Recent experimental data [4] indicate that the cosmic rays energy spectrum is suppressed, but it is impossible to determine whether the observations are consistent with the mentioned theory. On the other hand, the nature of detected cosmic radiation particles with energies exceeding 10 20 eV is unknown, as there are no potential sources of such particles in the Earth vicinity. Thus, the investigation of ultra-high energy cosmic rays is significant.
The problem of the highest energy cosmic rays gave rise to new theoretical models and predictions (see e.g. Ref. [1] and the references therein), that can be justified or falsified only by detailed analysis of processes at the ultra-high energy regime. One of these predictions is the existence of the Super Heavy Dark Matter, which decay or annihilation may lead to the production of ultra-high energy cosmic ray flux composed mainly of photons (see Ref. [5], [6], [7] and [8]). Hence, the study of ultra-high energy photons can be considered as an indirect search for heavy dark matter particles. Similarly, the rest of the new physical theories predict ultra-high energy (UHE) photon creation. However, the expected UHE photons flux is much higher than the observed one. A possible explanation for lowering the UHE photons flux is the cosmic ray ensembles (CRE) phenomenon. Providing that a UHE photon is cascading already in space, then a group of secondary cosmic rays might reach the top of the atmosphere instead of an individual particle. Thus, instead of only one extensive air shower (EAS) created in the atmosphere, a group of (EAS) is generated. This correlation may involve arrival time relations and/or spatial distributions to provide a characteristic signature, although a detection of this phenomenon might require an extended detector array, even about the size of the Earth, as proposed by the Cosmic-Ray Extremely Distributed Observatory (CREDO) [1]. The project aims at analysing data from professional observatories, educational cosmic ray facilities, and also from the smartphones with the CREDO Detector [9,10] mobile application, where a CMOS camera is used as a particle detector. For proper collected data interpretation dedicated to identification of CRE, the process of UHE photon cascading, formation of correlated air showers, and detector response need to be correctly modelled. This is the reason why, among other researches, the CREDO project probes the scenario of UHE photon cascading in the solar magnetic field [11]. The specific methods dedicated to utilizing this phenomenon for planning novel indirect UHE photon searches are presented in this article.
The previous studies of the CRE effect [1,11] showed that it is possible to observe some photons from cascades generated by the primary UHE photons also in cases where they are not aimed at the Earth. Here, we use the PRESHOWER program [12,13] with modifications to include such cases in the simulations for the first time, and deliver complete distributions of CRE photons at the top of the atmosphere which point to relevant observational strategies based on the properties of these distributions.

Model of ultra-high energy photon interaction
We modelled the Sun's vicinity by a surrounding sphere with radius 6 R . It was also assumed that UHE photons flux reaching the start sphere is isotropic. However, in the case of the UHE photon direction, only photons that come from outside of the sphere are considered to avoid the case duplication. Limits on the UHE photons are applied in accordance with [14]. Then, the cascading of UHE photon in the solar magnetic field is considered. Physical processes that are responsible for cascading the UHE photon in space are magnetic pair production, the bending of the charged particle trajectory, and synchrotron photon emission. In the considered case, a primary UHE photon may convert to a positron and electron in a sufficiently strong magnetic field. Then, trajectories of the created pair are bending, due to the presence of a magnetic field. Hence, the positron and the electron emit synchrotron photons, that may be observable as photons cascade on the top of the atmosphere. The following mathematical description of the UHE photon cascading is based on Ref. [12] and the references therein. The mathematical formulations of quantum mechanics of magnetic pair production used in this work were taken from [15] and [16]. On the basis of [15], the number of pairs created by n photons photons in transverse magnetic field B is described by the following expression where dl is photon path length and α(χ) is the attenuation coefficient, that depends on parameter where hν -photon energy, m e -electron mass and B cr = 4.414 × 10 13 G is the critical magnetic field strength.The attenuation coefficient depends on the probability transition for magnetic pair production. Using the first-order perturbation method and assuming that B B cr , the approximated formula for α(χ) may be expressed as where α QED is the fine-structure constant,λ c is reduced Compton wavelength and auxiliary function T(χ) is approximated by . Whereas, the K 1/3 2 3χ in the expression (2) is modified Bessel function. According to equation (1), the probability of conversion within infinitesimal distance dl is and for larger distance L it should be calculated as The conversion rate is negligible unless the condition χ = 1 2 hν m e c 2 B B cr ≈ 0.1 is satisfied. Thus, magnetic pair production has not been observed so far. The parameter χ also determines the shape of the pair-member energy distribution, which is given by the expression where is a fractional pair-member energy. The derivation of this formula is presented in detail in [16]. The radial acceleration of electron and positron in the solar magnetic field is calculated using wherev is a velocity versor, q is particle charge, E is particle kinetic energy and B is an external magnetic field. For short time interval ∆t, altered particle directionv(t + ∆t) may be approximated with Taylor series expansion and in this casev In the PRESHOWER simulator an appropriate time interval ∆t is chosen according to the particle energy and the magnetic field variability. During the first half of the ∆t period, the particle propagates in the previous direction. Then the particle direction is recalculated in accordance with (5) and in the second half of the ∆t particle is moving in an altered direction. In the end, the particle direction is recalculated once again. The energy distribution of emitted synchrotron radiation for electron at ultra-relativistic regime, according to [17], is expressed by where E is electron/positron energy, and where hν is emitted photon energy. In turn, energy emitted by electron or positron within path length dl is equal to where P is radiation power and r 0 is classical electron radius. Therefore, one can obtain that the probability of emission of synchrotron photon within infinitesimal path length dl is where and dN -the number of photons with energy [hν, hν + d(hν)], emitted within the path length dl. According to [18], at ultra-relativistic regime, synchrotron photons are emitted in half opening angle about charge particle direction, where γ is the Lorentz factor. The azimuthal angle of emitted photon direction is random, consistent with uniform distribution. Therefore, at the ultrarelativistic regime the direction of the emitted synchrotron photon is nearly tangential to the electron or positron trajectory. In this way, synchrotron photons may reproduce the projection of a pair-member trajectory onto the atmosphere plane. For that reason the spatial distribution of secondary photons on the top of the atmosphere is expected to be a very thin curve. Due to natural variations and complexity of the solar magnetic field, it is impossible to apply a satisfying, well approximated model in the PRESHOWER simulator. Thus, the solar magnetic field is modelled by the simple dipole approximation. The magnetic dipole moment used in simulations equals 6.87 × 10 32 G cm 3 .

Materials and Methods
In order to model the UHE photon cascading in the solar magnetic field, the PRESHOWER program is used in this work. In this program, interactions are modelled by Monte Carlo simulations. However, the program had to be modified to enable the time cumulative distribution of secondary photons on the top of the atmosphere. The detailed description of the PRESHOWER program is presented in [13]. Unfortunately, the article refers to the simulator of the UHE photon interaction with the geomagnetic field, so the simulator version applied in this work is different. However, this section and the indicated article provide sufficient knowledge to understand research methods. The PRESHOWER simulator generates sample data for determined initial conditions. Simulation parameters that could be changed were the primary photon energy, the impact parameter, and the heliocentric latitude of the primary photon trajectory. Therefore, only two coordinates of the simulation start point could be adjusted. What is more, the primary photon direction was completely determined after the choice of parameters. For these reasons, it was necessary to update the PRESHOWER simulator, with the purpose of applying the more realistic model. The introduced modifications and program functioning are described in the following sections.
The coordinate system used in the simulator is shown in Figure 1. The applied coordinate system is similar to the heliocentric ecliptic coordinate system, but in contrast to this system, the Earth motion is not simulated. Therefore, the x axis is determined by the centre of the Earth and the centre of the Sun. The Earth centre is also the coordinate system origin. The xy plane contains the ecliptic. In fact, the simulator neglects the Earth and its atmosphere curvature. The atmosphere flat is modelled by a plane with equation x = 112.83 km. On the scheme, simulation parameters -the impact parameter and the heliocentric latitude of the primary photon trajectory are indicated. The impact parameter in this model is the distance between the primary photon trajectory and the Sun centre. Whereas, the heliocentric latitude is defined as the angular distance between the ecliptic and a simulation start point in the heliocentric coordinate system. In this way, the simulation start point is completely determined, because the x coordinate of the start point is a constant value for all simulations equals 156.9 × 10 6 km, which is implemented in the code and cannot be adjusted. Moreover, the primary photon is always aimed at the centre of the atmosphere plane. Therefore, it is assigned after a simulation parameters choice. In this case, the heliocentric latitude remains unchangeable throughout the primary photon propagation. For proper functioning of the program, many functions with different destinations need to be applied, among others, the CERN routine DBSKA for calculating the modified Bessel functions. The most important function of the simulator performs a loop on particles and simulates physical processes involving them and propagates them. The loop includes only the primary UHE photon and possibly the primary pair, because secondary photons' directions are not altered and they do not need to be propagated step-wise. Thus, their position at the top of the atmosphere is calculated just after their generation and potential secondary photons cascading is neglected. The function calculates a simulation start point for given input parameters. Then, it propagates UHE photon and simulates its conversion gradually. In each step, the probability for conversion is computed according to the formula (3).
For a stochastic sampling of the probability for the process occurrence, the long period pseudo-random number generator is used. Generated numbers are consistent with the uniform distribution. Therefore, thanks to imposing the condition that the process occurs if the generated number is less than the calculated probability, the simulation reflects the probability distribution for the considered process. When the electron-positron pair is created, the fractional energy is simulated similarly to the conversion probability, with respect to its probability distribution (4). Then, the electron and positron are propagated step-wise according to the formula (5). The distance ∆s chosen in each step, which refers to the time interval ∆t, depends on the magnetic The purple vector t is the primary photon direction, which is randomly generated. The blue vector R points out the simulation start point, its length R is the simulation input parameter, whereas the polar angle θ between the vector R and the x axis and the azimuthal angle ϕ in the y z plane are randomly generated. The x p2 point determines the top of the atmosphere plane x = x p2 , but it is no longer the primary photon land point. field strength for the particle location in space and the current particle energy. Secondary photons emission is simulated by analogy with the primary photon conversion. The probability for this process is computed on the basis of the formula (6). And the approximation (7) is adopted as the angle of synchrotron photon deflection from the electron or positron trajectory. In turn, the angle of photon direction in the plane horizontal to the emitting particle trajectory is sampled from the uniform distribution. The simulation ends when the primary photon or created particles reach the top of the atmosphere. All obtained particle data is stored in a matrix, which has a fixed size. Therefore, the maximum number of particles generated in the simulation is limited.
On the basis of the simulation results, the spatial distribution of secondary photons at the top of the atmosphere may be analysed in detail, in the case of the UHE photon aimed at the top of the atmosphere. The influence of simulation input parameters on the distribution has been the subject of the previous study [11].
Previous results indicated that it is justified to consider a more realistic case, in which UHE photons reaching the Sun's vicinity from every direction are modelled. Furthermore, it was assumed that the UHE photons flux is isotropic. Nevertheless, this approach required substantial changes in the program's functioning. The most crucial modification concerned the arbitrary direction of the UHE photon, which implies the need for new simulation parameters and a new method of calculating the simulation start point. On the basis of assumptions made about the UHE photons flux, it was found that the most suitable for simulations is the heliocentric coordinate system. Figure 2 illustrates the modified coordinate system. It is worth mentioning that the inclination angle θ in the coordinate system is defined as an angle between the x axis and the position vector. In the new version of the PRESHOWER, simulation start points should be isotropically distributed at the sphere surface. Therefore, the radius of this start sphere was chosen as the simulation parameter. Whereas, the position on the start sphere is generated in such a way that every solid angle is equally probable. This condition may be described as where dN is the number of UHE photons at the start sphere in the solid angle dΩ. Therefore, the probability of generating the primary photon at the position with the azimuthal angle from interval [ϕ, ϕ + dϕ] should be equal where N is the total number of UHE photons reaching the start sphere. The probability density function of the azimuthal angles is constant and the variable is described by the uniform distribution. Similarly, one may obtain that f (cos θ) = const (10) and cos θ is uniformly distributed. Therefore, in the applied function ϕ and cos θ are randomly generated from uniform distribution and then the simulation start point is calculated using the standard transformation between spherical coordinates and Cartesian coordinates. In the transformation, the redefinitions of the θ angle and the Sun's position are included. The isotropic primary photon direction is modelled analogically. Nevertheless, the transformations from (θ,φ) coordinates to coordinates used in the simulation are more complex. Due to the fact that (θ,φ) are defined in the local coordinate system related to the position on the start sphere. Generated angles are firstly transformed to local Cartesian system ex = − sin θ cos ϕ,êỹ = − sin θ sin ϕ,êz = − cos θ.
The minus sign results from considering primary photons coming from the outside of the start sphere. The local versors may be identified as spherical versorŝ ez =ê r ,êx =ê θ ,êỹ =ê φ (11) related to the position on the start sphere (θ S , ϕ S ). Therefore, they are expressed in coordinates used in the simulator (Fig. 1) using formulaŝ e x =êz cos θ S −êx sin θ Ŝ e y =êz sin θ S cos ϕ S −êỹ sin ϕ S +êx cos θ S cos ϕ Ŝ e z =êỹ cos ϕ S +êx cos θ S sin ϕ S +êz sin θ S sin ϕ S .
Additionally, the drawn value cosθ belongs to the interval [0, 1], since we consider only the range of polar angles [0, π/2], which responds to directions beneath the sphere surface.
The generalisation of the program for an arbitrary primary photon direction required also changes in the simulation end condition. The previous condition may be effectively applied only to primary photons aimed at the atmosphere plane. Thus, assuming the plane is infinite, the previous condition is placed on photons with the negative component of the direction versor. In other cases, the simulation end condition is the distance to the Sun's centre larger than 7 R . What is more, the propagation of the primary photon, electron, or positron is aborted in case the particle lands on the Sun. In turn, secondary photons aimed at the Sun are not even included in the particle table. In the previous version of the program, these additional conditions were not needed, because of the specific determination of the primary photon direction. This approach enables considering also the potential special cases of the UHE photon cascading. For instance, it takes into account the primary photons not aimed at the Earth, that may generate observable secondary photons after conversion, thanks to the change of electron or positron direction in the magnetic field.

Results
Introduced amendments allowed to simulate the cumulative distribution of secondary photons at the top of the atmosphere. However, the new version of the simulator required to be verified in detail beforehand. Therefore, both the validation of the simulator and the preliminary analysis of results are described in the following sections.

Verification of the PRESHOWER functioning
It was significant to analyse simulation results firstly with respect to the correctness. This verification was performed both for appended functions and for the simulation as a whole. In the case of added functions, we visualised the way that examples of results model the phenomena. For this purpose, we analysed the distribution of generated solid angles and then visualised calculated positions at the sphere and directions. On the basis of 100 000 generated cos θ and ϕ values, it was found that the distribution of solid angles is flat. Similarly, the visualisation of position and direction distributions were in line with our expectations. With the purpose of verifying the whole program and the interoperability of its functions, it was necessary to compare results from the present version of the PRESHOWER with those from the previous version. However, this comparison should concern an identical physical event. Therefore, the simulation start point, the primary photon direction, and energy must be the same in both used simulator versions. We considered the case when the primary photon with energy 100 EeV aimed at the Earth propagates in the Sun's vicinity with the impact parameter 3 R and the heliocentric latitude 0 • . This responds to following parameters in the current PRESHOWER version: the simulation start point (R, θ S , φ S ) = (10.44 R , 0.29, 0), where R is the start sphere radius, and the primary photon direction (θ,φ) = (0.28, π). In this simple case the primary photon, electron, and positron are propagating only in the xy plane. Their movement modelled by both simulations is illustrated in Figure 3. The graph shows particle positions for each thousandth step of the simulation. It was found that the particles' propagation is similar in both cases. Minor differences such as distances between the next steps are caused by a different conversion point. This, in turn, results from the stochastic nature of physical processes. Furthermore, the distribution of secondary photons at the top of the atmosphere is presented in Figure 4. The generated distribution has a line-like shape and a size 8.48 × 10 4 km, which is a reasonable outcome. In accordance with previous results [11], the distribution size should be approximately 10 5 km. Therefore, obtained results are correct and it was concluded that the simulator properly models the UHE photon cascading.

Optimization of the simulation
The inclusion of more cases in the simulation, unfortunately, extended the computation time. It was found that the simulator needs to be optimised. Due to this fact, we reduced the cases included by the program and adjusted its parameters.      Figure 5. The graph of the number of simulations that fulfilled the land point limit in 10 10 simulations, executed for each limit, as a function of the land point limit, for different start sphere radius: 5 R , 6 R and 7 R

Reducing considered cases
First simulations showed that only a few of them gave rise to observable effects on the Earth. Thus, most of the calculation time was used for cases that are not the subject of this study. With the purpose of the effective disposal of computer resources, we resigned from modelling the special cases of the UHE photon cascading. Additional conditions were placed on the primary photon direction and its land point at the top of the atmosphere. Only primary photons moving towards the negative x axis are simulated. In turn, the primary photon land point is usually related to the centre of secondary photons distribution. In this way, limiting its distance from the x axis ( Figure 2) increases the probability of the observation on the Earth. The limit on the distance between the point at the top of the Earth's atmosphere where the UHE photon would have landed if no conversion had occurred, and the x axis connecting Earth's and Sun's centres is hereinafter referred to as the land point limit. Conditions are checked at the beginning of the program and further procedures are performed only for primary photons that fulfil limits. What is more, the number of secondary photons saved to the outfiles was limited to particles that are not further than 6500 km from the x axis at the top of the atmosphere.

Adjustment of simulation parameters
The implemented limit on the centre of the secondary photons spatial distribution and the start sphere radius -the simulation input parameter should be also adjusted for the optimization of the simulation. The limit on the primary photon land point has an influence on the ratio of observable on the Earth events to all simulations that fulfilled the condition. The number of simulations giving observable effects increases with the decrease in the distance from the x axis to the centre of secondary photons signature. On the other hand, the limit obviously affects the fraction of simulations that fulfil the condition and in this way the computation time. The strengthening of the limit on the primary photon land point results in the decline in the ratio of simulation fulfilling the condition to all running simulations and considerably prolongs the calculation time. These relations are shown in Figure 5 and Figure  6. On the basis of the obtained data, we calculated the ratio of observable on the Earth events to all running simulations, which is illustrated in Figure 7. It was observed that the fraction of detectable events is greater for the relaxation of the condition. However, the statistics for  this range is based on few entries, so it may be burdened with high risk. Figure 5 shows that the number of simulations fulfilling the condition depends weakly on the start sphere radius. Thus, the fraction of simulations observed on the Earth was calculated only for the start sphere radius 6 R ( Figure 6). In order to consider the largest number of events observable on the Earth in a relatively short time, the land point limit was set at 10 4 km. The number of CRE effects in all effects detectable on the Earth, in turn, is dependent on the start sphere radius. Due to the fact that the start sphere radius has an influence on primary photons' impact parameters, it also affects the probability of conversion. On the other hand, the decrease in the start sphere radius may also eliminate some cases that can give observable effects at the top of the atmosphere. Therefore the start sphere radius 6 R was chosen as an optimum value for the simulation parameter. The resulting cumulative distribution of secondary photons at the top of the atmosphere is shown in Figure 8. We enlarged the considered observable area to take into account the cross section of the Earth and the surrounding atmosphere. The distribution includes 324 CRE effect signatures, which corresponds to one year of observation taking into account present limits on the UHE photon flux [14]. Observed line-like signatures are parts of large cascades with centres beyond the Earth's cross-section. They justified the consideration of the cumulative distribution of secondary photons. Due to the fact that cascades are produced by UHE photons in the vicinity of the Sun, many secondary photons have directions contained in the solid angle determined by the Sun. Therefore CRE effect cascades should distinguish themselves from the rest of the observations, since secondary photons may reach higher energies than photons emitted by the Sun. The anomaly in the very-high-energy photons flux was already observed (see [19] and [20]). However, the UHE photon cascading was not taken into account as a possible source to date. Therefore, it is significant to analyse the cumulative distribution of secondary photons with regard to unique features that could help to classify them. The possible anisotropy may be observed in the distribution of the orientation of signatures at the top of the atmosphere. The horizontal or nearly horizontal orientation could be more frequent. It would be also theoretically motivated since signatures oriented in this way are generated by UHE photons that have the heliocentric latitude around 0 • or 90 • at the conversion point [1]. And for the heliocentric latitude 90 • , the conversion occurs more often [11]. In order to describe the orientation of signatures, the rotation angle parameter was introduced. It is given by the equation where points (y 1 , z 1 ) and (y 2 , z 2 ) are related to secondary photons that are the farthest, but still located in the observable area. The rotation angle is in the range of (−π/2, π/2). Obtained distribution is shown in the Figure 9.

1
The previous study ( [1,11]) has shown that the UHE photon interaction with the solar the interactions of the isotropic UHE photon flux with the solar magnetic field in order to 10 obtain for the first time particle distributions containing secondary photons from both central 11 and peripheral parts of the CRE cascades. The method is based on using the PRESHOWER 12 program appropriately modified and optimized. We verified the functioning of the simulator 13 and checked that obtained results are perfectly consistent with the previous analyses [11]. 14 In order to investigate the cumulative distribution of secondary photons at the top of the 15 atmosphere, we discussed the impact of simulation parameters on the modelled physical 16 phenomena and found optimal values for these parameters. The results obtained in this 17 way indicate that secondary photons produced by the UHE photons not aimed at the Earth 18 contribute significantly to the cumulative distribution of secondary photons at the top of the 19 atmosphere and to the distributions of CRE footprint orientations which might serve as a 20 valuable observational hint. In relation to the methods described here it is important to note 21 that recent observations ( [19,20]) demonstrate yet unexplained gamma-ray emission from 22 the solar disk. In the light of the simulations of CRE resulting from the interactions of UHE 23 photons nearby the Sun one might speculate that these observations may be a signature of 24 the CRE effect which motivates a dedicated analysis in this direction. We conclude that it is  The following abbreviations are used in this manuscript: