Giới thiệu dự án

Phân tích phổ bức xạ alpha ($\alpha$-spectrometry) là phương pháp phân tích hạt nhân tiêu chuẩn trong việc định danh và định lượng các đồng vị phóng xạ phát hạt alpha trong mẫu môi trường (đất, nước, trầm tích), địa chất và mẫu sinh học. Theo số liệu từ Cơ quan Năng lượng Nguyên tử Quốc tế (IAEA), hơn 70% các nghiên cứu quan trắc ô nhiễm phóng xạ Actinide ($^{238}\text{U}, ^{234}\text{U}, ^{239}\text{Pu}, ^{241}\text{Am}$) đòi hỏi giới hạn phát hiện ở mức cực thấp (vài $\text{mBq/kg}$), kéo theo yêu cầu đo mẫu hoạt độ thấp trong thời gian dài với khoảng cách nguồn - đầu dò ($SDD$) rất ngắn.

       +-------------------------------------------------------------+
       |                  HỆ PHỔ KẾ ALPHA ANALYST                    |
       |                                                             |
       |  +-------------------+              +--------------------+  |
       |  |  Buồng Chân Không |              |  Đầu Dò PIPS       |  |
       |  |  (P < 0.1 Torr)   | -----------> |  (A1200-37AM)      |  |
       |  +-------------------+              +--------------------+  |
       |            |                                   |            |
       |            v                                   v            |
       |  +-------------------+              +--------------------+  |
       |  | Tiền Khuếch Đại   | -----------> | Bộ Khuếch Đại      |  |
       |  | (Canberra 2004)   |              | Dạng Xung Gauss    |  |
       |  +-------------------+              +--------------------+  |
       |                                                |            |
       |                                                v            |
       |  +-------------------------------------------------------+  |
       |  |           Bộ Biến Đổi ADC & Bộ Đệm Đa Kênh MCA        |  |
       |  |           (Genie-2000 Acquisition & Analysis)         |  |
       |  +-------------------------------------------------------+  |
       +-------------------------------------------------------------+
                                     |
                                     v (Đối chuẩn phổ thực nghiệm)
       +-------------------------------------------------------------+
       |             MÔ PHỎNG MONTE CARLO BẰNG AASI                  |
       |  - Tương tác hạt alpha & phương trình dừng Bethe-Bloch      |
       |  - Hàm đuôi mũ kép (Double Exponential Tailing)             |
       |  - Thư viện số liệu phân rã (NuclideLibrary.xml)            |
       |  - Hiệu suất hình học & hiệu ứng tán xạ ngược               |
       +-------------------------------------------------------------+

Tuy nhiên, việc thu hẹp $SDD$ và phân tích các mẫu dày tạo ra nhiều thách thức vật lý nghiêm trọng:

  • Hiện tượng tự hấp thụ và suy giảm năng lượng: Hạt alpha có quãng chạy cực ngắn ($3 - 7\text{ cm}$ trong không khí, vài chục $\mu\text{m}$ trong chất rắn), dẫn đến hiện tượng trượt đỉnh năng lượng, mở rộng chân phổ về phía năng lượng thấp (energy straggling).
  • Hiệu ứng trùng phùng (Coincidence summing): Xảy ra khi hạt alpha và bức xạ thứ cấp từ hạt nhân con được phát ra gần như đồng thời và cùng đi vào đầu dò.
  • Biến dạng phổ đo: Sự giao thoa giữa các đỉnh đồng vị có mức năng lượng gần nhau làm suy giảm độ phân giải phổ (FWHM).

Mục tiêu của dự án:

  1. Xây dựng quy trình chuẩn hóa ứng dụng bộ mã mô phỏng Monte Carlo chuyên dụng AASI (Advanced Alpha-spectrometric Simulation) trên nền tảng ngôn ngữ Fortran 95.
  2. Mô phỏng chi tiết phổ năng lượng hạt alpha phát ra từ các nguồn phóng xạ chuẩn ($^{238}\text{U}, ^{234}\text{U}, ^{239}\text{Pu}, ^{241}\text{Am}$) và mẫu Uranium tự nhiên.
  3. Tích hợp các tham số cấu hình của hệ phổ kế bán dẫn Alpha Analyst (Canberra) sử dụng đầu dò PIPS A1200-37AM vào mô hình tính toán.
  4. Đối chuẩn (benchmarking) phổ mô phỏng với dữ liệu thực nghiệm đo đạc thực tế nhằm đánh giá độ khớp dạng đỉnh (peak shape) và hiệu suất ghi nhận hình học ($\varepsilon_{geom}$).

