155Tb production by cyclotrons: what level of 155Gd enrichment allows clinical applications?

Background 155Tb represents a potentially useful radionuclide for diagnostic medical applications, but its production remains a challenging problem, in spite of the fact that many production routes have been already investigated and tested. A recent experimental campaign, conducted with low-energy proton beams impinging on a 155Gd target with 91.9% enrichment, demonstrated a significant co-production of 156gTb, a contaminant of great concern since its half-life is comparable to that of 155Tb and its high-energy γ emissions severely impact on the dose released and on the quality of the SPECT images. In the present investigation, the isotopic purity of the enriched 155Gd target necessary to minimize the co-production of contaminant radioisotopes, in particular 156gTb, was explored using various computational simulations. Results Starting from the recent experimental data obtained with a 91.9% 155Gd-enriched target, the co-production of other Tb radioisotopes besides 155Tb has been theoretically evaluated using the Talys code. It was found that 156Gd, with an isotopic content of 5.87%, was the principal contributor to the co-production of 156gTb. The analysis also demonstrated that the maximum amount of 156Gd admissible for 155Tb production with a radionuclidic purity higher than 99% was 1%. A less stringent condition was obtained through computational dosimetry analysis, suggesting that a 2% content of 156Gd in the target can be tolerated to limit the dose increase to the patient below the 10% limit. Moreover, it has been demonstrated that the imaging properties of the produced 155Tb are not severely affected by this level of impurity in the target. Conclusions 155Tb can be produced with a quality suitable for medical applications using low-energy proton beams and 155Gd-enriched targets, if the 156Gd impurity content does not exceed 2%. Under these conditions, the dose increase due to the presence of contaminant radioisotopes remains below the 10% limit and good quality images, comparable to those of 111In, are guaranteed.


