Tổng quan luận án

Vật liệu composite giữa graphene và titan dioxide ($\text{TiO}_2$) thu hút nhiều sự quan tâm trong lĩnh vực quang xúc tác và quang điện nhờ khả năng xúc tác phân hủy các chất hữu cơ độc hại trong nước thải (từ các nhà máy chế biến nhựa, chất tẩy rửa, thuốc trừ sâu, phụ gia thực phẩm), quang phân hủy nước để sản xuất nhiên liệu sạch, ứng dụng trong pin Mặt Trời chất màu nhạy quang (DSSC), anode pin lithium, cảm biến khí và màng phủ bảo vệ.

Khoảng trống nghiên cứu được tác giả xác định xuất phát từ ba vấn đề chính trong các công bố khoa học trước đó:

  1. Về mô hình composite với rutile $\text{TiO}_2$: Các nghiên cứu lý thuyết trước đây chủ yếu tập trung vào composite của graphene với bề mặt anatase $\text{TiO}_2$ (101) do sự tương đồng về ô đơn vị dạng lục giác. Ngược lại, bề mặt rutile $\text{TiO}_2$ (110) có ô đơn vị hình chữ nhật, khiến việc ghép nối với ô đơn vị lục giác truyền thống của graphene gây ra sai lệch mạng (lattice mismatch) rất lớn về góc ($90^\circ$ so với $120^\circ$). Một số nghiên cứu trước đó vẫn gượng ép ghép ô lục giác của graphene lên mặt rutile (110) hoặc chỉ khảo sát các cụm cluster $\text{TiO}_2$ phân tán rời rạc trên bề mặt graphene.
  2. Về cấu trúc của dẫn xuất graphene: Trong thực nghiệm điều chế, sản phẩm thu được thường là graphene oxide (GO) hoặc dạng khử graphene oxide (RGO) với các nhóm chức chính trên mặt đáy là hydroxyl (-OH) và epoxy (-O-). Các nghiên cứu lý thuyết trước đây chủ yếu gắn ngẫu nhiên các nhóm chức hoặc chỉ dừng lại ở việc khảo sát một đến hai nhóm -OH trên cùng một phía của mặt phẳng graphene, chưa giải thích được quy luật phân bố không gian khi số lượng nhóm -OH tăng lên.
  3. Về phương pháp tính toán hệ rutile $\text{TiO}_2$: Phương pháp phiếm hàm mật độ (DFT) chuẩn với các phép gần đúng LDA/GGA thường đánh giá độ rộng vùng cấm ($E_g$) của rutile $\text{TiO}_2$ thấp hơn nhiều so với thực nghiệm do sai số tự tương tác (SIE). Việc áp dụng hiệu chỉnh DFT+U trước đây chủ yếu thực hiện trên orbital $3d$ của Ti ($U_d$) với các giá trị phân tán và chưa khảo sát đầy đủ sự hiệu chỉnh đồng thời trên orbital $2p$ của O ($U_p$).

Mục tiêu nghiên cứu của luận án được xác định gồm 4 nội dung:

  1. Tìm cấu trúc và tính chất electron của các dẫn xuất graphene với hydroxyl, ký hiệu là $\text{G}_n\text{OH}$ ($n \in \mathbb{N}^*$).
  2. Tìm cấu trúc và tính chất electron của các dẫn xuất graphene với một nhóm chức ($\text{G}_1\text{F}$), gồm các nhóm: $-\text{NH}_2$, $-\text{CH}_3$, $-\text{OCH}_3$, $-\text{CHO}$, $-\text{COOH}$ và epoxy.
  3. Xác định cấu trúc, tính chất electron của tinh thể rutile $\text{TiO}_2$ và bề mặt rutile $\text{TiO}_2$ (110).
  4. Xây dựng mô hình tính toán hai chiều cho composite giữa graphene với bề mặt rutile $\text{TiO}_2$ (110), và composite giữa dẫn xuất hydroxyl graphene với bề mặt rutile $\text{TiO}_2$ (110).

Đối tượng nghiên cứu của luận án là cấu trúc hình học, năng lượng liên kết và các tính chất electron (mật độ trạng thái DOS, mật độ trạng thái riêng phần PDOS, cấu trúc dải năng lượng Band structure, sự chênh lệch mật độ electron EDD) của graphene, các dẫn xuất chức hóa của graphene, tinh thể khối rutile $\text{TiO}_2$, bề mặt cắt rutile $\text{TiO}_2$ (110) và các cấu trúc composite hai chiều hình thành giữa chúng.

