Giới thiệu dự án

Công nghệ ghi nhận và phân tích phổ gamma đóng vai trò hạt nhân trong nghiên cứu vật lý hạt nhân, quan trắc phóng xạ môi trường, y học hạt nhân và kiểm định công nghiệp. Trong các hệ phổ kế gamma hiện đại, detector bán dẫn siêu tinh khiết HPGe (High-Purity Germanium) – đặc biệt là dòng detector năng lượng thấp LEGe (Low Energy Germanium) – được đánh giá là tiêu chuẩn vàng nhờ độ phân giải năng lượng vượt trội (lên đến 0,1%), vượt xa các dòng detector nhấp nháy truyền thống như NaI(Tl). Tuy nhiên, việc xác định hàm đáp ứng phổ và hiệu suất ghi nhận thực nghiệm đòi hỏi chi phí chuẩn bị nguồn chuẩn phóng xạ đắt đỏ, tốn thời gian căn chỉnh hình học và tiềm ẩn rủi ro chiếu xạ.

Phương pháp mô phỏng ngẫu nhiên Monte Carlo thông qua bộ công cụ chuẩn quốc tế MCNP4C2 (Monte Carlo N-Particle v4C2) từ Phòng thí nghiệm Quốc gia Los Alamos (LANL) mở ra giải pháp số hóa mạnh mẽ để giải quyết bài toán vận chuyển bức xạ phức tạp. Dự án tập trung nghiên cứu ảnh hưởng của các tham số vật lý chuyên sâu (NOCOH, IDES, MODE P E) và mật độ detector đến hàm đáp ứng phổ gamma, đồng thời đánh giá hiệu năng che chắn phông phóng xạ của ba cấu hình buồng chì chuẩn công nghiệp (737, 747, 777).

                      +------------------------------------------+
                      |       Nguon Buc Xa (Diem / Dia)          |
                      +------------------------------------------+
                                           |
                                           v [Buc xa Gamma]
                      +------------------------------------------+
                      |        Cua so Be (Beryllium 1.5mm)       |
                      +------------------------------------------+
                                           |
                                           v
                      +------------------------------------------+
                      |        Lop chet Ge (Dead Layer 100um)    |
                      +------------------------------------------+
                                           |
                                           v
                      +------------------------------------------+
                      |    Vung hoat dong HPGe (Vung Ngheo P-I-N)|
                      +------------------------------------------+
                                           |
                                           v
                      +------------------------------------------+
                      |  Tally F8 + GEB -> Pho Bien Do Xung Gauss |
                      +------------------------------------------+

Vấn đề kỹ thuật (Problem Statement)

Hiệu chỉnh phổ kế gamma thực nghiệm gặp phải các nút thắt kỹ thuật lớn:

  1. Thiếu khả năng khảo sát vi mô: Thực nghiệm không thể bóc tách tác động độc lập của từng cơ chế tương tác (tán xạ kết hợp Thomson, tán xạ không kết hợp Compton, tạo cặp, bức xạ hãm).
  2. Hiện tượng suy giảm năng lượng ở lớp chết (Dead layer): Lớp bán dẫn không nhạy dày $100,\mu\text{m}$ trên bề mặt tinh thể làm suy giảm nghiêm trọng hiệu suất ở vùng năng lượng thấp ($< 50\text{ keV}$).
  3. Nhiễu phông môi trường và bức xạ huỳnh quang chì: Bức xạ tán xạ và tia X đặc trưng của chì ($75 - 85\text{ keV}$) làm biến dạng vùng phổ quan tâm nếu không có cấu hình che chắn tối ưu.

Mục tiêu dự án

  1. Xây dựng mô hình hình học 3D chi tiết của hệ phổ kế LEGe kết hợp cấu trúc buồng chì trong tệp đầu vào MCNP4C2.
  2. Khảo sát định lượng sai số hiệu suất đỉnh khi kích hoạt/vô hiệu hóa các tham số mô phỏng: NOCOH (tán xạ kết hợp), IDES (vận chuyển electron/bức xạ hãm TTB), và chuyển đổi giữa MODE P (chỉ photon) và MODE P E (kết hợp photon-electron).
  3. Đánh giá sự phụ thuộc của hiệu suất đỉnh ($\varepsilon_p$) vào mật độ tinh thể Germanium trong dải biến thiên từ $5,0\text{ g/cm}^3$ đến $5,5\text{ g/cm}^3$ (chuẩn $5,323\text{ g/cm}^3$).
  4. So sánh khả năng suy giảm phông môi trường của 3 cấu hình buồng chì thương mại: Model 737, Model 747 và Model 777 (Ultra-low background).

