Tổng quan về luận án

Nghiên cứu về vật liệu composite dị thể hai chiều giữa graphene và titanium dioxide ($\text{TiO}_2$) đại diện cho một trong những hướng tiếp cận tiên phong của khoa học vật liệu và hóa học tính toán hiện đại. Hệ composite này giữ vai trò then chốt trong các phản ứng quang xúc tác phân hủy hợp chất hữu cơ độc hại trong nước thải công nghiệp, phân tách nước sản xuất hydro sạch nhằm chuyển hóa quang năng thành hóa năng, nâng cao hiệu suất pin mặt trời chất màu nhạy quang (DSSC), cũng như ứng dụng làm vật liệu anode cho pin lithium-ion và cảm biến khí. Tuy nhiên, một khoảng trống nghiên cứu sâu sắc (research gap) tồn tại dai dẳng trong y văn quốc tế: phần lớn các nghiên cứu lý thuyết trước đây chỉ tập trung vào bề mặt anatase $\text{TiO}_2$ (101) kết hợp với ô đơn vị lục giác của graphene nhờ sự tương đồng tương đối về đối xứng mạng [Long et al., 2012; Du et al., 2011]. Ngược lại, pha rutile $\text{TiO}_2$ — dạng thù hình bền vững nhất về mặt nhiệt động học với thế năng biến thiên enthalpy chuyển pha anatase–rutile là $-0{,}78 \pm 0{,}2\text{ kcal/mol}$ [Mitsuhashi et al., 1974] — lại bị bỏ ngỏ do sự bất tương thích hình học giữa ô đơn vị lục giác ($\gamma = 120^\circ$) của graphene và ô đơn vị hình chữ nhật ($\gamma = 90^\circ$) của bề mặt rutile (110). Thêm vào đó, trong thực nghiệm tổng hợp nano composite, sản phẩm thu được chủ yếu là graphene oxide (GO) hoặc reduced graphene oxide (RGO) mang các nhóm chức oxy hóa phức tạp, nhưng cấu trúc vi mô và cơ chế sắp xếp của các nhóm hydroxyl (-OH) đa chức trên mặt đáy graphene vẫn chưa từng được sáng tỏ một cách hệ thống.

Luận án tiến sĩ Hóa lý thuyết và Hóa lý của tác giả Trần Thị Thoa, thực hiện tại Trường Đại học Sư phạm Hà Nội dưới sự hướng dẫn của PGS. Nguyễn Thị Minh Huệ và PGS. Hoàng Văn Hùng (Mã số: 9440119, năm 2022), đã giải quyết triệt để các thách thức cốt lõi trên thông qua ba câu hỏi nghiên cứu và giả thuyết khoa học:

  1. Câu hỏi nghiên cứu 1 (RQ1): Quy luật sắp xếp hình học, độ bền liên kết và cấu trúc electron của các dẫn xuất hydroxyl graphene ($G_n\text{OH}$, với $n = 1 \div 6$) cùng các nhóm chức đơn ($G_1F$: -$\text{NH}_2$, -$\text{CH}_3$, -$\text{OCH}_3$, -$\text{CHO}$, -$\text{COOH}$, epoxy) trên mặt phẳng graphene diễn ra như thế nào?
    Giả thuyết 1 (H1): Các nhóm -OH tự tổ chức ưu tiên 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 hoàn hảo, được làm bền bởi mạng lưới liên kết hydrogen chuyển dời đỏ $\text{O-H}\cdots\text{O}$ và chuyển dời xanh $\text{O-H}\cdots\pi$.
  2. Câu hỏi nghiên cứu 2 (RQ2): Làm thế nào để triệt tiêu sai số tự tương tác (Self-Interaction Error - SIE) cố hữu của lý thuyết phiếm hàm mật độ (DFT) nhằm mô tả chính xác thông số mạng và độ rộng vùng cấm ($E_g$) của rutile $\text{TiO}_2$ khối và bề mặt rutile (110)?
    Giả thuyết 2 (H2): Phương pháp $\text{DFT}+U$ với sự hiệu chỉnh đồng thời tham số Hubbard $U_d$ trên orbital $\text{Ti } 3d$ và $U_p$ trên orbital $\text{O } 2p$ với giá trị tối ưu $U_p = 10\text{ eV}$ sẽ loại bỏ hiện tượng dao động chẵn–lẻ theo số lớp slab và tái lập chính xác thực nghiệm $E_g = 3{,}03\text{ eV}$.
  3. Câu hỏi nghiên cứu 3 (RQ3): Mô hình tinh thể học nào cho phép ghép nối giao diện dị thể hai chiều giữa dẫn xuất graphene và bề mặt rutile $\text{TiO}_2$ (110) mà không gây biến dạng mạng tinh thể quá mức?
    Giả thuyết 3 (H3): Việc chuyển đổi mạng graphene sang ô đơn vị hình chữ nhật 4 nguyên tử ($a = 2{,}463\text{ \AA}$, $b = 4{,}266\text{ \AA}$) cho phép khớp nối trực giao với ô đơn vị chữ nhật của rutile (110) ($a = 4{,}595\text{ \AA}$, $b = 2{,}959\text{ \AA}$), tối thiểu hóa sai lệch mạng (lattice mismatch).

