Design study of a novel geometrical arrangement for an in-beam small animal positron emission tomography scanner

Objective.We designed a geometrical solution for a small animal in-beam positron emission tomography (PET) scanner to be used in the project SIRMIO (Small animal proton irradiator for research in molecular image-guided radiation-oncology). The system is based on 56 scintillator blocks of pixelated LYSO crystals. The crystals are arranged providing a pyramidal-step shape to optimize the geometrical coverage in a spherical configuration.Approach.Different arrangements have been simulated and compared in terms of spatial resolution and sensitivity. The chosen setup enables us to reach a good trade-off between a solid angle coverage and sufficient available space for the integration of additional components of the first design prototype of the SIRMIO platform. The possibility of moving the mouse holder inside the PET scanner furthermore allows for achieving the optimum placement of the irradiation area for all the possible tumor positions in the body of the mouse. The work also includes a study of the scintillator material where LYSO and GAGG are compared with a focus on the random coincidence noise due to the natural radioactivity of Lutetium in LYSO, justifying the choice of LYSO for the development of the final system.Main results.The best imaging performance can be achieved with a sub-millimeter spatial resolution and sensitivity of 10% in the center of the scanner, as verified in thorough simulations of point sources. The simulation of realistic irradiation scenarios of proton beams in PMMA targets with/without air gaps indicates the ability of the proposed PET system to detect range shifts down to 0.2 mm.Significance.The presented results support the choice of the identified optimal design for a novel spherical in-beam PET scanner which is currently under commissioning for application to small animal proton and light ion irradiation, and which might find also application, e.g. for biological image-guidance in x-ray irradiation.

Observation of the radiative decay of the 229Th nuclear clock isomer

The radionuclide thorium-229 features an isomer with an exceptionally low excitation energy that enables direct laser manipulation of nuclear states. It constitutes one of the leading candidates for use in next-generation optical clocks^1-3. This nuclear clock will be a unique tool for precise tests of fundamental physics^4-9. Whereas indirect experimental evidence for the existence of such an extraordinary nuclear state is substantially older¹⁰, the proof of existence has been delivered only recently by observing the isomer's electron conversion decay¹¹. The isomer's excitation energy, nuclear spin and electromagnetic moments, the electron conversion lifetime and a refined energy of the isomer have been measured^12-16. In spite of recent progress, the isomer's radiative decay, a key ingredient for the development of a nuclear clock, remained unobserved. Here, we report the detection of the radiative decay of this low-energy isomer in thorium-229 (^229mTh). By performing vacuum-ultraviolet spectroscopy of ^229mTh incorporated into large-bandgap CaF₂ and MgF₂ crystals at the ISOLDE facility at CERN, photons of 8.338(24) eV are measured, in agreement with recent measurements^14-16 and the uncertainty is decreased by a factor of seven. The half-life of ^229mTh embedded in MgF₂ is determined to be 670(102) s. The observation of the radiative decay in a large-bandgap crystal has important consequences for the design of a future nuclear clock and the improved uncertainty of the energy eases the search for direct laser excitation of the atomic nucleus.

Sub-millimeter precise photon interaction position determination in large monolithic scintillators via convolutional neural network algorithms

In this work, we present the development and application of a convolutional neural network (CNN)-based algorithm to precisely determine the interaction position ofγ-quanta in large monolithic scintillators. Those are used as an absorber component of a Compton camera (CC) system under development for ion beam range verification via prompt-gamma imaging. We examined two scintillation crystals: LaBr₃:Ce and CeBr₃. Each crystal had dimensions of 50.8 mm × 50.8 mm × 30 mm and was coupled to a 64-fold segmented multi-anode photomultiplier tube (PMT) with an 8 × 8 pixel arrangement. We determined the spatial resolution for three photon energies of 662, 1.17 and 1.33 MeV obtained from 2D detector scans with tightly collimated¹³⁷Cs and⁶⁰Co photon sources. With the new algorithm we achieved a spatial resolution for the CeBr3 crystal below 1.11(8) mm and below 0.98(7) mm for the LaBr3:Ce detector for all investigated energies between 662 keV and 1.33 MeV. We thereby improved the performance by more than a factor of 2.5 compared to the previously used categorical average pattern algorithm, which is a variation of the well-established k-nearest neighbor algorithm. The trained CNN has a low memory footprint and enables the reconstruction of up to 10⁴events per second with only one GPU. Those improvements are crucial on the way to future clinicalin vivoapplicability of the CC for ion beam range verification.

