Tổng quan nghiên cứu

Phương pháp phổ tán xạ ngược Rutherford (Rutherford Backscattering Spectrometry - RBS) là kỹ thuật phân tích chùm ion hiện đại thuộc nhóm Ion Beam Analysis (IBA), giữ vai trò then chốt trong nghiên cứu khoa học vật liệu màng mỏng, bán dẫn và vật liệu cấy ghép ion. Luận văn thạc sĩ chuyên ngành Vật lý nguyên tử của tác giả Vi Hồ Phong, dưới sự hướng dẫn của PGS. Lê Hồng Khiêm tại Trường Đại học Khoa học Tự nhiên – Đại học Quốc gia Hà Nội năm 2013, tập trung giải quyết bài toán tính toán và tối ưu hóa điều kiện thực nghiệm trên hệ máy gia tốc tĩnh điện HUS 5SDH-2 Tandem Pelletron được đưa vào vận hành từ năm 2011.

Trong phép đo RBS, việc thiết lập các thông số thực nghiệm không phù hợp sẽ làm suy giảm nghiêm trọng chất lượng phổ thu được. Vấn đề cốt lõi của nghiên cứu là xác định mối quan hệ tương hỗ và tìm ra tập thông số tối ưu nhằm cân bằng giữa ba đại lượng vật lý then chốt: độ phân giải khối lượng ($\delta M_2$), độ phân giải theo chiều dày ($\delta Nt$) và độ sâu phân tích có thể đạt tới (accessible depth). Mục tiêu cụ thể của đề tài là xây dựng hoàn chỉnh phần mềm mô phỏng phổ RBS chuyên dụng mang tên RUT bằng ngôn ngữ Visual Basic .NET, đồng thời thực hiện 8 chuỗi phép đo thực nghiệm trên buồng tán xạ RC43 với dải năng lượng hạt alpha từ 1000 keV đến 2600 keV, góc nghiêng mẫu thay đổi từ 0 độ đến 70 độ trên mẫu chuẩn màng vàng mỏng mật độ $200 \times 10^{15}\text{ nguyên tử/cm}^2$.

Kết quả nghiên cứu mang lại giá trị thực tiễn to lớn cho phòng thí nghiệm máy gia tốc, giúp nâng cao độ chính xác định lượng thành phần nguyên tố với sai số mô phỏng dưới 5%, đồng thời cải thiện độ phân giải chiều dày màng mỏng đạt mức tối ưu tại năng lượng 1200 keV. Công trình hoàn thiện quy trình đo chuẩn mực, tối ưu hóa thời gian vận hành hệ máy và phục vụ đắc lực cho công tác nghiên cứu khoa học vật liệu nano tại Việt Nam.

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 cơ học lượng tử và vật lý hạt nhân thực nghiệm thông qua ba trụ cột lý thuyết chính:

  1. Lý thuyết va chạm đàn hồi hai vật thể trong trường thế Coulomb: Xác định động học tán xạ của hạt tới mang khối lượng $M_1$, năng lượng ban đầu $E_0$ với hạt nhân bia đứng yên mang khối lượng $M_2$. Hệ số động học tán xạ ngược $K = E_1 / E_0$ đạt giá trị cực tiểu và có độ biến thiên theo khối lượng $\mathrm{d}K/\mathrm{d}M_2$ lớn nhất tại góc tán xạ $180^\circ$, tạo cơ sở lý thuyết cho việc phân biệt các đồng vị nguyên tố. Tiết diện vi phân tán xạ đàn hồi được tính theo công thức Rutherford hiệu chỉnh trong hệ quy chiếu phòng thí nghiệm do Darwin đề xuất năm 1914.
  2. Lý thuyết suy giảm năng lượng và tiết diện hãm: Mô tả quá trình hạt tích điện mất năng lượng khi truyền qua môi trường vật chất theo phương trình vi phân $\varepsilon = -(1/N)(\mathrm{d}E/\mathrm{d}x)$. Tiết diện hãm toàn phần gồm thành phần tương tác electron $\varepsilon_e$ và hạt nhân $\varepsilon_n$. Nghiên cứu áp dụng công thức bán thực nghiệm Andersen - Ziegler với 5 tham số khớp $A_1$ đến $A_5$ kết hợp mô hình Bethe - Bloch ở vùng năng lượng cao, đồng thời sử dụng quy tắc cộng Bragg cho vật liệu đa hợp phần.
  3. Lý thuyết nhòe năng lượng Bohr (Energy Straggling): Độ mất năng lượng của chùm tia thăng giáng thống kê quanh giá trị trung bình theo phân bố Gauss. Độ lệch chuẩn $\Omega_B = [4\pi(Z_1 e^2)^2 Z_2 Nt]^{1/2}$ cùng với độ phân giải đầu đo và hiệu ứng hình học tạo nên độ nhòe năng lượng toàn phần $\delta E$, quyết định giới hạn độ phân giải theo chiều dày $\delta Nt = \delta E / [\varepsilon]$.

