Table of Contents Author Guidelines Submit a Manuscript
Science and Technology of Nuclear Installations
Volume 2013, Article ID 972381, 14 pages
http://dx.doi.org/10.1155/2013/972381
Research Article

Evaluation of and Productions Based on a High-Performance Cyclotron

1INFN, Laboratori Nazionali di Legnaro (LNL), Via dell’Università 2, 35020 Legnaro, Italy
2Dipartimento di Fisica, Università di Ferrara and INFN, Sezione di Ferrara, Via Saragat 1, 44122 Ferrara, Italy
3Dipartimento di Morfologia, Chirurgia e Medicina Sperimentale, Università di Ferrara e INFN, Sezione di Ferrara, Via Luigi Borsari 46, 44121 Ferrara, Italy

Received 3 June 2013; Revised 26 August 2013; Accepted 27 August 2013

Academic Editor: George Vandegrift

Copyright © 2013 J. Esposito et al. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

Abstract

Following preliminary feasibility studies which started at Legnaro National Laboratories (LNL) in 2011, the Italian National Institute for Nuclear Physics (INFN) research activities are underway aiming at the alternative, accelerator-driven, production routes. One of the most promising approaches is to use 100Mo-enriched (i.e., >99%) molybdenum metallic targets, bombarded with high-beam-current, high-energy proton cyclotrons. In order to get a comprehensive map of radionuclides expected, a detailed theoretical investigation has been carried out using the TALYS-TENDL 2012 excitation functions extended up to (p,6n), (p,p5n), and (p,2p4n) levels. A series of quality parameters have thus been calculated both at the end of beam (EOB) and at longer times. Results point out that accelerator-99Mo is of limited interest for a possible massive production because of the quite low specific activity with respect to reactor-99Mo. Accelerator- quality parameters (i.e., radionuclidic purity (RNP), isotopic purity (IP), and specific activities) calculated are instead quite close to the generator-Tc. Calculations at 15, 20, and 25 MeV have thus been performed to assess the best operative irradiation condition for production while minimizing both the short-lived and long-lived Tc contaminant radionuclides. Although present in minimum quantities, Tc contaminants may indeed have an impact either on the pharmaceutical labeling procedures or on contributing to patient radiation dose during the diagnostic procedures.

1. Introduction

Nearly 80 of the radioisotopes used in nuclear medicine are currently produced by nuclear reactors. Among them, (6.01 h half-life), a decay product of (65.94 h half-life), is the most important and widely used one. However, such a vital radionuclide is mainly supplied by two production facilities only, which cover about 65 of the world demand: the research reactor NRU at Chalk River, Ontario (Canada), and the HFR reactor at Petten (The Netherlands). The current mass production is indeed based upon the isotope separation from fission products by using the highly enriched uranium, weapons-grade (HEU-WG) material (235U enrichment 80). Such a production is therefore subject to the strict international regulations and control actions against the proliferation of nuclear weapons. Both reactors have experienced in recent years (2009-2010) some long-scheduled/unscheduled shutdowns, which caused a temporary shortage of radionuclides on the international market in 2010. Such an event has thus pushed new ideas about alternative arrangements, all based on accelerators, as well as making use of nonstrategic materials.

Since the early 1970s, production through accelerator systems was demonstrated feasible, as proposed first by Beaver and Hupf [1], mainly by the 100Mo(p,xn) reaction route. However, both excitation functions and production yields, due to the different Tc isotopes and isomeric states generated, still need to be made clearer. That is because of an unexpected spread in results, obtained over the past 40 years in different experimental campaigns. Results provided in the last 20 years, using nonstandard experimental procedures, and reported in different works [213] may give an idea about that and explain why further experimental investigations are demanded.

Moreover, labelling efficiency and biodistribution studies of some interesting radiopharmaceuticals have been recently carried out, using both cyclotron- and generator-produced , as reported by Urbano et al. [14], Guérin et al. [15], and Targholizadeh et al. [16]. Such encouraging results have driven the first study aimed at the assessment of irradiation conditions for accelerator-based Tc production, as discussed by Celler et al. [17]. Based on the excitation functions from theoretical nuclear models, both the production yields and a reduced map for other radioactive and long-lived generated isotopes (considering some simplified assumptions and constraints) have been determined.

In the framework of the (Laboratory of Radionuclides for Medicine) (LARAMED) project, funded by the Istituto Nazionale di Fisica Nucleare (INFN), a possible future supply of a set of radioisotopes in Italy, aimed at both research and medical application purposes, has been taken into account. This project derives from the installation at LNL in the near future of a high-performance cyclotron (proton output energy tunable in the range 35–70 MeV, maximum beam current 500 A).

Taking into account the wide range of irradiation conditions which will be provided by such an accelerator, a theoretical assessment study has been preliminary carried out for both and productions. Assuming a -enriched molybdenum metallic target, as currently available on the isotopes and pure materials market, a beam power areal density of 500 W/cm2 and a beam current up to 500 A, a series of parameters have been calculated for.(i)Optimal Molybdenum Layer Thickness. Two basic configurations have been taken into account; in option 1, the incident beam energy down to the threshold for the reaction is considered, thus avoiding the useless heating due to the Bragg peak (i.e., thick target configuration). In option 2, a recommended thickness based on the estimation of yield distribution versus beam penetration depth is instead analyzed. In such a way, an optimized radionuclide production (i.e., avoiding the drop-off region because of cross section lowering) is achieved, with the minimum required thickness for the isotope-enriched material (not inherently thick target configuration).(ii)Different Proton Energies. Taking advantage from the cyclotron directly tunable beam output energy, the reference cases at 40 and 70 MeV have been taken into account for production. direct production at 15, 20, and 25 MeV has instead been investigated, in order to limit the amount of contaminant isotopes. Supposing, as reference, an energy degrader device made of 3.67 mm carbon layer, beam energy may be shifted from 35 MeV (i.e., the minimum cyclotron output energy) down to 20 MeV. In such a case, an approximate Gaussian energy spread () of about 1.25 MeV occurs. However, 90 of beam current (i.e., 450 A) is within an energy window of just 550 KeV which, through a Wien filter (if necessary), may be driven towards the target station at the accelerator beam exit. (iii)Different Irradiation Times. 12 h, 24 h, and full saturation ( d) for production, while several short irradiation times up to one half-life (i.e., 1 h, 2 h, 3 h, and 6 h) have been taken into consideration in case of production. Irradiation times longer than are not convenient from the point of view of either the product final cost or the accumulation of useless longer-lived isotopes.