Performance evaluation of a staggered three-layer DOI PET detector using a 1 mm LYSO pitch with PETsys TOFPET2 ASIC: comparison of HAMAMATSU and KETEK SiPMs

In this study, we propose a staggered three-layer depth-of-interaction (DOI) detector with a 1 mm crystal pitch and 19.8 mm total crystal thickness for a high-resolution and high-sensitivity small animal in-beam PET scanner. A three-layered stacked LYSO scintillation array (0.9 × 0.9 × 6.6 mm³crystals, 23 × 22 mm²surface area) read out by a SiPM array (8 × 8 channels, 3 × 3 mm²active area/channel and 50μm microcell size) with data acquisition, signal processing and digitization performed using the PETsys Electronics Evaluations kit (based on the TOFPET v2c ASIC) builds a DOI LYSO detector block. The performance of the DOI detector was evaluated in terms of crystal resolvability, energy resolution, and coincidence resolving time (CRT). A comparative performance evaluation of the staggered three-layer LYSO block was conducted with two different SiPM arrays from KETEK and HAMAMATSU. 100% (KETEK) and 99.8% (HAMAMATSU) of the crystals were identified, by using a flood irradiation the front- and back-side. The average energy resolutions for the 1st, 2nd, and 3rd layers were 16.5 (±2.3)%, 20.9(±4.0)%, and 32.7 (±21.0)% (KETEK) and 19.3 (±3.5)%, 21.2 (±4.1)%, and 26.6 (±10.3)% (HAMAMATSU) for the used SiPM arrays. The measured CRTs (FWHM) for the 1st, 2nd, and 3rd layers were 532 (±111) ps, 463 (±108) ps, and 447 (±111) ps (KETEK) and 402 (±46) ps, 392 (±54) ps, and 408 (±196) ps (HAMAMATSU). In conclusion, the performance of the staggered three-layer DOI detector with 1 mm LYSO pitch and 19.8 mm total crystal thickness was fully characterized. The feasibility of a highly performing readout of a high resolution DOI PET detector via SiPM arrays from KETEK and HAMAMATSU employing the PETsys TOFPET v2c ASIC could be demonstrated.

Radiation protection modelling for 2.5 Petawatt-laser production of ultrashort x-ray, proton and ion bunches: Monte Carlo model of the Munich CALA facility

The 'Centre for Advanced Laser Applications' (CALA) is a new research institute for laser-based acceleration of electron beams for brilliant x-ray generation, laser-driven sub-nanosecond bunches of protons and heavy ions for biomedical applications like imaging and tumour therapy as well as for nuclear and high-field physics.The radiation sources emerging from experiments using the up to 2.5 petawatt laser pulses with 25 femtosecond duration will be mixed particle-species of high intensity, high energy and pulsed, thus posing new challenges compared to conventional radiation protection. Such worldwide pioneering laser experiments result in source characteristics that require careful a-priori radiation safety simulations.The FLUKA Monte-Carlo code was used to model the five CALA experimental caves, including the corridors, halls and air spaces surrounding the caves. Beams of electrons (<5 GeV), protons (<200 MeV),¹²C (<400MeV/u) and¹⁹⁷Au (<10MeV/u) ions were simulated using spectra, divergences and bunch-charges based on expectations from recent scientific progress.Simulated dose rates locally can exceed 1.5 kSv h^-1inside beam dumps. Vacuum pipes in the cave walls for laser transport and extraction channels for the generated x-rays result in small dose leakage to neighboring areas. Secondary neutrons contribute to most of the prompt dose rate outside caves into which the beam is delivered. This secondary radiation component causes non-negligible dose rates to occur behind walls to which large fluences of secondary particles are directed.By employing adequate beam dumps matched to beam-divergence, magnets, passive shielding and laser pulse repetition limits, average dose rates in- and outside the experimental building stay below design specifications (<0.5μSv h-1) for unclassified areas,<2.5μSv h-1for supervised areas,<7.5μSv h-1maximum local dose rate) and regulatory limits (<1mSv a-1for unclassified areas).

Electronic Bridge Excitation in Highly Charged ^{229}Th Ions

The excitation of the 8 eV ^{229m}Th isomer through the electronic bridge mechanism in highly charged ions is investigated theoretically. By exploiting the rich level scheme of open 4f orbitals and the robustness of highly charged ions against photoionization, a pulsed high-intensity optical laser can be used to efficiently drive the nuclear transition by coupling it to the electronic shell. We show how to implement a promising electronic bridge scheme in an electron beam ion trap starting from a metastable electronic state. This setup would avoid the need for a tunable vacuum ultraviolet laser. Based on our theoretical predictions, determining the isomer energy with an uncertainty of 10^{-5}  eV could be achieved in one day of measurement time using realistic laser parameters.