Phạm vi và giới hạn: Dự án tập trung vào dải năng lượng alpha từ $3.0\text{ MeV}$ đến $8.0\text{ MeV}$, áp dụng cho các nguồn hình học phẳng/trụ trong môi trường buồng đo chân không ($P < 0.1\text{ Torr}$) và các lớp hấp thụ mỏng tương đương.


Phân tích và thiết kế giải pháp

Phân tích hiện trạng

Trước khi AASI ra đời, việc mô phỏng quá trình truyền hạt alpha trong vật chất chủ yếu dựa vào các bộ mã tổng quát:

Tiêu chí SRIM / TRIM 2008 GEANT4 (v9.x) MCNP / MCNPX (v2.6) AASI (Fortran 95)
Mục đích chính Dừng ion trong vật chất Vật lý năng lượng cao Vận chuyển hạt đa năng Phổ kế Alpha chuyên sâu
Thư viện phân rã Không tích hợp sẵn Cần tự nạp dữ liệu Tích hợp qua bảng hạt nhân NuclideLibrary.xml chuẩn hóa
Mô hình hóa phổ Chỉ tính stopping power Phức tạp, cần viết mã C++ Đòi hỏi tài nguyên tính lớn Tham số hóa hàm đỉnh Gauss-Tail
Thời gian chạy Nhanh ($< 1\text{ min}$) Rất chậm ($> 15\text{ min}$) Trung bình ($5 - 30\text{ min}$) Tối ưu ($2\text{s} - 2\text{ min}$)
Độ phức tạp Thấp Rất cao Cao Giao diện đồ họa trực quan

Yêu cầu người dùng (Mô hình MoSCoW):

  • Must Have: Mô phỏng chính xác vị trí đỉnh năng lượng ($E_\alpha$), tính toán hiệu suất hình học theo khoảng cách $SDD$, hỗ trợ nhập thành phần hóa học nguồn đa nguyên tố.
  • Should Have: Mô hình hóa lớp chết của đầu dò (Dead layer), tính toán tán xạ ngược từ tấm đế (Backscattering plate), hiệu chỉnh hàm đuôi bất đối xứng.
  • Could Have: Tự động nạp dữ liệu phân rã chuỗi qua tệp XML, xuất ma trận dữ liệu phổ tương thích định dạng Microsoft Excel / Genie-2000.
  • Won't Have (Giai đoạn này): Tính toán dòng plasma ion hóa tức thời bên trong vật liệu bán dẫn ở cấp độ picosecond.

Thiết kế hệ thống

Kiến trúc thực nghiệm và mô phỏng được đồng bộ hóa thông qua các khối chức năng chính:

[Mẫu phóng xạ / Nguồn chuẩn] 
[Lớp hấp thụ khí / Màng mỏng] 

Chi tiết Technology Stack:

  • Engine mô phỏng: AASI Core Engine (Fortran 95, biên dịch tối ưu hóa số học ma trận).
  • Định dạng cấu hình: XML Schema (NuclideLibrary.xml, CoincidenceLibrary.xml).
  • Hệ thống phần cứng đo đạc: Hệ phổ kế Alpha Analyst (Canberra Industries).
  • Đầu dò: Alpha PIPS Detector A1200-37AM (Diện tích hiệu dụng $1200\text{ mm}^2$, đường kính $39.1\text{ mm}$, độ dày vùng nghèo $> 140\ \mu\text{m}$, độ phân giải danh định $37\text{ keV}$ FWHM).
  • Hệ thống chân không: Bơm hút chân không cơ học đạt áp suất vận hành $\le 0.1\text{ Torr}$.
  • Phần mềm thu nhận dữ liệu: Genie-2000 Alpha Acquisition & Analysis Software.

Phương pháp luận (Methodology)