The isotopic composition of the enriched metallic molybdenum material, taken as reference in this study, is the one provided by the ISOFLEX company [18]: (99.05), (0.54), (0.07), (0.11), (0.10), (0.05), and (0.08). The reaction yields for each Mo isotope have been estimated from the theoretical excitation functions leading to nuclei both in isomeric and ground states. For all nuclides, either directly produced following different reaction routes or generated as a decay product in any of the decay chains, the cumulative yields have been calculated.

Considering the direct production, particular attention must indeed be paid to Tc contaminant radionuclides by reaction routes opened up by all Mo isotopes because of the same chemical species. Although present in rather low quantities, their contributions (except ) are however not negligible in the overall yield estimation of accelerator-produced Tc. Optimal irradiation conditions (i.e., energy intervals, beam current, and irradiation time) are therefore needed in order to minimize their production. Finally, the long-lived , , and , having large cross sections and mean lives longer than 105 years, need to be produced in as low quantities as possible, having a direct impact on the specific activity of the final product and thus on the labelling procedures of the radiopharmaceutical preparation. Care has also to be taken of the productions of shorter-lived Tc isotopes (mainly Tc) because of their contribution to extra radiation dose. However, such an issue is strictly related to the isotopic material composition made available by the supplier, rather than an optimal combination of irradiation time-beam energy parameters.

2. 99Mo and Production Routes

The mass production using the commercial proton-driven accelerators currently available, or next to be put into operation (as the one which will be installed at LNL), is known to be hardly achieved through neutron reaction routes, that is, using either the spallation reactions on heavy target systems or the relatively low energy beams driven by high intensity accelerators on light targets, as argued by Froment et al. [19] and Abbas et al. [20]. The radiative capture reaction (n,) on -enriched samples would provide production with good purity level. However, the resulting specific activities estimated and experimented are too low for a feasible massive production, even if 90 of neutrons are slowed down in the resonance energy interval (1–100 eV) [19, 20]. Moreover, exploiting the other possible (n,2n) inelastic scattering reaction on -enriched samples, the limited neutron flux makes the expected production level poor. Some works are under way, as reported by Nagai and Hatsukawa [21] and Minato [22], although the approach followed to work out the low specific activity is to get the through a thermal-induced sublimation of technetium from the irradiated molybdenum bulk targets.

On the other hand, the 235U(n,f) fission reactions, that is, the subcritical version of the standard reactor-based route using accelerator-driven neutron multiplier systems, have also to be considered as long-term solutions because of the lack of information on production yields and economics, as discussed in NEA report [23]. For such reasons, only the proton-induced reactions (p,pn) and (p,2n) are taken into account in the present study, being the best short-term solution to the possible future massive supply. The current cyclotrons have indeed the required beam current to ensure a high production level of these vital isotopes.

2.1. Production Route: The (p,pn) Reaction on

production may derive from two reaction routes: the main (p,pn) reaction and as a decay product of ( s) and ( m), due to the additional (p,2p) reaction, when operating at proton beam energies above 10 MeV. Figure 1 shows a collection of the (p,pn) experimental excitation functions measured in the last two decades.

972381.fig.001
Figure 1: The collection of (p,pn) excitation functions experimentally measured in the last two decades. Error bars are shown for the most recent measurement only performed by Tárkányi (2012) [13]. The most recent theoretical excitation function by the TENDL 2012 library is also shown [24].

As can be observed, results are consistent up to 25 MeV, while at higher proton energies an unusual spread (i.e., uncertainty band) up to 100 mb (including error bars extension) may be noted among results obtained in different experimental campaigns. Moreover, even considering % of -enriched samples, additional reaction routes for the production of different Tc, Mo, and Nb isotopes (including Zr in case of a real ≥99% -enriched material) are already opened at proton energies higher than 10–15 MeV. The theoretical excitation functions available from the TENDL 2012 library [24] about the other Mo and Nb isotopes, yielded by other open reaction routes on , are shown in Figures 2 and 3, up to the (p,p5n) and (p,2p4n) levels, respectively. When available, experimental measurements are plotted as a benchmark.

972381.fig.002
Figure 2: The collection of (p,pxn) theoretical excitation functions from the TENDL 2012 library [24] up to the (p,p5n) level. The experimental measurements of by Levkovskij [2] and Takacs et al. [5] only are included for comparison.
972381.fig.003
Figure 3: The collection of (p,2pxn) theoretical excitation functions from the TENDL 2012 library [24] up to the (p,2p4n) level. The available experimental measurements by Levkovskij [2] for , , and nuclides are included for comparison.

Although predictions about the excitation functions for some reaction routes still need to be experimentally validated in the entire energy range, we are aware that a production of these nuclides is indeed energetically feasible. Nb and Tc isotopes expected to be yielded by open reaction routes on isotope only, are listed in Table 1 as reference, including decay modes. Similar considerations have been performed for the remaining Mo isotopes included in the enriched molybdenum material as well, in order to determine the whole list of nuclides (both stable and radioactive) generated by open reaction routes. The long-lived radioisotopes created, like and , , , , , and , are however not considered to be a concern in the following chemical processing aimed at the preparation of generator. The Nb, Tc, and even Rb elements may indeed be easily chemically separated from Mo for later reuse.

tab1
Table 1: Niobium and technetium isotopes expected to be produced by the Mo(p, x) reactions.

2.2. Production Route: The (p,2n) Reaction on Mo

