Khóa luận tốt nghiệp: Nghiên cứu khả năng ứng dụng của mô phỏng quang phổ alpha AASI

Khóa luận nghiên cứu ứng dụng của mô phỏng quang phổ alpha AASI trong vật lý, mở ra hướng đi mới cho nghiên cứu và phát triển công nghệ.

• Tác giả: Lâm Thu Văn
• Người hướng dẫn: ThS. Lê Công Hảo
• 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ý
• Đề tài: Khả Năng Ứng Dụng Của AASI Advanced Alpha Spectrometric Simulation
• Loại tài liệu: khóa luận tốt nghiệp
• Năm xuất bản: 2010
• Địa điểm: Thành Phố Hồ Chí Minh

Mục lục chi tiết

• LỜI CẢM ƠN
• LỜI NÓI ĐẦU
• 1. CHƯƠNG 1: TỔNG QUAN VỀ HẠT ALPHA
• 1.1. Khái niệm hạt alpha
• 1.2. Sự phân rã alpha
• 1.3. Các đặc trưng chính của phân rã alpha
• 1.4. Cơ chế phân rã alpha
• 1.5. Sự tương tác của hạt alpha với vật chất
• 1.5.1. Hiệu ứng ion hóa
• 1.5.2. Sự mất năng lượng của hạt alpha do quá trình ion hóa
• 2. CHƯƠNG 2: GIỚI THIỆU VỀ HỆ ĐO ALPHA ANALYST
• 3. CHƯƠNG 3: GIỚI THIỆU VỀ PHẦN MỀM MÔ PHỎNG AASI
• 4. CHƯƠNG 4: CÁC KẾT QUẢ THỰC NGHIỆM MÔ PHỎNG VÀ SỰ SO SÁNH

Tóm tắt

I. Tổng quan về khóa luận tốt nghiệp vật lý ứng dụng AASI

Khóa luận tốt nghiệp vật lý nghiên cứu khả năng ứng dụng của AASI Advanced Alpha Spectrometric Simulation là một chủ đề quan trọng trong lĩnh vực vật lý hạt nhân. Nghiên cứu này không chỉ giúp sinh viên nắm vững kiến thức lý thuyết mà còn phát triển kỹ năng thực hành trong việc mô phỏng và phân tích dữ liệu. Khóa luận tốt nghiệp vật lý này sẽ trình bày chi tiết về các phương pháp mô phỏng và ứng dụng của phần mềm AASI trong nghiên cứu hạt alpha.

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

AASI là phần mềm mô phỏng tiên tiến, được thiết kế để nghiên cứu và phân tích các hạt alpha. Phần mềm này cho phép người dùng mô phỏng các quá trình phức tạp trong vật lý hạt nhân, từ đó giúp nâng cao độ chính xác trong việc đo đạc và phân tích dữ liệu.

1.2. Tầm quan trọng của nghiên cứu hạt alpha

Hạt alpha là một trong những dạng phóng xạ quan trọng trong vật lý hạt nhân. Nghiên cứu về hạt alpha không chỉ giúp hiểu rõ hơn về cấu trúc hạt nhân mà còn có ứng dụng trong nhiều lĩnh vực như y học, năng lượng hạt nhân và bảo vệ môi trường.

II. Vấn đề và thách thức trong nghiên cứu hạt alpha

Nghiên cứu hạt alpha gặp phải nhiều thách thức, đặc biệt là trong việc mô phỏng và phân tích dữ liệu. Các vấn đề như hiện tượng trùng phùng và sự suy giảm năng lượng của hạt alpha trong quá trình truyền qua vật chất là những yếu tố cần được xem xét kỹ lưỡng. Việc hiểu rõ các vấn đề này sẽ giúp cải thiện độ chính xác của các phép đo và mô phỏng.

2.1. Hiện tượng trùng phùng trong đo đạc

Hiện tượng trùng phùng xảy ra khi hai hạt alpha phát ra gần như đồng thời, gây khó khăn trong việc phân tích dữ liệu. Điều này đòi hỏi các nhà nghiên cứu phải phát triển các phương pháp mới để tách biệt và phân tích các tín hiệu này.

2.2. Suy giảm năng lượng của hạt alpha

Khi hạt alpha truyền qua vật chất, nó sẽ mất năng lượng do ion hóa và va chạm với các nguyên tử. Sự suy giảm năng lượng này ảnh hưởng đến độ chính xác của các phép đo và cần được tính toán cẩn thận trong các mô phỏng.

III. Phương pháp mô phỏng alpha spectrometry hiệu quả

Phương pháp mô phỏng alpha spectrometry sử dụng phần mềm AASI đã được chứng minh là hiệu quả trong việc nghiên cứu và phân tích dữ liệu hạt alpha. Các phương pháp này không chỉ giúp cải thiện độ chính xác mà còn tiết kiệm thời gian và chi phí trong nghiên cứu.

3.1. Sử dụng phần mềm AASI trong mô phỏng

Phần mềm AASI cho phép mô phỏng các quá trình phức tạp liên quan đến hạt alpha, từ đó giúp người dùng có cái nhìn sâu sắc hơn về các hiện tượng vật lý. Việc sử dụng AASI giúp tối ưu hóa quy trình nghiên cứu và nâng cao độ chính xác của kết quả.

3.2. So sánh kết quả mô phỏng và thực nghiệm

Việc so sánh kết quả từ mô phỏng với kết quả thực nghiệm là một bước quan trọng trong nghiên cứu. Điều này không chỉ giúp xác nhận độ chính xác của mô phỏng mà còn cung cấp thông tin quý giá để cải thiện các phương pháp nghiên cứu.

IV. Ứng dụng thực tiễn của AASI trong nghiên cứu vật lý

AASI không chỉ được sử dụng trong nghiên cứu hạt alpha mà còn có nhiều ứng dụng thực tiễn khác trong lĩnh vực vật lý. Các ứng dụng này bao gồm việc phát triển các thiết bị đo đạc chính xác hơn và cải thiện quy trình phân tích dữ liệu trong các nghiên cứu khoa học.

