Khóa Luận Tốt Nghiệp Về Mô Phỏng Phổ Tia Gamma Sử Dụng Phương Pháp Monte Carlo

Luận văn tốt nghiệp vật lý nghiên cứu tốt nghiệp vật lý mô phỏng phổ tia gamma bằng phương pháp monte carlo, điều tra thực trạng, phân tích số liệu, đề xuất biện pháp cải tiến

• Tác giả: Nguyễn Thị Bích Loan
• Người hướng dẫn: ThS. Tạ Hưng Quý
• Trường học: Trường Đại Học Khoa Học Tự Nhiên
• Chuyên ngành: Vật Lý
• Đề tài: Khóa Luận Tốt Nghiệp Vật Lý Mô Phỏng Phổ Tia Gamma 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: 2000 - 2004
• Địa điểm: Thành Phố Hồ Chí Minh

Mục lục chi tiết

• LỜI MỞ ĐẦU
• 1. CHƯƠNG 1: PHƯƠNG PHÁP MONTE - CARLO
• 1.1. Tổng Quan Về Phương Pháp Mô Phỏng Monte — Carlo
• 1.2. Số ngẫu nhiên - giả ngẫu nhiên trên máy tính
• 1.3. Phương pháp tạo đồng dư tuyến tính (Linear Congruential)
• 1.4. Phương pháp hỗn hợp (Mixed Generators)
• 1.5. Phương pháp nhân tử (Multiplicative Generators)
• 1.6. Các phương pháp tạo số ngẫu nhiên trong mô phỏng Monte - Carlo
• 1.7. Phương pháp dùng hàm ngược
• 1.8. Phương pháp kết hợp
• 1.9. Phương pháp xoắn lại (Convolution)
• 1.10. Phương pháp chấp nhận loại bỏ (Acceptance-Rejection)
• 1.11. Một số ứng dụng của phương pháp mô phỏng Monte — Carlo
• 2. CHƯƠNG 2: SỰ TRUYỀN LƯỢNG TIA GAMMA QUA MÔI TRƯỜNG VẬT CHẤT
• 2.1. Tính chất của bức xạ Gamma
• 2.2. Cơ cấu dịch chuyển và sự hấp thụ bức xạ Gamma
• 2.3. Tán xạ Gamma
• 2.4. Hệ số hấp thụ (Hệ số suy giảm tuyến tính)
• 3. CHƯƠNG 3: PHỔ NĂNG LƯỢNG BỨC XẠ GAMMA DO DETECTOR NHẤP NHÁY GHI NHẬN
• 3.1. Nguyên tắc hoạt động của Detector nhấp nháy
• 3.2. Hiệu suất của detector nhấp nháy
• 3.3. Chất nhấp nháy, tính chất chung của chất nhấp nháy và nguyên tắc ghi nhận bức xạ nhờ chất nhấp nháy
• 3.4. Detector nhấp nháy dùng tinh thể NaI(Tl)
• 3.5. Phổ do Detector Nhấp Nháy Ghi Nhận
• 3.5.1. Sự tạo thành phổ biên độ xung
• 3.5.2. Nguồn gốc các đỉnh phổ
• 3.5.3. Một số dạng phổ minh họa
• 4. PHẦN II: BÀI TOÁN MÔ PHỎNG BỨC XẠ GAMMA QUA VẬT CHẤT
• 4.1. Lựa chọn ngôn ngữ lập trình mô phỏng
• 4.2. Trình tự giải quyết
• 4.3. Xây dựng thuật toán cho chương trình
• 4.4. Chuyển sang ngôn ngữ lập trình Matlab

Tóm tắt

I. Tổng quan về khóa luận tốt nghiệp vật lý mô phỏng tia gamma

Khóa luận tốt nghiệp về vật lý mô phỏng tia gamma bằng phương pháp Monte Carlo là một chủ đề nghiên cứu quan trọng trong lĩnh vực vật lý hạt nhân. Mục tiêu chính của nghiên cứu này là áp dụng phương pháp mô phỏng để phân tích và hiểu rõ hơn về các quá trình tương tác của tia gamma với vật chất. Việc sử dụng phương pháp Monte Carlo giúp giảm thiểu chi phí thí nghiệm và tăng cường độ chính xác trong việc dự đoán các kết quả thực nghiệm.

1.1. Khái niệm về tia gamma và ứng dụng trong vật lý

Tia gamma là loại bức xạ điện từ có năng lượng cao, thường được phát ra từ các quá trình hạt nhân. Chúng có ứng dụng rộng rãi trong y học, công nghiệp và nghiên cứu khoa học. Việc hiểu rõ về tia gamma là cần thiết để phát triển các thiết bị an toàn và hiệu quả.

1.2. Phương pháp Monte Carlo trong mô phỏng vật lý

Phương pháp Monte Carlo là một kỹ thuật mô phỏng dựa trên việc sử dụng các số ngẫu nhiên để giải quyết các bài toán phức tạp. Kỹ thuật này cho phép mô phỏng các quá trình vật lý một cách chính xác và hiệu quả, đặc biệt trong việc phân tích tương tác của tia gamma với vật chất.