The excitation function (p,2n) has been repeatedly measured over the past 40 years in different experimental campaigns (part of them shown in Figure 4), proving that within the energy range from 5 up to 70 MeV a single peak is present, centered around 15 MeV and estimated to be about 300 mb by Levkovskij [1] twenty years ago. The measurements repeated later by Scholten [4], Takacs et al. [5], and Khandaker et al. [7] have reduced such a value to about 200 mb. More recent evaluations point out higher peak values, about 250 mb, as reported by Lebeda and Pruszyński [9] and Tárkányi [13], later increased up to 300 mb by Gagnon et al. [12]. Even higher cross section peak values (i.e. >350 mb), originally measured by Lagunas-Solar et al. [11], more recently repeated by Challan et al. [8] are now considered not credible. Analyzing the different works performed, many reasons may be addressed to meet such an unusual spread observed, without a clear trend (i.e., different experimental set-ups, target material purity, contaminant levels, isotopic compositions, and ad hoc procedures followed during sample irradiations and detection methods). Moreover, different correction methods are applied in calculating the activity, as deduced by the 140 keV photon peak emissions through gamma spectrometry measurements.

972381.fig.004
Figure 4: The collection of (p,2n) excitation functions experimentally measured in the last two decades. Error bars are shown for the most recent measurement only performed by Tárkányi (2012) [13]. The most recent theoretical excitation function by the TENDL 2012 library is also shown [24].

Considering such a production route, particular care has indeed to be taken of the Tc contaminants (see in Table 1 those yielded by ) since they would remain in the final Tc product. In order to minimize them, irradiation energies at 15, 20, and 25 MeV only have been investigated in this work. In Figure 5, the theoretical excitation functions from the TENDL 2012 library [24], concerning the 100Mo(p,X) reaction routes for additional Tc isotopes production, both in ground and isomeric states, are plotted up to the (p,6n) level. The experimental excitation functions measured for by Takacs et al. [5] and for by Skakun [25] are shown for comparison as well. In Figure 5, one can notice that only the long-lived contaminants , , and will be produced by isotope up to 25 MeV.

972381.fig.005
Figure 5: The collection of (p,pxn) theoretical excitation functions from the TENDL 2012 library [24] up to the (p,6n) level. The experimental measurements for (Gagnon et al. [12]), (Skakun et al. [25]), and (Takacs et al. [5]) productions are included for comparison.

Figure 5 also shows the first experimental evaluation of the (p,2n) excitation function, performed in the energy range 8–18 MeV and reported by Gagnon et al. [12], from which a general agreement may be found with the theoretical one. Being the largest Tc contaminant produced, further studies and experimental campaigns are however requested, mainly to validate the excitation function in the entire energy range of interest.

Moreover, an additional issue to be assessed is the necessary reuse of the -enriched material (quite expensive) in case of direct production. The Mo isotopic composition of the original enriched target material is indeed expected to be modified, if multiple reuse is planned due to economic reasons. In this regard, a preliminary investigation has been carried out on irradiated -enriched metal samples, but considering only one recycling stage, as reported by Gagnon et al. [26]. Such a test basically demonstrated the feasibility of molybdenum recovery from the chemical point of view but is however not able to provide a final answer to such a question. The issue about how many times the molybdenum material recovered might be reused, before an impact on the accelerator produced Tc quality becomes critical, still needs to be defined.

3. Materials and Methods

3.1. Decay Chain Calculations

When molybdenum targets are irradiated, some of the reaction routes lead to the production of different radioactive isotopes, thus resulting in different radioactive decay chains. Since the same radionuclide may be created through a number of production-decay routes, all of them must be accounted for when estimating the final number of atoms and activity (relative to radioactive ones) available inside the target. The total number of any atom species present at time , , and the activity may be obtained by summing up values from the contributing chains. The general production-decay scheme of a single radioactive chain is therefore the following: where are the number of radioactive products of a given chain, is the corresponding decay constants, is the number of the last (stable) nuclide, and is the yield rate for the formation of the first nuclide of each chain directly produced. Since the (p,X) cross sections are relatively low (hundred mb as order of magnitude), only nuclear reactions on target nuclei were considered in calculations. Any secondary reaction that may occur at each level of the newly produced isotopes has been ignored. Likewise, the change in the number of target nuclei during the irradiation process has been considered negligible (i.e., are time constants).

For each decay chain, the variation in the number of nuclei versus time is described by the general set of differential equations. In case of more than one decay mode of a given radionuclide, every branch of the chain has been considered as a different, independent chain according to the branching ratio . Consider with boundary conditions for each equation given by

Such initial conditions are the most general possible, since at starting irradiation time , a certain amount of each (radioactive) nuclide species might be present inside the target material (i.e., in case of multiple reuse of target material). During calculation runs, however, all initial amounts for radioactive species have been set to zero, since only stable molybdenum isotopes were supposed to be included inside the target material, and a unique irradiation run was considered.

In the first one of (2), is the number of nuclei of a given species yielded as a result of the nuclear reaction. If such a nuclide species is radioactive, the number of nuclei available at any time is thus a balance between the production rate and the radioactive decay with constant . The solution for can be easily obtained by multiplying both sides of the equation by and then integrating the following which leads to

On the other hand, if the nuclei are stable, then increases linearly with time and is given by

Likewise, the remaining equations may be analytically solved by direct integration, once the solution for the previous nuclide of the chain has been determined. Consider the following for the generic and nuclides.

After analyzing both stable and unstable solutions for decay chains as long as 6 elements, two general and compact solutions have been obtained, which represent the number of nuclei of, respectively, the unstable and the stable elements of a given decay chains. Consider

The estimation of every nuclide produced by any decay chain has been thus calculated at any given irradiation time. For each decay chain, the nuclide number and the related radioactivity have been calculated, adding up the different contributions to a given nuclide species. Once the irradiation time is over (), the nuclide yields by nuclear reactions stop, while the decays still continue. All for each decay chain are therefore set to zero, and the same set of equations may be solved to calculate both the nuclide number and activity change against the decay time .