Background
Terbium is one of the few rare earth radiometals that could be used in nuclear medicine for tumor diagnosis and treatment, due to the favorable physical decay properties such as half-lives, type and energy of emissions.The four Tb radionuclides with higher clinical interest are 152 Tb and 155 Tb, positron and γ emitters respectively, relevant for diagnostic purposes, and 149 Tb and 161 Tb, α and β − emitters respectively, suitable for therapeutic applications [1].The matched pairs therapeutic/diagnostic radioisotopes of this element could be ideal candidates for theranostic applications, with no differences in chemical and pharmacokinetics behaviours [2].Indeed, significant differences could be observed when using pairs of diverse elements for diagnostics and therapy, because it is well known that each metal ion has specific chemical demands arising from its fundamental characteristics, giving rise to different coordination and geometry numbers [3], which influence the pharmacological properties of the labelled molecules.
Positron emission tomography (PET) imaging is the clinical nuclear imaging technique with higher sensitivity and resolution.Although, single photon emission computed tomography (SPECT) is still the predominant technology worldwide, due to the easy availability of γ-emitter radionuclides and the low cost of SPECT camera compared to PET scanner [4].Therefore, the production of new γ-emitting radionuclides that could be detected with the large number of available SPECT scanners is encouraged.Among Tb isotopes, 155 Tb is the only suitable for SPECT imaging, thanks to its two γ emissions at 87 keV (32%) and 105 keV (25%).Preclinical studies in nude mice bearing tumor xenograft demonstrated that 155 Tb-labelled biomolecules were able to visualize tumor sites using a small-animal SPECT/CT scanner even 7 days after radiocomplexes administration [5]. 155Tb has been produced by different methods, although most of them generated a large number of isotopes.First investigations on the production of 155 Tb via proton-induced reactions on Gd targets date back to 1989 with the work by Dmitriev et al. [6].In this case, no cross-section measurements were performed, but rather irradiation experiments of nat Gd thick targets for yields and activities, with incident proton energies between 11 and 22 MeV.Using nat Gd targets, cross-section measurements were performed much later by Vermeulen et al. [7].Besides the detailed experimental work that included the measurement of cross sections for the formation of a variety of terbium radionuclides, theoretical calculations were also performed for all contributing reactions.Since nat Gd has seven stable isotopes, it is not suited as target material for radioisotope production with high radionuclidic purity (RNP).Therefore, comparisons were made with accurate theoretical simulations in order to identify the most promising reactions that could be effectively employed in studies with enriched targets [7].More recently Formento-Cavaier et al. [8] extended the measurements with nat Gd targets up to 70 MeV.Experimental investigations with enriched targets of 155 Gd and 156 Gd have been undertaken by Favaretto et al. [9], with the result that 156 Gd(p,2n) 155 Tb provides high production yields, but implies the use of higher-energy cyclotrons and a significant contamination by the co-production of 156 Tb.The 155 Gd(p,n) 155 Tb reaction, on the other hand, can be performed with medical cyclotrons, with incident protons up to around 18 or 20 MeV, and has the potential to provide significant yields with high purity.In a subsequent publication, Dellepiane et al. [10] used the same enriched gadolinium oxide ( 155 Gd 91.9% enrichment or 156 Gd 93.3% enrichment) as target materials, measuring the 155 Tb production cross sections as well as a variety of contaminants produced for such highly enriched targets.In the most favorable case of the 155 Gd target, using an input energy of 10.5 MeV, contamination from 156 Tb was still significant, leading to a maximum RNP not greater than 93% after a decay time of about 96 h.Thus, it was suggested, as possible solution, to purify the final product through an off-line mass separation technique at the expense of a lower production yield, due the low efficiency of the current mass separation approaches.The production of 155 Tb and other Tb-isotopes ( 149 Tb, 152 Tb and 161 Tb) has been explored also at CERN-ISOLDE using spallation of high-energy proton beams on Tantalum targets, followed by ionization and mass separation [1,11].The drawback stands in the limited quantity of 155 Tb that can be produced with this method, due to the observed low efficiency of the accumulation procedure [12,13].It has been suggested also to irradiate a 159 Tb target by intermediate energy (60, 70 MeV) proton beams [14,15].This would open the possibility to have a 155 Dy/ 155 Tb generator system, similarly to the renowned 99 Mo/ 99m Tc one.However, the double separation chemistry among lanthanides represents a crucial step that still needs to be solved and the possible co-production of very long-lived terbium contaminants could imply a too low specific activity.
In this work Talys calculations of 155 Tb cross section have been benchmarked with data measured by Dellepiane et al. [10] with an enriched 155 Gd target.The contribution of each isotopic component has been disentangled to expose the effects deriving from the impurities of the target.Once the modeling has been tested on the specific isotopic 155 Gd abundance of the enriched target used by Dellepiane et al., the level of enrichment of the 155 Gd target necessary to produce 155 Tb with high purity has been investigated.However, the assessment of RNP is not enough to evaluate the suitability of the produced 155 Tb for clinical purposes.In fact, different produced contaminants can have a different impact on the dose delivered to the patient and on the quality of the SPECT images.Detailed biodistribution data of the DOTA-folate conjugate 161 Tb-cm09 in IGROV-1 tumour-bearing mice have been reported [11].As the biodistribution is independent of the radioisotope used to label the molecule, these data can be used to estimate the dosimetric properties of the same radiopharmaceutical when labeled with different Tb radioisotopes.In this way, it was possible to assess the dose increase (DI) [16] to the patient following 155 Tb-cm09 injection, due to the presence of contaminants in the 155 Tb produced supposing different enrichment of the 155 Gd target.Moreover, the effect of the contaminants with high-energy γ emissions on the SPECT images quality was evaluated.

Cross sections and thick target irradiation
The study of the nuclear reaction routes implies the adoption of different models to consider both the compound nucleus formation/decay and pre-equilibrium dynamics.To this purpose the Talys code [17] has been used, specifically version 1.95 which includes the geometry dependent hybrid (GDH) model [18,19].To describe the nuclear reaction mechanisms this Talys version provides 5 pre-equilibrium (PE) and 6 level-density (LD) models, for a total of 30 possible combinations of models.The description of the PE processes is mainly based on the exciton model.Transition matrices are the basic building blocks of the exciton formalism and they are described (1) analytically; (2) numerically; or (3) derived from the imaginary component of the optical potential.A fully quantummechanical approach, based upon the Feshbach-Kerman-Koonin theory and alternative to the exciton model, is included and labelled as (4).The already mentioned GDH model, which considers nuclear surface effects in the exciton model, is denoted as option (5).Another important ingredient for the compound nucleus formation is the nuclear LD, which affects the cross section through the Hauser-Feshbach formalism.Talys considers three phenomenological LD models, namely (1) the Fermi gas with constant temperature; (2) the back-shifted Fermi gas; (3) the generalized superfluid model; and three microscopic models, that is (4) the Goriely's tabulated Hartree-Fock densities; (5) the Hilaire's tabulation based upon the combinatorial model, and (6) the Hartree-Fock-Bogoliubov temperature-dependent formalism.The variability of these models can be described statistically introducing an interquartile band, measuring the model dispersion between the lower Q1 and the upper Q3 quartile, as has been discussed in a recent paper [20].In addition, the modeling has been tested on the specific isotopic abundances of the enriched 155 Gd target used by Dellepiane et al.
The modeled cross sections have been employed to evaluate the radionuclides produced under specific irradiation conditions, including the level of enrichment of the target.The computational approach discussed by Canton et al. [21] has been used for the evaluation of rates, activities, yields, and purities.
The rate R of production of a radionuclide from a beam colliding on a thick target can be derived from the expression where I 0 is the beam current, z proj is the projectile charge (1 for a proton beam), e the electron charge, N a the Avogadro number, A the atomic mass of the target element, E in and E out the energy of the beam hitting the target and the one leaving the target after traveling through its thickness, respectively.The production cross section for the nuclide is σ(E), while the target density ρ t and the stopping power of the projectile in the target dE/dx, described by the Bethe-Bloch formula (22).
Once the rates given in Eq. ( 1) for all the Tb radionuclides of interest are determined, the time evolution of the number of produced radionuclides, and the evolution of the activities during and after an irradiation can be calculated from the Bateman equations.The time dependence of the number of Tb radionuclides produced and their corresponding activities are used to evaluate the isotopic and radionuclidic purities for the production of 155 Tb.The main decay data of interest for the present analysis are reported in Table 1.

