Nghiên Cứu Hiệu Chỉnh Thông Số Kỹ Thuật Của Đầu Đo Nal(Tl) Bằng Phương Pháp Monte Carlo

Khóa luận trình bày phương pháp Monte Carlo để hiệu chỉnh thông số kỹ thuật đầu dò nait1 trong vật lý, mang lại kết quả chính xác và đáng tin cậy.

• Tác giả: Nguyễn Thị Tiên
• Người hướng dẫn: TS. Hoàng Đức Tâm
• Trường học: Trường Đại Học Sư Phạm Thành Phố Hồ Chí Minh
• Chuyên ngành: Vật Lý Học
• Đề tài: Nghiên Cứu Hiệu Chỉnh Thông Số Kỹ Thuật Của Đầu Đo Nal(Tl) Bằng Phương Pháp Monte Carlo
• Loại tài liệu: Luận Văn Tốt Nghiệp
• Năm xuất bản: 2016
• Địa điểm: Thành Phố Hồ Chí Minh

Mục lục chi tiết

• LỜI CẢM ƠN
• MỤC LỤC
• 1. CHƯƠNG 1: TƯƠNG TÁC GIỮA BỨC XẠ GAMMA VÀ VẬT CHẤT
• 1.1. Một số tính chất của bức xạ gamma
• 1.2. Các cơ chế tương tác của bức xạ gamma với vật chất
• 1.3. Hiệu ứng quang điện
• 1.4. Hiệu ứng Compton
• 2. CHƯƠNG 2: PHƯƠNG PHÁP MONTE CARLO VÀ CHƯƠNG TRÌNH MÔ PHỎNG MCNP
• 3. CHƯƠNG 3: KHẢO SÁT VÀ HIỆU CHỈNH CÁC THÔNG SỐ KỸ THUẬT CỦA ĐẦU ĐO NHẤP NHÁY Nal(Tl)
• 3.1. Xây dựng bộ ân
• 3.2. Kiểm tra khả năng mô phỏng của tệp đầu vào
• 3.3. Khảo sát, đánh giá và hiệu chỉnh các thông số kỹ thuật của đầu dò
• TÀI LIỆU THAM KHẢO

Tóm tắt

I. Tổng Quan Về Nghiên Cứu Hiệu Chỉnh Thông Số Kỹ Thuật Đầu Đo Nal Tl

Nghiên cứu hiệu chỉnh thông số kỹ thuật của đầu đo Nal(Tl) bằng phương pháp Monte Carlo là một lĩnh vực quan trọng trong vật lý hạt nhân. Đầu đo Nal(Tl) được sử dụng rộng rãi trong việc ghi nhận bức xạ gamma. Tuy nhiên, sau một thời gian hoạt động, các thông số kỹ thuật của đầu đo có thể bị thay đổi, ảnh hưởng đến độ chính xác trong đo lường. Việc hiệu chỉnh các thông số này là cần thiết để đảm bảo độ chính xác và độ tin cậy của các kết quả thu được.

1.1. Tầm Quan Trọng Của Đầu Đo Nal Tl Trong Nghiên Cứu

Đầu đo Nal(Tl) đóng vai trò quan trọng trong việc ghi nhận bức xạ gamma. Nó giúp xác định cường độ và đặc trưng của bức xạ, từ đó phục vụ cho nhiều nghiên cứu trong lĩnh vực vật lý hạt nhân.

1.2. Các Thông Số Kỹ Thuật Cần Hiệu Chỉnh

Các thông số kỹ thuật như độ nhạy, hiệu suất ghi nhận và độ chính xác trong đo lường cần được hiệu chỉnh để đảm bảo kết quả nghiên cứu chính xác và đáng tin cậy.

II. Vấn Đề Thách Thức Trong Hiệu Chỉnh Đầu Đo Nal Tl

Việc hiệu chỉnh đầu đo Nal(Tl) gặp nhiều thách thức do sự thay đổi của các thông số kỹ thuật theo thời gian. Các yếu tố như độ dày lớp vật liệu, cấu trúc của đầu đo và môi trường hoạt động có thể ảnh hưởng đến hiệu suất ghi nhận. Do đó, việc xác định các thông số tối ưu là rất quan trọng.

2.1. Ảnh Hưởng Của Môi Trường Đến Hiệu Suất Đầu Đo

Môi trường hoạt động có thể ảnh hưởng đến độ nhạy và hiệu suất ghi nhận của đầu đo. Các yếu tố như nhiệt độ, độ ẩm và bức xạ nền cần được xem xét trong quá trình hiệu chỉnh.

2.2. Khó Khăn Trong Việc Xác Định Thông Số Tối Ưu

Việc xác định các thông số tối ưu cho đầu đo Nal(Tl) là một thách thức lớn. Cần có các phương pháp chính xác để đo lường và hiệu chỉnh các thông số này.