3.2. The Estimations of Mo and Yields Expected inside Irradiated Samples

As a starting approach, the 99Mo production level expected at EOB on 100Mo-enriched molybdenum metallic targets has been estimated at 40 MeV and 70 MeV proton energies, while for the direct production 15, 20, and 25 MeV have been instead considered. In both cases, the maximum of 500 A beam current has been assumed as a reference. Yields expected at different beam currents anyway may be easily inferred through simple scaling calculations.

For the thick-target configuration, the required thickness for the molybdenum material has to be primarily assessed both at the minimum and the maximum proton energies, in order to limit the amount of 100Mo and, consequently, the material cost. This can be done assuming, as a starting step, that the proton beam hits the molybdenum sample in normal direction and using tabulated data on stopping power values and projected ranges calculated by SRIM 2012 code [27] or taken by the PSTAR database of the National Institute for Standards and Technology (NIST) website (NIST-STAR database [28]). The estimated metallic molybdenum thickness to fully stop 70 MeV protons turns out to be about 7.1 mm, while 2.6 mm is instead needed when using 40 MeV protons. When the incident energy optimized for production is instead considered 0.48 mm, 0.78 mm, and 1.1 mm are enough to fully stop 15, 20, and 25 MeV proton beams, respectively.

Such penetration thicknesses for protons are actually needed also in other irradiation conditions (tilted targets). The proton energy loss per unit path, as well as the corresponding energy drop off versus penetration depth, has been calculated. Since the sample thickness required to fully stop the beam is much lower than the maximum linear size (radius) of the target material in the hypotheses of cylindrical configuration, an analytical calculation approach, based on a slab geometry model, may be used in such a case with good approximation (see Figure 6).

972381.fig.006
Figure 6: The main parameters in the slab approximation geometry for the calculation of yield distributions inside molybdenum sample thickness.

The estimation of local yield contribution , that is, the nuclide species produced inside the given infinitesimal thickness at depth of the target material, is described by the following equation: where is the atomic density of the enriched molybdenum target material (i.e., ) with the target density, the weight fraction of the isotope considered, the Avogadro number, and the molar mass of the target. Then, is the (p,x) cross section for the production of the given nuclide species, as a function of the proton energy at each thickness (estimated by an iterative calculation process once known the proton stopping power and the incident proton energy ). The proton beam current divided by the electric charge unit is the number of protons per unit time, hitting the target ().

The product expression within brackets is the fraction of incident protons, available at the generic depth inside the sample, because of the contribution of the exponential attenuations in any single infinitesimal slab. The removal cross section takes into account all the reaction routes that remove protons from the main stream.

The resulting contribution of the exponential, for each infinitesimal thickness inside samples, however turns out to be basically constant and quite close to unity, mainly because of the order of magnitude of the cross section involved (i.e., 10–100 mb). Expression (9) may therefore be simplified to a very good approximation in the next one. The overall nuclide production , normalized for incident proton, is obtained by integrating (9) over the whole target thickness and dividing by as follows:

Equation (10) was used for later calculations of the in-target production of any Mo, Nb, Zr and Tc nuclides. The theoretical excitation functions for the reaction routes of interest, involving up to 6 emitted particles (i.e., protons and/or neutrons), have been taken into account as provided by the TENDL 2012 library and available as output by TALYS code. The excitation functions for both isomeric and ground states may therefore be calculated taking into account a complete set of nuclear models developed (i.e., optical model, direct reactions, compound reactions, etc.). Further information may be found in [24].

4. Results and Discussion

and yields versus sample thicknesses calculated, normalized per incident proton, are plotted in Figure 7 for the energies concerned. As can be noticed, values are within one order of tenth, up to the point where the yield drop off starts down to the corresponding energy threshold for the reaction (i.e., energy exit without the useless heat power deposition corresponding to the Bragg peak). Figure 7 also shows the recommended thicknesses for an optimized production, saving the target material.

fig7
Figure 7: The distributions of the nuclide yields, normalized per incident proton, versus penetration depth, calculated inside 99.05  -enriched samples: (a) 70 MeV and 40 MeV are considered for production; (b) the same plot for 15, 20, and 25 MeV beam, in case of direct production. The eye-guide vertical lines show the recommended sample thicknesses to get an optimized production.

In Table 2, the isotope productions expected by both thick and optimized target thicknesses, normalized per incident proton, are reported for the selected proton energies. As can be seen, the accelerator production of is interesting at energies larger than 40 MeV, in order to operate in the most useful energy range for the the excitation function (see Figure 2).

tab2
Table 2: The integral production yields (normalized per incident proton) estimated on Mo-enriched sample thicknesses both for 99Mo and Tc productions. Additional production quality parameters for Tc are listed as well.

As for the direct production, the quality parameters (i.e., the and the nuclei ratios) have been calculated and listed. Such parameters, which directly relate the production quality expected to the proton beam energy selected, have to be as high as possible. Almost all of the Tc nuclides directly yielded at 15 MeV, caused by and mainly, with a contribution of about 10 from the short-lived , are shown in Table 2, in case of the optimized target configuration. At 20 and 25 MeV, significant contributions are instead from other Tc contaminants (i.e., mainly from ,   , and at lesser extent by   ) produced by the reaction routes due to the other Mo target isotopes, which basically cause a decrease of the ratio. The contribution from Tc nuclides other than , compared with , indeed increases to about 17 and 37, respectively.

On the basis of these considerations, it may be inferred that the best irradiation energy which should be considered for the accelerator production is 15 MeV. However, the production rate is about 2.0 times lower at this energy than expected at 20 MeV. A weighted balance, between the yield ratio and production yield expected, is therefore needed at such an energy range. The production at 25 MeV should instead be avoided because of the low ratio (basically the half with respect to 15 MeV) which is available just from direct reactions only. It has anyway to be kept in mind that once produced, the compositions of different Tc nuclides are continuously modified by the decay chain contributions during the irradiation time, up to the EOB and for times longer, as discussed later in a more detailed way. In other words, the yield ratio numbers listed in Table 2 are the upper values for the isotopic purity (IP) parameter that may be achieved in the theoretical approach concerned.