Khung lý thuyết của công trình tích hợp lý thuyết phiếm hàm mật độ mở rộng hiệu chỉnh lực phân tán van der Waals ($\text{optPBE-vdW}$), phương pháp $\text{DFT}+U$ hiệu chỉnh tương quan mạnh, và lý thuyết liên kết hydrogen phi cổ điển. Phạm vi nghiên cứu bao quát từ mức độ đơn lớp, tinh thể khối, các mô hình bề mặt slab (FR, FIL, F2B từ 3 đến 9 lớp) cho đến siêu mạng composite 2D, tạo ra bước đột phá định lượng hóa cơ chế chuyển dịch điện tích quang sinh tại giao diện tiếp xúc.

Literature Review và Positioning

Bức tranh tổng quan y văn về vật liệu graphene và oxide kim loại chuyển tiếp được định hình bởi ba dòng nghiên cứu chính:

  1. Dòng nghiên cứu cấu trúc Graphene Oxide: Các mô hình thực nghiệm cổ điển của Lerf-Klinowski [1998] và Dékány [2006] cùng phân tích cộng hưởng từ hạt nhân trạng thái rắn ($^{13}\text{C SSNMR}$) của Cai et al. [2008] xác nhận sự tồn tại của ba pic chính: nhóm epoxy ở $60\text{ ppm}$, nhóm $\text{C-OH}$ ở $70\text{ ppm}$, và carbon $\text{sp}^2$ ở $130\text{ ppm}$, chứng minh các nhóm chức phân bố sát nhau trên mặt đáy. Tuy nhiên, các tính toán lý thuyết ban đầu của Lahaye et al. [2009] và Yan et al. [2009] chỉ dừng lại ở cấu hình đơn giản gồm 1 hoặc 2 nhóm -OH, khẳng định hai nhóm -OH không thể đính vào hai carbon liền kề mà có xu hướng định vị para, nhưng hoàn toàn bế tắc khi số lượng nhóm -OH gia tăng lên $n \ge 3$.
  2. Dòng nghiên cứu sai số SIE trong $\text{TiO}_2$: Các phiếm hàm LDA và GGA truyền thống luôn đánh giá thấp nghiêm trọng độ rộng vùng cấm của rutile $\text{TiO}_2$ do sự không định xứ quá mức của các electron $d$ và $p$ [Perdew et al., 1996; Anisimov et al., 1991]. Các công trình của Morgan & Watson [2007], Nolan et al. [2008], và Janotti et al. [2010] đã áp dụng hiệu chỉnh Hubbard $U_d$ lên orbital $\text{Ti } 3d$. Mặc dù $U_d \approx 10\text{ eV}$ giúp mở rộng vùng cấm tiệm cận giá trị thực nghiệm $3{,}03\text{ eV}$ [Tang et al., 1995], nó lại gây ra sự biến dạng sai lệch nghiêm trọng về hằng số mạng tinh thể $a$ và $c$.
  3. Dòng nghiên cứu mô hình hóa giao diện dị thể Composite: Yamashita et al. [2011] và Martsinovich et al. [2017] lần đầu tiên đề xuất việc sử dụng ô đơn vị hình chữ nhật của graphene cho composite với anatase (001) và rutile (110). Tuy nhiên, các công trình quốc tế quy mô của Long et al. [2012] và Du et al. [2011] khi mô phỏng graphene/rutile (110) vẫn áp dụng ô đơn vị lục giác nguyên bản, dẫn đến sai lệch góc quay $30^\circ$ và gây ra ứng suất nhân tạo (artificial strain) làm sai lệch cấu trúc dải năng lượng.