Towards a novel small animal proton irradiation platform: the SIRMIO project

Background: Precision small animal radiotherapy research is a young emerging field aiming to provide new experimental insights into tumor and normal tissue models in different microenvironments, to unravel complex mechanisms of radiation damage in target and non-target tissues and assess efficacy of novel therapeutic strategies. For photon therapy, modern small animal radiotherapy research platforms have been developed over the last years and are meanwhile commercially available. Conversely, for proton therapy, which holds potential for an even superior outcome than photon therapy, no commercial system exists yet. Material and methods: The project SIRMIO (Small Animal Proton Irradiator for Research in Molecular Image-guided Radiation-Oncology) aims at realizing and demonstrating an innovative portable prototype system for precision image-guided small animal proton irradiation, suitable for installation at existing clinical treatment facilities. The proposed design combines precise dose application with in situ multi-modal anatomical image guidance and in vivo verification of the actual treatment delivery. Results and conclusions: This manuscript describes the status of the different components under development, featuring a dedicated beamline for degradation and focusing of clinical proton beams, along with novel detector systems for in situimaging and range verification. The foreseen workflow includes pre-treatment proton transmission imaging, complemented by ultrasonic tumor localization, for treatment planning and position verification, followed by image-guided delivery with on site range verification by means of ionoacoustics (for pulsed beams) and positron-emission-tomography (PET, for continuous beams). The proposed compact and cost-effective system promises to open a new era in small animal proton therapy research, contributing to the basic understanding of in vivo radiation action to identify areas of potential breakthroughs for future translation into innovative clinical strategies.

Preparing an Isotopically Pure 229Th Ion Beam for Studies of 229mTh

A methodology is described to generate an isotopically pure ²²⁹Th ion beam in the 2+ and 3+ charge states. This ion beam enables one to investigate the low-lying isomeric first excited state of ²²⁹Th at an excitation energy of about 7.8(5) eV and a radiative lifetime of up to 10⁴ seconds. The presented method allowed for a first direct identification of the decay of the thorium isomer, laying the foundations to study its decay properties as prerequisite for an optical control of this nuclear transition. High energy ²²⁹Th ions are produced in the α decay of a radioactive ²³³U source. The ions are thermalized in a buffer-gas stopping cell, extracted and subsequently an ion beam is formed. This ion beam is mass purified by a quadrupole-mass separator to generate a pure ion beam. In order to detect the isomeric decay, the ions are collected on the surface of a micro-channel plate detector, where electrons, as emitted in the internal conversion decay of the isomeric state, are observed.

A novel approach to electron data background treatment in an online wide-angle spectrometer for laser-accelerated ion and electron bunches

Laser-based ion acceleration is driven by electrical fields emerging when target electrons absorb laser energy and consecutively leave the target material. A direct correlation between these electrons and the accelerated ions is thus to be expected and predicted by theoretical models. We report on a modified wide-angle spectrometer, allowing the simultaneous characterization of angularly resolved energy distributions of both ions and electrons. Equipped with online pixel detectors, the RadEye1 detectors, the investigation of this correlation gets attainable on a single shot basis. In addition to first insights, we present a novel approach for reliably extracting the primary electron energy distribution from the interfering secondary radiation background. This proves vitally important for quantitative extraction of average electron energies (temperatures) and emitted total charge.

A Laser Excitation Scheme for ^{229m}Th

Direct laser excitation of the lowest known nuclear excited state in ^{229}Th has been a long-standing objective. It is generally assumed that reaching this goal would require a considerably reduced uncertainty of the isomer's excitation energy compared to the presently adopted value of (7.8±0.5)  eV. Here we present a direct laser excitation scheme for ^{229m}Th, which circumvents this requirement. The proposed excitation scheme makes use of already existing laser technology and therefore paves the way for nuclear laser spectroscopy. In this concept, the recently experimentally observed internal-conversion decay channel of the isomeric state is used for probing the isomeric population. A signal-to-background ratio of better than 10^{4} and a total measurement time of less than three days for laser scanning appear to be achievable.

Submillimeter ionoacoustic range determination for protons in water at a clinical synchrocyclotron