Assessment of organ absorbed doses and effective dose due to xxx Tb-cm09 injection
Dosimetric assessment were carried out using biodistribution data of the DOTA-folate conjugate 161 Tb-cm09 in IGROV-1 tumour-bearing mice [11].The cm09 is composed by a targeting vector (which selectively binds to the folate receptor expressed on a variety of tumour types) conjugated to small-molecular-weight albumin (which improves the blood circulation time and tissue distribution profile of folate conjugates) and to the DOTA chelating agent.Dosimetric evaluation has been performed supposing the cm09 (1) labelled with 155 Tb and also with other Tb-radioisotopes expected to be produced by proton irradiation of 155 Gd-enriched targets.Biodistribution data of 161 Tb-cm09, acquired in a time window of 7 days post injection in IGROV-1 tumour-bearing female nude mice [11], were used to estimate absorbed doses in humans due to the various Tb radioisotopes produced during the irradiation of all the considered enriched 155 Gd target.This evaluation was done through the relative mass scaling method, which takes into account the differences in human (H) and animal (A) organ masses compared to the total body masses [25].The activity concentrations in the different animal source organs (blood, lung, spleen, kidneys, stomach, intestines, liver, salivary glands, muscle and bone), reported as per cent of injected activity per gram of tissue ([%IA/g] A ), were scaled from mice to adult male humans to obtain the decay-corrected per cent of injected activity for each human source organ ([%IA/ organ] H ) through the following formula: with OW H the weight of human organ, TBW A and TBW H the total body weight for animal and human, respectively.OW H and TBW H values were obtained from the adult male phantom implemented in the Organ Level Internal Dose Assessment (OLINDA) software code [26].Biodistribution data were then plotted as a function of post injection time and fitted with CoKiMo software [26] by a tri-exponential equation, representing the phase of accumulation and the possibility of both a fast and a slow elimination of the radiopharmaceutical [27].At last, the number of disintegrations per unit of administered activity in the source organs was obtained by integration of organ activity curves for each xxx Tb-radioisotope, considering its physical half-life.These data were then used as input values in a human adult male phantom to perform dosimetric calculations with the OLINDA software code (version 2.2.3), as reported before [28], obtaining the organ absorbed doses per unit of administered activity.For each xxx Tb radioisotope, the specific effective dose, xxxTb ED corresponds to the sum of the product of the organ equivalent dose per unit of administered activity, xxxTb D org , and the respective tissue-weighting factor, w org , recommended by ICRP 103 [29], (2) The irradiation simulations with different target enrichment leads to a distinct distribution of xxx Tb activities, described by the fraction of total activity: The total effective dose (ED tot ) of Tb-cm09 was calculated at different times after the end of bombardment (EoB) by summing all the contributions of the produced Tb radioisotopes: Finally, the evaluation of the DI generated by the co-produced Tb-impurities is determined according to the following equation: and represents the ratio between the total effective dose and the effective dose due to an ideal injection of pure 155 Tb compound.The doses discussed above refer to the total doses of an injection delivered at time t post-production.