III. Phương Pháp Monte Carlo Trong Hiệu Chỉnh Đầu Đo Nal Tl

Phương pháp Monte Carlo là một công cụ mạnh mẽ trong việc mô phỏng và hiệu chỉnh các thông số kỹ thuật của đầu đo Nal(Tl). Phương pháp này cho phép mô phỏng các tương tác giữa bức xạ gamma và vật chất, từ đó xác định hiệu suất ghi nhận một cách chính xác.

3.1. Nguyên Tắc Hoạt Động Của Phương Pháp Monte Carlo

Phương pháp Monte Carlo dựa trên việc sử dụng các mô phỏng ngẫu nhiên để tính toán các kết quả. Điều này giúp xác định các thông số kỹ thuật một cách chính xác và hiệu quả.

3.2. Ứng Dụng Của MCNP Trong Nghiên Cứu

Chương trình mô phỏng MCNP được sử dụng để mô phỏng các tương tác giữa bức xạ gamma và đầu đo Nal(Tl). Kết quả từ mô phỏng sẽ được so sánh với thực nghiệm để hiệu chỉnh các thông số kỹ thuật.

IV. Kết Quả Nghiên Cứu Về Hiệu Chỉnh Đầu Đo Nal Tl

Kết quả nghiên cứu cho thấy rằng việc hiệu chỉnh các thông số kỹ thuật của đầu đo Nal(Tl) bằng phương pháp Monte Carlo mang lại độ chính xác cao. Các thông số như độ dày lớp vật liệu và cấu trúc đầu đo đã được tối ưu hóa để đạt được hiệu suất ghi nhận tốt nhất.

4.1. So Sánh Kết Quả Mô Phỏng Và Thực Nghiệm

Kết quả từ mô phỏng Monte Carlo được so sánh với dữ liệu thực nghiệm. Sự phù hợp giữa hai kết quả này cho thấy tính chính xác của phương pháp hiệu chỉnh.

4.2. Đánh Giá Hiệu Suất Ghi Nhận Sau Hiệu Chỉnh

Sau khi hiệu chỉnh, hiệu suất ghi nhận của đầu đo Nal(Tl) đã được cải thiện đáng kể. Điều này cho thấy tầm quan trọng của việc hiệu chỉnh các thông số kỹ thuật.

V. Kết Luận Về Nghiên Cứu Hiệu Chỉnh Đầu Đo Nal Tl

Nghiên cứu hiệu chỉnh thông số kỹ thuật của đầu đo Nal(Tl) bằng phương pháp Monte Carlo đã chứng minh được tính hiệu quả và độ chính xác cao. Việc này không chỉ giúp cải thiện hiệu suất ghi nhận mà còn mở ra hướng đi mới cho các nghiên cứu trong lĩnh vực vật lý hạt nhân.

5.1. Tương Lai Của Nghiên Cứu Hiệu Chỉnh

Nghiên cứu trong tương lai có thể mở rộng ra các loại đầu đo khác và áp dụng các phương pháp mô phỏng tiên tiến hơn để nâng cao độ chính xác trong đo lường.

5.2. Ứng Dụng Thực Tiễn Của Kết Quả Nghiên Cứu

Kết quả nghiên cứu có thể được áp dụng trong nhiều lĩnh vực khác nhau, từ y học đến công nghiệp, nơi mà việc ghi nhận bức xạ gamma là cần thiết.

Nội dung chính

Tổng quan nghiên cứu

Trong vật lý hạt nhân thực nghiệm và quan trắc phóng xạ môi trường, hệ phổ kế gamma sử dụng đầu dò nhấp nháy NaI(Tl) kích thước 76x76 mm chiếm hơn 60% các ứng dụng thực tế nhờ hiệu suất ghi cao và tính kinh tế vượt trội. Tuy nhiên, sau thời gian dài vận hành, các yếu tố lão hóa quang học và sự thiếu đồng nhất về thông số chế tạo thực tế so với thiết kế lý thuyết dẫn đến sai số hiệu suất ghi nhận từ 5% đến 18%. Vấn đề đặt ra là các thông số hình học do nhà sản xuất công bố chưa phản ánh đầy đủ độ dày thực của các lớp bao bọc tinh thể, khiến việc mô phỏng tính toán hiệu suất trên máy tính gặp sai lệch lớn so với phép đo thực tế.

Mục tiêu cốt lõi của nghiên cứu là khảo sát, đánh giá định lượng và hiệu chỉnh chi tiết các thông số kỹ thuật của đầu dò nhấp nháy NaI(Tl) Gamma-Rad5 bằng phương pháp Monte Carlo thông qua chương trình MCNP5. Đề tài tập trung xác định hiệu suất đỉnh năng lượng toàn phần trên dải năng lượng từ 88,03 keV đến 1332,49 keV nhằm đưa độ lệch tương đối giữa mô phỏng và thực nghiệm về mức dưới 2%. Phạm vi nghiên cứu được triển khai tại phòng thí nghiệm chuyên ngành với khoảng cách chuẩn hóa nguồn - đầu dò là 20 cm, sử dụng bộ 7 nguồn chuẩn phóng xạ dạng điểm. Kết quả nghiên cứu mang ý nghĩa thực tiễn sâu sắc, giúp các cơ sở hạt nhân tối ưu hóa quy trình hiệu chuẩn, tiết kiệm hơn 70% thời gian thực nghiệm và nâng cao độ tin cậy trong định lượng hoạt độ phóng xạ môi trường.