4.1. Ứng dụng trong y học

Trong y học, AASI có thể được sử dụng để nghiên cứu các phương pháp điều trị phóng xạ, giúp cải thiện hiệu quả điều trị cho bệnh nhân. Việc mô phỏng hạt alpha trong môi trường sinh học là một lĩnh vực nghiên cứu đầy tiềm năng.

4.2. Ứng dụng trong bảo vệ môi trường

AASI cũng có thể được áp dụng trong việc theo dõi và phân tích ô nhiễm phóng xạ trong môi trường. Điều này giúp các nhà khoa học đánh giá được mức độ ô nhiễm và đưa ra các biện pháp khắc phục kịp thời.

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

Nghiên cứu khả năng ứng dụng của AASI trong lĩnh vực vật lý hạt nhân mở ra nhiều cơ hội mới cho các nhà nghiên cứu. 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 việc hiểu biết về hạt alpha và các ứng dụng của nó trong thực tiễn.

5.1. Triển vọng phát triển phần mềm AASI

Phần mềm AASI có tiềm năng phát triển mạnh mẽ trong tương lai, với nhiều tính năng mới được cập nhật để đáp ứng nhu cầu nghiên cứu ngày càng cao. Việc cải tiến phần mềm sẽ giúp nâng cao độ chính xác và hiệu quả trong nghiên cứu.

5.2. 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 mở rộng khả năng ứng dụng của AASI trong các lĩnh vực khác nhau, từ y học đến công nghiệp. Điều này sẽ giúp nâng cao giá trị của phần mềm trong cộng đồng khoa học.

Nội dung chính

Bối cảnh và vấn đề nghiên cứu

Trong lĩnh vực vật lý hạt nhân, phương pháp mô phỏng Monte Carlo đóng vai trò quan trọng trong việc mô hình hóa các thiết bị đo đạc và tái hiện các quá trình tương tác bức xạ phức tạp như tán xạ Rutherford, tán xạ Compton, hấp thụ quang điện và sự mất năng lượng liên tục của hạt mang điện. Phép đo phổ kế alpha là kỹ thuật phân tích được sử dụng phổ biến nhằm xác định định tính và định lượng hoạt độ phóng xạ trong các mẫu môi trường. Tuy nhiên, do hạt alpha có khả năng ion hóa mạnh và quãng chạy ngắn, quá trình đo đạc thực nghiệm gặp nhiều khó khăn: các mẫu có hoạt độ thấp đòi hỏi thời gian đo kéo dài và khoảng cách giữa nguồn và đầu dò (Source to Detector Distance - SDD) phải thu hẹp. Khoảng cách SDD nhỏ dẫn đến hiện tượng trùng phùng (sự phát hạt alpha xảy ra gần như đồng thời với bức xạ tiếp theo từ hạt nhân con cháu) và làm gia tăng hiệu ứng chồng chập phổ do sự suy giảm phẩm chất phổ khi hạt mất năng lượng qua vật chất.

Trước thời điểm nghiên cứu (năm 2010), các bộ chương trình mô phỏng Monte Carlo kinh điển như TRIM, GEANT hay MCNP đã được ứng dụng để mô phỏng tương tác của hạt alpha, nhưng việc cấu hình thường phức tạp. Bộ mã AASI (Advanced Alpha-spectrometric Simulation) được phát triển trên ngôn ngữ Fortran 95, tích hợp cơ sở dữ liệu phân rã dạng XML nhằm mô phỏng toàn diện phổ năng lượng alpha cho nhiều dạng hình học nguồn khác nhau. Dù có thời gian tính toán nhanh (từ vài giây đến hai phút trên cấu hình máy tính Pentium 1.6 GHz), phần mềm AASI khi đó vẫn chưa được phổ biến tại Việt Nam. Khóa luận tốt nghiệp của tác giả Lâm Thu Văn, thực hiện tại Bộ môn Vật lý Hạt nhân, Khoa Vật lý, Trường Đại học Sư phạm Thành phố Hồ Chí Minh (dưới sự hướng dẫn của ThS. Lê Công Hảo - Trường Đại học Khoa học Tự nhiên, Đại học Quốc gia TP. Hồ Chí Minh), tập trung giải quyết vấn đề ứng dụng công cụ tính toán này vào điều kiện thực nghiệm trong nước.

Nhiệm vụ nghiên cứu của đề tài được xác định cụ thể gồm hai nội dung chính:

1. Nghiên cứu cấu trúc giao diện và hướng dẫn chi tiết cách thức thiết lập các thông số vận hành của phần mềm mô phỏng phổ kế alpha AASI.
2. Thiết lập mô phỏng tương ứng với các điều kiện đo đạc thực tế trên hệ phổ kế alpha chuyên dụng, tiến hành so sánh đối chiếu giữa phổ năng lượng alpha mô phỏng và phổ thực nghiệm nhằm đánh giá mức độ tương thích của phần mềm.

Đối tượng nghiên cứu của khóa luận là phổ năng lượng của bức xạ alpha phát ra từ các nguồn phóng xạ chuẩn và mẫu Uranium. Phạm vi thực nghiệm được thực hiện trên hệ đo phổ kế Canberra Alpha Analyst đặt tại Bộ môn Vật lý Hạt nhân thuộc Trường Đại học Khoa học Tự nhiên TP. Hồ Chí Minh, kết hợp xử lý và so sánh số liệu phổ thông qua phần mềm Microsoft Excel.

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

Khóa luận xây dựng nội dung dựa trên nền tảng lý thuyết phân rã phóng xạ alpha và cơ chế tương tác bức xạ ion hóa với vật chất. Phân rã alpha là quá trình biến đổi hạt nhân mẹ $^{A}{Z}X$ thành hạt nhân con $^{A-4}{Z-2}Y$ và phát ra hạt nhân $^{4}_{2}\text{He}$. Quá trình này tuân theo điều kiện cân bằng khối lượng và bảo toàn động lượng, trong đó phần lớn động năng phân rã $\Delta E$ được truyền cho hạt alpha theo tỷ lệ khối lượng:

$$E_{\alpha} = \Delta E \cdot \frac{m_c}{m_c + m_{\alpha}}$$

Về mặt động học, sự phụ thuộc của chu kỳ bán rã $T_{1/2}$ vào động năng hạt alpha $E$ được mô tả qua định luật bán thực nghiệm Geiger-Nuttall:

$$\lg T_{1/2} = C + \frac{D}{\sqrt{E}}$$

Trong đó $C$ và $D$ là các hằng số phụ thuộc vào điện tích hạt nhân $Z$. Về cơ chế lượng tử, quá trình hạt alpha thoát ra khỏi hạt nhân được giải thích bằng hiệu ứng xuyên ngầm qua hàng rào thế Coulomb và thế ly tâm dựa trên phương trình Schrödinger, giải thích nghịch lý hạt có năng lượng từ 4 MeV đến 9 MeV có thể vượt qua bờ thế cao khoảng 30 MeV.

Khi truyền qua môi trường, hạt alpha mất năng lượng chủ yếu qua hai quá trình: ion hóa và kích thích nguyên tử. Năng lượng trung bình cần thiết để tạo một cặp ion ($w$) trong chất khí hoặc nước vào khoảng 30 eV, trong khi đối với chất bán dẫn chỉ khoảng vài electron-volt. Mức độ mất năng lượng riêng theo quãng đường (Stopping Power) được tính toán dựa trên lý thuyết va chạm cổ điển của Bohr kèm theo các hiệu chỉnh hiệu ứng mật độ và lớp vỏ liên kết. Quãng chạy của hạt alpha trong không khí được xác định gần đúng theo công thức $R_{kk} \approx 0,318 E^{3/2}$ (đối với $3 < R_{kk} < 7\text{ cm}$), và trong môi trường có khối lượng nguyên tử $A$ theo hệ thức $R = 0,56 R_{kk} A^{1/3}$.

Phương pháp nghiên cứu kết hợp giữa mô phỏng số và thực nghiệm kiểm chứng:

• Phương pháp mô phỏng Monte Carlo: Sử dụng bộ mã AASI phiên bản Fortran 95 để theo dõi lịch sử chuyển động của các hạt alpha từ nguồn qua các lớp vật liệu đến detector. Nguồn dữ liệu hạt nhân đầu vào được trích xuất từ các tệp cơ sở dữ liệu NuclideLibrary.xml (thông số phân rã) và CoincidenceLibrary.xml (thông số trùng phùng).
• Phương pháp đo phổ thực nghiệm: Sử dụng hệ đo Alpha Analyst tích hợp buồng chân không (áp suất đạt dưới 0,1 torr nhằm loại trừ sự mất năng lượng trong không khí) và đầu dò bán dẫn PIPS (Passivated Implanted Planar Silicon) model A1200-37AM. Tín hiệu xung từ đầu dò được xử lý qua bộ tiền khuếch đại nhạy điện tích Canberra 2004, bộ khuếch đại tạo xung chuẩn Gauss, bộ biến đổi tương tự - số ADC và hiển thị trên phần mềm Genie-2000.
• Phương pháp xử lý số liệu: Chuyển đổi dữ liệu kênh - số đếm từ phần mềm Genie-2000 và kết quả mô phỏng từ thẻ đồ thị của AASI sang bảng tính Microsoft Excel để vẽ đường cong phân bố năng lượng và so sánh.

Nội dung chính theo từng chương

Chương 1: Tổng quan về hạt alpha

Chương 1 hệ thống hóa các đặc trưng vật lý của hạt alpha và cơ chế phân rã hạt nhân. Văn bản chỉ ra rằng hiện có hơn 200 hạt nhân phân rã alpha, chủ yếu tập trung ở các nguyên tố nặng cuối bảng tuần hoàn ($Z > 83$) và một nhóm nhỏ trong vùng đất hiếm ($A = 140 - 160$). Điều kiện để xảy ra phân rã alpha là năng lượng liên kết riêng của hạt nhân mẹ phải nhỏ hơn 7 MeV (năng lượng liên kết riêng của hạt alpha là 7 MeV, tương ứng năng lượng liên kết toàn phần 28 MeV).

Văn bản trình bày chi tiết hằng số định luật Geiger-Nuttall cho hai trường hợp cụ thể: đối với $Z=84$ ($C=-50,15$; $D=128,8$) và $Z=90$ ($C=-51,94$; $D=139,4$). Tác giả phân tích cơ chế vượt rào thế bằng hiệu ứng đường ngầm lượng tử qua việc giải phương trình Schrödinger với thế Coulomb $U_c = 2Ze^2/r$ và thế ly tâm bổ sung $U_l = \hbar^2 l(l+1) / (2mr^2)$. Phần cuối chương tập trung vào cơ chế mất năng lượng qua va chạm không đàn hồi với electron, định nghĩa về electron delta ($\delta$-electron), sự biến thiên của độ ion hóa riêng dọc theo quãng đường tạo nên đường cong Bragg, và phân biệt giữa khái niệm quãng chạy trung bình với quãng chạy ngoại suy.

Chương 2: Giới thiệu hệ đo Alpha Analyst

Chương 2 mô tả chi tiết kiến trúc phần cứng và chức năng của từng khối trong hệ phổ kế Alpha Analyst của hãng Canberra. Hệ thống bao gồm:

• Buồng chân không: Áp suất vận hành đạt dưới 0,1 torr, trang bị hệ thống rãnh kẹp mẫu đa tầng cách đều nhau giúp thay đổi linh hoạt khoảng cách SDD, tích hợp màng chắn ngăn chặn các hạt bị giật lùi để giảm thiểu nhiễu.
• Detector bán dẫn PIPS: Mẫu A1200-37AM có diện tích bề mặt hoạt động $1200\text{ mm}^2$ (bán kính $19,55\text{ mm}$), độ phân giải năng lượng alpha là $37\text{ keV}$, điện thế phân cực khuyến nghị $+40\text{ V}$ (giới hạn tối đa khoảng 100 V), dòng rò ở $20^\circ\text{C}$ là $12\text{ nA}$, độ sâu vùng nghèo tối thiểu $> 140\ \mu\text{m}$ (đủ để hấp thụ hoàn toàn hạt alpha năng lượng tới 15 MeV), và tốc độ đếm phông nền cực thấp ($< 0,05\text{ counts/hour/cm}^2$).
• Bộ tiền khuếch đại model 2004: Là loại nhạy điện tích đặt sát đầu dò để giảm tạp âm; độ trôi hệ số khuếch đại $< \pm 0,02%$ cho đầu ra 10 V; độ nhạy năng lượng $9\text{ mV/MeV}$ hoặc $45\text{ mV/MeV}$.
• Bộ khuếch đại định hình xung: Sử dụng mạch lọc tạo xung chuẩn Gauss, mang lại tỷ số tín hiệu trên tạp âm (S/N) tốt hơn từ $17%$ đến $19%$ so với mạch lọc CR-RC thông thường, đồng thời giảm độ rộng xung từ $22%$ đến $52%$ tại mức $0,1%$ biên độ.
• Hệ phân tích đa kênh (MCA/ADC): Chuyển đổi biên độ xung thành kênh số liệu và kết nối điều khiển hoàn toàn bằng máy tính qua phần mềm Genie-2000 Alpha Acquisition & Analysis.

Chương 3: Giới thiệu về phần mềm AASI

Chương 3 trình bày cấu trúc vận hành của phần mềm AASI qua 7 thẻ giao diện chức năng chính:

• Welcome: Cung cấp giao diện khởi động và 4 tệp ví dụ mẫu (Example 1 đến Example 4) có sẵn để làm quen.
• Alpha Detector: Cấu hình các đặc tính đầu dò gồm đường kính cửa sổ, độ phân giải FWHM, bề dày vùng hoạt động, bề dày lớp chết (Dead layer thickness) bằng vật liệu Silic, cùng 4 tham số hình dạng đỉnh phổ (Peak shape): Đuôi toàn phần (Total tailing từ 0 đến 1), thành phần nhanh (Fast component), thành phần chậm (Slow component), và tỷ số diện tích nhanh/chậm (Area ratio fast/slow).
• Source: Thiết lập khoảng cách SDD và lựa chọn 1 trong 3 dạng cấu trúc nguồn: Nguồn đồng nhất (Homogeneous source), hạt phóng xạ rời rạc (Radioactive particles), hoặc hạt phóng xạ phân bố trong nền ma trận đồng nhất (Radioactive particle in homogeneous source matrix). Thẻ này cho phép định nghĩa hình học nguồn (đường kính, bề dày tâm, bề dày mép, độ mấp mô, mật độ khối), thành phần nguyên tố cấu tạo (tỷ lệ phần trăm khối lượng), và lựa chọn đồng vị từ NuclideLibrary.xml.
• The absorbing layers: Cấu hình các lớp vật liệu hấp thụ nhân tạo đặt giữa nguồn và đầu dò, chỉ định mật độ, bề dày và độ biến thiên bề dày.
• Coincidences and backscattering: Tùy chọn kích hoạt tính toán hiệu ứng trùng phùng (sử dụng thư viện CoincidenceLibrary.xml) và tán xạ ngược trên tấm nền giá đỡ (Backingplate) với các thông số mật độ, bề dày vật liệu nền.
• Simulation: Cung cấp công cụ tính toán hiệu suất ghi nhận hình học (Geometric efficiency calculation) và thiết lập thông số chạy phổ (Run spectrum simulation) bao gồm kích thước gói năng lượng (Energy bin size), bước tích phân (Integration step), và hiệu ứng mở rộng tán xạ Gauss (Include Gaussian straggling).
• Simulation plot: Hiển thị phổ mô phỏng, cho phép chỉnh thang đo trục tuyến tính/logarit, tùy chỉnh màu sắc cho từng đỉnh đồng vị, và hỗ trợ lệnh trích xuất số liệu trực tiếp sang định dạng Microsoft Excel (Copy selected graph in Excel format).

Chương 4: Các kết quả thực nghiệm, mô phỏng và sự so sánh

Chương 4 là phần thực nghiệm trọng tâm, tiến hành nhập các tham số vật lý thực tế vào AASI để mô phỏng và đối chiếu với phổ đo từ hệ Canberra Alpha Analyst. Tác giả khảo sát 3 đối tượng mẫu:

1. Mẫu chuẩn gồm 4 đồng vị: $^{238}\text{U}$, $^{234}\text{U}$, $^{239}\text{Pu}$, $^{241}\text{Am}$.
2. Mẫu Uranium chứa 2 đồng vị: $^{238}\text{U}$ và $^{234}\text{U}$.
3. Mẫu hỗn hợp đồng vị: $^{241}\text{Am}$ và $^{239}\text{Pu}$.

Quá trình mô phỏng được thiết lập với kích thước gói năng lượng $1,0512\text{ keV}$, bước tích phân $0,2\ \mu\text{m}$, bỏ qua lớp hấp thụ trung gian và không xét hiệu ứng tán xạ ngược/trùng phùng. Kết quả tính toán hiệu suất ghi hình học trên AASI cho mẫu chuẩn ở khoảng cách $SDD = 13,0\text{ mm}$ đạt $20,393%$.