Phạm vi nghiên cứu được giới hạn trong khuôn khổ hóa học tính toán và vật lý chất rắn lượng tử, sử dụng lý thuyết phiếm hàm mật độ kết hợp các phép hiệu chỉnh tương tác phân tán van der Waals (optPBE-vdW) và hiệu chỉnh tương tác mạnh Hubbard (DFT+U), được thực hiện và hoàn thành tại Trường Đại học Sư phạm Hà Nội vào năm 2022.

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

Luận án đã tổng quan các hướng nghiên cứu lý thuyết và thực nghiệm liên quan đến hệ vật liệu nghiên cứu:

  • Nghiên cứu về graphene và ô đơn vị hình chữ nhật: K. Yamashita và cộng sự (2014) là những người đầu tiên đề xuất chuyển ô đơn vị lục giác 2 nguyên tử của graphene thành ô đơn vị hình chữ nhật 4 nguyên tử khi nghiên cứu composite anatase $\text{TiO}_2$ (001)/graphene để giảm sai lệch mạng. Năm 2017, N. Martsinovich và cộng sự sử dụng ô đơn vị hình chữ nhật để nghiên cứu cấu trúc electron và sự truyền điện tích của graphene/rutile $\text{TiO}_2$ (110), chỉ ra ảnh hưởng của ứng suất biến dạng. Cùng năm 2017, A. Suzuki và cộng sự áp dụng phương pháp liên kết chặt (tight-binding) cho ô đơn vị hình chữ nhật 4 nguyên tử và thu được kết quả tương đương ô lục giác.
  • Nghiên cứu về dẫn xuất graphene oxide (GO, RGO): Các công trình thực nghiệm của Weiwei Cai và cộng sự (2008), Lerf (1998), Dékány (1998) và các phân tích phổ cộng hưởng từ hạt nhân trạng thái rắn ($^{13}\text{C}$ SSNMR) chỉ ra sự tồn tại chủ yếu của nhóm epoxy (60 ppm), nhóm C-OH (70 ppm) và carbon lai hóa $sp^2$ (130 ppm) trên mặt đáy graphene. Về mặt tính toán, Lahaye và cộng sự (2009) chỉ ra hai nhóm -OH cạnh nhau trên cùng một phía rất kém bền; Yan và cộng sự (2011) nhận thấy hai nhóm -OH có xu hướng định vị ở vị trí para. Tuy nhiên, cấu hình sắp xếp của nhiều hơn 2 nhóm -OH chưa được giải quyết.
  • Nghiên cứu về rutile $\text{TiO}_2$ và bề mặt (110): Nghiên cứu nhiệt động của Mitsuhashi và cộng sự xác định biến thiên enthalpy chuyển pha anatase-rutile là $-0{,}78 \pm 0{,}2\text{ kcal/mol}$, khẳng định rutile là pha bền nhiệt động nhất. Đối với phương pháp tính toán DFT+U, các tác giả trước đây khảo sát tham số $U_d$ cho orbital $3d$ của Ti trong khoảng rộng; giá trị $U_d \approx 10\text{ eV}$ cho độ rộng vùng cấm sát thực nghiệm nhưng gây biến dạng lớn hằng số mạng.
  • Nghiên cứu về liên kết hydrogen mở rộng: Khái niệm liên kết hydrogen ban đầu của L. Pauling (1931) và bản chất chuyển điện tích của Coulson (1957) được bổ sung bởi các phát hiện về liên kết $\text{C}-\text{H}\cdots\text{O}$ của J. Sutor (1962), liên kết $\text{C}-\text{H}\cdots\text{Anion}$ của Pedzisa và cộng sự, liên kết $\text{A}-\text{H}\cdots\pi$, và liên kết dihydrogen $\text{A}-\text{H}\cdots\text{H}-\text{M}$ của Jackson và cộng sự (2005).