Tranh luận học thuật sâu sắc nhất trong lĩnh vực nằm ở hai điểm đối kháng:

  • Tranh luận 1 (Bản chất liên kết hydrogen chuyển dời xanh): Giữa trường phái ủng hộ cơ chế tái tổ hợp lai hóa orbital $s$ làm rút ngắn liên kết $\text{A-H}$ [Hobza et al., 1999] và trường phái lý giải bằng tương tác siêu liên hợp nội/ngoại phân tử chuyển dịch điện tích [Alabugin et al., 2003; Weinhold, 2001] kết hợp hiệu ứng phân cực trường điện [Hermansson, 2002].
  • Tranh luận 2 (Phương pháp xử lý tương quan điện tử trong rutile): Giữa việc sử dụng phiếm hàm lai hóa đắt đỏ ($B_3LYP$, $HSE06$, $PBE0$) và phương pháp $\text{DFT}+U$ tiết kiệm chi phí nhưng thiếu nhất quán về việc lựa chọn chỉ số $U$.

Luận án của NCS. Trần Thị Thoa đã định vị xuất sắc ở giao điểm này, tiến hành so sánh đối chuẩn trực tiếp với hai nghiên cứu quốc tế của Long et al. [2012] và Martsinovich et al. [2017], mở rộng biên giới hiểu biết về sự tự tổ chức của nhóm chức oxy hóa và cơ chế truyền điện tích dị thể.

Đóng góp lý thuyết và khung phân tích

Đóng góp cho lý thuyết

Luận án mang lại bốn đóng góp đột phá làm thay đổi nhận thức lý thuyết trong hóa lý chất rắn và vật liệu nano:

  1. Mở rộng lý thuyết liên kết Hydrogen vào thế giới vật liệu 2D: Mở rộng khuôn khổ kinh điển của Pauling [1931] và Coulson [1957] về liên kết hydrogen chuyển dời đỏ cổ điển $\text{A-H}\cdots\text{B}$ sang liên kết phi cổ điển chuyển dời xanh $\text{O-H}\cdots\pi$. Công trình chứng minh trực tiếp rằng sự hình thành các vòng hexa-hydroxyl bền vững trên bề mặt graphene là kết quả cộng hưởng của cả hai loại liên kết: liên kết hydrogen chuyển dời đỏ $\text{O-H}\cdots\text{O}$ nội vòng và liên kết chuyển dời xanh $\text{O-H}\cdots\pi$ giữa proton nhóm -OH với đám mây electron $\pi$ của mạng carbon.
  2. Khái quát hóa quy luật tự tổ chức para của dẫn xuất $G_n\text{OH}$: Bác bỏ quan điểm cho rằng các nhóm hydroxyl phân bố ngẫu nhiên. Bằng chứng tính toán lượng tử khẳng định mệnh đề: "Trong dẫn xuất $G_n\text{OH}$, 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" [Trần Thị Thoa, 2022, tr. 4].
  3. Hoàn thiện phương pháp $\text{DFT}+U$ đồng thời ($U_d - U_p$): Vượt qua hạn chế của mô hình Dudarev đơn biến số, luận án chứng minh rằng việc áp dụng đồng thời tham số Hubbard cho cả vùng dẫn (orbital $\text{Ti } 3d$, $U_d = 7\text{ eV}$) và vùng hóa trị (orbital $\text{O } 2p$, $U_p = 10\text{ eV}$) là điều kiện tiên quyết để triệt tiêu sai số tự tương tác SIE, đưa độ rộng vùng cấm rutile về chính xác giá trị thực nghiệm mà vẫn bảo toàn cấu trúc mạng tinh thể.
  4. Xác lập trật tự độ bền nhiệt động của các dẫn xuất đơn chức ($G_1F$): Thiết lập dãy năng lượng liên kết ($E_b$): $G_{\text{epo}} > G_{\text{CH}3} > G{\text{NH}2} > G{\text{OCH}3} > G{\text{CHO}} \gg G_{\text{COOH}}$, cung cấp cơ sở nhiệt động học sáng tỏ vì sao nhóm -COOH kém bền nhất và chỉ xuất hiện ở các cạnh khuyết tật trong thực nghiệm tổng hợp GO.