II. Thách thức trong mô phỏng tia gamma bằng phương pháp Monte Carlo

Mặc dù phương pháp Monte Carlo mang lại nhiều lợi ích, nhưng cũng tồn tại một số thách thức trong việc mô phỏng tia gamma. Các vấn đề như độ chính xác của số ngẫu nhiên, thời gian tính toán và khả năng xử lý dữ liệu lớn là những yếu tố cần được xem xét kỹ lưỡng.

2.1. Độ chính xác của số ngẫu nhiên trong mô phỏng

Độ chính xác của các số ngẫu nhiên được sử dụng trong mô phỏng có ảnh hưởng lớn đến kết quả cuối cùng. Việc lựa chọn phương pháp tạo số ngẫu nhiên phù hợp là rất quan trọng để đảm bảo tính chính xác của mô phỏng.

2.2. Thời gian tính toán và hiệu suất mô phỏng

Thời gian tính toán là một yếu tố quan trọng trong mô phỏng Monte Carlo. Việc tối ưu hóa thuật toán và sử dụng các công nghệ tính toán hiện đại có thể giúp giảm thiểu thời gian và nâng cao hiệu suất mô phỏng.

III. Phương pháp mô phỏng tia gamma bằng Monte Carlo

Để thực hiện mô phỏng tia gamma, cần xây dựng một mô hình toán học chính xác và áp dụng các thuật toán phù hợp. Phương pháp Monte Carlo cho phép mô phỏng các quá trình tương tác phức tạp giữa tia gamma và vật chất một cách hiệu quả.

3.1. Xây dựng mô hình toán học cho mô phỏng

Mô hình toán học là nền tảng cho việc mô phỏng. Cần xác định các tham số vật lý và điều kiện ban đầu để đảm bảo tính chính xác của mô phỏng. Việc xây dựng mô hình cần dựa trên các lý thuyết vật lý hiện có.

3.2. Thuật toán mô phỏng và triển khai trên máy tính

Sau khi xây dựng mô hình, cần phát triển thuật toán mô phỏng và triển khai trên máy tính. Việc sử dụng ngôn ngữ lập trình như Matlab giúp thực hiện các phép tính phức tạp một cách nhanh chóng và hiệu quả.

IV. Ứng dụng thực tiễn của mô phỏng tia gamma

Mô phỏng tia gamma bằng phương pháp Monte Carlo có nhiều ứng dụng thực tiễn trong các lĩnh vực như y học, an toàn bức xạ và nghiên cứu vật lý. Kết quả từ mô phỏng có thể giúp cải thiện thiết kế thiết bị và quy trình làm việc trong các lĩnh vực này.

4.1. Ứng dụng trong y học và điều trị ung thư

Mô phỏng tia gamma có thể được sử dụng để tối ưu hóa liều lượng bức xạ trong điều trị ung thư. Việc hiểu rõ về tương tác của tia gamma với mô tế bào giúp cải thiện hiệu quả điều trị và giảm thiểu tác dụng phụ.

4.2. An toàn bức xạ trong công nghiệp

Trong công nghiệp, mô phỏng tia gamma giúp đánh giá rủi ro và thiết kế các biện pháp an toàn cho người lao động. Việc sử dụng mô phỏng giúp giảm thiểu tiếp xúc với bức xạ và bảo vệ sức khỏe cho nhân viên.

V. Kết luận và triển vọng tương lai của nghiên cứu

Khóa luận tốt nghiệp về mô phỏng tia gamma bằng phương pháp Monte Carlo không chỉ cung cấp kiến thức lý thuyết mà còn mở ra nhiều cơ hội nghiên cứu mới. Tương lai của lĩnh vực này hứa hẹn sẽ có nhiều tiến bộ nhờ vào sự phát triển của công nghệ và phương pháp mô phỏng.

5.1. Tương lai của phương pháp Monte Carlo trong vật lý

Phương pháp Monte Carlo sẽ tiếp tục phát triển và được áp dụng rộng rãi trong nhiều lĩnh vực khác nhau. Sự kết hợp giữa công nghệ máy tính và lý thuyết vật lý sẽ tạo ra những bước tiến mới trong nghiên cứu.

5.2. Khuyến nghị cho nghiên cứu tiếp theo

Nghiên cứu tiếp theo nên tập trung vào việc cải thiện độ chính xác của mô phỏng và mở rộng ứng dụng của phương pháp Monte Carlo trong các lĩnh vực khác nhau. Việc hợp tác giữa các nhà nghiên cứu và ngành công nghiệp cũng là điều cần thiết để thúc đẩy sự phát triển.

Nội dung chính

Tổng quan nghiên cứu

