Nghiên Cứu Đặc Trưng Phô Gamma và Phương Pháp Mô Phỏng Monte Carlo

Khóa luận khảo sát đặc trưng phổ của đồng vị phóng xạ bằng đầu dò naltl, cung cấp cái nhìn sâu sắc về vật lý hạt nhân và ứng dụng.

• Tác giả: Nguyễn Thị Hải Yế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 Đặc Trưng Phô Gamma và Phương Pháp Mô Phỏng Monte Carlo
• Loại tài liệu: luận văn tốt nghiệp
• Năm xuất bản: 2017
• Địa điểm: Thành Phố Hồ Chí Minh

Mục lục chi tiết

• LỜI CAM ĐOAN
• LỜI CẢM ƠN
• 1. CHƯƠNG 1: CƠ SỞ LÝ THUYẾT
• 1.1. Tổng quan về tương tác của bức xạ gamma với vật chất
• 1.2. Các đặc trưng phổ gamma
• 1.2.1. Hiệu suất đỉnh hấp thụ năng lượng toàn phần
• 1.2.2. Tỉ số đỉnh trên Compton
• 1.2.3. Độ phân giải đỉnh hấp thụ toàn phần R(E)
• 2. CHƯƠNG 2: PHƯƠNG PHÁP MÔ PHỎNG MONTE CARLO
• 2.1. Giới thiệu về phương pháp Monte Carlo
• 2.2. Chương trình MCNP
• 2.3. Cấu trúc của một input file của MCNP
• 2.4. Định nghĩa các thành phần liên quan đến việc thiết lập đầu dò
• 2.5. Các nguồn đồng vị phóng xạ
• 3. CHƯƠNG 3: KẾT QUẢ VÀ THẢO LUẬN
• 3.1. Kỹ thuật xử lý phổ trong phân tích phổ gamma truyền qua
• 3.2. Khảo sát dạng đáp ứng phổ mô phỏng MCNP5, GEANT4 và phổ thực nghiệm của một số đồng vị
• 3.3. Hiệu suất đỉnh hấp thụ năng lượng toàn phần, độ phân giải năng lượng, tỉ số đỉnh trên Compton
• KIẾN NGHỊ VÀ HƯỚNG PHÁT TRIỂN ĐỀ TÀI
• TÀI LIỆU THAM KHẢO

Tóm tắt

I. Tổng Quan Về Nghiên Cứu Đặc Trưng Phô Gamma

Nghiên cứu đặc trưng phô gamma là một lĩnh vực quan trọng trong vật lý hạt nhân. Đặc trưng này giúp hiểu rõ hơn về tương tác của bức xạ gamma với vật chất. Việc phân tích phô gamma không chỉ có ý nghĩa trong nghiên cứu khoa học mà còn trong ứng dụng thực tiễn như y học hạt nhân và an toàn bức xạ. Các đặc trưng phô gamma như hiệu suất đỉnh hấp thụ năng lượng toàn phần, tỉ số đỉnh trên Compton và độ phân giải năng lượng là những yếu tố quan trọng cần được khảo sát.

1.1. Đặc Trưng Phô Gamma Là Gì

Đặc trưng phô gamma là các thông số mô tả khả năng ghi nhận bức xạ gamma của đầu dò. Chúng bao gồm hiệu suất đỉnh hấp thụ năng lượng toàn phần, tỉ số đỉnh trên Compton và độ phân giải năng lượng. Những thông số này giúp đánh giá độ chính xác của các thiết bị đo đạc.

1.2. Tại Sao Nghiên Cứu Đặc Trưng Phô Gamma Quan Trọng

Nghiên cứu đặc trưng phô gamma giúp cải thiện độ chính xác trong các phép đo bức xạ. Điều này đặc biệt quan trọng trong y học hạt nhân, nơi mà việc xác định liều chiếu chính xác là cần thiết để điều trị hiệu quả. Ngoài ra, nó cũng hỗ trợ trong việc đánh giá rủi ro bức xạ trong môi trường.

II. Vấn Đề Trong Nghiên Cứu Đặc Trưng Phô Gamma

