Search

Purpose: Beam angle selection is critical in proton therapy treatment planning, yet automated approaches remain underexplored. This study presents and evaluates GAMBAS, a novel, fast machine learning model for automatic beam angle selection. Methods: The model extracts a predefined set of anatomical features from a patient's CT and structure contours. Using these features, it identifies the most similar patient from a training database and suggests that patient's beam arrangement. A retrospective study with 19 patients was conducted, comparing this model's suggestions to human planners' choices and randomly selected beam arrangements from the training dataset. An expert treatment planner evaluated the plans on quality (scale 1-5), ranked them, and guessed the method used. Results: The number of acceptable (score 4 or 5) plans was comparable between human-chosen 17 (89%) and model-selected 16(84%) beam arrangements. The fully automatic treatment planning took between 4 - 7 min (mean 5 min). Conclusion: The model produces beam arrangements of comparable quality to those chosen by human planners, demonstrating its potential as a fast tool for quality assurance and patient selection, although it is not yet ready for clinical use., Comment: 8 pages, 2 figures, 1 table

The importance of computers is continually increasing in radiotherapy. Efficient algorithms, implementations and the ability to leverage advancements in computer science are crucial to improve cancer care even further and deliver the best treatment to each patient. Yet, the software landscape for radiotherapy is fragmented into proprietary systems that do not share a common interface. Further, the radiotherapy community does not have access to the vast possibilities offered by modern programming languages and their ecosystem of libraries yet. We present JulianA.jl, a novel Julia package for radiotherapy. It aims to provide a modular and flexible foundation for the development and efficient implementation of algorithms and workflows for radiotherapy researchers and clinicians. JulianA.jl can be interfaced with any scriptable treatment planning system, be it commercial, open source or in-house developed. This article highlights our design choices and showcases the package's simplicity and powerful automatic treatment planning capabilities.

We derive the generalized second law (GSL) for arbitrary cuts of Killing horizons from the perspective of crossed-product gravitational algebras, making use of a recent proposal by one of us for the construction of local gravitational algebras. This construction relies on the existence of a state whose modular flow is geometric on the horizon. In both free and interacting quantum field theories, such states are guaranteed to exist by the properties of half-sided translations on the horizon. Using geometric identities derived from the canonical analysis of general relativity on null surfaces, we show that the crossed product entropy agrees with the generalized entropy of the horizon cut in a semiclassical limit, and further reproduce Wall's result relating the GSL to monotonicity of relative entropy of the quantum field algebras. We also give a novel generalization of the GSL for interacting theories in asymptotically flat spacetimes involving the concept of an algebra at infinity for a half-sided translation, which accounts for triviality of the algebra of fields smeared only on the horizon. Going beyond the semiclassical limit, we compute subleading corrections to the crossed product entropy, but are unable to determine if the GSL continues to hold after accounting for these. We speculate that an improved GSL could follow from a hidden subalgebra structure of the crossed products, assuming the existence of an operator-valued weight between horizon cut algebras., Comment: 62 pages + appendices and references, 4 figures; v2 minor edits, references added

Background: In MR-guided proton therapy planning, aligning MR and CT images is key for MR-based CT synthesis, especially in mobile regions like the head-and-neck. Misalignments here can lead to less accurate synthetic CT (sCT) images, impacting treatment precision. Purpose: This study introduces a novel network that cohesively unifies image generation and registration processes to enhance the quality and anatomical fidelity of sCTs derived from better-aligned MR images. Methods: The approach synergizes a generation network (G) with a deformable registration network (R), optimizing them jointly in MR-to-CT synthesis. This goal is achieved by alternately minimizing the discrepancies between the generated/registered CT images and their corresponding reference CT counterparts. The generation network employs a UNet architecture, while the registration network leverages an implicit neural representation of the Deformable Vector Fields (DVFs). We validated this method on a dataset comprising 60 Head-and-Neck patients, reserving 12 cases for holdout testing. Results: Compared to the baseline Pix2Pix method with MAE 124.95\pm 30.74 HU, the proposed technique demonstrated 80.98\pm 7.55 HU. The unified translation-registration network produced sharper and more anatomically congruent outputs, showing superior efficacy in converting MR images to sCTs. Additionally, from a dosimetric perspective, the plan recalculated on the resulting sCTs resulted in a remarkably reduced discrepancy to the reference proton plans. Conclusions: This study conclusively demonstrates that a holistic MR-based CT synthesis approach, integrating both image-to-image translation and deformable registration, significantly improves the precision and quality of sCT generation, particularly for the challenging body area with varied anatomic changes between corresponding MR and CT.