Khoảng trống học thuật mà luận án lựa chọn giải quyết là:

  1. Thiết lập quy luật sắp xếp không gian và bản chất liên kết hydrogen nội lớp của các dẫn xuất hydroxyl graphene $\text{G}_n\text{OH}$ với số lượng nhóm chức tăng dần ($n = 1 \div 6$).
  2. Xác định ảnh hưởng của các nhóm chức hữu cơ khác nhau ($\text{G}_1\text{F}$) lên cấu trúc dải và mật độ trạng thái của graphene.
  3. Tối ưu hóa đồng thời các tham số Hubbard ($U_d, U_p$) cho rutile $\text{TiO}_2$ để vừa mô tả chính xác độ rộng vùng cấm vừa bảo đảm độ chính xác của hằng số mạng tinh thể.
  4. Xây dựng mô hình ô mạng giao diện hai chiều tương thích giữa ô đơn vị hình chữ nhật của dẫn xuất graphene và bề mặt rutile $\text{TiO}_2$ (110).

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

Luận án sử dụng nền tảng hóa học lượng tử và vật lý chất rắn tính toán:

  • Lý thuyết phiếm hàm mật độ (DFT): Dựa trên hai định lý nền tảng của Hohenberg-Kohn (1964) về tính duy nhất của phiếm hàm mật độ electron trạng thái cơ bản và nguyên lý biến phân năng lượng, kết hợp phương trình Kohn-Sham (1965) đưa vào hệ giả tưởng các electron không tương tác di chuyển trong trường thế hiệu dụng $V_S$.
  • Phiếm hàm tương quan - trao đổi và hiệu chỉnh tán xạ: Sử dụng phép gần đúng gradient tổng quát GGA với phiếm hàm PBE (Perdew, Burke, Ernzerhof, 1996) chứa hai tham số tăng cường trao đổi cục bộ $\kappa = 0{,}804$ và gradient tương quan hiệu dụng $\mu = 0{,}21951$. Để mô tả tương tác yếu, luận án áp dụng phiếm hàm van der Waals không định xứ optPBE-vdW của Klimeš và cộng sự (kết hợp $94{,}5268%$ PBE và $5{,}4732%$ revPBE với $\mu = 0$).
  • Phương pháp DFT+U: Sử dụng phương pháp hiệu chỉnh sai số tự tương tác SIE theo mô hình Dudarev với tham số hiệu dụng $U_{\text{eff}} = U - J$, áp dụng độc lập và đồng thời cho orbital Ti-$3d$ ($U_d$) và O-$2p$ ($U_p$).
  • Bộ hàm cơ sở và thế giả: Mô tả hệ tuần hoàn bằng bộ cơ sở sóng phẳng mở rộng theo định lý Bloch với các điều kiện biên tuần hoàn (PBCs) trên các ô siêu mạng (supercell). Tương tác giữa lõi và electron hóa trị được xử lý thông qua phương pháp sóng phẳng tăng cường xạ ảnh (PAW) của Blöchl, kết hợp thuật toán tối ưu hóa hình học và giải lặp tự nhất quán (SCF) được triển khai trong các gói phần mềm lượng tử chuyên dụng (như VASP).

Tư liệu và dữ liệu nghiên cứu bao gồm các thông số thực nghiệm đối chiếu: hằng số mạng tinh thể của graphene ($a = 2{,}463\text{ \AA}$, $r_{\text{C}-\text{C}} = 1{,}422\text{ \AA}$), hằng số mạng và độ rộng vùng cấm thực nghiệm của các pha $\text{TiO}_2$, phổ nhiễu xạ tia X (XRD), cùng các dữ liệu phổ cộng hưởng từ hạt nhân $^{13}\text{C}$ SSNMR của graphene oxide.

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

Chương 1: Phương pháp phiếm hàm mật độ (DFT)

Chương 1 trình bày cơ sở toán học và vật lý của các phương pháp tính toán từ mô hình Thomas-Fermi (1927) biểu diễn động năng qua mật độ electron $\rho(r)$, chỉ ra các hạn chế trong việc bỏ qua năng lượng trao đổi - tương quan. Tiếp đó, chương phân tích chi tiết hệ phương trình Kohn-Sham và sơ đồ giải lặp tự nhất quán.

