|
Pallis, C. (2014). Linking Starobinsky-type inflation in no-scale supergravity to MSSM. J. Cosmol. Astropart. Phys., 04(4), 024–31pp.
Abstract: A novel realization of the Starobinsky inflationary model within a moderate extension of the Minimal Supersymmetric Standard Model (MSSM) is presented. The proposed superpotential is uniquely determined by applying a continuous R and a Z2 discrete symmetry, whereas the Kahler potential is associated with a no-scale-type SU(54, 1)/ SU(54) x U(1) R X Z2 Kahler manifold. The inflaton is identified with a Higgs-like modulus whose the vacuum expectation value controls the gravitational strength. Thanks to a strong enough coupling (with a parameter CT involved) between the inflaton and the Ricci scalar curvature, inflation can be attained even for subplanckian values of the inflaton with CT >= 76 and the corresponding effective theory being valid up to the Planck scale. The inflationary observables turn out to be in agreement with the current data and the inflaton mass is predicted to be 3 10(3) GeV. At the cost of a relatively small superpotential coupling constant, the model offers also a resolution of the f,t problem of MSSM for CT <= 4500 and gravitino heavier than about 10(4) GeV. Supplementing MSSM by three right-handed neutrinos we show that spontaneously arising couplings between the inflaton and the particle content of MSSM not only ensure a sufficiently low reheating temperature but also support a scenario of non-thermal leptogenesis consistently with the neutrino oscillation parameters.
|
|
|
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.
|
|
|
Oliver, S., Gimenez-Alventosa, V., Berumen, F., Gimenez, V., Beaulieu, L., Ballester, F., et al. (2023). Benchmark of the PenRed Monte Carlo framework for HDR brachytherapy. Z. Med. Phys., 33(4), 511–528.
Abstract: Purpose: The purpose of this study is to validate the PenRed Monte Carlo framework for clinical applications in brachytherapy. PenRed is a C++ version of Penelope Monte Carlo code with additional tallies and utilities. Methods and materials: Six benchmarking scenarios are explored to validate the use of PenRed and its improved bachytherapy-oriented capabilities for HDR brachytherapy. A new tally allowing the evaluation of collisional kerma for any material using the track length kerma estimator and the possibility to obtain the seed positions, weights and directions processing directly the DICOM file are now implemented in the PenRed distribution. The four non-clinical test cases developed by the Joint AAPM-ESTRO-ABG-ABS WG-DCAB were evaluated by comparing local and global absorbed dose differences with respect to established reference datasets. A prostate and a palliative lung cases, were also studied. For them, absorbed dose ratios, global absorbed dose differences, and cumulative dose-volume histograms were obtained and discussed. Results: The air-kerma strength and the dose rate constant corresponding to the two sources agree with the reference datatests within 0.3% (Sk) and 0.1% (K). With respect to the first three WG-DCAB test cases, more than 99.8% of the voxels present local (global) differences within +/- 1%(+/- 0.1%) of the reference datasets. For test Case 4 reference dataset, more than 94.9%(97.5%) of voxels show an agreement within +/- 1%(+/- 0.1%), better than similar benchmarking calculations in the literature. The track length kerma estimator scorer implemented increases the numerical efficiency of brachytherapy calculations two orders of magnitude, while the specific brachytherapy source allows the user to avoid the use of error-prone intermediate steps to translate the DICOM information into the simulation. In both clinical cases, only minor absorbed dose differences arise in the low-dose isodoses. 99.8% and 100% of the voxels have a global absorbed dose difference ratio within +/- 0.2%for the prostate and lung cases, respectively. The role played by the different segmentation and composition material in the bone structures was discussed, obtaining negligible absorbed dose differ-ences. Dose-volume histograms were in agreement with the reference data.Conclusions: PenRed incorporates new tallies and utilities and has been validated for its use for detailed and precise high-dose-rate brachytherapy simulations.
|
|
|
Oliveira, C. A. B., Sorel, M., Martin-Albo, J., Gomez-Cadenas, J. J., Ferreira, A. L., & Veloso, J. F. C. A. (2011). Energy resolution studies for NEXT. J. Instrum., 6, P05007–13pp.
Abstract: This work aims to present the current state of simulations of electroluminescence (EL) produced in gas-based detectors with special interest for NEXT – Neutrino Experiment with a Xenon TPC. NEXT is a neutrinoless double beta decay experiment, thus needs outstanding energy resolution which can be achieved by using electroluminescence. The process of light production is reviewed and properties such as EL yield and associated fluctuations, excitation and electroluminescence efficiencies, and energy resolution, are calculated. An EL production region with a 5 mm width gap between two infinite parallel planes is considered, where a uniform electric field is produced. The pressure and temperature considered are 10 bar and 293 K, respectively. The results show that, even for low values of VUV photon detection efficiency, good energy resolution can be achieved: below 0.4% (FWHM) at Q(beta beta) = 2.458 MeV.
|
|
|
Mosbech, M. R., Boehm, C., Hannestad, S., Mena, O., Stadler, J., & Wong, Y. Y. Y. (2021). The full Boltzmann hierarchy for dark matter-massive neutrino interactions. J. Cosmol. Astropart. Phys., 03(3), 066–31pp.
Abstract: The impact of dark matter-neutrino interactions on the measurement of the cosmological parameters has been investigated in the past in the context of massless neutrinos exclusively. Here we revisit the role of a neutrino-dark matter coupling in light of ongoing cosmological tensions by implementing the full Boltzmann hierarchy for three massive neutrinos. Our tightest 95% CL upper limit on the strength of the interactions, parameterized via u(chi) = sigma(0)/sigma(Th) (m(chi)/100GeV)(-1), is u(chi) <= 3.34 . 10(-4), arising from a combination of Planck TTTEEE data, Planck lensing data and SDSS BAO data. This upper bound is, as expected, slightly higher than previous results for interacting massless neutrinos, due to the correction factor associated with neutrino masses. We find that these interactions significantly relax the lower bounds on the value of sigma 8 that is inferred in the context of Lambda CDM from the Planck data, leading to agreement within 1-2 sigma with weak lensing estimates of sigma 8, as those from KiDS1000. However, the presence of these interactions barely affects the value of the Hubble constant H-0.
|
|