Assessment of the imaging properties of 155 Tb in the presence of other contaminant radioisotopes
The main problem in the quality of SPECT images when 155 Tb contains γ-emitting contaminants is the noise produced by the Compton scattering of high energy γ rays.Even using energy windows to acquire only the γ rays of interest, higher energy γ rays can be Compton scattered and registered by the acquisition system.The effect of such signals results in blurry and bad contrast images.Therefore, a study of the noise contribution of high-energy γ rays emitted by contaminant terbium radioisotopes has been carried out using a preclinical PET/SPECT/CT system (VECTor 5, Milabs) as a model of the acquisition system.
To simulate the γ-ray spectrum of each Tb radionuclide in VECTor 5, a homemade menu-driven spreadsheet software Visual Gamma was used.It makes use of specific initial system parameters such as efficiency and peak shape and simulates the spectrum of the radionuclides, using as input data the energy and the intensity of the γ ray, taken from the NUDAT3 database [24].Parameters such as γ-ray efficiency, the peak width and the shape of the Compton area as a function of the γ-rays energies have been determined using standard calibration point-like sources ( 57 Co, 22 Na, 60 Co and 137 Cs) in air.The process of assessment of the image quality, obtained using a particular γ-ray peak, was based on two main factors: the intensity of the γ-ray emission and the presence of Compton-scattered γ rays in the energy window of the peak selected for imaging.An energy window of 10% interval around the peak barycenter was chosen.To evaluate the influence of Compton-scattered γ rays on the final quality of the image, the ratio δ = N extr /N intr was used, where N extr is the number of all higher-energy γ rays that (5)

ED
underwent Compton scattering and fell into the selected energy window, and N intr is the number of γ rays from the selected photo-peak that fell into the energy window, i.e. not including Compton scattering contribution.It is convenient to use such Compton-topeak ratio since in the ideal case (lack of other noises) its value directly corresponds to the amount of the noises in the image.Imaging qualities were evaluated for 155 Tb produced from the proton irradiation of 155 Gd targets with enrichment of 100, 99 and 98%, at the EoB and 96 h later.

Cross sections
The same composition of the enriched 155 Gd targets employed by Dellepiane [10] for cross-section measurements ( 154 Gd 0.5%, 155 Gd 91.90%, 156 Gd 5.87%, 157 Gd 0.81%, 158 Gd 0.65%, 160 Gd 0.27%) was used in the irradiation simulations In Fig. 1 the experimental cross sections for 155 Tb production are compared with theoretical calculations.The gray band represents the models variability between Q1 and Q3, the dashed lines are the minimum/maximum values of the cross sections obtained with all models, and the blue solid line is the Talys default option (PE2-LD1), commonly referred to as the standard simulation in the literature.For the 155 Tb case all models reproduce equally the cross section up to 10 MeV, and for higher energy the band is relatively thin and this corresponds to a limited model variability, as expected from a typical (p,n) reaction.On account of this, and because the Talys default reproduces the very recent data measured by Dellepiane [10] quite satisfactorily, the analysis will be presented considering the Talys default as the benchmark calculation.
Figure 2 shows the contribution of each isotopic component of the target to the 155 Tb production.The main contribution to the 155 Tb cross section comes, as should be expected, from the 155 Gd component of the target.For energies higher than 10 MeV, also the contribution from the 156 Gd component becomes significant, and this explains the increase of the experimental data with respect to an ideal target with 100% 155 Gd enrichment.The contribution from 157 Gd is quite small and that from other Gadolinium components of the targets, such as 154 Gd, is negligible.This figure includes cross-section data for 100%-enriched targets assessed by Dellepiane et al., 2023 [30]; they were published after the submission of this study and added in the final revision of the text.
Next, the 156t Tb total cross section, which refers to the cumulative cross sections of 156g,156m1,156m2 Tb (ground and the first two metastable states), is presented in Fig. 3A.
The curve including all contributions is in agreement with the data measured by Dellepiane et al.It is evident that the main contribution belongs to the 156 Gd component of the target, followed by the 157 Gd component at slightly higher energies.At even higher energies, about 19 MeV, the contribution from 158 Gd also becomes significant.Instead, the contribution to the 156 Tb cross section from the main 155 Gd component of the target remains quite small in the entire range of considered energies.
Figure 3B, C and D describe the 156 Tb production cross sections for ground and the two metastable states, respectively.For the ground state, the model calculations are in fair agreement with the measurements, although at higher energies there is a tendency to underestimate the data.For the first metastable state there are no measured production data, while for the second one the curves underestimate the measurements.In all cases, by separating the individual isotopic components in the Gd target, it is evident that the main contributions to the production of the contaminants derive from the Gd isotopes heavier than 155.
The remaining contaminants that may have an impact on the quality of the produced 155 Tb are 153 Tb and 154 Tb (separated in ground and the first two metastable states) and their cross sections are given in Fig. 4. The theoretical curves are in clear agreement with the measurements, with the exception of 154m2 Tb where slight discrepancies can be observed.In all figures it is evident that the contamination at Fig. 2 155 Tb cross sections from enriched 155 Gd target: Talys calculations and contributions from different Gd isotopes of the target Fig. 3 156g, 156m1, 156m2, 156t Tb (where t in panel A is the cumulative cross section of 156g, 156m1, 156m2 Tb) cross sections from enriched 155 Gd target obtained from Talys calculations: contributions from different Gd isotopes of the target Fig. 4 153, 154g, 154m1,154m2 Tb cross sections from enriched 155 Gd target obtained from Talys calculations: contributions from different Gd isotopes of the target energies lower that 10 MeV can be entirely ascribed to the small presence of 154 Gd in the target, responsible of a characteristic bump seen in the lower energy data.At higher energies, 155 Gd becomes the principal contributor to the production of these contaminants.