Khung phân tích độc đáo

Khung phân tích của luận án được xây dựng dựa trên sự tích hợp chặt chẽ của ba trụ cột lý thuyết: (1) Lý thuyết phiếm hàm mật độ hiệu chỉnh phân tán $\text{optPBE-vdW}$, (2) Mô hình hiệu chỉnh tương quan Hubbard đa orbital $\text{DFT}+U_d+U_p$, và (3) Phương pháp hàm sóng tăng cường sóng phẳng (PAW).

Đóng góp khái niệm trung tâm bao gồm:

  • Khái niệm "Vòng Hexa-hydroxyl Para": Cấu trúc liên kết gồm 6 nhóm -OH đính xen kẽ tại các vị trí carbon đối diện (para) trong vòng sáu cạnh, tối ưu hóa năng lượng liên kết $E_b$ thông qua mạng lưới cho–nhận proton khép kín.
  • Điều kiện biên phân tích: Tất cả các mô hình được tính toán với điều kiện biên tuần hoàn 3 chiều (PBCs), lớp chân không dày tối thiểu $15\text{ \AA}$ dọc trục $z$ nhằm loại bỏ tương tác giả giữa các phiến tuần hoàn, và nhiệt độ quy chuẩn tại $0\text{ K}$.

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 theo hệ hình triết học thực chứng tính toán (computational positivism) kết hợp với hiện thực phản biện (critical realism), nơi các đại lượng cơ lượng tử (mật độ electron, cấu trúc dải năng lượng, ten-sơ ứng suất) được tính toán từ các nguyên lý ban đầu (ab initio) để phản ánh bản chất khách quan của cấu trúc vật chất. Thiết kế nghiên cứu đa cấp độ (multi-level design) được phân tầng cụ thể:

  • Cấp độ 1 (Orbital & Phân tử): Tương tác xen phủ orbital $\text{C } 2p_z$, $\text{O } 2p$, $\text{Ti } 3d$, dịch chuyển mật độ electron (EDD).
  • Cấp độ 2 (Bề mặt & Màng đơn lớp): Nghiên cứu các dẫn xuất $G_n\text{OH}$, $G_1F$, và bề mặt slab rutile (110) theo các mô hình thư giãn khác nhau.
  • Cấp độ 3 (Giao diện dị thể 2D Composite): Khớp nối siêu ô mạng (supercell) giữa graphene và rutile $\text{TiO}_2$ (110).

Kích thước mẫu tính toán bao gồm 6 họ dẫn xuất hydroxyl ($n = 1 \div 6$ với hàng chục đồng phân hình học), 6 loại nhóm chức đơn $G_1F$, 3 mô hình slab bề mặt rutile (FR - Full Relax, FIL - Fix Inner Layer, F2B - Fix Two Bottom) với độ dày từ 3 đến 9 lớp nguyên tử, và hàng trăm cấu hình ghép nối dị thể $mA \times pA'$ và $nB \times qB'$.

Quy trình nghiên cứu rigorous

