Tổng quan luận án

Luận án tiến sĩ với đề tài "Nghiên cứu áp dụng một số phương pháp hạt nhân nguyên tử trong phân tích vật liệu TiO2/SiO2 sử dụng chùm ion từ máy gia tốc" (tên tiếng Anh: "Study on application of the nuclear methods for analysis TiO2/SiO2 material using accelerated ion beam") do nghiên cứu sinh Trần Văn Phúc thực hiện dưới sự hướng dẫn của GS. Lê Hồng Khiêm và TS. Mirosław Kulik tại Học viện Khoa học và Công nghệ thuộc Viện Hàn lâm Khoa học và Công nghệ Việt Nam. Luận án thuộc chuyên ngành Vật lý nguyên tử và hạt nhân (mã số: 9440106), được hoàn thành và bảo vệ năm 2023.

Tính cấp thiết và khoảng trống nghiên cứu

Công nghệ chùm ion được ứng dụng trong việc biến tính và phân tích các đặc trưng vật liệu màng mỏng, bề mặt và mặt phân cách. Trong các hệ màng kép chống phản xạ (antireflection coatings - ARC) tự làm sạch cho pin mặt trời silic, hệ hai lớp TiO2/SiO2 được sử dụng rộng rãi nhờ khả năng truyền qua quang học và độ bám dính. Lớp SiO2 có chiết suất thấp đóng vai trò làm giảm tổn hao phản xạ, trong khi lớp TiO2 cung cấp khả năng tự làm sạch qua cơ chế quang xúc tác và tính siêu ưa nước cảm ứng quang. Hiệu suất chống phản xạ phụ thuộc trực tiếp vào việc kiểm soát chiết suất và bề dày của các lớp vật liệu cũng như vùng mặt phân cách. Phương pháp trộn lẫn chùm ion (Ion Beam Mixing - IBM) cho phép điều chỉnh cấu trúc và bề dày lớp chuyển tiếp tại mặt phân cách giữa các lớp vật liệu.

Tuy nhiên, phần lớn các nghiên cứu IBM trước đây tập trung vào các hệ kim loại/kim loại, kim loại/silic hoặc kim loại/chất điện môi. Các dữ liệu thực nghiệm và cơ chế vi mô về hiện tượng trộn lẫn chùm ion trên hệ màng oxit/oxit, đặc biệt là hệ TiO2/SiO2 trong dải năng lượng ion thấp từ 100 keV đến 250 keV, còn rất hạn chế. Mối quan hệ định lượng giữa các thông số chùm ion (khối lượng, năng lượng) với mức độ trộn lẫn nguyên tử, sự biến đổi thành phần hóa học và các hằng số quang học của lớp chuyển tiếp chưa được khảo sát đầy đủ.

Mục tiêu nghiên cứu

Luận án xác định 5 mục tiêu nghiên cứu cụ thể:

  1. Xác định các biến đổi về cấu trúc độ dày của hệ vật liệu TiO2/SiO2/Si, bao gồm cả lớp chuyển tiếp giữa TiO2 và SiO2 do chiếu xạ ion khí hiếm trong dải năng lượng 100 - 250 keV bằng phương pháp phổ tán xạ ngược Rutherford (Rutherford Backscattering Spectrometry - RBS).
  2. Khảo sát sự phụ thuộc của mức độ trộn lẫn cảm ứng ion tại mặt phân cách TiO2/SiO2 vào năng lượng và khối lượng của các ion tới đối với các hệ màng có độ dày khác nhau.
  3. Giải thích cơ chế trộn lẫn nguyên tử theo động học truyền hạt bằng mô phỏng chương trình Dừng và Quãng chạy của Ion trong Vật chất (Stopping and Range of Ions in Matter - SRIM).
  4. Khảo sát ảnh hưởng của các biến đổi thành phần hóa học do chiếu xạ ion lên mức độ trộn lẫn theo năng lượng ion bằng phương pháp phổ quang điện tử tia X (X-ray Photoelectron Spectroscopy - XPS).
  5. Khảo sát sự biến đổi các thông số quang học (chiết suất, hệ số tắt, khe năng lượng quang học) của vùng chuyển tiếp TiO2/SiO2 trước và sau chiếu xạ theo năng lượng ion bằng phương pháp quang phổ elip (Ellipsometry Spectroscopy - ES).

Đối tượng và phạm vi nghiên cứu

  • Đối tượng nghiên cứu: Cấu trúc màng hai lớp TiO2/SiO2 trên đế Si (TiO2/SiO2/Si) được chế tạo theo hai nhóm độ dày khác nhau: Nhóm 1 có cấu trúc 20 nm TiO2 / 40 nm SiO2 / đế Si; Nhóm 2 có cấu trúc 30 nm TiO2 / 130 nm SiO2 / đế Si. Các chùm ion khí hiếm đơn hóa trị gồm $\text{Ne}^+$, $\text{Ar}^+$, $\text{Kr}^+$, $\text{Xe}^+$ ở 4 mức năng lượng 100, 150, 200 và 250 keV với liều chiếu xạ cố định $3 \times 10^{16}\text{ ions/cm}^2$.
  • Phạm vi nghiên cứu:
    • Thực nghiệm cấy ion được thực hiện trên máy gia tốc cấy ion UNIMAS tại Đại học Maria Curie-Skłodowska (MCSU), Lublin, Ba Lan ở nhiệt độ phòng.
    • Phân tích cấu trúc và thành phần theo chiều sâu bằng phương pháp RBS sử dụng chùm ion $\text{He}^+$ năng lượng 1,5 MeV từ máy gia tốc Van de Graaff tại nhóm EG-5, Phòng thí nghiệm Vật lý Neutron Frank, Viện Liên hiệp Nghiên cứu Hạt nhân (JINR), Dubna, Nga.
    • Phân tích quang phổ elip và phổ quang điện tử tia X được thực hiện tại các phòng thí nghiệm chuyên ngành của Đại học Maria Curie-Skłodowska (MCSU), Ba Lan.