Yields, isotopic, and radionuclidic purity
Starting from the cross section analysis of the effects of the isotopic components of the target used for 155 Tb production, yields, isotopic purity and RNP were assessed for different target compositions.Four different targets were considered in the simulations to evaluate the minimum enrichment required for a 155 Tb production with the purity needed for medical applications: the first target with the exact isotopic composition considered in measurements by Favaretto et al. [9] and Dellepiane et al. [10] (in particular with 91.90% of 155 Gd and 5.87% of 156 Gd), two highly-enriched targets (99% and 98% of 155 Gd and 156 Gd the only contaminant), and the ideal target with 100% of 155 Gd.The irradiation conditions, for all cases, were set to 1 µA current, 1 h irradiation time and 10.5-8 MeV energy interval, corresponding to the optimal energy selection defined by Dellepiane.
Table 2 reports the yields of the main Tb radionuclides involved in the production.The 156 Tb contamination grows proportionally with the fraction of 156 Gd in the target.At the energies considered, the production of 153 Tb is negligible in all cases.The 154g, 154m1, 154m2 Tb contamination remains small and stable when varying the 156 Gd component in the target.It increases only if the target contains a fraction of 154 Gd, as in the case of the target employed by Dellepiane, with a 0.5% contribution.Table 2 also exhibits the yields 72 h and 96 h after EoB, when the activities of 154g Tb and all metastable states are considerably reduced, since their half-lives are about 1 d or less.Figure 5 shows the time evolution of the 155 Tb RNP considering the four different target enrichment.The solid green line ( 156 Gd contamination 5.87%) levels at 93.5%, in agreement with the RNP measured [10] after 96 h from EoB. Significantly higher values are reached in the other three cases.In particular, 97.8, 98.8, and 99.8% RNP is obtained, after 96 h, with a target enrichment Gd directly affects the 155 Tb RNP, so it is crucial to limit it as much as possible.
Figure 6 shows the fraction of total activity of 154g Tb, 156g Tb, 156m1 Tb, and 156m2 Tb.Clearly, the amount of 156 Gd in the target proportionally influences the activity of all three states of 156 Tb.However, the main problem is represented by 156g Tb because its long half-life is comparable to the one of 155 Tb.Conversely, the production of 154g Tb represents a minor issue, because of its shorter half-life.Moreover, in the selected energy region, any residual amount of 156 Gd does not produce 154g Tb.This contaminant is mainly produced from 154 Gd, which appears, with a residue of 0.5%, only in the less enriched target.In addition, we have evaluated the isotopic purities and they correspond to 93.9, 97.9, 98.9, and 99.9% for the target enrichments of 91.90, 98, 99, and 100%, respectively, 96 h after the EoB..This implies that the production route is essentially carrier-free, without the presence of contaminants, including stable or long-lived ones.

