TW201711053A - Method of measurement of physical irradiation for high purity zirconium-89 product obtaining an optimum plating thickness value by plotting a minimum attenuation position according to an interval irradiation power absorption range - Google Patents

Abstract

A method of measurement of physical irradiation for a high purity zirconium-89 product includes steps of: plotting a function graph of incident power of yttrium-89(p,n)zirconium-89 and related nuclear species with respect to reaction cross-sectional areas; selecting a reaction probability or barn value as a horizontal line that intersects the Zr-89 curve at two intersection points; drawing a vertical line that intersects the X-axis to obtain two irradiation power values and taking the values into respective function formula of each curve; obtaining a set of reaction cross-sectional areas by integration; repeating the foregoing steps to generate at least another set of reaction cross-sectional areas; comparing the increase and decrease of each set of reaction cross-sectional areas; selecting a reaction cross-sectional area of a maximum Zr-89 and an tolerable Zr-88 average reaction cross-sectional area ; obtaining two corresponding intersection points by reverse tracking to be optimum irradiation power values; calculating an interval irradiation power absorption range thereof; plotting a function graph of a thickness of the Y-89(p,n)zirconium-89 with respect to attenuation of incident power; selecting an attenuation curve corresponding to the optimum power Eb; and plotting a minimum attenuation position according to the interval irradiation power absorption range to obtain an optimum plating thickness value.

Description

High-purity zirconium (Zr)-89 production physical exposure measurement method

The invention relates to a method for measuring physical irradiation of high-purity zirconium (Zr)-89, in particular to a physical irradiation measuring method capable of obtaining the best generating probability of main nuclear species and avoiding other nuclear species reactions as much as possible.

Generally, in the process of producing high-purity zirconium (Zr)-89, the metal ruthenium (Y)-89 metal ruthenium is solid-plated on a solid target, and the oxidized yttrium (Y)-89 is solid target or After packaging the yttrium (Y)-89 film (foil) after the solid target, since the relationship between the irradiation energy and the thickness of the yttrium (Y)-89 metal plating is not considered, it can only be applied by try and error. The solid target is irradiated with different irradiation energy (MeV), and the activity is measured by a radioactivity measuring instrument, and used to calculate the yield after irradiation.

After the irradiated solid target, if the radioactive nucleus (Y)-89 is washed away from the target of the solid target by using a mineral acid (such as hydrochloric acid, HCl), the activity is measured by a radioactivity measuring instrument and then Many impurities (ie, other nuclear species) can be found by directly adsorbing organic and inorganic adsorbents. When the solid target is irradiated with different irradiation energy (MeV), other nuclear reactions other than the main nuclear reaction are generated in parallel, and a plurality of impurities are generated; and the impurity is caused by the half-life and the half-life of the main nucleus, resulting in a pseudo-radioactive dose (Dose). Value, and cause the generator to wash the ytterbium (Y)-89 after the zirconium (Zr)-89 metal ion is used for labeling the drug, the metal ions in the impurity interfere with the pretreatment efficiency and decrease The yield of the drug label.

In view of the above-mentioned shortcomings in the production method of high-purity zirconium (Zr)-89, the inventors have made research on the above-mentioned shortcomings, and finally the present invention has been produced.