Table 3 reports the 99Mo in-target yields estimated at EOB for different irradiation conditions, taking into account the contribution from the decay of parent nuclides . The proton beam spot size and energy, resulting from the mean areal power density, uniquely defines the size and mass of the requested samples (thick or optimized), thus allowing to estimate the specific activity as well.

tab3
Table 3: The Mo production yields expected at EOB for different irradiation conditions, based on 99.05Mo-enriched metallic molybdenum (thick and optimized target configurations), 500 A proton current, and 500 W/cm2 mean areal power density on target.

Table 4 reports the production yields estimated at EOB at 15 and 20 MeV for the optimized sample configuration at the irradiation times as long as (1 h, 2 h, 3 h, and 6 h). The specific activity (in isotopic carrier condition) is estimated as well. Moreover, all technetium radionuclides, as well as activity, versus activities of all radioactive species, are calculated. The specific activity, the radionuclidic purity (i.e., ), and the isotopic purity (i.e., ) expected at EOB are also listed, taking into account the contributions from all the radioactive chains available.

tab4
Table 4: The Tc production yields estimated at EOB for different irradiation times at 15 and 20 MeV proton beams on 99.05Mo-enriched metallic molybdenum and optimized target configuration. 500 A proton current and 500 W/cm2 mean areal power density are considered on target. A series of production quality parameters are calculated and listed (see body text).

As may be observed, almost all of the in-target radioactivity produced at EOB is due to Tc radioisotopes (i.e., about 98 at 15 MeV and 95 at 20 MeV, resp.). The radionuclidic purity (RNP) is however relatively low (around 20–30%), even considering the best irradiation combination (e.g., 15 MeV proton beam and the short irradiation times in order to get the highest IP). The reason for that is basically the contribution of the quite short-lived 100Tc, the presence of which may not be avoided. At both energies considered, the 100Tc activity contribution ranges from 5 times larger to values almost equal to that provided by , at increasing irradiation times.

As a reference, in Table 5, the calculated integral yields (mCi/A) at 3 and 6 hours irradiation, from all of the Tc radioactive products expected at EOB inside the thin target configuration, are listed at the three bombarding energies of 15, 20, and 25 MeV.

tab5
Table 5: The calculated integral yields (mCi/A) at EOB after 3 and 6 hours irradiations for the expected Tc radioactive species inside the target. 99.05Mo-enriched metallic molybdenum in the optimized target configuration irradiated at 15, 20, and 25 MeV is considered. The whole list of all Mo, Nb, and Zr radioactive yields is not reported here for the sake of brevity.

Based on all the calculations performed, the resulting activity expected after 24 h irradiation inside sample at EOB is 109.6 Ci, considering an optimized target configuration under 70 MeV beam, while the specific activity is being 256.9 mCi/g. When 40 MeV protons are instead taken into account, 32.6 Ci activity is estimated at EOB, with a specific activity increase to 386.0 mCi/g.

Although these in-target activity levels seem to be enough to cover a regional demand, it should be however noted that the resulting specific activities are a factor of 104 lower than the one present inside the current Mo/Tc generators, achieving levels of about 104 Ci/g. If the current industrial method for Mo/Tc generator manufacturing has to be maintained, the accelerator production is therefore not a favourable option of practical interest, considering the high cost of the target material, the very large alumina column needed, and the resulting large elution volumes, as discussed by Morley et al. [29].

When the direct production is instead taken into account in the range 15–25 MeV, the estimation of the resulting specific activity (see Table 4) is similar to the one provided by the current Mo/Tc generators, which turns out to be about 1.5  Ci/g at any standard 24 h elution. Table 5 shows that the major contribution to the Tc-produced activity at EOB is from the short-lived and the production energy range of which is only partly included in the optimized target configuration. Although the exit proton energy, at the indicated target thicknesses reported in Table 2, turns out to be around 9.6 MeV, its production may not be avoided. Because of the quite short mean life, does not play a disturbing role in the final Tc contaminants. Just a few minutes after EOB, the activity indeed drops off, transmuting in the stable which may be later separated in the chemical process of the target dissolution. The contributions from other main Tc contaminants to the overall Tc activity which, instead, remain in the final Tc solution are, in order of decreasing activity (reference case at 20 MeV protons and 3 hrs irradiation), , , , , , , and . On the other hand, the activity contribution from the long-lived , , and , which are produced in larger amounts because of the highest values of the excitation functions, is minimal and therefore irrelevant for dosimetry considerations. Nevertheless, their impact on the isotopic purity level is not negligible and may affect the radiochemical quality of the accelerator-Tc labelled pharmaceuticals.

As recalled, a key parameter to be carefully assessed considering the accelerator production is the ground state , which is useless for the diagnostic procedures. In the report by the european commission [30], a limiting purity of the end product, approximately composed by 25%   and 75% (i.e., ratio equal to 0.25), is reported to interfere with the function of some labelled radiopharmaceuticals, thus reducing the effectiveness of Tc-based scans. Unfortunately, it is not clear in [30] how such an IP level should affect the diagnostic procedures, as neither the radiolabelling processes nor the resulting SPECT images quality are mentioned.

In order to get a reference, when technetium (), under pertechnetate form (), is eluted as a sterile solution from a standard Mo/Tc generator after intervals of 24 h, the maximum activity is reached and the corresponding isotopic purity IP turns out to be about 0.26, as may be inferred by the considerations reported by Alfassi et al. [31]. On the other hand, on the basis of our theoretical approach it can be noted that such a ratio is never higher than 0.22 (see Table 4). It has to be reminded that the excitation function from the TENDL 2012 library shows the peak located at the bottom of the uncertainty band from the experimental measurements, as shown in Figure 4. Therefore, it has to be reasonably expected that an IP value similar to Tc generators might be obtained at both short irradiation times and low bombarding energies (i.e., 1-2 hr, 15 MeV).

