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.

Abstract: We study the energy distribution and equation of state of the universe between the end of inflation and the onset of radiation domination (RD), considering observationally consistent single-field inflationary scenarios, with a potential 'flattening' at large field values, and a monomial shape V(phi) proportional to vertical bar phi vertical bar(p) around the origin. As a proxy for (p)reheating, we include a quadratic interaction g(2)phi X-2(2) between the inflaton phi and a light scalar 'daughter' field X, with g(2) > 0. We capture the non-perturbative and non-linear nature of the system dynamics with lattice simulations, obtaining that: i) the final energy transferred to X depends only on p, not on g(2); ii) the final transfer of energy is always negligible for 2 <= p < 4, and of order similar to 50% for p >= 4; iii) the system goes at late times to matter-domination for p = 2, and always to RD for p > 2. In the latter case we calculate the number of e-folds until RD, significantly reducing the uncertainty in the inflationary observables Tl-s and r.

Abstract: We calculate neutrino emissivities from self-annihilating dark matter (DM) (chi) in the dense and hot stellar interior of a (proto)neutron star. Using a model where DM interacts with nucleons in the stellar core through a pseudoscalar boson (a) we find that the neutrino production rates from the dominant reaction channels chi -> nu(nu) over bar or chi chi -> aa, with subsequent decay of the mediator a -> nu(nu) over bar, could locally match and even surpass those of the standard neutrinos from the modified nuclear URCA processes at early ages. We find that the emitting region can be localized in a tiny fraction of the star (less than a few percent of the core volume) and the process can last its entire lifetime for some cases under study. We discuss the possible consequences of our results for stellar cooling in light of existing DM constraints.

Abstract: In this paper, we explore the impact of flavor violating neutral current non-standard interaction (NSI) parameter epsilon(mu tau) in the oscillation of atmospheric neutrinos and antineutrinos separately using the 50 kt magnetized ICAL detector at INO. We find that due to non-zero epsilon(mu tau), nu(mu) -> nu(mu) and (nu) over bar (mu) -> (nu) over bar (mu) transition probabilities get modified substantially at higher energies and longer baselines, where vacuum oscillation dominates. We demonstrate for the first time that by adding the hadron energy information along with the muon energy and muon direction in each event, the sensitivity of ICAL to the NSI parameter epsilon(mu tau) can be enhanced significantly. The most optimistic bound on epsilon(mu tau) that we obtain is – 0.01 < epsilon(mu tau) < 0.01 at 90% C.L. using 500 kt.yr exposure and considering E-mu, cos theta(mu), and E-had' as observables in their ranges of [1, 21] GeV, [- 1, 1], and [0, 25] GeV, respectively. We discuss for the first time the importance of the charge identification capability of the ICAL detector to have better constraints on epsilon(mu t). We also study the impact of non-zero epsilon(mu tau) on mass hierarchy determination and precision measurement of oscillation parameters.

Abstract: We examine the relationship between three approaches (Hadamard, DeWitt-Schwinger, and adiabatic) to the renormalization of expectation values of field operators acting on a charged quantum scalar field. First, we demonstrate that the DeWitt-Schwinger representation of the Feynman Green's function is a particular case of the Hadamard representation. Next, we restrict attention to a spatially flat Friedmann-Lemaitre-Robertson-Walker universe with time-dependent, purely electric, background electromagnetic field, considering two-, three-, and four-dimensional space-times. Working to the order required for the renormalization of the stress-energy tensor, we find the adiabatic and DeWitt-Schwinger expansions of the Green's function when the space-time points are spatially separated. In two and four dimensions, the resulting DeWitt-Schwinger and adiabatic expansions are identical. In three dimensions, the DeWittSchwinger expansion contains terms of adiabatic order 4 that are not necessary for the renormalization of the stress-energy tensor and hence absent in the adiabatic expansion. The equivalence of the DeWittSchwinger and adiabatic approaches to renormalization in the scenario considered is thereby demonstrated in even dimensions. In odd dimensions the situation is less clear and further investigation is required in order to determine whether adiabatic renormalization is a locally covariant renormalization prescription.