Tổng quan tài liệu và vị trí của luận án

Các hướng nghiên cứu tiền đề

Hiện tượng trộn lẫn chùm ion được phát hiện lần đầu vào năm 1973 bởi Lee và cộng sự khi cấy ion P vào hệ Pd/Si, cùng công trình của Van der Weg và cộng sự khi bắn phá ion Ar lên hệ Pd/Si tạo thành pha palladium silicide. Trong các thập kỷ tiếp theo, nghiên cứu IBM được mở rộng trên nhiều hệ vật liệu:

  • Sự phụ thuộc vào liều ion ($\Phi$) và năng lượng truyền ($F_D$): Mô hình va chạm của Haff và Switkowski (1978) chỉ ra mức độ trộn lẫn phụ thuộc tuyến tính vào mật độ năng lượng tổn hao hạt nhân $F_D$. Tsaur và cộng sự (1979) cấy ion Xe 300 keV vào marker Pt trong Si; Matteson và cộng sự (1981) khảo sát các marker Ni, Ge, Pd, Sn, Sb, Au trong Si được chiếu bởi ion Ne, Ar, Kr, Xe đều ghi nhận sự tăng tuyến tính của độ trộn lẫn theo $F_D$. Datta (2003) nghiên cứu hệ màng hai lớp Au-Ni dưới chùm ion $\text{Kr}^+$ 300 keV; Bibic (2002) khảo sát hệ Ta/Si dưới chùm ion Ar; Dhar (2005) khảo sát hệ Pt/Co cấy ion $\text{Pt}^+$ 600 keV đều xác nhận độ rộng vùng trộn lẫn tăng tuyến tính với căn bậc hai của liều chiếu.
  • Ảnh hưởng của nhiệt độ: Phân tích của Nagel (2001) trên các hệ Ni/SiC, Cu/SiC, Pt/SiC dưới chùm ion $\text{Ar}^+$ 150 keV ở nhiệt độ từ 77 K đến 773 K và Qing Su (2016) trên composite SiOC/Fe chỉ ra sự tồn tại của nhiệt độ tới hạn, phân tách vùng trộn lẫn đạn đạo (ballistic mixing) ở nhiệt độ thấp và khuếch tán tăng cường bức xạ (radiation-enhanced diffusion) ở nhiệt độ cao.
  • Ảnh hưởng của khối lượng ion và các hiệu ứng phi tuyến: Averback (1986) chiếu xạ hệ Pt/Si bằng các ion từ He (4 amu) đến Xe (131 amu); J. Conrad (1993, 1996) cấy ion Ar, Kr, Xe, Pb vào hệ Pt/Ni và Ta/Si; F. Shi (1995) cấy ion vào hệ Sb/Ni. Các kết quả này chứng minh rằng ion nặng làm tăng mạnh hiệu suất trộn lẫn vượt qua dự đoán của mô hình va chạm đạn đạo thuần túy, do sự hình thành các đột biến nhiệt cục bộ (local thermal spikes) hoặc đột biến nhiệt toàn phần (global thermal spikes).
  • Nghiên cứu trên hệ oxit/oxit: Galuska và cộng sự (1991) nghiên cứu hệ $\text{Al}_2\text{O}_3/\text{SiO}_2$ chiếu ion $\text{Xe}^+$ và $\text{Kr}^+$, ghi nhận sự tăng cường liên kết Si-O-Al tại mặt phân cách. Ngược lại, Ogale và cộng sự (1987) cấy ion $\text{Kr}^{3+}$ 300 keV vào hệ $\text{Fe}_2\text{O}_3/\text{Al}_2\text{O}_3$ nhưng không quan sát thấy sự trộn lẫn trên phổ RBS. Kraft và cộng sự (2014) khảo sát hệ màng mỏng $\text{ZnO}/\text{SiO}_2$ dưới tác xạ ion nặng ở năng lượng 300 - 900 keV trong chế độ dừng hạt nhân, ghi nhận mức độ trộn lẫn nhỏ và phụ thuộc tuyến tính vào mật độ năng lượng lắng đọng.

Khoảng trống nghiên cứu luận án lựa chọn giải quyết

Luận án tập trung giải quyết khoảng trống thực nghiệm và lý thuyết đối với hệ màng hai lớp oxit/oxit $\text{TiO}_2/\text{SiO}_2$ trong dải năng lượng ion thấp 100 - 250 keV:

  1. Chưa có công trình nào khảo sát đầy đủ ảnh hưởng đồng thời của cả 4 loại ion khí hiếm ($\text{Ne}^+, \text{Ar}^+, \text{Kr}^+, \text{Xe}^+$) và 4 mức năng lượng (100, 150, 200, 250 keV) lên cấu trúc độ dày lớp chuyển tiếp $\text{TiO}_2/\text{SiO}_2$.
  2. Chưa làm rõ vai trò của độ dày màng ban đầu đối với số lượng sai hỏng dịch chuyển nguyên tử (Displacements Per Atom - DPA) và hiệu suất trộn lẫn.
  3. Thiếu dữ liệu thực nghiệm về mối liên hệ giữa sự biến đổi thành phần hóa học (các trạng thái liên kết của Ti) tại lớp bề mặt với động học truyền năng lượng ion.
  4. Chưa có nghiên cứu chi tiết về sự thay đổi các hằng số quang học ($n, k, E_g$) của lớp chuyển tiếp $\text{TiO}_2/\text{SiO}_2$ sau khi trộn lẫn bằng chùm ion.

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

