|
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.
|
|
|
Caron, S., Eckner, C., Hendriks, L., Johannesson, G., Ruiz de Austri, R., & Zaharijas, G. (2023). Mind the gap: the discrepancy between simulation and reality drives interpretations of the Galactic Center Excess. J. Cosmol. Astropart. Phys., 06(6), 013–56pp.
Abstract: The Galactic Center Excess (GCE) in GeV gamma rays has been debated for over a decade, with the possibility that it might be due to dark matter annihilation or undetected point sources such as millisecond pulsars (MSPs). This study investigates how the gamma-ray emission model (-yEM) used in Galactic center analyses affects the interpretation of the GCE's nature. To address this issue, we construct an ultra-fast and powerful inference pipeline based on convolutional Deep Ensemble Networks. We explore the two main competing hypotheses for the GCE using a set of-yEMs with increasing parametric freedom. We calculate the fractional contribution (fsrc) of a dim population of MSPs to the total luminosity of the GCE and analyze its dependence on the complexity of the ryEM. For the simplest ryEM, we obtain fsrc = 0.10 f 0.07, while the most complex model yields fsrc = 0.79 f 0.24. In conclusion, we find that the statement about the nature of the GCE (dark matter or not) strongly depends on the assumed ryEM. The quoted results for fsrc do not account for the additional uncertainty arising from the fact that the observed gamma-ray sky is out-of-distribution concerning the investigated ryEM iterations. We quantify the reality gap between our ryEMs using deep-learning-based One-Class Deep Support Vector Data Description networks, revealing that all employed ryEMs have gaps to reality. Our study casts doubt on the validity of previous conclusions regarding the GCE and dark matter, and underscores the urgent need to account for the reality gap and consider previously overlooked “out of domain” uncertainties in future interpretations.
|
|
|
ATLAS Collaboration(Aaboud, M. et al), Alvarez Piqueras, D., Aparisi Pozo, J. A., Bailey, A. J., Barranco Navarro, L., Cabrera Urban, S., et al. (2019). Modelling radiation damage to pixel sensors in the ATLAS detector. J. Instrum., 14, P06012–52pp.
Abstract: Silicon pixel detectors are at the core of the current and planned upgrade of the ATLAS experiment at the LHC. Given their close proximity to the interaction point, these detectors will be exposed to an unprecedented amount of radiation over their lifetime. The current pixel detector will receive damage from non-ionizing radiation in excess of 10(15) 1 MeV n(eq)/cm(2), while the pixel detector designed for the high-luminosity LHC must cope with an order of magnitude larger fluence. This paper presents a digitization model incorporating effects of radiation damage to the pixel sensors. The model is described in detail and predictions for the charge collection efficiency and Lorentz angle are compared with collision data collected between 2015 and 2017 (<= 10(15) 1 MeV n(eq)/cm(2)).
|
|
|
ATLAS Collaboration(Aad, G. et al), Alvarez Piqueras, D., Cabrera Urban, S., Castillo Gimenez, V., Costa, M. J., Fernandez Martinez, P., et al. (2015). Modelling Z -> ττ processes in ATLAS with τ-embedded Z -> μμ data. J. Instrum., 10, P09018–41pp.
Abstract: This paper describes the concept, technical realisation and validation of a largely data-driven method to model events with Z -> tau tau decays. In Z -> μμevents selected from proton-proton collision data recorded at root s = 8 TeV with the ATLAS experiment at the LHC in 2012, the Z decay muons are replaced by tau leptons from simulated Z -> tau tau decays at the level of reconstructed tracks and calorimeter cells. The tau lepton kinematics are derived from the kinematics of the original muons. Thus, only the well-understood decays of the Z boson and tau leptons as well as the detector response to the tau decay products are obtained from simulation. All other aspects of the event, such as the Z boson and jet kinematics as well as effects from multiple interactions, are given by the actual data. This so-called tau-embedding method is particularly relevant for Higgs boson searches and analyses in tau tau final states, where Z -> tau tau decays constitute a large irreducible background that cannot be obtained directly from data control samples. In this paper, the relevant concepts are discussed based on the implementation used in the ATLAS Standard Model H -> tau tau analysis of the full datataset recorded during 2011 and 2012.
|
|
|
Jaworski, G., Palacz, M., Nyberg, J., de Angelis, G., de France, G., Di Nitto, A., et al. (2012). Monte Carlo simulation of a single detector unit for the neutron detector array NEDA. Nucl. Instrum. Methods Phys. Res. A, 673, 64–72.
Abstract: A study of the dimensions and performance of a single detector of the future neutron detector array NEDA was performed by means of Monte Carlo simulations, using GEANT4. Two different liquid scintillators were evaluated: the hydrogen based BC501A and the deuterated BC537. The efficiency and the probability that one neutron will trigger a signal in more than one detector were investigated as a function of the detector size. The simulations were validated comparing the results to experimental measurements performed with two existing neutron detectors, with different geometries, based on the liquid scintillator BC501.
|
|
|
Guerrero, C., Cano-Ott, D., Mendoza, E., Tain, J. L., Algora, A., Berthoumieux, E., et al. (2012). Monte Carlo simulation of the n_TOF Total Absorption Calorimeter. Nucl. Instrum. Methods Phys. Res. A, 671, 108–117.
Abstract: The n_TOF Total Absorption Calorimeter (TAC) is a 4 pi BaF2 segmented detector used at CERN for measuring neutron capture cross-sections of importance for the design of advanced nuclear reactors. This work presents the simulation code that has been developed in GEANT4 for the accurate determination of the detection efficiency of the TAC for neutron capture events. The code allows to calculate the efficiency of the TAC for every neutron capture state, as a function of energy, crystal multiplicity, and counting rate. The code includes all instrumental effects such as the single crystal detection threshold and energy resolution, finite size of the coincidence time window, and signal pile-up. The results from the simulation have been validated with experimental data for a large set of electromagnetic de-excitation patterns: beta-decay of well known calibration sources, neutron capture reactions in light nuclei with well known level schemes like Ti-nat, reference samples used in (n,gamma) measurements like Au-197 and experimental data from an actinide sample like Pu-240. The systematic uncertainty in the determination of the detection efficiency has been estimated for all the cases. As a representative example, the accuracy reached for the case of Au-197(n,gamma) ranges between 0.5% and 2%, depending on the experimental and analysis conditions. Such a value matches the high accuracy required for the nuclear cross-section data needed in advanced reactor design.
|
|
|
Hornillos, M. B. G., Gorlychev, V., Caballero, R., Cortes, G., Poch, A., Pretel, C., et al. (2011). Monte Carlo Simulations for the Study of a Moderated Neutron Detector. J. Korean Phys. Soc., 59(2), 1573–1576.
Abstract: This work presents the Monte Carlo simulations performed with the MCNPX and GEANT4 codes for the design of a BEta deLayEd Neutron detector, BELEN-20. This detector will be used for the study of beta delayed neutron emission and consists of a block of polyethylene with dimensions 90 x 90 x 80 cm(3) and 20 cylindrical (3)He gas counters. The results of these simulations have been validated experimentally with a (252)Cf source in the laboratory at UPC, Barcelona. Also the first experiment with this detector has been carried out in November 2009 in JYFL, Finland. In this experiment the neutron emission probability after beta decay of the fission products (88)Br, (94,95)Rb, and (138)I has been measured; this data is still under analysis. Simulations with MCNPX and GEANT4 have been performed in order to obtain the efficiency of the BELEN-20 detector for each of the above nuclei using the neutron energy distribution corresponding to each nucleus.
|
|
|
Zhang, X., Xiao, Y. T., & Gimeno, B. (2020). Multipactor Suppression by a Resonant Static Magnetic Field on a Dielectric Surface. IEEE Trans. Electron Devices, 67(12), 5723–5728.
Abstract: In this article, we study the suppression of the multipactor phenomenon on a dielectric surface by a resonant static magnetic field. A homemade Monte Carlo algorithm is developed for multipactor simulations on a dielectric surface driven by two orthogonal radio frequency (RF) electric field components. When the static magnetic field is perpendicular to the tangential and normal RF electric fields, it is shown that if the normal electric field lags the tangential electric field by pi/2, the superposition of the normal and tangential electric fields will trigger a gyro-acceleration of the electron cloud and restrain the multipactor discharge effectively. By contrast, when the normal electric field is in advance of the tangential electric field by pi/2, the difference between the normal and tangential electric fields drives gyro-motion of the electron cloud. Consequently, two enhanced discharge zones are inevitable. The suppression effects of the resonant static magnetic field that is parallel to the tangential RF electric field or to the normal RF electric field are also presented.
|
|
|
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.
|
|
|
Balaudo, A., Calore, F., De Romeri, V., & Donato, F. (2024). NAJADS: a self-contained framework for the direct determination of astrophysical J-factors. J. Cosmol. Astropart. Phys., 02(2), 001–33pp.
Abstract: Cosmological simulations play a pivotal role in understanding the properties of the dark matter (DM) distribution in both galactic and galaxy -cluster environments. The characterization of DM structures is crucial for informing indirect DM searches, aiming at the detection of the annihilation (or decay) products of DM particles. A fundamental quantity in these analyses is the astrophysical J -factor. In the DM phenomenology community, J -factors are typically computed through the semi -analytical modelling of the DM mass distribution, which is affected by large uncertainties. With the scope of addressing and possibly reducing these uncertainties, we present NAJADS, a self-contained framework to derive the DM J -factor directly from the raw simulations data. We show how this framework can be used to compute all -sky maps of the J -factor, automatically accounting for the complex 3D structure of the simulated halos and for the boosting of the signal due to the density fluctuations along the line of sight. After validating our code, we present a proof -of -concept application of NAJADS to a realistic halo from the IllustrisTNG suite, and exploit it to make a thorough comparison between our numerical approach and traditional semi -analytical methods. JCAP02(2024)001
|
|