Laser isotope separation of 176Lu through offtheshelf lasers  Scientific Reports

Abstract . We propose a novel and simple method for the laser isotope separation of 176 Lu a precursor for the production of 177 Lu medical isotope. The physics of the laseratom interaction has been studied through the dynamics of the atomic level populations using the density matrix formalism. It has been shown that a combination of cw excitation lasers and pulsed ionization laser can be used for the laser isotope separation of 176 Lu. The optimum conditions for the efficient and selective separation of 176 Lu have been derived by studying the time evolution of level population under laser excitation. It has also been shown that, it might be possible to produce?~?100% enriched 176 Lu isotope at a rate of 5?mg/h, which is higher than all previously reported methods so far. The isotope separation process proposed can be easily adopted using offtheshelf lasers, for similar atomic systems. Download PDF Introduction . Atomic Vapor Laser Isotope Separation (AVLIS) of actinides (particularly Uranium) has applications in nuclear industry and defence 1 ; therefore, a lot of research has been carried out to achieve the above objective. The technology thus developed is largely limited to a handful of technologically advanced nations and strongly controlled due to proliferation risks. Most of the AVLIS activity is centred on utilization of high repetition rate Copper Vapor Laser (CVL) or SolidState Laser (SSL) pumped dye laser systems for the excitation and ionization of the desired isotope(s). The high pulse repetition frequency (PRF) (10?kHz or higher) of these lasers ensures high interaction efficiency with the atomic vapor enabling the production of the desired isotopes in large scales. It has been realised that AVLIS technology has the ability to play a pivotal role in nuclear medicine 2 . 177 Lu isotope has gained importance in the targeted radionuclide therapy of small tumors and metastatic lesions 3 , 4 . 177 Lu can be produced in a nuclear reactor by irradiation of its precursor (or parent) 176 Lu isotope. The crosssection? 4 of this nuclear reaction is 2090 b. 177 Lu having a halflife of 6.65 d decays to 177 Hf emitting βparticles with energies 497?keV (76%), 384?keV (9.7%) and 176?keV (12%). The 177 Hf which is formed in nuclear excited states decays to the ground state emitting low energy γradiation with energies 208?keV (11%) and 113?keV (6.6%). With mean penetration depths of βparticles in submm, the radionuclide 177 Lu is capable of delivering energy in small volumes, thus having localized cytotoxic effect on the disease cells. The γradiation is useful for imaging and studies of biodistribution and excretion kinetics of the infused radionuclide. Due to these advantages, 177 Lu has evolved as the most effective theranostic (therapeutic?+?diagnostic) isotope in nuclear medicine. However, natural Lu cannot be used for the production of 177 Lu due to the low natural abundance of its precursor isotope 176 Lu (2.59%). It is required to be enriched to?≥?50% level 5 . The typical administered quantity 4 of 177 Lu per course for a patient is?~?7.4?GBq (1.8??g). Therefore, after taking into account the conversion efficiency, 1 mg of enriched 176 Lu can be used for the treatment of?>?500 patients after irradiation. Despite the requirements in small quantities, currently only a very few laboratories in the world have the ability to produce enriched 176 Lu isotope leading to the worldwide shortage of this isotope. This is primarily due to the nonavailability of the suitable commercial tunable narrow band high repetition rate lasers. AVLIS exploits the small resonance frequency shifts between the constituent isotopes of an element (called isotope shifts). The lasers in the stepwise ionization are resonantly tuned to the desired isotope causing preferential ionization enabling the separation from the rest of the isotopes. However, laser isotope separation of Lutetium is extremely complex. The nuclear volume changes 6 in adjacent stable Lu isotopes are small (0.041 fm 2 ) resulting in small field shift. Since the mass shift is also small for the heavy element like Lu, the resultant isotope shifts are small. Fortunately, in case of Lutetium, both stable isotopes of Lu ( 175 Lu and 176 Lu) have widely varying nuclear spin, magnetic moment and quadrupole moment values. Consequently, 176 Lu has a widely spread hyperfine spectrum than the interfering 175 Lu isotope, which can be exploited for the enrichment of 176 Lu. Kurchatov Institute 7 , Moscow has reported successful enrichment of 176 Lu to?>?68% with production rate of?~?4?mg/h. They have used the following photoionization scheme for the enrichment of 176 Lu. $$5d6{s}^{2}\,\,^{2}{D}_{3/2} (0.00 {\mathrm{cm}}^{1}) \stackrel{540.4068 \mathrm{nm}}{\longrightarrow}\,5d6s6p\,\, ^{4}{F}_{5/2}^{o}\left(18504.58 {\mathrm{cm}}^{1}\right)\stackrel{535.0626 nm}{\longrightarrow }$$ $$5d6s7s \,\,^{4}{D}_{3/2}\left({37,193.98} {\mathrm{cm}}^{1}\right)\stackrel{618.0061\,\mathrm{nm}}{\longrightarrow}\,{53,375} {\mathrm{cm}}^{1}Autoionization State \stackrel{}{\longrightarrow }{Lu}^{+}$$ We have recently studied 8 the isotope selective photoionization of 176 Lu and optimization of the process conditions for the enrichment of 176 Lu through the above photoionization scheme using density matrix formalism. The obtained results were in good agreement with the reported experimental data 7 . Further, we have recently reported a new possible photoionization pathway for the enrichment of 176 Lu using broad band dye lasers 9 . However, most of these schemes are accessible by high repetition rate (10–30?kHz) CVL or SSL pumped dye laser systems. Since these laser systems have been categorised as potential dual use technology equipment, hence they have export restrictions to several countries. Therefore, isotope separation using the above scheme(s) cannot be easily adopted by all the producers of nuclear medicine. This is severely restricting the availability of enriched 176 Lu isotope resulting in global shortage of the 177 Lu medical isotope. In order to enhance the production of enriched 176 Lu isotope, it is required to look into the photoionization schemes which are easily accessible by offtheshelf lasers. In the present work, we have theoretically studied the possibility of using offtheshelf lasers for the enrichment of 176 Lu isotope. Photoionization of Lu . Lutetium has two stable isotopes. They are 175 Lu (97.41%) and 176 Lu (2.59%). It has a ground state configuration of 5d6s 2 2 D 3/2 (0.00?cm ?1 ) and has an ionization potential of 43,762.60?cm ?1 . From the known energy levels 10 of Lu, it is possible to formulate innumerable multistep photoionization schemes. Apart from the wellknown first excitation transitions 540.4068?nm and 573.8130?nm accessible by dye lasers, several strong first excitation transitions 11 exist between the wavelength range 272–495?nm. Among them, the transitions between the wavelength range of 350–495?nm are accessible by continuous wave (cw) diode laser or second harmonic of the Ti:Sapphire lasers. The hyperfine structure 13 of Lu isotopes of several photoionization schemes originating from these transitions have been studied for the relative separation of the hyperfine components of the constituent isotopes. A list of possible photoionization schemes which showed good separation of the most intense hyperfine component of 176 Lu from the hyperfine structure of 175 Lu have been tabulated in Table 1 . The requisite wavelengths for these schemes are accessible by commercially available diode lasers or Ti:Sapphire lasers (fundamental frequency or through second harmonic generation). Among them, the following excitation transition has been considered for further discussion. Table 1 Table of two step photoionization schemes for the laser isotope separation of Lu isotopes. Energy level data obtained from Ref. 10 . Full size table 465 nm–1179?nm scheme . $$5d6{s}^{2}\,^{2}{D}_\frac{3}{2}\left(0.00 \,\,{\mathrm{cm}}^{1}\right)\stackrel{465.9316\,\,\mathrm{nm}}{\longrightarrow} 5d6s6p \,\,^{2}{D}_\frac{5}{2}^{o}\left(21462.38 {\mathrm{cm}}^{1}\right)\stackrel{1179.87 \mathrm{nm}}{\longrightarrow}$$ $${5d}^{2}6s\,^{2}{P}_\frac{3}{2}\left({29,937.9} {\mathrm{cm}}^{1}\right)\stackrel{}{\longrightarrow}\,\,Autoionization\,\,State\stackrel{}{\longrightarrow}{Lu}^{+}$$ Photoionization of Lu through pulsed blue dye lasers 14 for the first step excitation transition was reported way back in 1981. Later, photoionization schemes originating from UVblue wavelength region 15 , 16 have been studied for the photoionization of Lu using second and third harmonics of the pulsed Ti:Sapphire lasers to access the requisite wavelengths. However, such lasers have linewidths in several GHz, hence, cannot be easily adopted for the isotope separation of Lu. Utilization of diode lasers for isotope separation was first reported by Olivares et al. 17 , wherein twostep photoionization method was used for the isotope separation of Li isotopes. Recently Matsuoka et al. 18 have studied isotope separation of 48 Ca using diode lasers. In both cases ionization has been carried out by a nonresonance ionization process using pulsed lasers. Further, applicability of such a method for large scale separations has not been stated. The isotope shifts for the first and second excitation transitions of Lu are typically in the range of 400?MHz and 100?MHz respectively. From the hyperfine structure constant data 12 , 13 (Table 2 ), the frequency positions of the twostep hyperfine excitation pathways have been calculated and tabulated in Table 3 . As it can be seen from the Table 3 , the most intense hyperfine excitation pathway 17/2–19/2–17/2 at the frequency position (15,150?MHz, 25,160?MHz) of 176 Lu is well separated from the nearest and the most intense hyperfine excitation pathway 5–65 which is at the frequency position (11,100?MHz, ?17,720?MHz)of 175 Lu. Table 2 Table of HFS coupling constants of lutetium energy levels for the 465–1179?nm photoionization scheme. Full size table Table 3 Table of the hyperfine excitation pathways of 175 Lu and 176 Lu for the 465–1179?nm photoionization pathway (the frequency position is with reference to the centre of gravity of 176 Lu isotope). Full size table The schematic of the photoionization scheme is shown in Fig. 1 . The atoms in the ground state 5d6s 2 2 D 3/2 (0.00?cm ?1 ) are excited to the 5d6s6p 2 D 0 5/2 (21,462.38?cm ?1 ) level using the 465.9316?nm laser and atoms from this level are further excited to the 5d 2 6s 2 P 3/2 (29,937.9?cm ?1 ) level using the 1179.87?nm laser, which are ionized using an ionization laser through either autoionization level or nonresonant ionization. The ionization process is considered as an incoherent process. During the excitation process, due to the finite lifetime of the resonant states, the atoms decay either back to their original states or into the trapped states and lost from ionization process. The excitation rates, decay 11 rates of atoms and the ionization rates are also given in Fig. 1 . Figure 1 Schematic of the photionization pathway of 176 Lu (not to scale). The excitation rate (Rabi frequency), decay rate to lower level and decay rate to trapped levels are represented by Ω, Γ and γ respectively 11 . The radiative decay rate for the 1179.87?nm is taken as 10?MHz. Please see main section for the discussion. Full size image The radiative decay rate of the transition 1179.87?nm has not been reported so far. Several transitions originating from the lower levels ranging between 20,000 cm 1 to 25,000?cm ?1 have shown 11 decay rates of about 1?×?10 7 ?Hz, therefore; the decay rate is taken as 1?×?10 7 ?Hz for the 1179.87?nm transition. Several articles have been published on the even parity autoionization states 19 , 20 , 21 , 22 , 23 wherein two step laser excitation method has been used. However, very little work has been on the odd parity autoionization states. D’yachkov 7 et al. have found several oddparity autoionization states between the energy range 53,067?cm ?1 –53,861?cm ?1 among which, autoionization transition to 53,375?cm ?1 was the strongest. Another work 24 recently reported nineteen new odd parity autoionization levels in the energy range 50,650–51,650?cm ?1 . Using Ti:Sapphire laser for the third excitation step, it is possible to access Rydberg or Autoionization state in the range between 39,938?cm ?1 and 49,560?cm ?1 from the second excitation level 29,937.9?cm ?1 . Since this is a large energy range, it is very likely to find a suitable autoionization state, nonetheless, can be confirmed only by experiments. In a laser isotope separation process, the atomic parameters such as isotope shifts, hyperfine structures of the constituent isotopes, radiative lifetimes of levels, branching ratios and ionization crosssection; laser parameters such as power, spectral bandwidth, frequency, pulse width (or interaction time) and delay between the pulses; effusive atom source parameters such as atomic velocity distribution and angular divergence (collimation ratio), all have a complex interplay on the ionization efficiency and selectivity of the photoionization process. The coherent laser–atom interactions can be accurately described by the density matrix theory 25 . The density matrix equations relevant for the pulsed twostep excitation process have been published in our recent article 8 . When, the atoms traverse perpendicular to the spatially overlapped copropagating laser beams, for all CW lasers case, the laseratom interaction time is determined by the atomic velocity and the laser beam diameter. While in case of CW laser excitation followed by pulsed laser ionization, the laseratom interaction time is determined by the temporal width of the pulsed laser. The coupled density matrix elements have been integrated for the entire laseratom interaction time determined as described above. The effect of linewidth of excitation lasers has been incorporated according to the phase diffusion model as discussed in our recent article 8 . In an AVLIS process, when the first and second excitation lasers are tuned to the frequencies ν 1 and ν 2 respectively, isotope selectivity is defined as $${\text{Isotope Selectivity S }}({\upnu }_{1},{\upnu }_{2})=\left(\frac{{\upeta }_{{176}_{\mathrm{Lu}}}}{{\upeta }_{{175}_{\mathrm{Lu}}}}\right)$$ (1) where, η is the ionization efficiency of the isotope. When the abundance of the isotopes is not equal, one needs to normalize the selectivity with the abundance of the constituent isotopes; thus, isotope ratio enhancement factor is defined as $${\text{Isotope Ratio Enhancement Factor IRE }}({\upnu }_{1},{\upnu }_{2})=\left(\frac{{\upeta }_{{176}_{\mathrm{Lu}}}}{{\upeta }_{{175}_{\mathrm{Lu}}}}\right)\times \left(\frac{{\mathrm{f}}_{{176}_{\mathrm{Lu}}}}{{\mathrm{f}}_{{175}_{\mathrm{Lu}}}}\right)$$ (2) where, f is the fractional abundance of the isotope. Degree of enrichment is calculated using the expression $${\text{Degree of enrichment }}(\mathrm{\%}) =\left\{\frac{{\upeta }_{{176}_{\mathrm{Lu}}}{\times \mathrm{f}}_{{176}_{\mathrm{Lu}}}}{\left({\upeta }_{{175}_{\mathrm{Lu}}}{\times \mathrm{f}}_{{175}_{\mathrm{Lu}}}+{\upeta }_{{176}_{\mathrm{Lu}}}{\times \mathrm{f}}_{{176}_{\mathrm{Lu}}}\right)}\right\}\times 100$$ (3) Results and discussion . For a source temperature of 1600 °C, the most probable velocity of 176 Lu is 421.8?m/s. At this temperature the Doppler broadening of the 465.9316?nm transition is?~?1500?MHz. Let us consider that the atom source to have a lengthtodiameter (L/D) ratio of 5 corresponding to the fullangle divergence (θ) of 22.6°. For the orthogonal laseratom interaction, Doppler broadening along the laser propagation axis reduces to?~?565?MHz. Let us also consider that the laser atom interaction region 26 to have dimensions of 10?mm diameter and 270?mm long. First, we look at the possibility of using cw lasers for excitation and ionization of 176 Lu isotope. Excitation and ionization using cw lasers . As mentioned earlier, the primary aim of the present work is to investigate the enrichment of 176 Lu through offtheshelf lasers. A detailed survey has been carried out on the availability of laser systems in the wavelength ranges indicated in Table 1 . Blue diode laser systems within the wavelength range of the first excitation transition are available with a typical cw power of 30 mW; while diode lasers with accessible wavelengths of the second excitation transition have a typical power of 120 mW. Therefore, we set the limit to the powers of first and second excitation lasers as 30 mW and 120 mW corresponding to the power densities of 38.2 mW/cm 2 and 153 mW/cm 2 respectively. The saturation intensities of the transitions 465.9316?nm and 1179.87?nm are 13 mW/cm 2 and 3 mW/cm 2 respectively which can be easily achieved using the available cw lasers. The frequency of the excitation lasers can be controlled to?~?1?MHz using the commercially available wavelength meters. Nonresonance ionization can be carried out using the high power 10?W cw lasers which are available in wavelengths of 532?nm (Nd:YAG laser) or 527?nm laser (Nd:YVO 4 laser). Therefore, we set the limit to the power of ionization laser to 10?W and the corresponding power density is 12.73?W/cm 2 . Calculation of ionization efficiencies of Lu isotopes has been carried out varying the frequencies of the first and second excitation lasers around the resonance of Lu isotopes and the results are plotted as twodimensional contour plots in Fig. 2 . The contours showed resonances in agreement with expected positions (Table 3 ). Diagonal ridges corresponding to the coherent two photoionization have been observed even at longer detunings due to the strong first and second excitation transitions. The resonance frequency position of the most intense hyperfine excitation pathway 17/2–19/2–17/2 (15,150?MHz, 25,160?MHz) of 176 Lu lies far away from the most intense hyperfine excitation pathway 5–6–5 (11,100?MHz, ?17,720?MHz) of 175 Lu. Figure 2 Two dimensional contour plots of ionization efficiencies of (A) 175 Lu and (B) 176 Lu isotopes under Doppler free conditions. The excitation and ionization laser power densities are 0.038?W/cm 2 , 0.153?W/cm 2 and 12.73?W/cm 2 respectively. Full angle divergence of the atomic beam is 22.6° and both the excitation lasers are counter propagating having a frequency jitter of 1?MHz. Resonance frequency positions of the hyperfine excitation pathways have been numbered as per Table 3 . Resonance frequency position of 17/2–19/2–17/2 (15,150?MHz, 25,160?MHz) of 176 Lu isotope is shown as red circle in (A). Full size image Since the most probable velocity of Lu atoms at 1600 °C is 421.8?m/s, for the cw laser beam diameter of 1?cm, the interaction time of atoms with spatially overlapped laser beams is 23.7??s. The cw lasers typically operate in the fundamental mode TEM 00 having a Gaussian intensity distribution across the diameter; nevertheless it is possible to convert them into a flattop intensity distribution using beam shaping optics. For the present work, we have considered lasers to have flattop intensity distribution across its diameter. First, evolution of level populations with interaction time (Fig.? 3 A,B) have been calculated setting the excitation lasers to the frequencies corresponding to the hyperfine excitation pathway of 17/2–19/2–17/2 (15,150?MHz, 25,160?MHz) of 176 Lu and setting the power of the ionization laser to zero. As expected in the case of nonresonant 175 Lu isotope (Fig.? 3 A), a very small fraction of atoms which are excited into the upper levels decay into the trapped states. At the end of the interaction, about 0.56% atoms can be found in the 1993.92?cm ?1 trapped level and rest of the atoms in the ground state 0.00?cm ?1 . On the other hand, for the case of resonant 176 Lu isotope (Fig.? 3 B), the atoms from the ground F?=?17/2 hyperfine level are excited through the 17/2–19/2–17/2 hyperfine excitation pathway. Since the ionization laser power is set to zero, the atoms excited into the upper hyperfine levels decay to the 1993.92?cm ?1 , 4136.13?cm ?1 and 7476.39?cm ?1 trapped levels. Due to the repetitive excitation and decay processes, at the end of the interaction, the entire population of F?=?17/2 (30%) is transferred to the trapped levels mentioned above. It is therefore expected that the introduction of cw ionization laser transfers atoms from the excited levels into the ionization continuum. The dynamics of the photoionization process with the introduction of a cw ionization laser is discussed below. Figure 3 Time evolution of level populations of the ground, excited, ion and trapped states of 175 Lu (A,C) and 176 Lu (B,D) under cw laser excitation. For graphs (A) and (B) excitation and ionization laser power densities are 0.038?W/cm 2 , 0.153?W/cm 2 and 0.0?W/cm 2 respectively. For graphs (C,D) excitation and ionization laser power densities are 0.038?W/cm 2 , 0.153?W/cm 2 and 12.73?W/cm 2 respectively. Full angle divergence of the atomic beam is 22.6° and both the excitation lasers are counter propagating having frequency jitter of 1?MHz. Excitation lasers are set to the frequencies corresponding to the hyperfine excitation pathway 17/2–19/2–17/2 (15,150?MHz, 25,160?MHz) of 176 Lu isotope. Full size image When a cw ionization laser with a power density of 12.73?W/cm 2 along with the cw excitation lasers interacts with the atomic beam; the ionization efficiency of the nonresonant 175 Lu isotope (Fig.? 3 C) reaches to the value of 2.3?×?10 –10 and the ionization efficiency of resonant 176 Lu isotope (Fig.? 3 D) reaches to 2.5?×?10 –5 implying that most of the population excited into the higher levels decays to the trapped states. The low ionization efficiency of the resonant 176 Lu isotope can be attributed to multiple factors such as low power density of the cw ionization laser and low ionization crosssection?(1?×?10 –16 cm 2 ) of the nonresonant ionization process. From the Eqs.?( 1 )–( 3 ), one can calculate the degree of enrichment, which is?~?100%. Considering the Lu number density of 1?×?10 13 atoms/cm 3 and the dimensions 7 , 26 of the laseratom interaction region of 10?mm diameter and 270?mm length, the production rate of enriched 176 Lu is calculated to be 4.9??g/h. In comparison to the value of 3.7?mg/h reported by D’yachkov et al. 7 , the production rate in the present case is three orders lower. Therefore, enrichment of 176 Lu through “ all cw” laser process is very inefficient. Hence, we have adopted a combination of cw lasers for the excitation process and pulsed laser for the ionization process which is discussed in the next section. cw laser excitation and pulsed laser ionization . As discussed in the previous section, the low ionization efficiency of 176 Lu under cw laser excitation and ionization is due to the low nonresonant ionization crosssection and low peak power of the cw ionization laser. This can be overcome by employing the tunable high repetition rate pulsed laser for the ionization. Theoretical modelling ?of cw laser excitation followed by pulsed laser ionization for the isotope selective photoionization of Sr has been studied previously 27 . Currently high repetition rate (up to 10?kHz) Ti:Sapphire lasers are available with the average power of 1?W. The pulse duration of these lasers is 50?ns. These lasers are also having?a wavelength tunability between 700 and 1000?nm. The lasers can produce peak power density of 2500?W/cm 2 . Wavelength tunability of these lasers allows atoms to be excited either into high lying Rydberg states followed by electric field ionization or directly into the autoionization levels. Crosssection 28 for such a process can be in the range of 1?×?10 –14 cm 2 which is twoorders higher than the crosssection for the nonresonant ionization process. As a result of high peak power density and the high resonance ionization crosssection, ionization efficiency is expected to be enhanced significantly. If one carefully observes the time evolution of level populations of the first (21,462.38?cm ?1 ) and the second (29,937.9?cm ?1 ) excited states of 176 Lu under cw laser excitation (Fig.? 3 A,B); one can find that the atomic population from these states starts decaying into the trapped states after?~?20?ns of its interaction with the excitation lasers. Therefore, for efficient ionization, it is important for the pulsed ionization laser to interact with the “nascent” atomic beam within 50?ns of its interaction with the cw lasers. Experimentally, this can be achieved by passing the diode lasers through mechanical choppers and synchronising the pulsed ionization laser with the sync signal of the mechanical chopper (Supplementary material Fig. S1 ). Alternatively, it is also possible to use highfrequency intensity modulation of narrowlinewidth laser light using fastswitching electrooptical modulators 29 . Evolution of level populations have been calculated for the Lu isotopes by setting the pulsed ionization laser peak power density to 2500?W/cm 2 with a pulse width of 50?ns. The results are plotted in Fig.? 4 . It can be seen from the Fig.? 4 B that, in this case the ionization efficiency of 176 Lu isotope reaches a value of 4.7?×?10 –2 , as compared to the 175 Lu ionization efficiency of 5.7?×?10 –9 (Fig.? 4 A) corresponding to the degree of enrichment of?~?100% of the 176 Lu isotope. Considering the atom number density to be 1?×?10 13 atoms/cm 3 and the laseratom interaction region having dimensions of 10?mm diameter and 270?mm length, the production rate of enriched 176 Lu is calculated to be 2.7?mg/h. This is comparable to the value of 3.7?mg/h reported by D’yachokov et al. 7 . However, in our case the degree of enrichment is significantly higher (100%) than that has been previously reported (68.4%)?value. Figure 4 Time evolution of level populations of the ground, excited, ion and trapped states of (A) 175 Lu and (B) 176 Lu under cw laser excitation and pulsed laser ionization. Excitation laser power densities are 0.038?W/cm 2 , 0.153?W/cm 2 respectively and frequency jitter of excitation lasers is 1?MHz. Peak power density of the ionization laser is 2500?W/cm 2 witha pulsed width of 50?ns. Full angle divergence of the atomic beam is 22.6° and both the excitation lasers are counter propagating. Excitation lasers are set to the frequencies corresponding to the hyperfine excitation pathway 17/2–19/2–17/2 (15,150?MHz, 25,160?MHz) of 176 Lu isotope. Full size image In the “nascent” atomic beam, the population of the ground state 0.00?cm ?1 is distributed among the four hyperfine components F?=?11/2 to F?=?17/2 based on the 2F?+?1 rule. Therefore, the F?=?17/2 hyperfine level possesses 0.30 (30%) of the population of which 5?×?10 –2 is ionized. This implies ionization efficiency of 17% of the targeted F?=?17/2 hyperfine level population. The ionization efficiency can be further increased by increasing the power of the ionization laser; however, it is not practically viable as tunable high repetition rate lasers with average power?>?1?W are not commercially available. Another reason for the low ionization efficiency is due to the large Doppler broadening of the atomic ensemble. As mentioned earlier, for selfcollimating atomic beam effusing out of the atom source having L/D ratio of 5 (full angle divergence of 22.6°), the Doppler broadening along the laser propagation axis will be 565?MHz for the 465.9316?nm transition. The diode laser which can be considered as the single frequency system during short interaction times (50?ns), causes a “hole burning” in the velocity groups of the Doppler broadened atomic ensemble. The group of atoms having velocities nearly perpendicular to the laser propagation axis induce a very small proportion of velocity induced Doppler shifts. When the excitation laser is tuned to the resonance of a transition, this group of atoms of the resonant isotope are excited and ionized. The remaining atoms experiencing the velocity induced Doppler shifts are less likely to get excited and ionized. This causes low ionization efficiency. One of the ways to overcome such a problem is by increasing the linewidth of the excitation lasers; however, this is associated with the loss in selectivity. Therefore, we have adopted an alternative approach. Angular divergence of the atomic beam can be reduced by increasing the L/D ratio of the atom source which causes reduction in the angular divergence and hence velocity induced Doppler shifts. Using this approach, most of the atoms in the atomic beam can be brought into the resonance by adjusting the L/D ratio of the atom source. Calculation of ionization efficiency and degree of enrichment has been carried out varying the L/D ratio of the atom source. It has been found that for a L/D ratio of 10 corresponding to the full angular divergence of 11.5° enhances the ionization efficiency to 8.8?×?10 –2 while the degree of enrichment remains at?~?100%. Such source geometry allows production of 176 Lu at a rate of 5?mg/h. Effect of delay of ionization laser . Effect of delay and delay jitter of the ionization laser pulse with reference to the interaction time has been studied (Supplementary material Fig. S2 ). Ionization efficiency and the degree of enrichment have been calculated varying delay of the ionization laser with reference to the start time of interaction of cw excitation lasers with “nascent” atomic beam. It has been observed that a delay of 50?ns results in highest ionization efficiency. Therefore, it is necessary to set the delay of the ionization laser to 50 ?±?10?ns to achieve the ionization efficiency within 90% of its value. This requirement can be easily met with the currently available technological advances in the synchronization electronics. Conclusion . We have proposed a novel and simple method for the enrichment of 176 Lu isotope a precursor to 177 Lu medical isotope. The physics of the laseratom interaction has been studied through the dynamics of the atomic level populations using density matrix formalism. It has been shown that a combination of cw excitation lasers and pulsed ionization laser can be used for the laser isotope separation of 176 Lu. The optimum conditions for the efficient and selective separation of 176 Lu have been derived by studying the time evolution of level populations during laser excitation. It has been shown that, it might be possible to produce?~?100% enriched 176 Lu isotope at a rate of 5?mg/h, which is higher than all previously reported methods so far. The enrichment process proposed can be easily adopted using offtheshelf lasers for similar atomic systems. Method . The isotope shift and hyperfine structure data of natural Lu isotopes have been taken from the literature and the frequency positions of all hyperfine pathways have been derived with reference to centre of gravity of 176 Lu. A density matrix code for threestep photoionization has been developed. Separate codes have been used for all cw lasers case; cw laser excitation followed by pulsed laser ionization case. Ionization efficiencies of natural Lu isotopes have been calculated for various conditions. From the ionization efficiencies, degree of enrichment of 176 Lu has been derived. Data availability . All the relevant data is provided with this paper. Any additional data relevant to the work in this paper are available upon reasonable request to the corresponding author. Code availability . The code developed is not available for other users. References . 1. Cohen, K. & Murphy, G. M. (eds) The Theory of Isotope Separation as Applied to the Large Scale Production of U235 (McgrawHill Book Company Inc., 1951). Google Scholar ? 2. Eerkens, J.W., Kunze, J.F., & Bond, L.J.. Laser isotope enrichment for medical and industrial applications. in Proceedings of ICONE14 International Conference on Nuclear Engineering, July 17–20, 2006, Miami, USA . 3. de Jong, M. et al. Somatostatin receptortargeted radionuclide therapy of tumors: Preclinical and clinical findings. Semin. Nuclear Med. 32 , 133–140 (2002). Article ? Google Scholar ? 4. Dash, A., Raghavan, M., Pillai, A. & Knapp Jr, F. F. Production of 177 Lu for targeted radionuclide therapy: Available options. Nucl. Med. Mol. Imaging 49 , 85–107 (2015). CAS ? Article ? Google Scholar ? 5. Dvorakova, Z. Production and chemical processing of 177 Lu for nuclear medicine at the Munich research reactor FRMII. Doctoral Thesis, Technische Universiat Munchen (2007). 6. Georg, U. et al. Laser spectroscopy investigation of the nuclear moments and radii of lutetium isotopes. Eur. Phys. J. A 3 , 225–235 (1998). ADS ? CAS ? Article ? Google Scholar ? 7. D’yachkov, A. B. et al. Selective photoionisation of lutetium isotopes. Quantum Electron. 42 , 953–956 (2012). Article ? Google Scholar ? 8. Suryanarayana, M. V. Isotope selective threestep photoionization of 176 Lu. J. Opt. Soc. Am. 38B , 353–370 (2021). ADS ? Article ? Google Scholar ? 9. Suryanarayana, M. V. Isotope separation of 176 Lu a precursor to 177 Lu medical isotope using broadband lasers. Sci. Rep. 11 , 6118 (2021). ADS ? CAS ? Article ? Google Scholar ? 10. Kramida, A., Ralchenko, Yu., Reader, J. &NIST ASD Team (2020). NIST Atomic Spectra Database (version 5.8) (online). https://physics.nist.gov/asd . https://doi.org/10.18434/T4W30F (National Institute of Standards and Technology, 2021). Accessed 12 Jul 2021. 11. Atomic spectral line database from CDROM 23 of R. L. Kurucz. https://www.cfa.harvard.edu/amp/ampdata/kurucz23/sekur.html . 12. Brenner, T., Buttgenbach, S., Rupprecht, W. & Traber, F. Nuclear moments of the low abundant natural isotope 176 Lu and hyperfine anomalies in the lutetium isotopes. Nucl. Phys. A 440 , 407–423 (1985). ADS ? CAS ? Article ? Google Scholar ? 13. Vergès, J. & Wyart, J.F. Infrared emission spectrum of lutecium and extended analysis of Lu I. Phys. Scr. 17 , 495 (1978). ADS ? Article ? Google Scholar ? 14. Bekov, G. I. & VidolovaAngelova, E. P. Optimal scheme for multistage photoionization of lutetium atoms by laser radiation. Sov. J. Quantum Electron. 11 (1), 137–139 (1981). ADS ? Article ? Google Scholar ? 15. Gadelshin, V. et al. Laser resonance ionization spectroscopy on lutetium for the MEDICIS project. Hyperfine Interact 238 , 28 (2017). ADS ? Article ? Google Scholar ? 16. Gadelshin, V. M. et al. Measurement of the laser resonance ionization efficiency for lutetium. Radiochim. Acta 107 (7), 653–661 (2019). CAS ? Article ? Google Scholar ? 17. Olivares, I. E., Duarte, A. E., Saravia, E. A. & Duarte, F. J. Lithium isotope separation with tunable diode lasers. Appl. Opt. 41 , 2973–2977 (2002). ADS ? CAS ? Article ? Google Scholar ? 18. Matsuoka, K. et al. The laser Isotope separation (LIS) methods for the enrichment of 48Ca. J. Phys. Conf. Ser. 1468 , 012199 (2020). CAS ? Article ? Google Scholar ? 19. Bekov, G. I. & VidolovaAngelova, E. P. Optimal scheme for multistage photoionization of lutetium atoms by laser radiation. Sov. J. Quantum Electron. 11 , 137 (1981). ADS ? Article ? Google Scholar ? 20. Xu, C. B. et al. The study of autoionizing states of lutetium atoms by resonance ionization spectroscopy. J. Phys. B At. Mol. Opt. Phys. 26 , 2821–2835 (1993). ADS ? Article ? Google Scholar ? 21. Kujirai, O. & Ogawa, Y. Observation of evenparity autoionization states of lutetium atom by optogalvanic spectroscopy. J. Phys. Soc. Jpn. 67 , 1056–1057 (1998). ADS ? CAS ? Article ? Google Scholar ? 22. Ogawa, Y. & Kujirai, O. Study of evenparity autoionization states of lutetium atom by atomic beam laser resonance ionization spectroscopy. J. Phys. Soc. Japan. 68 , 428–433 (1999). ADS ? CAS ? Article ? Google Scholar ? 23. Li, R. et al. Evenparity Rydberg and autoionizing states of lutetium by laser resonanceionization spectroscopy. Phys. Rev. A 95 , 052501 (2017). ADS ? Article ? Google Scholar ? 24. Rath, A. D., Biswal, D. & Kundu, S. Odd parity autoionization levels of Lu I by threestep photoionization. J. Quant. Spectrosc. Radiat. Transfer 270 , 107696 (2021). CAS ? Article ? Google Scholar ? 25. Shore, B. W. The Theory of Coherent Atomic Excitation: Simple Atoms and Fields Vol. 1 (Wiley, 1990). Google Scholar ? 26. Babichev, A. P. et al. Development of the laser isotope separation method (AVLIS) for obtaining weight amounts of highly enriched 150 Nd isotope. Quantum Electron. 35 , 879–890 (2005). ADS ? CAS ? Article ? Google Scholar ? 27. Jones, O. R., Perks, R. M., Telle, H. H. & Morgan, T. G. Isotopespecific resonance ionization mass spectrometry using continuouswave diode and pulsed dye lasers for trace detection and enrichment of elements of environmental and medical importance, exemplified for strontium (Sr)Part I: Theoretical modelling. Rapid Commun. Mass Spectrom. 10 , 1725–1738 (1996). ADS ? CAS ? Article ? Google Scholar ? 28. D’yachkov, A. B. et al. A study of the kinetic parameters of Lu laser photoionization scheme. Opt. Spectrosc. 128 , 289–296 (2020). ADS ? Article ? Google Scholar ? 29. Voss, A. et al. First use of highfrequency intensity modulation of narrowlinewidth laser light and its application in determination of 206,205,204 Fr groundstate properties. Phys. Rev. Lett. 111 , 122501 (2013). ADS ? Article ? Google Scholar ? Download references Acknowledgements . Authors acknowledge the support of Computational Analysis Division, Bhabha Atomic Research Centre, Visakhapatnam for providing the Super Computer Facility for this work. Author information . Affiliations . Bhabha Atomic Research Centre, Visakhapatnam, Andhra Pradesh, 530 011, India M. V. Suryanarayana?&?M. Sankari Authors M. V. Suryanarayana View author publications You can also search for this author in PubMed ? Google Scholar M. Sankari View author publications You can also search for this author in PubMed ? Google Scholar Contributions . M.V.S. conceived the idea and developed the requisite codes. Both the authors contributed towards computations, data interpretation and manuscript preparation. Corresponding author . Correspondence to M. V. Suryanarayana. Ethics declarations . Competing interests . The authors declare no competing interests. Additional information . Publisher's note . Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations. Supplementary Information . Supplementary Figures. . Rights and permissions . Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http://creativecommons.org/licenses/by/4.0/ . Reprints and Permissions About this article . Cite this article . Suryanarayana, M.V., Sankari, M. Laser isotope separation of 176 Lu through offtheshelf lasers. Sci Rep 11, 18292 (2021). https://doi.org/10.1038/s41598021977738 Download citation Received : 11 March 2021 Accepted : 31 August 2021 Published : 14 September 2021 DOI : https://doi.org/10.1038/s41598021977738 Comments . By submitting a comment you agree to abide by our Terms and Community Guidelines . If you find something abusive or that does not comply with our terms or guidelines please flag it as inappropriate. .