Phương pháp tiếp cận và kết quả kỳ vọng

  • Phương pháp: Ứng dụng kỹ thuật truy xuất phân bố độ cao xung (Tally F8) kết hợp hàm mở rộng Gauss (GEB - Gaussian Energy Broadening) với số lượng lịch sử hạt $NPS \ge 10^6$ để đạt sai số thống kê $< 1%$.
  • Kết quả kỳ vọng: Xác lập bộ tham số mô phỏng tối ưu cho hệ phổ kế LEGe, định lượng độ lệch phổ so với benchmark quốc tế (Quados 2003), và đưa ra khuyến nghị lựa chọn buồng chì che chắn cho các phòng thí nghiệm đo phóng xạ năng lượng thấp.

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

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

Trong đo lường bức xạ gamma, việc xác định hàm đáp ứng của detector có thể thực hiện qua ba hướng tiếp cận chính:

Tiêu chí Đo thực nghiệm nguồn chuẩn Mô hình giải tích bán thực nghiệm Mô phỏng Monte Carlo (MCNP)
Độ chính xác Rất cao tại các điểm năng lượng chuẩn Trung bình (sai số lớn ở biên hình học) Rất cao ($R^2 > 0,99$), mô phỏng toàn dải năng lượng
Chi phí vận hành Đắt (mua nguồn phóng xạ, bảo dưỡng LN2) Thấp (chỉ tính toán công thức) Rất thấp (không cần nguồn mẫu thực tế)
Khả năng tùy biến Cố định theo hình học phòng thí nghiệm Hạn chế ở các hình học bất đối xứng Không giới hạn (thay đổi detector, buồng chì tùy ý)
Phân tích vi mô Không thể tách rời cơ chế tán xạ Không hỗ trợ Phân tích chi tiết từng cơ chế tương tác vật lý

Yêu cầu hệ thống theo mô hình MoSCoW

  • Must-have (Bắt buộc): Mô hình hóa chính xác 9 ô vật liệu (cell) và 19 mặt biên (surface); kích hoạt xử lý mở rộng đỉnh Gauss (GEB) theo hàm thực nghiệm; áp dụng Tally F8 cho phổ năng lượng 8192 kênh.
  • Should-have (Nên có): Thiết lập thẻ vật lý PHYS:P với giới hạn năng lượng chi tiết EMCPE = 100\text{ MeV}; khảo sát tương tác photon-electron kết hợp (MODE P E).
  • Could-have (Có thể có): Mô phỏng chi tiết 3 cấu hình buồng chì phức hợp nhiều lớp chì/đồng/thiếc (Pb/Cu/Sn).
  • Won't-have (Chưa thực hiện đợt này): Hiệu ứng giãn nở Doppler trong tinh thể (do giới hạn phiên bản MCNP4C2) và tán xạ hạt nhân photonuclear (ISPN = 0).

Thiết kế hệ thống mô phỏng

Hệ thống phổ kế HPGe LEGe được cấu trúc thành 9 ô phân vùng vật liệu đồng nhất trong không gian pha 3D:

+-------------------------------------------------------------------------+
| [O 9: Khong gian vo han ben ngoai - Importance = 0]                    |
|  +-------------------------------------------------------------------+  |
|  | [O 8: Khong gian khong khi bao quanh he do - Importance = 1]      |  |
|  |   +-------------------------------------------------------------+ |  |
|  |   | [O 1: Nguon buc xa dia r=2.5mm, dat cach cua so Be 2.0 cm]   | |  |
|  |   +-------------------------------------------------------------+ |  |
|  |                                 |                                 |  |
|  |   +-------------------------------------------------------------+ |  |
|  |   | [O 7: Cua so Beryllium (Be) r=26.0mm, day 1.5mm]            | |  |
|  |   | [O 5: Buong chan khong cach tinh the 5.0mm]                 | |  |
|  |   | [O 6: Vo nhom ben ngoai Al day 1.5mm, r=41.5mm, h=58.0mm]   | |  |
|  |   | [O 4: Vanh kep Al giu tinh the (phan sau & phan ben)]        | |  |
|  |   | [O 3: Lop chet Ge Dead Layer day 100 um bao quanh]          | |  |
|  |   | [O 2: Vung hoat dong Ge sieu tinh khiet r=19.9mm, h=14.8mm]  | |  |
|  |   +-------------------------------------------------------------+ |  |
|  +-------------------------------------------------------------------+  |
+-------------------------------------------------------------------------+

