|
|
Oliver, S., Rodriguez Bosca, S., & Gimenez-Alventosa, V. (2024). Enabling particle transport on CAD-based geometries for radiation simulations with penRed. Comput. Phys. Commun., 298, 109091–11pp.
Abstract: Geometry construction is a fundamental aspect of any radiation transport simulation, regardless of the Monte Carlo code being used. Typically, this process is tedious, time-consuming, and error-prone. The conventional approach involves defining geometries using mathematical objects or surfaces. However, this method comes with several limitations, especially when dealing with complex models, particularly those with organic shapes. Furthermore, since each code employs its own format and methodology for defining geometries, sharing and reproducing simulations among researchers becomes a challenging task. Consequently, many codes have implemented support for simulating over geometries constructed via Computer-Aided Design (CAD) tools. Unfortunately, this feature is lacking in penRed and other PENELOPE physics-based codes. Therefore, the objective of this work is to implement such support within the penRed framework. New version program summary Program Title: Parallel Engine for Radiation Energy Deposition (penRed) CPC Library link to program files: https://doi.org/10.17632/rkw6tvtngy.2 Developer's repository link: https://github.com/PenRed/PenRed Code Ocean capsule: https://codeocean.com/capsule/1041417/tree Licensing provisions: GNU Affero General Public License v3 Programming language: C++ standard 2011. Journal reference of previous version: V. Gimenez-Alventosa, V. Gimenez Gomez, S. Oliver, PenRed: An extensible and parallel Monte-Carlo framework for radiation transport based on PENELOPE, Computer Physics Communications 267 (2021) 108065. doi:https://doi.org/10.1016/j.cpc.2021.108065. Does the new version supersede the previous version?: Yes Reasons for the new version: Implements the capability to simulate on CAD constructed geometries, among many other features and fixes. Summary of revisions: All changes applied through the code versions are summarized in the file CHANGELOG.md in the repository package. Nature of problem: While Monte Carlo codes have proven valuable in simulating complex radiation scenarios, they rely heavily on accurate geometrical representations. In the same way as many other Monte Carlo codes, penRed employs simple geometric quadric surfaces like planes, spheres and cylinders to define geometries. However, since these geometric models offer a certain level of flexibility, these representations have limitations when it comes to simulating highly intricate and irregular shapes. Anatomic structures, for example, require detailed representations of organs, tissues and bones, which are difficult to achieve using basic geometric objects. Similarly, complex devices or intricate mechanical systems may have designs that cannot be accurately represented within the constraints of such geometric models. Moreover, when the complexity of the model increases, geometry construction process becomes more difficult, tedious, time-consuming and error-prone [2]. Also, as each Monte Carlo geometry library uses its own format and construction method, reproducing the same geometry among different codes is a challenging task. Solution method: To face the problems stated above, the objective of this work is to implement the capability to simulate using irregular and adaptable meshed geometries in the penRed framework. This kind of meshes can be constructed using Computer-Aided Design (CAD) tools, the use of which is very widespread and streamline the design process. This feature has been implemented in a new geometry module named “MESH_BODY” specific for this kind of geometries. This one is freely available and usable within the official penRed package1. It can be used since penRed version 1.9.3b and above.
|
|
|
|
Gimenez-Alventosa, V., Gimenez, V., & Oliver, S. (2021). PenRed: An extensible and parallel Monte-Carlo framework for radiation transport based on PENELOPE. Comput. Phys. Commun., 267, 108065–12pp.
Abstract: Monte Carlo methods provide detailed and accurate results for radiation transport simulations. Unfortunately, the high computational cost of these methods limits its usage in real-time applications. Moreover, existing computer codes do not provide a methodology for adapting these kinds of simulations to specific problems without advanced knowledge of the corresponding code system, and this restricts their applicability. To help solve these current limitations, we present PenRed, a general-purpose, standalone, extensible and modular framework code based on PENELOPE for parallel Monte Carlo simulations of electron-photon transport through matter. It has been implemented in C++ programming language and takes advantage of modern object-oriented technologies. In addition, PenRed offers the capability to read and process DICOM images as well as to construct and simulate image-based voxelized geometries, so as to facilitate its usage in medical applications. Our framework has been successfully verified against the original PENELOPE Fortran code. Furthermore, the implemented parallelism has been tested showing a significant improvement in the simulation time without any loss in precision of results. Program summary Program title: PenRed: Parallel Engine for Radiation Energy Deposition. CPC Library link to program files: https://doi .org /10 .17632/rkw6tvtngy.1 Licensing provision: GNU Affero General Public License (AGPL). Programming language: C++ standard 2011. Nature of problem: Monte Carlo simulations usually require a huge amount of computation time to achieve low statistical uncertainties. In addition, many applications necessitate particular characteristics or the extraction of specific quantities from the simulation. However, most available Monte Carlo codes do not provide an efficient parallel and truly modular structure which allows users to easily customise their code to suit their needs without an in-depth knowledge of the code system. Solution method: PenRed is a fully parallel, modular and customizable framework for Monte Carlo simulations of the passage of radiation through matter. It is based on the PENELOPE [1] code system, from which inherits its unique physics models and tracking algorithms for charged particles. PenRed has been coded in C++ following an object-oriented programming paradigm restricted to the C++11 standard. Our engine implements parallelism via a double approach: on the one hand, by using standard C++ threads for shared memory, improving the access and usage of the memory, and, on the other hand, via the MPI standard for distributed memory infrastructures. Notice that both kinds of parallelism can be combined together in the same simulation. Moreover, both threads and MPI processes, can be balanced using the builtin load balance system (RUPER-LB [30]) to maximise the performance on heterogeneous infrastructures. In addition, PenRed provides a modular structure with methods designed to easily extend its functionality. Thus, users can create their own independent modules to adapt our engine to their needs without changing the original modules. Furthermore, user extensions will take advantage of the builtin parallelism without any extra effort or knowledge of parallel programming. Additional comments including restrictions and unusual features: PenRed has been compiled in linux systems withg++ of GCC versions 4.8.5, 7.3.1, 8.3.1 and 9; clang version 3.4.2 and intel C++ compiler (icc) version 19.0.5.281. Since it is a C++11-standard compliant code, PenRed should be able to compile with any compiler with C++11 support. In addition, if the code is compiled without MPI support, it does not require any non standard library. To enable MPI capabilities, the user needs to install whatever available MPI implementation, such as openMPI [24] or mpich [25], which can be found in the repositories of any linux distribution. Finally, to provide DICOM processing support, PenRed can be optionally compiled using the dicom toolkit (dcmtk) [32] library. Thus, PenRed has only two optional dependencies, an MPI implementation and the dcmtk library.
|
|
|
|
ANTARES Collaboration(Bhandari, S. et al), Barrios-Marti, J., Coleiro, A., Hernandez-Rey, J. J., Illuminati, G., Tönnis, C., et al. (2018). The SUrvey for Pulsars and Extragalactic Radio Bursts – II. New FRB discoveries and their follow-up. Mon. Not. Roy. Astron. Soc., 475(2), 1427–1446.
Abstract: We report the discovery of four Fast Radio Bursts (FRBs) in the ongoing SUrvey for Pulsars and Extragalactic Radio Bursts at the Parkes Radio Telescope: FRBs 150610, 151206, 151230 and 160102. Our real-time discoveries have enabled us to conduct extensive, rapid multimessenger follow-up at 12 major facilities sensitive to radio, optical, X-ray, gamma-ray photons and neutrinos on time-scales ranging from an hour to a few months post-burst. No counterparts to the FRBs were found and we provide upper limits on afterglow luminosities. None of the FRBs were seen to repeat. Formal fits to all FRBs show hints of scattering while their intrinsic widths are unresolved in time. FRB 151206 is at low Galactic latitude, FRB 151230 shows a sharp spectral cut-off, and FRB 160102 has the highest dispersion measure (DM = 2596.1 +/- 0.3 pc cm(-3)) detected to date. Three of the FRBs have high dispersion measures (DM > 1500 pc cm(-3)), favouring a scenario where the DMis dominated by contributions from the intergalactic medium. The slope of the Parkes FRB source counts distribution with fluences > 2 Jy ms is alpha = – 2.2(-1.2)(+0.6) and still consistent with a Euclidean distribution (alpha = -3/2). We also find that the all-sky rate is 1.7(-0.9)(+1.5) x 10(3)FRBs/(4 pi sr)/day above similar to 2 Jy ms and there is currently no strong evidence for a latitude- dependent FRB sky rate.
|
|
|
|
ATLAS Collaboration(Aad, G. et al), Aparisi Pozo, J. A., Bailey, A. J., Cabrera Urban, S., Cardillo, F., Castillo, F. L., et al. (2021). Measurements of sensor radiation damage in the ATLAS inner detector using leakage currents. J. Instrum., 16(8), P08025–46pp.
Abstract: Non-ionizing energy loss causes bulk damage to the silicon sensors of the ATLAS pixel and strip detectors. This damage has important implications for data-taking operations, charged-particle track reconstruction, detector simulations, and physics analysis. This paper presents simulations and measurements of the leakage current in the ATLAS pixel detector and semiconductor tracker as a function of location in the detector and time, using data collected in Run 1 (2010-2012) and Run 2 (2015-2018) of the Large Hadron Collider. The extracted fluence shows a much stronger vertical bar z vertical bar-dependence in the innermost layers than is seen in simulation. Furthermore, the overall fluence on the second innermost layer is significantly higher than in simulation, with better agreement in layers at higher radii. These measurements are important for validating the simulation models and can be used in part to justify safety factors for future detector designs and interventions.
|
|
|
|
NEXT Collaboration(Alvarez, V. et al), Carcel, S., Cervera-Villanueva, A., Diaz, J., Ferrario, P., Gil, A., et al. (2013). Radiopurity control in the NEXT-100 double beta decay experiment: procedures and initial measurements. J. Instrum., 8, T01002–19pp.
Abstract: The “Neutrino Experiment with a Xenon Time-Projection Chamber” (NEXT) is intended to investigate the neutrinoless double beta decay of Xe-136, which requires a severe suppression of potential backgrounds. An extensive screening and material selection process is underway for NEXT since the control of the radiopurity levels of the materials to be used in the experimental set-up is a must for rare event searches. First measurements based on Glow Discharge Mass Spectrometry and gamma-ray spectroscopy using ultra-low background germanium detectors at the Laboratorio Subterraneo de Canfranc (Spain) are described here. Activity results for natural radioactive chains and other common radionuclides are summarized, being the values obtained for some materials like copper and stainless steel very competitive. The implications of these results for the NEXT experiment are also discussed.
|
|
|
|
NEXT Collaboration(Cebrian, S. et al), Alvarez, V., Carcel, S., Cervera-Villanueva, A., Diaz, J., Ferrario, P., et al. (2015). Radiopurity assessment of the tracking readout for the NEXT double beta decay experiment. J. Instrum., 10, P05006–16pp.
Abstract: The “Neutrino Experiment with a Xenon Time-Projection Chamber” (NEXT) is intended to investigate the neutrinoless double beta decay of Xe-136, which requires a severe suppression of potential backgrounds; therefore, an extensive screening and selection process is underway to control the radiopurity levels of the materials to be used in the experimental set-up of NEXT. The detector design combines the measurement of the topological signature of the event for background discrimination with the energy resolution optimization. Separate energy and tracking readout planes are based on different sensors: photomultiplier tubes for calorimetry and silicon multi-pixel photon counters for tracking. The design of a radiopure tracking plane, in direct contact with the gas detector medium, was specially challenging since the needed components like printed circuit boards, connectors, sensors or capacitors have typically, according to available information in databases and in the literature, activities too large for experiments requiring ultra-low background conditions. Here, the radiopurity assessment of tracking readout components based on gamma-ray spectroscopy using ultra-low background germanium detectors at the Laboratorio Subterraneo de Canfranc (Spain) is described. According to the obtained results, radiopure enough printed circuit boards made of kapton and copper, silicon photomultipliers and other required components, fulfilling the requirement of an overall background level in the region of interest of at most 8 x 10(-4) counts keV(-1) kg(-1) y(-1), have been identified.
|
|
|
|
Angles-Castillo, A., Perez, A., & Roldan, E. (2024). Bright and dark solitons in a photonic nonlinear quantum walk: lessons from the continuum. New J. Phys., 26(2), 023004–16pp.
Abstract: We propose a nonlinear quantum walk model inspired in a photonic implementation in which the polarization state of the light field plays the role of the coin-qubit. In particular, we take profit of the nonlinear polarization rotation occurring in optical media with Kerr nonlinearity, which allows to implement a nonlinear coin operator, one that depends on the state of the coin-qubit. We consider the space-time continuum limit of the evolution equation, which takes the form of a nonlinear Dirac equation. The analysis of this continuum limit allows us to gain some insight into the existence of different solitonic structures, such as bright and dark solitons. We illustrate several properties of these solitons with numerical calculations, including the effect on them of an additional phase simulating an external electric field.
|
|
|
|
Jay, G., Arnault, P., & Debbasch, F. (2021). Dirac quantum walks with conserved angular momentum. Quantum Stud. Math. Found., 8, 419–430.
Abstract: A quantum walk (QW) simulating the flat (1+2)D Dirac equation on a spatial polar grid is constructed. Because fermions are represented by spinors, which do not constitute a representation of the rotation group SO(3), but rather of its double cover SU(2), the QW can only be defined globally on an extended spacetime where the polar angle extends from 0 to 4 pi. The coupling of the QW with arbitrary electromagnetic fields is also presented. Finally, the cylindrical relativistic Landau levels of the Dirac equation are computed explicitly and simulated by the QW.
|
|
|
|
Nzongani, U., Zylberman, J., Doncecchi, C. E., Perez, A., Debbasch, F., & Arnault, P. (2023). Quantum circuits for discrete-time quantum walks with position-dependent coin operator. Quantum Inf. Process., 22(7), 270–46pp.
Abstract: The aim of this paper is to build quantum circuits that implement discrete-time quantum walks having an arbitrary position-dependent coin operator. The position of the walker is encoded in base 2: with n wires, each corresponding to one qubit, we encode 2(n) position states. The data necessary to define an arbitrary position-dependent coin operator is therefore exponential in n. Hence, the exponentiality will necessarily appear somewhere in our circuits. We first propose a circuit implementing the position-dependent coin operator, that is naive, in the sense that it has exponential depth and implements sequentially all appropriate position-dependent coin operators. We then propose a circuit that “transfers” all the depth into ancillae, yielding a final depth that is linear in n at the cost of an exponential number of ancillae. Themain idea of this linear-depth circuit is to implement in parallel all coin operators at the different positions. Reducing the depth exponentially at the cost of having an exponential number of ancillae is a goal which has already been achieved for the problem of loading classical data on a quantum circuit (Araujo in Sci Rep 11:6329, 2021) (notice that such a circuit can be used to load the initial state of the walker). Here, we achieve this goal for the problem of applying a position-dependent coin operator in a discrete-time quantum walk. Finally, we extend the result of Welch (New J Phys 16:033040, 2014) from position-dependent unitaries which are diagonal in the position basis to position-dependent 2 x 2-block-diagonal unitaries: indeed, we show that for a position dependence of the coin operator (the block-diagonal unitary) which is smooth enough, one can find an efficient quantum-circuit implementation approximating the coin operator up to an error epsilon (in terms of the spectral norm), the depth and size of which scale as O(1/epsilon). A typical application of the efficient implementation would be the quantum simulation of a relativistic spin-1/2 particle on a lattice, coupled to a smooth external gauge field; notice that recently, quantum spatial-search schemes have been developed which use gauge fields as the oracle, to mark the vertex to be found (Zylberman in Entropy 23:1441, 2021), (Fredon arXiv:2210.13920). A typical application of the linear-depth circuit would be when there is spatial noise on the coin operator (and hence a non-smooth dependence in the position).
|
|
|
|
Hinarejos, M., Bañuls, M. C., & Perez, A. (2013). A Study of Wigner Functions for Discrete-Time Quantum Walks. J. Comput. Theor. Nanosci., 10(7), 1626–1633.
Abstract: We perform a systematic study of the discrete time Quantum Walk on one dimension using Wigner functions, which are generalized to include the chirality (or coin) degree of freedom. In particular, we analyze the evolution of the negative volume in phase space, as a function of time, for different initial states. This negativity can be used to quantify the degree of departure of the system from a classical state. We also relate this quantity to the entanglement between the coin and walker subspaces.
|
|