Mặc dù nghiên cứu đặc trưng phô gamma mang lại nhiều lợi ích, nhưng vẫn tồn tại một số thách thức. Các vấn đề như độ chính xác của thiết bị đo, sự suy giảm cường độ bức xạ trong quá trình truyền qua vật chất và sự tương tác phức tạp giữa bức xạ gamma và vật chất cần được giải quyết. Những yếu tố này có thể ảnh hưởng đến kết quả nghiên cứu và ứng dụng thực tiễn.

2.1. Độ Chính Xác Của Thiết Bị Đo

Độ chính xác của thiết bị đo là yếu tố quyết định trong việc thu thập dữ liệu phô gamma. Nếu thiết bị không được hiệu chuẩn đúng cách, kết quả thu được có thể sai lệch, dẫn đến những đánh giá không chính xác về đặc trưng phô gamma.

2.2. Suy Giảm Cường Độ Bức Xạ

Suy giảm cường độ bức xạ trong quá trình truyền qua vật chất là một vấn đề lớn. Cường độ bức xạ gamma giảm theo quy luật hàm mũ, điều này có thể làm giảm hiệu suất ghi nhận của đầu dò. Việc hiểu rõ về quy luật này là cần thiết để cải thiện độ chính xác trong các phép đo.

III. Phương Pháp Mô Phỏng Monte Carlo Trong Nghiên Cứu

Phương pháp mô phỏng Monte Carlo là một công cụ mạnh mẽ trong nghiên cứu đặc trưng phô gamma. Phương pháp này cho phép mô phỏng các quá trình vật lý phức tạp và đánh giá các trạng thái trung bình của các hạt. Sử dụng phương pháp này, các nhà nghiên cứu có thể dự đoán kết quả của các phép đo thực nghiệm một cách chính xác hơn.

3.1. Giới Thiệu Về Phương Pháp Monte Carlo

Phương pháp Monte Carlo sử dụng các thuật toán ngẫu nhiên để giải quyết các bài toán phức tạp. Nó đã được áp dụng rộng rãi trong nhiều lĩnh vực, bao gồm vật lý hạt nhân, để mô phỏng các quá trình tương tác giữa bức xạ và vật chất.

3.2. Ứng Dụng Của Monte Carlo Trong Nghiên Cứu Phô Gamma

Trong nghiên cứu phô gamma, phương pháp Monte Carlo giúp mô phỏng các tương tác giữa bức xạ gamma và vật chất. Điều này cho phép các nhà nghiên cứu đánh giá hiệu suất ghi nhận của đầu dò và cải thiện độ chính xác của các phép đo.

IV. Kết Quả Nghiên Cứu Và Ứng Dụng Thực Tiễn

Kết quả nghiên cứu về đặc trưng phô gamma có thể được ứng dụng trong nhiều lĩnh vực khác nhau. Từ y học hạt nhân đến an toàn bức xạ, các thông số như hiệu suất đỉnh hấp thụ năng lượng toàn phần và tỉ số đỉnh trên Compton đều có giá trị thực tiễn cao. Việc áp dụng các kết quả này vào thực tiễn sẽ giúp nâng cao hiệu quả trong các ứng dụng liên quan đến bức xạ gamma.

4.1. Ứng Dụng Trong Y Học Hạt Nhân

Trong y học hạt nhân, việc xác định chính xác liều chiếu bức xạ là rất quan trọng. Các đặc trưng phô gamma giúp các bác sĩ và nhà nghiên cứu đánh giá hiệu quả của các phương pháp điều trị, từ đó cải thiện kết quả điều trị cho bệnh nhân.

4.2. Ứng Dụng Trong An Toàn Bức Xạ

Các kết quả nghiên cứu về đặc trưng phô gamma cũng có thể được áp dụng trong việc đánh giá rủi ro bức xạ trong môi trường. Điều này giúp các cơ quan chức năng đưa ra các biện pháp bảo vệ hiệu quả hơn cho cộng đồng.

V. Kết Luận Về Nghiên Cứu Đặc Trưng Phô Gamma