Bảng cấu hình công nghệ và thông số vật liệu

  • Mã nguồn tính toán: MCNP4C2 (Monte Carlo N-Particle Transport Code System, Version 4C2).
  • Thư viện tiết diện vi phân: Thư viện mở rộng photon/electron chuẩn LANL (MCPLIB02, EL03).
  • Môi trường vận hành: Hệ thống máy trạm tính toán hiệu năng cao chạy nhân Linux/Unix.
  • Hệ thống buồng chì giảm phông:
    • Cấu hình 737: Vỏ thép carbon $0,95\text{ cm}$, lớp chì $10\text{ cm}$, lớp lót hấp thụ tia X thứ cấp gồm $\text{Cu } 1,5\text{ mm}$ và $\text{Sn } 1,0\text{ mm}$. Thiết kế có cổng ngang kết nối bình nitơ lỏng ($\text{LN}_2$).
    • Cấu hình 747: Vỏ thép $0,95\text{ cm}$, lớp chì $10\text{ cm}$, $\text{Cu } 1,6\text{ mm}$, $\text{Sn } 1,0\text{ mm}$, tối ưu diện tích sàn ($0,4\text{ m}^2$).
    • Cấu hình 777: Cấu hình phông siêu thấp (Ultra-low background), vỏ thép $0,95\text{ cm}$, lớp chì siêu dày $15\text{ cm}$ (tăng 50% bề dày chì so với bản 737), $\text{Cu } 1,5\text{ mm}$, $\text{Sn } 1,0\text{ mm}$.

Implementation và kết quả

Cơ sở động lực học và thuật toán tương tác photon

MCNP4C2 mô hình hóa 4 quá trình tương tác cơ bản của photon với vật chất:

1. Tán xạ Compton (Incoherent Scattering)