Quy trình tính toán lượng tử tuân thủ nghiêm ngặt các tiêu chuẩn kiểm soát chất lượng quốc tế thông qua phần mềm VASP:

  • Bộ hàm cơ sở & Thế giả: Sử dụng phương pháp sóng tăng cường chiếu (Projector Augmented Wave - PAW) của Blöchl [1994], năng lượng cắt sóng phẳng (plane-wave cutoff energy) được thiết lập ở mức cao $E_{\text{cut}} = 500\text{ eV}$.
  • Lấy mẫu vùng Brillouin: Lưới điểm $k$-point Monkhörst-Pack được kiểm tra hội tụ nghiêm ngặt: $15 \times 15 \times 1$ cho ô đơn vị graphene nguyên bản, $9 \times 9 \times 1$ cho supercell $6 \times 3$ của graphene, và $6 \times 6 \times 6$ cho tinh thể rutile $\text{TiO}_2$ khối. Các điểm đối xứng cao được quét chính xác: $\Gamma - \text{K} - \text{M}$ cho ô lục giác và $\Gamma - \text{X} - \text{Y} - \text{W}$ cho ô chữ nhật.
  • Tiêu chuẩn hội tụ: Lực tác dụng lên mỗi nguyên tử sau khi tối ưu hóa hình học nhỏ hơn $0{,}01\text{ eV/\AA}$, và tổng năng lượng điện tử hội tụ ở mức sai số dưới $10^{-6}\text{ eV}$.
  • Độ tin cậy và kiểm chuẩn: Sự biến thiên năng lượng liên kết $E_b$ và khoảng cách liên kết được đối chiếu chéo giữa các phiếm hàm LDA, GGA-PBE, và $\text{optPBE-vdW}$.

Data và phân tích

Phân tích dữ liệu định lượng sâu sắc dựa trên các thông số vật lý then chốt:

  • Hằng số mạng tinh thể: Ô đơn vị hình chữ nhật 4 nguyên tử của graphene đạt $a_3 = 2{,}463\text{ \AA}$, $a_4 = 4{,}266\text{ \AA}$, trùng khớp tuyệt đối với quan hệ hình học $a_3 = \sqrt{3} r_{\text{CC}}$ và $a_4 = 3 r_{\text{CC}}$ ($r_{\text{CC}} = 1{,}422\text{ \AA}$).
  • Sai số tương quan Hubbard: Khi chỉ áp dụng $U_d = 7\text{ eV}$, $E_g$ của rutile tăng lên nhưng hằng số mạng $a$ lệch $+1{,}2%$. Khi đồng hiệu chỉnh $U_d = 7\text{ eV}$ và $U_p = 10\text{ eV}$, độ rộng vùng cấm đạt $3{,}01\text{ eV}$ (sai số $< 0{,}7%$ so với thực nghiệm $3{,}03\text{ eV}$ ở $10\text{ K}$), đồng thời hằng số mạng $a = 4{,}594\text{ \AA}$ và $c = 2{,}961\text{ \AA}$ đạt độ chính xác gần như tuyệt đối.

Phát hiện đột phá và implications

Những phát hiện then chốt