Cơ sở lý thuyết và phương pháp nghiên cứu

Khung lý thuyết áp dụng

Nghiên cứu được xây dựng trên nền tảng lý thuyết tương tác giữa bức xạ gamma và vật chất trong dải năng lượng từ vài chục keV đến dưới 2 MeV. Quá trình này được chi phối bởi ba cơ chế vật lý cơ bản:

Thứ nhất, hiệu ứng quang điện là quá trình lượng tử gamma truyền toàn bộ năng lượng cho electron liên kết ở các lớp trong nguyên tử. Tiết diện quang điện tỷ lệ thuận với lũy thừa bậc cao của nguyên tử số Z và giảm nhanh theo hàm bậc 3,5 của năng lượng gamma tới, chiếm ưu thế tuyệt đối ở dải năng lượng thấp dưới 200 keV.

Thứ hai, tán xạ Compton xảy ra khi photon va chạm với electron liên kết yếu, truyền một phần năng lượng cho electron giật lùi và tán xạ lệch hướng theo công thức Klein-Nishina. Đây là kênh tương tác chủ đạo trong khoảng năng lượng trung bình từ 0,2 MeV đến 5,0 MeV.

Thứ ba, hiệu ứng tạo cặp electron - positron xảy ra khi photon có năng lượng lớn hơn 1,022 MeV tương tác trong điện trường hạt nhân. Sau khi hạt positron giảm chậm sẽ hủy cặp với electron môi trường, giải phóng 2 photon có năng lượng 0,511 MeV bay theo hai hướng ngược nhau.

Bên cạnh đó, mô hình đánh giá phổ kế áp dụng hai khái niệm then chốt là hiệu suất đỉnh năng lượng toàn phần (FEPE) và độ phân giải năng lượng (FWHM). Độ mở rộng phổ dạng Gauss được chuẩn hóa qua phương trình thực nghiệm liên hệ giữa bề rộng nửa chiều cao đỉnh với căn bậc hai của năng lượng gamma theo hệ số chuẩn hóa 2,355 lần độ lệch chuẩn.

Phương pháp nghiên cứu

Nghiên cứu sử dụng phương pháp thống kê ngẫu nhiên Monte Carlo được hiện thực hóa qua chương trình tính toán MCNP5. Dữ liệu thực nghiệm được thu thập từ bộ 7 nguồn chuẩn phóng xạ gồm Cadimi-109 (88,03 keV), Natri-22 (511 keV và 1274,54 keV), Cesi-137 (661,66 keV), Mangan-54 (834,84 keV), Kẽm-65 (1115,54 keV), Coban-60 (1173,23 keV và 1332,49 keV) cùng nguồn đa năng Europi-152 đặt cách bề mặt đầu dò 20 cm.

Mô hình không gian 3 chiều trong MCNP5 được phân chia thành 14 ô cơ bản (cell) mô phỏng chính xác cấu trúc tinh thể NaI(Tl), lớp bột phản xạ oxit nhôm ($Al_2O_3$), lớp đệm silicon, vỏ bọc nhôm và ống nhân quang điện (PMT). Cỡ mẫu mô phỏng được thiết lập ở mức 1.000.000 lịch sử hạt (NPS = $10^6$) cho mỗi cấu hình để đảm bảo sai số thống kê tương đối R luôn nhỏ hơn 0,05, đạt độ tin cậy chuẩn cho các phép đo phổ kế hạt nhân.

Lý do lựa chọn chương trình MCNP5 kết hợp thẻ phân bố độ cao xung F8 (Pulse Height Tally) chia thành 8192 bin năng lượng là vì thuật toán này cho phép theo dõi toàn bộ vết chuyển động và sự truyền năng lượng của từng photon trong vật liệu phức tạp. Để tái hiện chính xác phổ thực tế, tham số làm rộng đỉnh dạng Gauss (GEB) được làm khớp theo hàm đa thức bậc hai với các hệ số thực nghiệm a = -0,0137257 MeV, b = 0,0684 MeV$^{1/2}$ và c = 0,0012 MeV$^{-1}$, tạo nền tảng vững chắc cho quá trình phân tích và đối sánh dữ liệu.

Kết quả nghiên cứu và thảo luận

Những phát hiện chính

Quá trình khảo sát và mô phỏng trên chương trình MCNP5 đã đem lại 3 phát hiện khoa học mang tính đột phá:

Thứ nhất, việc mô phỏng đầu dò chỉ với khối tinh thể NaI(Tl) trần mà bỏ qua các lớp bao bọc tạo ra sai số rất lớn so với thực nghiệm. Cụ thể, tại mức năng lượng 661,66 keV của Cesi-137, sai lệch hiệu suất lên tới 17,76%, và tại đỉnh 1115,54 keV của Kẽm-65, độ lệch là 10,08%. Điều này chứng minh các lớp vật liệu bảo vệ đóng vai trò quyết định đến việc hấp thụ và tán xạ bức xạ trước khi tới vùng nhạy của tinh thể.

Thứ hai, việc hiệu chỉnh lớp phản xạ oxit nhôm ($Al_2O_3$) mặt trước dày 1,0 mm, mặt bên dày 2,0 mm kết hợp lớp silicon 2,0 mm và vỏ nhôm 1,5 mm đã cải thiện vượt bậc độ chính xác của mô hình. Tỷ số hiệu suất đỉnh toàn phần giữa mô phỏng và thực nghiệm đạt mức lý tưởng, đưa độ lệch tương đối của cả 8 mức năng lượng khảo sát xuống dưới 2% (từ 0,05% đến 1,86%).

Thứ ba, khảo sát độc lập về bề dày của ống nhân quang điện (PMT) mô phỏng bằng khối nhôm đặc từ 0 cm đến 9 cm cho thấy hiệu suất đỉnh toàn phần hầu như không thay đổi. Biến thiên hiệu suất chỉ dao động trong khoảng cực hẹp từ 0,01% đến 0,11%, chứng minh bề dày PMT ít ảnh hưởng trực tiếp đến hiệu suất đỉnh năng lượng toàn phần.

Thảo luận kết quả

Nguyên nhân dẫn đến sai lệch lớn khi không mô phỏng các lớp vật liệu bao quanh là do các photon gamma năng lượng thấp dễ bị tán xạ ngược hoặc suy giảm bởi vỏ bọc kim loại trước khi vào tinh thể. Khi bổ sung đầy đủ lớp oxit nhôm, silicon và vỏ nhôm theo đúng cấu trúc hiệu chỉnh, sự mất mát năng lượng do tán xạ ngoài được tái hiện hoàn toàn trùng khớp với điều kiện thực tế.

Đối với ống nhân quang PMT, do nằm ở mặt đáy phía sau tinh thể, các photon tới từ hướng mặt trước đã bị tinh thể NaI(Tl) kích thước 76x76 mm hấp thụ phần lớn năng lượng thông qua hiệu ứng quang điện và tán xạ Compton nhiều lần. Lượng photon xuyên qua đáy để tán xạ ngược trở lại từ PMT là rất nhỏ, giải thích tại sao thay đổi bề dày khối nhôm PMT từ 0 cm lên 9 cm không gây tác động đáng kể đến hiệu suất ghi nhận.

Dữ liệu nghiên cứu được biểu diễn trực quan qua các bảng đối chiếu hiệu suất FEPE và đồ thị phổ xung 8192 kênh. Trên biểu đồ phổ gamma của Coban-60 và Cesi-137, đường cong phổ mô phỏng MCNP5 sau khi áp dụng hiệu chỉnh GEB hoàn toàn trùng khít với đỉnh thực nghiệm cả về vị trí lẫn độ rộng chân phổ. So sánh với các công bố của Shi và cộng sự (vốn cố định bề dày nhôm PMT ở mức 3,0 cm), nghiên cứu này chứng minh rằng việc tối giản hóa khối PMT trong tệp đầu vào không làm suy giảm độ chính xác mà còn giúp tiết kiệm đáng kể tài nguyên tính toán.

Đề xuất và khuyến nghị

Nhằm nâng cao hiệu quả ứng dụng phương pháp Monte Carlo trong phân tích hạt nhân, tác giả đưa ra 4 khuyến nghị hành động cụ thể:

1. Chuẩn hóa bộ thông số hình học đầu vào MCNP5: Các viện nghiên cứu và phòng thí nghiệm hạt nhân cần áp dụng bộ thông số đã hiệu chỉnh gồm lớp oxit nhôm mặt trước 1,0 mm, mặt bên 2,0 mm, lớp silicon 2,0 mm và vỏ nhôm 1,5 mm cho đầu dò NaI(Tl) 76x76 mm. Giải pháp này giúp kiểm soát sai số mô phỏng dưới mức 1,5% ngay trong giai đoạn thiết lập tệp đầu vào trong vòng 1 đến 3 tháng tới.

2. Tối giản hóa mô hình ống nhân quang PMT: Các kỹ sư tính toán nên thay thế cấu trúc PMT phức tạp bằng một khối nhôm hình trụ đồng trục đơn giản có bề dày từ 2 cm đến 3 cm. Việc này giúp giảm hơn 25% thời gian xử lý của bộ xử lý trung tâm (CPU) trong các bài toán mô phỏng lớn với số hạt trên $10^7$ mà vẫn đảm bảo độ tin cậy.