The main object of the present invention is to provide a high-purity zirconium (Zr)-89 production physical irradiation measurement method, which is a 89Y (p, n) 89Zr nuclear reaction incident energy and reaction cross-sectional area function diagram, and 89Y (p, n The application of the 89Zr target thickness and the fundamental physical principle of the incident energy attenuation function diagram can calculate the irradiation energy 參數 of the radioactive nucleus 89-89 (Y-89) solid target process, and the process 产生-89 ( Y-89) The quality of the nuclear species is stable and uniform, and the impurity content can be predicted and controlled by scientific methods and meets its physical and chemical properties.

In order to achieve the above objects and effects, the technical means adopted by the present invention include:

"Drawing y(Y)-89(p,n) zirconium (Zr)-89 and related nuclear species incident energy and reaction cross-sectional area function map" step S11, according to the atomic physical properties, draw 钇(Y)- 89(p,n) zirconium (Zr)-89 and related zirconium (Zr)-88, zirconium (Zr)-87 and other nuclear species, the incident energy and the reaction cross-sectional area function graph, and calculate the function equation of each curve;

"Select a target (Barn) value, make a horizontal line and the zirconium (Zr)-89 incident energy and reaction cross-sectional area curve at the intersection point, draw a vertical line from the two intersection points to the X-axis, The second irradiation energy value (E1, E2) step S12 is to select a target (Barn) value of the reaction cross-sectional area, and make a horizontal line at the target (Barn) value, which is intersected with the radioactive nucleus zirconium (Zr)-89. The energy and the cross-sectional area of the reaction are two intersection points, and a vertical line is drawn from the two intersection points to the X-axis to obtain an irradiation energy value (E1, E2) of the two intersection points;

"The above two illuminating energy values (E1, E2) are substituted into the function formula of each incident energy and the reaction cross-sectional area curve, and the integral can obtain a set of curves corresponding to the two illuminating energy values (E1, E2). Reaction cross-sectional area" step S13;

a "repetition S12 ~ S13 step, another target (Barn) value, and generate a set of reaction cross-sectional area corresponding to each curve" step S14;

A determination step is to determine whether the number of groups available for reference cross-sectional area has reached the required number. If it is negative, then repeat the steps of S12~S13 a target value (Barn), and a set of reaction cross-sectional areas corresponding to each curve" step S15, if affirmative, the next step is performed;

"Comparing the growth and decline of the cross-sectional area of each group; selecting the maximum zirconium (Zr)-89) reaction cross-sectional area Aa-Zr89 and the tolerable zirconium (Zr)-88 average reaction cross-sectional area Bb-Zr88, and deducting the group The two intersection points corresponding to the cross-sectional area of the reaction are taken as the optimal irradiation energy value (Ea, Eb), and the interval irradiation energy absorption range ΔEi of the two intersection points can be calculated. Step S16, the cross-sectional area of each group is compared, and the maximum zirconium is selected. A set of data of the cross-sectional area Aa-Zr89 of Zr)-89 and the average cross-sectional area Bb-Zr88 of the zirconium (Zr)-88 which can withstand (as small as possible), and reversely obtain the corresponding two intersection points as the irradiation energy The value (Ea, Eb), then, the interval energy absorption range ΔEi(MeV)=Eb(MeV)-Ea(MeV) of the two intersection points can be calculated;

"Draw a graph of the yttrium (Y)-89(p,n) zirconium (Zr)-89 target thickness and incident energy attenuation function, select the attenuation curve corresponding to the optimal energy Eb, and illuminate the energy absorption range ΔEi according to the interval Draw the lowest attenuation position to obtain the optimum plating thickness (d) value. Step S17, according to the atomic physical properties, draw the ytterbium (Y)-89(p,n) zirconium (Zr)-89 at different incident energies with respect to The attenuation curve of different target thickness is selected, and the attenuation curve corresponding to the optimal energy Eb is selected. According to the interval, the energy absorption range ΔEi can draw a horizontal line at Ea and the attenuation curve of the optimal energy Eb intersects at the intersection point, and The lower intersection will draw a vertical line to the X-axis to obtain the optimum plating thickness (d) value.

In order to obtain a more specific understanding of the above objects, functions and features of the present invention, the following figures are illustrated as follows:

Referring to Figures 1 to 5, it can be seen that the method of the present invention comprises, in order, the step S11, which is characterized by atomic physical properties of yttrium (Y)-89 (p, n) zirconium (Zr)-89 and related zirconium ( Zr)-88, zirconium (Zr)-87 and other nuclear species in different reaction cross-sectional areas corresponding to different incident energy corresponding points, and then draw the corresponding curves A, B, C (as shown in Figure 2), and calculate The equation for each curve is as follows: Curve A: Zr-89ε(4,20)=-0.0115x ² +0.2829x-1.0542 R ² =0.9216 Curve B: Zr-88ε(10,40)=-0.0038x ² + 0.1899x-1.6021 R ² =0.7915 Curve C: Zr-87ε(26,54)=-0.0016x ² +0.1275x-2.2649 R ² =0.8192

In the second figure, the curve B is located between the curve A and the curve C, and the curve B is partially interlaced with the curve A and the curve C, and the curve A and the curve C are separated from each other.

In step S12, a target (Barn) value of a reaction cross-sectional area is selected, and a horizontal line is formed at the target (Barn) value, and the horizontal line is intersected with a radioactive nucleus zirconium (Zr)-89 incident energy and a reaction cross-sectional area curve (ie, a curve). A) The upper two intersection points, a vertical line is drawn from the two intersection points to the X axis, and the irradiation energy values (E1, E2) of the two intersection points are obtained.