Công trình công bố 5 phát hiện cốt lõi mang tính bước ngoặt:

  1. Sự hình thành vòng lặp Hexa-Hydroxyl bền vững: Trong các dẫn xuất $G_n\text{OH}$, khi số lượng nhóm -OH tăng từ $1$ đến $6$, năng lượng liên kết trung bình trên mỗi nhóm -OH ($E_b$) giảm dần (tức độ bền tăng lên), đạt cực đại bền vững tại $G_6\text{OH}$ với cấu hình 6 nhóm -OH định vị para khép kín trên cùng một phía mặt phẳng carbon.
  2. Cơ chế ổn định lưỡng tính của liên kết Hydrogen: Phát hiện sự cùng tồn tại của hai hiệu ứng liên kết: liên kết hydrogen chuyển dời đỏ $\text{O-H}\cdots\text{O}$ với khoảng cách $\text{H}\cdots\text{O}$ ngắn ($1{,}85 \div 2{,}05\text{ \AA}$) làm kéo dài liên kết cộng hóa trị $\text{O-H}$, và liên kết chuyển dời xanh $\text{O-H}\cdots\pi$ tương tác trực tiếp với hệ electron phân tán của vòng benzen lân cận.
  3. Quy luật mật độ trạng thái quanh mức Fermi của $G_1F$: Đối với các dẫn xuất $G_1F$, các trạng thái điện tử cục bộ nằm sát mức Fermi ($E_F$) được đóng góp chủ yếu bởi các orbital $2p_z$ của nguyên tử $X$ thuộc nhóm chức liên kết trực tiếp với carbon graphene, chứng minh khả năng tinh chỉnh vùng cấm graphene thông qua hóa học nhóm chức.
  4. Giải mã nghịch lý mô hình Slab rutile (110): Khảo sát ba mô hình slab cho thấy: "Sự dao động chẵn-lẻ theo số lớp của các tính chất bề mặt rutile (110) có thể được hạn chế bằng phương pháp $\text{DFT}+U$. Mô hình F2B mặc dù nhanh chóng hội tụ các tính chất của bề mặt theo số lớp, nhưng lại gây ra sai số độ rộng vùng cấm ($E_g$) lớn nhất" [Trần Thị Thoa, 2022, tr. 4]. Mô hình FR (Full Relax) với số lớp $\ge 5$ là mô hình mô tả chính xác nhất bản chất vật lý bề mặt.
  5. Thiết lập mô hình 2D Composite không ứng suất: Xây dựng thành công ma trận ghép nối giữa siêu ô mạng graphene hình chữ nhật ($m \times n$) và rutile (110) ($p \times q$), chỉ ra cặp ghép $2A \times 1A'$ và $3B \times 4B'$ triệt tiêu độ sai lệch mạng xuống dưới $1{,}8%$, tạo điều kiện lý tưởng cho sự chuyển dịch điện tích quang sinh từ rutile sang graphene.

Implications đa chiều

  • Về mặt lý thuyết: Cung cấp cơ sở vi mô vững chắc cho các mô hình cấu trúc RGO/GO thực nghiệm, lấp đầy khoảng trống cấu trúc mà các kỹ thuật nhiễu xạ tia X (XRD) và phổ NMR chưa thể giải mã chi tiết.
  • Về mặt phương pháp luận: Chuẩn hóa quy trình tính toán $\text{DFT}+U_d+U_p$ cho các oxide kim loại chuyển tiếp có vùng hóa trị tạo bởi orbital $p$ và vùng dẫn tạo bởi orbital $d$, có thể chuyển giao trực tiếp cho các hệ vật liệu tương tự như $\text{ZnO}$, $\text{Fe}_2\text{O}_3$, $\text{WO}_3$.
  • Về mặt thực tiễn và công nghệ: Định hướng chính xác cho các nhà thực nghiệm trong việc kiểm soát mức độ khử GO; đề xuất giải pháp chức hóa bề mặt bằng nhóm hydroxyl và amine nhằm tối đa hóa hiệu suất phân tách cặp electron-lỗ trống quang sinh, nâng cao hiệu suất quang xúc tác xử lý nước thải dệt nhuộm và sản xuất pin nhiên liệu hydro.
  • Khả năng tổng quát hóa: Khung tính toán ô đơn vị hình chữ nhật có thể mở rộng cho toàn bộ các cấu trúc dị thể giữa vật liệu 2D (như $\text{MoS}_2$, $\text{WS}_2$, phosphorene, MXenes) và các bề mặt oxide có đối xứng trực giao hoặc tứ giác.

Limitations và Future Research

Nhằm duy trì tính trung thực và chuẩn mực học thuật khắt khe, luận án thẳng thắn chỉ rõ 4 giới hạn nội tại:

  1. Giới hạn nhiệt động học tại $0\text{ K}$: Các tính toán $\text{DFT}$ tĩnh được thực hiện tại điểm không tuyệt đối ($0\text{ K}$), chưa tính đến dao động mạng nhiệt động (phonon dispersions) và năng lượng tự do Gibbs ở nhiệt độ phòng ($298\text{ K}$).
  2. Hiệu ứng môi trường dung môi: Mô hình mô phỏng trong chân không, chưa bao gồm các phân tử nước hoặc dung môi phân cực xung quanh — yếu tố có thể cạnh tranh tạo liên kết hydrogen với các nhóm -OH trên bề mặt graphene.
  3. Khuyết tật mạng thực tế: Nghiên cứu tập trung vào bề mặt lý tưởng, chưa khảo sát các khuyết tật điểm thường gặp trong thực nghiệm như khuyết tật vacancy oxy trên mặt rutile (110) hoặc khuyết tật Stone-Wales trên graphene.
  4. Quy mô siêu mạng composite: Do hạn chế về tài nguyên tính toán hiệu năng cao, mô hình composite mới dừng lại ở bước khảo sát hình học ô mạng tiếp xúc tối ưu, chưa thực hiện tính toán cấu trúc dải năng lượng toàn phần cho supercell composite chứa hàng trăm nguyên tử.