3. Thiết lập quy trình tái hiệu chuẩn định kỳ 12 tháng: Trung tâm an toàn bức xạ cần ban hành quy định đánh giá lại hệ số làm rộng đỉnh GEB (a, b, c) hàng năm cho các hệ phổ kế NaI(Tl) nhằm cập nhật kịp thời sự trôi dạt thông số do lão hóa linh kiện điện tử và tinh thể.

4. Ứng dụng nguồn đa năng Europi-152 trong kiểm soát chất lượng phổ: Bộ phận quan trắc phóng xạ môi trường nên sử dụng nguồn chuẩn Europi-152 với hơn 10 vạch phổ đặc trưng để kiểm tra đồng thời tính tuyến tính năng lượng và hình dạng cạnh Compton trước mỗi đợt đo mẫu định kỳ.

Đối tượng nên tham khảo luận văn

Nội dung luận văn mang lại giá trị ứng dụng thiết thực cho 4 nhóm đối tượng chuyên môn:

Nhóm thứ nhất là các nhà nghiên cứu và kỹ sư vật lý hạt nhân tại các viện chuyên ngành. Luận văn cung cấp phương pháp chi tiết để xây dựng tệp đầu vào MCNP5 chuẩn xác, giúp tối ưu hóa thiết kế che chắn và hiệu chuẩn đầu dò mà không cần tiến hành quá nhiều phép đo thực nghiệm tốn kém.

Nhóm thứ hai là cán bộ quản lý và quan trắc an toàn bức xạ môi trường. Tài liệu giúp nâng cao độ chính xác khi xác định hoạt độ của các đồng vị phóng xạ nhân tạo như Cesi-137 và Coban-60 trong mẫu đất, nước và không khí xung quanh các cơ sở bức xạ.

Nhóm thứ ba là giảng viên, học viên cao học và sinh viên ngành Vật lý học, Kỹ thuật hạt nhân. Luận văn là tài liệu tham khảo mẫu mực về cấu trúc khai báo Cell Card, Surface Card, Material Card và kỹ thuật làm khớp hàm Gauss GEB trên 8192 kênh phổ kế.

Nhóm thứ tư là các đơn vị phát triển thiết bị đo lường hạt nhân. Các dữ liệu khảo sát độ dày lớp phản xạ $Al_2O_3$ và lớp đệm silicon cung cấp chỉ dẫn kỹ thuật giá trị để tối ưu hóa quy trình đóng gói tinh thể và nâng cao hiệu suất truyền sáng tới photocathode.

Câu hỏi thường gặp

1. Tại sao cần phải hiệu chỉnh các thông số kỹ thuật của đầu dò thay vì dùng trực tiếp số liệu nhà sản xuất? Thông số nhà sản xuất thường là kích thước danh định và có dung sai chế tạo. Ngoài ra, sau thời gian vận hành, các lớp đệm quang học có thể bị biến tính. Nghiên cứu chứng minh rằng nếu chỉ dùng mô hình tinh thể trần hoặc thông số chưa hiệu chỉnh, sai số hiệu suất tại đỉnh 661,66 keV lên tới 17,76%, gây sai lệch nghiêm trọng trong định lượng hoạt độ phóng xạ.

2. Bề dày của ống nhân quang điện (PMT) ảnh hưởng như thế nào đến hiệu suất đỉnh năng lượng toàn phần? Khảo sát mô phỏng với bề dày PMT thay đổi từ 0 cm đến 9 cm cho thấy hiệu suất ghi nhận chỉ biến thiên cực nhỏ từ 0,01% đến 0,11%. Độ lệch tương đối giữa mô phỏng và thực nghiệm luôn duy trì ổn định dưới 2%, chứng minh thông số này không tác động đáng kể đến hiệu suất đỉnh toàn phần.

3. Phương pháp Monte Carlo có ưu thế gì vượt trội so với các phép đo thực nghiệm truyền thống? Phương pháp Monte Carlo cho phép mô phỏng vận chuyển hạt ngẫu nhiên trong không gian 3 chiều phức tạp với độ chính xác cao. Chạy mô phỏng $10^6$ hạt giúp thu được sai số thống kê R dưới 0,05, giúp khảo sát sâu các chi tiết cấu trúc bên trong đầu dò mà mắt thường và thực nghiệm không thể can thiệp trực tiếp.

4. Thông số hiệu chỉnh làm rộng Gauss (GEB) có ý nghĩa gì trong bài toán mô phỏng MCNP5? Trong mô phỏng lý thuyết, các vạch phổ thu được là các đường thẳng đơn sắc. Việc khai báo thẻ FT8 GEB với các tham số a = -0,0137257 MeV, b = 0,0684 MeV$^{1/2}$ giúp chuyển đổi phổ lý thuyết thành các đỉnh phân bố Gauss mở rộng, tương thích hoàn toàn với 8192 kênh của hệ phổ kế thực nghiệm.