Về phía thực nghiệm, phép đo mẫu chuẩn thực hiện trong thời gian $15559\text{ giây}$, thu được tổng số đếm 24192 hạt, tốc độ đếm $1,5549\text{ hạt/s}$, hoạt độ của $^{241}\text{Am}$ xác định được là $1,89453\text{ Bq}$. Phép đo mẫu Uranium thực hiện trong $1500\text{ giây}$ với phương trình đường chuẩn năng lượng $E = 1,019\text{ MeV} + 6,825 \cdot 10^{-4} \times \text{kênh}$, hoạt độ $^{238}\text{U}$ là $1,983\text{ Bq}$ và $^{234}\text{U}$ là $1,750\text{ Bq}$. Khi đưa số liệu mô phỏng và thực nghiệm lên cùng một đồ thị Excel, tác giả nhận thấy vị trí đỉnh năng lượng và độ rộng bán phần của các đồng vị trùng khớp với nhau, chứng minh tính khả thi của việc dùng AASI để mô phỏng hệ phổ kế.

Kết quả và đóng góp

Khóa luận đã tổng hợp và lượng hóa các thông số kỹ thuật thực nghiệm và mô phỏng thông qua các số liệu đo đạc cụ thể:

Thông số	Mẫu chuẩn ($^{238}\text{U}, ^{234}\text{U}, ^{239}\text{Pu}, ^{241}\text{Am}$)	Mẫu Uranium ($^{238}\text{U}, ^{234}\text{U}$)	Mẫu ($^{241}\text{Am}, ^{239}\text{Pu}$)
Loại đầu dò PIPS	1200PIPS (Canberra)	1200PIPS (Canberra)	1200PIPS (Canberra)
Đường kính đầu dò	$39,1\text{ mm}$	$39,1\text{ mm}$	$39,1\text{ mm}$
Độ phân giải FWHM	$37\text{ keV}$	$37\text{ keV}$	$37\text{ keV}$
Bề dày lớp chết (Silic)	$0,05\ \mu\text{m}$	$0,05\ \mu\text{m}$	$0,05\ \mu\text{m}$
Đuôi toàn phần (Total tailing)	0,007666	0,0099551	0,999900
Thành phần nhanh / chậm	$10\text{ keV}$ / $66\text{ keV}$	$10\text{ keV}$ / $66\text{ keV}$	$10\text{ keV}$ / $66\text{ keV}$
Tỷ số diện tích nhanh/chậm	1,2	1,2	1,2
Khoảng cách nguồn - đầu dò (SDD)	$13,0\text{ mm}$	$5,0\text{ mm}$	$9,0\text{ mm}$ (hoặc $30,0\text{ mm}$)
Đường kính nguồn	$24,1\text{ mm}$	$24,1\text{ mm}$	$10,0\text{ mm}$
Bề dày tâm và mép nguồn	$0,01\ \mu\text{m}$	$0,01\ \mu\text{m}$	$0,01\ \mu\text{m}$
Thành phần vật liệu nguồn	N (61%), O (33%), He (3%), Ar (1%)	Hợp chất $\text{UO}_2$ (Uranium, Oxygen)	Hỗn hợp khí tương tự mẫu chuẩn
Mật độ nguồn	$0,0012\text{ g/cm}^3$	$1,88\text{ g/cm}^3$	$0,0213\text{ g/cm}^3$
Thời gian đo thực nghiệm	$15559\text{ s}$	$1500\text{ s}$	—
Tổng số đếm thực nghiệm	24192 hạt	—	—
Hiệu suất ghi nhận hình học tính toán	$20,393%$	Phụ thuộc SDD 5 mm	Phụ thuộc SDD 9 mm
Tổng số đếm mô phỏng mẫu chuẩn	15721 hạt	—	—

Đóng góp của khóa luận bao gồm:

• Thiết lập tài liệu hướng dẫn kỹ thuật chi tiết bằng tiếng Việt đầu tiên về cách sử dụng phần mềm mô phỏng phổ alpha AASI, giúp người nghiên cứu nắm bắt các bước cấu hình thông số nguồn, đầu dò và môi trường.
• Xác định được bộ thông số vật lý chuẩn xác (đặc biệt là các tham số hình dạng đỉnh phổ: Total tailing, Fast/Slow component) tương ứng với đầu dò PIPS A1200-37AM trên hệ Canberra Alpha Analyst.
• Chứng minh được sự phù hợp về mặt hình dạng phổ và vị trí các đỉnh năng lượng giữa phổ tính toán bằng thuật toán Monte Carlo của AASI với phổ thực nghiệm từ hệ đo chuyên dụng.

Hạn chế và hướng nghiên cứu tiếp

Khóa luận phản ánh một số giới hạn kỹ thuật trong phạm vi thực hiện:

• Phép mô phỏng trong nghiên cứu đã bỏ qua ảnh hưởng của các lớp vật liệu hấp thụ trung gian (chọn No absorbing layers) và chưa kích hoạt module tính toán hiện tượng trùng phùng cũng như hiệu ứng tán xạ ngược từ tấm nền (Coincidences & Backscattering).
• Số lượng mẫu thực nghiệm khảo sát chỉ dừng lại ở mẫu chuẩn 4 đỉnh năng lượng và mẫu Uranium tự nhiên; chưa tiến hành khảo sát các loại mẫu môi trường có ma trận phức tạp hoặc mẫu sol khí (aerosol).
• Sự sai khác giữa tổng số đếm thực nghiệm (24192 hạt trong 15559 giây) và số đếm mô phỏng (15721 hạt) bắt nguồn từ việc số phân rã gán cho từng hạt nhân trong tệp mô phỏng được nhập theo các giá trị giả định độc lập thay vì chuẩn hóa tuyệt đối theo hoạt độ thực nghiệm của nguồn.

Hướng mở rộng được xác định từ nội dung nghiên cứu bao gồm:

• Khai thác đầy đủ các tệp dữ liệu phân rã trong CoincidenceLibrary.xml để mô phỏng chi tiết hiện tượng trùng phùng cho các chuỗi phóng xạ có thời gian sống ngắn.
• Áp dụng AASI để tối ưu hóa khoảng cách đo SDD cho các mẫu môi trường thực tế nhằm giảm thiểu thời gian đo mà vẫn kiểm soát được độ mở rộng chân đỉnh phổ do hiện tượng tự hấp thụ trong nguồn dày.