In actual measurement, the target (Barn) value is limited to between 0.5 and 1, in order to obtain data that is closer to practicality, so as to shorten the measurement time.

In an S13 step, the two irradiation energy values (E1, E2) are respectively substituted into a function formula of zirconium (Zr)-89 and Zr)-88 incident energy and a reaction cross-sectional area curve, and the integral can obtain a two-reaction cross-sectional area A1-Zr89. With B1-Zr88, the reaction cross-sectional area A1-Zr89 is the area covered by the curve A between the two irradiation energy values (E1, E2), and the reaction cross-sectional area B1-Zr88 is the curve B for the second irradiation. The area covered by the energy values (E1, E2) (as shown in Figure 3), which is separated from curve A, so there is no area covered.

An S14 step; generating a set of reaction cross-sectional areas AX-Zr89 and BX-Zr88 corresponding to another target (Barn) value in the same manner (as shown in Figure 4).

An S15 step, if yes, the next step is performed.

In step S16, the cross-sectional area of each group is compared, and one of the cross-sectional area Ab-Zr89 having the largest zirconium (Zr)-89 reaction cross-sectional area and the tolerable (as small as possible) zirconium (Zr)-88 average reaction cross-sectional area Bb-Zr88 is selected. Group data, and reversely obtain the corresponding two intersection points as the irradiation energy value (Ea, Eb) (as shown in Fig. 4), and then calculate the interval of the two intersection points, the irradiation energy absorption range ΔEi(MeV)= Eb(MeV)-Ea(MeV);

In step S17, the attenuation curve of the ytterbium (Y)-89(p,n) zirconium (Zr)-89 at different incident energies with respect to different target thicknesses is plotted according to the atomic physical properties, and the optimum energy Eb is selected. The attenuation curve, according to the interval, the energy absorption range ΔEi can draw a horizontal line at Ea and the optimal energy Eb attenuation curve intersects at the intersection point P, and draw a vertical line at the lower intersection point P to the X-axis. Get the best plating thickness (d) value.

In a possible embodiment:

1. Select a target cross-sectional area (Barn) value of 0.5 as a horizontal line, which is intersected at the intersection point of the radioactive nucleus zirconium (Zr)-89 incident energy and the reaction cross-sectional area curve (ie curve A). The second intersection point draws a vertical line to the X axis to obtain the irradiation energy value (E1, E2) of the two intersection points. As shown in Fig. 3, the irradiation energy value (E1, E2) = (6.5, 17.5) .

2. Substituting the two irradiation energy values (E1, E2) = (6.5, 17.5) into the function formula of the incident energy and the cross-sectional area of the zirconium (Zr)-89 and (Zr)-88, respectively, and integrating the curve A The reaction cross-sectional area A1-Zr89 between the two irradiation energy values (E1, E2) and the reaction cross-sectional area B1-Zr88 between the two irradiation energy values (E1, E2).

3. Alternatively, the target (Barn) value is 0.6, and the above steps are repeated to generate another set of reaction cross-sectional areas AX-Zr89 and BX-Zr88 in the same manner, as shown in FIG.

4. Determine whether the number of groups available for reference reaction cross-sectional area has reached the required quantity. If it is insufficient, perform the above steps to generate another set of reaction cross-sectional area. If it is satisfied, compare the cross-sectional areas of each group and select a set of data. Has a maximum zirconium (Zr)-89 reaction cross-sectional area Aa-Zr89 and a tolerable (as small as possible) zirconium (Zr)-88 average reaction cross-sectional area Bb-Zr88, assuming that Aa-Zr89=AX-Zr89 and Bb-Zr88 =BX-Zr88 (i.e., the AX-Zr89 is the maximum zirconium-89 reaction cross-sectional area Aa-Zr89, and the BX-Zr88 is the tolerable zirconium-88 average reaction cross-sectional area Bb-Zr88).

Then, the corresponding two intersection points can be obtained by the Aa-Zr89 and Bb-Zr88 as the irradiation energy value (Ea, Eb) = (8, 16) (as shown in Fig. 4), and the second meeting can be calculated. The interval irradiation energy absorption range of the point ΔEi=Eb-Ea=16-8=8 (MeV).

5. In Fig. 5, select the attenuation curve corresponding to the optimal energy Eb of 16 (MeV). According to the interval, the energy absorption range ΔEi=8 (MeV) can draw a horizontal line at Ea=8 (MeV). The optimal energy Eb is a 16 (MeV) attenuation curve intersecting at the intersection point P, and a vertical line is drawn at the lower intersection point P to the X axis (if the illumination energy is between the displayed values, the inside can be utilized) Insertion calculation), and finally obtain the best plating thickness d = 700mm.