Tiết diện tán xạ vi phân được tính thông qua biểu thức Klein-Nishina kết hợp thừa số hiệu chỉnh dạng tán xạ không kết hợp $I(Z, v)$: $$\frac{d\sigma_c(Z, \alpha)}{d\mu} = \pi r_0^2 \left( \frac{\alpha'}{\alpha} \right)^2 \left[ \frac{\alpha'}{\alpha} + \frac{\alpha}{\alpha'} + \mu^2 - 1 \right] I(Z, v)$$ Trong đó: $r_0 = 2,8179 \times 10^{-13}\text{ cm}$ là bán kính cổ điển của electron; $\alpha = \frac{E}{m_e c^2}$ là năng lượng photon tới theo đơn vị khối lượng nghỉ electron; $\mu = \cos\theta$; $I(Z, v)$ suy giảm tiết diện tán xạ về phía trước ở năng lượng thấp đối với vật liệu $Z$ cao.

2. Tán xạ Thomson (Coherent Scattering)

Tán xạ kết hợp không làm thay đổi năng lượng photon mà chỉ đổi hướng bay. Tiết diện tán xạ vi phân được hiệu chỉnh qua hàm $C^(Z, v)$: $$\sigma_T(Z, \alpha, \mu) d\mu = C^(Z, v) \pi r_0^2 (1 + \mu^2) d\mu$$ Khi thiết lập NOCOH = 1, thuật toán MCNP bỏ qua tán xạ kết hợp, gây ảnh hưởng trực tiếp đến phân bố phổ ở vùng năng lượng thấp ($< 50\text{ keV}$).

3. Hiệu ứng quang điện (Photoelectric Effect)

Năng lượng photon tới $E$ bị hấp thụ hoàn toàn bởi electron liên kết có năng lượng $e < E$, giải phóng quang electron với động năng $E_k = E - e$. Lỗ trống điện tử được lấp đầy kèm theo phát xạ tia X huỳnh quang đặc trưng ($K_\alpha, K_\beta$) hoặc electron Auger, được mô hình hóa chi tiết qua chế độ MODE P E hoặc xấp xỉ bức xạ hãm bia dày TTB (Thick-Target Bremsstrahlung).

4. Hiệu ứng tạo cặp (Pair Production)

Xảy ra trong trường Coulomb hạt nhân khi $E_\gamma \ge 2 m_e c^2 = 1,022\text{ MeV}$. Cặp electron-positron ($\beta^-, \beta^+$) sinh ra; positron hủy hạt tạo ra 2 photon $0,511\text{ MeV}$ bay ngược chiều góc $180^\circ$.

5. Mở rộng phổ thực nghiệm Gauss (GEB Card)

Để chuyển đổi phân bố xung lý thuyết thành phổ Gauss thực tế, độ mở rộng xung $\Delta E$ (FWHM) được khớp thực nghiệm theo hàm: $$\text{FWHM}(E) = a + b\sqrt{E} + cE^2 = 3,592\times 10^{-4} + 1,832\times 10^{-2}\sqrt{E} + 0,01918E$$ Hệ số bán độ rộng Gauss $A$ liên hệ với FWHM qua công thức: $$A = \frac{\text{FWHM}}{2\sqrt{\ln 2}} \approx 0,600561 \times \text{FWHM}$$

c =================================================================
c FILE DAU VAO MCNP4C2: MO PHONG PHO KE DETECTOR HPGe LEGe
c =================================================================
c 1. KHOI DINH NGHIA O MANG (CELL CARDS)
1    1  -0.0012   -1             imp:p,e=1  $ O nguon buc xa dia
2    2  -5.3230    1 -2 -3       imp:p,e=1  $ Vung hoat dong Ge
3    2  -5.3230    2 -4 -3       imp:p,e=1  $ Lop chet Dead Layer 100um
4    3  -2.7000    4 -5 -6       imp:p,e=1  $ Vanh kep nhom Al
5    0             5 -7 -8       imp:p,e=1  $ Buong chan khong
6    3  -2.7000    7 -9 -10      imp:p,e=1  $ Vo boc ngoai Al
7    4  -1.8500    9 -11 -12     imp:p,e=1  $ Cua so Beryllium Be
8    1  -0.0012   11 -13 -14     imp:p,e=1  $ Khong khi bao quanh
9    0            13 : 14        imp:p,e=0  $ Khong gian ngoai (Kill cell)

c 2. KHOI THE NANG LUONG VA VAT LY (DATA CARDS)
MODE P E
PHYS:P 100 0 0 $ EMCPE=100MeV, IDES=0, NOCOH=0
SDEF POS=0 0 2.0 AXS=0 0 1 RAD=d1 ERG=0.030 PAR=2
SI1  0 0.25
SP1 -21 1
F8:P 2
FT8  GEB 3.592E-4 1.832E-2 0.01918
E8   0.001 8191i 0.040
NPS  2000000

Kết quả đo lường và kiểm chứng

1. Ảnh hưởng của tham số NOCOH và IDES (Mật độ detector chuẩn $5,323\text{ g/cm}^3$)

Do Sai Biet Hieu Suat Dinh (%)
  8.0% |                                      [x] MODE P E vs P (15 keV: 7.30%)
  7.0% |
  6.0% |
  5.0% |  [*] NOCOH=1 (15 keV: 4.69%)
  4.0% |
  3.0% |                                            [#] IDES=1 (1 MeV: 3.02%)
  2.0% |
  1.0% |        [*] NOCOH=1 (30 keV: 0.61%)
  0.0% +--------------------------------------------------------------------
             15 keV            30 keV            750 keV           1 MeV
Mức năng lượng ($E_\gamma$) TH 00 (Mặc định) $\varepsilon_p$ TH 01 (NOCOH=1) Sai biệt (%) TH 10 (IDES=1) Sai biệt (%) TH 11 (NOCOH=1, IDES=1) Sai biệt (%) So sánh Quados 2003 (NOCOH=1)
15 keV $4,862 \times 10^{-3}$ 4,69% 0,00% 4,69% 3,40% (Lệch 37,94%)
30 keV $5,871 \times 10^{-3}$ 0,61% 0,00% 0,61% 0,70% (Lệch 12,86%)
100 keV $4,320 \times 10^{-3}$ 0,05% 0,12% 0,16% 0,10% (Tương thích tốt)
750 keV $8,412 \times 10^{-4}$ 0,00% 2,52% 2,52% 1,50% (Lệch 68,00%)
1000 keV (1 MeV) $6,215 \times 10^{-4}$ 0,00% 3,02% 3,02% 2,20% (Lệch 37,27%)

Phân tích: Ở năng lượng thấp ($15 - 30\text{ keV}$), việc tắt tán xạ kết hợp (NOCOH = 1) làm sai lệch hiệu suất đỉnh tới $4,69%$ do tán xạ Thomson bẻ cong hướng photon đi vào vùng nhạy. Ở vùng năng lượng cao ($750 - 1000\text{ keV}$), tham số IDES = 1 (tắt sinh electron thứ cấp) chiếm ưu thế gây sai lệch $3,02%$ do không tính đến sự thất thoát electron động năng lớn ra khỏi thể tích nhạy.

2. Sai biệt giữa MODE P E và MODE P

  • Tại 15 keV và 30 keV, MODE P E cho hiệu suất đỉnh cao hơn MODE P lần lượt là 7,30% và 1,52% (so sánh với Quados 2003 đạt $5,50%$ và $1,40%$, mức sai khác tương ứng $32,73%$ và $8,57%$).
  • Tại 750 keV và 1 MeV, hiệu suất đỉnh của MODE P E thấp hơn MODE P lần lượt là -2,52% và -3,02% do hiện tượng rò rỉ electron thứ cấp khỏi thể tích nhạy của tinh thể.

3. Tương quan giữa mật độ tinh thể và hiệu suất đỉnh

Khảo sát biến thiên mật độ Ge từ $5,000\text{ g/cm}^3$ đến $5,500\text{ g/cm}^3$:

  • Tại dải năng lượng trung bình và cao ($> 50\text{ keV}$): Hiệu suất đỉnh $\varepsilon_p$ tỷ lệ thuận tuyến tính với mật độ detector (mật độ hạt nhân tương tác trong thể tích nhạy tăng).
  • Tại dải năng lượng thấp ($15\text{ keV}$ và $30\text{ keV}$): Hiệu suất đỉnh tỷ lệ nghịch với mật độ detector. Nguyên nhân vật lý: Mật độ tinh thể tăng kéo theo mật độ khối của lớp chết ($100,\mu\text{m}$) tăng, làm tăng hấp thụ quang điện ngay tại lớp chết trước khi photon tiếp cận được vùng hoạt động $I$.

4. Đánh giá định lượng hiệu quả buồng chì (Độ lệch phổ giảm phông)

Cấu hình buồng chì Bề dày lớp chì (Pb) Lớp lót chống tia X chì Độ lệch phổ tại 1 MeV ($\rho = 5,0\text{ g/cm}^3$) Khả năng giảm phông môi trường
Model 747 $10\text{ cm}$ $\text{Cu } 1,6\text{ mm} + \text{Sn } 1,0\text{ mm}$ 0,7799% Khá (Tối ưu không gian hẹp)
Model 737 $10\text{ cm}$ $\text{Cu } 1,5\text{ mm} + \text{Sn } 1,0\text{ mm}$ 0,8230% Tốt (Có cổng dewar $\text{LN}_2$)
Model 777 15 cm (+50%) $\text{Cu } 1,5\text{ mm} + \text{Sn } 1,0\text{ mm}$ 1,0814% Xuất sắc (Ultra-low background)

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

  1. Chuẩn hóa quy trình mô phỏng vi mô cho detector LEGe: Xây dựng thành công tệp cấu hình tham số MCNP4C2 chi tiết cho dòng detector LEGe với đầy đủ lớp bọc Al, cửa sổ Be, chân không và lớp chết Ge.
  2. Làm rõ cơ chế đối nghịch của lớp chết ở năng lượng thấp: Chứng minh định lượng hiện tượng nghịch đảo hiệu suất theo mật độ tinh thể tại vùng $15 - 30\text{ keV}$, giải thích cơ chế suy giảm quang điện trong lớp chết $100,\mu\text{m}$.
  3. Định lượng tác động của các cờ vật lý (Physics Cards): Xác định ngưỡng năng lượng mà tại đó các tham số NOCOH và IDES có thể bỏ qua (dưới 5%) để tối ưu hóa thời gian tính toán máy tính mà vẫn duy trì độ tin cậy khoa học.
  4. Đánh giá hiệu năng che chắn của 3 dòng buồng chì công nghiệp: Cung cấp dữ liệu định lượng về mức độ suy giảm phông của buồng chì Model 777 (với bề dày chì $15\text{ cm}$) vượt trội $31,4%$ so với Model 737 tại năng lượng $1\text{ MeV}$.

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

Trường hợp ứng dụng (Use Cases)

  • Quan trắc phóng xạ môi trường biển & đất: Xác định hàm lượng các đồng vị phát gamma năng lượng thấp như $^{210}\text{Pb}$ ($46,5\text{ keV}$), $^{241}\text{Am}$ ($59,5\text{ keV}$) với giới hạn phát hiện thấp (LOD) tối ưu.
  • Kiểm tra không phá hủy (NDT) & Giám định hạt nhân: Phân tích đồng vị Uranium/Plutonium trong mẫu nhiên liệu hạt nhân sử dụng buồng chì Model 777 siêu phông thấp.
  • Hiệu chuẩn ảo cho phòng thí nghiệm hạt nhân: Giảm thiểu 80% thời gian thực nghiệm căn chỉnh đường cong hiệu suất nguồn chuẩn đắt tiền.
+--------------------------------------------------------------------------+
|                       LO TRINH TRIEN KHAI 4 GIAI DOAN                    |
+--------------------------------------------------------------------------+
|  Giai doan 1 (Tuan 1 - 2): Chuan hoa hinh hoc detector & Buong chi      |
|  Giai doan 2 (Tuan 3 - 5): Chay mo phong Monte Carlo MCNP4C2 (NPS=2M)    |
|  Giai doan 3 (Tuan 6 - 7): Tich hop ham GEB & Chuan hoa du lieu thuc do |
|  Giai doan 4 (Tuan 8):     Xuat duong cong hieu suat dinh & Trien khai   |
+--------------------------------------------------------------------------+

Phân tích hiệu quả kinh tế (Cost-Benefit Analysis)

  • Chi phí thực nghiệm truyền thống: Bộ nguồn chuẩn đa đỉnh chuẩn hóa quốc tế ($\approx $15.000 - $25.000$) kèm chi phí suy hao phóng xạ định kỳ 2 năm/lần.
  • Giải pháp mô phỏng Monte Carlo MCNP: Chi phí đầu tư máy trạm tính toán một lần ($\approx $3.000$), tiết kiệm hơn 85% ngân sách hiệu chuẩn ban đầu, thời gian xác lập đường cong hiệu suất rút ngắn từ 3 tuần xuống còn 48 giờ tính toán.

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

Hạn chế kỹ thuật

  • Giới hạn phiên bản MCNP4C2: Chưa hỗ trợ tính năng làm rộng phổ do chuyển động nhiệt Doppler của electron liên kết (NODOP = 0), ảnh hưởng nhẹ đến hình dạng đáy pic tán xạ Compton.
  • Sai biệt với dữ liệu tham chiếu Quados (2003): Độ lệch hiệu suất đỉnh ở mức $750\text{ keV}$ và $1\text{ MeV}$ giữa tính toán mô phỏng và Quados còn cao ($37,27% - 68%$), cần hiệu chỉnh thêm về mô hình hình học chi tiết của nguồn phóng xạ.

Định hướng mở rộng

  1. Nâng cấp mô hình lên phiên bản MCNP5 / MCNP6 hoặc GEANT4 để tích hợp hiệu ứng Doppler đầy đủ và khảo sát tương tác photonuclear.
  2. Khảo sát hàm đáp ứng với các nguồn dạng hình học thể tích phức tạp (mẫu đất Marinelli beaker, mẫu khí dạng trụ).
  3. Ứng dụng thuật toán tối ưu hóa bề dày lớp chết Ge thông qua so sánh phổ mô phỏng với phổ thực nghiệm từ nguồn chuẩn $^{133}\text{Ba}$, $^{152}\text{Eu}$.

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

+-----------------------------------------------------------------------------+
| Sinh vien & Hoc vien cao hoc: Tiep can giao trinh so & Bo input card chuan  |
| Ky su & Chuyen gia do luong:  Toi uu cau hinh buong chi & Thoi gian do      |
| Phong thi nghiem hat nhan:    Tiet kiem 85% chi phi mua nguon chuan quoc te |
| Nha nghien cuu vat ly ly thuyet: Co so du lieu benchmark tuong tac vi mo   |
+-----------------------------------------------------------------------------+
  • Sinh viên & Học viên chuyên ngành Vật lý hạt nhân: Nắm vững phương pháp xây dựng input card MCNP, kỹ thuật chia cell, tally F8 và hàm khớp GEB chuẩn mực.
  • Kỹ sư vận hành hệ phổ kế: Lựa chọn chính xác cấu hình buồng chì (737/747/777) phù hợp với diện tích phòng thí nghiệm và yêu cầu giới hạn phông nền.
  • Các viện nghiên cứu & Trung tâm quan trắc: Giảm thiểu chi phí mua sắm nguồn chuẩn phóng xạ, chuẩn hóa đường cong hiệu suất tự động cho các mẫu đo môi trường.

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

1. Cần cấu hình phần cứng tối thiểu nào để chạy mô phỏng MCNP4C2 cho hệ phổ kế?

Hệ thống cần máy trạm CPU tối thiểu 4 nhân (khuyến nghị 8 nhân trở lên), 8GB RAM, hệ điều hành Linux (Ubuntu/CentOS) hoặc Windows có trình biên dịch Fortran/C. Thời gian chạy 1 triệu lịch sử hạt ($NPS = 10^6$) mất khoảng 15 - 30 phút tùy thuộc vào chế độ MODE P hay MODE P E.

2. Khi nào có thể bỏ qua tham số NOCOH và IDES trong mô phỏng?

Đối với các bài toán đo bức xạ gamma thông thường có mức năng lượng trong khoảng $100\text{ keV} - 600\text{ keV}$, ảnh hưởng của NOCOH và IDES đến hiệu suất đỉnh đều nhỏ hơn $0,2%$. Do đó, có thể bỏ qua để tăng tốc độ tính toán mà không làm suy giảm độ chính xác của phổ.

3. Tại sao cấu hình buồng chì 777 lại cho khả năng giảm phông tốt nhất?

Cấu hình 777 sở hữu lớp chì nguyên chất dày $15\text{ cm}$ (dày hơn 50% so với cấu hình 737 và 747 có lớp chì dày $10\text{ cm}$). Bề dày tăng thêm này làm tăng xác suất hấp thụ quang điện và tán xạ Compton đối với bức xạ gamma phông năng lượng cao từ môi trường (như $^{40}\text{K}$, $^{208}\text{Tl}$).

4. Lớp lót Cu/Sn bên trong buồng chì có vai trò gì?

Khi tia gamma tương tác với thành chì sẽ phát ra tia X huỳnh quang chì đặc trưng ($75 - 85\text{ keV}$). Lớp thiếc ($\text{Sn } 1,0\text{ mm}$) hấp thụ tia X của chì và phát ra tia X thiếc ($25\text{ keV}$). Tiếp theo, lớp đồng ($\text{Cu } 1,5\text{ mm}$) hấp thụ tia X của thiếc và chỉ phát ra tia X đồng năng lượng rất thấp ($8\text{ keV}$), mức năng lượng này bị hấp thụ hoàn toàn bởi không khí và cửa sổ Beryllium trước khi tới detector.

5. Sự khác biệt căn bản giữa MODE P và MODE P E là gì?

MODE P chỉ vận chuyển hạt photon; electron sinh ra từ tương tác quang điện/Compton được giả định dừng lại và truyền toàn bộ năng lượng tại điểm tương tác (hoặc dùng xấp xỉ TTB). MODE P E vận chuyển chi tiết cả quỹ đạo của electron và positron, cho độ chính xác cao nhất nhưng thời gian tính toán tăng từ 3 đến 5 lần.

Kết luận

Nghiên cứu đã ứng dụng thành công chương trình mô phỏng Monte Carlo MCNP4C2 để phân tích chuyên sâu các yếu tố ảnh hưởng đến hàm đáp ứng phổ gamma của detector bán dẫn siêu tinh khiết HPGe LEGe. Kết quả chỉ rõ tán xạ kết hợp (NOCOH) chi phối chủ yếu ở vùng năng lượng thấp ($< 50\text{ keV}$), trong khi cơ chế vận chuyển electron (IDES) đóng vai trò quyết định ở vùng năng lượng cao ($> 700\text{ keV}$). Hiện tượng nghịch đảo hiệu suất theo mật độ tinh thể tại dải năng lượng thấp đã làm sáng tỏ cơ chế suy giảm quang điện tại lớp chết $100,\mu\text{m}$. Đồng thời, việc đánh giá thực nghiệm mô phỏng trên 3 cấu hình buồng chì khẳng định Model 777 là lựa chọn tối ưu cho các phép đo phóng xạ mức siêu thấp. Công trình tạo tiền đề vững chắc cho việc thiết lập hệ hiệu chuẩn phổ kế ảo độ chính xác cao trong nghiên cứu hạt nhân ứng dụng.