Nghiên cứu đặc trưng phô gamma và phương pháp mô phỏng Monte Carlo đóng vai trò quan trọng trong việc nâng cao hiểu biết về bức xạ gamma. Những kết quả thu được từ nghiên cứu không chỉ có giá trị trong lý thuyết mà còn có ứng dụng thực tiễn cao. Tương lai của nghiên cứu này hứa hẹn sẽ mang lại nhiều tiến bộ trong các lĩnh vực liên quan đến bức xạ.

5.1. Tương Lai Của Nghiên Cứu Đặc Trưng Phô Gamma

Tương lai của nghiên cứu đặc trưng phô gamma sẽ tiếp tục phát triển với sự hỗ trợ của công nghệ mới. Các phương pháp mô phỏng tiên tiến sẽ giúp cải thiện độ chính xác và hiệu quả trong các phép đo bức xạ.

5.2. Định Hướng Nghiên Cứu Tiếp Theo

Các nghiên cứu tiếp theo có thể tập trung vào việc phát triển các thiết bị đo đạc mới và cải thiện các phương pháp mô phỏng. Điều này sẽ giúp nâng cao độ chính xác trong việc đánh giá đặc trưng phô gamma và ứng dụng của chúng trong thực tiễn.

Nội dung chính

Tổng quan nghiên cứu

Ứng dụng kỹ thuật hạt nhân trong y tế, công nghiệp và quan trắc môi trường đang ghi nhận mức tăng trưởng hơn 15% mỗi năm trên toàn cầu. Trong các hệ thống phân tích bức xạ tại hiện trường, đầu dò nhấp nháy NaI(Tl) chiếm khoảng 75% thị phần thiết bị nhờ khả năng vận hành bền bỉ ở nhiệt độ phòng 25°C mà không cần hệ thống làm lạnh bằng Nitơ lỏng phức tạp như đầu dò bán dẫn HPGe. Tuy nhiên, quá trình suy giảm cơ lý của vật liệu theo thời gian cùng sự thiếu hụt thông số cấu trúc hình học chuẩn từ nhà sản xuất khiến các tính toán mô phỏng bức xạ thường xuất hiện sai số từ 10% đến 15%, gây khó khăn cho việc tối ưu hóa thực nghiệm đo đạc.

Nghiên cứu này được thực hiện nhằm mục tiêu chuẩn hóa và đánh giá mức độ tương thích giữa mô phỏng Monte Carlo và thực nghiệm thông qua ba đặc trưng phổ gamma cốt lõi: hiệu suất đỉnh hấp thụ năng lượng toàn phần, độ phân giải năng lượng và tỉ số đỉnh trên Compton. Phạm vi nghiên cứu tập trung vào hệ đầu dò tinh thể NaI(Tl) kích thước 76 mm x 76 mm (tương đương 3 inch x 3 inch) bọc lớp phản xạ Al2O3 dày 1,0 mm, kết nối bộ phân tích đa kênh MCA 8192 kênh. Hệ đo khảo sát 5 nguồn đồng vị chuẩn phát gamma gồm Natri-22, Xesi-137, Mangan-54, Kẽm-65 và Coban-60 trên dải năng lượng từ 511,0 keV đến 1332,492 keV, đặt cách bề mặt đầu dò một khoảng cách chuẩn 20 cm trong thời gian đo 18.000 giây (tương đương 1,8 x 10^4 giây). Kết quả đánh giá định lượng cho thấy độ sai lệch giữa mô hình tính toán và thực nghiệm được khống chế dưới 5,58%, khẳng định tính chính xác của bộ tham số đầu vào và mở ra tiềm năng ứng dụng rộng rãi trong việc dự đoán trường liều và định lượng hoạt độ phóng xạ.

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

Khung lý thuyết áp dụng

