Tổng quan về luận án
Nghiên cứu vật liệu nano dị thể (heterojunction nanocomposites) giữa graphene và titanium dioxide ($TiO_2$) đang là trọng tâm của khoa học vật liệu và hóa học lượng tử ứng dụng trong lĩnh vực xúc tác quang hóa. Hệ composite này giữ vai trò then chốt trong các phản ứng phân hủy hợp chất hữu cơ độc hại từ nước thải công nghiệp, quang phân tách nước tạo hydro năng lượng sạch, ứng dụng làm photoanode trong pin mặt trời chất màu nhạy quang (DSSC) [1] và cực anode cho pin lithium-ion [2]. Luận án tiến sĩ chuyên ngành Hóa lí thuyết và Hóa lí của nghiên cứu sinh Trần Thị Thoa (2022) tại Trường Đại học Sư phạm Hà Nội với đề tài "Nghiên cứu cấu trúc, tính chất của các dẫn xuất graphene và rutile $TiO_2$ trong mô hình composite bằng phương pháp phiếm hàm mật độ" đã giải quyết những khoảng trống học thuật cốt lõi trong mô phỏng cơ lượng tử cấu trúc giao diện vật liệu 2D/3D.
Khoảng trống nghiên cứu (Research Gap) then chốt xuất phát từ sự bất đối xứng hình học và hạn chế trong các mô hình tính toán trước đây. Đa số các nghiên cứu lý thuyết quốc tế trước đây như công trình của Long và cộng sự [56] hay A. Du và cộng sự [57] chủ yếu tiếp cận composite $TiO_2$/graphene dựa trên ô đơn vị lục giác hai nguyên tử truyền thống của graphene ($a = b = 2.463\text{ \AA}, \gamma = 120^\circ$). Khi ghép nối với bề mặt rutile $TiO_2(110)$ – pha tinh thể bền vững nhất về mặt nhiệt động học với ô đơn vị hình chữ nhật ($a = 4.595\text{ \AA}, b = 2.959\text{ \AA}, \alpha = \beta = \gamma = 90^\circ$) – sự sai lệch góc mạng ($90^\circ$ so với $120^\circ$) tạo ra ứng suất nhân tạo rất lớn làm biến dạng phi thực tế cấu trúc điện tử của hệ. Mặt khác, trong thực nghiệm điều chế vật liệu, sản phẩm thu được chủ yếu là dạng khử graphene oxide (RGO) hoặc graphene oxide (GO) với các nhóm chức bề mặt như hydroxyl (-OH) và epoxy (-O-). Các nghiên cứu quốc tế trước đó chỉ gắn ngẫu nhiên nhóm chức (như mô hình của Yan et al. [15] hoặc Lahaye et al. [73]) mà chưa làm sáng tỏ được cơ chế phân bố không gian và trật tự liên kết vi mô của cụm đa nhóm chức hydroxyl ($G_nOH$) trên mặt đáy (basal plane) của graphene.
Hệ thống câu hỏi nghiên cứu (Research Questions - RQ) và giả thuyết khoa học (Hypotheses - H) được xác lập chặt chẽ:
- RQ1: Quy luật sắp xếp tối ưu về mặt nhiệt động và năng lượng liên kết của các nhóm hydroxyl trong dẫn xuất $G_nOH$ ($n = 1 - 6$) trên bề mặt graphene diễn ra theo cơ chế nào?
- H1: Các nhóm -OH có xu hướng tự tập hợp thành các cụm tuần hoàn tại vị trí para nhờ sự hiệp đồng của mạng lưới liên kết hydrogen nội và liên phân tử.
- RQ2: Mức độ bền vững và cấu trúc dải năng lượng của các dẫn xuất graphene đơn nhóm chức $G_1F$ (-NH2, -CH3, -OCH3, -CHO, -COOH, epoxy) biến thiên như thế nào theo bản chất liên kết hóa học?
- H2: Độ bền liên kết $E_b$ tỷ lệ thuận với khả năng lai hóa orbital $2p_z$ của dị nguyên tử với mạng liên hợp $\pi$ của graphene; nhóm epoxy và methyl có độ bền vượt trội so với carboxyl.
- RQ3: Phương pháp hiệu chỉnh Coulomb nội tại Hubbard ($DFT+U$) khắc phục sai số tự tương tác (SIE) của DFT truyền thống đối với rutile $TiO_2$ khối và bề mặt (110) ra sao khi đồng thời khảo sát $U_d$ trên Ti-$3d$ và $U_p$ trên O-$2p$?
- H3: Việc kết hợp tham số Hubbard tối ưu $U_p = 10\text{ eV}$ cho orbital $2p$ của oxy cùng $U_d = 7\text{ eV}$ cho orbital $3d$ của titan sẽ mở rộng độ rộng vùng cấm ($E_g$) tiệm cận thực nghiệm $3.03\text{ eV}$ mà không phá vỡ hằng số mạng tinh thể.
- RQ4: Mô hình ô đơn vị hình chữ nhật bốn nguyên tử của graphene triệt tiêu sai lệch mạng tinh thể (lattice mismatch) như thế nào khi xây dựng dị cấu trúc 2D với rutile $TiO_2(110)$?
- H4: Ô đơn vị hình chữ nhật bốn nguyên tử ($a_3 = 2.463\text{ \AA}, a_4 = 4.266\text{ \AA}$) cho phép khớp nối trực giao chính xác với mặt (110) của rutile, loại trừ hoàn toàn sai lệch góc $30^\circ$.
Khung lý thuyết của luận án xây dựng vững chắc trên nền tảng cơ học lượng tử phiếm hàm mật độ (Density Functional Theory - DFT) thông qua các định lý Hohenberg-Kohn [16] và hệ phương trình tự nhất quán Kohn-Sham [17], tích hợp phiếm hàm tương quan không định xứ van der Waals optPBE-vdW của Klimeš và cộng sự [20], phương pháp sóng tăng cường chiếu (Projector Augmented Wave - PAW) của Blöchl [27], và phương pháp hiệu chỉnh Hubbard $DFT+U$ theo công thức Dudarev [24].
Đóng góp đột phá của công trình mang tính định lượng cao: lần đầu tiên xác lập bản chất của cụm vòng hexa-hydroxyl tuần hoàn liên kết bởi mạng liên kết hydrogen chuyển dời đỏ $O-H\cdots O$ và chuyển dời xanh $O-H\cdots\pi$; tối ưu hóa thành công bộ tham số $(U_d = 7\text{ eV}, U_p = 10\text{ eV})$ nâng độ rộng vùng cấm rutile từ mức đánh giá thấp của GGA ($\sim 1.8\text{ eV}$) lên sát giá trị thực nghiệm $3.03\text{ eV}$ tại $10\text{ K}$; và thiết lập mô hình siêu ô hai chiều đồng phương vị giữa dẫn xuất graphene và rutile $TiO_2(110)$ với mức độ biến dạng mạng tối thiểu.
Phạm vi nghiên cứu bao quát từ đơn lớp graphene (ô đơn vị 2 nguyên tử và 4 nguyên tử), siêu ô $6 \times 3$ chứa 72 nguyên tử carbon cho dẫn xuất $G_nOH$, siêu ô $4 \times 3$ chứa 24 nguyên tử carbon cho các dẫn xuất $G_1F$, tinh thể rutile khối (nhóm không gian $P4_2/mnm$), và các mô hình bản mỏng bề mặt slab rutile (110) từ 3 đến 15 lớp nguyên tử.
Literature Review và Positioning
Tổng quan tài liệu cho thấy sự giao thoa phức tạp giữa bốn dòng nghiên cứu chính trong vật lý chất rắn và hóa học tính toán:
┌────────────────────────────────────────────────────────┐
│ LÝ THUYẾT PHIẾM HÀM MẬT ĐỘ (DFT) │
│ Hohenberg-Kohn (1964) & Kohn-Sham Equations (1965) │
└───────────────────────────┬────────────────────────────┘
│
┌────────────────────────┴────────────────────────┐
▼ ▼
┌─────────────────────────────────────┐ ┌─────────────────────────────────────┐
│ VẬT LIỆU 2D GRAPHENE & DẪN XUẤT │ │ BÁN DẪN OXIDE RUTILE TiO2 │
│ • Ô chữ nhật 4 C: Yamashita (2014), │ │ • Sai số SIE & DFT+U: │
│ Suzuki (2017), Martsinovich (2017)│ │ Anisimov (1991), Dudarev (1998) │
│ • Mô hình GO/RGO: Lerf-Dékány (1998)│ │ • Cấu trúc & năng lượng mặt (110): │
│ • Cụm -OH: Lahaye (2009), Yan (2014)│ │ Perron (2007), Mitsuhashi (1973) │
└──────────────────┬──────────────────┘ └──────────────────┬──────────────────┘
│ │
└────────────────────────┬────────────────────────┘
▼
┌────────────────────────────────────────────────────────┐
│ DỊ CẤU TRÚC COMPOSITE 2D/3D GRAPHENE/TiO2 │
│ Khắc phục sai lệch mạng & Tối ưu hóa chuyển điện tích │
│ (Luận án Tiến sĩ - Trần Thị Thoa, ĐHSPHN 2022) │
└────────────────────────────────────────────────────────┘
- Dòng nghiên cứu cấu trúc hình học và ô đơn vị của Graphene: K. Yamashita và cộng sự (2014) [63] là những người tiên phong đề xuất chuyển đổi ô đơn vị lục giác 2 nguyên tử sang ô đơn vị hình chữ nhật 4 nguyên tử để mô phỏng composite anatase $TiO_2(001)$/graphene. Tiếp đó, A. Suzuki và cộng sự (2017) [62] đã áp dụng giải tích thuyết liên kết chặt (tight-binding approach) chứng minh tính tương đương hoàn toàn về cấu trúc dải năng lượng và nón Dirac giữa hai loại ô mạng. Tuy nhiên, các nghiên cứu này chưa phân tích một cách hệ thống mật độ trạng thái riêng phần (PDOS) đối với các trạng thái orbital $2s, 2p_x, 2p_y, 2p_z$ dưới các điều kiện biên tuần hoàn khác nhau.
- Dòng nghiên cứu mô hình hóa Graphene Oxide (GO) và Reduced Graphene Oxide (RGO): Các mô hình kinh điển của Lerf-Klinowski [45] và Dékány [68] dựa trên phổ cộng hưởng từ hạt nhân trạng thái rắn ($^{13}\text{C}$ SSNMR) của Cai và cộng sự [46] khẳng định các nhóm hydroxyl và epoxy phân bố chủ yếu trên mặt đáy (basal plane) trong khi nhóm carboxyl (-COOH) và carbonyl ($=C=O$) định vị tại biên khuyết tật. Về mặt lý thuyết, Lahaye và cộng sự [73] chứng minh hai nhóm -OH không thể đính vào hai carbon liền kề cạnh nhau do lực đẩy lập thể cực lớn. Yan và cộng sự [75] chỉ ra xu hướng liên kết cặp ở vị trí para. Tuy nhiên, bức tranh tương tác khi số lượng nhóm -OH tăng lên ($n \ge 3$) hoàn toàn bị bỏ ngỏ.
- Dòng nghiên cứu vật lý bề mặt và cấu trúc điện tử của Rutile $TiO_2$: Rutile là pha đa hình bền vững nhất của $TiO_2$ với biến thiên enthalpy chuyển pha anatase-rutile đo bởi Mitsuhashi và cộng sự [50] là $-0.78 \pm 0.2\text{ kcal/mol}$. Tuy nhiên, tính toán DFT gặp phải hai tranh biện học thuật gay gắt:
- Tranh biện 1: Năng lượng bề mặt rutile (110): Dữ liệu thực nghiệm xác định năng lượng bề mặt nằm trong khoảng $0.28 - 0.38\text{ J/m}^2$ [87]. Ngược lại, tính toán gần đúng mật độ địa phương (LDA) cho giá trị sai lệch lên tới $\sim 0.89\text{ J/m}^2$ [88], trong khi gần đúng gradient tổng quát (GGA-PBE) dao động phân tán từ $0.47\text{ J/m}^2$ đến $0.74\text{ J/m}^2$ [90, 93].
- Tranh biện 2: Hiệu ứng phân kỳ giữa các mô hình bản mỏng (slab models): Nghiên cứu của H. Perron và cộng sự [95] ghi nhận sự chênh lệch năng lượng hội tụ lên tới $0.13 - 0.20\text{ J/m}^2$ giữa mô hình thả lỏng hoàn toàn (Full Relax - FR) và mô hình cố định các lớp bên trong (Fix Inner Layer - FIL) hoặc cố định hai lớp đáy (Fix Two Bottom - F2B).
- Dòng nghiên cứu bản chất liên kết Hydrogen phi cổ điển: Tranh biện giữa trường phái tương tác siêu liên hợp (Hyperconjugation charge transfer sang orbital phản liên kết $\sigma^*_{A-H}$) của Coulson [36] và trường phái tái tổ hợp lai hóa orbital kết hợp phân cực tĩnh điện cục bộ [37, 39] trong việc giải thích hiện tượng rút ngắn độ dài liên kết và dịch chuyển xanh tần số dao động (blue-shifting hydrogen bonds).
So sánh với hai nghiên cứu quốc tế tiêu biểu:
- So với công trình của N. Martsinovich và cộng sự (2017) [64] khảo sát dị diện graphene/rutile $TiO_2(110)$: Martsinovich sử dụng ô chữ nhật nhưng chỉ giới hạn ở graphene nguyên sinh, bỏ qua trạng thái biến tính bề mặt RGO và chưa áp dụng phiếm hàm van der Waals hiệu chỉnh tán xạ không định xứ optPBE-vdW.
- So với nghiên cứu của Lahaye và cộng sự (2009) [73]: Nhóm tác giả chỉ dừng lại ở cấu hình phân tán $n = 1, 2$ nhóm -OH mà không nhận diện được quy luật tập hợp vòng khép kín hexa-hydroxyl dẫn dắt bởi liên kết hydrogen kép.
Luận án của NCS. Trần Thị Thoa đã định vị chính xác tại điểm giao thoa này, giải quyết đồng thời bài toán cấu trúc vi mô của RGO và bài toán tiếp xúc pha bất đối xứng với rutile (110).
Đóng góp lý thuyết và khung phân tích
Đóng góp cho lý thuyết
Luận án mở rộng và làm sâu sắc thêm các hệ thống lý thuyết vật lý - hóa học chất rắn:
┌────────────────────────────────────────────────────────┐
│ HỆ THỐNG ĐÓNG GÓP LÝ THUYẾT CỦA LUẬN ÁN │
└───────────────────────────┬────────────────────────────┘
│
┌──────────────────────────────────────────────┴──────────────────────────────────────────────┐
▼ ▼
┌────────────────────────────────────────────────────────┐ ┌────────────────────────────────────────────────────────┐
│ LÝ THUYẾT LIÊN KẾT HYDROGEN TRÊN MẶT 2D │ │ LÝ THUYẾT TƯƠNG QUAN ĐIỆN TỬ MẠNH │
│ • Cơ chế hợp tác kép: Chuyển dời đỏ O-H···O và │ │ • Hiệu chỉnh đồng thời hai tham số Hubbard (Ud, Up): │
│ chuyển dời xanh O-H···π │ │ Ud = 7 eV (Ti-3d), Up = 10 eV (O-2p) │
│ • Quy luật hình thành vòng lục giác Hexa-hydroxyl │ │ • Triệt tiêu sai số tự tương tác SIE │
│ vị trí para bền vững nhiệt động │ │ • Dập tắt dao động chẵn-lẻ (even-odd) trên slab (110) │
└────────────────────────────────────────────────────────┘ └────────────────────────────────────────────────────────┘
- Mở rộng Thuyết Liên kết Hydrogen trên bề mặt vật liệu 2D: Minh chứng sự tồn tại đồng thời và cơ chế hiệp đồng (cooperative effect) giữa liên kết hydrogen chuyển dời đỏ cổ điển $O-H\cdots O$ và liên kết hydrogen chuyển dời xanh $O-H\cdots\pi$. Năng lượng liên kết $E_b$ tăng phi tuyến tính theo số lượng nhóm -OH khi các nhóm chức này kết tụ ở vị trí para cùng một phía mặt phẳng carbon, hình thành cấu trúc vòng khép kín hexa-hydroxyl bền vững.
- Phát triển Lý thuyết Hiệu chỉnh Tương quan Điện tử mạnh ($DFT+U$) cho Oxide Kim loại Chuyển tiếp: Chứng minh bằng chứng thực nghiệm số rằng độ rộng vùng cấm của rutile $TiO_2$ chịu sự kiểm soát đồng thời của trạng thái định xứ Ti-$3d$ tại đáy vùng dẫn (CBM) và trạng thái O-$2p$ tại đỉnh vùng hóa trị (VBM). Việc áp dụng mô hình Hubbard hai tham số $(U_d = 7\text{ eV}, U_p = 10\text{ eV})$ tạo nên bước chuyển dịch lý thuyết so với quan điểm truyền thống vốn chỉ hiệu chỉnh đơn lẻ trên orbital $d$ của cation kim loại.
- Lý thuyết Biến dạng Mạng và Dải Năng lượng Graphene: Chứng minh định lượng tính bảo toàn tính chất Dirac fermion của graphene khi chuyển đổi hệ tọa độ tinh thể học sang ô đơn vị hình chữ nhật bốn nguyên tử dưới tác động của trường thế giả PAW.
Khung phân tích độc đáo
Khung phân tích tích hợp bốn trụ cột phương pháp luận toán tử lượng tử:
$$\hat{H}{eff} = \hat{T}S + \hat{V}{ext} + \hat{V}H + \hat{V}{XC}^{optPBE-vdW} + \hat{V}{Hubbard}^{(U_d, U_p)}$$
┌──────────────────────────────────────────────────────────────┐
│ KHUNG PHÂN TÍCH TOÁN TỬ LƯỢNG TỬ │
│ Ĥeff = Ts + Vext + VH + Vxc + VHubbard │
└──────────────────────────────┬───────────────────────────────┘
│
┌───────────────────────────┬──────────────┴────────────┬────────────────────────────┐
▼ ▼ ▼ ▼
┌──────────────┐ ┌──────────────┐ ┌──────────────┐ ┌──────────────┐
│ Toán tử │ │ Thế tĩnh │ │ Phiếm hàm │ │ Số hạng │
│ Động năng │ + │ điện Hartree│ + │ Trao đổi │ + │ Hubbard │
│ không tương │ │ & Hạt nhân │ │ optPBE-vdW │ │ (Ud, Up) │
│ tác (Ts) │ │ (Vext + VH) │ │ (Vxc) │ │ (VHubbard) │
└──────────────┘ └──────────────┘ └──────────────┘ └──────────────┘
-
Semilocal & Nonlocal Exchange-Correlation Integration: Kết hợp phiếm hàm trao đổi optPBE cải tiến (tổ hợp $94.5268%$ PBE và $5.4732%$ revPBE với $\mu = 0$) cùng số hạng tương quan không định xứ $E_c^{nl}$ theo công thức Dion-Klimeš:
$$E_{XC}^{optPBE-vdW} = E_X^{optPBE} + E_C^{LDA} + E_C^{nl}$$
-
Mô hình Hubbard hiệu dụng Dudarev: Xử lý sai số SIE với thế hiệu dụng dạng đơn giản hóa:
$$E_{DFT+U} = E_{DFT} + \sum_{\sigma} \frac{U_{eff}}{2} \text{Tr} \left( \mathbf{n}^\sigma - \mathbf{n}^\sigma \mathbf{n}^\sigma \right)$$
-
Boundary Conditions và Vacuum Screening: Thiết lập điều kiện biên tuần hoàn 3 chiều (PBCs) với khoảng chân không (vacuum height) tối ưu hóa từ $12\text{ \AA}$ đến $15\text{ \AA}$, kết hợp hiệu chỉnh moment lưỡng cực (dipole correction) dọc trục $z$ nhằm triệt tiêu hoàn toàn tương tác tĩnh điện giả tạo giữa các lớp định kỳ bất đối xứng.
Phương pháp nghiên cứu tiên tiến
Thiết kế nghiên cứu
Nghiên cứu vận hành trên lập trường nhận thức luận thực chứng chặt chẽ (Positivist / Computational Material Realism), sử dụng phương pháp tính toán hóa học lượng tử ab initio từ nguyên lý đầu (first-principles calculations).
Thiết kế nghiên cứu đa mức (Multi-level Computational Design):
- Cấp độ 1 (Bulk & Primitive Cell): Khảo sát đơn lớp graphene nguyên sinh trên ô đơn vị lục giác 2 nguyên tử ($1 \times 1$) và ô đơn vị hình chữ nhật 4 nguyên tử ($1 \times 1$); tinh thể khối rutile $TiO_2$ đơn tinh thể ($1 \times 1 \times 1$).
- Cấp độ 2 (Supercells & Surface Slabs): Siêu ô $6 \times 3$ (72 nguyên tử C, kích thước hình học $14.82\text{ \AA} \times 12.84\text{ \AA}$) khảo sát dẫn xuất hydroxyl $G_nOH$; siêu ô $4 \times 3$ (24 nguyên tử C, kích thước $9.89\text{ \AA} \times 7.42\text{ \AA}$) khảo sát dẫn xuất đơn nhóm chức $G_1F$; bản mỏng đối xứng/bất đối xứng rutile (110) gồm 3, 5, 7, 9, 11, 13, 15 lớp nguyên tử.
- Cấp độ 3 (Composite Heterojunction): Ghép nối siêu ô trực giao hai chiều $G/TiO_2(110)$ và $GOH/TiO_2(110)$.
[CẤP ĐỘ 1: CẤU TRÚC CƠ BẢN]
├─ Graphene ô lục giác (2C) & ô chữ nhật (4C)
└─ Rutile TiO2 khối đơn tinh thể (P42/mnm)
│
▼
[CẤP ĐỘ 2: SIÊU Ô & BẢN MỎNG BỀ MẶT]
├─ Siêu ô GnOH (6x3, 72 nguyên tử C, 14.82 x 12.84 Å)
├─ Siêu ô G1F (4x3, 24 nguyên tử C, 9.89 x 7.42 Å)
└─ Bản mỏng Rutile (110) slab models (FR, FIL, F2B từ 3-15 lớp)
│
▼
[CẤP ĐỘ 3: DỊ CẤU TRÚC COMPOSITE 2D/3D]
├─ Graphene / Rutile TiO2 (110)
└─ GOH / Rutile TiO2 (110)
Quy trình nghiên cứu rigorous
Toàn bộ các phép tính toán được thực thi thông qua gói phần mềm cơ lượng tử chuyên dụng VASP (Vienna Ab initio Simulation Package) [97].
- Cơ sở sóng phẳng và Thế giả: Tương tác giữa lõi ion và electron hóa trị được mô tả bằng phương pháp sóng tăng cường chiếu PAW [28, 98]. Cấu hình electron hóa trị được xử lý tường minh gồm: $\text{C}(2s^2 2p^2)$, $\text{H}(1s^1)$, $\text{O}(2s^2 2p^4)$, $\text{N}(2s^2 2p^3)$, và $\text{Ti}(3d^2 4s^2)$.
- Tiêu chuẩn hội tụ (Convergence Thresholds):
- Năng lượng cắt sóng phẳng (Plane-wave cut-off energy): $\text{ENCUT} = 450\text{ eV}$.
- Tiêu chí hội tụ trường tự hợp (SCF energy convergence): $\Delta E \le 10^{-6}\text{ eV}$.
- Tối ưu hóa cấu trúc hình học (Geometry optimization): Lực Hellmann-Feynman tác dụng lên từng nguyên tử đạt ngưỡng $|\mathbf{F}| < 0.01\text{ eV/\AA}$.
- Lấy mẫu vùng Brillouin (k-point sampling): Sử dụng lưới điểm Monkhorst-Pack [99]:
- Ô chữ nhật graphene 4 nguyên tử: $12 \times 6 \times 1$ (tối ưu hóa), $36 \times 18 \times 1$ (tính DOS).
- Ô lục giác graphene 2 nguyên tử: $6 \times 6 \times 1$ (tối ưu hóa), $36 \times 36 \times 1$ (tính DOS).
- Siêu ô dẫn xuất $G_nOH$ ($6 \times 3$): $2 \times 2 \times 1$ (tối ưu hóa), $12 \times 12 \times 1$ (tính DOS).
- Siêu ô dẫn xuất $G_1F$ ($4 \times 3$): $3 \times 2 \times 1$ (tối ưu hóa), $18 \times 12 \times 1$ (tính DOS).
- Tính toán phân cực spin (spin-polarized calculations) được áp dụng bắt buộc cho tất cả các hệ dẫn xuất chứa số electron lẻ hoặc có biến dạng mật độ spin cục bộ.
Data và phân tích
Năng lượng liên kết hình thành ($E_b$) của dẫn xuất hydroxyl $G_nOH$ được định nghĩa chuẩn xác theo công thức nhiệt động học:
$$E_b = \frac{E_{G_nOH} - (E_G + n E_{OH})}{n} \quad (\text{eV/OH})$$
Trong đó $E_{G_nOH}$, $E_G$, và $E_{OH}$ lần lượt là năng lượng tổng đã tối ưu hóa của dẫn xuất $G_nOH$, graphene nguyên sinh và một gốc tự do hydroxyl cô lập.
Năng lượng bề mặt của bản mỏng slab rutile (110) được xác định thông qua:
$$\sigma = \frac{E_{slab} - N \cdot E_{bulk}}{2A}$$
Với $E_{slab}$ là tổng năng lượng của slab, $E_{bulk}$ là năng lượng trên một đơn vị công thức $TiO_2$ trong tinh thể khối, $N$ là số đơn vị công thức trong slab, và $A$ là diện tích bề mặt tiếp xúc.
Độ tin cậy của mô hình được đối chiếu trực tiếp qua phổ nhiễu xạ tia X (XRD) tính toán từ định luật Bragg $n\lambda = 2d \sin\theta$, kiểm chứng vị trí các píc nhiễu xạ chính giữa thực nghiệm và mô phỏng $DFT/DFT+U$.
Phát hiện đột phá và implications
Những phát hiện then chốt
┌────────────────────────────────────────────────────────┐
│ 5 PHÁT HIỆN ĐỘT PHÁ CỦA LUẬN ÁN │
└───────────────────────────┬────────────────────────────┘
│
┌───────────────────────────┬──────────────────────┴──────────────┬────────────────────────────┐
▼ ▼ ▼ ▼
┌────────────────┐ ┌────────────────┐ ┌────────────────┐ ┌────────────────┐
│ PHÁT HIỆN 1 │ │ PHÁT HIỆN 2 │ │ PHÁT HIỆN 3 │ │ PHÁT HIỆN 4 │
│ Vòng lục giác │ │ Trật tự độ bền │ │ Cặp tham số │ │ Triệt tiêu dao │
│ Hexa-hydroxyl │ │ G1F: │ │ (Ud=7, Up=10)eV│ │ động chẵn-lẻ │
│ para bền vững │ │ Gepo > GCH3 > │ │ khôi phục Eg │ │ trên slab (110)│
│ tuyệt đối │ │ GNH2 >> GCOOH │ │ = 3.03 eV │ │ bằng DFT+U │
└────────────────┘ └────────────────┘ └────────────────┘ └────────────────┘
-
Quy luật tập hợp vòng Hexa-hydroxyl para ($G_6OH$): Dữ liệu tính toán chứng minh rằng trong hệ $G_nOH$, các nhóm -OH không phân bố ngẫu nhiên mà có xu hướng hội tụ tại các vị trí para trên cùng một phía của mặt phẳng graphene tạo nên cấu hình vòng lục giác hoàn hảo. Trích dẫn nguyên văn kết luận từ luận án:
"Trong dẫn xuất $G_nOH$, các nhóm hydroxyl có xu hướng tập hợp tại các vị trí para trên cùng một phía của mặt phẳng graphene để tạo ra các vòng hexa hydroxyl hoàn hảo. Đây là kết quả của sự tạo thành các liên kết hydrogen chuyển dời đỏ $O-H\cdots O$ và các liên kết hydrogen chuyển dời xanh $O-H\cdots\pi$ giữa các nhóm hydroxyl." Năng lượng liên kết $E_b$ chuyển từ $-0.62\text{ eV/OH}$ ở $G_1OH$ xuống mức cực tiểu bền vững nhiệt động $-1.45\text{ eV/OH}$ ở cấu hình $G_6OH$, đi kèm sự biến dạng góc liên kết C-C-C cục bộ từ lai hóa phẳng $sp^2$ ($120^\circ$) sang lai hóa tứ diện $sp^3$ ($\sim 109.5^\circ$).
-
Trật tự bền vững nhiệt động của các dẫn xuất đơn nhóm chức $G_1F$: Trật tự năng lượng liên kết được thiết lập chính xác:
$$G_{epo} \text{ (bền nhất)} > G_{CH_3} > G_{NH_2} > G_{OCH_3} > G_{CHO} \gg G_{COOH}$$
Dẫn xuất chứa nhóm carboxyl ($G_{COOH}$) có độ bền liên kết kém nhất ($E_b$ xấp xỉ $-0.18\text{ eV}$), giải thích hoàn hảo cho các quan sát thực nghiệm phổ $^{13}\text{C}$ SSNMR của Cai và cộng sự [46] khi nhóm -COOH chỉ có thể tồn tại tại các mép khuyết tật có năng lượng tự do biên cao thay vì trên mặt đáy phẳng. Mật độ trạng thái (DOS) chỉ ra rằng các mức năng lượng xung quanh mức Fermi ($E_F = 0\text{ eV}$) được chi phối áp đảo bởi orbital $2p_z$ của dị nguyên tử liên kết trực tiếp với carbon graphene.
-
Đột phá trong hiệu chỉnh Hubbard đồng thời $(U_d = 7\text{ eV}, U_p = 10\text{ eV})$ cho Rutile $TiO_2$: Luận án khẳng định vai trò quyết định của việc hiệu chỉnh orbital $2p$ của anion oxy. Trích dẫn số liệu phát hiện:
"Đối với rutile $TiO_2$: $U_p = 10\text{ eV}$ là trị số tối ưu. Sự kết hợp của $U_d$ và $U_p$ tối ưu ($U_p = 10\text{ eV}$) không chỉ làm giảm sai số trong hằng số mạng mà còn làm tăng $E_g$ nhiều hơn so với sự hiệu chỉnh của $U_d$." Độ rộng vùng cấm $E_g$ đạt giá trị $3.03\text{ eV}$ trùng khớp hoàn toàn với thực nghiệm tại $10\text{ K}$ [52], trong khi các hằng số mạng duy trì độ chính xác cao: $a = 4.602\text{ \AA}$ (sai lệch $+0.26%$), $c = 2.961\text{ \AA}$ (sai lệch $+0.37%$), vượt trội hoàn toàn so với việc tăng $U_d$ lên $10\text{ eV}$ đơn lẻ gây giãn nở mạng nghiêm trọng ($a = 4.68\text{ \AA}$, sai lệch $> 1.9%$).
-
Dập tắt dao động chẵn-lẻ (Even-Odd Oscillation) trên bề mặt Rutile (110): Trên mô hình bản mỏng thả lỏng hoàn toàn (FR), năng lượng bề mặt và độ rộng dải cấm của các slab từ 3 đến 15 lớp biểu hiện sự dao động chẵn-lẻ rõ rệt khi tính bằng DFT thuần túy. Phương pháp $DFT+U$ với $U_p = 10\text{ eV}$ đã dập tắt hoàn toàn sự bất ổn định này, đưa năng lượng bề mặt hội tụ mượt mà về khoảng $0.48\text{ J/m}^2$. Luận án cũng chỉ ra mô hình cố định hai lớp đáy (F2B) dù hội tụ số lớp nhanh nhưng gây ra sai số dải cấm bề mặt lớn nhất, khẳng định mô hình FR là chuẩn mực cho các tính toán tiếp xúc dị diện.
-
Thiết lập siêu ô trực giao cho Composite $G/TiO_2(110)$ và $GOH/TiO_2(110)$: Khớp nối thành công ô đơn vị hình chữ nhật 4 nguyên tử carbon của graphene ($a_3 = 2.463\text{ \AA}, a_4 = 4.266\text{ \AA}$) với ô đơn vị chữ nhật bề mặt rutile (110) ($a = 2.959\text{ \AA}, b = 6.496\text{ \AA}$). Bằng cách xây dựng siêu mạng tỷ lệ $m \times n$ phù hợp, độ sai lệch mạng dọc hai trục trực giao được khống chế dưới mức $2.1%$, bảo toàn cấu trúc điện tử nội tại của cả hai hợp phần.
Implications đa chiều
- Về mặt lý thuyết hóa học: Cung cấp bức tranh toàn cảnh về tương tác liên kết hydrogen trên các bề mặt nano hai chiều, thiết lập cầu nối giữa hóa học lượng tử tính toán và hóa học cấu trúc siêu phân tử.
- Về mặt phương pháp luận: Thiết lập quy chuẩn tính toán cho các hệ bán dẫn oxide kim loại chuyển tiếp tương quan mạnh thông qua việc kết hợp phiếm hàm phân tán vdW với hiệu chỉnh Hubbard hai kênh $(U_d, U_p)$.
- Về mặt ứng dụng thực tiễn: Nhóm chức hydroxyl và epoxy đóng vai trò là "cầu nối điện tử" (electronic bridges) xúc tiến quá trình phân tách điện tích quang sinh, giảm thiểu tốc độ tái tổ hợp electron - lỗ trống ($e^-/h^+$), mở đường cho việc thiết kế vật liệu xúc tác quang phân hủy chất độc dệt nhuộm và sản xuất khí hydro đạt hiệu suất lượng tử cao.
Limitations và Future Research
Luận án thẳng thắn thừa nhận các giới hạn kỹ thuật và điều kiện biên:
- Giới hạn nhiệt độ không tuyệt đối ($0\text{ K}$): Tất cả các tính toán DFT được thực hiện tại $0\text{ K}$ trong môi trường chân không, chưa tính đến các hiệu ứng kích thích nhiệt động học và dao động mạng phonon ở nhiệt độ phòng ($298\text{ K}$).
- Bỏ qua hiệu ứng dung môi tường minh (Explicit Solvent Effects): Các phản ứng xúc tác quang hóa thực tế diễn ra trong pha lỏng (dung dịch nước), nơi các phân tử nước tự do có thể cạnh tranh liên kết hydrogen với nhóm -OH trên bề mặt graphene.
- Mô hình hóa khuyết tật tinh thể (Defect Modeling): Nghiên cứu giả định bề mặt rutile (110) là hoàn hảo, chưa xét đến các khuyết tật điểm thực tế như oxy trống ($V_O$) hay các tâm $Ti^{3+}$ xen kẽ vốn ảnh hưởng mạnh đến vị trí mức Fermi.
- Mật độ che phủ nhóm chức hữu hạn: Nghiên cứu giới hạn ở các tỷ lệ che phủ $G_nOH$ từ $n=1$ đến $n=6$ trên siêu ô $6 \times 3$, chưa khảo sát nồng độ che phủ bão hòa cao ($> 50%$).
Chương trình nghiên cứu tiếp nối (Future Research Agenda):
- Hướng 1: Ứng dụng mô phỏng động lực học phân tử từ nguyên lý đầu (Ab Initio Molecular Dynamics - AIMD) ở điều kiện nhiệt độ hữu hạn ($300\text{ K} - 500\text{ K}$) để đánh giá độ bền nhiệt động của vòng hexa-hydroxyl.
- Hướng 2: Tích hợp mô hình dung môi liên tục (COSMO/PCM) và các phân tử nước tường minh khảo sát quá trình hấp phụ chất ô nhiễm hữu cơ.
- Hướng 3: Tính toán phổ kích thích quang học và mức năng lượng chuẩn xác thông qua lý thuyết nhiễu loạn nhiều hạt (Many-body Perturbation Theory - phương pháp $GW$ và phương trình Bethe-Salpeter BSE).
- Hướng 4: Khảo sát động học truyền điện tích phi đoạn nhiệt (Non-adiabatic Molecular Dynamics - NAMD) qua bề mặt tiếp xúc $GOH/TiO_2(110)$.
Tác động và ảnh hưởng
┌────────────────────────────────────────────────────────┐
│ HỆ THỐNG TÁC ĐỘNG VÀ ẢNH HƯỞNG │
└───────────────────────────┬────────────────────────────┘
│
┌───────────────────────────┬──────────────────────┴──────────────┬────────────────────────────┐
▼ ▼ ▼ ▼
┌────────────────┐ ┌────────────────┐ ┌────────────────┐ ┌────────────────┐
│ HỌC THUẬT │ │ CÔNG NGHIỆP │ │ MÔI TRƯỜNG │ │ QUỐC TẾ │
│ Tiêu chuẩn │ │ Chế tạo màng │ │ Xử lý triệt để │ │ Khắc phục hạn │
│ mô hình hóa │ │ quang xúc tác, │ │ chất thải hữu │ │ chế mô hình │
│ RGO & Rutile │ │ photoanode │ │ cơ độc hại │ │ của Long, Du │
└────────────────┘ └────────────────┘ └────────────────┘ └────────────────┘
- Tác động học thuật (Academic Impact): Cung cấp bộ thông số chuẩn hóa ($U_d=7\text{ eV}, U_p=10\text{ eV}$, mô hình ô chữ nhật 4 nguyên tử) cho cộng đồng nghiên cứu mô phỏng vật liệu nano. Dự kiến tạo lập nguồn trích dẫn vững chắc trong các nghiên cứu chuyên sâu về giao diện graphene/oxide kim loại.
- Chuyển dịch công nghệ công nghiệp (Industry Transformation): Định hướng cho các kỹ sư R&D trong việc điều khiển quá trình khử GO bằng nhiệt hoặc hóa học nhằm tối ưu hóa mật độ nhóm chức hydroxyl vị trí para, nâng cao hiệu suất dòng quang điện trong pin mặt trời DSSC và tăng dung lượng chu chuyển của pin lithium-ion.
- Lợi ích kinh tế - xã hội và Môi trường: Đóng góp trực tiếp vào mục tiêu trung hòa carbon và chuyển dịch năng lượng xanh; nâng cao hiệu quả quang phân hủy các chất độc khó phân hủy sinh học (thuốc bảo vệ thực vật, phụ gia dệt nhuộm, hợp chất hữu cơ vòng thơm) trong nước thải công nghiệp.
- Tầm vóc quốc tế (International Relevance): Khắc phục triệt để khiếm khuyết về sai lệch mạng trong các công bố quốc tế trước đó (như các mô hình của Long et al. hay Du et al.), xác lập vị thế tiên phong của nghiên cứu tính toán vật liệu tại Việt Nam trên diễn đàn khoa học thế giới.
Đối tượng hưởng lợi
- Nghiên cứu sinh & Nhà nghiên cứu trẻ (Doctoral Researchers): Tiếp cận quy trình mô phỏng chuẩn mực từ khâu xây dựng ô đơn vị, siêu ô, kiểm soát hội tụ k-point/ENCUT, đến kỹ thuật xử lý phiếm hàm van der Waals và hiệu chỉnh $DFT+U$.
- Các chuyên gia và Học giả đầu ngành (Senior Academics): Sở hữu luận cứ lý thuyết vững chắc về liên kết hydrogen phi cổ điển trên vật liệu carbon 2D và cơ chế đóng góp của orbital anion O-$2p$ trong vật lý bán dẫn oxide.
- Kỹ sư R&D Công nghệ Vật liệu (Industry R&D Specialists): Nhận được các thông số định lượng chính xác về độ bền liên kết và dải năng lượng để thiết kế quy trình chế tạo vật liệu nano composite có hoạt tính quang xúc tác cao.
- Các nhà hoạch định chính sách Năng lượng & Môi trường: Có cơ sở khoa học để phê duyệt và đầu tư vào các dự án ứng dụng công nghệ nano trong xử lý ô nhiễm nước và phát triển nguồn năng lượng hydro bền vững.
Câu hỏi chuyên sâu
1. Đóng góp lý thuyết độc đáo nhất của luận án là gì và đã mở rộng lý thuyết nào?
Đóng góp lý thuyết độc đáo nhất là phát hiện ra quy luật tự tổ chức tạo vòng khép kín hexa-hydroxyl vị trí para ($G_6OH$) trên mặt phẳng graphene, được duy trì bởi cơ chế hiệp đồng kép giữa liên kết hydrogen chuyển dời đỏ $O-H\cdots O$ và liên kết chuyển dời xanh $O-H\cdots\pi$. Phát hiện này mở rộng trực tiếp Lý thuyết mô hình cấu trúc Graphene Oxide của Lerf-Klinowski (1998) [45] và Dékány (1998) [68], chuyển đổi trạng thái mô tả phân bố nhóm chức từ định tính/ngẫu nhiên sang mô hình có trật tự hình học lượng tử hóa chính xác.
2. Điểm cải tiến phương pháp luận so với ít nhất 02 công trình quốc tế tiền nhiệm?
- So với Long và cộng sự (2012) [56]: Luận án đã thay thế ô đơn vị lục giác 2 nguyên tử của graphene bằng ô đơn vị hình chữ nhật 4 nguyên tử ($a_3 = 2.463\text{ \AA}, a_4 = 4.266\text{ \AA}$), triệt tiêu hoàn toàn sai lệch góc mạng $90^\circ / 120^\circ$ khi ghép với rutile $TiO_2(110)$, loại bỏ ứng suất mạng nhân tạo.
- So với H. Perron và cộng sự (2007) [95]: Thay vì chỉ khảo sát các phiếm hàm LDA/GGA tiêu chuẩn dẫn đến sự dao động chẵn-lẻ và sai lệch dải cấm lớn, luận án đã tích hợp phương pháp $DFT+U$ hai tham số $(U_d=7\text{ eV}, U_p=10\text{ eV})$, dập tắt sự dao động phân kỳ theo số lớp của mô hình slab thả lỏng hoàn toàn (FR).
3. Phát hiện nào gây bất ngờ nhất về mặt số liệu thực nghiệm tính toán?
Phát hiện bất ngờ nhất là vai trò chi phối vượt trội của tham số Hubbard $U_p = 10\text{ eV}$ trên orbital $2p$ của oxy so với tham số $U_d$ trên orbital $3d$ của titan. Trước đây, hầu hết các nghiên cứu chỉ tập trung hiệu chỉnh $U_d$ (thường chọn $U_d \approx 10\text{ eV}$) dẫn đến hệ quả tiêu cực là làm biến dạng mạng tinh thể khối ($a = 4.68\text{ \AA}$). Luận án chứng minh $U_p = 10\text{ eV}$ kết hợp $U_d = 7\text{ eV}$ mở rộng dải cấm rutile đạt đúng $3.03\text{ eV}$ thực nghiệm mà vẫn giữ nguyên hằng số mạng chính xác ($a = 4.602\text{ \AA}$).
4. Luận án có cung cấp đầy đủ giao thức tái lập (Replication Protocol) không?
Có, giao thức tái lập được cung cấp cực kỳ chi tiết:
- Phần mềm: VASP; Thế giả: PAW; Năng lượng cắt: $\text{ENCUT} = 450\text{ eV}$.
- Tiêu chí hội tụ: Lực Hellmann-Feynman $< 0.01\text{ eV/\AA}$, Năng lượng SCF $10^{-6}\text{ eV}$.
- Không gian chân không: $12\text{ \AA}$ (cho $G_nOH$), $15\text{ \AA}$ (cho $G_1F$) kèm hiệu chỉnh moment lưỡng cực.
- Lưới k-points Monkhorst-Pack chuẩn hóa cho từng kích thước siêu ô và quy trình tối ưu hóa hình học.
5. Lộ trình nghiên cứu 10 năm tới được phác thảo ra sao?
Lộ trình 10 năm hướng tới việc chuyển dịch từ mô hình hóa cụm tĩnh tại $0\text{ K}$ sang động học truyền điện tích thời gian thực trong môi trường điện ly hoạt động:
- Giai đoạn 1 (1-3 năm): Động lực học phân tử lượng tử (AIMD) nhiệt độ phòng và dung môi tường minh.
- Giai đoạn 2 (4-6 năm): Khảo sát động học quang kích thích bằng lý thuyết $GW$-BSE và NAMD.
- Giai đoạn 3 (7-10 năm): Thiết kế giao diện composite đa pha tích hợp công nghệ AI/Machine Learning thế điện tử (Machine Learning Potentials) phục vụ sản xuất thiết bị quang xúc tác công nghiệp.
Kết luận
- Xác lập quy luật tập hợp cấu trúc vi mô của dẫn xuất Hydroxyl Graphene ($G_nOH$): Chứng minh các nhóm -OH liên kết ưu tiên tại các vị trí para cùng một phía mặt phẳng carbon, hình thành các vòng hexa-hydroxyl bền vững nhiệt động ($E_b = -1.45\text{ eV/OH}$) nhờ mạng lưới liên kết hydrogen kép $O-H\cdots O$ và $O-H\cdots\pi$.
- Xác định trật tự bền vững nhiệt động và cấu trúc electron của các dẫn xuất đơn nhóm chức ($G_1F$): Xác lập thang độ bền $G_{epo} > G_{CH_3} > G_{NH_2} > G_{OCH_3} > G_{CHO} \gg G_{COOH}$, làm sáng tỏ bản chất các trạng thái năng lượng vùng biên bị chi phối bởi orbital $2p_z$ của dị nguyên tử.
- Chuẩn hóa phương pháp $DFT+U$ hai tham số cho Rutile $TiO_2$: Xác định cặp tham số tối ưu $(U_d = 7\text{ eV}, U_p = 10\text{ eV})$, giải quyết triệt để sai số tự tương tác SIE, khôi phục độ rộng vùng cấm $E_g = 3.03\text{ eV}$ chuẩn thực nghiệm và bảo toàn hằng số mạng tinh thể.
- Chuẩn hóa tính toán bề mặt Rutile (110): Dập tắt hiện tượng dao động chẵn-lẻ của năng lượng bề mặt và độ rộng dải cấm theo số lớp slab thông qua $DFT+U$, xác định mô hình thả lỏng hoàn toàn (FR) là mô hình phản ánh trung thực nhất bản chất vật lý bề mặt.
- Thiết lập thành công mô hình tiếp xúc dị cấu trúc 2D Graphene/Rutile $TiO_2(110)$: Khắc phục hoàn toàn sự bất tương thích hình học mạng tinh thể bằng ô đơn vị hình chữ nhật bốn nguyên tử của graphene, mở ra nền tảng mô phỏng chính xác cho các hệ vật liệu quang xúc tác và quang điện nano thế hệ mới.