5. Bề dày lớp oxit nhôm ($Al_2O_3$) và silicon được khuyến nghị là bao nhiêu để đạt kết quả tối ưu? Kết quả nghiên cứu xác nhận bề dày tối ưu của lớp phản xạ $Al_2O_3$ mặt trước là 1,0 mm, mặt bên là 2,0 mm kết hợp với lớp đệm silicon 2,0 mm. Bộ thông số này đảm bảo độ lệch tương đối giữa tính toán MCNP5 và số liệu thực nghiệm luôn dưới 2% trên toàn bộ dải năng lượng từ 88,03 keV đến 1332,49 keV.

Kết luận

• Thiết lập thành công mô hình hình học 3 chiều chi tiết gồm 14 ô cơ bản cho hệ phổ kế NaI(Tl) kích thước 76x76 mm trên chương trình MCNP5.
• Khẳng định vai trò quyết định của lớp phản xạ $Al_2O_3$ (1,0 mm), silicon (2,0 mm) và vỏ nhôm (1,5 mm) trong việc giảm sai số tính toán từ 17,76% xuống dưới 2%.
• Chứng minh tính bất biến thực tiễn của bề dày ống nhân quang PMT (biến thiên dưới 0,11%), cho phép tối giản hóa mô hình tính toán để tiết kiệm thời gian chạy mô phỏng.
• Xác lập bộ thông số GEB chuẩn hóa giúp tái hiện hoàn hảo hình dạng phổ thực nghiệm trên hệ thống phân tích 8192 kênh.
• Hoàn thiện quy trình hiệu chỉnh thông số kỹ thuật cho 8 mức năng lượng gamma chuẩn từ 88,03 keV đến 1332,49 keV đạt độ tin cậy khoa học cao.

Công trình là đóng góp quan trọng cho việc chuẩn hóa phương pháp mô phỏng phổ kế hạt nhân tại Việt Nam. Trong lộ trình 6 đến 12 tháng tới, hướng nghiên cứu cần mở rộng sang các dòng đầu dò nhấp nháy thế hệ mới và dải năng lượng cao trên 3 MeV. Hãy tải toàn văn luận văn và áp dụng ngay bộ thông số hiệu chỉnh MCNP5 vào hệ thống nghiên cứu của bạn để nâng cao tối đa độ chính xác của các phép đo hạt nhân.

Trích đoạn nội dung tài liệu

CHƯƠNG 1. Tương tác giữa bức xa gamma và vật chat 1. Một số tính chất của bức xạ gamma Bức xạ gamma lả các photon có năng lượng cao, với bản chất là sóng điện từ với bước sóng nhỏ hơn nhiều so với kích thước nguyên tử. Năng lượng của nó được tính theo công thức: (1.1) Bức xạ gamma có khả năng xuyên sâu rat lớn và khi đi qua vật chất chúng bị hap thụ do tương tác điện từ.

Tuy nhiên, cơ chế của quá trình này khác với các hạt tích điện. Các hạt tích điện va chạm nhiều với các electron và hạt nhân nên phương bay bị lệch nhiều so với ban đầu, trong khi bức xạ gamma khi qua môi trường vật chất xác suất xảy ra va chạm với các electron và hạt nhân là thấp hơn so với các hạt tích điện. Nguyên nhân có sự khác nhau này là do lượng tử gamma không mang điện nên không chịu tác dụng của lực Coulomb tương tác xa. tương tác với electron trong miền với bán kính cỡ 103m.

Các hạt tích điện va chạm nhiều với các electron nguyên tử nên dé dang bị làm chậm trong môi trường, có quãng chạy hữu hạn trong vật chất nên có thê bị hấp thụ hoàn toàn. Khi tăng bè dày vật chat thì chùm tia gamma chỉ suy giảm về cường độ mà không bị hap thụ hoàn toàn nên không có khái niệm quãng chạy. Các cơ chế tương tác của bức xạ gamma với vật chất Bức xạ gamma tương tác với nguyên tử vật chất thông qua nhiều cơ chế khác nhau, các tương tác này không gây ra hiện tượng ion hóa trực tiếp như hạt tích điện nhưng có kha năng làm but electron quỹ đạo ra khỏi nguyên tử. Cac electron tự do này bị làm chậm trong vật chất và gây ion hóa tao ra các cặp electron-ion, electron-l6 trong và các tia gamma thứ cấp.