Organ absorbed doses and effective dose due to XXX Tb-cm09 injection
Biokinetics curves were obtained plotting the radiopharmaceutical concentration corrected by the radioactive decay vs. time for each source organ of the ICRP 89 male phantom of 73 kg [31].The total volume of the blood (5110 ml) was obtained using the specific volume value of 70 ml/kg.Figure 7 shows a fast blood clearance followed by a quick radiopharmaceutical uptake by the main organs, with a slow wash-out.Liver and kidneys were the organs with slower clearance.
The number of disintegrations in the source organs, calculated by assuming that the injected radiopharmaceutical was labelled with only one of the radioisotopes 154g Tb, 155 Tb, 156g Tb, 156m1 Tb and 156m2 Tb, are reported in Table 3.These are the main radionuclides expected to be produced by proton irradiation of 155 Gd-enriched targets.The dosimetric properties of 154m1 Tb and 154m2 Tb were not assessed because these two metastable states are not included in the OLINDA software.However, at an irradiation energy of 10.5 MeV, 154m1 Tb and 154m2 Tb production is essentially due to the presence of 154 Gd in the target, not considered in the case of 155 Gd target enrichment ≥ 98%.Besides, since the energy of these two metastable states is very close to that of the ground state, their contribution to the absorbed dose occurs mainly through the decay of the ground state, properly taken into account through the application of the Bateman equations.For the dosimetric assessment, the total activity in the intestines was distributed in left colon, small intestine, right colon and rectum, according to their mass ratio with respect to the total intestine mass.The mean maximum volume of blood that can be contained in the four chambers of the heart, two atria and two ventricles, of an adult man is about 505 ml [32], therefore, just the 10% of the blood activity was assigned to the "heart contents" and the rest to "Remaining".The activity in mouse muscle was extrapolated to humans considering that the human muscle is 40% of total body weight.Muscle activity was also assigned to the"Remaining" organs when performing calculations with OLINDA 2.2.3, because this tissue is not included in the source organs of the ICRP 89 phantom model.
For each xxx Tb-radiopharmaceutical, the organs with the highest number of disintegrations were the kidneys, followed by the liver and the cortical bone.Comparing the different radioisotopes, it was found that the radionuclides with the highest number of disintegrations were 155 Tb and 156g Tb, due to their long half-life.156m2 Tb showed the lowest number of disintegrations because of its relatively short half-life.
The absorbed doses per unit of administered activity in the main human male organs, calculated for different xxx Tb-cm09, are reported in Table 4.For all the radioisotopes, the organs receiving the highest absorbed dose are the kidneys, followed by the osteogenic cells for 155 Tb, 156m1 Tb and 156m2 Tb, and by the adrenals for 154g Tb and 156g Tb. 156g Tb is the radioisotope giving the highest values of absorbed doses, due to its long half-life and the large amount of energy emitted per decay (see Table 1).The absorbed doses due to 154g Tb administration are also higher than those of 155 Tb, because, even if the 154g Tb half-life is much shorter than that of 155 Tb, the energy emitted per decay is about ten times higher.The absorbed doses due to 156m1 Tb and 156m2 Tb are about one order of magnitude lower than those of 155 Tb.The ED after administration of 156m1 Tb-cm09 and 156m2 Tb-cm09 are lower than the one of 155 Tb-cm09, consequently their presence as contaminants in the produced 155 Tb will not increase the dose to the patient.In contrast, when the radiopharmaceutical is labelled with 156g Tb or 154g Tb, the ED is 5.9 and 2.4 times the ED values of 155 Tb-cm09.Therefore, to guarantee a dose increment per unit of activity administered lower than 10%, the presence of 156g Tb or 154g Tb as contaminants must be lower than about 2% or 7%, respectively.

The dose increase due to 155 Tb contaminants obtained with different levels of target enrichment
The DI due to 155 Tb-cm09 obtained with different levels of target enrichment is plotted in Fig. 8 vs. time t post-production.The DI is well above 25% for the 91.9% target enrichment and this confirms that the isotopic contamination of that target is inadequate for medical purposes.The ideal case of a 155 Gd target with 100% enrichment is described by the solid red line.Here, the RNP is very close to 100% and the DI is negligible.Quite interesting are the dashed and dot-dashed curves obtained with 99% and 98% enrichment.The former reaches a 5% DI, while the latter remains within the 10% limit.

