|
Jordan, D., Algora, A., & Tain, J. L. (2016). An event generator for simulations of complex beta-decay experiments. Nucl. Instrum. Methods Phys. Res. A, 828, 52–57.
Abstract: This article describes a Monte Carlo event generator for the design, optimization and performance characterization of beta decay spectroscopy experimental set-ups. The event generator has been developed within the Geant4 simulation architecture and provides new features and greater flexibility in comparison with the current available decay generator.
|
|
|
de los Rios, M., Petac, M., Zaldivar, B., Bonaventura, N. R., Calore, F., & Iocco, F. (2023). Determining the dark matter distribution in simulated galaxies with deep learning. Mon. Not. Roy. Astron. Soc., 525(4), 6015–6035.
Abstract: We present a novel method of inferring the dark matter (DM) content and spatial distribution within galaxies, using convolutional neural networks (CNNs) trained within state-of-the-art hydrodynamical simulations (Illustris-TNG100). Within the controlled environment of the simulation, the framework we have developed is capable of inferring the DM mass distribution within galaxies of mass similar to 10(11)-10(13)M(circle dot) from the gravitationally baryon-dominated internal regions to the DM-rich, baryon-depleted outskirts of the galaxies, with a mean absolute error always below approximate to 0.25 when using photometrical and spectroscopic information. With respect to traditional methods, the one presented here also possesses the advantages of not relying on a pre-assigned shape for the DM distribution, to be applicable to galaxies not necessarily in isolation, and to perform very well even in the absence of spectroscopic observations.
|
|
|
Athenodorou, A., Binosi, D., Boucaud, P., De Soto, F., Papavassiliou, J., Rodriguez-Quintero, J., et al. (2016). On the zero crossing of the three-gluon vertex. Phys. Lett. B, 761, 444–449.
Abstract: We report on new results on the infrared behavior of the three-gluon vertex in quenched Quantum Chromodynamics, obtained from large-volume lattice simulations. The main focus of our study is the appearance of the characteristic infrared feature known as 'zero crossing', the origin of which is intimately connected with the nonperturbative masslessness of the Faddeev-Popov ghost. The appearance of this effect is clearly visible in one of the two kinematic configurations analyzed, and its theoretical origin is discussed in the framework of Schwinger-Dyson equations. The effective coupling in the momentum subtraction scheme that corresponds to the three-gluon vertex is constructed, revealing the vanishing of the effective interaction at the exact location of the zero crossing.
|
|
|
n_TOF Collaboration(Zugec, P. et al), Domingo-Pardo, C., Giubrone, G., & Tain, J. L. (2014). GEANT4 simulation of the neutron background of the C6D6 set-up for capture studies at n_TOF. Nucl. Instrum. Methods Phys. Res. A, 760, 57–67.
Abstract: The neutron sensitivity of the Cr6D6 detector setup used at nTOF facility for capture measurements has been studied by means of detailed GEANT4 simulations. A realistic software replica of the entire nTOF experimental hall, including the neutron beam line, sample, detector supports and the walls of the experimental area has been implemented in the simulations. The simulations have been analyzed in the same manner as experimental data, in particular by applying the Pulse Height Weighting Technique. The simulations have been validated against a measurement of the neutron background performed with a(nat)-C sample, showing an excellent agreement above 1 keV. At lower energies, an additional component in the measured C-nat yield has been discovered, which prevents the use of C-nat data for neutron background estimates at neutron energies below a few hundred eV. The origin and time structure of the neutron background have been derived from the simulations. Examples of the neutron background for two different samples are demonstrating the important role of accurate simulations of the neutron background in capture cross-section measurements.
|
|
|
AGATA Collaboration(Soderstrom, P. A. et al), & Gadea, A. (2011). Interaction position resolution simulations and in-beam measurements of the AGATA HPGe detectors. Nucl. Instrum. Methods Phys. Res. A, 638(1), 96–109.
Abstract: The interaction position resolution of the segmented HPGe detectors of an AGATA triple cluster detector has been studied through Monte Carlo simulations and in an in-beam experiment. A new method based on measuring the energy resolution of Doppler-corrected gamma-ray spectra at two different target to detector distances is described. This gives the two-dimensional position resolution in the plane perpendicular to the direction of the emitted gamma-ray. The gamma-ray tracking was used to determine the full energy of the gamma-rays and the first interaction point, which is needed for the Doppler correction. Five different heavy-ion induced fusion-evaporation reactions and a reference reaction were selected for the simulations. The results of the simulations show that the method works very well and gives a systematic deviation of <1 mm in the FVVHM of the interaction position resolution for the gamma-ray energy range from 60 keV to 5 MeV. The method was tested with real data from an in-beam measurement using a (30)5i beam at 64 MeV on a thin C-12 target. Pulse-shape analysis of the digitized detector waveforms and gamma-ray tracking was performed to determine the position of the first interaction point, which was used for the Doppler corrections. Results of the dependency of the interaction position resolution on the gamma-ray energy and on the energy, axial location and type of the first interaction point, are presented. The FVVHM of the interaction position resolution varies roughly linearly as a function of gamma-ray energy from 8.5 mm at 250 key to 4 mm at 1.5 MeV, and has an approximately constant value of about 4 mm in the gamma-ray energy range from 1.5 to 4 MeV.
|
|
|
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.
|
|
|
Doncel, M., Cederwall, B., Martin, S., Quintana, B., Gadea, A., Farnea, E., et al. (2015). Conceptual design of a high resolution Ge array with tracking and imaging capabilities for the DESPEC (FAIR) experiment. J. Instrum., 10, P06010–15pp.
Abstract: We present results of Monte Carlo simulations for the conceptual design of the high-resolution DESPEC Germanium Array Spectrometer (DEGAS) proposed for the Facility for Ion and Antiproton Research (FAIR) under construction at Darmstadt, Germany. The project is carried out in three phases, although only results for the two first phases will be addressed in this work. The first phase will consist of a re-arrangement of the EUROBALL cluster detectors previously used in the RISING campaign at GSI. The second phase is based on coupling AGATA-type triple-cluster detectors with EUROBALL cluster detectors in a compact geometry around the active ion implantation target of DESPEC.
|
|
|
de Putter, R., Wagner, C., Mena, O., Verde, L., & Percival, W. J. (2012). Thinking outside the box: effects of modes larger than the survey on matter power spectrum covariance. J. Cosmol. Astropart. Phys., 04(4), 019–31pp.
Abstract: Accurate power spectrum (or correlation function) covariance matrices are a crucial requirement for cosmological parameter estimation from large scale structure surveys. In order to minimize reliance on computationally expensive mock catalogs, it is important to have a solid analytic understanding of the different components that make up a covariance matrix. Considering the matter power spectrum covariance matrix, it has recently been found that there is a potentially dominant effect on mildly non-linear scales due to power in modes of size equal to and larger than the survey volume. This beat coupling effect has been derived analytically in perturbation theory and while it has been tested with simulations, some questions remain unanswered. Moreover, there is an additional effect of these large modes, which has so far not been included in analytic studies, namely the effect on the estimated average density which enters the power spectrum estimate. In this article, we work out analytic, perturbation theory based expressions including both the beat coupling and this local average effect and we show that while, when isolated, beat coupling indeed causes large excess covariance in agreement with the literature, in a realistic scenario this is compensated almost entirely by the local average effect, leaving only similar to 10% of the excess. We test our analytic expressions by comparison to a suite of large N-body simulations, using both full simulation boxes and subboxes thereof to study cases without beat coupling, with beat coupling and with both beat coupling and the local average effect. For the variances, we find excellent agreement with the analytic expressions for k < 0.2 hMpc(-1) at z = 0.5, while the correlation coefficients agree to beyond k = 0.4 hMpc(-1). As expected, the range of agreement increases towards higher redshift and decreases slightly towards z = 0. We finish by including the large-mode effects in a full covariance matrix description for arbitrary survey geometry and confirming its validity using simulations. This may be useful as a stepping stone towards building an actual galaxy (or other tracer's) power spectrum covariance matrix.
|
|
|
Figueroa, D. G., Florio, A., Torrenti, F., & Valkenburg, W. (2023). CosmoLattice: A modern code for lattice simulations of scalar and gauge field dynamics in an expanding universe. Comput. Phys. Commun., 283, 108586–13pp.
Abstract: This paper describes CosmoGattice, a modern package for lattice simulations of the dynamics of interacting scalar and gauge fields in an expanding universe. CosmoGattice incorporates a series of features that makes it very versatile and powerful: i) it is written in C++ fully exploiting the object oriented programming paradigm, with a modular structure and a clear separation between the physics and the technical details, ii) it is MPI-based and uses a discrete Fourier transform parallelized in multiple spatial dimensions, which makes it specially appropriate for probing scenarios with well -separated scales, running very high resolution simulations, or simply very long ones, iii) it introduces its own symbolic language, defining field variables and operations over them, so that one can introduce differential equations and operators in a manner as close as possible to the continuum, iv) it includes a library of numerical algorithms, ranging from O(delta t(2)) to O(delta t(10)) methods, suitable for simulating global and gauge theories in an expanding grid, including the case of 'self-consistent' expansion sourced by the fields themselves. Relevant observables are provided for each algorithm (e.g. energy densities, field spectra, lattice snapshots) and we note that, remarkably, all our algorithms for gauge theories (Abelian or non-Abelian) always respect the Gauss constraint to machine precision. Program summary Program Title:: CosmoGattice CPC Library link to program files: https://doi .org /10 .17632 /44vr5xssc6 .1 Developer's repository link: http://github .com /cosmolattice /cosmolattice Licensing provisions: MIT Programming language: C++, MPI Nature of problem: The phenomenology of high energy physics in the early universe is typically characterized by non-linear dynamics, which cannot be captured accurately with analytical techniques. In order to fully understand the non-linearities developed in a given scenario, one needs to carry out lattice simulations. A number of public packages for lattice simulations have appeared over the years, but most of them are only capable of simulating scalar fields. However, realistic models of particle physics do contain other kind of field species, such as (Abelian or non-Abelian) gauge fields, whose non-linear dynamics can also play a relevant role in the early universe. Tensor modes representing gravitational waves are also naturally expected in many scenarios. Solution method: CosmoGattice represents a modern code for lattice simulations of scalar-gauge field theories in an expanding universe. It allows for the simulation of the evolution of interacting (singlet) scalar fields, charged scalar fields under U(1) and/or SU(2) gauge groups, and the corresponding associated Abelian and/or non-Abelian gauge fields. From version 1.1 onward, CosmoGattice also allows to simulate the production of gravitational waves. Simulations can be done either in a flat space-time background, or in a homogeneous and isotropic (spatially flat) expanding FLRW background. CosmoGattice provides symplectic integrators, with accuracy ranging from O (delta t(2)) up to O(delta t(10)), to simuate the non-linear dynamics of the appropriate fields in comoving three-dimensional lattices. The code is parallelized with MPI, and uses a discrete Fourier Transform parallelized in multiple spatial dimensions, which makes it a very powerful code for probing physical problems with well-separated scales. Moreover, the code has been designed as a `platform' to implement any system of dynamical equations suitable for discretization on a lattice.
|
|
|
DUNE Collaboration(Abud, A. A. et al), Amedo, P., Antonova, M., Barenboim, G., Cervera-Villanueva, A., De Romeri, V., et al. (2023). Highly-parallelized simulation of a pixelated LArTPC on a GPU. J. Instrum., 18(4), P04034–35pp.
Abstract: The rapid development of general-purpose computing on graphics processing units (GPGPU) is allowing the implementation of highly-parallelized Monte Carlo simulation chains for particle physics experiments. This technique is particularly suitable for the simulation of a pixelated charge readout for time projection chambers, given the large number of channels that this technology employs. Here we present the first implementation of a full microphysical simulator of a liquid argon time projection chamber (LArTPC) equipped with light readout and pixelated charge readout, developed for the DUNE Near Detector. The software is implemented with an end-to-end set of GPU-optimized algorithms. The algorithms have been written in Python and translated into CUDA kernels using Numba, a just-in-time compiler for a subset of Python and NumPy instructions. The GPU implementation achieves a speed up of four orders of magnitude compared with the equivalent CPU version. The simulation of the current induced on 103 pixels takes around 1 ms on the GPU, compared with approximately 10 s on the CPU. The results of the simulation are compared against data from a pixel-readout LArTPC prototype.
|
|