Cơ sở lý thuyết và khung phân tích

Luận án xây dựng trên nền tảng lý thuyết tương tác giữa ion mang năng lượng với chất rắn:

  1. Động học va chạm đàn hồi hai vật: Bảo toàn năng lượng và động lượng trong hệ quy chiếu phòng thí nghiệm. Thừa số động học $K = E_1 / E_0$ được xác định bởi: $$K = \left[ \frac{\sqrt{M_2^2 - M_1^2 \sin^2\theta} + M_1 \cos\theta}{M_1 + M_2} \right]^2$$ trong đó $M_1$ là khối lượng ion tới, $M_2$ là khối lượng hạt nhân bia, $\theta$ là góc tán xạ.
  2. Tiết diện tán xạ vi phân Coulomb (Rutherford): Xác định xác suất tán xạ ion vào một góc khối $d\Omega$, tỷ lệ nghịch với bình phương năng lượng ion ($1/E^2$) và tăng mạnh khi góc tán xạ giảm.
  3. Quá trình tổn hao năng lượng ($dE/dx$): Gồm tổn hao hạt nhân $(dE/dx)_n$ do va chạm đàn hồi trực tiếp với hạt nhân bia và tổn hao điện tử $(dE/dx)_e$ do kích thích, ion hóa đám mây điện tử: $$\frac{dE}{dx} = \left.\frac{dE}{dx}\right|_n + \left.\frac{dE}{dx}\right|_e$$ Ở dải năng lượng thấp (< 250 keV), tổn hao hạt nhân chiếm ưu thế, trực tiếp tạo ra sự dịch chuyển nguyên tử khỏi vị trí nút mạng và gây ra hiện tượng trộn lẫn nguyên tử.
  4. Cơ chế trộn lẫn chùm ion (IBM):
    • Trộn lẫn giật lùi (Recoil mixing): Do va chạm trực diện đơn lẻ giữa ion tới và nguyên tử bia, đẩy nguyên tử bia xuyên sâu qua mặt phân cách.
    • Trộn lẫn dòng thác (Cascade mixing): Do chuỗi các va chạm thứ cấp từ nguyên tử bị giật lùi ban đầu (Primary Knock-on Atom - PKA), tạo ra sự chuyển dời ngẫu nhiên, đẳng hướng của nhiều nguyên tử qua mặt phân cách.
    • Hệ số khuếch tán hiệu dụng dòng thác: Theo mô hình Andersen (1979) và mô hình Sigmund - Gras-Marti (1981): $$D_{cas} t = \frac{\Gamma_m F_D(x) \langle r^2 \rangle \Phi}{6 E_d N}$$ với $\Gamma_m = 0{,}608 \left[ \frac{4 M_1 M_2}{(M_1 + M_2)^2} \right]^{1/2}$, $F_D$ là năng lượng tạo sai hỏng, $E_d$ là năng lượng ngưỡng dịch chuyển nguyên tử, $N$ là mật độ nguyên tử bia, $\Phi$ là liều chiếu xạ, và $\langle r^2 \rangle$ là quãng chạy bình phương trung bình của nguyên tử bị dịch chuyển.

Phương pháp thực nghiệm và kỹ thuật phân tích

  • Chế tạo mẫu và chiếu xạ cấy ion:
    • Tổng số 40 mẫu cấu trúc $\text{TiO}_2/\text{SiO}_2/\text{Si}$ được chia thành 2 nhóm độ dày và chiếu xạ tại máy gia tốc cấy ion UNIMAS (Đại học Maria Curie-Skłodowska, Ba Lan).
    • Nguồn ion phóng điện hồ quang (ARC discharge), điện áp trích 25 kV, hệ nam châm điện phân tách khối lượng góc $90^\circ$.
    • Mức chân không buồng chiếu $10^{-6}\text{ mbar}$, mật độ dòng ion duy trì $\sim 1\ \mu\text{A/cm}^2$ ở nhiệt độ phòng.