Chương tổng kết các mức gần đúng phiếm hàm năng lượng trao đổi - tương quan:

  • Xấp xỉ mật độ địa phương (LDA, LSDA) mô tả tốt khí electron đồng nhất nhưng đánh giá quá cao năng lượng liên kết.
  • Xấp xỉ gradient tổng quát (GGA) như PW91, PBE, BLYP, BP86 bổ sung gradient mật độ electron, mô tả tốt hơn các liên kết hóa học và liên kết hydrogen.
  • Phiếm hàm lai hóa (B3LYP, PBE0, HSE06) kết hợp tỷ lệ năng lượng trao đổi Hartree-Fock chính xác.
  • Phiếm hàm van der Waals (vdW-DF của Dion và optPBE-vdW của Klimeš) giải quyết tương tác phân tán London không định xứ.
  • Phương pháp DFT+U của Anisimov và Dudarev nhằm xử lý sai số tự tương tác SIE đối với các electron $d, f$ tương quan mạnh.
  • Lý thuyết sóng phẳng, điều kiện biên tuần hoàn, phép biến đổi mạng đảo vùng Brillouin, cùng sự phát triển của các phương pháp thế giả từ HSC bảo toàn định mức, thế siêu mềm USPP của Vanderbilt đến phương pháp PAW của Blöchl, Kresse và Joubert.

Chương 2: Tổng quan hệ chất và phương pháp tính toán

Chương 2 hệ thống hóa các khái niệm tinh thể học, mạng Bravais, chỉ số Miller $(hkl)$, khoảng cách giữa các mặt phẳng song song $d_{hkl}$, định luật Bragg cho nhiễu xạ tia X ($2d\sin\theta = n\lambda$), mật độ trạng thái (DOS, PDOS), cấu trúc dải năng lượng và mặt Fermi.

Chương phân tích bản chất liên kết hydrogen cổ điển (chuyển dời đỏ: làm dãn dài liên kết $\text{A}-\text{H}$, giảm tần số dao động hóa trị $\nu_{\text{A}-\text{H}}$, tăng cường độ IR, chuyển dịch điện tích $0{,}01 \div 0{,}03\text{ e}$) và liên kết hydrogen chuyển dời xanh (rút ngắn liên kết $\text{A}-\text{H}$, tăng tần số dao động, năng lượng từ vài đến vài chục $\text{kJ/mol}$).

Chương tổng hợp dữ liệu cấu trúc và tính chất của 3 dạng thù hình phổ biến của $\text{TiO}_2$:

Thuộc tính Rutile Anatase Brookite
Cấu trúc tinh thể Hệ tứ giác Hệ tứ giác Hệ trực giao
Hằng số mạng ($\text{\AA}$) $a = 4{,}59;\ c = 2{,}95$ $a = 3{,}78;\ c = 9{,}52$ $a = 9{,}18;\ b = 5{,}45;\ c = 5{,}15$
Nhóm không gian $P4_2/mnm$ $I4_1/amd$ $Pbca$
Độ rộng vùng cấm $E_g$ ở 10 K (eV) $3{,}03$ $3{,}46$ —

Đối với graphene, chương khảo sát hai mô hình ô đơn vị:

  • Ô đơn vị lục giác: Gồm 2 nguyên tử C, vector cơ sở có độ lớn $a = 2{,}463\text{ \AA}$, các điểm đối xứng cao trong vùng Brillouin là $\Gamma, \text{K}, \text{M}$.
  • Ô đơn vị hình chữ nhật: Gồm 4 nguyên tử C, các hằng số mạng tương ứng là $a_3 = 2{,}463\text{ \AA}$ và $a_4 = 4{,}266\text{ \AA}$, các điểm đối xứng cao là $\Gamma, \text{X}, \text{Y}, \text{W}$. Các electron quanh điểm Dirac có vận tốc Fermi $v_F \approx 1 \times 10^6\text{ m/s}$.

Chương 3: Kết quả và thảo luận