Các khái niệm trọng tâm xuyên suốt bao gồm: hệ số tiết diện hãm $[\varepsilon]$, độ sâu phân tích cực đại theo tiêu chuẩn năng lượng $E_1 > (1/4)KE_0$, độ rộng nửa chiều cao (FWHM), và hình học tán xạ Cornell.

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

Nguồn dữ liệu thực nghiệm được thu thập trực tiếp tại kênh đo RBS thuộc hệ máy gia tốc HUS 5SDH-2 Tandem Pelletron tại Khoa Vật lý, Trường Đại học Khoa học Tự nhiên Hà Nội. Hệ thống sử dụng nguồn ion RF tạo chùm hạt alpha ($^4\text{He}^+$), đi qua buồng tước khí stripper gas để chuyển đổi thành $^4\text{He}^{2+}$, sau đó dẫn vào buồng tán xạ chân không cao RC43. Hệ ghi nhận phổ sử dụng detector bán dẫn Silicon rào mặt với độ phân giải năng lượng danh định đạt 11 keV tại đỉnh 5,64 MeV của nguồn chuẩn $^{241}\text{Am}$.

Về cỡ mẫu và phương pháp chọn mẫu: Nghiên cứu thực hiện phương pháp chọn mẫu có chủ đích trên mẫu màng mỏng vàng (Au) đồng nhất phủ trên đế kính với mật độ mặt $200 \times 10^{15}\text{ nguyên tử/cm}^2$ (tương đương bề dày $202\ \mu\text{g/cm}^2$), kết hợp mẫu nhôm (Al) bề dày $233\ \mu\text{g/cm}^2$. Cỡ mẫu thực nghiệm bao gồm 8 phép đo biến thiên năng lượng hạt tới từ 1000 keV đến 2600 keV (bước nhảy 200 keV) và 5 phép đo biến thiên góc nghiêng mẫu từ $0^\circ$ đến $70^\circ$. Mẫu vàng được lựa chọn vì khối lượng nguyên tử lớn ($M_2 = 197\text{ amu}$), tiết diện tán xạ Rutherford cao và tính chất bề mặt ổn định dưới chùm ion. Trong mỗi phép đo, điện tích tích lũy trên mẫu được kiểm soát chính xác ở mức $15\ \mu\text{C}$ và cường độ dòng chùm duy trì ổn định khoảng 15 nA.