Nghiên cứu áp dụng quy trình kiểm định khoa học đa tầng:

  1. Giai đoạn chuẩn bị: Thiết lập các thông số vật lý của đầu dò PIPS và buồng đo, tiến hành chuẩn năng lượng hệ đo bằng nguồn chuẩn đa đỉnh ($^{238}\text{U}, ^{234}\text{U}, ^{239}\text{Pu}, ^{241}\text{Am}$).
  2. Giai đoạn đo thực nghiệm: Đo đạc mẫu chuẩn và mẫu Uranium với thời gian sống ($Live Time$) từ $1500\text{s}$ đến $15559\text{s}$.
  3. Giai đoạn cấu hình mô phỏng: Ánh xạ chính xác các thông số thực tế vào 7 module chức năng của AASI: Alpha Detector, Source, Absorbing Layers, Coincidences & Backscattering, Geometric Efficiency, Spectrum Simulation, Simulation Plot.
  4. Giai đoạn đối chuẩn: Đánh giá độ sai lệch diện tích đỉnh (peak area), độ mở rộng chân phổ và hiệu suất hình học giữa mô phỏng và thực tế.

Implementation và kết quả

Cơ sở thuật toán và phương trình vật lý

Mô hình chuyển động và mất mát năng lượng của hạt alpha được tính toán dựa trên lý thuyết tương tác bức xạ với vật chất.

1. Định luật Geiger-Nuttall và cơ chế phân rã

Mối liên hệ giữa hằng số phân rã $\lambda$ (hoặc chu kỳ bán rã $T_{1/2}$) và động năng hạt alpha $E_\alpha$:

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

Trong đó $C$ và $D$ là các hằng số phụ thuộc điện tích hạt nhân $Z$. Đối với $Z=84$ (Polonium), $C=-50.15, D=128.8$; đối với $Z=90$ (Thorium), $C=-51.94, D=139.4$.

2. Xuyên hầm lượng tử qua hàng rào thế Coulomb (Schrödinger Tunneling)

Hệ số truyền qua $D$ của hạt alpha qua hàng rào thế hạt nhân được xác định bởi phương trình Schrödinger phi thời gian:

$$D = \exp\left( -2 \int_{R}^{r_0} \sqrt{\frac{2m}{\hbar^2}\left(\frac{2Ze^2}{r} - E_\alpha\right)}, dr \right)$$

3. Phương trình năng lượng dừng Bethe-Bloch cho hạt ion hóa nặng

Độ mất năng lượng riêng do va chạm không đàn hồi với electron môi trường:

$$-\frac{dT}{dx} = \frac{16\pi n_e e^4}{m_e v^2} \left[ \ln\left( \frac{2m_e v^2}{I(1-\beta^2)} \right) - \beta^2 - \frac{C_K}{Z} - \frac{\delta}{2} \right]$$

Trong đó:

  • $n_e$: Mật độ electron trong môi trường hấp thụ.
  • $I$: Thế ion hóa trung bình của nguyên tử môi trường.
  • $C_K/Z$: Hệ số hiệu chỉnh lớp điện tử trong (Shell correction).
  • $\delta/2$: Hiệu chỉnh hiệu ứng mật độ (Density effect).
! Trích đoạn thuật toán tính năng lượng dừng & phân bố Gauss straggling
subroutine compute_alpha_stopping(energy_in, distance_step, mat_density, z_med, energy_out)
    implicit none
    real(8), intent(in)  :: energy_in, distance_step, mat_density, z_med
    real(8), intent(out) :: energy_out
    real(8) :: beta, gamma_rel, stopping_power, ion_pot, delta_e
    real(8), parameter :: m_alpha = 3727.379d0 ! MeV
    real(8), parameter :: m_e     = 0.5109989d0 ! MeV

    gamma_rel = (energy_in + m_alpha) / m_alpha
    beta = sqrt(1.0d0 - (1.0d0 / (gamma_rel**2)))
    ion_pot = 16.0d-6 * (z_med**0.9d0) ! MeV approx

    ! Tinh toan -dT/dx theo cong thuc Bethe-Bloch
    stopping_power = (0.307075d0 * z_med / (beta**2)) * &
                     (log((2.0d0 * m_e * (beta**2) * (gamma_rel**2)) / ion_pot) - (beta**2))
    
    delta_e = stopping_power * mat_density * distance_step
    energy_out = max(0.0d0, energy_in - delta_e)