Chương 3 trình bày các kết quả tính toán chi tiết:

  1. So sánh hai loại ô đơn vị của graphene: Tính toán cấu trúc dải năng lượng, DOS và PDOS cho thấy ô đơn vị hình chữ nhật 4 nguyên tử phản ánh hoàn toàn chính xác các tính chất electron của ô đơn vị lục giác truyền thống, bảo toàn điểm Dirac và độ rộng vùng cấm bằng 0.
  2. Dẫn xuất hydroxyl graphene ($\text{G}_n\text{OH}$): Khảo sát các cấu hình từ $\text{G}_1\text{OH}$ đến $\text{G}_6\text{OH}$. Tính toán năng lượng liên kết $E_b$, độ dài liên kết và góc liên kết cho thấy các nhóm -OH 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 đáy graphene, tạo thành các vòng hexa-hydroxyl khép kín. Cấu trúc này được làm bền nhờ sự hình thành đồng thời của liên kết hydrogen chuyển dời đỏ $\text{O}-\text{H}\cdots\text{O}$ và liên kết hydrogen chuyển dời xanh $\text{O}-\text{H}\cdots\pi$.
  3. Dẫn xuất một nhóm chức ($\text{G}_1\text{F}$): So sánh năng lượng liên kết chỉ ra thứ tự độ bền giảm dần: $\text{G}{\text{epo}} > \text{G}{\text{CH}3} > \text{G}{\text{NH}2} > \text{G}{\text{COOH}}$. Dẫn xuất $\text{G}_{\text{COOH}}$ có độ bền liên kết rất thấp trên mặt phẳng graphene. Mật độ trạng thái quanh mức Fermi của các dẫn xuất $\text{G}_1\text{F}$ được đóng góp chủ yếu bởi orbital $2p_z$ của nguyên tử thuộc nhóm chức liên kết trực tiếp với nguyên tử carbon của mạng graphene.
  4. Khảo sát tinh thể rutile $\text{TiO}_2$ bằng DFT+U: Nghiên cứu ảnh hưởng của $U_d$ và $U_p$ độc lập và đồng thời. Kết quả xác định giá trị $U_p = 10\text{ eV}$ là tối ưu; việc kết hợp $U_d = 7\text{ eV}$ và $U_p = 10\text{ eV}$ vừa điều chỉnh độ rộng vùng cấm $E_g$ tiệm cận giá trị thực nghiệm ($3{,}03\text{ eV}$), vừa kiểm soát sai lệch hằng số mạng $a, c, u$ ở mức tối thiểu so với chỉ dùng $U_d$.
  5. Bề mặt rutile $\text{TiO}_2$ (110): Xây dựng và đánh giá 3 mô hình phiến (slab): FR (thả lỏng toàn bộ - full relax), FIL (cố định các lớp bên trong - fix inner layer), và F2B (cố định hai lớp đáy - fix two bottom). Kết quả cho thấy phương pháp DFT+U khắc phục được hiện tượng dao động chẵn - lẻ của các tính chất bề mặt theo số lớp. Mô hình F2B giúp tính chất bề mặt nhanh hội tụ theo độ dày lớp nhưng gây ra sai số độ rộng vùng cấm $E_g$ lớn nhất.
  6. Mô hình composite hai chiều: Xây dựng ô mạng giao diện giữa ô đơn vị hình chữ nhật của graphene/dẫn xuất graphene với mặt rutile $\text{TiO}_2$ (110) thông qua việc khảo sát các tỷ lệ ghép vector mạng $mA$ với $pA'$, $nB$ với $qB'$, $mA$ với $qB'$, và $nB$ với $pA'$, thiết lập mô hình tính toán giao diện tuần hoàn hai chiều với độ bất đối xứng mạng ở mức thấp nhất.

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