6. Through the above measurement, an optimal illumination energy parameter of 16 MeV can be obtained, and the optimal plating thickness is 700 mm. The other actual illumination parameters of the accelerator can be adjusted according to the empirical value. Finally, the actual illumination parameters are as follows: a. :16 MeV b. Accelerating particles: proton (fixed conditions for cyclotron irradiation) c. Beam current: 200 μA (fixed conditions for cyclotron irradiation) d. Irradiation time: 60 hr (fixed conditions for cyclotron irradiation) e. Irradiation angle: 7 degrees (radiation fixing condition of cyclotron)

The best yield and the lowest other nucleus production can be obtained by irradiation with a cyclotron.

The invention can calculate the irradiation energy parameters of the radioactive nuclear species 钇-89 (Y-89) solid target process by using the above-mentioned parameter measurement method, and the quality of the 钇-89 (Y-89) nuclear species produced by the irradiation is stable and uniform, and the impurity content is It can be predicted and controlled scientifically and conforms to its physical and chemical properties.

In summary, the high-purity zirconium (Zr)-89 physical irradiation measurement method of the present invention can achieve the improvement of the main nuclear species production probability while avoiding the effects of other nuclear species reactions, which is a novelty and progress. The invention of the present invention is filed in accordance with the law; however, the above description is only for the preferred embodiment of the invention, and variations, modifications, alterations or equivalent substitutions of the technical means and scope of the invention are It is also within the scope of the patent application of the present invention.

A‧‧‧Zirconium (Zr)-89 incident energy and reaction cross-sectional area function curve
B‧‧‧Zirconium (Zr)-88 incident energy and reaction cross-sectional area function curve
C‧‧‧Zirconium (Zr)-87 incident energy and reaction cross-sectional area function curve
Reaction cross-sectional area of A1-Zr89, AX-Zr89‧‧‧ zirconium (Zr)-89
Reaction cross-sectional area of B1-Zr88, BX-Zr88‧‧‧ zirconium (Zr)-88
D‧‧‧ plating thickness value
E1, E2, Ea, Eb‧‧‧. Irradiation energy value ΔEi‧‧‧ interval irradiation energy absorption range
P‧‧‧
S11‧‧‧Draw a 數(Y)-89(p,n) zirconium (Zr)-89 and related nuclear species incident energy and reaction cross-sectional area
S12‧‧‧ Select a target (Barn) value, make a horizontal line and the zirconium (Zr)-89 incident energy and reaction cross-sectional area curve at the intersection point, draw a vertical line from the two intersection points to the X-axis , get the two irradiation energy values (E1, E2)
S13‧‧‧ The two illumination energy values (E1, E2) are respectively substituted into the function formulas of the incident energy and the reaction cross-sectional area curve, and the integral can obtain a set of curves between the two illumination energy values (E1, E2). Corresponding reaction cross-sectional area
S14‧‧‧ Repeat steps S12~S13, select another target (Barn) value, and generate a set of reaction cross-sectional areas corresponding to each curve
Does S15‧‧‧ generate another set of reaction cross-sectional areas?
S16‧‧‧Compare the growth and cross-sectional area of each group; select the maximum zirconium (Zr)-89 reaction cross-sectional area Aa-Zr89 and the tolerable zirconium (Zr)-88 average reaction cross-sectional area Bb-Zr88, and reverse the The two intersection points corresponding to the cross-sectional area of the reaction group are taken as the optimal irradiation energy enthalpy (Ea, Eb), and the interval irradiation energy absorption range ΔEi of the two intersection points can be calculated.
S17‧‧‧Draw a 數(Y)-89(p,n) zirconium (Zr)-89 target thickness and incident energy attenuation function diagram, select the attenuation curve corresponding to the optimal energy Eb, and illuminate the energy absorption according to the interval The range ΔEi draws the lowest attenuation position for the best plating thickness (d) value

Figure 1 is a flow chart of the main method of the present invention.

Fig. 2 is a graph showing the incident energy and reaction cross-sectional area of ytterbium (Y)-89(p,n) zirconium (Zr)-89 and related nucleus of the present invention.

Figure 3 is a diagram showing the corresponding identification of a set of relative illumination energy values by taking a target (Barn) value in Figure 2.

