|
|
Wagner, C., Verde, L., & Boubekeur, L. (2010). N-body simulations with generic non-Gaussian initial conditions I: power spectrum and halo mass function. J. Cosmol. Astropart. Phys., 10(10), 022–24pp.
Abstract: We address the issue of setting up generic non-Gaussian initial conditions for N-body simulations. We consider inflationary-motivated primordial non-Gaussianity where the perturbations in the Bardeen potential are given by a dominant Gaussian part plus a non-Gaussian part specified by its bispectrum. The approach we explore here is suitable for any bispectrum, i.e. it does not have to be of the so-called separable or factorizable form. The procedure of generating a non-Gaussian field with a given bispectrum (and a given power spectrum for the Gaussian component) is not univocal, and care must be taken so that higher-order corrections do not leave a too large signature on the power spectrum. This is so far a limiting factor of our approach. We then run N-body simulations for the most popular inflationary-motivated non-Gaussian shapes. The halo mass function and the non-linear power spectrum agree with theoretical analytical approximations proposed in the literature, even if they were so far developed and tested only for a particular shape (the local one). We plan to make the simulations outputs available to the community via the non-Gaussian simulations comparison project web site http://icc.ub.edu/similar to liciaverde/NGSCP.html.
|
|
|
|
ATLAS Collaboration(Aaboud, M. et al), Alvarez Piqueras, D., Bailey, A. J., Barranco Navarro, L., Cabrera Urban, S., Castillo Gimenez, V., et al. (2018). Comparison between simulated and observed LHC beam backgrounds in the ATLAS experiment at E-beam=4 TeV. J. Instrum., 13, P12006–41pp.
Abstract: Results of dedicated Monte Carlo simulations of beam-induced background (BIB) in the ATLAS experiment at the Large Hadron Collider (LHC) are presented and compared with data recorded in 2012. During normal physics operation this background arises mainly from scattering of the 4 TeV protons on residual gas in the beam pipe. Methods of reconstructing the BIB signals in the ATLAS detector, developed and implemented in the simulation chain based on the FLUKA Monte Carlo simulation package, are described. The interaction rates are determined from the residual gas pressure distribution in the LHC ring in order to set an absolute scale on the predicted rates of BIB so that they can be compared quantitatively with data. Through these comparisons the origins of the BIB leading to different observables in the ATLAS detectors are analysed. The level of agreement between simulation results and BIB measurements by ATLAS in 2012 demonstrates that a good understanding of the origin of BIB has been reached.
|
|
|
|
Guadilla, V. et al, Tain, J. L., Algora, A., Agramunt, J., Gelletly, W., Jordan, D., et al. (2018). Characterization and performance of the DTAS detector. Nucl. Instrum. Methods Phys. Res. A, 910, 79–89.
Abstract: DTAS is a segmented total absorption y-ray spectrometer developed for the DESPEC experiment at FAIR. It is composed of up to eighteen NaI(Tl) crystals. In this work we study the performance of this detector with laboratory sources and also under real experimental conditions. We present a procedure to reconstruct offline the sum of the energy deposited in all the crystals of the spectrometer, which is complicated by the effect of NaI(Tl) light-yield non-proportionality. The use of a system to correct for time variations of the gain in individual detector modules, based on a light pulse generator, is demonstrated. We describe also an event-based method to evaluate the summing-pileup electronic distortion in segmented spectrometers. All of this allows a careful characterization of the detector with Monte Carlo simulations that is needed to calculate the response function for the analysis of total absorption gamma-ray spectroscopy data. Special attention was paid to the interaction of neutrons with the spectrometer, since they are a source of contamination in studies of beta-delayed neutron emitting nuclei.
|
|
|
|
Kuo, J. L., Lattanzi, M., Cheung, K., & Valle, J. W. F. (2018). Decaying warm dark matter and structure formation. J. Cosmol. Astropart. Phys., 12(12), 026–24pp.
Abstract: We examine the cosmology of warm dark matter (WDM), both stable and decaying, from the point of view of structure formation. We compare the matter power spectrum associated to WDM masses of 1.5 keV and 0.158 keV, with that expected for the stable cold dark matter ACDM Xi SCDM paradigm, taken as our reference model. We scrutinize the effects associated to the warm nature of dark matter, as well as the fact that it decays. The decaying warm dark matter (DWDM) scenario is well-motivated, emerging in a broad class of particle physics theories where neutrino masses arise from the spontaneous breaking of a continuous global lepton number symmetry. The majoron arises as a Nambu-Goldstone boson, and picks up a mass from gravitational effects, that explicitly violate global symmetries. The majoron necessarily decays to neutrinos, with an amplitude proportional to their tiny mass, which typically gives it cosmologically long lifetimes. Using N-body simulations we show that our DWDM picture leads to a viable alternative to the ACDM scenario, with predictions that can differ substantially on small scales.
|
|
|
|
Olleros, P., Caballero, L., Domingo-Pardo, C., Babiano, V., Ladarescu, I., Calvo, D., et al. (2018). On the performance of large monolithic LaCl3(Ce) crystals coupled to pixelated silicon photosensors. J. Instrum., 13, P03014–17pp.
Abstract: We investigate the performance of large area radiation detectors, with high energy-and spatial-resolution, intended for the development of a Total Energy Detector with gamma-ray imaging capability, so-called i-TED. This new development aims for an enhancement in detection sensitivity in time-of-flight neutron capture measurements, versus the commonly used C6D6 liquid scintillation total-energy detectors. In this work, we study in detail the impact of the readout photosensor on the energy response of large area (50 x 50 mm(2)) monolithic LaCl3(Ce) crystals, in particular when replacing a conventional mono-cathode photomultiplier tube by an 8 x 8 pixelated silicon photomultiplier. Using the largest commercially available monolithic SiPM array (25 cm(2)), with a pixel size of 6 x 6 mm(2), we have measured an average energy resolution of 3.92% FWHM at 662 keV for crystal thick-nesses of 10, 20 and 30 mm. The results are confronted with detailed Monte Carlo (MC) calculations, where optical processes and properties have been included for the reliable tracking of the scintillation photons. After the experimental validation of the MC model, we use our MC code to explore the impact of a smaller photosensor segmentation on the energy resolution. Our optical MC simulations predict only a marginal deterioration of the spectroscopic performance for pixels of 3 x 3 mm(2).
|
|
|
|
Caron, S., Gomez-Vargas, G. A., Hendriks, L., & Ruiz de Austri, R. (2018). Analyzing gamma rays of the Galactic Center with deep learning. J. Cosmol. Astropart. Phys., 05(5), 058–24pp.
Abstract: We present the application of convolutional neural networks to a particular problem in gamma ray astronomy. Explicitly, we use this method to investigate the origin of an excess emission of GeV gamma rays in the direction of the Galactic Center, reported by several groups by analyzing Fermi-LAT data. Interpretations of this excess include gamma rays created by the annihilation of dark matter particles and gamma rays originating from a collection of unresolved point sources, such as millisecond pulsars. We train and test convolutional neural networks with simulated Fermi-LAT images based on point and diffuse emission models of the Galactic Center tuned to measured gamma ray data. Our new method allows precise measurements of the contribution and properties of an unresolved population of gamma ray point sources in the interstellar diffuse emission model. The current model predicts the fraction of unresolved point sources with an error of up to 10% and this is expected to decrease with future work.
|
|
|
|
ATLAS Collaboration(Aaboud, M. et al), Alvarez Piqueras, D., Barranco Navarro, L., Cabrera Urban, S., Castillo Gimenez, V., Cerda Alberich, L., et al. (2017). Study of the material of the ATLAS inner detector for Run 2 of the LHC. J. Instrum., 12, P12009–59pp.
Abstract: The ATLAS inner detector comprises three different sub-detectors: the pixel detector, the silicon strip tracker, and the transition-radiation drift-tube tracker. The Insertable B-Layer, a new innermost pixel layer, was installed during the shutdown period in 2014, together with modifications to the layout of the cables and support structures of the existing pixel detector. The material in the inner detector is studied with several methods, using a low-luminosity root s = 13 TeV pp collision sample corresponding to around 2.0 nb(-1) collected in 2015 with the ATLAS experiment at the LHC. In this paper, the material within the innermost barrel region is studied using reconstructed hadronic interaction and photon conversion vertices. For the forward rapidity region, the material is probed by a measurement of the efficiency with which single tracks reconstructed from pixel detector hits alone can be extended with hits on the track in the strip layers. The results of these studies have been taken into account in an improved description of the material in the ATLAS inner detector simulation, resulting in a reduction in the uncertainties associated with the charged-particle reconstruction efficiency determined from simulation.
|
|
|
|
de Salas, P. F., Gariazzo, S., Lesgourgues, J., & Pastor, S. (2017). Calculation of the local density of relic neutrinos. J. Cosmol. Astropart. Phys., 09(9), 034–24pp.
Abstract: Nonzero neutrino masses are required by the existence of flavour oscillations, with values of the order of at least 50 meV. We consider the gravitational clustering of relic neutrinos within the Milky Way, and used the N – one-body simulation technique to compute their density enhancement factor in the neighbourhood of the Earth with respect to the average cosmic density. Compared to previous similar studies, we pushed the simulation down to smaller neutrino masses, and included an improved treatment of the baryonic and dark matter distributions in the Milky Way. Our results are important for future experiments aiming at detecting the cosmic neutrino background, such as the Princeton Tritium Observatory for Light, Early-universe, Massive-neutrino Yield (PTOLEMY) proposal. We calculate the impact of neutrino clustering in the Milky Way on the expected event rate for a PTOLEMY-like experiment. We find that the effect of clustering remains negligible for the minimal normal hierarchy scenario, while it enhances the event rate by 10 to 20% (resp. a factor 1.7 to 2.5) for the minimal inverted hierarchy scenario (resp. a degenerate scenario with 150 meV masses). Finally we compute the impact on the event rate of a possible fourth sterile neutrino with a mass of 1.3 eV.
|
|
|
|
NEXT Collaboration(Azevedo, C. D. R. et al), Gomez-Cadenas, J. J., Alvarez, V., Benlloch-Rodriguez, J. M., Botas, A., Carcel, S., et al. (2018). Microscopic simulation of xenon-based optical TPCs in the presence of molecular additives. Nucl. Instrum. Methods Phys. Res. A, 877, 157–172.
Abstract: We introduce a simulation framework for the transport of high and low energy electrons in xenon-based optical time projection chambers (OTPCs). The simulation relies on elementary cross sections (electron-atom and electron-molecule) and incorporates, in order to compute the gas scintillation, the reaction/quenching rates (atom-atom and atom-molecule) of the first 41 excited states of xenon and the relevant associated excimers, together with their radiative cascade. The results compare positively with observations made in pure xenon and its mixtures with CO2 and CF4 in a range of pressures from 0.1 to 10 bar. This work sheds some light on the elementary processes responsible for the primary and secondary xenon-scintillation mechanisms in the presence of additives, that are of interest to the OTPC technology.
|
|
|
|
Guadilla, V. et al, Algora, A., Tain, J. L., Agramunt, J., Jordan, D., Monserrate, M., et al. (2017). Characterization of a cylindrical plastic beta-detector with Monte Carlo simulations of optical photons. Nucl. Instrum. Methods Phys. Res. A, 854, 134–138.
Abstract: In this work we report on the Monte Carlo study performed to understand and reproduce experimental measurements of a new plastic beta-detector with cylindrical geometry. Since energy deposition simulations differ from the experimental measurements for such a geometry, we show how the simulation of production and transport of optical photons does allow one to obtain the shapes of the experimental spectra. Moreover, taking into account the computational effort associated with this kind of simulation, we develop a method to convert the simulations of energy deposited into light collected, depending only on the interaction point in the detector. This method represents a useful solution when extensive simulations have to be done, as in the case of the calculation of the response function of the spectrometer in a total absorption gamma-ray spectroscopy analysis.
|
|