Luận án đã đạt được các đóng góp mới về mặt lý luận và khoa học tính toán:

  1. Về cơ chế phân bố nhóm chức trên graphene: Chỉ ra quy luật tập hợp của các nhóm hydroxyl trong dẫn xuất $\text{G}_n\text{OH}$ tại các vị trí para trên cùng một phía của mặt phẳng graphene để tạo thành các vòng hexa-hydroxyl hoàn hảo, làm sáng tỏ vai trò đồng thời của liên kết hydrogen chuyển dời đỏ $\text{O}-\text{H}\cdots\text{O}$ và liên kết hydrogen chuyển dời xanh $\text{O}-\text{H}\cdots\pi$.
  2. Về độ bền tương đối của các dẫn xuất graphene $\text{G}_1\text{F}$: Xác lập trật tự độ bền liên kết $\text{G}{\text{epo}} > \text{G}{\text{CH}3} > \text{G}{\text{NH}2} > \text{G}{\text{COOH}}$, cung cấp bằng chứng lý thuyết giải thích vì sao nhóm -COOH kém bền trên mặt đáy graphene và thực tế chủ yếu tồn tại ở mép tấm hoặc vị trí khuyết tật. Đồng thời chỉ ra các trạng thái electron lân cận mức Fermi bị chi phối bởi orbital $2p_z$ của nguyên tử liên kết trực tiếp.
  3. Về phương pháp luận tính toán rutile $\text{TiO}_2$ và bề mặt (110): Xác định bộ tham số tối ưu $U_p = 10\text{ eV}$ kết hợp với $U_d$ trong phương pháp DFT+U, hạn chế được sai số hằng số mạng và cải thiện độ chính xác của vùng cấm $E_g$. Đánh giá định lượng sai số của 3 mô hình slab (FR, FIL, F2B), chỉ ra phương pháp DFT+U làm triệt tiêu sự dao động chẵn - lẻ theo số lớp nguyên tử.
  4. Về mô hình hóa cấu trúc composite: Đề xuất và xây dựng thành công mô hình ô mạng tuần hoàn hai chiều ghép nối giữa ô đơn vị hình chữ nhật 4 nguyên tử của graphene (và dẫn xuất hydroxyl graphene) với mặt rutile $\text{TiO}_2$ (110), giải quyết bài toán sai lệch góc mạng tinh thể.

Về mặt kiến nghị, tác giả đề xuất tiếp tục áp dụng mô hình ô đơn vị hình chữ nhật và bộ tham số DFT+U tối ưu để khảo sát sâu hơn về động học chuyển dịch điện tích, phổ hấp thụ quang học và cơ chế phản ứng quang xúc tác chi tiết trên giao diện composite $\text{G}_n\text{OH}/\text{rutile TiO}_2\text{ (110)}$.

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

Văn bản luận án nêu rõ các giới hạn nghiên cứu:

  • Nghiên cứu tập trung vào mô hình hóa lý thuyết ở trạng thái cơ bản ($0\text{ K}$) dưới điều kiện chân không, chưa khảo sát ảnh hưởng động học nhiệt độ hữu hạn và môi trường dung môi xung quanh.
  • Khảo sát nhóm chức hydroxyl mới thực hiện đến số lượng $n = 6$ trên cùng một phía của mặt phẳng graphene; chưa mở rộng cho các cấu hình phân bố đồng thời ở cả hai phía mặt phẳng hoặc tổ hợp đa nhóm chức phức tạp (như đồng thời epoxy và hydroxyl với mật độ cao).
  • Hướng nghiên cứu tiếp theo được tác giả xác định là hoàn thiện tính toán tương tác phân tán và chuyển dịch điện tích trên toàn bộ mô hình composite hai chiều đã xây dựng, mở rộng khảo sát các trạng thái kích thích quang hóa của hệ composite.

Giá trị tham khảo

Luận án là tài liệu tham khảo trong các lĩnh vực:

  • Nghiên cứu sinh và học viên cao học ngành Hóa lý thuyết, Hóa lý, Vật lý chất rắn, Khoa học vật liệu: Tham khảo quy trình xây dựng ô đơn vị hình chữ nhật cho vật liệu hai chiều dạng tổ ong, phương pháp xác định tham số hiệu chỉnh $U_d, U_p$ trong DFT+U, và phương pháp tính toán năng lượng bề mặt trên các mô hình slab khác nhau.
  • Giảng viên và các nhà nghiên cứu quang xúc tác: Tham khảo các dữ liệu định lượng về độ bền liên kết của các nhóm chức trên graphene, cơ chế hình thành liên kết hydrogen $\text{O}-\text{H}\cdots\pi$, và cấu trúc vùng năng lượng giao diện vật liệu composite carbon/oxide kim loại.

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

1. Tại sao luận án sử dụng ô đơn vị hình chữ nhật 4 nguyên tử thay vì ô đơn vị lục giác 2 nguyên tử truyền thống của graphene khi mô hình hóa composite với rutile $\text{TiO}_2$ (110)?
Bề mặt rutile $\text{TiO}_2$ (110) có cấu trúc mạng ô đơn vị hình chữ nhật (góc mạng $90^\circ$). Nếu sử dụng ô đơn vị lục giác truyền thống của graphene (góc mạng $120^\circ$), sự sai lệch góc không thể bù trừ bằng cách mở rộng supercell, gây biến dạng mạng rất lớn trong mô hình tính toán. Ô đơn vị hình chữ nhật 4 nguyên tử của graphene có tính đối xứng tương thích ($90^\circ$), giúp ghép nối chính xác với mặt rutile (110) mà vẫn bảo toàn đầy đủ các đặc tính electron của graphene.

