To help understand our universe better, researchers and scientists currently run extreme-scale cosmology simulations on leadership supercomputers. However, such simulations can generate large amounts of scientific data, which often result in expensive costs in data associated with data movement and storage. Lossy compression techniques have become attractive because they significantly reduce data size and can maintain high data fidelity for post-analysis. In this paper, we propose to use GPU-based lossy compression for extreme-scale cosmological simulations. Our contributions are threefold: (1) we implement multiple GPU-based lossy compressors to our opensource compression benchmark and analysis framework named Foresight; (2) we use Foresight to comprehensively evaluate the practicality of using GPU-based lossy compression on two real-world extreme-scale cosmology simulations, namely HACC and Nyx, based on a series of assessment metrics; and (3) we develop a general optimization guideline on how to determine the best-fit configurations for different lossy compressors and cosmological simulations. Experiments show that GPU-based lossy compression can provide necessary accuracy on post-analysis for cosmological simulations and high compression ratio of 5~15x on the tested datasets, as well as much higher compression and decompression throughput than CPU-based compressors.

The advent of large cosmological sky surveys - ushering in the era of precision cosmology - has been accompanied by ever larger cosmological simulations. The analysis of these simulations, which currently encompass tens of billions of particles and up to trillion particles in the near future, is often as daunting as carrying out the simulations in the first place. Therefore, the development of very efficient analysis tools combining qualitative and quantitative capabilities is a matter of some urgency. In this paper we introduce new analysis features implemented within ParaView, a parallel, open-source visualization toolkit, to analyze large N-body simulations. The new features include particle readers and a very efficient halo finder which identifies friends-of-friends halos and determines common halo properties. In combination with many other functionalities already existing within ParaView, such as histogram routines or interfaces to Python, this enhanced version enables fast, interactive, and convenient analyses of large cosmological simulations. In addition, development paths are available for future extensions.

Current and upcoming cosmological observations allow us to probe structures on smaller and smaller scales, entering highly nonlinear regimes. In order to obtain theoretical predictions in these regimes, large cosmological simulations have to be carried out. The promised high accuracy from observations make the simulation task very demanding: the simulations have to be at least as accurate as the observations. This requirement can only be fulfilled by carrying out an extensive code validation program. The first step of such a program is the comparison of different cosmology codes including gravitation interactions only. In this paper we extend a recently carried out code comparison project to include five more simulation codes. We restrict our analysis to a small cosmological volume which allows us to investigate properties of halos. For the matter power spectrum and the mass function, the previous results hold, with the codes agreeing at the 10% level over wide dynamic ranges. We extend our analysis to the comparison of halo profiles and investigate the halo count as a function of local density. We introduce and discuss ParaView as a flexible analysis tool for cosmological simulations, the use of which immensely simplifies the code comparison task.

The Antarctic Muon And Neutrino Detector Array (AMANDA) is a high-energy neutrino telescope operating at the geographic South Pole. It is a lattice of photo-multiplier tubes buried deep in the polar ice between 1500m and 2000m. The primary goal of this detector is to discover astrophysical sources of high energy neutrinos. A high-energy muon neutrino coming through the earth from the Northern Hemisphere can be identified by the secondary muon moving upward through the detector. The muon tracks are reconstructed with a maximum likelihood method. It models the arrival times and amplitudes of Cherenkov photons registered by the photo-multipliers. This paper describes the different methods of reconstruction, which have been successfully implemented within AMANDA. Strategies for optimizing the reconstruction performance and rejecting background are presented. For a typical analysis procedure the direction of tracks are reconstructed with about 2 degree accuracy.

We present the results of a search for point sources of high energy neutrinos in the northern hemisphere using AMANDA-II data collected in the year 2000. Included are flux limits on several AGN blazars, microquasars, magnetars and other candidate neutrino sources. A search for excesses above a random background of cosmic-ray-induced atmospheric neutrinos and misreconstructed downgoing cosmic-ray muons reveals no statistically significant neutrino point sources. We show that AMANDA-II has achieved the sensitivity required to probe known TeV gamma-ray sources such as the blazar Markarian 501 in its 1997 flaring state at a level where neutrino and gamma-ray fluxes are equal.