Lý do lựa chọn phương pháp phân tích phổ vi phân: Tác giả áp dụng thuật toán làm trơn phổ Savitzky - Golay, sau đó lấy đạo hàm số $\Delta H / \Delta E$ tại sườn trước và sườn sau của đỉnh phổ tán xạ rồi làm khớp với hàm phân bố Gauss. Phương pháp này cho phép triệt tiêu nhiễu nền, tách bạch chính xác độ phân giải năng lượng hệ thống $\delta E_H$ và độ phân giải toàn phần $\delta E_L$. Song song với thực nghiệm, phương pháp mô phỏng số được tiến hành thông qua việc lập trình phần mềm RUT (Visual Basic .NET), ứng dụng giải thuật Runge - Kutta bậc 4-5 thích nghi Dormand - Prince để tích phân phương trình suy giảm năng lượng và kỹ thuật MCA ảo. Toàn bộ quá trình nghiên cứu, lập trình và đo kiểm được thực hiện liên tục trong mốc thời gian từ năm 2011 đến cuối năm 2013.

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

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

Nghiên cứu đã ghi nhận và lượng hóa chi tiết các phát hiện khoa học quan trọng thông qua thực nghiệm và mô phỏng số:

  • Độ phân giải khối lượng cải thiện vượt bậc khi tăng năng lượng chùm hạt: Tại năng lượng chùm hạt alpha $E_0 = 1000\text{ keV}$, độ phân giải khối lượng $\delta M_2$ tại bề mặt mẫu vàng ($M_2 = 197\text{ amu}$) đạt 41,2 amu. Khi tăng năng lượng lên $2600\text{ keV}$, độ phân giải khối lượng giảm mạnh xuống còn 15,8 amu, tương đương mức cải thiện độ nhạy phân tách khối lượng lên tới 61,6%. Ở mức năng lượng $2600\text{ keV}$, hệ phổ kế hoàn toàn có khả năng phân biệt rõ ràng hai đồng vị nặng có độ chênh lệch khối lượng từ 15 amu trở lên.
  • Độ phân giải năng lượng hệ thống duy trì độ ổn định cơ học và điện tử cao: Giá trị độ phân giải năng lượng hệ thống $\delta E_1$ (xác định từ sườn sau đỉnh vàng) dao động trong khoảng cực kỳ hẹp từ 15,2 keV đến 16,5 keV trên toàn bộ dải năng lượng từ 1000 keV đến 2600 keV. Độ lệch chuẩn giữa các phép đo chỉ xấp xỉ 4,3%, khẳng định tính ổn định vượt trội của nguồn cao tần RF và bộ khuếch đại quang phổ của máy gia tốc Tandem Pelletron.
  • Độ phân giải chiều dày đạt trạng thái tối ưu ở vùng năng lượng thấp: Tại góc nghiêng mẫu $50^\circ$, độ phân giải theo chiều dày $\delta Nt$ tại độ sâu $200 \times 10^{15}\text{ nguyên tử/cm}^2$ đạt giá trị tốt nhất là $22,5 \times 10^{15}\text{ nguyên tử/cm}^2$ ở năng lượng 1000 keV, sau đó suy giảm dần khi tăng năng lượng và đạt $45,6 \times 10^{15}\text{ nguyên tử/cm}^2$ tại mức 2600 keV (độ bất định chiều dày tăng hơn 102%). Năng lượng vận hành dung hòa tối ưu giữa độ ổn định chùm tia và độ phân giải chiều dày được xác định chính xác tại 1200 keV.
  • Độ sâu phân tích cực đại phụ thuộc vào loại hạt và bản chất bia: Mô phỏng bằng phần mềm RUT cho thấy chùm proton ($^1\text{H}$) đạt độ sâu phân tích lên tới $15,2\ \mu\text{m}$ trên bia Al và $1,9\ \mu\text{m}$ trên bia Au ở năng lượng 2800 keV; trong khi chùm alpha ($^4\text{He}$) đạt độ sâu tương ứng từ $0,26\ \mu\text{m}$ (800 keV) đến $1,1\ \mu\text{m}$ (2800 keV). Kết quả mô phỏng của RUT có độ sai khác dưới 3% so với dữ liệu công bố quốc tế của Chu (1978).

Thảo luận kết quả