Fig. 4 is a diagram showing the corresponding identification value of another set of relative irradiation energy values by taking another target (Barn) value in Fig. 2.

Figure 5 is a graph of the yttrium (Y)-89 (p,n) zirconium (Zr)-89 target thickness and incident energy decay function of the present invention.

S11‧‧‧Draw a graph of incident energy and reaction cross-sectional area of yttrium (Y)-89(p,n) zirconium (Zr)-89 and related nuclear species

S12‧‧‧ Select a target (Barn) value, make a horizontal line and the zirconium (Zr)-89 incident energy and reaction cross-sectional area curve at the intersection point, draw a vertical line from the two intersection points to the X-axis , get the two irradiation energy values (E1, E2)

S13‧‧‧Substituting the above two illumination energy values (E1, E2) into a function formula of each incident energy and the reaction cross-sectional area curve, the integral can obtain a set of curves between the two illumination energy values (E1, E2) Corresponding reaction cross-sectional area

S14‧‧‧ Repeat steps S12~S13, select another target (Barn) value, and generate a set of reaction cross-sectional areas corresponding to each curve

Does S15‧‧‧ generate another set of reaction cross-sectional areas?

S16‧‧‧Compare the growth and cross-sectional area of each group; select the maximum zirconium (Zr)-89 reaction cross-sectional area Aa- Zr89 and the tolerable zirconium (Zr)-88 average reaction cross-sectional area Bb-Zr88, and reversely obtain the two intersection points corresponding to the cross-sectional area of the reaction as the optimal irradiation energy value (Ea, Eb), and can calculate two Intersection illumination energy absorption range ΔEi

S17‧‧‧Draw a graph of the yttrium (Y)-89(p,n) zirconium (Zr)-89 target thickness and incident energy attenuation function, select the attenuation curve corresponding to the optimal energy Eb, and illuminate the energy absorption according to the interval The range ΔEi draws the lowest attenuation position for the best plating thickness (d) value

Claims (2)

A method for measuring physical irradiation of high-purity zirconium (Zr)-89, comprising: an S11 step, according to atomic physical properties, drawing yttrium (Y)-89 (p, n) zirconium (Zr)-89 and Correlation diagram of the incident energy and reaction cross-sectional area of the nuclear species such as zirconium (Zr)-88 and zirconium (Zr)-87, and calculate the functional equation of each curve. One step S12 is to select a target with a cross-sectional area (Barn). The value is a horizontal line at the target (Barn) value, and the horizontal line is intersected with the intersection point of the radioactive nucleus zirconium (Zr)-89 incident energy and the reaction cross-sectional area curve, and a vertical line is drawn from the two intersection points to the X. The axis obtains the irradiation energy value (E1, E2) of the two intersection points; and in the step S13, the two irradiation energy values (E1, E2) are respectively substituted into a function formula of each incident energy and the reaction sectional area curve, and the integral is obtained. a set of reaction curves corresponding to the two irradiation energy values (E1, E2); an S14 step, repeating steps S12 to S13, selecting a target (Barn) value, and generating a set of corresponding curves Reaction cross-sectional area; an S15 judgment step, confirming the group available for reference reaction cross-sectional area Whether the quantity has been met, if it is negative, then execute the step of “repeat S12~S13, select another target (Barn) value, and generate a set of reaction cross-sectional areas corresponding to each curve”. If yes, execute The next step; an S16 step, comparing the cross-sectional areas of each group, selecting the cross-sectional area Bb of the zirconium (Zr)-88 with the largest zirconium (Zr)-89 reaction cross-sectional area Aa-Zr89 and the tolerable (as small as possible) a set of data of -Zr88, and reversely obtain the corresponding two intersection points as the irradiation energy value (Ea, Eb), and then calculate the interval of the second intersection point of the irradiation energy absorption range ΔEi(MeV)=Eb(MeV) -Ea(MeV); A step S17, based on the atomic physical properties, draws the attenuation curve of the ytterbium (Y)-89(p,n) zirconium (Zr)-89 at different incident energies with respect to different target thicknesses. The attenuation curve of the best energy Eb should be based on the energy absorption range ΔEi of the interval, and a horizontal line can be drawn at Ea and the attenuation curve of the best energy Eb intersects at the intersection point, and a vertical line is drawn at the lower intersection point. The optimum plating thickness (d) value can be obtained with the X axis. The high-purity zirconium (Zr)-89 according to Item 1 is a physical irradiation measuring method in which a target (Barn) value is limited to between 0.5 and 1.