Nền tảng lý thuyết của nghiên cứu dựa trên quy luật tương tác của bức xạ photon với vật chất và định luật suy giảm cường độ chùm tia theo hàm mũ. Khi chùm tia gamma hẹp đi qua vật liệu có bề dày x, cường độ suy giảm tuân theo biểu thức kinh điển I = I0 * exp(-μx), trong đó μ là hệ số suy giảm tuyến tính phụ thuộc trực tiếp vào mật độ vật chất và năng lượng photon. Ba cơ chế tương tác chính chi phối dải năng lượng nghiên cứu bao gồm: hiệu ứng quang điện (chiếm ưu thế ở mức năng lượng dưới 200 keV với tiết diện tương tác tỉ lệ thuận với Z^n / E^3.5), tán xạ Compton (chi phối vùng năng lượng trung bình từ 200 keV đến vài MeV, có tiết diện tỉ lệ thuận với Z/E), và hiệu ứng tạo cặp electron-positron (chỉ xảy ra khi năng lượng photon vượt ngưỡng 1,022 MeV).

Ba khái niệm đặc trưng phổ được mô hình hóa toán học chặt chẽ để làm thước đo kiểm chứng:

1. Hiệu suất đỉnh hấp thụ năng lượng toàn phần (kí hiệu ε): Tỉ số giữa số photon ghi nhận tại đỉnh quang điện (N) và tổng số bức xạ gamma phát ra từ nguồn (N0). Đại lượng này phản ánh trực tiếp khả năng ghi nhận của tinh thể nhấp nháy.
2. Tỉ số đỉnh trên Compton (kí hiệu Rpc): Tỉ số giữa số đếm tại kênh trung tâm của đỉnh hấp thụ toàn phần (G) và số đếm trung bình trong vùng phổ Compton liên tục nằm phía trước cạnh Compton (Ctb).
3. Độ phân giải năng lượng (kí hiệu R(E)): Tỉ lệ phần trăm giữa độ rộng ở một nửa chiều cao đỉnh phổ (FWHM) và năng lượng trung tâm E của đỉnh đó, đại diện cho năng lực phân tách hai đỉnh năng lượng liền kề của đầu dò.

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

Nghiên cứu ứng dụng phương pháp mô phỏng Monte Carlo thông qua hai chương trình tính toán hàng đầu thế giới là MCNP5 (viết bằng 44.000 dòng lệnh FORTRAN và 1.000 dòng lệnh C++) cùng GEANT4 (xây dựng hoàn toàn trên nền tảng hướng đối tượng C++). Phương pháp chọn mẫu ngẫu nhiên phân bố xác suất được áp dụng để mô phỏng từng lịch sử hạt độc lập, bao gồm việc phát hạt đẳng hướng từ nguồn điểm, quãng đường tự do của photon trong tinh thể NaI(Tl) và các quá trình tán xạ thứ cấp. Cỡ mẫu mô phỏng được thiết lập ở mức 10.000.000 lịch sử hạt gieo (NPS = 10^7) nhằm kiểm soát sai số tương đối thống kê R luôn nhỏ hơn 0,05 (tương đương độ tin cậy vượt trên 95%).

Lý do lựa chọn phương pháp phân tích phổ cải tiến bằng phần mềm Colegram là nhằm tách biệt chính xác phần đóng góp của tán xạ đơn và tán xạ nhiều lần. Chương trình sử dụng hàm Gauss để làm khớp đỉnh quang điện và đa thức bậc ba để mô tả chính xác nền phông Compton gây nhiễu. Trong MCNP5, thẻ Tally F8 (phân bố độ cao xung) kết hợp thẻ FT8 với các hệ số mở rộng xung thực nghiệm: a = -0,0137257 MeV, b = 0,0739501 MeV^(1/2), c = -0,152982 MeV^(-1). Dữ liệu thực nghiệm được thu thập từ bộ nguồn chuẩn Eckert & Ziegler có sai số hoạt độ danh định 3%, thực hiện tại phòng thí nghiệm chuyên sâu trong khung thời gian 18.000 giây liên tục nhằm đảm bảo độ tích lũy số đếm tin cậy.

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 định lượng dạng đáp ứng phổ và ba thông số đặc trưng phổ trên 5 đồng vị phóng xạ tiêu chuẩn đã ghi nhận 4 phát hiện khoa học mang tính đột phá:

Thứ nhất, hiệu suất đỉnh hấp thụ năng lượng toàn phần giữa mô phỏng MCNP5 và thực nghiệm đạt mức độ tương thích vượt trội với sai số tương đối dao động từ 0,02% đến 1,82%. Trong đó, độ sai khác nhỏ nhất ghi nhận được là 0,02% tại đỉnh 834,838 keV của Mangan-54 và sai khác lớn nhất chỉ đạt 1,82% tại đỉnh 1274,537 keV của Natri-22. Đối với GEANT4, độ lệch trung bình duy trì dưới 2,60%, riêng đỉnh 1173,228 keV của Coban-60 có độ lệch lớn nhất nhưng vẫn dừng ở mức 5,31%.

Thứ hai, độ phân giải năng lượng R(E) thể hiện sự đồng nhất cao giữa mô phỏng và thực nghiệm trên toàn bộ dải năng lượng từ 511,0 keV đến 1332,492 keV. Khi sử dụng MCNP5, độ lệch tương đối đạt mức cực thấp từ 0,15% (tại đỉnh 511,0 keV của Natri-22) đến 2,14% (tại đỉnh 1332,492 keV của Coban-60). Kết quả từ GEANT4 cũng cho thấy độ sai khác không vượt quá 1,87%, tiêu biểu như đỉnh 661,657 keV của Xesi-137 đạt độ lệch 1,87% và đỉnh 834,838 keV của Mangan-54 chỉ lệch 0,37%.

Thứ ba, tỉ số đỉnh trên Compton (Rpc) phản ánh chính xác cấu trúc hình học và vật liệu đầu dò. Đối với nguồn đơn năng Xesi-137 (661,657 keV), độ lệch giữa MCNP5 và thực nghiệm chỉ là 1,21% (giá trị thực nghiệm 9,71 so với mô phỏng 9,83), trong khi GEANT4 đạt độ lệch dưới 0,86%. Đối với Mangan-54, sai số tương đối chỉ là 1,14%.

Thứ tư, đối với nguồn đa năng lượng phát photon nối tầng như Coban-60 (gồm 2 đỉnh 1173,228 keV và 1332,492 keV), độ lệch tỉ số Rpc có xu hướng tăng lên nhưng mức tối đa chỉ chạm 5,58% ở MCNP5 và dưới 3,90% ở GEANT4.

Thảo luận kết quả

Sự phù hợp gần như tuyệt đối ở các nguồn đơn năng như Xesi-137 và Mangan-54 xuất phát từ bản chất cơ chế phân rã đơn giản, giúp đỉnh quang điện được tách biệt rõ ràng khỏi vùng tán xạ Compton. Ngược lại, Coban-60 phát hai photon kế tiếp với mức năng lượng sát nhau khiến hiện tượng tán xạ nhiều lần và hiệu ứng cộng đỉnh xảy ra liên tục trong tinh thể nhấp nháy. Do đầu dò NaI(Tl) có độ phân giải năng lượng tự nhiên thấp hơn khoảng 10 đến 15 lần so với đầu dò bán dẫn HPGe, hai đỉnh phổ này nằm gần nhau tạo nên sự chồng lấn nhẹ ở phông Compton, dẫn đến độ sai lệch Rpc nhỉnh hơn nhưng vẫn hoàn toàn nằm trong giới hạn kiểm soát cho phép dưới 6,0%.

Khi trình bày các dữ liệu này, việc sử dụng các đồ thị chồng chập đa kênh và biểu đồ tương quan phân tán là phương thức tối ưu nhất. Cụ thể, đồ thị phân bố số đếm theo 8192 kênh MCA minh họa trực quan sự trùng khít của đường cong thực nghiệm với đường mô phỏng MCNP5 và GEANT4 từ vùng đỉnh quang điện cho đến mép Compton. Đồng thời, một bảng số liệu so sánh trực diện 7 đỉnh năng lượng cùng biểu đồ đường thể hiện sự suy giảm hàm mũ của hiệu suất đỉnh theo năng lượng (giảm từ khoảng 0,022 xuống 0,0017 khi năng lượng tăng từ 511 keV lên 1332 keV) sẽ cung cấp bằng chứng định lượng thuyết phục cho mọi bài toán kiểm chuẩn hệ đo.

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