From：


监测目标主题

(1)

(1)

系统抽取对象

机构

(1)

(1)

(2)

(6)

(13)

(4)

(1)

(5)

(1)

(1)

(1)

(1)

(1)

(2)

(1)

(9)

(1)

(2)

活动

(1)

法案

(1)

出版物

(1)

(1)

(1)

地理

(1)

(1)

(1)

(1)

(1)

(2)

(1)

(3)

(1)

(3)

(2)

(1)

(1)

(2)

(3)

(1)

(1)

(1)

(1)

(1)

(1)

人物

(1)

(1)

(1)

(1)

(1)

(2)

(2)

(13)

(1)

系统抽取主题

(1)

(1)

(1)

(1)

(1)

(1)

(1)

(1)

(1)

(1)

(1)

(1)

(1)

(1)

(1)

(1)

(1)

(1)

(1)

(1)

(1)

(1)

(1)

(1)

(1)

(1)

(1)

(1)

(1)

(1)

(1)

(1)

(1)

(1)

(1)

(1)

(1)

(1)

(1)

(1)

(1)

(1)

(1)

(1)

(1)

(1)

(1)

(1)

(1)

(1)

(1)

(1)

(1)

(1)

(1)

(1)

(1)

(1)

(1)

(1)

(1)

(1)