Nghiên cứu tương tác giữa bức xạ gamma và môi trường vật chất đóng vai trò then chốt trong vật lý hạt nhân hiện đại, công nghệ y học phóng xạ và thiết kế an toàn bức xạ công nghiệp. Bức xạ gamma là sóng điện từ mang năng lượng cao, dao động từ khoảng 10 keV đến hơn 20 GeV, có khả năng đâm xuyên mạnh và gây nguy cơ phơi nhiễm bức xạ nghiêm trọng cho con người nếu không được kiểm soát chặt chẽ. Trong thực tế, việc thực hiện các phép đo thực nghiệm phổ kế gamma đòi hỏi chi phí đầu tư trang thiết bị đo lường đắt đỏ, đồng thời tiềm ẩn rủi ro chiếu xạ đối với nhân viên phòng thí nghiệm. Phương pháp giải tích truyền thống khó có thể giải quyết trọn vẹn các bài toán vận chuyển bức xạ đa chiều do tính chất tương tác ngẫu nhiên và phức tạp của photon trong vật chất.

Mục tiêu trọng tâm của nghiên cứu là ứng dụng phương pháp mô phỏng Monte-Carlo để mô hình hóa toàn diện quá trình truyền bức xạ gamma qua môi trường vật chất và tái lập phổ năng lượng ghi nhận bởi detector nhấp nháy tinh thể NaI(Tl). Đề tài tập trung khảo sát các nguồn phóng xạ chuẩn điển hình gồm Cs-137 phát photon 0,662 MeV và Co-60 phát hai mức năng lượng 1,173 MeV và 1,332 MeV. Toàn bộ chương trình tính toán được phát triển trên nền tảng ngôn ngữ lập trình khoa học và đối chiếu trực tiếp với hệ đo thực nghiệm tại Phòng thí nghiệm Bộ môn Vật lý Hạt nhân, Trường Đại học Khoa học Tự nhiên Thành phố Hồ Chí Minh. Kết quả nghiên cứu mang ý nghĩa thực tiễn to lớn khi giúp giảm thiểu hơn 70% chi phí thử nghiệm thực tế, đạt độ tương thích phổ với sai số tương đối dưới 5%, đồng thời cung cấp giải pháp an toàn tuyệt đối phục vụ công tác đào tạo và nghiên cứu chuyên sâu.

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 điện động lực học lượng tử và vật lý bức xạ hạt nhân, tập trung vào bốn cơ chế tương tác cơ bản giữa photon gamma và vật chất:

1. Hiệu ứng quang điện: Quá trình photon bị hấp thụ hoàn toàn bởi electron liên kết trong nguyên tử, làm bứt electron ra ngoài với động năng bằng hiệu giữa năng lượng photon tới và năng lượng liên kết. Tiết diện vi mô của hiệu ứng tỉ lệ thuận với bậc số nguyên tử $Z^5$ và giảm nhanh theo hàm bậc $E_\gamma^{-3.5}$, đóng vai trò chủ đạo ở vùng năng lượng dưới 100 keV đối với các vật liệu có số khối lớn.
2. Tán xạ Compton: Quá trình photon va chạm đàn hồi với electron tự do hoặc electron liên kết yếu, làm photon bị lệch hướng một góc $\theta$ từ 0 độ đến 180 độ và truyền một phần động năng cho electron giật lùi. Tiết diện vi phân tán xạ được xác định chính xác theo công thức Klein-Nishina, chi phối dải năng lượng trung bình từ 0,2 MeV đến vài MeV.
3. Hiệu ứng tạo cặp: Hiện tượng photon năng lượng cao triệt tiêu trong điện trường hạt nhân để tạo ra cặp electron và positron. Quá trình này chỉ xảy ra khi năng lượng photon vượt qua ngưỡng năng lượng tối thiểu 1,022 MeV (gấp đôi năng lượng nghỉ của electron $2m_0c^2$).
4. Tán xạ Rayleigh: Hiện tượng tán xạ đàn hồi trên toàn bộ nguyên tử, trong đó photon giữ nguyên tần số và năng lượng mà chỉ thay đổi hướng truyền.

Cơ chế ghi nhận bức xạ sử dụng detector nhấp nháy NaI(Tl) kết hợp ống nhân quang điện (PMT). Tinh thể NaI pha tạp Thallium với nồng độ tối ưu khoảng 1% đến 2%, sở hữu khối lượng riêng $3,67\text{ g/cm}^3$, số nguyên tử hiệu dụng trung bình xấp xỉ 50, thời gian phát quang cực ngắn chỉ $0,25,\mu\text{s}$ ($2,5 \times 10^{-7}\text{ s}$) và sản lượng ánh sáng kỹ thuật đạt 20% đến 25%.

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

Nghiên cứu ứng dụng phương pháp thống kê Monte-Carlo để mô phỏng ngẫu nhiên từng lịch sử chuyển động của hạt photon. Thuật toán phát sinh chuỗi số giả ngẫu nhiên dựa trên phương pháp tạo đồng dư tuyến tính theo công thức Lehmer, đảm bảo chu kỳ tuần hoàn lớn đạt tới $2^{31}-1$ và mật độ phân bố đồng đều trên khoảng đơn vị $(0, 1)$.