Nhóm mẫu Cấu trúc lớp ban đầu Ion chiếu xạ Mức năng lượng (keV) Liều chiếu ($\text{ions/cm}^2$)
Nhóm 1 $20\text{ nm TiO}_2\ /\ 40\text{ nm SiO}_2\ /\ \text{Si}$ $\text{Ne}^+, \text{Ar}^+, \text{Kr}^+, \text{Xe}^+$ 100, 150, 200, 250 $3 \times 10^{16}$
Nhóm 2 $30\text{ nm TiO}_2\ /\ 130\text{ nm SiO}_2\ /\ \text{Si}$ $\text{Ne}^+, \text{Ar}^+, \text{Kr}^+, \text{Xe}^+$ 100, 150, 200, 250 $3 \times 10^{16}$
  • Phương pháp phổ tán xạ ngược Rutherford (RBS):
    • Thực hiện trên máy gia tốc Van de Graaff tại JINR Dubna.
    • Sử dụng chùm ion $\text{He}^+$ năng lượng 1,5 MeV, độ phân kỳ chùm nhỏ, đường kính vết chùm $\sim 5\text{ mm}$, độ phân giải năng lượng của hệ đo là 15 keV với 1024 kênh phân tích.
    • Cấu hình hình học IBM: góc tới $\alpha = 60^\circ$, góc tán xạ $\theta = 170^\circ$, góc thoát $\beta = 50^\circ$.
    • Xử lý và khớp phổ thực nghiệm bằng phần mềm SIMNRA để xác định bề dày các lớp (tính theo $\text{nguyên tử/cm}^2$ và chuyển đổi sang nm) và phân bố nồng độ các nguyên tố Ti, Si, O theo độ sâu.
  • Phương pháp mô phỏng SRIM/TRIM:
    • Chế độ "Detailed calculation with full damage cascades" sử dụng 99.999 ion tới.
    • Tính toán quãng chạy hình chiếu ($R_p$), độ tán xạ quãng chạy ($\Delta R_p$), mật độ năng lượng tổn hao hạt nhân và điện tử, số lượng sai hỏng trên mỗi nguyên tử (DPA), và phân bố lỗ trống tạo ra tại vùng chuyển tiếp.
  • Phương pháp quang phổ elip (ES / SE):
    • Đo các góc elip $\Psi$ (tỷ số biên độ) và $\Delta$ (độ lệch pha) trong dải bước sóng quang học.
    • Áp dụng mô hình quang học nhiều lớp kết hợp mô hình xấp xỉ môi trường hiệu dụng (Effective Medium Approximation - EMA) để xác định bề dày lớp chuyển tiếp, chiết suất $n(\lambda)$, hệ số tắt $k(\lambda)$ và khe năng lượng quang học $E_g$ thông qua hệ thức Tauc.
  • Phương pháp phổ quang điện tử tia X (XPS):
    • Hệ chân không siêu cao UHV tại Khoa Hóa học, MCSU.
    • Thu nhận phổ vạch Ti 2p để phân tích tỷ lệ nồng độ tương đối của các trạng thái hóa học $\text{Ti}, \text{TiO}, \text{Ti}_2\text{O}_3, \text{TiO}_2$ tại lớp gần bề mặt trước và sau khi cấy ion.

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

Chương 1: Cơ sở lý thuyết (Theoretical Background)

Chương 1 trình bày tổng quan các nghiên cứu thực nghiệm và mô hình lý thuyết về quá trình trộn lẫn chùm ion, cùng các nguyên lý vật lý chi phối tương tác ion - chất rắn ở dải năng lượng thấp:

  • Cấu hình mẫu nghiên cứu IBM: Phân tích ba dạng cấu hình chuẩn gồm lớp đánh dấu (marker), màng hai lớp (bilayer) và màng đa lớp (multilayer). Luận án lựa chọn cấu hình màng hai lớp $\text{TiO}_2/\text{SiO}_2$ trên đế Si.
  • Động học va chạm và năng lượng dừng: Trình bày chi tiết động học va chạm đàn hồi trong hệ quy chiếu khối tâm và hệ quy chiếu phòng thí nghiệm; chứng minh góc tán xạ cực đại $\theta_m$ phụ thuộc tỷ số khối lượng $M_1/M_2$. Thiết lập đường cong năng lượng dừng ion $S(v)$ với ba vùng vận tốc ion khác nhau: vùng I (năng lượng thấp, dừng hạt nhân chiếm ưu thế), vùng II (dừng điện tử theo lý thuyết Bohr), vùng III (va chạm nhanh).
  • Cơ chế biến tính chất rắn bằng ion năng lượng thấp: Phân tích 4 quá trình cơ bản xảy ra đồng thời gồm: cấy ion (implantation), tổn hại bức xạ (radiation damage), trộn lẫn chùm ion (IBM), và phún xạ bề mặt (sputtering).
  • Phân tích cơ chế trộn lẫn đạn đạo:
    • Trộn lẫn giật lùi: Mô tả quá trình một ion truyền động năng lớn cho một nguyên tử bia qua va chạm trực diện ($\theta = 0$). Phân tích chỉ ra xác suất va chạm trực diện rất thấp, phần lớn là va chạm mềm ($\theta > 0$) với năng lượng truyền nhỏ, do đó trộn lẫn giật lùi sơ cấp không thể giải thích toàn bộ mức độ trộn lẫn quan sát được.
    • Trộn lẫn dòng thác: Mô tả quá trình va chạm đa dịch chuyển liên tiếp tạo thành dòng thác va chạm. Các nguyên tử chuyển động ngẫu nhiên đẳng hướng với năng lượng xấp xỉ năng lượng ngưỡng dịch chuyển $E_d$. Hệ số khuếch tán dòng thác $D_{cas}t$ tỷ lệ thuận với mật độ năng lượng tổn hại $F_D$ và liều ion $\Phi$, độc lập với nhiệt độ ở vùng nhiệt độ thấp.

Chương 2: Phương pháp thực nghiệm (The Experimental Methods)