Đôi với hệ do gamma, các cặp mang điện tạo ra do quá trình ion hóa được dùng đẻ đo lượng điện tích do tương tác tạo ra nhằm xác định năng lượng của gamma tới. Trong luận văn này dé cập đến việc ghi nhận các bức xạ gamma có năng lượng nằm trong dãy từ vài chục keV đến dưới 2 MeV chủ yếu xảy ra ba loại tương tác cơ bản là hap thụ quang điện, tán xa Compton và hiệu ứng tạo cặp. Khi đi qua vật chất, bức xạ gamma bị mắt năng lượng do hiệu ứng quang điện, tán xạ Compton bởi sự va chạm với electron trở nên vượt trội ở vùng năng lượng thấp và trung bình. vùng năng lượng cao là chiếm ưu thé khi đó bức xạ có thẻ tương tác với trường hạt nhân gây ra hiệu ứng tạo cặp.

Ngoài ra, các quá trình khác ít xảy ra hoặc không quan trọng trong vùng năng lượng được dé cập như tin xạ Thomson xảy ra với electron tự do, tán xa Rayleigh xảy ra với electron các lớp ngoài cùng, cả hai quá trình đều làm thay đôi hướng bay của photon nhưng không làm mất năng lượng, một số trường hợp photon tương tác với hạt nhân cũng sẽ không được quan tâm. Hiệu ứng quang điện Hiệu ứng quang điện là quá trình gamma tương tac với electron các lớp trong của nguyên tử và truyền toàn bộ năng lượng E = hv cho các electron này. Một phần năng lượng giúp electron thắng lực liên kết trên lớp thứ i (i= K, L, M.) với năng lượng liên kết là £, (£„ >£, >£„ >.) và một phần năng lượng trở thành động năng E, của electron bay ra khỏi lớp vỏ nguyên tử, các electron này gọi là các quang electron. Electron quang điện Hình 1.

Minh họa hiệu ứng quang điện 2 Theo định luật bảo toàn năng lượng và xung lượng, một phần nhỏ năng lượng cũng được truyền cho nguyên tử giật lùi và xem như không đáng kê nên electron thoát ra khỏi nguyên tử với năng lượng xấp xi theo công thức: E.2) Từ công thức (1.2) thi động năng cua electron là xác định và được hệ phô kế ghi nhận hình thành đính năng lượng toàn phan trong phd gamma. Các clectron trong nguyên tử hấp thụ hoàn toàn năng lượng gamma tới đề thoát ra ngoài đồng thời dé lại các lỗ trồng. Các electron từ các lớp khác trong nguyên tử sẽ chuyên vẻ lap day các lỗ trong và phát ra tia X đặc trưng. Ngoài ra, hiệu ứng Auger có thê xảy ra khi năng lượng tia X chuyên cho một electron khác trong cùng nguyên tử, làm birt electron này ra khỏi lớp vỏ nguyên tử và được gọi là electron Auger.

Electron Auger Hinh 1. Minh hoa sy phat electron Auger Do nang lượng liên kết thay đôi theo số nguyên tử Z nên tiết diện tương tác quang điện 6, không những phụ thuộc vào năng lượng gamma tới mà còn phụ thuộc vào số Z của vật liệu theo công thức [1]: 6, =const .8) 3 Theo hai công thức (1.8), đối với những vật liệu có Z lớn thì xác xuất xảy ra hiệu ứng quang điện là lớn ngay cả đối với những tia gamma có năng lượng cao. Điều này có ý nghĩa rất quan trọng trong việc tính toán khả năng che chắn và hấp thụ các tia gamma năng lượng cao. Đối với những vật liệu nhẹ thì hiệu ứng quang điện chủ yếu chỉ xây ra với những gamma có năng lượng thấp dong thời sự suy giảm mạnh tiết diện quang điện theo năng lượng gamma là lí do vì sao hiệu ứng quang điện là kênh trội của tương tác giữa gamma với vật chất ở vùng năng lượng tương đối thấp.

Tiết diện hiệu ứng quang điện phụ thuộc năng lượng gamma [1] Từ hình 1.3 cho ta thấy ở miền năng lượng gamma lớn tiết điện quang điện rất bé bởi khi đó gamma tới coi clectron có liên kết rất yếu trong nguyên tử. Khi năng lượng gamma bắt đầu giảm thì tiết diện tăng theo quy luật E† và khi E giảm dan đến £„ thì tiết điện tăng theo hàm E** va đạt cực đại tại E=e,. Nếu năng lượng gamma giảm xuống dưới giá trị ¢, thì hiệu ứng quang điện không thé xảy ra với các electron lớp K nên tiết điện giảm đột ngột. Tiết diện tăng trở lại khi năng lượng gamma tiếp tục giảm bởi lúc này xảy ra hiệu ứng quang điện đối với electron lớp L và đạt cực đại tại E=¢,, tiết điện lại giảm đột ngột khi E giảm xuống thấp hơn e,.