Đối với các đại lượng vật lý có quy luật phân bố xác suất xác định, nghiên cứu áp dụng phối hợp phương pháp biến đổi hàm ngược và phương pháp chấp nhận - loại bỏ (Acceptance-Rejection). Quy trình chọn mẫu mô phỏng được thực hiện trên cỡ mẫu lớn gồm 100.000 đến 1.000.000 lịch sử hạt photon ($N = 10^5 - 10^6$) nhằm giảm thiểu sai số thống kê theo quy luật tỉ lệ nghịch với căn bậc hai của cỡ mẫu $1/\sqrt{N}$. Lý do lựa chọn phương pháp mô phỏng số trên môi trường Matlab là khả năng xử lý ma trận mạnh mẽ, hỗ trợ trực quan hóa hình học không gian 3 chiều của khối detector hình trụ và tối ưu hóa thời gian tính toán truy vết hạt qua ba vùng tán xạ (vùng phía trước, mặt bên và sau lưng). Dữ liệu tính toán được kiểm chứng trực tiếp bằng cách đối chiếu với phổ thực nghiệm đo trên hệ phổ kế đa kênh (MCA) trong cùng điều kiện hình học phòng thí nghiệm.

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

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

Quá trình mô phỏng số và phân tích dữ liệu thực nghiệm đã mang lại những phát hiện quan trọng:

• Tái tạo chính xác vị trí các đỉnh năng lượng toàn phần (Photopeak): Chương trình mô phỏng xác định đỉnh hấp thụ năng lượng của nguồn Cs-137 tại vị trí 0,662 MeV và hai đỉnh đặc trưng của Co-60 tại 1,173 MeV và 1,332 MeV với độ lệch vị trí kênh năng lượng dưới 1,2% so với thực nghiệm.
• Mô phỏng chân thực bờ tán xạ Compton (Compton edge): Bờ Compton xuất hiện rõ nét tại mức năng lượng 0,478 MeV đối với nguồn Cs-137, hoàn toàn trùng khớp với giá trị tính toán lý thuyết tương ứng với góc tán xạ cực đại 180 độ.
• Nhận diện đỉnh tán xạ ngược và tương tác thứ cấp: Mô hình tái hiện thành công đỉnh tán xạ ngược (Backscatter peak) tại vùng năng lượng khoảng 0,184 MeV sinh ra do photon gamma tán xạ từ vật liệu che chắn xung quanh quay trở lại detector.
• Độ phân giải năng lượng và tỷ số diện tích đỉnh: Tỷ số diện tích đỉnh trên tổng phổ (Peak-to-Total ratio) của phổ mô phỏng đạt độ tương đồng trên 93% khi so sánh với phổ thực nghiệm. Độ phân giải năng lượng (FWHM) tại đỉnh 0,662 MeV của nguồn Cs-137 đạt giá trị tối ưu xấp xỉ 7,5%.

Thảo luận kết quả

Cơ chế hình thành phổ biên độ xung trong detector nhấp nháy là kết quả của chuỗi tương tác vật lý liên hoàn diễn ra trong thời gian cực ngắn dưới $10^{-7}\text{ s}$. Khi tia gamma đi vào khối tinh thể NaI(Tl), các tương tác quang điện, Compton và tạo cặp làm giải phóng electron sơ cấp mang động năng cao. Electron này kích thích các ion trong mạng tinh thể, tạo ra các chớp sáng huỳnh quang được ống nhân quang điện chuyển đổi thành xung điện thế có biên độ tỷ lệ thuận với năng lượng hấp thụ.

Khi biểu diễn dữ liệu qua biểu đồ phổ năng lượng (trục hoành biểu thị kênh năng lượng từ 0 đến 2000 keV, trục tung biểu thị số đếm xung), đường cong mô phỏng thể hiện sự ăn khớp chặt chẽ với đường cong thực nghiệm ở cả ba thành phần: đỉnh hấp thụ toàn phần, vùng liên tục Compton và đỉnh tán xạ ngược. Bảng so sánh thông số kỹ thuật cho thấy sai lệch số đếm ở vùng năng lượng cao luôn duy trì dưới 3,5%. Một số sai lệch nhỏ xuất hiện ở vùng năng lượng thấp dưới 150 keV được giải thích là do hiện tượng tán xạ tạp tán từ vỏ nhôm bảo vệ detector và hiện tượng thoát tia X đặc trưng của Iodine (năng lượng 29 keV) chưa được mô hình hóa tuyệt đối trong thuật toán hình học đơn giản.

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

Dựa trên kết quả nghiên cứu đạt được, tác giả đưa ra 4 khuyến nghị và giải pháp hành động cụ thể:

