Tính toán lý thuyết trong điều chế YTTRI 90 và hỗn hợp YTTRI 90 – LUTETI 177

TRƯỜNG ĐẠI HỌC ĐÀ LẠT KHOA KỸ THUẬT HẠT NHÂN PHAN HÀ PHƯƠNG TÍNH TỐN LÝ THUYẾT TRONG ĐIỀU CHẾ YTTRI-90 VÀ HỖN HỢP YTTRI-90 – LUTETI-177 KHOÁ LUẬN TỐT NGHIỆP KỸ SƯ KỸ THUẬT HẠT NHÂN LÂM ĐỒNG, NĂM 2018 TRƯỜNG ĐẠI HỌC ĐÀ LẠT KHOA KỸ THUẬT HẠT NHÂN PHAN HÀ PHƯƠNG – 1410714 TÍNH TỐN LÝ THUYẾT TRONG ĐIỀU CHẾ YTTRI-90 VÀ HỖN HỢP YTTRI-90 – LUTETI-177 KHOÁ LUẬN TỐT NGHIỆP KỸ SƯ KỸ THUẬT HẠT NHÂN GIÁO VIÊN HƯỚNG DẪN THẠC SĨ DƯƠNG VĂN ĐƠNG KHỐ 2014-2019 NHẬN XÉT CỦA GIÁO VIÊN HƯỚNG DẪN ………………………………………………………………………… ………………………………………………………………………… ………………………………………………………………………… ………………………………………………………………………… ………………………………………………………………………… ………………………………………………………………………… ………………………………………………………………………… ………………………………………………………………………… ………………………………………………………………………… ………………………………………………………………………… ………………………………………………………………………… ………………………………………………………………………… ………………………………………………………………………… ………………………………………………………………………… ………………………………………………………………………… ………………………………………………………………………… ………………………………………………………………………… ………………………………………………………………………… ………………………………………………………………………… ………………………………………………………………………… ………………………………………………………………………… ………………………………………………………………………… ………………………………………………………………………… ………………………………………………………………………… NHẬN XÉT CỦA GIÁO VIÊN PHẢN BIỆN ………………………………………………………………………… ………………………………………………………………………… ………………………………………………………………………… ………………………………………………………………………… ………………………………………………………………………… ………………………………………………………………………… ………………………………………………………………………… ………………………………………………………………………… ………………………………………………………………………… ………………………………………………………………………… ………………………………………………………………………… ………………………………………………………………………… ………………………………………………………………………… ………………………………………………………………………… ………………………………………………………………………… ………………………………………………………………………… ………………………………………………………………………… ………………………………………………………………………… ………………………………………………………………………… ………………………………………………………………………… ………………………………………………………………………… ………………………………………………………………………… ………………………………………………………………………… ………………………………………………………………………… LỜI CẢM ƠN Lời đầu tiên, em xin dành lời cảm ơn sâu sắc, tình cảm q mến, kính trọng đến Thạc sĩ Dương Văn Đông anh, chị làm việc Trung tâm nghiên cứu điều chế đồng vị phóng xạ, người tận tình dẫn, giúp đỡ suốt trình nghiên cứu thực khóa luận Viện nghiên cứu hạt nhân Em xin gửi lời cám ơn chân thành tới toàn thể quý thầy cô Trường Đại học Đà Lạt, quý thầy cô khoa Kỹ thuật hạt nhân dạy dỗ, truyền đạt kiến thức quý báu cho em suốt thời gian học tập rèn luyện trường Xin chân thành cảm ơn bạn sinh viên lớp HNK38, người đồng hành học tập, nghiên cứu trao đổi kiến thức Con chân thành cảm ơn gia đình ln đồng hành, hỗ trợ, động viên tin tưởng suốt thời gian qua Sau cùng, dù cố gắng chỉnh sửa luận cách hồn thiện với vốn kiến thức cịn hạn hẹp thời gian thực khóa luận hạn chế nên sai sót điều khơng thể tránh khỏi Vì em mong nhận đóng góp ý kiến q báu từ q thầy Đó hành trang quý giá giúp em hoàn thiện kiến thức sau Em xin chân thành cảm ơn Lâm Đồng, ngày 10 tháng 12 năm 2018 Sinh viên i DANH MỤC CHỮ VIẾT TẮT ĐVPX Đồng vị phóng xạ YHHN Y học hạt nhân MCNP Monte Carlo N-Particle ii DANH MỤC HÌNH ẢNH Hình Sơ đồ hoạt động bình sinh phóng xạ 99mTc từ 99Mo Hình Cấu hình electron yttri (Z = 39) .9 Hình Sơ đồ phân rã 90Sr 10 Hình Cấu hình electron luteti (Z = 71) 12 Hình Hai phương pháp sản xuất 177Lu lò phản ứng .13 Hình Nguyên tắc hoạt động phương pháp Monte Carlo 30 Hình Đồ thị hoạt độ 90Y theo thời gian chiếu .37 Hình Đồ thị hoạt độ 177Lu theo thời gian chiếu 40 Hình Số liệu đặc trưng phản ứng 89Y(n, γ)90Y 44 Hình 10 Kết tra cứu thời gian chiếu cho 177Lu 45 iii DANH MỤC BẢNG BIỂU Bảng Một số hạt nhân phóng xạ thơng dụng dùng y học hạt nhân .5 Bảng Các nguồn sinh đồng vị phóng xạ thường dùng y học hạt nhân Bảng Thông số đầu vào cho tính tốn hoạt độ 90Y .25 Bảng Thơng số đầu vào cho tính tốn hoạt độ 177Lu 26 Bảng Thơng số đầu vào cho tính tốn hoạt độ hỗn hợp 90Y 177Lu 26 Bảng Tổng hợp kết tính với hai mức thơng lượng khác 90Y .33 Bảng Tổng hợp kết tính với hai mức thơng lượng khác 177Lu 34 Bảng Kết tính tốn hoạt độ 90Y chương trình Excel dựng sẵn lị có thơng lượng chuẩn 35 Bảng Kết tính tốn hoạt độ 90Y chương trình Excel dựng sẵn lị có thơng lượng 2.3×1013n/cm2.s lị Đà Lạt 36 Bảng 10 Kết tính tốn hoạt độ 177Lu chương trình Excel dựng sẵn lị có thơng lượng chuẩn 38 Bảng 11 Kết tính tốn hoạt độ 177Lu chương trình Excel dựng sẵn lị có thơng lượng 2.3×1013n/cm2.s lò Đà Lạt .39 Bảng 12 Kết tính tốn hoạt độ tạo thành hỗn hợp 90Y 177Lu chương trình Excel phản ứng 89Y(n,) 90Y 176Lu(n,) 177Lu 41 Bảng 13 Bảng tổng hợp kết tính hoạt độ lý thuyết 90Y 46 Bảng 14 Bảng tổng hợp kết tính hoạt độ lý thuyết 177Lu 46 Bảng 15 Bảng tổng hợp kết tính hoạt độ lý thuyết 46 iv MỤC LỤC MỞ ĐẦU .1 CHƯƠNG – TỔNG QUAN LÝ THUYẾT .3 1.1 Sản xuất đồng vị phóng xạ .3 1.2 Nguyên tắc sản xuất đồng vị phóng xạ .4 1.2.1 Điều chế từ lò phản ứng hạt nhân 1.2.2 Điều chế hạt nhân phóng xạ từ máy gia tốc hạt 1.2.3 Điều chế Generator (nguồn sinh đồng vị phóng xạ) .6 1.3 Sản xuất đồng vị 90Y .8 1.3.1 Khái quát yttri 1.3.2 Sản xuất đồng vị 90Y 1.3.3 Trực tiếp kích hoạt 89Y phản ứng (n, γ) 10 1.3.4 Hệ sinh đồng vị phóng xạ 90Sr/90Y 10 1.4 Sản xuất đồng vị 177Lu 12 1.4.1 Khái quát luteti 12 1.4.2 Sản xuất đồng vị 177Lu 13 1.4.3 Phương pháp kích hoạt trực tiếp qua phản ứng (n,γ) 14 1.4.4 Phương pháp sản xuất gián tiếp từ 176Yb .15 1.5 Phương pháp sản xuất hỗn hợp 90Y 177Lu 16 CHƯƠNG – PHƯƠNG PHÁP NGHIÊN CỨU 17 2.1 Tính tốn hoạt độ tạo thành phương trình kích hoạt phản ứng (n, gamma) 19 2.1.1 Tính tốn hoạt độ 90Y tạo thành từ phản ứng 89Y(n,γ) 90Y 22 2.1.2 Tính tốn hoạt độ 177Lu tạo thành phản ứng 176Lu(n,γ) 177Lu 23 2.1.3 Tính toán hoạt độ hỗn hợp 89 90 Y Lu tạo thành phản ứng 177 Y(n,γ)90Y 176Lu(n,γ) 177Lu .23 2.2 Tính tốn hoạt độ tạo thành chương trình Excel phản ứng (n, gamma) 24 v 2.2.1 Tính tốn hoạt độ tạo thành 90Y chương trình Excel phản ứng 89Y(n,γ) 90Y 25 2.2.2 Tính tốn hoạt độ tạo thành 177Lu chương trình Excel phản ứng 176Lu(n,γ) 177Lu .25 2.2.3 Tính tốn hoạt độ tạo thành hỗn hợp 90Y 177Lu chương trình Excel phản ứng 176Lu(n,γ) 177Lu 89Y(n,γ) 90Y .26 2.3 Tính tốn hoạt độ tạo thành phương pháp Monte Carlo phản ứng (n, gamma) .27 2.3.1 Giới thiệu phương pháp Monte Carlo .27 2.3.2 Giới thiệu chương trình MCNP 30 CHƯƠNG – KẾT QUẢ VÀ THẢO LUẬN 32 3.1 Kết tính tốn hoạt độ phóng xạ 32 3.1.1 Tính tốn hoạt độ 90Y 32 3.1.2 Kết tính tốn hoạt độ 177Lu .33 3.1.3 Kết tính tốn hoạt độ hỗn hợp 90Y 177Lu 34 3.2 Kết tính tốn hoạt độ phóng xạ tạo thành chương trình Excel 35 3.2.1 Kết tính tốn hoạt độ phóng xạ 90Y chương trình Excel phản ứng 89Y(n,) 90Y 35 3.2.2 Kết tính tốn hoạt độ phóng xạ 177Lu chương trình Excel phản ứng 176Lu(n,γ) 177Lu 38 3.2.3 Kết tính toán hoạt độ tạo thành hỗn hợp 90 Y 177 Lu chương trình Excel phản ứng 89Y(n,γ)90Y 176Lu(n,γ)177Lu .41 3.3 Kết kiểm tra hoạt độ tạo thành bảng tra cứu số liệu kích hoạt .43 3.4 Tính tốn hoạt độ phóng xạ phần mềm MCNP 46 3.5 Nhận xét 46 KẾT LUẬN .48 TÀI LIỆU THAM KHẢO .49 PHỤ LỤC 50 vi 2.50000E+05 2.00000E+05 1.50000E+05 1.00000E+05 5.00000E+04 0.00000E+00 20 30 45 55 70 80 95 105 120 135 150 170 192 240 500 3000 Hình Đồ thị hoạt độ 90Y theo thời gian chiếu Nhận xét: Từ kết tính tốn ta nhận biết thêm thơng số quan trọng thời gian chiếu bão hoà, nghĩa cân lượng đồng vị tạo thành phân rã Trong trường hợp 90Y cân đạt chiếu xạ với thời gian 500 37 3.2.2 Kết tính tốn hoạt độ phóng xạ 177 Lu chương trình Excel phản ứng 176Lu(n,γ) 177Lu Bảng 10 Kết tính tốn hoạt độ 177Lu chương trình Excel dựng sẵn lị có thông lượng chuẩn Calculated yields of 177Lu from enriched 176Lu (2.59%) target irradiated in IVV9 Neutron flux (cm2･sec) 1.00E+13 Enrichment of 176Lu 0.0259 Cross section of formation of 176Lu 2050b t1/2 (h) 177Lu 160.8 Number of atoms of 176 Lu 176 Lu in 1µg 2.05E-21 cm2 8.86192E+13 T (hour) A(177Lu) (Bq/µg) 10 15 20 25 30 35 90 95 100 105 110 115 120 125 130 135 140 3.87361E+04 7.66463E+04 1.13748E+05 1.50059E+05 1.78549E+05 1.85595E+05 2.20374E+05 5.84168E+05 6.10448E+05 6.36168E+05 6.61340E+05 6.85974E+05 7.10084E+05 7.33679E+05 7.56772E+05 7.79372E+05 8.01490E+05 8.23136E+05 Theo chương trình ta thấy điều kiện chiếu mẫu hot thu c cho 1àg Lu l: 6.36168ì105 Bq = 0.017 mCi (bia tự nhiên 2,59%) 38 Bảng 11 Kết tính tốn hoạt độ 177Lu chương trình Excel dựng sẵn lị có thơng lượng 2.3×1013n/cm2.s lò Đà Lạt Calculated yields of 177Lu from enriched 176Lu (2.59%) target irradiated in IVV9 Neutron flux (cm2･sec) 2.30E+13 Enrichment of 176Lu 0.0259 Cross section of formation of 176Lu 2050b t1/2 (h) 177Lu 160.8 Number of atoms of 176 Lu 176 Lu in 1µg 8.86192E+13 T (hour) A(177Lu) (Bq/µg) 10 15 20 25 30 35 90 95 100 105 110 115 120 125 130 135 140 8.90931E+04 1.76286E+05 2.61621E+05 3.45135E+05 4.26869E+05 5.06861E+05 5.85146E+05 1.34359E+06 1.40403E+06 1.46319E+06 1.52108E+06 1.57774E+06 1.63319E+06 1.68746E+06 1.74058E+06 1.79255E+06 1.84343E+06 1.89321E+06 39 2.05E-21 cm2 2.00000E+06 1.80000E+06 1.60000E+06 1.40000E+06 1.20000E+06 1.00000E+06 8.00000E+05 6.00000E+05 4.00000E+05 2.00000E+05 0.00000E+00 20 30 45 55 70 80 95 105 120 135 150 170 192 240 500 3000 Hình Đồ thị hoạt độ 177Lu theo thời gian chiếu Nhận xét: Tương tự kết tính tốn trường hợp 90Y Từ kết tính tốn ta nhận thấy thơng số thời gian chiếu bão hoà đạt chiếu xạ với thời gian 500 40 177 Lu, cân 3.2.3 Kết tính tốn hoạt độ tạo thành hỗn hợp 90Y 177 Lu chương trình Excel phản ứng 89Y(n,γ)90Y 176Lu(n,γ)177Lu Bảng 12 Kết tính tốn hoạt độ tạo thành hỗn hợp 90 Y 177 Lu chương trình Excel phản ứng 89Y(n,) 90Y 176Lu(n,) 177Lu Calculated yields of Y-90 and Lu-177 from enriched Y-89 (100%) and Lu-176 (2.59%) target Neutron flux(cm2･sec) 1.0E+13 Enrichment of Y-89 Enrichment of Lu-89 0.0259 Cross section of formation of Y-90 1.3b 1.30E-24 cm2 Cross section of formation of Lu-177 2050b 2.05E-21 cm2 Number of atoms of Y-89, in 1µg Y 89 6.76629E+15 Number of atoms of Lu-176 in 1µg Lu 176 64h t1/2 (h) Y-90 t1/2 (h) Lu-177 160.8h Y Lu 8.86192E+13 A(Y-90) Bq/ug A(Lu-177) Bq/ug 4.63661E+03 3.87361E+04 10 9.02881E+03 7.66463E+04 15 1.31895E+04 1.13748E+05 20 1.71309E+04 1.50059E+05 24 2.01338E+04 1.78549E+05 25 2.08645E+04 1.85595E+05 30 2.44013E+04 2.20374E+05 35 2.77517E+04 2.54411E+05 40 3.09254E+04 2.87723E+05 45 3.39319E+04 3.20324E+05 48 3.56592E+04 3.39550E+05 50 3.67799E+04 3.52230E+05 55 3.94778E+04 3.83456E+05 60 4.20335E+04 4.14016E+05 65 4.44544E+04 4.43924E+05 70 4.67478E+04 4.73195E+05 T(hour) 41 72 4.76309E+04 4.84727E+05 75 4.89202E+04 5.01841E+05 80 5.09782E+04 5.29877E+05 85 5.29276E+04 5.57315E+05 90 5.47743E+04 5.84168E+05 95 5.65237E+04 6.10448E+05 96 5.68624E+04 6.15636E+05 100 5.81809E+04 6.36168E+05 105 5.97507E+04 6.61340E+05 110 6.12377E+04 6.85974E+05 115 6.26464E+04 7.10084E+05 120 6.39808E+04 7.33679E+05 125 6.52449E+04 7.56772E+05 130 6.64423E+04 7.79372E+05 135 6.75766E+04 8.01490E+05 140 6.86512E+04 8.23136E+05 144 6.94699E+04 8.40121E+05 150 7.06333E+04 8.65054E+05 160 7.24120E+04 9.05204E+05 168 7.37026E+04 9.36101E+05 170 7.40081E+04 9.43660E+05 180 7.54404E+04 9.80493E+05 190 7.67256E+04 1.01577E+06 192 7.69664E+04 1.02265E+06 200 7.78790E+04 1.04956E+06 Kết tính hoạt độ Lị Đà Lạt có thơng lượng 2.3×1013 hồn tồn tương tự trường hợp tính phương trình kích hoạt Hoạt độ tạo thành tỷ lệ thuận với thông lượng neutron lị phản ứng hoạt độ hỗn hợp cộng hợp hoạt độ hai đồng vị 90Y 177Lu 42 3.3 Kết kiểm tra hoạt độ tạo thành bảng tra cứu số liệu kích hoạt Bảng tra cứu tài liệu số liệu kích hoạt phản ứng hạt nhân xây dựng sở phương trình kích hoạt có thêm thời gian làm nguội t, phương trình nhân thêm hệ số biểu thị trình phân rã A1 (τ, t) = A1 (τ) D1 (t) (1) Trong (1) phương trình sau: A1 (τ) = Nϕσact Gg 100M [1 − e − 0.693 τ T1 ] Phương trình chung để thiết lập bảng tra cứu là: A1 (x, ϕ, τ, t) = xϕ 107 A1 (τ) D1 (t) dps (2) Từ công thức đưa số liệu chuẩn để tính tốn như: g = 1àg G = 100% = 1013 n.cm-2.s-1 (2.3ì1013 n.cm-2.s-1 Lị Đà Lạt) σact = 1.3×10-24 cm2 (90Y); σact = 2050×10-24 (177Lu) M = 88.905g (89Y); 175.9 (176Lu) T1/2 = 64 (90Y); 160.8 (177Lu) τ: thời gian chiếu xạ neutron biến đổi 43 Hình Số liệu đặc trưng phản ứng 89Y(n, γ)90Y Theo kết tra cứu bảng ta thấy với thời gian chiếu 100 thơng lượng 1013 n/cm2.s, hoạt độ to thnh trờn àg Y l: 5.8246ì104 dps/àg Y tính đơn vị mCi chia cho hệ số 3.7×107 kết 1.57×10-3 mCi 44 Hình 10 Kết tra cứu thời gian chiếu cho 177Lu Theo kết tra cứu bảng ta thấy với thời gian chiếu 100 = 4.17 ngày, thông lượng 1013 n/cm2.s, hoạt độ tạo thành µg Lu khong bng l: 6.5011ì105 dps/àg Lu (giỏ tr ti 4.25 ngày) Để tính đơn vị mCi chia cho hệ số 3.7×107 kết 0.018 mCi Như cơng thức tính cơng thức dùng để tính tốn cho bảng tra cứu liệt kê Điều nói lên hoạt độ tạo thành bia chiếu xạ tỷ lệ thuận với thơng lượng neutron lị phản ứng Bên cạnh hoạt độ tạo thành bão hòa sau thời gian chiếu định phụ thuộc vào thời gian bán rã đồng vị tạo thành 45 Bảng 13 Bảng tổng hợp kết tính hoạt độ lý thuyết 90Y Hoạt độ tính cơng Kết thức lý thuyết Hoạt độ 1.57×10-3 mCi Hoạt độ tra cứu Hoạt độ tính Exel 1.57×10-3 mCi 1.57×10-3 mCi Bảng 14 Bảng tổng hợp kết tính hoạt độ lý thuyết 177Lu Hoạt độ tính cơng Kết thức lý thuyết Hoạt độ 1.71×10-2 mCi Hoạt độ tra cứu Hoạt độ tính 1.8×10-2 mCi Exel 1.71×10-2 mCi 3.4 Tính tốn hoạt độ phóng xạ phần mềm MCNP Khi áp dụng chương trình phần mềm MCNP vướng mắc nhiều vấn đề thông số đầu vào chưa có tài liệu chuẩn để tham khảo chương trình thao tác ln bị gián đoạn, nên kết cuối chưa xuất thành cơng, việc tính tốn cần phải có thời gian để xây dựng lại từ đầu 3.5 Nhận xét Các kết thu phương pháp tính toán lý thuyết tổng hợp Bảng sau: Bảng 15 Bảng tổng hợp kết tính hoạt độ lý thuyết Hoạt độ tính cơng thức Hoạt độ tra cứu Hoạt độ tính excel Hoạt độ tính MCNP Y 1.57×10-3 mCi 1.57×10-3 mCi 1.57×10-3 mCi - Lu 1.71×10-2 mCi 1.8×10-2 mCi 1.71×10-2 mCi - 90 177 Về kết tính tốn từ cơng thức kích hoạt chung từ chương trình Exel dựng sẵng hồn tồn khớp với tài liệu tra cứu kích hoạt, điều nói lên cách tính chương trình tính phù hợp cho việc ứng dụng đánh giá sản xuất đồng vị phóng xạ lị phản ứng hạt nhân Đà Lạt nói riêng hệ thống lị sản xuất nói chung áp dụng cho việc tính tốn hiệu suất sản xuất sản phẩm phản 46 ứng (n, gamma), nhiên kết lý thuyết thực nghiệm chênh lệch nhiều, nhà sản xuất cần có chương trình nghiên cứu tiếp tục tìm hiểu sâu chất điều kiện tối ưu nhằm tăng suất cách hợp lý Từ kết tính tốn ta nhận biết thêm thơng số quan trọng thời gian chiếu bão hoà, nghĩa cân lượng đồng vị tạo thành phân rã Trong trường hợp 90Y ta cân đạt chiếu xạ với thời gian 500 Về chương trình MCNP, lúc đầu khảo sát, quy hoạch, đặt vấn đề lựa chọn đưa thêm chương trình tính tốn với tham vọng tạo tốn đột phá mơ chung nhằm tối ưu trực tiếp cho điều kiện thực nghiệm, song chương trình khơng thể thực q trình hiệu chỉnh thông số đầu vào chưa đủ (không có tài liệu tham khảo cho trường hợp riêng biệt), cần phải có thời gian dài để khảo sát thông số liên quan, khuôn khổ thời gian thực khố luận khơng cho phép Mặc dù chưa có kết việc ứng dụng chương trình MCNP cho mơ hình tính tốn hiệu suất phản ứng hạt nhân để điều chế đồng vị phóng xạ song bước đầu góp phần khởi động ý tưởng cho tính tốn mơ liên quan đến q trình hố học hạt nhân 47 KẾT LUẬN Với lò phản ứng hạt nhân có cơng suất thấp lị phản ứng hạt nhân Đà Lạt nay, khả sản xuất đồng vị phóng xạ 90Y 177Lu hồn tồn khả thi thực cách ổn định để đạt sản phẩm có đủ số lượng chất lượng mong muốn Sự chênh lệch nhiều so với tính tốn lý thuyết điều ngẫu nhiên, dễ hiểu, trách nhiệm nhà sản xuất phải điều chỉnh thực nghiệm để đem lại kết cao điều quan tâm 48 TÀI LIỆU THAM KHẢO Manual for reactor produced radioisotopes, IAEA, Vienna (2003) http://www.benhvien103.vn/vietnamese/bai-giang-chuyen-nganh/yhoc-hat-nhan/duoc-chat-phong-xa/837/ Dương Văn Đơng, Giáo trình Sản xuất đồng vị phóng xạ hợp chất đánh dấu, Viện nghiên cứu hạt nhân (2017) Phan Sỹ An, Y học hạt nhân, Trường Đại học Y Hà Nội (2005) https://vi.wikipedia.org/wiki/Yttri Pillai M R A., Chapter Yttrium-90 production and radiopharmaceuticals development https://vi.wikipedia.org/wiki/Luteti Ashutosh Dash, Production of 177Lu for targeted radionuclide therapy: Available options Đặng Nguyên Phương, 2012, Hướng dẩn sử dụng MCNP cho hệ điều hành windows 10 https://vi.wikipedia.org/wiki/Ph%C6%B0%C6%A1ng_ph%C3%A1p_ Monte_Carlo 11 Hoàng Sỹ Minh Phương, Mơ Monte Carlo chương trình MCNP kiểm chứng thực nghiệm phép đo chiều dày vật liệu hệ chuyên dụng MYO-101, Tạp chí phát triển KH&CN, Tập 13, Số T2 2010 49 PHỤ LỤC Kết tính đối chứng hoạt độ chương trình Neutron Activation Calculator thuộc Dự án WISE Uranium WISE (World Information Service on Energy) tiến hành Hình Giao diện thơng số nhập vào kết tính cho 90Y 50 Hình Giao diện thơng số nhập vào kết tính cho 177Lu 51 ...TRƯỜNG ĐẠI HỌC ĐÀ LẠT KHOA KỸ THUẬT HẠT NHÂN PHAN HÀ PHƯƠNG – 1410714 TÍNH TỐN LÝ THUYẾT TRONG ĐIỀU CHẾ YTTRI- 90 VÀ HỖN HỢP YTTRI- 90 – LUTETI- 177 KHOÁ LUẬN TỐT NGHIỆP KỸ SƯ KỸ THUẬT HẠT NHÂN GIÁO VIÊN... 90Y 176Lu(n,) 177Lu 41 Bảng 13 Bảng tổng hợp kết tính hoạt độ lý thuyết 90Y 46 Bảng 14 Bảng tổng hợp kết tính hoạt độ lý thuyết 177Lu 46 Bảng 15 Bảng tổng hợp kết tính hoạt độ lý thuyết. .. số đầu vào cho tính tốn hoạt độ 90Y .25 Bảng Thơng số đầu vào cho tính tốn hoạt độ 177Lu 26 Bảng Thông số đầu vào cho tính tốn hoạt độ hỗn hợp 90Y 177Lu 26 Bảng Tổng hợp kết tính với