Once produced, and eventually chemically separated from other radioactive and stable isotopes of other chemical species, the expected evolution of both the isotopic purity (IP) and the radionuclidic purity (RNP) for the accelerator-produced , against the decay time after EOB, is finally shown in the three plots of Figure 8. The sharp increase of RNP, with respect to the ones enlisted in Table 4 up to values higher than 95, occurs in just a few minutes, mainly because of the quick decay of . The decays of other short-lived Tc isotopes lead the RNP increase to values as high as 99 in about 1 hour after EOB, considering the reference case of 3 h irradiation at 20 MeV energy. RNP values above 99 for accelerator-produced Tc are quite important, considering that the same parameter from the generator-produced Tc is about 99.99. Such a limit is basically approached (i.e., 99.58) if the irradiations at proton energies as low as 15 MeV are performed, whatever the irradiation time chosen. Moreover, the accelerator produced , having RNP values never higher than 99 and IP ones always lower than 10 1 h after EOB, clearly points out that irradiations at 25 MeV have to be avoided, even using highly enriched molybdenum material.

fig8
Figure 8: The evolution of isotopic purity (IP) and radionuclidic purity (RNP) expected inside the target (or supposing a quick chemical separation process of accelerator Tc produced) versus the decay time after EOB. 99.05  -enriched metallic molybdenum in optimized sample configuration bombarded at 15, 20, and 25 MeV protons and 1, 3, and 6 h irradiation times has been considered.

Some experimental investigations have been performed in recent times to assess the influence of the IP parameter on radiolabelling procedures, using from standard Mo/Tc eluates at times as long as after 3 days from the previous elution. In a recent study performed by Urbano et al. [14], the and oxidizing impurities amounts present in different eluates were accurately determined and their influence in radiochemical purity determination of different radiopharmaceuticals, for example, -DTPA (diethylenetriaminepenta-acetic acid), -MIBI (methoxyisobutylisonitrile), and -HMPAO (hexamethylpropylene amine oxime), was carefully evaluated. It was found that the radiolabelling procedures of these commercial kits were not affected when using eluates at 72 h and at longer times from the previous one (i.e., solutions with an IP value as low as about 0.08). Moreover, the USP requirements were fulfilled in all the cases considered, even with those radiopharmaceuticals whose stannous contents were low (i.e., MIBI and HMAPO). Such a consideration suggests that, as shown in Figure 8, the useful time-window (1–10 h) after EOB for labelling procedures may be fulfilled by accelerator-produced Tc at energies below 20 MeV and irradiation times within 3 h.

It has to be noted that some interesting imaging in vivo tests have also been performed on healthy rats, to get a comparison between the generator-produced against the cyclotron-produced , as reported in a newsline of JNM by Guérin et al. [15]. After short irradiations (1.5–3 h) at 17 MeV using 99.5% 100Mo-enrichment, the resulting labeling, the resulting labelling efficiencies of -MDP and -MIBI were above the USP requirements (>90%), and identical patterns were found for the cyclotron- and the generator-produced radiopharmaceuticals, within normal individual variations between each pair of animals.

Moreover, in vivo biodistribution studies in rats were performed, using the cyclotron-produced -BRIDA pharmaceutical, whose radiochemical purity after the technetium labelling procedure was above 97, finding agreement with the expected time behaviour of the kit in different organs (Targholizadeh et al. [16]). Therefore, it may be concluded that the long-lived isotopes and 98Tc, certainly present in the radiopharmaceuticals, seem to induce a marginal effect (if any) on the results of the diagnostic procedures. Additional and more accurate investigations are however needed in order to have clearer indications.

It should at last be reminded that the chemical separation process described by Chattopadhyay [32, 33] may be completed very quickly, just within 1 h from EOB, as confirmed by Guérin et al. [15]. In case of a direct accelerator production, it must be assumed that a fully dedicated radiochemical plant for separation and purification processes, up to the final radionuclides in GMP chemical form ready for use, should be available in the same regional area of the production site. In such a way, the already established distribution network for short-lived radionuclides may be exploited.

5. Conclusions

In this work, a comprehensive feasibility study, aimed at the alternative production through high-performance cyclotron-based proton accelerators up to 70 MeV, has been presented. For such a purpose, a wide set of TENDL 2012 theoretical excitation functions have been used, available for the reaction routes on all Mo isotopes present inside -enriched molybdenum metallic material, up to (p,6n), (p,p5n), and (p,2p4n) levels. No simplified constraints, about either the magnitude of cross section values for the various isotope productions (both ground or isomeric states) or assumptions about jumps on the different decay chains (based on the short mean life of radionuclide concerned) have been adopted, thus following the most general theoretical approach.

The estimated theoretical yields obtained in this study are, on the whole, consistent with the previous ones performed by Celler et al. [17], based on a different theoretical approach. An agreement with the experimental results reported in the works by Scholten [4], Takacs et al. [5], and Gagnon et al. [12], mainly about the and and yields by direct reactions on Mo, has been found. For all of the theoretical excitation functions taken into account about the Mo isotopes present inside the target material, a benchmark study has been performed as well in the EXFOR database [34], collecting all of the experimental nuclear reaction data available. A general agreement (although with some deviations) has been found, as may be seen in Figures 2, 3, and 5 for the case of Mo isotope, thus confirming the general reliability of the TENDL 2012 theoretical excitation functions.

This study proves that the expected in-target specific activities for the production (mainly in the 40–70 MeV energy range), even using the incoming high-performance cyclotron at LNL, are about a factor of 104 lower than the one currently available inside the standard Mo/Tc generators. Although the in-target yields estimated for the irradiation parameters investigated may be considered enough, such a route is however not interested as a production way, if maintaining the current standard Mo/Tc generators manufacturing because of the high cost of the target material.