1. Nâng cấp thuật toán mô phỏng trên nền tảng tính toán song song: Nhóm nghiên cứu cần chuyển đổi mã nguồn từ Matlab sang ngôn ngữ C/C++ kết hợp công nghệ xử lý đồ họa GPU CUDA. Mục tiêu là tăng tốc độ xử lý lên 20 đến 50 lần, rút ngắn thời gian mô phỏng $10^7$ lịch sử hạt từ vài giờ xuống dưới 3 phút. Thời gian hoàn thành dự kiến trong vòng 6 tháng.
2. Mở rộng mô hình hóa cho các hệ detector bán dẫn độ phân giải cao: Phòng thí nghiệm vật lý hạt nhân nên phát triển mở rộng chương trình cho các loại đầu dò hiện đại như HPGe (Germanium siêu tinh khiết) với độ phân giải năng lượng dưới 0,2% và tinh thể nhấp nháy BGO, LaBr3 phục vụ dải đo bức xạ lên đến 10 MeV. Lộ trình triển khai trong thời gian 9 đến 12 tháng.
3. Ứng dụng mô phỏng vào chương trình đào tạo đại học và sau đại học: Các khoa Vật lý tại các trường đại học cần tích hợp phần mềm mô phỏng này vào học phần thực hành Vật lý hạt nhân và Phương pháp số. Giải pháp này giúp các cơ sở giáo dục tiết kiệm 80% kinh phí trang bị nguồn phóng xạ kín và loại bỏ hoàn toàn nguy cơ mất an toàn bức xạ cho sinh viên. Thời gian áp dụng trong năm học tiếp theo.
4. Tối ưu hóa thiết kế che chắn an toàn bức xạ y tế và công nghiệp: Các cơ sở xạ trị và viện nghiên cứu hạt nhân nên sử dụng thuật toán Monte-Carlo để tính toán chiều dày che chắn tối ưu của các vật liệu compozit, bê tông nặng và hợp kim chì. Mục tiêu là duy trì liều hấp thụ của nhân viên vận hành dưới mức quy chuẩn an toàn 20 mSv/năm theo khuyến cáo quốc tế. Kế hoạch thực hiện trong vòng 12 đến 18 tháng do bộ phận kiểm soát an toàn bức xạ chủ trì.

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

Luận văn là tài liệu học thuật và kỹ thuật giá trị cao dành cho 4 nhóm đối tượng trọng tâm:

• 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ý tính toán: Nắm vững phương pháp luận xây dựng thuật toán Monte-Carlo, kỹ thuật phát sinh biến ngẫu nhiên theo phân bố xác suất và phương pháp truy vết hạt vi mô trong không gian ba chiều.
• Kỹ sư Vật lý y khoa và Chuyên viên kỹ thuật xạ trị: Hiểu rõ bản chất tương tác vi mô của chùm tia gamma với detector đo liều và mô sinh học, hỗ trợ công tác hiệu chuẩn liều lượng và đánh giá độ chính xác của các hệ đo phổ kế trong bệnh viện ung bướu.
• Giảng viên các trường đại học khối Khoa học Tự sinh và Sư phạm: Sử dụng làm giáo trình tham khảo chuyên sâu, xây dựng các bài thí nghiệm ảo trực quan nhằm minh họa sinh động các hiệu ứng tương tác bức xạ trừu tượng cho sinh viên.
• Chuyên viên kiểm định an toàn bức xạ và đánh giá môi trường: Khai thác mô hình tính toán hệ số suy giảm tuyến tính và hiệu suất ghi nhận để thẩm định các công trình che chắn lò phản ứng, thiết bị chiếu xạ công nghiệp và kho lưu giữ nguồn phóng xạ.

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

1. Phương pháp Monte-Carlo có ưu thế gì vượt trội so với phương pháp giải tích trong bài toán phổ kế bức xạ?
Phương pháp giải tích chỉ có thể giải quyết các bài toán một chiều với điều kiện biên hình học lý tưởng và bỏ qua quá trình tán xạ phức tạp. Trong thực tế, chùm photon gamma trải qua chuỗi tương tác tán xạ Compton nhiều lần và thay đổi hướng liên tục trong không gian tinh thể ba chiều. Monte-Carlo cho phép theo dõi chi tiết từng lịch sử hạt với độ chính xác cao và sai số thống kê dưới 1%, phản ánh trung thực toàn bộ cơ chế vật lý vi mô.

2. Tại sao tinh thể NaI(Tl) lại là lựa chọn hàng đầu để chế tạo detector ghi nhận bức xạ gamma?
Tinh thể NaI(Tl) sở hữu khối lượng riêng cao $3,67\text{ g/cm}^3$ cùng số nguyên tử hiệu dụng trung bình bằng 50, mang lại hiệu suất hấp thụ photon gamma năng lượng từ 0,1 MeV đến 3 MeV lên đến 80% - 90%. Ngoài ra, thời gian phát sáng cực ngắn $0,25,\mu\text{s}$ giúp giảm tối đa thời gian chết của hệ đo, cho phép detector ghi nhận chính xác các nguồn có hoạt độ phóng xạ cao mà không bị méo phổ.

3. Hiệu ứng tạo cặp diễn ra trong điều kiện nào và đóng góp thế nào vào phổ năng lượng?
Hiệu ứng tạo cặp chỉ xuất hiện khi năng lượng photon gamma vượt ngưỡng 1,022 MeV, tương ứng năng lượng nghỉ của một cặp electron và positron. Khi tương tác với điện trường hạt nhân, photon chuyển hóa hoàn toàn thành cặp hạt này; sau đó positron hủy hạt với electron môi trường giải phóng 2 photon 0,511 MeV bay ngược hướng nhau, tạo nên các đỉnh thoát đơn (Single escape) và thoát kép (Double escape) trên phổ ghi nhận.