Giá trị tham khảo

Tài liệu này hữu ích cho:

• Sinh viên, học viên chuyên ngành Vật lý hạt nhân, Kỹ thuật hạt nhân và Kỹ thuật bức xạ đang tìm hiểu về nguyên lý phổ kế alpha và phương pháp mô phỏng Monte Carlo.
• Kỹ sư, cán bộ an toàn bức xạ và phân tích phóng xạ môi trường cần tài liệu hướng dẫn vận hành hệ đo Canberra Alpha Analyst kết hợp phần mềm Genie-2000.
• Các nhà nghiên cứu cần tham khảo quy trình tối ưu hóa các tham số hình học và tham số đuôi phổ (Peak tailing) khi giải phổ alpha bằng máy tính. Phần đáng tham khảo nhất là Chương 3 (hướng dẫn cấu hình chi tiết từng thẻ chức năng của AASI) và Chương 4 (bộ số liệu thiết lập thực nghiệm trên đầu dò PIPS A1200-37AM).

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

1. Phần mềm AASI sử dụng thuật toán gì và được lập trình bằng ngôn ngữ nào?

Phần mềm AASI ứng dụng phương pháp mô phỏng thống kê Monte Carlo để mô hình hóa quá trình vận chuyển của hạt alpha qua nguồn, các lớp hấp thụ và đầu dò. Bộ mã nguồn của chương trình được viết bằng ngôn ngữ lập trình Fortran 95, sử dụng các tệp cấu trúc định dạng XML để lưu trữ dữ liệu hạt nhân.

2. Đầu dò được sử dụng trong hệ phổ kế Alpha Analyst của nghiên cứu có đặc điểm gì?

Đầu dò được trang bị là loại bán dẫn PIPS (Passivated Implanted Planar Silicon) mã hiệu A1200-37AM do hãng Canberra chế tạo. Đầu dò có diện tích hoạt động $1200\text{ mm}^2$, bán kính $19,55\text{ mm}$, độ phân giải alpha danh định $37\text{ keV}$, bề dày lớp chết Silic $0,05\ \mu\text{m}$, độ sâu vùng nghèo trên $140\ \mu\text{m}$ và hoạt động ở điện thế phân cực $+40\text{ V}$.

3. Tại sao buồng đo phổ alpha cần phải hút chân không tới áp suất dưới 0,1 torr?

Hạt alpha là hạt mang điện nặng, có độ ion hóa riêng rất cao nên bị mất năng lượng nhanh chóng và có quãng chạy rất ngắn trong không khí (chỉ từ 3 đến 7 cm đối với dải năng lượng 4 - 9 MeV). Việc hút chân không xuống dưới 0,1 torr giúp loại bỏ môi trường vật chất giữa nguồn và đầu dò, ngăn ngừa hiện tượng mất năng lượng do va chạm với phân tử khí, từ đó tránh làm suy giảm độ phân giải và méo mó hình dạng đỉnh phổ.

4. Hiệu suất ghi nhận hình học của mẫu chuẩn được tính toán trong luận văn đạt giá trị bao nhiêu?

Với thiết lập khoảng cách từ nguồn tới đầu dò là $SDD = 13,0\text{ mm}$, đường kính đầu dò $39,1\text{ mm}$ và đường kính nguồn $24,1\text{ mm}$, chức năng tính toán hiệu suất hình học của phần mềm AASI xác định hiệu suất ghi nhận của mẫu chuẩn đạt $20,393%$.

5. Dữ liệu phân rã và hiện tượng trùng phùng được AASI truy xuất từ các tệp tin nào?

AASI lưu trữ và truy xuất các đặc trưng phân rã phóng xạ alpha thông qua tệp NuclideLibrary.xml. Đối với việc kích hoạt tính toán hiện tượng trùng phùng hạt nhân và tán xạ giật lùi, phần mềm truy xuất từ cơ sở dữ liệu bổ sung có tên CoincidenceLibrary.xml.

Kết luận

Khóa luận của tác giả Lâm Thu Văn đã hệ thống hóa cơ sở lý thuyết phân rã alpha, cấu trúc phần cứng hệ đo Canberra Alpha Analyst và hoàn thiện quy trình thiết lập tham số trên phần mềm mô phỏng Monte Carlo AASI. Kết quả đối sánh trên đồ thị năng lượng cho thấy phổ mô phỏng bằng AASI có sự tương thích tốt về vị trí đỉnh và hình dạng phổ so với phổ thực nghiệm thu được từ đầu dò PIPS A1200-37AM đối với các mẫu chuẩn và mẫu Uranium. Công trình khẳng định khả năng ứng dụng hiệu quả của bộ mã AASI trong việc hỗ trợ phân tích phổ kế alpha, tính toán hiệu suất hình học và tối ưu hóa điều kiện đo đạc mẫu phóng xạ môi trường tại Việt Nam.

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

CHƯƠNG 1 : TONG QUAN VE HAT ALPHA [1], [2].1 Khái niệm hạt alpha Hat Alpha hay tia alpha là một dạng của phóng xạ. Đó là hạt ion hóa cao và khả năng đâm xuyên thấp. Hat alpha gồm hai proton và hai neutron liên kết với nhau thành một hạt giống hệt hạt nhân nguyên tử hellium, do đó, hat alpha có thê được viết là He””. Hat alpha xuất hiện trong phân rã của hạt nhân phóng xạ như là uranium, radium,.