On the contrary, the direct proton-driven production is feasible and quite interesting. Although the in-target yields are improved moving towards beam energies higher than 20 MeV, the production of other Tc contaminants starts increasing to levels that are not negligible. Some impacts may therefore occur on the quality of accelerator-Tc radiopharmaceuticals, as well as on dosimetric aspects. The calculation results clearly show that the energy region of 15–20 MeV is the best operative solution for the accelerator production of . In such a, case the main quality-related parameters, that is, the isotopic purity (IP) as well as the resulting radionuclidic purity (RNP) evolution some hours after EOB, are indeed closer to the generator-produced if short irradiations (i.e., not longer than 3 h) are chosen. Moreover, results also show that the expected specific activity at EOB is similar to the one available from the current Mo/Tc generators, that is, around  Ci/g for 15–20 MeV proton beam. A good balance between production yields and isotopic/radionuclidic purities concern is therefore likely to be found within such constraints. A series of short irradiations are thus preferable instead of a unique per day, in order to minimize the amount of Tc-contaminants. Moreover, a key parameter is the production of the ground state that has a strong influence on the IP parameter. Although recent studies have shown that labelling procedures performed on some commercial kits were not affected using eluates from the standard Mo/Tc generator at 72 h and longer from the previous one, such a consideration suggests that a time window, ranging from 2 h and up to 10 h after EOB, may be exploited for medical procedures with the accelerator- pharmaceuticals. Further studies and experimental campaigns are however requested to validate the (p,2n) excitation function in the entire energy range and to find out the acceptable IP and RNP limits preferably through imaging in vivo tests. These basic aspects will be the subject of future experimental investigations as well as the evaluation of the amount of other technetium isotopes produced (e.g., ) from the irradiation of available enriched targets.

Acknowledgments

The authors would like to thank Professor Adriano Duatti and Professor Giovanni Fiorentini for many helpful suggestions and comments for the paper. The work has been promoted and supported by the INFN 5th National Scientific Commission in the framework of the funded research project APOTEMA (2012–2014).