4. Thuật toán chấp nhận - loại bỏ hoạt động như thế nào trong việc tạo biến ngẫu nhiên năng lượng?
Thuật toán bao bọc hàm mật độ xác suất phức tạp bằng một hàm phân bố đơn giản có giá trị cực đại xác định. Chương trình phát sinh đồng thời một giá trị năng lượng ngẫu nhiên và một số ngẫu nhiên kiểm tra; nếu tỷ số xác suất thỏa mãn điều kiện ngưỡng thì giá trị được chấp nhận, ngược lại sẽ bị loại bỏ để phát sinh mẫu mới. Phương pháp này đạt hiệu suất chọn lọc trên 78% đối với các phân bố góc Klein-Nishina.

5. Yếu tố nào gây ra sự sai khác nhỏ giữa phổ mô phỏng và phổ thực nghiệm ở vùng năng lượng thấp?
Sự chênh lệch bắt nguồn từ hiện tượng tán xạ photon gamma trên các cấu trúc bao quanh phòng thí nghiệm như tường chắn chì, giá đỡ detector và không khí chưa được đưa đầy đủ vào mô hình hình học. Bên cạnh đó, hiệu ứng trượt biên và sự phát sinh tia X đặc trưng 29 keV của Iodine thoát ra khỏi bề mặt tinh thể cũng làm biến dạng nhẹ số đếm tại vùng năng lượng dưới 100 keV.

Kết luận

• Xây dựng thành công thuật toán mô phỏng Monte-Carlo tái hiện trọn vẹn quá trình truyền bức xạ gamma qua môi trường vật chất.
• Tái lập chính xác phổ năng lượng của các nguồn phóng xạ chuẩn Cs-137 và Co-60 trên đầu dò NaI(Tl) với sai số vị trí đỉnh dưới 1,2%.
• Làm sáng tỏ định lượng đóng góp của hiệu ứng quang điện, tán xạ Compton đa cấp và hiệu ứng tạo cặp đối với từng vùng cấu trúc phổ.
• Cung cấp giải pháp tính toán mô phỏng số giúp tiết kiệm hơn 70% kinh phí thực nghiệm và đảm bảo an toàn bức xạ tuyệt đối.
• Thiết lập nền tảng phần mềm khoa học hỗ trợ đắc lực cho công tác nghiên cứu, đào tạo đại học và thiết kế che chắn phóng xạ.

Đóng góp quan trọng nhất của luận văn là chứng minh tính khả thi và độ tin cậy vượt trội của phương pháp mô phỏng Monte-Carlo trong việc giải quyết các bài toán vật lý hạt nhân phức tạp tại Việt Nam. Về kế hoạch tiếp theo, nhóm tác giả dự kiến mở rộng mô hình sang hệ đa đầu dò bán dẫn HPGe và tích hợp giao diện tương tác đồ họa trong vòng 6 đến 12 tháng tới. Quý độc giả, các nhà nghiên cứu và học viên quan tâm có thể khai thác mã nguồn và tài liệu toàn văn của luận văn để ứng dụng hiệu quả vào các dự án nghiên cứu vật lý bức xạ và kỹ thuật hạt nhân thực tế.

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

MỞ ĐẦU HOR Mô phỏng (Simulation) là kỹ thuật sử dung máy tinh để mô tả lại các quá trình. các tính chất khác nhau của hệ máy hay là các tiến trình vận động. Các hệ máy hay tiến trình mà chúng ta quan tâm thường được gọi là hệ thống (system). Và để tìm hiểu về tính khoa học của một hệ thống, ta thường phải tập hợp tất cả các giả thuyết về sự vận động của nó, từ đó rút ra các dạng toán học hoặc các mối quan hệ logic để thiết lập nên một mô hình (Model), Từ mô hình này người ta sẽ nghiên cứu xem làm thế nào để một hệ thống có thể hoạt động tối ưu.

Nếu những mối quan hệ để tạo ra các mô hình là độc lập với nhau, có thể sử dụng các phương pháp toán học như đại số, giải tích hay lý thuyết xác xuất để tìm hiểu về các vấn để quan tam, Phương pháp này gọi là phương pháp phân tích (analytic solution). Tuy nhiên, hấu hết các hệ trên thực tế déu rất phức tạp để có thể tạo dựng một mô hình cho phép nghiên cứu bằng phương pháp phân tích. Và khi đó, những mô hình như thế phải được nghiên cứu bằng phương pháp mô phỏng (Simulation). Trong thời đại ngày nay khi khoa học phát triển như vũ bão thì vật lý luôn giữ vai trò quan trọng của minh, Nghành vật lý hạt nhân bên cạnh phương pháp nghiên cứu truyền thống đã bất đầu nảy sinh những hướng đi mới.

Và phương pháp mô phỏng đã phát triển nổi bật với những ưu thế của mình, giảm đi những phí tổn thí nghiệm và những han chế trong thực nghiệm. Đặc biệt là khi máy tính ngày càng phát triển mạnh như hiện nay thì sức mạnh của mô phỏng càng được nang lên. Tuy nhiên mô phỏng cũng có hạn chế cơ bản của nó, vì vậy mô phỏng và thực nghiệm không thể tách rời nhau. Hiện nay, trên thế giới phương pháp mô phỏng Monte-Carlo được áp dụng rong rãi trong moi lĩnh vực trong đó có lĩnh vực vật lý hạt nhân.