Proton ranges in water between 145 MeV to 227 MeV initial energy have been measured at a clinical superconducting synchrocyclotron using the acoustic signal induced by the ion dose deposition (ionoacoustic effect). Detection of ultrasound waves was performed by a very sensitive hydrophone and signals were stored in a digital oscilloscope triggered by secondary prompt gammas. The ionoacoustic range measurements were compared to existing range data from a calibrated range detector setup on-site and agreement of better than 1 mm was found at a Bragg peak dose of about 10 Gy for 220 MeV initial proton energy, compatible with the experimental errors. Ionoacoustics has thus the potential to measure the Bragg peak position with submillimeter accuracy during proton therapy, possibly correlated with ultrasound tissue imaging.

Lifetime Measurement of the ^{229}Th Nuclear Isomer

The first excited isomeric state of ^{229}Th possesses the lowest energy among all known excited nuclear states. The expected energy is accessible with today's laser technology and in principle allows for a direct optical laser excitation of the nucleus. The isomer decays via three channels to its ground state (internal conversion, γ decay, and bound internal conversion), whose strengths depend on the charge state of ^{229m}Th. We report on the measurement of the internal-conversion decay half-life of neutral ^{229m}Th. A half-life of 7±1  μs has been measured, which is in the range of theoretical predictions and, based on the theoretically expected lifetime of ≈10^{4}  s of the photonic decay channel, gives further support for an internal conversion coefficient of ≈10^{9}, thus constraining the strength of a radiative branch in the presence of internal conversion.

Comparative Characterization Study of a LaBr3(Ce) Scintillation Crystal in Two Surface Wrapping Scenarios: Absorptive and Reflective

The properties of a 50 mm × 50 mm × 30 mm monolithic LaBr₃:Ce scintillator crystal coupled to a position-sensitive multi-anode photomultiplier (PMT, Hamamatsu H9500), representing the absorbing detector of a Compton camera under study for online ion (proton) beam range verification in hadron therapy, was evaluated in combination with either absorptive or reflective crystal surface coating. This study covered an assessment of the energy and position-dependent energy resolution, exhibiting a factor of 2.5–3.5 improvement for the reflectively wrapped crystal at 662 keV. The spatial dependency was investigated using a collimated ¹³⁷Cs source, showing a steep degradation of the energy resolution at the edges and corners of the absorptively wrapped crystal. Furthermore, the time resolution was determined to be 273 ps (FWHM) and 536 ps (FWHM) with reflective and absorptive coating, respectively, using a ⁶⁰Co source. In contrast, the light spread function (LSF) of the light amplitude distribution on the PMT segments improved for the absorptively wrapped detector. Both wrapping modalities showed almost no differences in the energy-dependent photopeak detection efficiency.

Ionoacoustic characterization of the proton Bragg peak with submillimeter accuracy

PURPOSE

Range verification in ion beam therapy relies to date on nuclear imaging techniques which require complex and costly detector systems. A different approach is the detection of thermoacoustic signals that are generated due to localized energy loss of ion beams in tissue (ionoacoustics). Aim of this work was to study experimentally the achievable position resolution of ionoacoustics under idealized conditions using high frequency ultrasonic transducers and a specifically selected probing beam.

METHODS

A water phantom was irradiated by a pulsed 20 MeV proton beam with varying pulse intensity and length. The acoustic signal of single proton pulses was measured by different PZT-based ultrasound detectors (3.5 and 10 MHz central frequencies). The proton dose distribution in water was calculated by Geant4 and used as input for simulation of the generated acoustic wave by the matlab toolbox k-WAVE.

RESULTS

In measurements from this study, a clear signal of the Bragg peak was observed for an energy deposition as low as 10(12) eV. The signal amplitude showed a linear increase with particle number per pulse and thus, dose. Bragg peak position measurements were reproducible within ±30 μm and agreed with Geant4 simulations to better than 100 μm. The ionoacoustic signal pattern allowed for a detailed analysis of the Bragg peak and could be well reproduced by k-WAVE simulations.

CONCLUSIONS

The authors have studied the ionoacoustic signal of the Bragg peak in experiments using a 20 MeV proton beam with its correspondingly localized energy deposition, demonstrating submillimeter position resolution and providing a deep insight in the correlation between the acoustic signal and Bragg peak shape. These results, together with earlier experiments and new simulations (including the results in this study) at higher energies, suggest ionoacoustics as a technique for range verification in particle therapy at locations, where the tumor can be localized by ultrasound imaging. This acoustic range verification approach could offer the possibility of combining anatomical ultrasound and Bragg peak imaging, but further studies are required for translation of these findings to clinical application.

Direct mapping of nuclear shell effects in the heaviest elements