trong một quá trình gọi là phân rã alpha. Sự phân rã làm hạt nhân ở trạng thái kích thích và trở vẻ trạng thái cơ ban bằng cách phát tia gamma.2 Sự phân rã alpha Xay ra khi hạt nhân phóng xạ có tỉ số > quá thấp. Khi phan rã alpha, hạt nhân ban dau 4X chuyển thành hạt nhân {4Y và phát ra hạt alpha. AX —*> AY + tHe (1.1) Vẻ quan hệ khôi lượng, phân rã alpha thỏa mãn điền kiện sau đây: Af„ =À(_ +m, +2m +@Q (1.2) trong đó M,, Af_,m, vam, tương ứng là khối lượng các nguyên tử mẹ, nguyên tử con, hạt nhân alpha và hat clcctron.

Ø là khỗi lượng tương đương với năng lượng tong cộng giải phóng khi phan ra, bằng tông động năng của hạt nhân con và hat alpha. Hai hat electron quỹ đạo bị mat đi khi hạt nhân mẹ phân rã ra hạt nhân con có số nguyên tử thấp hơn. Hạt alpha phát ra với năng lượng cô định và suất ra cô định. Ví dụ như quá trình phân rã của 28 Ra —*-» Rn + Ÿ//c, gồm hai nhánh phát alpha, nhánh thứ nhất với hạt alpha năng lượng 4,602 MeV, suất ra 5,5% và nhánh thứ hai với hạt alpha năng lượng 4,785 MeV, suất ra 94,4%.

Hạt nhân Rn sau phân rã theo nhánh thứ nhất nằm ở trạng thái kích thích và tiếp tục phân rã gamma dé chuyên về trạng thái cơ bản. Hạt nhân ?°#z sau phân rã theo nhánh thứ hai nằm ở trạng thái cơ bản. Sơ đồ phân rã của * 2 [9] Hiện nay người ta đã biết hơn 200 hạt nhân phân rã alpha, Phân rã alpha chủ yeu xảy ra với các hạt nhân nặng ở cuối bảng tuần hoàn Mendeleev với Z>83. Chang hạn thallium (Z=81) không có đồng vị nao phat alpha, chi (Z=82) có hai đồng vị phát alpha, con polonium (Z=§4) có không ít hơn hai mươi mốt đông vị phát alpha,.

Ngoài ra có một nhóm nhỏ các hạt nhân trong vùng đất hiếm cũng phân rã alpha với A=140 đến 160, hạt nhân nhẹ nhất là ‘Ce chứa 84 neutron. Một trường hợp đặc biệt là hạt nhân nhẹ {Be cũng phat alpha với thời gian bán rã 3.3 Các đặc trưng chính của phân rã alpha Có ba đặc trưng chính là thời gian bán rã Tị¿; của hạt nhân trước phân rã, động năng E và quãng chạy R của hạt alpha. Thời gian bán rã của các hạt nhân phân rã alpha thay đổi trong một dãy rất rộng. Chang hạn “Pb có T;2=1,4.10” năm, “Rn có T,;=10 s, Năng lượng các hạt bay ra chỉ thay đôi trong một dãy hẹp.

Đối với các hạt nhân nặng thì năng lượng các hat alpha thay đổi từ 4 MeV đến 9 MeV, còn doi với nhóm đất hiểm từ 2 MeV đến 4,5 MeV. Quang chạy của hat alpha rat ngắn đo nó bj hap thụ rat mạnh trong vật chat.4 Sự liên hệ giữa thời gian bán rã T,2, động năng E và quãng chạy R 1.1 Thời gian bán rã T;„ và động năng E Có sự phụ thuộc rất mạnh của thời gian bán rã T;; vào động năng E của hat alpha, đây là tinh chat quan trọng nhất của các hạt nhân phân rã alpha. Chang hạn nếu giảm 1% năng lượng thì có the làm tăng thời gian bán rã lên một bậc. Nếu giảm 10% năng lượng thì thời gian bán rã thay đỗi từ hai đến ba bậc.

Sự phụ thuộc nay tuần theo định luật Geiger-Nuttal như sau: D lg7,, =C+—— (1.3) VE trong đó C va D là các hing số không phụ thuộc vào số khối A ma chỉ phy thuộc vao điện tích Z. Vi dụ nếu E tính theo MeV thi: Với Z=84 ta có C=-50,15 và D=128,8 Với Z=90 ta có C=-51,94 và D=139,4 1.2 Quang chạy R và động lượng E Trong phép gan đúng đầu tiên có thé sử dụng công thức liên hệ giữa quãng chạy của hat alpha trong không khí (tính theo cm) va năng lượng của no (tính theo MeV) trong miễn nang lượng 4-9 MeV hay quãng chạy trong không khí 3-7 cm như sau: Ry, ~0.318E*2 Trong môi trường với khối lượng hạt nhân A thì quãng chạy hạt alpha đưếỳ tính theo công thức: R=0,56R,, 4° (1.5) Nhờ hai công thức nay có thé tinh được quãng chạy khi đo được năng lượng hay ngược lại đối với hạt alpha.5 Điều kiện về năng lượng đối với phân rã alpha Xét quá trình phân rã: ⁄X —*> PUY + ‡/fe Điều kiện dé xảy ra phân rã alpha là: AE = Ey.6) Tức là năng lượng liên kết của hạt nhân mẹ E, ,, phải nhỏ hơn tong năng lượng liên kết của hạt nhân con #„ „,„ „ vả hạt alpha £, „. Năng lượng liên kết của hạt alpha là 28 MeV nên năng lượng liên kết riêng trên một nucleon 1a 7 MeV. Như vậy, dé phân rã alpha xảy ra thì năng lượng liên kết riêng của hạt nhân mẹ phải nhỏ hơn 7 MeV.

Điều kiện nay không cho phép các hạt nhân nhẹ phân rã alpha vì năng lượng liên kết riêng của chúng vào cỡ 8 MeV. Còn năng lượng liên kết riêng đối với các hạt nhân nặng giảm xuống do năng lượng đây Coulomb của các proton. Năng lượng phân rã AE phân bố cho động năng hạt nhân con #, ,„ ; và động năng hat alpha AE=E.7) Gọi động lượng của hạt nhân mẹ là p,„, hạt nhân con là p,,; „ vả hat alpha là p,, giả thiết là hạt nhân mẹ đứng yên khi phân rã p„; =0, thi theo định luật bảo toàn động lượng ta có P¿:; ; = P„- Từ đó nhận được: mn Eusza EB, = (1.8) Pu; trong đó m, và m,,,, là khối lượng của hat alpha và hạt nhân con. Kết hợp hai biểu thức (1.8) ta được: E = AE MN 4.5 7-2 (1:9) Megyg TM, Do khối lượng hạt nhân con rất lớn so với khối lượng hạt alpha nên theo công thức (1.

tức là hạt alpha mang phan lớn năng lượng tỏa ra khi phân rã. Vi dụ, năng lượng phân rã alpha “Bi(ThC)—— ` THThC )+ø là AE =6,203 MeV. Năng lượng này phân bố cho hạt alpha theo công thức (1.9) là £ =6,086 MeV và cho hạt nhân con “STI(TAC ) theo công thức (1.6 Cơ chế phân rã alpha Trong cơ chế phân rã alpha cần tính đến ba yếu tổ là trường thé Coulomb quanh hạt nhân, lực ly tâm và cấu trúc hạt nhân.1 Trường thé Coulomb và hiệu ứng đường ngầm Đề giải thích sự phụ thuộc rất mạnh của thời gian ban rã T,,; của hạt nhân vào năng lượng hạt alpha, người ta xem xét cơ chế đẻ hat alpha thoát ra khỏi hạt nhân. Giả thuyết gần đúng nhất là xem hạt alpha hình thành vả ton tại trong hạt nhân trước khi thoát ra khỏi hạt nhân.

Hạt alpha mang điện tích dương +2e nên ngoài lực tương tác hạt nhân, nó còn chịu tác dụng của lực Coulomb. Đề giái bai toán tương tác nay ta gia sử hat alpha từ ngoài đi vào hạt nhân. Thế Coulomb do hạt nhân tương tác lên nó tăng tỉ lệ nghịch với khoảng cách r theo biểu thức: 2Ze? Ú,„.10) Hàm này được the hiện trên hình sau: Hình 1. Thé tương tac hạt nhân và thé Coulomb của hat alpha Thế tương tác nảy tăng dan trong miễn ngoài bán kính hạt nhân, tại đó lực hạt nhân bang không.

Đến biên hạt nhân + = # thi lực hạt nhân đóng vai trò quan trọng và đường biểu diễn giảm đột ngột theo đường thăng đứng. dạng thế bên trong hạt nhân z < R chưa được biết tường tận, ở đây giả thuyết nó có dang hỗ chữ nhật với thế không đôi bên trong hạt nhân. Chiều cao bờ thế Coulomb tại r= # = 10” em và với Z = 100 là: U„ = #2“ =30 MeV r (111) Như đã trình bày trên day, hat alpha phân rã từ các hạt nhân nặng có năng lượng từ 4 MeV đến 9 MeV, tức là nhỏ hơn chiêu cao hang rảo thé. Theo cơ học cô điện thi hat alpha không thé vượt qua hàng rào thể để ra ngoài, tức là không thẻ xảy ra quá trình phân rã alpha.

Tuy nhiên, trong thế giới vi mô, theo cơ học lượng tử, hat alpha có thể xuyên qua hang rao Coulomb theo cơ chế đường ngâm. Đề giải bài toán về hiệu ứng đường ngắm của hạt alpha, một lần nữa ta đơn giản hóa bang cách tính xác suất hat alpha có năng lượng E vượt qua bở thé hình chữ nhật có độ cao U và độ rộng d. Trạng thái của hạt alpha được mô tả bởi hàm (r) thỏa mãn phương trình Schrédinger: (1.12) - ? - - Ay(r) +2 )) =0 ] E~U0wr h trong đó: eo FF & fos ¬. F : R : A ===++==*==-: toán từ đạo ham riêng bậc hai theo không gian Ox” Gy" ae” m: khối lượng hat alpha fh: hằng số Planck Giải phương trình Schrodinger với rao thế có dang chữ nhật ta sẽ tìm được hệ số truyền qua =ep|-[ 28 - Ey (1.13) Nếu hat alpha trong hạt nhân có vận tốc v thì nó đi đến bờ thé trung bình v/ & lần trong | giây.

Nhu vậy, hằng số phan ra alpha bằng: w v [2m 3— vis Vận tốc y có thé xác định từ hệ thức bất định Heisenberg giữa động lượng vả tọa độ, tức lả h vz=—. Vậ mR ¥ “aifPtw- bar | 2m (1.15) Thời gian bán rã 7;,, = a „ đo đó nó phụ thuộc rat mạnh vào bán kính hạt nhân R. Ta thử đánh giá bậc của thời gian bán rã, với U,-E=20 MeV, «2. Khi đó D=e" ~10'5, Hệ số của hảm =—— = 10".== 10x z 10 năm a Thời gian nay là hợp lí vì nó vào cỡ thời gian bán rã của 2”.2 Vai trò của bờ thé li tâm Trong các tính toán trên ta coi hạt alpha bay ra với momen quỹ đạo £ =0.

Nếu hạt alpha bay ra với £ +0 thì nó phải vượt qua bờ thé li tâm bổ sung ngoài thé Coulomb: h/(+l) U, i = : (1.16) 2mr- Bờ thé li tâm nay không lớn do nó giảm theo ham = trong lúc bờ thé Coulomb giảm chậm F hơn theo ham 5 nhưng do độ thay đổi này con chia cho hang số Planck trong ham số mũ nên nó ' làm tăng đáng ké thời gian bán rã của hạt alpha.1 dưới đây sẽ dan ra hệ số K suy giảm của hằng số phân ra 4 với các giá trị khác nhau của ¢ đổi với trường hợp £=5 MeV và R=9,6.