Chương 2 mô tả chi tiết quy trình chế tạo, điều kiện chiếu xạ và các phương pháp phân tích thực nghiệm:

  • Hệ màng chống phản xạ $\text{TiO}_2/\text{SiO}_2$ cho pin mặt trời: Phân tích cơ chế giao thoa triệt tiêu ánh sáng phản xạ của hệ màng kép, giúp giảm hệ số phản xạ của tấm pin silic từ 36% xuống 15% (đối với đơn lớp $\text{SiO}_2$) và xuống 7% (đối với màng kép $\text{TiO}_2/\text{SiO}_2$).
  • Hệ thống cấy ion UNIMAS: Mô tả sơ đồ cấu tạo và nguyên lý hoạt động của máy gia tốc cấy ion UNIMAS tại MCSU (Ba Lan). Chùm ion sau khi trích ở điện áp 25 kV được hội tụ bằng bộ ba thấu kính tĩnh điện, phân tách khối lượng bằng từ trường quét góc $90^\circ$, sau đó gia tốc đến năng lượng yêu cầu (100 - 250 keV) với độ chính xác năng lượng tốt hơn 1%. Hệ thống quét chùm 2 chiều bảo đảm độ đồng đều không gian trên bề mặt mẫu $> 99%$.
  • Mô phỏng SRIM/TRIM: Thiết lập bài toán mô phỏng chuyển dời ion $\text{Ne}^+, \text{Ar}^+, \text{Kr}^+, \text{Xe}^+$ trong cấu trúc đa lớp $\text{TiO}_2/\text{SiO}_2/\text{Si}$ với 99.999 hạt; trích xuất các thông số phân bố độ sâu, số nguyên tử Ti và Si bị dịch chuyển, mật độ lỗ trống và năng lượng mất mát $dE/dx$.
  • Hệ đo phổ RBS tại JINR Dubna: Mô tả máy gia tốc Van de Graaff, hệ nam châm lái chùm, buồng bia chân không $10^{-6}\text{ Torr}$, đầu dò rào cản bề mặt bán dẫn với độ phân giải 15 keV nối với bộ phân tích đa kênh (MCA) 1024 kênh. Thiết lập cấu hình đo nghiêng IBM ($\alpha = 60^\circ, \theta = 170^\circ$) để tăng độ phân giải độ sâu tại vùng mặt phân cách. Quy trình khớp phổ bằng mã SIMNRA được chuẩn hóa để xác định độ dày lớp chuyển tiếp.
  • Phương pháp quang phổ elip (ES) và XPS: Trình bày nguyên lý đo phân cực elip (RAE - Rotating Analyzer Ellipsometer), mô hình môi trường hiệu dụng EMA để tính toán thành phần thể tích các pha trong lớp chuyển tiếp; nguyên lý đo quang điện tử tia X để xác định trạng thái hóa trị của Ti tại vùng bề mặt.

Chương 3: Kết quả và thảo luận (Results and Discussion)

Chương 3 trình bày các kết quả thực nghiệm và phân tích định lượng về quá trình trộn lẫn ion cùng sự biến đổi tính chất lý hóa:

1. Đặc trưng phổ RBS và sự hình thành lớp chuyển tiếp

  • Phổ RBS của các mẫu trước khi chiếu xạ (virgin) thể hiện các đỉnh phổ rõ nét, bờ sau của đỉnh Ti và bờ trước của nền $\text{SiO}_2/\text{Si}$ có độ dốc đứng, phản ánh mặt phân cách sắc nét.
  • Sau khi chiếu xạ ion ($\text{Ne}^+, \text{Ar}^+, \text{Kr}^+, \text{Xe}^+$), bờ tín hiệu Ti và Si bị choãi rộng ra, biên độ đỉnh Ti giảm xuống. Hiện tượng này chứng minh sự khuếch tán tương hỗ giữa các nguyên tử Ti vào lớp $\text{SiO}_2$ và Si/O vào lớp $\text{TiO}_2$, hình thành lớp chuyển tiếp hỗn hợp ($\text{TiO}_2/\text{SiO}_2$ mixed layer).
  • Độ mở rộng độ rộng nửa chiều cao ($\Delta\text{FWHM}$) của đỉnh phổ Ti tăng tuyến tính theo năng lượng ion tới, thể hiện mức độ trộn lẫn tăng dần.

2. Sự phụ thuộc của mức độ trộn lẫn vào năng lượng ion tới

  • Bề dày thực tế của lớp chuyển tiếp ($t$) được xác định từ phần mềm SIMNRA (đơn vị $\text{nguyên tử/cm}^2$) và quy đổi sang nanomet (nm).
  • Tỷ số bề dày tương đối $r_t = t_{irradiated} / t_{un-irradiated}$ tăng tuyến tính theo năng lượng ion chiếu xạ từ 100 keV đến 250 keV đối với tất cả các loại ion.
  • Kết quả mô phỏng SRIM cho thấy khi năng lượng ion tăng, quãng chạy hình chiếu $R_p$ tăng sâu hơn vào vùng mặt phân cách, năng lượng truyền cho các nguyên tử giật lùi tăng lên, đẩy các nguyên tử Ti và Si đi xa hơn vị trí ban đầu.

3. Sự phụ thuộc của mức độ trộn lẫn vào khối lượng ion tới

  • Khảo sát 4 loại ion với khối lượng nguyên tử tăng dần: Ne (20 amu), Ar (40 amu), Kr (84 amu), Xe (131 amu).
  • Mức độ trộn lẫn tăng mạnh theo khối lượng của ion tới. Độ dốc của hàm khớp tuyến tính bề dày tương đối theo năng lượng tăng rõ rệt từ ion nhẹ (Ne) đến ion nặng (Xe).
  • Mô phỏng SRIM giải thích: Các ion nặng hơn có tiết diện va chạm hạt nhân lớn hơn nhiều, tạo ra mật độ năng lượng tổn hại hạt nhân $F_D$ và số lượng va chạm thứ cấp trong dòng thác lớn hơn tại vùng mặt phân cách, dẫn đến số lượng nguyên tử bị dịch chuyển trên mỗi ion tới cao hơn đáng kể.

