|
Balazs, C. et al, Mamuzic, J., & Ruiz de Austri, R. (2021). A comparison of optimisation algorithms for high-dimensional particle and astrophysics applications. J. High Energy Phys., 05(5), 108–46pp.
Abstract: Optimisation problems are ubiquitous in particle and astrophysics, and involve locating the optimum of a complicated function of many parameters that may be computationally expensive to evaluate. We describe a number of global optimisation algorithms that are not yet widely used in particle astrophysics, benchmark them against random sampling and existing techniques, and perform a detailed comparison of their performance on a range of test functions. These include four analytic test functions of varying dimensionality, and a realistic example derived from a recent global fit of weak-scale supersymmetry. Although the best algorithm to use depends on the function being investigated, we are able to present general conclusions about the relative merits of random sampling, Differential Evolution, Particle Swarm Optimisation, the Covariance Matrix Adaptation Evolution Strategy, Bayesian Optimisation, Grey Wolf Optimisation, and the PyGMO Artificial Bee Colony, Gaussian Particle Filter and Adaptive Memory Programming for Global Optimisation algorithms.
|
|
|
Balbinot, R., & Fabbri, A. (2023). The Hawking Effect in the Particles-Partners Correlations. Physics, 5(4), 968–982.
Abstract: We analyze the correlations functions across the horizon in Hawking black hole radiation to reveal the correlations between Hawking particles and their partners. The effects of the underlying space-time on this are shown in various examples ranging from acoustic black holes to regular black holes.
|
|
|
Balibrea-Correa, J., Lerendegui-Marco, J., Babiano-Suarez, V., Caballero, L., Calvo, D., Ladarescu, I., et al. (2021). Machine Learning aided 3D-position reconstruction in large LaCl3 crystals. Nucl. Instrum. Methods Phys. Res. A, 1001, 165249–17pp.
Abstract: We investigate five different models to reconstruct the 3D gamma-ray hit coordinates in five large LaCl3(Ce) monolithic crystals optically coupled to pixelated silicon photomultipliers. These scintillators have a base surface of 50 x 50 mm(2) and five different thicknesses, from 10 mm to 30 mm. Four of these models are analytical prescriptions and one is based on a Convolutional Neural Network. Average resolutions close to 1-2 mm fwhm are obtained in the transverse crystal plane for crystal thicknesses between 10 mm and 20 mm using analytical models. For thicker crystals average resolutions of about 3-5 mm fwhm are obtained. Depth of interaction resolutions between 1 mm and 4 mm are achieved depending on the distance of the interaction point to the photosensor surface. We propose a Machine Learning algorithm to correct for linearity distortions and pin-cushion effects. The latter allows one to keep a large field of view of about 70%-80% of the crystal surface, regardless of crystal thickness. This work is aimed at optimizing the performance of the so-called Total Energy Detector with Compton imaging capability (i-TED) for time-of-flight neutron capture cross-section measurements.
|
|
|
Ballester, F., Tedgren, A. C., Granero, D., Haworth, A., Mourtada, F., Fonseca, G. P., et al. (2015). A generic high-dose rate Ir-192 brachytherapy source for evaluation of model-based dose calculations beyond the TG-43 formalism. Med. Phys., 42(6), 3048–3062.
Abstract: Purpose: In order to facilitate a smooth transition for brachytherapy dose calculations from the American Association of Physicists in Medicine (AAPM) Task Group No. 43 (TG-43) formalism to model-based dose calculation algorithms (MBDCAs), treatment planning systems (TPSs) using a MBDCA require a set of well-defined test case plans characterized by Monte Carlo (MC) methods. This also permits direct dose comparison to TG-43 reference data. Such test case plans should be made available for use in the software commissioning process performed by clinical end users. To this end, a hypothetical, generic high-dose rate (HDR) Ir-192 source and a virtual water phantom were designed, which can be imported into a TPS. Methods: A hypothetical, generic HDR Ir-192 source was designed based on commercially available sources as well as a virtual, cubic water phantom that can be imported into any TPS in DICOM format. The dose distribution of the generic Ir-192 source when placed at the center of the cubic phantom, and away from the center under altered scatter conditions, was evaluated using two commercial MBDCAs [Oncentra (R) Brachy with advanced collapsed-cone engine (ACE) and BrachyVision AcuRos (TM)]. Dose comparisons were performed using state-of-the-art MC codes for radiation transport, including ALGEBRA, BrachyDose, GEANT4, MCNP5, MCNP6, and pENELopE2008. The methodologies adhered to recommendations in the AAPM TG-229 report on high-energy brachytherapy source dosimetry. TG-43 dosimetry parameters, an along-away dose-rate table, and primary and scatter separated (PSS) data were obtained. The virtual water phantom of (201)(3) voxels (1 mm sides) was used to evaluate the calculated dose distributions. Two test case plans involving a single position of the generic HDR Ir-192 source in this phantom were prepared: (i) source centered in the phantom and (ii) source displaced 7 cm laterally from the center. Datasets were independently produced by different investigators. MC results were then compared against dose calculated using TG-43 and MBDCA methods. Results: TG-43 and PSS datasets were generated for the generic source, the PSS data for use with the ACE algorithm. The dose-rate constant values obtained from seven MC simulations, performed independently using different codes, were in excellent agreement, yielding an average of 1.1109 +/- 0.0004 cGy/(h U) (k = 1, Type A uncertainty). MC calculated dose-rate distributions for the two plans were also found to be in excellent agreement, with differences within type A uncertainties. Differences between commercial MBDCA and MC results were test, position, and calculation parameter dependent. On average, however, these differences were within 1% for ACUROS and 2% for ACE at clinically relevant distances. Conclusions: A hypothetical, generic HDR Ir-192 source was designed and implemented in two commercially available TPSs employing different MBDCAs. Reference dose distributions for this source were benchmarked and used for the evaluation of MBDCA calculations employing a virtual, cubic water phantom in the form of a CT DICOM image series. The implementation of a generic source of identical design in all TPSs using MBDCAs is an important step toward supporting univocal commissioning procedures and direct comparisons between TPSs.
|
|
|
Ballester, F., Granero, D., Perez-Calatayud, J., Venselaar, J. L. M., & Rivard, M. J. (2010). Study of encapsulated Tm-170 sources for their potential use in brachytherapy. Med. Phys., 37(4), 1629–1637.
Abstract: Methods: The authors have assumed a theoretical Tm-170 cylindrical source encapsulated with stainless steel and typical dimensions taken from the currently available HDR Ir-192 brachytherapy sources. The dose-rate distribution was calculated for this source using the GEANT4 Monte Carlo (MC) code considering both photon and electron Tm-170 spectra. The AAPM TG-43 U1 brachytherapy dosimetry parameters were derived. To study general properties of Tm-170 encapsulated sources, spherical sources encapsulated with stainless steel and platinum were also studied. Moreover, the influence of small variations in the active core and capsule dimensions on the dosimetric characteristics was assessed. Treatment times required for a Tm-170 source were compared to those for Ir-192 and Yb-169 for the same contained activity. Results: Due to the energetic beta spectrum and the large electron yield, the bremsstrahlung contribution to the dose was of the same order of magnitude as from the emitted gammas and characteristic x rays. Moreover, the electron spectrum contribution to the dose was significant up to 4 mm from the source center compared to the photon contribution. The dose-rate constant Lambda of the cylindrical source was 1.23 cGy h(-1) U-1. The behavior of the radial dose function showed promise for applications in brachytherapy. Due to the electron spectrum, the anisotropy was large for r < 6 mm. Variations in manufacturing tolerances did not significantly influence the final dosimetry data when expressed in cGy h(-1) U-1. For typical capsule dimensions, maximum reference dose rates of about 0.2, 10, and 2 Gy min(-1) would then be obtained for Tm-170, Ir-192, and Yb-169, respectively, resulting in treatment times greater than those for HDR Ir-192 brachytherapy. Conclusions: The dosimetric characteristics of source designs exploiting the low photon energy of Tm-170 were studied for potential application in HDR-brachytherapy. Dose-rate distributions were obtained for cylindrical and simplified spherical Tm-170 source designs (stainless steel and platinum capsule materials) using MC calculations. Despite the high activity of Tm-170, calculated treatment times were much longer than for Ir-192.
|
|
|
Bandos, I. A., de Azcarraga, J. A., & Meliveo, C. (2012). Conformal higher spin theory in extended tensorial superspace. Fortschritte Phys.-Prog. Phys., 60(7-8), 861–867.
Abstract: We discuss the formulation of free conformal higher spin theories with extended N = 2, 4, 8 supersymmetry in N-extended tensorial superspaces. The superfield higher spin equations can be obtained by quantizing a superparticle model in N-extended tensorial superspace. The N-extended higher spin supermultiplets just contain scalar and spinor fields in tensorial space so that, in contrast with the standard (super)space approach, no nontrivial generalizations of the Maxwell or Einstein equations to tensorial space appear when N > 2. For N = 4, 8, the higher spin-tensorial components of the extended tensorial superfields are expressed through additional scalar and spinor fields in tensorial space which obey the same free higher spin equations, but that are axion-like in the sense that they possess Peccei-Quinn-like symmetries.
|
|
|
Baran, J. et al, & Brzezinski, K. (2024). Feasibility of the J-PET to monitor the range of therapeutic proton beams. Phys. Medica, 118, 103301–9pp.
Abstract: Purpose: The aim of this work is to investigate the feasibility of the Jagiellonian Positron Emission Tomography (J -PET) scanner for intra-treatment proton beam range monitoring. Methods: The Monte Carlo simulation studies with GATE and PET image reconstruction with CASToR were performed in order to compare six J -PET scanner geometries. We simulated proton irradiation of a PMMA phantom with a Single Pencil Beam (SPB) and Spread -Out Bragg Peak (SOBP) of various ranges. The sensitivity and precision of each scanner were calculated, and considering the setup's cost-effectiveness, we indicated potentially optimal geometries for the J -PET scanner prototype dedicated to the proton beam range assessment. Results: The investigations indicate that the double -layer cylindrical and triple -layer double -head configurations are the most promising for clinical application. We found that the scanner sensitivity is of the order of 10-5 coincidences per primary proton, while the precision of the range assessment for both SPB and SOBP irradiation plans was found below 1 mm. Among the scanners with the same number of detector modules, the best results are found for the triple -layer dual -head geometry. The results indicate that the double -layer cylindrical and triple -layer double -head configurations are the most promising for the clinical application, Conclusions: We performed simulation studies demonstrating that the feasibility of the J -PET detector for PET -based proton beam therapy range monitoring is possible with reasonable sensitivity and precision enabling its pre -clinical tests in the clinical proton therapy environment. Considering the sensitivity, precision and cost-effectiveness, the double -layer cylindrical and triple -layer dual -head J -PET geometry configurations seem promising for future clinical application.
|
|
|
n_TOF Collaboration(Barbagallo, M. et al), Domingo-Pardo, C., & Tain, J. L. (2018). Experimental setup and procedure for the measurement of the Be-7(n,p)Li-7 reaction at n_TOF. Nucl. Instrum. Methods Phys. Res. A, 887, 27–33.
Abstract: Following the completion of the second neutron beam line and the related experimental area (EAR2) at the n_TOF spallation neutron source at CERN, several experiments were planned and performed. The high instantaneous neutron flux available in EAR2 allows to investigate neutron induced reactions with charged particles in the exit channel even employing targets made out of small amounts of short-lived radioactive isotopes. After the successful measurement of the Be-7(n,alpha)alpha cross section, the Be-7(n,p)Li-7 reaction was studied in order to provide still missing cross section data of relevance for Big Bang Nucleosynthesis (BBN), in an attempt to find a solution to the cosmological Lithium abundance problem. This paper describes the experimental setup employed in such a measurement and its characterization.
|
|
|
Barenboim, G., Chen, J. Z., Hannestad, S., Oldengott, I. M., Tram, T., & Wong, Y. Y. Y. (2021). Invisible neutrino decay in precision cosmology. J. Cosmol. Astropart. Phys., 03(3), 087–53pp.
Abstract: We revisit the topic of invisible neutrino decay in the precision cosmological context, via a first-principles approach to understanding the cosmic microwave background and large-scale structure phenomenology of such a non-standard physics scenario. Assuming an effective Lagrangian in which a heavier standard-model neutrino nu(H) couples to a lighter one nu(l) and a massless scalar particle phi via a Yukawa interaction, we derive from first principles the complete set of Boltzmann equations, at both the spatially homogeneous and the firstorder inhomogeneous levels, for the phase space densities of nu(H), nu(l), and phi in the presence of the relevant decay and inverse decay processes. With this set of equations in hand, we perform a critical survey of recent works on cosmological invisible neutrino decay in both limits of decay while nu(H) is ultra-relativistic and non-relativistic. Our two main findings are: (i) in the non-relativistic limit, the effective equations of motion used to describe perturbations in the neutrino-scalar system in the existing literature formally violate momentum conservation and gauge invariance, and (ii) in the ultra-relativistic limit, exponential damping of the anisotropic stress does not occur at the commonly-used rate Gamma(T) = (1/tau(0))( m(nu H)/E-nu H)(3), but at a rate similar to (1/ tau(0))(m(nu H)/E-nu H)(5). Both results are model-independent. The impact of the former finding on the cosmology of invisible neutrino decay is likely small. The latter, however, implies a significant revision of the cosmological limit on the neutrino lifetime tau(0) from tau(old)(0) greater than or similar to 1.2 x 10(9) s (m(nu H)/50 meV)(3) to tau(0) greater than or similar to (4 x 10(5) -> 4 x 10(6)) s (m(nu H)/50 meV)(5).
|
|
|
Barenboim, G., & Park, W. I. (2016). New- vs. chaotic- inflations. J. Cosmol. Astropart. Phys., 02(2), 061–20pp.
Abstract: We show that “spiralized” models of new-inflation can be experimentally identified mostly by their positive spectral running in direct contrast with most chaotic-inflation models which have negative runnings typically in the range of O(10(-4)-10(-3)).
|
|