While von Neumann entropies for subregions in quantum field theory universally contain ultraviolet divergences, differences between von Neumann entropies are finite and well-defined in many physically relevant scenarios. We demonstrate that such a notion of entropy differences can be rigorously defined in quantum field theory in a general curved spacetime by introducing a novel, covariant regulator for the entropy based on the modular crossed product. This regulator associates a type II von Neumann algebra to each spacetime subregion, resulting in well-defined renormalized entropies. This prescription reproduces formulas for entropy differences that coincide with heuristic formulas widely used in the literature, and we prove that it satisfies desirable properties such as unitary invariance and concavity. As an application, we provide proofs of the Bekenstein bound and the quantum null energy condition, formulated directly in terms of vacuum-subtracted von Neumann entropies., Comment: 6 pages + appendices, 2 figures

Creating high quality treatment plans is crucial for a successful radiotherapy treatment. However, it demands substantial effort and special training for dosimetrists. Existing automated treatment planning systems typically require either an explicit prioritization of planning objectives, human-assigned objective weights, large amounts of historic plans to train an artificial intelligence or long planning times. Many of the existing auto-planning tools are difficult to extend to new planning goals. A new spot weight optimisation algorithm, called JulianA, was developed. The algorithm minimises a scalar loss function that is built only based on the prescribed dose to the tumour and organs at risk (OARs), but does not rely on historic plans. The objective weights in the loss function have default values that do not need to be changed for the patients in our dataset. The system is a versatile tool for researchers and clinicians without specialised programming skills. Extending it is as easy as adding an additional term to the loss function. JulianA was validated on a dataset of 19 patients with intra- and extracerebral neoplasms within the cranial region that had been treated at our institute. For each patient, a reference plan which was delivered to the cancer patient, was exported from our treatment database. Then JulianA created the auto plan using the same beam arrangement. The reference and auto plans were given to a blinded independent reviewer who assessed the acceptability of each plan, ranked the plans and assigned the human-/machine-made labels. The auto plans were considered acceptable in 16 out of 19 patients and at least as good as the reference plan for 11 patients. Whether a plan was crafted by a dosimetrist or JulianA was only recognised for 9 cases. The median time for the spot weight optimisation is approx. 2 min (range: 0.5 min - 7 min).

Objective: 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. Approach: The Monte Carlo simulation studies with GATE and PET image reconstruction with CASToR were performed in order to compare six J-PET scanner geometries (three dual-heads and three cylindrical). 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. Main 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. Significance: 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 the future clinical application. Experimental tests are needed to confirm these findings., Comment: 19 pages, 7 pages

Gravitational subsystems with boundaries carry the action of an infinite-dimensional symmetry algebra, with potentially profound implications for the quantum theory of gravity. We initiate an investigation into the quantization of this corner symmetry algebra for the phase space of gravity localized to a region bounded by a 2-dimensional sphere. Starting with the observation that the algebra $\mathfrak{sdiff}(S^2)$ of area-preserving diffeomorphisms of the 2-sphere admits a deformation to the finite-dimensional algebra $\mathfrak{su}(N)$, we derive novel finite-$N$ deformations for two important subalgebras of the gravitational corner symmetry algebra. Specifically, we find that the area-preserving hydrodynamical algebra $\mathfrak{sdiff}(S^2)\oplus_{\mathcal{L}}\mathbb{R}^{S^2}$ arises as the large-$N$ limit of $\mathfrak{sl}(N,\mathbb C)\oplus\mathbb{R}$ and that the full area-preserving corner symmetry algebra $\mathfrak{sdiff}(S^2)\oplus_{\mathcal{L}}\mathfrak{sl}(2,\mathbb{R})^{S^2}$ is the large-$N$ limit of the pseudo-unitary group $\mathfrak{su}(N,N)$. We find matching conditions for the Casimir elements of the deformed and continuum algebras and show how these determine the value of the deformation parameter $N$ as well as the representation of the deformed algebra associated with a quantization of the local gravitational phase space. Additionally, we present a number of novel results related to the various algebras appearing, including a detailed analysis of the asymptotic expansion of the $\mathfrak{su}(N)$ structure constants, as well as an explicit computation of the full $\mathfrak{diff}(S^2)$ structure constants in the spherical harmonic basis. A consequence of our work is the definition of an area operator which is compatible with the deformation of the area-preserving corner symmetry at finite $N$., Comment: v3: 57 pages + Appendices + References = 99 pages; Typos fixed; Published in JHEP