(1)

(1)

(1)

(1)

(1)

(1)

(1)

(1)

(1)

(1)

(1)

(1)

(1)

(1)

(1)

(1)

(1)

(1)

(1)

(1)

(1)

(1)

(1)

(1)

(1)

(1)

(1)

(1)

(1)

(1)

(1)

(1)

(1)

(1)

(1)

(1)

(1)

(1)

(1)

(1)

(1)

(1)

(1)

(1)

(1)

(1)

(1)

(1)

(1)

(1)

(1)

(1)

(1)

(1)

(1)

(1)

(1)

(1)

(1)

(1)

(1)

(1)

(1)

(1)

(1)

(1)

(1)

(1)

(1)

(1)

(1)

(1)

(1)

(1)

(1)

(1)

(1)

(1)

(1)

(1)

(1)

(1)

(1)

(1)

(1)

(1)

(1)

(1)

(1)

(1)

(1)

(1)

(1)

(1)

(1)

(1)

(1)

(1)

(1)

(1)

(1)

(1)

(1)

(1)

(1)

(1)

(1)

(1)

(1)

(1)

(1)

(1)

(1)

(1)

(1)

(1)

(1)

(1)

(1)

(1)

(1)

(1)

(1)

(1)

(1)

(1)

(1)

(1)

(1)

(1)

(1)

(1)

(1)

(1)

(1)