Chương trình nghiên cứu 10 năm tới (Future Research Agenda)

  1. Thực hiện mô phỏng động lực học phân tử lượng tử (Ab initio Molecular Dynamics - AIMD) ở $300\text{ K}$ và $500\text{ K}$ để khảo sát độ bền nhiệt của vòng hexa-hydroxyl.
  2. Ứng dụng mô hình dung môi hiện (explicit solvent) kết hợp mô hình điện môi liên tục (implicit solvent - VASPsol) để mô phỏng chính xác môi trường quang xúc tác trong pha lỏng.
  3. Tính toán trạng thái kích thích quang học bằng phương pháp đắt tiền hơn như lý thuyết nhiễu loạn nhiều hạt $GW$ và phương trình Bethe-Salpeter (BSE) để thu được phổ hấp thụ UV-Vis chính xác tuyệt đối.
  4. Mở rộng khảo sát tương tác dị thể với sự xuất hiện của các khuyết tật mạng điểm ($\text{Ti}^{3+}$ interstitials và $\text{O}$ vacancies) nhằm mô phỏng sát thực nghiệm hơn nữa.

Tác động và ảnh hưởng

  • Tác động học thuật: Luận án mở ra một hướng tiếp cận mẫu mực trong hóa học tính toán tại Việt Nam, dự kiến thu hút trích dẫn cao trong các tạp chí ISI/Scopus chuyên ngành vật liệu quang xúc tác và nano carbon.
  • Chuyển đổi công nghiệp: Cung cấp dữ liệu nền tảng cho các tập đoàn hóa chất, năng lượng tái tạo và xử lý môi trường trong việc thiết kế các màng lọc nano composite $\text{TiO}_2/\text{RGO}$ có tuổi thọ cao và hoạt tính quang hóa vượt trội.
  • Ý nghĩa xã hội & môi trường: Đóng góp giải pháp công nghệ vật liệu xanh nhằm giải quyết bài toán ô nhiễm nguồn nước công nghiệp và thúc đẩy lộ trình phát triển kinh tế hydro (Hydrogen Economy), giảm thiểu phát thải khí nhà kính.
  • Tầm vóc quốc tế: Kết quả nghiên cứu thiết lập chuẩn mực mới về mô hình hóa giao diện dị thể 2D trực giao, được công bố và chia sẻ học thuật với cộng đồng nghiên cứu vật lý chất rắn toàn cầu.

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

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 độc đáo nhất là việc phát hiện và chứng minh quy luật tự tổ chức tạo vòng hexa-hydroxyl para trên mặt đáy graphene, mở rộng thuyết liên kết hydrogen của Pauling [1931] và Coulson [1957] sang hệ vật liệu 2D thông qua việc tích hợp đồng thời liên kết hydrogen chuyển dời đỏ $\text{O-H}\cdots\text{O}$ và liên kết chuyển dời xanh $\text{O-H}\cdots\pi$.

2. Cải tiến phương pháp luận nổi bật nhất so với ít nhất hai nghiên cứu quốc tế trước đây?

So với nghiên cứu của Long et al. [2012] và Du et al. [2011] (sử dụng ô đơn vị lục giác gây méo mạng góc $30^\circ$), luận án đã chuẩn hóa việc sử dụng ô đơn vị hình chữ nhật 4 nguyên tử của graphene ($a = 2{,}463\text{ \AA}$, $b = 4{,}266\text{ \AA}$), kết hợp kỹ thuật $\text{DFT}+U_d+U_p$ ($U_p = 10\text{ eV}$) loại bỏ hoàn toàn sai số SIE và triệt tiêu biến dạng mạng nhân tạo.