4. Ảnh hưởng của độ dày lớp màng ban đầu

  • So sánh giữa Nhóm 1 ($20\text{ nm TiO}_2 / 40\text{ nm SiO}_2$) và Nhóm 2 ($30\text{ nm TiO}_2 / 130\text{ nm SiO}_2$):
  • Mức độ trộn lẫn tương đối ($r_t$) và tốc độ tăng bề dày lớp chuyển tiếp ở Nhóm 1 cao hơn so với Nhóm 2 ở cùng điều kiện chiếu xạ.
  • Giá trị DPA (Displacements Per Atom) tính toán tại độ sâu 20 nm dưới đáy lớp $\text{TiO}_2$ cho thấy ở màng mỏng hơn (Nhóm 1), ion tới mất ít năng lượng hơn ở lớp phủ bề mặt, do đó năng lượng khả dụng khi đến mặt phân cách $\text{TiO}_2/\text{SiO}_2$ lớn hơn, tạo ra số sai hỏng DPA cao hơn.
Nhóm mẫu Ion Năng lượng (keV) Bề dày lớp chuyển tiếp RBS (nm) $r_t$ (Bề dày tương đối) DPA tại vùng phân cách (SRIM)
Nhóm 1 Virgin 0 5,2 1,00 -
(20/40 nm) $\text{Ne}^+$ 100 12,1 2,33 14,2
$\text{Ne}^+$ 250 18,5 3,56 23,8
$\text{Xe}^+$ 100 24,6 4,73 48,6
$\text{Xe}^+$ 250 41,2 7,92 82,4
Nhóm 2 Virgin 0 6,1 1,00 -
(30/130 nm) $\text{Ne}^+$ 100 9,8 1,61 8,7
$\text{Ne}^+$ 250 14,2 2,33 15,1
$\text{Xe}^+$ 100 18,3 3,00 31,5
$\text{Xe}^+$ 250 32,5 5,33 56,2

5. Biến đổi thành phần hóa học lớp bề mặt bằng XPS

  • Phổ XPS vùng Ti 2p của các mẫu trước và sau khi cấy ion $\text{Ne}^+$ ở các mức năng lượng 100, 150, 200, 250 keV được phân tách thành 4 thành phần hóa học: $\text{Ti}$ kim loại, $\text{TiO}, \text{Ti}_2\text{O}_3$ và $\text{TiO}_2$.
  • Khi năng lượng ion tăng từ 100 keV lên 250 keV, tỷ lệ nồng độ tương đối của $\text{TiO}_2$ tăng dần, trong khi nồng độ của các hợp chất khử ($\text{Ti}, \text{TiO}, \text{Ti}_2\text{O}_3$) giảm xuống.
  • Kết quả mô phỏng SRIM cho thấy tổn hao năng lượng của ion $\text{Ne}^+$ trong lớp $\text{TiO}_2$ giảm khi năng lượng ion tới tăng (do phân bố đỉnh dừng hạt nhân chuyển dịch sâu hơn vào trong). Điều này đồng nghĩa với việc mức độ phá hủy liên kết hóa học tại lớp sát bề mặt giảm, ion truyền năng lượng tương tác lớn hơn xuống vùng phân cách sâu bên dưới, thúc đẩy mức độ trộn lẫn $\text{TiO}_2/\text{SiO}_2$.

6. Tính chất quang học của lớp chuyển tiếp bằng phổ elip (SE)

  • Dữ liệu phổ elip ($\Psi, \Delta$) được xử lý qua mô hình môi trường hiệu dụng EMA với cấu trúc 3 pha thành phần trong lớp chuyển tiếp: $\text{TiO}_2$, $\text{SiO}_2$ và lỗ xốp (voids).
  • Bề dày lớp chuyển tiếp xác định từ phổ elip có sự phù hợp tốt với kết quả thu được từ phổ RBS.
  • Chiết suất ($n$) và hệ số tắt ($k$) của lớp chuyển tiếp tăng khi năng lượng ion cấy tăng từ 100 keV đến 200 keV, sau đó giảm nhẹ ở mức 250 keV. Sự biến đổi này được giải thích bởi sự thay đổi hàm lượng phần trăm thể tích của pha $\text{TiO}_2$ (vốn có chiết suất cao $n \approx 2{,}4 - 2{,}6$) dịch chuyển vào lớp chuyển tiếp dưới tác động của các mức sai hỏng khác nhau.
  • Khe năng lượng quang học ($E_g$) của lớp chuyển tiếp được xác định từ đồ thị Tauc phụ thuộc vào mức năng lượng ion cấy, phản ánh sự tái cấu trúc mạng tinh thể và biến đổi mật độ trạng thái điện tử cục bộ sau chiếu xạ.

Kết quả và những đóng góp mới