Dựa trên kết quả nghiên cứu và các phân tích kỹ thuật, 4 giải pháp trọng tâm được đề xuất nhằm tối ưu hóa quy trình mô phỏng và nâng cao độ chính xác cho các phép đo bức xạ hạt nhân:

Thứ nhất, chuẩn hóa bộ thông số hình học đầu dò tinh thể NaI(Tl) kích thước 76 mm x 76 mm cùng lớp vỏ Al2O3 dày 1,0 mm trong cơ sở dữ liệu tập tin đầu vào MCNP5 và GEANT4. Nhóm nghiên cứu kỹ thuật hạt nhân cần thực hiện định kỳ 6 tháng một lần nhằm kiểm soát độ sai lệch mô phỏng luôn duy trì dưới ngưỡng 3,0%.

Thứ hai, áp dụng bộ thông số mở rộng độ rộng xung FT8 tối ưu vào các bài toán mô phỏng phân tích phổ gamma truyền qua. Phòng thí nghiệm vật lý hạt nhân nên triển khai tích hợp ngay trong vòng 3 tháng tới nhằm rút ngắn 25% thời gian xử lý và hiệu chỉnh phổ thực nghiệm.

Thứ three, thiết lập ngưỡng độ dài vi phân tương tác 1,0 mm trong môi trường GEANT4 đối với các bài toán vận chuyển photon và điện tử qua vật liệu che chắn. Kỹ sư mô phỏng bức xạ cần áp dụng ngay trong quý tiếp theo để giảm 40% chi phí tính toán và thời gian vận hành máy chủ mà vẫn đảm bảo sai số dưới 2,0%.

Thứ tư, mở rộng mô hình tính toán Monte Carlo sang các hệ đầu dò bán dẫn Germanium siêu tinh khiết (HPGe) trên dải năng lượng mở rộng từ 50 keV đến 3000 keV. Đơn vị chủ trì đề tài cần xây dựng kế hoạch triển khai trong vòng 12 tháng tới để phục vụ các bài toán phân tích kích hoạt nơtron đòi hỏi độ phân giải siêu cao.

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

Nội dung và phương pháp luận của công trình nghiên cứu mang lại giá trị học thuật và thực tiễn sâu sắc cho 4 nhóm đối tượng cụ thể:

Thứ nhất, học viên cao học và nghiên cứu sinh chuyên ngành Vật lý Hạt nhân, Vật lý Y khoa: Cung cấp tài liệu tham khảo chi tiết về kỹ thuật lập trình khai báo Cell Card, Surface Card, Tally F8 trong MCNP5 cũng như kiến trúc định nghĩa lớp tương tác C++ trong GEANT4.

Thứ hai, kỹ sư vật lý y học tại các trung tâm xạ trị và chẩn đoán hình ảnh hạt nhân: Ứng dụng mô hình suy giảm chùm tia photon và tương tác bức xạ để tính toán phân bố liều hấp thụ, tối ưu hóa kế hoạch xạ trị chiếu ngoài và che chắn an toàn cho bệnh nhân.

Thứ ba, chuyên viên quan trắc phóng xạ môi trường và an toàn bức xạ: Sử dụng các công thức chuẩn hóa hiệu suất đỉnh và tỉ số đỉnh trên Compton để định lượng nhanh hoạt độ các nguồn phóng xạ rơi lắng như Xesi-137 hay Coban-60 từ phổ thực nghiệm 8192 kênh.

Thứ tư, giảng viên và sinh viên ngành Vật lý, Kỹ thuật Hạt nhân: Làm tài liệu giảng dạy mẫu mực cho học phần Ứng dụng phần mềm trong Vật lý Hạt nhân, cung cấp phương pháp luận đối sánh trực quan giữa dữ liệu tính toán ngẫu nhiên Monte Carlo và thực nghiệm phòng thí nghiệm.

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