3. Phát hiện bất ngờ nhất từ dữ liệu tính toán lượng tử là gì?

Phát hiện bất ngờ nhất là mô hình slab F2B (Fix Two Bottom) dù giúp các đại lượng bề mặt hội tụ cực nhanh theo số lớp nhưng lại gây ra sai số vùng cấm $E_g$ lớn nhất; đồng thời giá trị hiệu chỉnh $U_p = 10\text{ eV}$ trên orbital $\text{O } 2p$ đóng vai trò quyết định trong việc mở rộng vùng cấm rutile $\text{TiO}_2$ về giá trị thực nghiệm $3{,}03\text{ eV}$ hơn là việc tăng tham số $U_d$ trên $\text{Ti } 3d$.

4. Luận án có cung cấp quy trình tái lập (Replication Protocol) hoàn chỉnh không?

Hoàn toàn đầy đủ và chi tiết: Luận án cung cấp tường minh tọa độ không gian thực và mạng đảo, các điểm đối xứng cao vùng Brillouin ($\Gamma, \text{K}, \text{M}$ và $\Gamma, \text{X}, \text{Y}, \text{W}$), năng lượng cắt $E_{\text{cut}} = 500\text{ eV}$, lưới $k$-mesh, tham số phiếm hàm $\text{optPBE-vdW}$ ($x = 94{,}5268% \text{ PBE}$, $\mu = 0$), và các ma trận ghép nối siêu ô mạng $mA \times pA'$, $nB \times qB'$.

5. Lộ trình nghiên cứu 10 năm được phác thảo ra sao?

Lộ trình định hướng phát triển từ mô phỏng tĩnh $0\text{ K}$ sang động lực học phân tử lượng tử nhiệt độ phòng (AIMD $300\text{ K}$), tích hợp hiệu ứng solvat hóa phân cực, mô hình hóa các khuyết tật vacancy quang hóa, và chế tạo thử nghiệm màng composite quang xúc tác dị thể nano quy mô pilot.

Kết luận

Công trình nghiên cứu của NCS. Trần Thị Thoa đã xác lập những dấu ấn khoa học xuất sắc với 5 đóng góp cốt lõi:

  1. Giải mã toàn diện cơ chế tự tổ chức không gian của dẫn xuất hydroxyl graphene ($G_n\text{OH}$), chứng minh cấu hình vòng hexa-hydroxyl định vị para là trạng thái bền vững nhiệt động nhất.
  2. Khám phá bản chất tương tác kép của mạng lưới liên kết hydrogen chuyển dời đỏ $\text{O-H}\cdots\text{O}$ và chuyển dời xanh $\text{O-H}\cdots\pi$ trên bề mặt graphene.
  3. Thiết lập chính xác trật tự năng lượng liên kết của 6 dẫn xuất đơn chức ($G_1F$), giải thích hoàn hảo dữ liệu thực nghiệm phổ $^{13}\text{C SSNMR}$.
  4. Tối ưu hóa thành công phương pháp $\text{DFT}+U_d+U_p$ với $U_p = 10\text{ eV}$, khắc phục triệt để sai số tự tương tác SIE và giải thích hiện tượng dao động chẵn–lẻ trên bề mặt rutile (110).
  5. Xây dựng thành công mô hình ghép nối composite 2D giữa dẫn xuất graphene và bề mặt rutile $\text{TiO}_2$ (110) dựa trên ô đơn vị trực giao 4 nguyên tử, triệt tiêu độ sai lệch mạng.

Công trình đã nâng tầm hệ hình nghiên cứu hóa lý tính toán, mở ra ít nhất 3 dòng nghiên cứu mới về vật liệu nano dị thể, khẳng định vị thế học thuật vững chắc và mang lại giá trị thực tiễn lâu dài cho nền khoa học vật liệu tiên tiến.