We present the results of a Monte-Carlo study of the sensitivity of the planned IceCube detector to predicted fluxes of muon neutrinos at TeV to PeV energies. A complete simulation of the detector and data analysis is used to study the detector's capability to search for muon neutrinos from sources such as active galaxies and gamma-ray bursts. We study the effective area and the angular resolution of the detector as a function of muon energy and angle of incidence. We present detailed calculations of the sensitivity of the detector to both diffuse and pointlike neutrino emissions, including an assessment of the sensitivity to neutrinos detected in coincidence with gamma-ray burst observations. After three years of datataking, IceCube will have been able to detect a point source flux of E^2*dN/dE = 7*10^-9 cm^-2s^-1GeV at a 5-sigma significance, or, in the absence of a signal, place a 90% c.l. limit at a level E^2*dN/dE = 2*10^-9 cm^-2s^-1GeV. A diffuse E-2 flux would be detectable at a minimum strength of E^2*dN/dE = 1*10^-8 cm^-2s^-1sr^-1GeV. A gamma-ray burst model following the formulation of Waxman and Bahcall would result in a 5-sigma effect after the observation of 200 bursts in coincidence with satellite observations of the gamma-rays.

This paper describes the search for astronomical sources of high-energy neutrinos using the AMANDA-B10 detector, an array of 302 photomultiplier tubes, used for the detection of Cherenkov light from upward traveling neutrino-induced muons, buried deep in ice at the South Pole. The absolute pointing accuracy and angular resolution were studied by using coincident events between the AMANDA detector and two independent telescopes on the surface, the GASP air Cherenkov telescope and the SPASE extensive air shower array. Using data collected from April to October of 1997 (130.1 days of livetime), a general survey of the northern hemisphere revealed no statistically significant excess of events from any direction. The sensitivity for a flux of muon neutrinos is based on the effective detection area for through-going muons. Averaged over the Northern sky, the effective detection area exceeds 10,000 m^2 for E_{mu} ~ 10 TeV. Neutrinos generated in the atmosphere by cosmic ray interactions were used to verify the predicted performance of the detector. For a source with a differential energy spectrum proportional to E_{nu}^{-2} and declination larger than +40 degrees, we obtain E^2(dN_{nu}/dE) <= 10^{-6}GeVcm^{-2}s^{-1} for an energy threshold of 10 GeV.

We report on a search for electromagnetic and/or hadronic showers (cascades) induced by a diffuse flux of neutrinos with energies between 5 TeV and 300 TeV from extraterrestrial sources. Cascades may be produced by matter interactions of all flavors of neutrinos, and contained cascades have better energy resolution and afford better background rejection than through-going \nu_\mu-induced muons. Data taken in 1997 with the AMANDA detector were searched for events with a high-energy cascade-like signature. The observed events are consistent with expected backgrounds from atmospheric neutrinos and catastrophic energy losses from atmospheric muons. Effective volumes for all flavors of neutrinos, which allow the calculation of limits for any neutrino flux model, are presented. The limit on cascades from a diffuse flux of \nu_e+\nu_\mu+\nu_\tau+\bar{\nu}_e+\bar{\nu}_\mu+\bar{\nu}_\tau is E^2 d\Phi/dE < 9.8 x 10^{-6} GeV/(cm^2 s sr), assuming a neutrino flavor flux ratio of 1:1:1 at the detector. The limit on cascades from a diffuse flux of \nu_e+\bar{\nu}_e is E^ d\Phi/dE < 6.5 x 10{-6} GeV/(cm^2 s sr), independent of the assumed neutrino flavor flux ratio.

A search for nearly vertical up-going muon-neutrinos from neutralino annihilations in the center of the Earth has been performed with the AMANDA-B10 neutrino detector. The data sample collected in 130.1 days of live-time in 1997, ~10^9 events, has been analyzed for this search. No excess over the expected atmospheric neutrino background is oberved. An upper limit at 90% confidence level on the annihilation rate of neutralinos in the center of the Earth is obtained as a function of the neutralino mass in the range 100 GeV-5000 GeV, as well as the corresponding muon flux limit.

The Antarctic Muon and Neutrino Detector Array (AMANDA) began collecting data with ten strings in 1997. Results from the first year of operation are presented. Neutrinos coming through the Earth from the Northern Hemisphere are identified by secondary muons moving upward through the array. Cosmic rays in the atmosphere generate a background of downward moving muons, which are about 10^6 times more abundant than the upward moving muons. Over 130 days of exposure, we observed a total of about 300 neutrino events. In the same period, a background of 1.05*10^9 cosmic ray muon events was recorded. The observed neutrino flux is consistent with atmospheric neutrino predictions. Monte Carlo simulations indicate that 90% of these events lie in the energy range 66 GeV to 3.4 TeV. The observation of atmospheric neutrinos consistent with expectations establishes AMANDA-B10 as a working neutrino telescope.