Và nếu năng lượng gamma tiếp tục giảm thì hiệu ứng quang điện sẽ xảy ra với các electron lớp M,. 4 * Hiệu ứng Compton Hiệu ứng Compton là quá trình tán xạ giữa gamma với các electron liên kết yếu trong nguyên tử trong đó gamma truyền một phan năng lượng cho electron và sau tán xạ phương bay của gamma bị lệch so với ban dau, tán xạ này gọi là tán xạ Compton. Sơ đô tán xạ gamma lên clectron tự do Theo định luật bảo toàn năng lượng và xung lượng ta tính được năng lượng gamma sau tán xạ, động năng giật lùi của electron và tìm được mối liên hệ giữa các góc tán xạ sau va chạm được trình bày theo các công thức (1. Năng lượng gamma sau tán xạ [12]: E E=—————— 1.12) với =——y; mic Từ công thức (1.9) cho thấy E: phụ thuộc vào góc tán xạ 0 nên trong phố gamma thu được nẻn liên tục với các xung mang năng lượng trải đài từ 0 đến năng lượng cực đại E-.

Ngoài ra E’ < E và phan năng lượng chênh lệch E- được truyền cho electron giật lùi. Nói cách khác thì năng lượng gamma tới có thê không được hấp thụ hoàn toàn trong lần tương tác dau tiên vì vậy dé theo đõi toàn bộ quá trình mat năng lượng của gamma cân xem xét các gamma thứ câp vả tương tác của chúng. Hiệu ứng Compton xảy ra mạnh ở vùng năng lượng trung bình từ vài chục keV đến vài MeV và phụ thuộc vào năng lượng liên kết của electron trong nguyên tử. Quá trình bat đầu xảy ra khi tăng năng lượng gamma đến giá trị lớn hơn nhiều so với năng lượng liên kết của các electron lớp trong nguyên tử khi đó có thể xem các electron này là các electron tự do và hiệu ứng quang điện không còn đáng kê.

Công thức tiết điện tin xa Compton toàn phan theo Klein-Nishina [12]: ; 2(1+ o = 2ar,” ive 2(1+4) a | 1l+2a _ 1n(1+2ø) a +L in(1+20)-—432.13) 2a (1+ 2a)" voi r= ` x= 2,82x10!Šm là bán kính electron cô điền. Phân bỏ tiết điện tán xa Compton theo góc tán xạ tính theo công thức Klein- Nishina với các giá trị năng lượng gamma tới từ IkeV đến 10MeV [7] 6 Xét 2 trường hợp giới hạn của tiết điện tán xạ Compton: Đối với các gamma có năng lượng thắp ( E<€ m,c”) tiết diện tán xa Compton 6, tăng tuyên tính khi năng lượng gamma giảm và dat xắp xi tiết điện tán xa Thomson được xác định theo công thức [1]: Gc = Gy {I~20+ Ea? +.14) : ÑT oy ue Q. trong đó Gy, = r, là tiết điện tấn xa Thomson. Đối với gamma có năng lượng cao (ơ >I hay E>>mc}) thi 6G, được xác định theo công thức { 1]: oc =m | 2 +In24, e a (1.15) Từ công thức (1.15) cho thay Ge biến thiên ti lệ nghịch với E ngoài ra trong nguyên tử có Z electron nên Ø, được xác định [1]: G,.16) Hiệu ứng tao cặp Hiệu ứng tạo cặp là quá trình gamma tương tac với trường điện từ quanh hạt nhân sinh ra cặp electron-positron và truyền toàn bộ năng lượng cho cặp electron-positron và nhân giật lùi.

Do đó, năng lượng gamma tới phải lớn hon hai lần năng lượng nghỉ của electron (E = 2m c”), trong thực tế xác suất dé hiệu ứng xảy đáng kê là vào khoảng vài MeV. Vậy sự tạo cặp chỉ chiếm ưu thé ở vùng năng lượng cao. Do năng lượng giật lùi là không đáng kề nên theo định luật bảo toàn năng lượng ta có tông động năng của electron và positron được xác định theo công thức: E,.17) Hai hat electron va positron chuyén dong cham dan trong vat liệu trong đó electron mat dan năng lượng đẻ ion hóa các nguyên tử môi trường, positron kết hợp với 7 electron khác gan đó xảy ra hiện tượng hủy cặp electron-positron và hai lượng tử gamma được phát ra gần như ngược chiêu nhau và mỗi lượng tử mang phan năng lượng là 0,511 MeV. y, OY \ Electron Positron : 0.511keV raed ee O.

Minh họa hiệu ứng tao cặp electron-positron Hiệu ứng tạo cặp là kênh trội ở vùng năng lượng cao nên khi xét trong miền năng lượng Sm.c’ <E<50m,c° thì tiết diện tạo cặp tỉ lệ với Z?. Với những vật liệu có số nguyên tử Z lớn thì tiết điện tạo cặp nằm trong miền năng lượng này là khá lớn. Ngoài ra nó còn phụ thuộc vào InE biéu diễn qua công thức [1]: Ø;„ = const .