We expand on our results in arXiv:2111.02429 to present a broad new class of gravitational observables in asymptotically Anti-de Sitter space living on general codimension-zero regions of the bulk spacetime. By taking distinct limits, these observables can reduce to well-studied holographic complexity proposals, e.g., the volume of the maximal slice and the action or spacetime volume of the Wheeler-DeWitt patch. As with the codimension-one family found in arXiv:2111.02429, these new observables display two key universal features for the thermofield double state: they grow linearly in time at late times and reproduce the switchback effect. Hence we argue that any member of this new class of observables is an equally viable candidate as a gravitational dual of complexity. Moreover, using the Peierls construction, we show that variations of the codimension-zero and codimension-one observables are encoded in the gravitational symplectic form on the semi-classical phase-space, which can then be mapped to the CFT., Comment: 53 pages + 6 appendices

Recently, Ciambelli, Leigh, and Pai (CLP) [arXiv:2111.13181] have shown that nonzero charges integrating Hamilton's equation can be defined for all diffeomorphisms acting near the boundary of a subregion in a gravitational theory. This is done by extending the phase space to include a set of embedding fields that parameterize the location of the boundary. Because their construction differs from previous works on extended phase spaces by a covariant phase space ambiguity, the question arises as to whether the resulting charges are unambiguously defined. Here, we demonstrate that ambiguity-free charges can be obtained by appealing to the variational principle for the subregion, following recent developments on dealing with boundaries in the covariant phase space. Resolving the ambiguity produces corrections to the diffeomorphism charges, and also generates additional obstructions to integrability of Hamilton's equation. We emphasize the fact that the CLP extended phase space produces nonzero diffeomorphism charges distinguishes it from previous constructions in which diffeomorphisms are pure gauge, since the embedding fields can always be eliminated from the latter by a choice of unitary gauge. Finally, we show that Wald-Zoupas charges, with their characteristic obstruction to integrability, are associated with a modified transformation in the extended phase space, clarifying the reason behind integrability of Hamilton's equation for standard diffeomorphisms., Comment: 22 pages; v2, updated discussion, references added

We develop a general framework for constructing charges associated with diffeomorphisms in gravitational theories using covariant phase space techniques. This framework encompasses both localized charges associated with spacetime subregions, as well as global conserved charges of the full spacetime. Expressions for the charges include contributions from the boundary and corner terms in the subregion action, and are rendered unambiguous by appealing to the variational principle for the subregion, which selects a preferred form of the symplectic flux through the boundaries. The Poisson brackets of the charges on the subregion phase space are shown to reproduce the bracket of Barnich and Troessaert for open subsystems, thereby giving a novel derivation of this bracket from first principles. In the context of asymptotic boundaries, we show that the procedure of holographic renormalization can be always applied to obtain finite charges and fluxes once suitable counterterms have been found to ensure a finite action. This enables the study of larger asymptotic symmetry groups by loosening the boundary conditions imposed at infinity. We further present an algorithm for explicitly computing the counterterms that renormalize the action and symplectic potential, and, as an application of our framework, demonstrate that it reproduces known expressions for the charges of the generalized Bondi-Metzner-Sachs algebra., Comment: 57 pages + 29 pages appendices and references, 1 figure; v2 added 3 tables, additional references, and new discussion of Casimir energy in AdS

We present a new infinite class of gravitational observables in asymptotically Anti-de Sitter space living on codimension-one slices of the geometry, the most famous of which is the volume of the maximal slice. We show that these observables display universal features for the thermofield-double state: they grow linearly in time at late times and reproduce the switch-back effect in shock wave geometries. We argue that any member of this class of observables is an equally viable candidate as the extremal volume for a gravitational dual of complexity., Comment: 6 pages + appendices, 3 figures; v2 minor edits and title change - version as in PRL