2. Quy luật sắp xếp không gian của các nhóm hydroxyl (-OH) trên mặt đáy graphene được luận án phát hiện như thế nào?
Các nhóm -OH có xu hướng định vị tại các vị trí para trên cùng một phía của mặt phẳng graphene để tạo thành các vòng hexa-hydroxyl hoàn chỉnh. Cấu hình này đạt năng lượng bền vững nhờ sự phối hợp đồng thời giữa liên kết hydrogen chuyển dời đỏ thông thường ($\text{O}-\text{H}\cdots\text{O}$) giữa các nhóm -OH và liên kết hydrogen chuyển dời xanh ($\text{O}-\text{H}\cdots\pi$) giữa nhóm -OH với hệ electron $\pi$ của khung carbon.

3. Thứ tự độ bền liên kết của các dẫn xuất graphene một nhóm chức ($\text{G}_1\text{F}$) được xác định ra sao?
Năng lượng liên kết $E_b$ cho thấy thứ tự độ bền giảm dần: $\text{G}{\text{epo}} > \text{G}{\text{CH}3} > \text{G}{\text{NH}2} > \text{G}{\text{COOH}}$. Dẫn xuất gắn nhóm carboxyl ($\text{G}_{\text{COOH}}$) trên mặt phẳng graphene có độ bền rất thấp, phù hợp với các quan sát thực nghiệm rằng nhóm -COOH chủ yếu chỉ liên kết tại mép hoặc các vị trí khuyết tật của tấm graphene.

4. Vai trò của việc hiệu chỉnh đồng thời tham số Hubbard $U_d$ và $U_p$ đối với rutile $\text{TiO}_2$ là gì?
Nếu chỉ tăng tham số $U_d$ trên orbital Ti-$3d$, dù mở rộng được độ rộng vùng cấm $E_g$ về gần thực nghiệm nhưng lại làm sai lệch lớn các hằng số mạng tinh thể ($a, c, u$). Việc tối ưu hóa đồng thời với $U_p = 10\text{ eV}$ trên orbital O-$2p$ kết hợp với $U_d = 7\text{ eV}$ giúp vừa đưa $E_g$ tiệm cận giá trị thực nghiệm ($3{,}03\text{ eV}$), vừa giảm thiểu sai số của các hằng số mạng mạng tinh thể.

5. Trong ba mô hình slab khảo sát bề mặt rutile (110) (FR, FIL, F2B), mô hình nào có sai số vùng cấm lớn nhất?
Mô hình F2B (cố định hai lớp nguyên tử ở đáy) mặc dù giúp các tính chất bề mặt nhanh chóng đạt trạng thái hội tụ theo số lớp, nhưng lại là mô hình gây ra sai số về độ rộng vùng cấm ($E_g$) lớn nhất trong ba mô hình được khảo sát.

Kết luận

Luận án của nghiên cứu sinh Trần Thị Thoa đã giải quyết bài toán mô hình hóa lý thuyết cấu trúc và tính chất electron của các dẫn xuất graphene, tinh thể khối và bề mặt rutile $\text{TiO}_2$ (110), cùng cấu trúc giao diện composite hai chiều tương ứng bằng phương pháp phiếm hàm mật độ kết hợp các phép hiệu chỉnh optPBE-vdW và DFT+U. Công trình đã làm sáng tỏ cơ chế liên kết hydrogen ổn định cấu trúc hexa-hydroxyl trên mặt graphene, xác lập trật tự độ bền các dẫn xuất chức hóa $\text{G}_1\text{F}$, xác định bộ tham số tối ưu $U_p = 10\text{ eV}$ cho rutile $\text{TiO}_2$, và thiết lập mô hình ghép mạng hình chữ nhật đồng nhất cho hệ composite graphene/rutile $\text{TiO}_2$ (110). Các kết quả tính toán cung cấp luận cứ định lượng tin cậy cho việc giải thích cấu trúc thực nghiệm và định hướng thiết kế vật liệu quang xúc tác mới.