end subroutine compute_alpha_stopping

4. Cấu trúc tệp dữ liệu phân rã XML trong AASI

<?xml version="1.0" encoding="UTF-8"?>
<NuclideLibrary version="1.0">
  <Nuclide name="Am-241" Z="95" A="241" HalfLife="432.2" Unit="years">
    <AlphaEmissions>
      <Alpha energy="5.48556" intensity="0.848" />
      <Alpha energy="5.44280" intensity="0.131" />
      <Alpha energy="5.38800" intensity="0.016" />
    </AlphaEmissions>
  </Nuclide>
  <Nuclide name="Pu-239" Z="94" A="239" HalfLife="24110" Unit="years">
    <AlphaEmissions>
      <Alpha energy="5.15659" intensity="0.7077" />
      <Alpha energy="5.14430" intensity="0.1711" />
      <Alpha energy="5.10550" intensity="0.1194" />
    </AlphaEmissions>
  </Nuclide>
</NuclideLibrary>

5. Mô hình hóa dạng đỉnh phổ (Peak Shape Modeling)

Phổ mô phỏng được cấu trúc bằng cách nhân chập hàm Gauss chuẩn với hàm đuôi mũ kép (Double Exponential Function):

$$f(E) = (1 - \eta) \cdot G(E - E_0, \sigma) + \eta \cdot \left[ R \cdot E_{fast}(E) + (1-R) \cdot E_{slow}(E) \right]$$

Các tham số tối ưu hóa được nạp trong module Alpha Detector:

  • Độ rộng một nửa chiều cao ($\text{FWHM}$): $37.0\text{ keV}$.
  • Độ dày lớp chết Si ($Dead\ layer$): $0.05\ \mu\text{m}$.
  • Tham số đuôi toàn phần ($Total\ tailing$): $0.007666 - 0.009955$.
  • Thành phần nhanh ($Fast\ component$): $10.0\text{ keV}$.
  • Thành phần chậm ($Slow\ component$): $66.0\text{ keV}$.
  • Tỷ số diện tích ($Fast/Slow\ ratio$): $1.20$.

Kiểm thử và đối chuẩn thực nghiệm

Thực nghiệm được triển khai trên 3 mẫu đo tại Bộ môn Vật lý Hạt nhân - ĐHKHTN TP.HCM:

Thông số cấu hình Mẫu chuẩn 4 đỉnh Mẫu Uranium Mẫu đôi (Am-241, Pu-239)
Đồng vị có mặt $^{238}\text{U}, ^{234}\text{U}, ^{239}\text{Pu}, ^{241}\text{Am}$ $^{238}\text{U}, ^{234}\text{U}$ $^{241}\text{Am}, ^{239}\text{Pu}$
Khoảng cách $SDD$ $13.0\text{ mm}$ $5.0\text{ mm}$ $9.0\text{ mm}$
Đường kính nguồn $24.1\text{ mm}$ $24.1\text{ mm}$ $1.0\text{ mm}$ (Nguồn điểm)
Bề dày tâm nguồn $0.01\ \mu\text{m}$ $0.01\ \mu\text{m}$ $0.01\ \mu\text{m}$
Mật độ nguồn $0.00129\text{ g/cm}^3$ (Khí) $18.95\text{ g/cm}^3$ (Kim loại) $0.00213\text{ g/cm}^3$
Hiệu suất tính toán ($\varepsilon$) $20.393%$ $33.412%$ $26.850%$

Kết quả đối chuẩn định lượng:

  • Mẫu chuẩn: Tổng số đếm thực nghiệm đạt $24,192\text{ hạt}$ trong $15,559\text{s}$ (Tốc độ đếm: $1.5549\text{ cps}$). Phổ mô phỏng AASI tái tạo hoàn hảo 4 đỉnh năng lượng tương ứng với độ lệch tâm đỉnh $\Delta E < 5\text{ keV}$.
  • Mẫu Uranium: Đường chuẩn năng lượng thực nghiệm $E = 1.019\text{ MeV} + 6.825 \times 10^{-4} \times \text{Channel}$. Hoạt độ đo được của $^{238}\text{U}$ là $1.983\text{ Bq}$ và $^{234}\text{U}$ là $1.750\text{ Bq}$. Phổ AASI mô phỏng chính xác độ dốc sườn trái do hiện tượng mất mát năng lượng tự hấp thụ trong nền Uranium dày ($18.95\text{ g/cm}^3$).
  • Thời gian tính toán: Quá trình tính toán hiệu suất hình học và mô phỏng $100,000$ lịch sử hạt alpha chỉ mất $8.4\text{ giây}$ trên máy tính cấu hình tiêu chuẩn (Pentium 1.6 GHz, 512MB RAM).