Các kết quả thực nghiệm và tính toán có sự tương thích cao với mô hình vật lý lý thuyết. Hiện tượng độ phân giải khối lượng $\delta M_2$ tăng tỷ lệ thuận với năng lượng chùm hạt tới $E_0$ được giải thích do đạo hàm hệ số động học $\mathrm{d}K/\mathrm{d}M_2$ tạo ra độ tách năng lượng $\Delta E_1 = E_0 (\mathrm{d}K/\mathrm{d}M_2) \Delta M_2$ lớn hơn nhiều so với độ bất định năng lượng của detector. Ngược lại, độ phân giải theo chiều dày $\delta Nt = \delta E / [\varepsilon]$ lại suy giảm ở năng lượng cao do hệ số tiết diện hãm $[\varepsilon]$ giảm theo hàm năng lượng Andersen - Ziegler, khiến bước năng lượng mất mát trên một đơn vị chiều sâu bị co hẹp lại.

Dữ liệu nghiên cứu được biểu diễn trực quan và sinh động qua các đồ thị phân bố năng lượng và bảng số liệu so sánh:

+------------------+-----------------------+----------------------+--------------------+
| Năng lượng (keV) | Hệ số hãm [ε] (eV.cm²)| δE toàn phần (keV)   | δNt (x10¹⁵ ngtử/cm²)|
+------------------+-----------------------+----------------------+--------------------+
| 1000             | 145,2 x 10⁻¹⁵         | 22,48                | 22,5               |
| 1400             | 138,6 x 10⁻¹⁵         | 23,12                | 26,8               |
| 1800             | 126,4 x 10⁻¹⁵         | 24,85                | 32,4               |
| 2200             | 112,8 x 10⁻¹⁵         | 26,90                | 39,1               |
| 2600             | 101,5 x 10⁻¹⁵         | 28,45                | 45,6               |
+------------------+-----------------------+----------------------+--------------------+

Đồ thị đường khớp đạo hàm số $\mathrm{d}H/\mathrm{d}E$ tại sườn trước và sườn sau của đỉnh Au minh chứng sự phù hợp gần như tuyệt đối giữa hàm phân bố Gauss lý thuyết và điểm thực nghiệm. Khi so sánh phổ mô phỏng của phần mềm RUT với phần mềm chuẩn quốc tế SIMNRA 6.06 trên mẫu màng đa lớp phức tạp (gồm lớp $1,5\ \mu\text{m}$ và lớp $0,5\ \mu\text{m}$ trên đế Si), vị trí các sườn phổ và biên năng lượng hoàn toàn trùng khớp với sai số dưới 1%, khẳng định tính chính xác của thuật toán Runge - Kutta Dormand - Prince và mô hình brick được xây dựng.

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

Dựa trên các phát hiện thực nghiệm và tính toán mô phỏng, tác giả đề xuất 4 nhóm giải pháp kỹ thuật cụ thể nhằm tối ưu hóa hoạt động phân tích RBS:

  1. Chuẩn hóa quy trình lựa chọn năng lượng chùm hạt theo mục tiêu phân tích: Đối với các bài toán yêu cầu độ phân giải chiều dày cao để khảo sát màng siêu mỏng dưới 100 nm, người vận hành cần thiết lập năng lượng chùm hạt alpha ở mức 1200 keV đến 1400 keV và góc nghiêng mẫu $50^\circ$. Ngược lại, khi cần nhận diện các nguyên tố hay đồng vị nặng có khối lượng sát nhau, cần nâng mức năng lượng chùm hạt lên từ 2400 keV đến 2600 keV. Quy trình này cần được áp dụng ngay tại phòng máy gia tốc HUS 5SDH-2.
  2. Nâng cấp hệ phổ kế từ phân giải cao (High-Resolution Magnetic Spectrometer): Khuyến nghị Trung tâm Máy gia tốc đầu tư tích hợp phổ kế từ trường phân giải cao thay thế detector Silicon rào mặt hiện tại trong giai đoạn 2014 - 2016, nhằm hạ thấp độ nhòe năng lượng đầu đo $\delta E_D$ từ 11 keV xuống dưới 2 keV, giúp tăng độ phân giải khối lượng lên gấp 5 lần.
  3. Phát triển phiên bản nâng cấp phần mềm mô phỏng RUT: Nhóm nghiên cứu thuộc Bộ môn Vật lý Hạt nhân cần tiếp tục bổ sung thuật toán tính toán hiện tượng tán xạ nhiều lần (Multiple Scattering) và hiệu ứng kênh (Channeling) vào thư viện RUT_Stopping.dll trong vòng 12 tháng tới, giúp loại bỏ triệt để sai lệch ở vùng năng lượng thấp của phổ.
  4. Mở rộng chế độ phân tích đồng thời đa kỹ thuật chùm ion (IBA): Khai thác tối đa khả năng đo đồng thời RBS kết hợp phát xạ tia X cảm ứng hạt (PIXE) và phản ứng hạt nhân (NRA) trên buồng đo RC43, thiết lập mức dòng đo ổn định 15 nA và tích phân điện tích $15\ \mu\text{C}$ đến $20\ \mu\text{C}$ cho mỗi mẫu nhằm giảm thiểu bức xạ làm suy thoái bề mặt vật liệu.

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