Time series data mining is an important field of research in the era of "Big Data". Next generation astronomical surveys will generate data at unprecedented rates, creating the need for automated methods of data analysis. We propose a method of light curve characterisation that employs a pipeline consisting of a neural network with a Long-Short Term Memory Variational Autoencoder architecture and a Gaussian mixture model. The pipeline performs extraction and aggregation of features from light curve segments into feature vectors of fixed length which we refer to as light curve "fingerprints". This representation can be readily used as input of down-stream machine learning algorithms. We demonstrate the proposed method on a data set of Rossi X-ray Timing Explorer observations of the galactic black hole X-ray binary GRS 1915+105, which was chosen because of its observed complex X-ray variability. We find that the proposed method can generate a representation that characterises the observations and reflects the presence of distinct classes of GRS 1915+105 X-ray flux variability. We find that this representation can be used to perform efficient classification of light curves. We also present how the representation can be used to quantify the similarity of different light curves, highlighting the problem of the popular classification system of GRS 1915+105 observations, which does not account for intermediate class behaviour.

The Brown-York stress tensor provides a means for defining quasilocal gravitational charges in subregions bounded by a timelike hypersurface. We consider the generalization of this stress tensor to null hypersurfaces. Such a stress tensor can be derived from the on-shell subregion action of general relativity associated with a Dirichlet variational principle, which fixes an induced Carroll structure on the null boundary. The formula for the mixed-index tensor $T{}^i{}_j$ takes a remarkably simple form that is manifestly independent of the choice of auxiliary null vector at the null surface, and we compare this expression to previous proposals for null Brown-York stress tensors. The stress tensor we obtain satisfies a covariant conservation equation with respect to any connection induced from a rigging vector at the hypersurface, as a result of the null constraint equations. For transformations that act covariantly on the boundary structures, the Brown-York charges coincide with canonical charges constructed from a version of the Wald-Zoupas procedure. For anomalous transformations, the charges differ by an intrinsic functional of the boundary geometry, which we explicity verify for a set of symmetries associated with finite null hypersurfaces. Applications of the null Brown-York stress tensor to symmetries of asymptotically flat spacetimes and celestial holography are discussed., Comment: 22 pages + appendix; v2: fixed sign conventions, references added; v3 references added, matches published version

Investigations that were carried out over the last two decades with novel and more sensitive instrumentation have dramatically improved our knowledge of the more violent physical processes taking place in galactic and extra-galactic Black-Holes, Neutron Stars, Supernova Remnants/Pulsar Wind Nebulae, and other regions of the Universe where relativistic acceleration processes are in place. Simultaneous and/or combined observations with gamma-ray satellites and ground based high-energy telescopes, have clarified the scenario of the mechanisms responsible for high energy photon emission by leptonic and hadronic accelerated particles in the presence of magnetic fields. Specifically, the European Space Agency INTEGRAL soft gamma-ray observatory has detected more than 1000 sources in the soft gamma-ray band, providing accurate positions, light curves and time resolved spectral data for them. Space observations with Fermi-LAT and observations that were carried out from the ground with H.E.S.S., MAGIC, VERITAS, and other telescopes sensitive in the GeV-TeV domain have, at the same time, provided evidence that a substantial fraction of the cosmic sources detected are emitting in the keV to TeV band via Synchrotron-Inverse Compton processes, in particular from stellar galactic BH systems as well as from distant black holes. In this work, employing a spatial cross correlation technique, we compare the INTEGRAL/IBIS and TeV all-sky data in search of secure or likely associations. Although this analysis is based on a subset of the INTEGRAL all-sky observations, we find that there is a significant correlation: 39 objects show emission in both soft gamma-ray and TeV wavebands. The full INTEGRAL database, now comprising almost 19 years of public data available, will represent an important legacy that will be useful for the Cherenkov Telescope Array (CTA) and other ground based large projects., Comment: 30 pages, 10 figures, accepted for publication in MDPI - Universe Special Issue "High-Energy Gamma-Ray Astronomy: Results on Fundamental Questions after 30 Years of Ground-Based Observations"