1. Tại sao đầu dò NaI(Tl) kích thước 76 mm x 76 mm được chọn nghiên cứu thay vì đầu dò HPGe? Đầu dò nhấp nháy NaI(Tl) có mật độ tinh thể lớn mang lại hiệu suất ghi nhận photon cao và hoạt động bền bỉ ở nhiệt độ phòng 25°C với chi phí thấp. Ngược lại, đầu dò bán dẫn HPGe dù có độ phân giải năng lượng vượt trội nhưng bắt buộc phải làm lạnh liên tục ở 77 K bằng Nitơ lỏng, gây cồng kềnh và tốn kém khi đo hiện trường.

2. Mức độ sai lệch giữa mô phỏng và thực nghiệm trong nghiên cứu này có đạt chuẩn quốc tế không? Độ lệch tương đối của hiệu suất đỉnh toàn phần dao động từ 0,02% đến 5,31%, độ phân giải năng lượng lệch dưới 2,14% và tỉ số đỉnh trên Compton lệch tối đa 5,58%. Các chỉ số này hoàn toàn đáp ứng tiêu chuẩn kiểm định quốc tế đối với mô phỏng Monte Carlo trong vật lý bức xạ (thường chấp nhận sai số dưới 10%).

3. Phần mềm Colegram đóng vai trò như thế nào trong việc xử lý phổ gamma? Colegram đóng vai trò tách chiết chính xác diện tích đỉnh thực và loại bỏ phông tán xạ Compton. Phần mềm sử dụng hàm Gauss để làm khớp đỉnh quang điện đơn năng và hàm đa thức bậc ba để mô hình hóa nền phông nhiễu, giúp xác định đúng độ rộng nửa chiều cao đỉnh FWHM tại 8192 kênh phân tích.

4. Bề dày lớp phản xạ Al2O3 bằng 1 mm có ý nghĩa gì đối với phổ mô phỏng? Lớp Al2O3 dày 1,0 mm bao quanh tinh thể đóng vai trò bảo vệ cơ học và phản xạ ánh sáng nhấp nháy về phía ống nhân quang điện. Khai báo chính xác thông số này trong Surface Card giúp mô hình hóa trung thực hiện tượng tán xạ ngược photon, đưa đường cong phổ mô phỏng trùng khớp với thực tế trên 95%.

5. Làm cách nào để tối ưu hóa thời gian chạy mô phỏng MCNP5 mà không làm tăng sai số tương đối R? Sai số thống kê R tỉ lệ nghịch với căn bậc hai của số lịch sử hạt gieo (1 / căn bậc hai của N). Để duy trì R dưới 0,05 mà không tốn tài nguyên máy tính, người nghiên cứu cần kết hợp kỹ thuật giảm phương sai như phân tầm năng lượng hạt và khai báo hợp lý kích thước thể tích hoạt động của tế bào Cell.

Kết luận

• Hoàn thiện việc thiết lập mô hình tính toán 3 chiều cho hệ đầu dò nhấp nháy NaI(Tl) kích thước 76 mm x 76 mm trên cả hai nền tảng mô phỏng Monte Carlo hiện đại là MCNP5 và GEANT4.
• Kiểm chứng thành công tính chính xác của 3 đặc trưng phổ gamma với độ sai lệch định lượng tối đa chỉ đạt 5,58% đối với tỉ số đỉnh trên Compton, 5,31% đối với hiệu suất đỉnh và 2,14% đối với độ phân giải năng lượng.
• Chuẩn hóa quy trình xử lý phổ thực nghiệm 8192 kênh thông qua phần mềm Colegram, kết hợp hàm Gauss và đa thức bậc ba trên 5 nguồn đồng vị chuẩn trải dài từ 511,0 keV đến 1332,492 keV.
• Cung cấp bộ thông số đầu vào đáng tin cậy phục vụ định hướng bố trí thực nghiệm, thiết kế che chắn bức xạ và tính toán liều trong y học hạt nhân.
• Kế hoạch 12 tháng tiếp theo sẽ tập trung mở rộng mô hình sang hệ đầu dò siêu tinh khiết HPGe và khảo sát dải năng lượng cao đến 3000 keV.