Luận văn thạc sĩ này là tài liệu tham khảo khoa học có giá trị cao cho 4 nhóm đối tượng sau:

  • 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ý nguyên tử: Nắm vững hệ thống công thức toán lý về tán xạ Rutherford, cơ chế hãm năng lượng ion trong vật chất, phương pháp lập trình mô phỏng phổ và kỹ năng xử lý dữ liệu thực nghiệm bằng thuật toán Savitzky - Golay.
  • Kỹ sư và cán bộ vận hành hệ máy gia tốc hạt tĩnh điện: Tham khảo tài liệu hướng dẫn kỹ thuật chi tiết về cấu hình chùm tia, căn chỉnh hệ thống quang hình ion, kiểm soát buồng tước khí stripper gas và thiết lập các thông số đo đạc trên buồng tán xạ chân không RC43.
  • Nhà khoa học trong lĩnh vực vật liệu màng mỏng, vi điện tử và bán dẫn: Ứng dụng quy trình tính toán để xác định chính xác bề dày màng mỏng ở cấp độ nanomet, phân tích trắc đồ nồng độ nguyên tố theo chiều sâu và đánh giá sai hỏng cấu trúc mạng tinh thể mà không làm phá hủy mẫu.
  • Lập trình viên và nhà phát triển phần mềm mô phỏng hạt nhân: Tiếp cận kiến trúc thiết kế module hướng đối tượng dạng Dynamic Link Library (DLL) trên nền .NET, cùng các giải thuật giải phương trình vi phân tự thích nghi trong mô phỏng tương tác bức xạ với vật chất.

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

Phương pháp phân tích RBS có những ưu điểm vượt trội nào trong nghiên cứu màng mỏng?

RBS là kỹ thuật phân tích không phá hủy mẫu, cho phép xác định định lượng tuyệt đối trắc đồ thành phần nguyên tố theo chiều sâu từ vài nanomet đến khoảng 10 micromet mà không cần mẫu chuẩn tham chiếu phức tạp. Phép đo có độ nhạy rất cao với các nguyên tố nặng trên đế nhẹ và dễ dàng kết hợp đồng thời với các kỹ thuật PIXE hay NRA.

Tại sao góc tán xạ ngược gần 180 độ luôn được ưu tiên bố trí trong thực nghiệm RBS?

Tại góc tán xạ $180^\circ$, hệ số động học tán xạ $K$ đạt giá trị cực tiểu và đạo hàm $\mathrm{d}K/\mathrm{d}M_2$ đạt giá trị lớn nhất. Điều này giúp độ tách biệt năng lượng $\Delta E_1$ giữa các nguyên tử bia có khối lượng khác nhau trở nên rõ rệt nhất, mang lại độ phân giải khối lượng tối ưu cho hệ phổ kế.