References

  1. J. E. Beaver and H. B. Hupf, “Production of 99mTc on a medical cyclotron: a feasibility study,” Journal of Nuclear Medicine, vol. 12, no. 11, pp. 739–741, 1971. View at Google Scholar · View at Scopus
  2. N. Levkovskij, “Middle mass nuclides (A=40±100) activation cross sections by medium energy (E=10±50 MeV) protons and α-particles (experiment and systematics),” (Moscow: Inter-Vesti) p. 155, 1991.
  3. P. Chodash, C. T. Angell, J. Benitez et al., “Measurement of excitation functions for the natMo(d,x) 99Mo and natMo(p,x) 99Mo reactions,” Applied Radiation and Isotopes, vol. 69, no. 10, pp. 1447–1452, 2011. View at Publisher · View at Google Scholar · View at Scopus
  4. B. Scholten, R. M. Lambrecht, M. Cogneau, H. Vera Ruiz, and S. M. Qaim, “Excitation functions for the cyclotron production of 99mTc and 99Mo,” Applied Radiation and Isotopes, vol. 51, no. 1, pp. 69–80, 1999. View at Publisher · View at Google Scholar · View at Scopus
  5. S. Takacs, Z. Szucs, F. Tarkanyi, A. Hermanne, and M. Sonck, “Evaluation of proton-induced reactions on 100Mo: new cross sections for the production of 99mTc and 99Mo,” Journal of Radioanalytical and Nuclear Chemistry, vol. 257, pp. 195–201, 2003. View at Google Scholar
  6. M. S. Uddin, M. Hagiwara, F. Tarkanyi, F. Ditroi, and M. Baba, “Experimental studies on the proton-induced activation reactions of molybdenum in the energy range 22–67 MeV,” Applied Radiation and Isotopes, vol. 60, no. 6, pp. 911–920, 2004. View at Publisher · View at Google Scholar · View at Scopus
  7. M. U. Khandaker, M. S. Uddin, K. S. Kim, Y. S. Lee, and G. N. Kim, “Measurement of cross-sections for the (p,xn) reactions in natural molybdenum,” Nuclear Instruments and Methods in Physics Research B, vol. 262, pp. 171–181, 2007. View at Google Scholar
  8. M. B. Challan, M. H. N. Comsan, and M. A. Abou-Zeid, “Thin Target yields and EMPIRE-II predictions on the accelerator production of Technetium-99m,” Journal of Nuclear and Radiation Physics, vol. 1-2, pp. 1–12, 2007. View at Google Scholar
  9. O. Lebeda and M. Pruszyński, “New measurement of excitation functions for (p,x) reactions on natMo with special regard to the formation of 95mTc, 96m+gTc, 99mTc and 99Mo,” Applied Radiation and Isotopes, vol. 68, no. 12, pp. 2355–2365, 2010. View at Publisher · View at Google Scholar · View at Scopus
  10. A. A. Alharbi, A. Azzam, M. McCleskey et al., “Medical radioisotopes production: a comprehensive cross-section study for the production of Mo and Tc radioisotopes via proton induced nuclear reactions on natMo, dioisotopes—applications,” in Bio-Medical Science, N. Singh, Ed., pp. 748–742, InTech, 2011. View at Google Scholar
  11. M. C. Lagunas-Solar, N. X. Zeng, I. Mirshad, and T. Grey-Morgan, “An Update on the Direct Production of 99mTc with Proton Beams and enriched 100Mo targets,” Transactions of the American Nuclear Society, vol. 74, p. 137, 1996. View at Google Scholar
  12. K. Gagnon, F. Bénard, M. Kovacs et al., “Cyclotron production of 99mTc: experimental measurement of the 100Mo(p,x) 99Mo, 99mTc and 99gTc excitation functions from 8 to 18 MeV,” Nuclear Medicine and Biology, vol. 38, no. 6, pp. 907–916, 2011. View at Publisher · View at Google Scholar · View at Scopus
  13. F. Tárkányi, F. Ditrói, A. Hermanne, S. Takács, and A. V. Ignatyuk, “Investigation of activation cross-sections of proton induced nuclear reactions on natMo up to 40 MeV: new data and evaluation,” Nuclear Instruments and Methods in Physics Research B, vol. 280, pp. 45–73, 2012. View at Publisher · View at Google Scholar · View at Scopus
  14. N. Urbano, S. Modoni, M. Guerra, and M. Chinol, “Evaluation of fresh and old eluate of 99Mo/99mTc generators used for labeling of different pharmaceutical kits,” Journal of Radioanalytical and Nuclear Chemistry, vol. 265, no. 1, pp. 7–10, 2005. View at Publisher · View at Google Scholar · View at Scopus
  15. B. Guérin, S. Tremblay, S. Rodrigue et al., “Cyclotron production of 99mTc: an approach to the medical isotope crisis,” Journal of Nuclear Medicine, vol. 51, no. 4, pp. N13–N16, 2010. View at Google Scholar · View at Scopus
  16. H. Targholizadeh, G. Raisali, A. R. Jalilian, N. Rostampour, M. Ensaf, and M. K. Dehghan, “Cyclotron production of technetium radionuclides using a natural metallic molybdenum thick target and consequent preparation of [Tc]-BRIDA as a radio-labelled kit sample,” Nukleonika, vol. 55, no. 1, pp. 113–118, 2010. View at Google Scholar · View at Scopus
  17. A. Celler, X. Hou, F. Bénard, and T. Ruth, “Theoretical modeling of yields for proton-induced reactions on natural and enriched molybdenum targets,” Physics in Medicine and Biology, vol. 56, no. 17, pp. 5469–5484, 2011. View at Publisher · View at Google Scholar · View at Scopus
  18. Isoflex Batch, “Certificate of Analyses Mo-100 enriched molybdenum metallic material,” (ISOFLEX-USA), 2012.
  19. P. Froment, I. Tilquin, M. Cogneau, T. Delbar, J. Vervier, and G. Ryckewaert, “The production of radioisotopes for medical applications by the adiabatic resonance crossing (ARC) technique,” Nuclear Instruments and Methods in Physics Research A, vol. 493, no. 3, pp. 165–175, 2002. View at Publisher · View at Google Scholar · View at Scopus
  20. K. Abbas, S. Buono, N. Burgio et al., “Development of an accelerator driven neutron activator for medical radioisotope production,” Nuclear Instruments and Methods in Physics Research A, vol. 601, no. 3, pp. 223–228, 2009. View at Publisher · View at Google Scholar · View at Scopus
  21. Y. Nagai and Y. Hatsukawa, “Production of 99Mo for nuclear medicine by 100Mo(n, 2n) 99Mo,” Journal of the Physical Society of Japan, vol. 78, no. 3, Article ID 033201, 2009. View at Publisher · View at Google Scholar · View at Scopus
  22. F. Minato and Y. Nagai, “Estimation of production yield of 99Mo for medical use using neutrons from natC(d,n) at Ed = 40 MeV,” Journal of the Physical Society of Japan, vol. 79, no. 9, Article ID 093201, 2010. View at Publisher · View at Google Scholar · View at Scopus
  23. Nuclear Energy Agency, Organisazion for economic co-operation and development 2010 The supply of medical isotopes—review of potential molybdenum-99/technetium-99m production technologies, http://www.oecd-nea.org/med-radio/reports/Med-Radio-99Mo-Prod-Tech.pdf.
  24. A. J. Koning, D. Rochman, S. van der Marck et al., TENDL-2012: TALYS-based evaluated nuclear data library, www.talys.eu/tendl-2012.html.
  25. E. A. Skakun, V. G. Batij, Y. N. Rakivnenko, and O. A. Rastrepin, “Excitation functions and isomer ratios for up-to-9 MeV proton interactions with Zr and Mo isotope nuclei,” Soviet Journal of Nuclear Physics, vol. 46, no. 17, 1987. View at Google Scholar
  26. K. Gagnon, J. S. Wilson, C. M. B. Holt et al., “Cyclotron production of 99mTc: recycling of 100Mo metal targets,” Applied Radiation and Isotopes, vol. 70, pp. 1685–1690, 2012. View at Google Scholar
  27. SRIM, “Particle interactions with matter,” 2012, http://www.srim.org.
  28. NIST Physical Laboratory, STAR database, http://www.nist.gov/pml/data/star/index.cfm.
  29. T. J. Morley, M. Dodd, K. Gagnon et al., “An automated module for the separation and purification of cyclotron-produced 99mTcO4-,” Nuclear Medicine and Biology, vol. 39, no. 4, pp. 551–559, 2012. View at Publisher · View at Google Scholar · View at Scopus
  30. European Commission, “Health and consumers directorate general preliminary report on supply of radioisotopes for medical use and current developments in nuclear medicine,” Tech. Rep., SANCO/C/3/HW, 2009.
  31. Z. Alfassi, F. Groppi, M. Bonardi, and J. Goeij, “On the artificial nature of Tc and the carrier-free nature of 99mTc from 99Mo/99mTc generators,” Applied Radiation and Isotopes, vol. 63, pp. 37–40, 2002. View at Google Scholar
  32. S. Chattopadhyay, S. S. Das, M. K. Das, and N. C. Goomer, “Recovery of 99mTc from Na2[99Mo]MoO4 solution obtained from reactor-produced (n,γ)99Mo using a tiny Dowex-1 column in tandem with a small alumina column,” Applied Radiation and Isotopes, vol. 66, no. 12, pp. 1814–1817, 2008. View at Publisher · View at Google Scholar · View at Scopus
  33. S. Chattopadhyay, S. S. Das, and L. Barua, “A simple and rapid technique for recovery of 99mTc from low specific activity (n,γ)99Mo based on solvent extraction and column chromatography,” Applied Radiation and Isotopes, vol. 68, no. 1, pp. 1–4, 2010. View at Publisher · View at Google Scholar · View at Scopus
  34. Experimental Nuclear Reaction Data (EXFOR), http://www-nds.iaea.org/exfor/exfor.htm.