Đóng góp mới về mặt lý luận và khoa học

  1. Xác lập cơ chế trộn lẫn đạn đạo cho hệ oxit/oxit: Luận án chứng minh quá trình trộn lẫn nguyên tử tại mặt phân cách $\text{TiO}_2/\text{SiO}_2$ dưới tác xạ ion khí hiếm năng lượng 100 - 250 keV xảy ra chủ yếu theo cơ chế dòng thác va chạm đạn đạo (cascade mixing), được điều khiển bởi tổn hại dừng hạt nhân $(dE/dx)_n$.
  2. Quy luật phụ thuộc năng lượng và khối lượng ion: Xác định mức độ trộn lẫn nguyên tử tăng tuyến tính với năng lượng ion tới và tăng phi tuyến mạnh theo khối lượng nguyên tử của ion chiếu xạ ($\text{Ne} < \text{Ar} < \text{Kr} < \text{Xe}$).
  3. Ảnh hưởng của độ dày lớp màng: Phát hiện và định lượng hiệu ứng độ dày màng ban đầu: cấu trúc màng mỏng ($20\text{ nm TiO}_2 / 40\text{ nm SiO}_2$) chịu mức độ tạo sai hỏng DPA tại mặt phân cách cao hơn, dẫn đến hiệu suất trộn lẫn nguyên tử vượt trội so với cấu trúc màng dày ($30\text{ nm TiO}_2 / 130\text{ nm SiO}_2$).
  4. Mối liên hệ giữa biến đổi hóa học, quang học và trộn lẫn: Lần đầu tiên thiết lập mối liên hệ định lượng giữa sự suy giảm mất mát năng lượng ion ở lớp bề mặt $\text{TiO}_2$ (khảo sát qua XPS và SRIM) với sự tăng cường mức độ trộn lẫn tại mặt phân cách, cũng như sự biến đổi phi đơn điệu của các hằng số quang học ($n, k$) của lớp chuyển tiếp theo năng lượng ion cấy.

Đóng góp thực tiễn và kiến nghị ứng dụng

  • Cung cấp bộ thông số thực nghiệm chuẩn xác về độ dày, thành phần pha, chiết suất và hệ số tắt của lớp chuyển tiếp $\text{TiO}_2/\text{SiO}_2$ được biến tính bằng chùm ion.
  • Đề xuất phương pháp sử dụng chùm ion khí hiếm năng lượng thấp như một công cụ kiểm soát chính xác bề dày vùng chuyển tiếp mặt phân cách nhằm tối ưu hóa độ bám dính cơ học và tính chất quang học cho các hệ màng phủ chống phản xạ đa lớp ứng dụng trong công nghệ pin mặt trời quang điện.

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

Hạn chế của luận án

  • Thực nghiệm cấy ion chỉ được tiến hành ở nhiệt độ phòng, chưa khảo sát ảnh hưởng của nhiệt độ đế (nhiệt độ thấp < 80 K hoặc nhiệt độ cao > 373 K) để phân tách định lượng đóng góp của hiệu ứng đột biến nhiệt (thermal spike) và khuếch tán tăng cường bức xạ.
  • Khảo sát thực nghiệm được cố định ở một mức liều chiếu duy nhất ($3 \times 10^{16}\text{ ions/cm}^2$), chưa khảo sát đường cong động học trộn lẫn phụ thuộc liều ion chi tiết trên hệ $\text{TiO}_2/\text{SiO}_2$.
  • Các ion sử dụng giới hạn trong nhóm ion khí hiếm trơ ($\text{Ne}^+, \text{Ar}^+, \text{Kr}^+, \text{Xe}^+$), chưa mở rộng sang các ion hoạt tính hóa học (như ion kim loại hoặc ion oxy).

Hướng nghiên cứu tiếp theo

  • Mở rộng khảo sát động học quá trình trộn lẫn ion theo dải liều chiếu biến thiên từ $10^{15}$ đến $10^{17}\text{ ions/cm}^2$.
  • Nghiên cứu ảnh hưởng của quá trình ủ nhiệt sau chiếu xạ (post-implantation annealing) đến sự phục hồi sai hỏng mạng, sự tái kết tinh các pha oxit và độ bền quang học của lớp chuyển tiếp.
  • Khảo sát trực tiếp hiệu suất biến đổi quang điện của pin mặt trời silic thực tế khi được phủ lớp màng $\text{TiO}_2/\text{SiO}_2$ đã qua xử lý trộn lẫn bằng chùm ion.

Giá trị tham khảo

Luận án là tài liệu tham khảo chuyên khảo có giá trị trong các lĩnh vực:

  • Nghiên cứu sinh, học viên cao học chuyên ngành Vật lý nguyên tử và hạt nhân, Vật lý chất rắn: Tham khảo phương pháp luận kết hợp giữa phân tích chùm ion hạt nhân (IBA/RBS), mô phỏng truyền hạt (SRIM/TRIM) và các phương pháp quang phổ bề mặt (XPS, SE).
  • Giảng viên và nhà nghiên cứu công nghệ vật liệu màng mỏng: Tham khảo các mô hình va chạm nguyên tử, cơ chế dừng hạt nhân, và phương pháp xác định cấu trúc mặt phân cách vật liệu oxit vô cơ.
  • Kỹ sư và chuyên gia công nghệ quang điện tử, pin mặt trời: Tham khảo dữ liệu thực nghiệm về điều khiển chiết suất $n$, hệ số tắt $k$ và bề dày lớp chuyển tiếp $\text{TiO}_2/\text{SiO}_2$ nhằm cải thiện độ bền bám dính và hiệu suất quang học của màng chống phản xạ.

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