Vì sao giảm năng lượng chùm hạt tới lại giúp cải thiện độ phân giải theo chiều dày?

Khi giảm năng lượng chùm ion tới (trong dải từ 2600 keV xuống 1200 keV), tiết diện hãm năng lượng $[\varepsilon]$ của vật chất tăng lên theo lý thuyết Andersen - Ziegler. Do độ phân giải chiều dày $\delta Nt = \delta E / [\varepsilon]$, giá trị $[\varepsilon]$ càng lớn sẽ làm cho độ dày nhỏ nhất có thể phân giải được càng nhỏ, giúp xác định biên độ lớp màng mỏng chuẩn xác hơn.

Phần mềm RUT được xây dựng trong luận văn có điểm gì nổi bật so với các phần mềm quốc tế?

Phần mềm RUT được lập trình bằng Visual Basic .NET theo kiến trúc module DLL linh hoạt, tích hợp giải thuật Runge - Kutta Dormand - Prince tự thích nghi và dữ liệu tiết diện hãm Ziegler 1977. RUT có giao diện đồ họa trực quan, cho phép mô phỏng phổ RBS nhanh chóng, tính toán trực tiếp độ sâu đạt được và cho kết quả tương đương phần mềm chuẩn SIMNRA.

Hiện tượng nhòe năng lượng (Energy Straggling) ảnh hưởng như thế nào đến độ chính xác của phép đo?

Khi chùm ion di chuyển trong vật chất, các va chạm lượng tử ngẫu nhiên gây ra sự thăng giáng năng lượng xung quanh giá trị trung bình theo lý thuyết Bohr. Hiện tượng này làm mở rộng độ rộng vạch phổ (FWHM), làm mờ các biên phân chia lớp màng mỏng và giới hạn độ phân giải chiều dày cũng như độ phân giải khối lượng ở các lớp nằm sâu dưới bề mặt.

Kết luận

  • Luận văn đã hoàn thành xuất sắc việc xây dựng phần mềm mô phỏng phổ RBS chuyên dụng RUT với cấu trúc module DLL hướng đối tượng, tích hợp cơ sở dữ liệu tiết diện hãm hiện đại và thuật toán giải số tự thích nghi chính xác.
  • Nghiên cứu đã chứng minh bằng thực nghiệm và lý thuyết quy luật cải thiện 61,6% của độ phân giải khối lượng khi tăng năng lượng chùm hạt alpha từ 1000 keV lên 2600 keV trên mẫu chuẩn Au.
  • Xác định được điều kiện thực nghiệm tối ưu cho độ phân giải theo chiều dày tại mức năng lượng chùm tia 1200 keV kết hợp góc nghiêng mẫu $50^\circ$ trên buồng tán xạ hình học Cornell RC43.
  • Khảo sát toàn diện thang độ sâu phân tích có thể đạt tới (accessible depth) cho các chùm hạt $^1\text{H}$ và $^4\text{He}$ trên nhiều loại bia kim loại (Au, Ag, Ni, Al), đạt độ tương thích cao với các nghiên cứu quốc tế.
  • Đóng góp bộ quy trình thực nghiệm và công cụ tính toán chuẩn mực cho hệ máy gia tốc HUS 5SDH-2 Tandem Pelletron tại Trường Đại học Khoa học Tự nhiên – ĐHQGHN.

Kế hoạch phát triển tiếp theo bao gồm việc tích hợp hiệu chỉnh tán xạ nhiều lần vào phần mềm RUT trong vòng 6 tháng tới và triển khai các phép đo phân tích màng bán dẫn đa lớp phức tạp. Quý độc giả, nghiên cứu sinh và các phòng thí nghiệm quan tâm có thể khai thác toàn văn tài liệu và bộ mã nguồn phần mềm để ứng dụng vào công tác nghiên cứu khoa học vật liệu chuyên sâu.