Impact on the image quality
The imaging quality of 155 Tb produced with the different 155 Gd-enriched targets was assessed at the EoB and 96 h later by considering the 156g Tb and 154g Tb yields reported in Table 2.The metastable states 156m1 Tb, 156m2 Tb and 154m1 Tb were disregarded in the simulations because of their lower energy γ-ray emission.Using the data in Table 2, a γ-ray spectrum produced at the EoB from the irradiation of 100% 155 Gd-enriched target was simulated, as shown in Fig. 9.The γ rays with the energies of 86.55 keV and 105.318 keV emitted from 155 Tb and the low intensity 88.97 keV γ ray emitted by 156g Tb are situated close with respect to the energy resolution of the imaging system and form a compound peak with a barycenter of 88.5 keV.Other four peaks from 155 Tb with energies of 148.64 keV, 161.29 keV, 163.28 keV and 180.08 keV form a less intense compound peak with a barycenter of 167 keV.A smaller 155 Tb peak at 262 keV (5.3%) is exposed to the higher energy γ rays from the isotopes of 154g Tb and 156g Tb only.For imaging purposes, the compound peak with the energy of 88.5 keV is most convenient, for its higher intensity, although the other two peaks at 167 keV and 262 keV can also be used for imaging.
Table 5 shows the calculated Compton-to-peak ratio for the three principal peaks of 155 Tb with respect to the γ-ray background, generated by the contaminant isotopes 156g Tb and 154g Tb.The analysis of the data reveals that in all cases the level of noise in the γ-ray images reconstructed using a 10% energy window around the barycenter of the peaks does not exceed 27%.For the intense peak at 88.5 keV the Compton-to-peak ratio remains within the interval of 19%-23% for all enriched targets immediately after the EoB and no significant improvement could be achieved 96 h later.Lower Comptonto-peak values, slightly below 20%, are obtained for the compound peak at 167 keV.This suggests its use as a single peak for image reconstruction in some particular cases.The lower intensity of the peak at 262 keV makes the ratio more affected by the presence of the contaminant nuclides 156g Tb and 154g Tb, although its maximum value still remains below 27%.However, it is worth noting that the present image quality estimation is made for point-like γ-ray sources in air, hence the values quoted refer to the imaging system only and do not take into account Compton scattering inside the tissues.To estimate this contribution for the case of small animals imaging, a water phantom was used in our laboratory, obtaining a noise increase level smaller than 5%.

Discussion
The use of low-energy proton beams on 155 Gd-enriched target represents a promising possibility for the production of 155 Tb, however the isotopic purity of the target strongly influences the amount of co-produced radioisotopes.The most problematic contaminant is 156g Tb, because its half-life is comparable to that of 155 Tb, and its γ emissions have severe impact on the dose released and on the quality of the SPECT images.For these reasons, its production should be kept as low as possible.The co-production of 156g Tb can be minimized by limiting the 156 Gd content in the target.In particular, a RNP higher than 98.5% can be obtained starting from about 30 h after the EoB with a 156 Gd content lower than 1%.For longer times, the RNP levels at 99%.However, the 99% value is not an established limit and a smaller one can be tolerated if the DI due to the presence of contaminant radioisotopes is low.By assuming a 10% limit in DI as an acceptable condition for the contamination of a production route [33], the maximum content of 156 Gd in the target could be increased to 2%.It is worth to note that in this case the DI does not exceed the 10% value in the entire time range shown in Fig. 8, namely from 0 to 30 days.Therefore the available activity can be utilized much earlier with respect to the timing when the RNP is close to 98%, and specifically, the product can be available right after the time needed to perform a radiochemical purification [9].The main limits of these results is that extreme care should always be taken when using biodistribution data from animals to predict absorbed doses to humans and that the DI due to contaminant radioisotopes must be determined for each specific radiopharmaceutical, because it directly depends on its biodistribution and kinetics characteristics [34].Taking into account these limitations, the dosimetric results obtained in this study should be interpreted as useful to have an idea of the contribution of Tb impurities to the absorbed doses imparted after injection of 155 Tb-cm09 and to select the 155 Gd target requirements for 155 Tb production, however a more precise estimation is possible only using biodistribution data collected on humans.
Enriched 155 Gd is currently commercially available from various companies (ISO-FLEX, CIL, AMT, etc.) with 155 Gd purity greater than 90%, but with a 156 Gd component larger than 2%, not yet sufficient for the production of 155 Tb with the necessary purity for medical application.Gd-isotopes separation is performed commercially using the "Calutron method", namely using a very large mass spectrometer for electro-magnetic separation.Since this technique is very energy inefficient, the prices for the enriched Gd material are high.However, it is possible to tune the Calutron to produce highly enriched 155 Gd with isotopic enrichment > 98.0%, or with 156 Gd content < 2%.To make a producer interested in changing the operational plan implies an economical return and profitability, which presently is not yet expected by the isotope-producer companies (Allan Pashkovski, Managing Director, Isoflex, Private communications, 6th July, 2023).A scenario of a diffuse worldwide production of 155 Tb by hospital cyclotrons may be much more favorable of economical return for the producer companies.
The targeted therapeutic agents currently used in nuclear medicine are based on biological structures such as peptides, monoclonal antibodies and their fragments, coupled to α and β -emitters radionuclides.These agents, due to their nature, exhibit maximum absorption time range from several hours to a few days after injection.To perform SPECT imaging of these radiocomplexes, a radiopharmaceutical labeled with a radionuclide having half-life comparable to this maximum absorption time is required [35].For this reason, 111 In is the only radionuclide currently used in the clinic to perform SPECT imaging of these radiocomplexes, since its half-life (2.8 d) is long enough for image acquisition even several days after the radiocomplex administration.Therefore, to evaluate the 155 Tb potential as matched pair for targeted theranostic agents labelled with 149 Tb and 161 Tb, it is interesting to compare its imaging and dosimetric properties with those of 111 In.
Imaging simulations were performed for pure 111 In using the main photon emissions at energies of 171 keV (91%) and 245 keV (94%) to compare its Compton-to-peak ratio with those of 155 Tb.The Compton-to-peak ratio of the peak at 171 keV, the closest in energy to the 88.5 keV and 167 keV peaks of 155 Tb, turned out to be 22.6%, which is slightly larger than the values obtained for 155 Tb, although they are comparable within 10% of uncertainty interval.These results are supported by the excellent imaging properties of 155 Tb reported by Favaretto [9], where a spatial resolution up to 1 mm was obtained from a Derenzo phantom.Moreover, Müller et al. [5] published a comparison between SPECT images of Derenzo phantoms filled with solutions of pure 155 Tb (2.6 MBq) and In (4 MBq), demonstrating comparable spatial resolution and the ability to produce equal images with lower activity of 155 Tb.
To evaluate the influence of the physical properties of the labeling radionuclide, the dosimetric properties of 111 In-labelled-cm09 were estimated with OLINDA by assuming the same biodistribution of Tb-cm09 compound.The calculated ED of 111 In-cm09 was 0.0217 mSv/MBq, very similar to that of 155 Tb-cm09 (0.0186 mSv/MBq).These results are a consequence of the fact that, despite the half-life of 155 Tb is almost 2 times longer than that of 111 In, its total energy emission per decay is about 2 times lower compared to 111 In (0.4409 MeV/nt).It should therefore be expected that the imaging and dosimetric properties of 155 Tb-radiocomplexes are comparable and even better than those of 111 In ones.