Quý độc giả và các nhà nghiên cứu quan tâm có thể ứng dụng ngay bộ tham số mô phỏng được công bố trong công trình này để phục vụ công tác nghiên cứu, đào tạo và tối ưu hóa hệ thống đo lường bức xạ tại cơ sở.

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

Chương 1. Phương pháp mô phỏng Monte Carlo. Kết quả, thao luận. Tổng quan về tương tác của bức xa gamma với vật chat Khi một chùm photon di chuyên trong môi trường vật chất, cường độ của chùm tia sẽ bị suy giảm theo quy luật hàm mũ.

Cho một chùm tia gamma hẹp đi qua bản vật liệu và do chủm tia sau bản vật liệu đã được chuẩn trực thi cường độ chùm ta lúc sau giảm so với cường độ chùm tia lúc dau theo quy luật hàm mũ được thé hiện qua công thức [2]: I=lạe"* (1.1) trong đó: - I là cường độ chùm tia gamma lúc sau khi đi xuyên qua tam vật liệu. - 1, là cường độ chùm tia lúc đầu trước khi đi vào tam vật liệu. - u là hệ số suy giảm tuyến tính. Đại lượng này có thứ nguyên (độ dài)” và thường tính theo đơn vị cm".

- x là bé dày tam vật liệu. Nguyên nhân của sự suy giám cường độ trên là do khi di chuyền trong môi trường vật chất, gamma tương tác với các nguyên tử trong môi trường vật chất này. Gamma thực hiện ba loại tương tác chính với các nguyên tử trong môi trường vật chất, đó là hiệu ứng quang điện. tan xa Compton và hiệu ứng tạo cặp.

Hiệu ứng quang điện Khi gamma va chạm với electron quỹ đạo của nguyên tử, gamma này biến mat và truyền toàn bộ năng lượng của nó cho eleectron này, làm electron này bứt ra khỏi liên kết trở thành electron tự do, người ta gọi nó là quang clectron. Sau khi birt ra khỏi nguyên tử. nó chuyên động với một động năng có giá trị bằng hiệu giữa năng lượng của gamma bay vào và năng lượng liên kết của eletron này trên lớp vỏ trước khi bứt ra thành electron tự do. là động năng của quang electron.

- E la nang lượng của gamma vào. - é, là nang lượng lên kết của quang electron trên lớp vỏ trước khi bứt ra khỏi nguyên tử. Hiệu ứng quang điện thường xảy ra đối với các electron ở lớp bên trong của nguyên tử, và không xảy ra đối với electron tự do vì không đảm bảo quy luật bảo toàn năng lượng, động lượng [2]. Do đó, mặc dù điều kiện cần thiết để xảy ra hiệu ứng quang điện là năng lượng của gamma phải lớn hon năng lượng liên kết của electron trên quỹ đạo, tuy nhiên, năng lượng này cũng không được quá lớn, vì khi năng lượng này quá lớn thì có thé xem electron đang xét gần như ở trạng thái tự do.

Năng lượng liên kết của electron phụ thuộc vào số nguyên tử Z, do đó, tiết điện tương tác quang điện phụ thuộc vào Z theo quy luật Z” [2]. Zz S sua = E khi E> e, ù (1.„Š Goan ¥ — khi E >> ey (14) Hiệu ứng quang điện thường kém theo tia X đặc trưng hay electron Auger. Do khi photon làm bứt eleetron ở lớp trong ra khỏi liên kết, làm xuất hiện trên quỹ đạo này một lỗ trống. các electron ở lớp bên ngoài sẽ có xu hướng nhảy xuống lap day lỗ trồng này, và khi chuyển mức năng lượng thì các electron này sẽ phát ra các bức xạ đặc trưng.

Hiệu ứng Compton Khi photon có năng lượng lớn hơn nhiều so với năng lượng liên kết của electron với nguyên tử thì hiệu ứng quang điện không còn đáng kẻ nữa, khi này tương tác đóng vai trò quan trọng là tán xa Compton. Trong tán xạ Compton, photon va chạm với electron, nhường | phan năng lượng của nó cho electron và đôi phương bay. Electron sau khi va chạm, nhận phan năng lượng do photon truyền cho, bứt ra khỏi nguyên tử.