Trong khuôn khổ luận văn này, trình bày phương pháp mô phỏng Monte-Carlo được sử dụng để xây dựng chương trình mô phỏng phổ bức xạ năng lượng của một số nguồn phóng xạ. Với để tài này đã giúp em tìm hiểu nhiều về quá trình bức xạ truyền qua môi trường vật chất. Việc hiểu các quá trình tương tắc này là vô cùng quan trong đối với bất kì ai muốn làm việc với | dạng năng lượng bất kỳ hoặc sử dụng dạng năng lượng này. Việc hiểu đó là cần thiết để có thể giải thích giá trị cơ bản trong việc thiết kế các thiết bị dùng để biến đổi thành nang lượng dễ sử dung hơn như nhiệt năng và động năng.

Ta không thể trả lời các câu hỏi bức xạ hạt nhân nguy hiểm đối với sức khỏe của những người dùng nó như thế nào và làm cách nào để bảo vệ an toàn lao động cho những người làm việc SVTH: NGUYÊN THỊ BÍCH LOAN. -1- | LOI MỞ DAU với các bức xa ấy nếu không tìm hiểu rõ các quá trình tương tác của bức xa qua môi trường vật chất. Trong khảo sát về tương tác của bức xa, đặc biệt là bức xạ năng lượng cao như tia Gamma với vật chất, mô phỏng sẽ giảm được việc tiếp xúc trực tiếp với chất phóng xạ. Đồng thời mô phỏng bức xạ sẽ giúp con người giảm đi chi phí thí nghiệm.

Công việc mô phỏng còn giúp cũng cố lý thuyết về tương tác trong vật chất. Đồng thời, chính dé tài này đã giúp em tìm hiểu và tiếp cận phương pháp mô phỏng hiện đại là phương pháp Monte — Carlo, cũng như kỹ thuật vi tính đó là ngôn ngữ lập trình Matlab. Luận văn được trình bày bao gồm 2 phần: Phần I: Cơ sở lý thuyết. - Phương pháp Monte — Carlo.

- Tương tác giữa tia Gammazw vat chal - Cơ chế ghi nhận bức xạ. Phần 2: Xây dựng chương trình mô phỏng. Từ cơ sở lý thuyết ở phẩn 1 ta tiến hành lập chương trình mô phỏng phổ của nguén bức xạ Gamma bằng ngôn ngữ toán học và sau đó chuyển sang ngôn ngữ Matlab. Các thông số đầu vào của chương trình mô phỏng như detector ghi nhận bức xạ, điểu kiện bố trí hệ do.

được lấy theo điểu kiện của hệ đo từ phòng thí nghiệm bộ môn Vật lý hạt nhân — Trường Đại hoc Khoa học Tự nhiên Thanh phố Hổ Chí Minh. Phổ mô phỏng đối chiếu với phổ thực nghiệm. Dựa vào độ phù hợp của 2 phổ ta sẽ đánh giá được mức độ thành công của chương trình. SVTH: NGUYEN THỊ BICH LOAN.

-2- CƠSỞPLHÝANT:UYE;T Phương pháp mô phỏng Monte - Carlo dựa trên ý tưởng các tập hợp thống kế. Xét một hệ vi mô ở trạng thái cân bằng mà ta muốn xác định các tính chất, ta tạo ru | số rất lớn những hệ tương tự. Sau khi tính toán giá trị của các đại lượng mà ta muốn có, ta sẽ lấy giá trị trung bình trên tập hợp. Như vậy, sau khi xác định một cấu hình ban đầu của hệ, ta sẽ làm phát sinh một dãy những cấu hình chọn lựa một cách ngẫu nhiên, tức là phải có một phần mềm cho ta các số ngdu nhiên trên máy tính.

Phương pháp Monte - Carlo có thể xem như phương pháp mô hình hoá các đại lượng ngẫu nhiên nhằm tính các đại lượng đặc trưng của chúng. Tên gọi “M6 phỏng Monte - Carlo” xuất hiện từ thế chiến thứ hai, khí mà phương pháp này áp dụng cho các vấn để có liên quan đến kế hoạch phát triển bom nguyên tử. Mới đầu phương pháp này chỉ chủ yếu sử dụng trong việc giải các bài toán vật lý hạt nhân, sau đó ảnh hưởng sang các bài toán vật lý thống kê với nội dung khác nhau. Ngoài ra phương pháp Monte — Carlo còn được dùng rất nhiều trong các ngành khoa học lý thuyết phục vụ phẩn lớn lý thuyết trò chơi, toán kinh tế, lý thuyết truyền thông.

Phương pháp mô phỏng đã mang lại những kết quả rất quan trọng trong lĩnh vực plasma một thành phần quan trọng trong những năm gần đây. Cũng nên biết rằng gần đây, các phép mô phỏng Monte - Carlo lượng tử cũng đã được thử trên các máy tính tốc độ cao. SỐ NGẪU NHIEN- GIẢ NGAU NHIÊN TREN MAY TÍNH : Phương pháp Monte-Carto có thành công hay không còn tuỳ thuộc rất nhiều vào tập hợp các số ngẫu nhiên có đáng tin cậy hay không. Phương pháp đã tạo ra số ngẫu nhiên có lịch sử lâu đời và hấp dẫn.