Conclusions 155
Tb can be produced with a quality suitable for medical applications using low-energy proton beams and 155 Gd-enriched targets, if the 156 Gd impurity content does not exceed 2%.Under these conditions, the dose increase due to the presence of contaminant radioisotopes remains below the 10% limit and good quality images, comparable to those of

Fig. 1
Fig. 1 Experimental 155 Tb cross sections from enriched 155 Gd target and theoretical curves expressing the variability of nuclear reaction models

Fig. 5
Fig. 5 Time-evolution of 155 Tb radionuclidic purity for the different target enrichments.The curves start 1 min after EoB to eliminate the transient effects of the rapidly decaying products

Fig. 7
Fig. 7 Time-activity curves of Tb-cm09 in the main source organs: symbols show experimental data obtained from biodistribution studies; the lines depict fitted time-activity curves

Fig. 8 Fig. 9
Fig. 8 Dose increase for 155 Tb-cm09 radiopharmaceutical labelled at different times after the end of irradiation of different 155 Gd-enriched targets

Table 1
[24] decay data of xxx Tb radionuclides.mode,half-life,emittedenergy per nuclear transformation (nt) in the form of electron, photon, and the total one of xxx Tb radionuclides were obtained by the ICRP 107 publication[23]and NUDAT3 database[24](emitted energy data not available for 154m1 Tb and 154m2 Tb).

enrichment 155 Tb 156g Tb 156m1 Tb 156m2 Tb 154g Tb 154m1 Tb 154m2 Tb
of 98, 99, and 100%, respectively.It is evident that the contamination of the target with 156

Table 3
Number of nuclear transitions (MBq × h/MBq) in source organs per unit administered activity of xxx Tb-cm09 for male ICRP 89 human phantom

Table 4
Organ absorbed doses (mGy/MBq) and ED values (mSv/MBq) per unit administered activity calculated for xxx Tb-cm09 for male ICRP 89 phantoms

Table 5
Compton-to-peak ratio calculated from the simulation of 155 Tb γ-ray spectra, obtained for 100%, 99% and 98% enrichment of155Gd targets, at the EoB and 96 h later