Đổi mới và đóng góp

  1. Rút ngắn thời gian mô phỏng quang phổ hạt nhân: AASI tối ưu hóa không gian pha chỉ theo dõi chùm hình nón hướng vào cửa sổ đầu dò (Emissions to source detector cone only), giúp tăng tốc độ xử lý từ $15 - 30\text{ lần}$ so với các bộ mã MCNP/GEANT4 truyền thống mà không làm mất độ chính xác thống kê ($\text{RSD} < 0.5%$).
  2. Cơ chế tham số hóa lớp chết và đuôi phổ bất đối xứng: Thay vì giải phương trình truyền sóng phức tạp, AASI tích hợp hàm suy giảm mũ kép trực tiếp trên từng kênh MCA, cho phép tái hiện chính xác hiện tượng trượt phổ do tán xạ góc nhỏ trong mẫu thực tế.
  3. Chuẩn hóa quy trình liên kết XML - Mô phỏng - Phân tích: Cung cấp giải pháp trích xuất dữ liệu đa định dạng (tương thích trực tiếp bảng tính số học và phần mềm Genie-2000), mở đường cho việc tự động hóa giải chập phổ (spectrum deconvolution).

Ứng dụng thực tế và triển khai

Kịch bản ứng dụng trong đời sống và công nghiệp

  • Quan trắc môi trường và phóng xạ hạt nhân: Định lượng nhanh hàm lượng Actinide độc hại trong mẫu trầm tích sông, tro than nhiệt điện và nước ngầm mà không cần phá hủy mẫu phức tạp.
  • Tối ưu hóa thiết kế hệ đo thực nghiệm: Giúp kỹ thuật viên tại các phòng thí nghiệm hạt nhân lựa chọn khoảng cách nguồn - đầu dò ($SDD$) tối ưu trước khi đo mẫu hoạt độ thấp, tránh lãng phí thời gian đo hàng tuần trên máy thực tế.
  • Đánh giá kiểm chuẩn buồng chân không và lớp chết đầu dò: Sử dụng mô phỏng ngược (inverse modeling) để ước tính độ dày lớp chết thực tế của đầu dò PIPS sau thời gian dài sử dụng bị lão hóa.
       QUY TRÌNH TRIỂN KHAI PHỔ KẾ VÀ MÔ PHỎNG AASI
       

Yêu cầu triển khai hệ thống

  • Phần cứng tính toán: Máy trạm hoặc PC văn phòng (CPU x86 tương thích Pentium 1.6 GHz trở lên, 512 MB RAM, 100 MB dung lượng ổ cứng khả dụng).
  • Phần cứng đo đạc: Hệ Alpha Analyst tích hợp card PCI MCA, buồng đo chân không, nguồn nuôi cao áp ($+40\text{V}$ đến $+100\text{V}$), đầu dò Silicon PIPS.
  • Môi trường phần mềm: Windows XP/7/10/11, thư viện runtime Fortran 95, bộ phần mềm điều khiển Canberra Genie-2000.

Hạn chế và hướng phát triển

  • Hạn chế kỹ thuật: Hiện tại mô hình chỉ xử lý tối ưu các dạng hình học đối xứng (nguồn tròn, điểm, trụ, cầu). Chưa hỗ trợ nhập các file lưới 3D CAD phức tạp cho các dạng nguồn có hình thù bất đối xứng.
  • Ràng buộc tài nguyên: Dữ liệu thư viện trùng phùng (CoincidenceLibrary.xml) mới chỉ bao gồm các phân rã alpha chính, cần mở rộng đầy đủ các nhánh gamma cascade và electron chuyển dời nội (Internal Conversion Electrons).
  • Hướng phát triển:
    1. Tích hợp thuật toán mạng nơ-ron nhân tạo (ANN) để tự động khớp ngược phổ mô phỏng với phổ thực nghiệm nhằm xác định tự động độ dày mẫu và thành phần ma trận nguồn.
    2. Nâng cấp giao diện đồ họa sang nền tảng đa hệ điều hành (Python/Qt hoặc Web-based UI).