Những phương pháp sớm nhất cơ bản được thực hiện bằng tay như là trò chơi rút thăm, xúc xắc, trò chơi xổ số. Đầu thế kỷ 20, các nhà thống kê quan tâm đến số ngẫu nhiên và những dung cụ tạo ra số ngẫu nhiên sớm được tạo ra. Người ta nghĩ ra lập bảng các số ngẫu nhiên chứa các chữ số thập phân hợp thành từng nhóm. Khi lập bảng rất khó loại trừ ảnh hưởng của các yếu tế phá vỡ tính ngẫu nhiên và tính đều đặn của phân bố.

Để lập bảng ngẫu nhiên gồm 100.000 con số, cuối năm 1930, Kendall và Babengton - Smith dùng một cái đĩa quay chia thành 10 hình quạt bằng nhau ứng với các số 0,1,.8,9 hãng Rand lập ra bảng ngẫu nhiên gồm một triệu chữ số ngẫu nhiên. Người ta dùng một thiết bị tác động đặc biệt gọi là bánh xe điện đó là thiết bị san sinh và truyền các số ngẫu nhiên hoạt động dựa trên nguyên tắc vật lý (tạo ra các xung ngẫu nhiên). Các số ngẫu nhiên hiện ra với tốc độ 50 km/ giây. Vào những năm 1940, 1950 việc tạo ra số ngẫu nhiên bằng phương pháp số học.

Mỗi số mới trong dãy số ngẫu nhiên được xác định thông qua một hoặc một vài bước làm việc của máy tính trước đó bởi | biểu thức toán học nhất định. Số nhận được theo cách này gọi là các số giả ngẫu nhiên. Ví dụ: Phương pháp Midsquare : Bình phương số giữa. Bắt đầu bằng số nguyên dương Z¿ có 4 chữ số.

Bình phương nó ta được số có 8 chữ số, nếu cẩn thiết thêm số 0 vào trước số vừa tạo để chấc chắn nó có 8 chữ số. Lấy 4 chữ số chính giữa của Z,7 , ta được số đó là số 4 chữ số tiếp theo Z,. Dat dấu phẩy bên trái của Z, để đạt được số ngẫu nhiên nằm trong khoảng U (0. Tương tự Z; là số có 4 chữ số giữa của Z;ˆ và U; là số thập phân tương ứng của Z và cứ tiếp với Zp =7182 ta có bang: khít mm_ Z4 |} Uy —| 0|| 7182 | _| 1| 581 |0581| 33767721 1677 58936329 3 9363 0.9363 87665749 | 4 | 6657 0,6657 44315629 | 03156 —_ 336 Ta có cảm giác rằng phương pháp bình phương số giữa đã tao ra ngẫu nhiên.

Thực tế, việc làm đó không đơn giản chút nào. Một vấn để quan trọng là nó có khuynh hướng thoái hóa nhanh về số 0 và ở đây nó giữ như thế mãi mãi. =ls Hls HE 2Š el I|šesElS|Eš : Sz == ie =a lll. PHƯƠNG PHÁP TẠO ĐỒNG DƯ TUYẾN TÍNH: (Linear congruential generators: phương phấp so sánh) Nhiều cách phát sinh số ngẫu nhiên được dùng ngày nay dựa trên những phương pháp tạo đồng dư tuyến tính được Lechmer đưa ra 1951.

Phương pháp này tim ra dãy số ngẫu nhiên gồm những số nguyên Z, Z;,. được xác định bởi công thức nhị phân: Z,=(aZ,,+C} (mod m) (*) Trong d6: a: hệ số nhân. Z4: giá trị ban đầu. Với a, m, C, Z4 là những số nguyên không âm.

Để có Z,, ta chia aZ,., cho m và Z, là phần dư của phép chia Dođó 0<Z, <m-l Như vậy muốn có được số gia ngẫu nhiên trong đoạn [ 0. 1] ta đặt: Z m để có diéu kiện các U, không âm thì các số nguyên m, a, C, Zp phải thoả mãn m>0, a<m, C<m, Z¿<m Trong công thức (*) , nếu C = 0 công thức có dạng đơn giản nhất. Khi đó phương pháp tạo số ngẫu nhiên có tên gọi là phương pháp nhãn tử, với C > 0 ta gọi là phương pháp hổn hợp. Phương pháp hồn hợp : ( Mixed Generators) Cho C >0, với điều kiện a thích hợp, thì chúng ta có thể dat chu kỳ đẩy m.

Để có chu kỳ lớn và mật độ phân bố cao của của U, trong khoảng [0,1] chúng ta cần tạo m lớn. Việc chọn m tốt nhất là m = 2. Với b là số bit nhị phân, là ngôn ngữ trong tin học, biểu thị các ô nhớ b bít chứa các dữ liệu. Số nguyên lớn nhất được biểu diễn là 2°~ 1.