Photon tracking with GPUs in IceCube Chirkin, Dmitry
Nuclear instruments & methods in physics research. Section A, Accelerators, spectrometers, detectors and associated equipment,
10/2013, Volume:
725
Journal Article
Peer reviewed
GPUs (graphics processing units) have become increasingly popular in the recent years for scientific calculations involving large numbers of similar steps. Photon propagation is a necessary part of ...simulating detector response to passing charged particles in IceCube that is an ideal application for use with GPUs. We discuss the principle ideas and practical issues of running such an application within the simulation chain used within our collaboration.
Airborne radar has detected ≈100 lakes under the Antarctic ice cap, the largest of which is Lake Vostok. International planning is underway to search in Lake Vostok for microbial life that may have ...evolved in isolation from surface life for millions of years. It is thought, however, that the lakes may be hydraulically interconnected. If so, unsterile drilling would contaminate not just one but many of them. Here we report measurements of temperature vs. depth down to 2,345 m in ice at the South Pole, within 10 km from a subglacial lake seen by airborne radar profiling. We infer a temperature at the 2,810-m deep base of the South Pole ice and at the lake of -9°C, which is 7°C below the pressure-induced melting temperature of freshwater ice. To produce the strong radar signal, the frozen lake must consist of a mix of sediment and ice in a flat bed, formed before permanent Antarctic glaciation. It may, like Siberian and Antarctic permafrost, be rich in microbial life. Because of its hydraulic isolation, proximity to South Pole Station infrastructure, and analog to a Martian polar cap, it is an ideal place to test a sterile drill before risking contamination of Lake Vostok. From the semiempirical expression for strain rate vs. shear stress, we estimate shear vs. depth and show that the IceCube neutrino observatory will be able to map the three-dimensional ice-flow field within a larger volume (0.5 km3) and at lower temperatures (-20°C to -35°C) than has heretofore been possible.
The IceCube Neutrino Observatory is a cubic kilometer neutrino detector located at the geographic South Pole designed to detect high-energy astrophysical neutrinos. To thoroughly understand the ...detected neutrinos and their properties, the detector response to signal and background has to be modeled using Monte Carlo techniques. An integral part of these studies are the optical properties of the ice the observatory is built into. The simulated propagation of individual photons from particles produced by neutrino interactions in the ice can be greatly accelerated using graphics processing units (GPUs). In this paper, we (a collaboration between NVIDIA and IceCube) reduced the propagation time per photon by a factor of up to 3 on the same GPU. We achieved this by porting the OpenCL parts of the program to CUDA and optimizing the performance. This involved careful analysis and multiple changes to the algorithm. We also ported the code to NVIDIA OptiX to handle the collision detection. The hand-tuned CUDA algorithm turned out to be faster than OptiX. It exploits detector geometry and only a small fraction of photons ever travel close to one of the detectors.
The IceCube Neutrino Observatory instruments about 1 km3 of deep, glacial ice at the geographic South Pole. It uses 5160 photomultipliers to detect Cherenkov light emitted by charged relativistic ...particles. An unexpected light propagation effect observed by the experiment is an anisotropic attenuation, which is aligned with the local flow direction of the ice. We examine birefringent light propagation through the polycrystalline ice microstructure as a possible explanation for this effect. The predictions of a first-principles model developed for this purpose, in particular curved light trajectories resulting from asymmetric diffusion, provide a qualitatively good match to the main features of the data. This in turn allows us to deduce ice crystal properties. Since the wavelength of the detected light is short compared to the crystal size, these crystal properties include not only the crystal orientation fabric, but also the average crystal size and shape, as a function of depth. By adding small empirical corrections to this first-principles model, a quantitatively accurate description of the optical properties of the IceCube glacial ice is obtained. In this paper, we present the experimental signature of ice optical anisotropy observed in IceCube light-emitting diode (LED) calibration data, the theory and parameterization of the birefringence effect, the fitting procedures of these parameterizations to experimental data, and the inferred crystal properties.
GPGPU computing offers extraordinary increases in pure processing power for parallelizable applications. In IceCube we use GPUs for ray-tracing of cherenkov photons in the Antarctic ice as part of ...detector simulation. We report on how we implemented the mixed simulation production chain to include the processing on the GPGPU cluster for the IceCube Monte-Carlo production. We also present ideas to include GPGPU accelerated reconstructions into the IceCube data processing.
The IceCube detector, planned to reach 1 km3 in the next 3 years, is now 50% complete with 40 strings deployed in the ice and 40 IceTop stations installed on the surface. To improve the sensitivity ...of IceCube to high-energy neutrinos, the spacing between strings deployed in the detector's last construction phase may be increased. We consider several final detector geometry possibilities, and demonstrate that a substantial increase in the expected number of neutrino events with energies from 0.5 to 500 PeV, and a similar improvement for the signal expectation from GZK neutrinos, can be achieved.
The optical sensors of the IceCube Neutrino Observatory are attached on vertical strings of cables. They were frozen into the ice in the deployment holes made by hot water drill. This hole ice, to ...the best of our knowledge, consists of a bubbly central column, with the remainder of the re-frozen volume being optically clear. The bubbly ice often blocks one or several of the calibration LEDs in every optical sensor and significantly distorts the angular profile of the calibration light pulses. It also affects the sensors' response to in-coming photons at different locations and directions. We present our modeling of the hole ice optical properties as well as optical sensor location and orientation within the hole ice. The shadowing effects of cable string and possible optical sensor tilt away from the nominal vertical alignment are also discussed.
It is often not possible to construct a probability density function that describes the data. This can happen if there is no analytic description, and the number of parameters is too large so that it ...is impossible to simulate and tabulate all combinations. In these situations it is still interesting to rank simulation sets performed with different parameters in how well they compare to data. We propose a solution that appears to be better suited to this task than some of the obvious alternatives.
The IceCube Neutrino Observatory instruments about 1 km\(^3\) of deep, glacial ice at the geographic South Pole using 5160 photomultipliers to detect Cherenkov light from relativistic, charged ...particles. Most IceCube science goals rely on precise understanding and modelling of the optical properties of the instrumented ice. A peculiar light propagation effect observed by IceCube is an anisotropic attenuation, which is aligned with the local flow of the ice. Recent efforts have shown this effect is most likely due to curved photon trajectories resulting from the asymmetric light diffusion in the birefringent polycrystalline microstructure of the ice. This new model can be optimized by adjusting the average orientation, size and shape of the ice crystals. We present the parametrization of the birefringence effect in our photon propagation simulation, the fitting procedures and results. The anticipated potential of calibration instrumentation in the upcoming IceCube Upgrade to improve on known shortcomings of the current ice modelling is also discussed.
The IceCube Neutrino Observatory instruments about 1 km\(^3\) of deep, glacial ice at the geographic South Pole using 5160 photomultipliers to detect Cherenkov light of charged relativistic ...particles. Most of IceCube's science goals rely heavily on an ever more precise understanding of the optical properties of the instrumented ice. A curious light propagation effect observed by the experiment is an anisotropic attenuation, which is aligned with the local flow of the ice. Having recently identified curved photon trajectories resulting from asymmetric light diffusion in the birefringent polycrystalline microstructure of the ice as the most likely underlying cause of this effect, work is now ongoing to optimize the model parameters (effectively deducing the average crystal size and shape in the detector). We present the parametrization of the birefringence effect in our photon propagation simulation, the fitting procedures and results as well as the impact of the new ice model on data-MC agreement.