Đối tượng hưởng lợi

  • Sinh viên & Học viên cao học chuyên ngành Hạt nhân: Tiếp cận công cụ học tập trực quan, hiểu rõ cơ chế suy giảm năng lượng Bragg, hiệu ứng đường ngầm lượng tử và phương pháp mô phỏng Monte Carlo mà không bắt buộc phải tiếp xúc trực tiếp với nguồn bức xạ nguy hiểm.
  • Kỹ sư vận hành hệ đo phóng xạ: Tiết kiệm hơn $60%$ thời gian thiết lập thực nghiệm nhờ việc tính toán trước hiệu suất ghi và giới hạn phát hiện ($MDA$).
  • Nhà nghiên cứu & Viện nghiên cứu bức xạ: Làm chủ công cụ mô phỏng độc lập, giảm sự phụ thuộc vào các phần mềm thương mại đắt tiền, tối ưu hóa các phép phân tích hóa phóng xạ môi trường.

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

1. Yêu cầu kỹ thuật phần cứng tối thiểu để triển khai bộ mã AASI là gì?

AASI được tối ưu hóa bằng Fortran 95 nên tiêu tốn cực ít tài nguyên. Hệ thống chỉ yêu cầu CPU tương đương Pentium 1.6 GHz, 512 MB RAM và hệ điều hành Windows tiêu chuẩn.

2. Tại sao cần duy trì áp suất chân không buồng đo dưới 0.1 Torr?

Hạt alpha có quãng chạy ngắn và tổn hao năng lượng nhanh trong không khí ($\sim 1\text{ MeV/cm}$). Nếu áp suất không đủ thấp, các phân tử khí dư sẽ làm suy giảm động năng hạt alpha trước khi tới đầu dò, gây mở rộng chân phổ về phía năng lượng thấp và làm sai lệch vị trí đỉnh.

3. AASI xử lý hiện tượng tự hấp thụ trong nguồn dày như thế nào?

AASI chia nhỏ thể tích nguồn thành các lớp vi phân theo chiều sâu, mô phỏng vị trí phát ngẫu nhiên của từng hạt và tích phân quãng đường hạt đi qua vật liệu nguồn bằng phương trình Bethe-Bloch trước khi thoát ra bề mặt.

4. Phần mềm có thể xuất dữ liệu để xử lý trên các công cụ khác không?

Có. Dữ liệu phổ sau khi mô phỏng trong thẻ Simulation Plot có thể sao chép trực tiếp sang định dạng bảng tính Microsoft Excel hoặc xuất thành mảng dữ liệu ASCII tương thích với phần mềm phân tích chuyên dụng Genie-2000.

5. Chi phí đầu tư hệ thống thực nghiệm và tính khả thi kinh tế?

Việc kết hợp mô phỏng AASI giúp giảm thiểu tần suất vận hành và hao mòn thiết bị buồng chân không của hệ Alpha Analyst, kéo dài tuổi thọ đầu dò bán dẫn PIPS (vốn có chi phí thay thế từ $2,000 - $5,000 USD/đầu dò), nâng cao hiệu quả kinh tế cho phòng phân tích.


Kết luận

Đồ án đã nghiên cứu và làm chủ thành công phần mềm mô phỏng phổ kế alpha AASI trên nền tảng Monte Carlo, đồng thời thực hiện quy trình kiểm chuẩn khoa học nghiêm ngặt với hệ đo Alpha Analyst thực tế tại Đại học Khoa học Tự nhiên TP.HCM. Kết quả cho thấy sự tương thích xuất sắc giữa lý thuyết mô phỏng và thực nghiệm đo phổ của các nguồn chuẩn $^{238}\text{U}, ^{234}\text{U}, ^{239}\text{Pu}, ^{241}\text{Am}$ cũng như mẫu Uranium tự nhiên. Công trình là tài liệu kỹ thuật và giải pháp ứng dụng giá trị cho công tác nghiên cứu, giảng dạy và quan trắc phóng xạ môi trường tại Việt Nam.