1. Cơ chế vật lý chủ đạo chi phối quá trình trộn lẫn nguyên tử tại mặt phân cách $\text{TiO}_2/\text{SiO}_2$ trong luận án là gì?

Cơ chế chủ đạo là quá trình trộn lẫn dòng thác va chạm đạn đạo (cascade mixing). Khi các ion khí hiếm mang năng lượng 100 - 250 keV đi vào vật chất, tương tác dừng hạt nhân $(dE/dx)_n$ chiếm ưu thế, tạo ra các nguyên tử giật lùi sơ cấp (PKA). Các PKA này tiếp tục va chạm đa dịch chuyển liên tiếp với các nguyên tử lân cận, tạo nên dòng thác va chạm đẩy các nguyên tử Ti vào lớp $\text{SiO}_2$ và Si/O vào lớp $\text{TiO}_2$ tại mặt phân cách.

2. Năng lượng và khối lượng ion ảnh hưởng như thế nào đến bề dày lớp chuyển tiếp?

Bề dày tương đối của lớp chuyển tiếp ($r_t$) tăng tuyến tính với năng lượng ion tới (từ 100 đến 250 keV) do ion năng lượng cao xâm nhập sâu hơn và truyền động năng lớn hơn cho các nguyên tử giật lùi. Đồng thời, mức độ trộn lẫn tăng mạnh phi tuyến theo khối lượng ion ($\text{Ne}^+ < \text{Ar}^+ < \text{Kr}^+ < \text{Xe}^+$) do ion nặng có tiết diện tán xạ hạt nhân lớn hơn, tạo ra mật độ năng lượng tổn hại $F_D$ và số sai hỏng dịch chuyển nguyên tử tại mặt phân cách cao hơn.

3. Tại sao hệ màng có độ dày mỏng hơn (Nhóm 1) lại có mức độ trộn lẫn lớn hơn hệ màng dày (Nhóm 2)?

Ở hệ màng mỏng ($20\text{ nm TiO}_2 / 40\text{ nm SiO}_2$), ion tới tiêu hao ít năng lượng hơn khi đi qua lớp phủ $\text{TiO}_2$ bề mặt so với lớp phủ 30 nm của Nhóm 2. Do đó, năng lượng khả dụng của ion khi tiếp cận mặt phân cách $\text{TiO}_2/\text{SiO}_2$ ở Nhóm 1 cao hơn, tạo ra số sai hỏng trên mỗi nguyên tử (DPA) lớn hơn (ví dụ với ion Xe 250 keV, DPA của Nhóm 1 đạt 82,4 so với 56,2 của Nhóm 2), dẫn đến hiệu suất trộn lẫn nguyên tử cao hơn.

4. Phổ XPS ghi nhận sự thay đổi thành phần hóa học tại lớp bề mặt $\text{TiO}_2$ như thế nào khi tăng năng lượng ion chiếu xạ?

Phổ XPS vùng Ti 2p cho thấy khi năng lượng ion $\text{Ne}^+$ tăng từ 100 keV lên 250 keV, nồng độ tương đối của hợp chất $\text{TiO}_2$ tăng dần, trong khi nồng độ các trạng thái khử ($\text{Ti}, \text{TiO}, \text{Ti}_2\text{O}_3$) giảm xuống. Điều này phù hợp với tính toán SRIM cho thấy tổn hao năng lượng của ion tại lớp bề mặt giảm khi năng lượng ion tăng, làm giảm sự phá hủy liên kết hóa học ở bề mặt và tập trung năng lượng tương tác xuống vùng phân cách sâu hơn.

5. Chiết suất ($n$) và hệ số tắt ($k$) của lớp chuyển tiếp biến đổi ra sao sau khi chiếu xạ ion?

Theo phân tích quang phổ elip sử dụng mô hình EMA, chiết suất $n$ và hệ số tắt $k$ của lớp chuyển tiếp tăng khi năng lượng ion cấy tăng từ 100 keV đến 200 keV, sau đó có xu hướng giảm nhẹ ở mức 250 keV. Sự biến đổi này bắt nguồn từ sự thay đổi hàm lượng nồng độ pha $\text{TiO}_2$ (pha có chiết suất cao) thâm nhập vào vùng chuyển tiếp dưới tác động của mật độ sai hỏng được tạo ra ở các mức năng lượng ion khác nhau.

Kết luận

Luận án của nghiên cứu sinh Trần Văn Phúc đã nghiên cứu thực nghiệm và mô phỏng cơ chế trộn lẫn chùm ion (IBM) trên hệ màng hai lớp oxit/oxit $\text{TiO}_2/\text{SiO}_2/\text{Si}$ sử dụng 4 loại ion khí hiếm ($\text{Ne}^+, \text{Ar}^+, \text{Kr}^+, \text{Xe}^+$) trong dải năng lượng 100 - 250 keV. Bằng sự kết hợp giữa phương pháp phổ tán xạ ngược Rutherford (RBS), mô phỏng động học SRIM, quang phổ elip (SE) và phổ quang điện tử tia X (XPS), công trình đã định lượng chính xác sự phụ thuộc của bề dày lớp chuyển tiếp, thành phần hóa học và các hằng số quang học vào năng lượng, khối lượng ion cũng như độ dày ban đầu của màng. Kết quả nghiên cứu cung cấp cơ sở dữ liệu thực nghiệm và lý thuyết cho việc ứng dụng công nghệ chùm ion hạt nhân trong chế tạo và kiểm soát chất lượng màng phủ chống phản xạ quang học đa chức năng.