Quantum-mechanical shell effects are expected to strongly enhance nuclear binding on an "island of stability" of superheavy elements. The predicted center at proton number Z = 114, 120, or 126 and neutron number N = 184 has been substantiated by the recent synthesis of new elements up to Z = 118. However, the location of the center and the extension of the island of stability remain vague. High-precision mass spectrometry allows the direct measurement of nuclear binding energies and thus the determination of the strength of shell effects. Here, we present such measurements for nobelium and lawrencium isotopes, which also pin down the deformed shell gap at N = 152.

Mass measurements of very neutron-deficient Mo and Tc isotopes and their impact on rp process nucleosynthesis

The masses of ten proton-rich nuclides, including the N=Z+1 nuclides ⁸⁵Mo and ⁸⁷Tc, were measured with the Penning trap mass spectrometer SHIPTRAP. Compared to the Atomic Mass Evaluation 2003 a systematic shift of the mass surface by up to 1.6 MeV is observed causing significant abundance changes of the ashes of astrophysical x-ray bursts. Surprisingly low α separation energies for neutron-deficient Mo and Tc are found, making the formation of a ZrNb cycle in the rp process possible. Such a cycle would impose an upper temperature limit for the synthesis of elements beyond Nb in the rp process.

Discovery of the shape coexisting 0+ state in 32 Mg by a two neutron transfer reaction

The "island of inversion" nucleus 32 Mg has been studied by a (t, p) two neutron transfer reaction in inverse kinematics at REX-ISOLDE. The shape coexistent excited 0+ state in 32 Mg has been identified by the characteristic angular distribution of the protons of the Δ L=0 transfer. The excitation energy of 1058 keV is much lower than predicted by any theoretical model. The low γ-ray intensity observed for the decay of this 0+ state indicates a lifetime of more than 10 ns. Deduced spectroscopic amplitudes are compared with occupation numbers from shell-model calculations.

Shape coexistence near neutron number N=20: first identification of the E0 decay from the deformed first excited Jpi=0+ state in 30Mg

The 1789 keV state in 30Mg was identified as the first excited 0+ state via its electric monopole (E0) transition to the ground state. The measured small value of rho2(E0,0(2)+-->0(1)+)=(26.2+/-7.5)x10(-3) implies within a two-level model a small mixing of competing configurations with largely different intrinsic quadrupole deformation near the neutron shell closure at N=20. Axially symmetric configuration mixing calculations identify the ground state of 30Mg to be based on neutron configurations below the N=20 shell closure, while the excited 0+ state mainly consists of two neutrons excited into the nu 1f7/2 orbital. The experimental result represents the first case where an E0 back decay from a strongly deformed second to the normal deformed first nuclear potential minimum well has been unambiguously identified, thus directly proving shape coexistence at the borderline of the much-debated "island of inversion."

First Penning trap mass measurements beyond the proton drip line

The masses of six neutron-deficient rare holmium and thulium isotopes close to the proton drip line were determined with the SHIPTRAP Penning trap mass spectrometer. For the first time the masses of the proton-unbound isotopes 144,145Ho and 147,148Tm were directly measured. The proton separation energies were derived from the measured mass values and compared to predictions from mass formulas. The new values of the proton separation energies are used to determine the location of the proton drip line for holmium and thulium more accurately.

Excitation energies of superdeformed States in 196Pb: towards a systematic study of the second well in Pb isotopes

The excitation energy of the lowest-energy superdeformed band in 196Pb is established using the techniques of time-correlated gamma-ray spectroscopy. Together with previous measurements on 192Pb and 194Pb, this result allows superdeformed excitation energies, binding energies, and two-proton and two-neutron separation energies to be studied systematically, providing stringent tests for current nuclear models. The results are examined for evidence of a "superdeformed shell gap."

"Safe" Coulomb excitation of 30Mg

We report on the first radioactive beam experiment performed at the recently commissioned REX-ISOLDE facility at CERN in conjunction with the highly efficient gamma spectrometer MINIBALL. Using 30Mg ions accelerated to an energy of 2.25 MeV/u together with a thin (nat)Ni target, Coulomb excitation of the first excited 2+ states of the projectile and target nuclei well below the Coulomb barrier was observed. From the measured relative deexcitation gamma-ray yields the B(E2;0(+)gs-->2(+)1) value of 30Mg was determined to be 241(31)e2 fm4. Our result is lower than values obtained at projectile fragmentation facilities using the intermediate-energy Coulomb excitation method, and confirms the theoretical conjecture that the neutron-rich magnesium isotope 30Mg resides outside the "island of inversion."