Physical Review D
covering particles, fields, gravitation, and cosmology
- Highlights
- Recent
- Accepted
- Collections
- Authors
- Referees
- Search
- Press
- About
- Editorial Team
Simulation of radio signals from cosmic-ray cascades in air and ice as observed by in-ice Askaryan radio detectors
Simon De Kockere, Dieder Van den Broeck, Uzair Abdul Latif, Krijn D. de Vries, Nick van Eijndhoven, Tim Huege, and Stijn Buitink
Phys. Rev. D 110, 023010 – Published 8 July 2024
- Article
- References
- No Citing Articles
PDFHTMLExport Citation
- Abstract
- Authors
- See Also
- Article Text
Abstract
A new generation of neutrino observatories will search for PeV-EeV neutrinos interacting in the ice by detecting radio pulses. Extended air showers propagating into the ice will form an important background and could be a valuable calibration signal. We present results from a Monte-Carlo simulation framework developed to fully simulate radio emission from cosmic-ray particle cascades as observed by in-ice radio detectors in the polar regions. The framework involves a modified version of coreas (a module of corsika7) to simulate in-air radio emission and a geant4-based framework for simulating in-ice radio emission from cosmic-ray showers as observed by in-ice antennas. The particles that reach the surface of the polar ice sheet at the end of the corsika7 simulation are injected into the geant4-based shower simulation code that takes the particles and propagates them further into the ice sheet, using an exponential density profile for the ice. The framework takes into account curved ray paths caused by the exponential refractive index profiles of air and ice. We present the framework and discuss some key features of the radio signal and radio shower footprint for in-ice observers.
11 More
- Received 25 March 2024
- Accepted 7 June 2024
DOI:https://doi.org/10.1103/PhysRevD.110.023010
© 2024 American Physical Society
Physics Subject Headings (PhySH)
- Research Areas
Particle astrophysicsParticle detection signatures
- Techniques
Cherenkov detectorsMulti-purpose particle detectorsNeutrino detection
Gravitation, Cosmology & AstrophysicsParticles & Fields
Authors & Affiliations
Simon De Kockere1,*, Dieder Van den Broeck2,†, Uzair Abdul Latif1, Krijn D. de Vries1,‡, Nick van Eijndhoven1, Tim Huege3,2, and Stijn Buitink2
- 1Vrije Universiteit Brussel, Dienst ELEM, IIHE, 1050 Brussels, Belgium
- 2Vrije Universiteit Brussel, Astrophysical Institute, IIHE, 1050 Brussels, Belgium
- 3Institut für Astroteilchenphysik, Karlsruher Institut für Technologie, Karlsruhe, Germany
- *Contact author: simon.de.kockere@vub.be
- †Contact author: dieder.jan.van.den.broeck@vub.be
- ‡Contact author: krijn.de.vries@vub.be
See Also
Simulation of in-ice cosmic ray air shower induced particle cascades
S. De Kockere, K. D. de Vries, N. van Eijndhoven, and U. A. Latif
Phys. Rev. D 106, 043023 (2022)
Article Text (Subscription Required)
Click to Expand
References (Subscription Required)
Click to Expand
Issue
Vol. 110, Iss. 2 — 15 July 2024
Access Options
- Buy Article »
- Log in with individual APS Journal Account »
- Log in with a username/password provided by your institution »
- Get access through a U.S. public or high school library »
![Simulation of radio signals from cosmic-ray cascades in air and ice as observed by in-ice Askaryan radio detectors (13) Simulation of radio signals from cosmic-ray cascades in air and ice as observed by in-ice Askaryan radio detectors (13)](https://i0.wp.com/cdn.journals.aps.org/development/journals/images/author-services-placard.png)
Images
Figure 1
A sketch of the three different types of rays: a direct ray (solid line), an indirect refracted ray (dotted line), and an indirect reflected ray (dashed line). In reality, only one of the two indirect rays will be a solution of the ray tracing.
Figure 2
A visualization of an interpolation table used for the in-ice ray tracing, for a given receiver position. Interpolation tables for the in-air ray tracing follow a similar structure, but instead cover a much larger area in air. A single in-air table is used for multiple receivers at the same depth in the ice.
Figure 3
Left: visualisation of several rays traced from the receiver to a line of emitters, including an air-ice boundary at . Direct rays are shown in blue. Reflected rays are only considered when total internal reflection occurs and are shown in red. Transmitted rays are shown in green. Right: the boost factor estimations along with the numerically computed value of as a function of depth. The diamonds correspond to the calculations using pointing directly from emission point to receiver point. The stars correspond to the values of the boost factor when evaluating at the emission point and using the launch direction for . The dots correspond to the numerically calculated values of the boost factor. The different regions are associated with direct (blue), reflected (red), and transmitted rays (green).
Figure 4
A visualization of the variables used in the end-point formula including ray tracing [Eq.(3)], in the specific case where the start point, end point, and receiver coincide in the same plane. The circles represent the start and end point of the emitter step; the square represents the receiver. All vectors shown in the illustration lie in the same plane.
Figure 5
A visualization of the relation between a global Cartesian coordinate system and a local spherical coordinate system, defined by the incoming direction of the ray.
Figure 6
The fluence footprint of the simulated cosmic-ray shower at a depth of 100m for the in-air emission (left), in-ice emission (middle), and the combined emission (right), using a primary energy of and zenith angle . The simulation was performed using 121 antennas in a star-shaped grid with eight arms and an antenna spacing of 10m, indicated by the white dots. For the interpolation of the star grid we used the code described in [57].
Figure 7
The electric field components as a function of time for four different antennas positioned along the west axis at a depth of 100m, using a primary energy of and zenith angle . The coordinate of each antenna is indicated at the top of the plots. The blue line indicates the in-air emission, and the red line indicates the in-ice emission. Note that the range on the y axis is different for every plot.
Figure 8
The electric field components as a function of time for four different antennas positioned along the north axis at a depth of 100m, using a primary energy of and zenith angle . The coordinate of each antenna is indicated at the top of the plots. The blue line indicates the in-air emission, and the red line indicates the in-ice emission. Note that the range on the y axis is different for every plot.
Figure 9
The fluence on the west axis of the simulated cosmic-ray shower at a depth of 100m for the in-ice emission only using a constant ice density of and a constant index of refraction , compared to the result using ray tracing with an exponential density and index of refraction profile. In both cases an eighth order digital Butterworth bandpass filter for a frequency band of 300–1000MHz was applied. The simulation was performed using 201 antennas using an antenna spacing of 1.5m, with a primary energy of and zenith angle .
Figure 10
The fluence on the north axis (left) and west axis (right) of the simulated cosmic-ray shower at a depth of 100m for the combined in-air and in-ice emission when applying an eighth order digital Butterworth bandpass filter, for the frequency bands 30–100MHz (upper), 100–300MHz (middle), and 300–1000MHz (lower). The simulation was performed using 201 antennas each on the north and west axis using an antenna spacing of 1.5m, with a primary energy of and zenith angle .
Figure 11
The fluence footprint of the simulated cosmic-ray shower at a depth of 100m for the in-air emission (left), in-ice emission (middle), and the combined emission (right), using a primary energy of and zenith angle . The simulation was performed using 121 antennas in a star-shaped grid with eight arms and an antenna spacing of 10m, indicated by the white dots. For the interpolation of the star grid we used the code described in [57].
Figure 12
The electric field components as a function of time for two different antennas positioned along the west axis at a depth of 100m, using a primary energy of and zenith angle . The coordinate of each antenna is indicated at the top of the plots. The blue line indicates the in-air emission, and the red line indicates the in-ice emission. Note that the range on the y axis is different for every plot.
Figure 13
The difference in arrival time of air and ice radio pulses as a function of distance to the shower axis for varying depths as indicated by the legend, using a primary energy and a zenith angle . The arrival time was defined as the time where the Hilbert envelope reaches 33% of its maximum value, after passing an eighth order digital Butterworth bandpass filter for the frequency band 30–1000MHz, averaged over the three components.
Figure 14
The number of particles (left) and the distribution of the energy (right) in the air shower as a function of depth with primary energy and zenith angle . The particle distributions are obtained over the full radial extent of the in-air particle cascade. Similar results were shown in [22].
Figure 15
The number of particles (left) and the distribution of the energy (right) in the air shower as a function of depth with primary energy and zenith angle . The particle distributions are obtained over the full radial extent of the in-air particle cascade. Similar results were shown in [22].
Figure 16
The total energy contained within a given radius from the shower core at the air-ice boundary up to 100m (left) and in more detail up to 5m (right), for the simulated shower with a primary energy and zenith angle and the simulated shower with a primary energy and zenith angle . The energy is expressed as a fraction of the primary energy . The air-ice boundary is located at an altitude of 2.835km, which corresponds to a vertical atmospheric depth of .
Figure 17
The energy deposited in ice by the in-ice particle cascade for the simulated shower with a primary energy and zenith angle (left) and the simulated shower with a primary energy and zenith angle (right). Shown here is the deposited energy density within a vertical 1-cm wide slice going through the center of the particle shower. The air-ice boundary is located at an altitude of 2.835km, which corresponds to a vertical atmospheric depth of . Similar results were shown in [22].
Figure 18
A distribution of the CPU time of 591 cores used for the simulation of the in-ice particle cascade and the corresponding radio emission for 121 antennas in the ice, using a primary energy and zenith angle .