Tổng quan về luận án

Nghiên cứu vật liệu nano tổ hợp (composite) 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 trong lĩnh vực quang xúc tác dị thể và chuyển hóa năng lượng mặt trời. Hệ vật liệu này đóng vai trò then chốt trong các phản ứng phân hủy chất hữu cơ độc hại từ nguồn nước thải công nghiệp (nhựa, chất tẩy rửa, thuốc trừ sâu), phản ứng quang phân tách nước sinh hydro sạch, các điện cực trong pin mặt trời chất màu nhạy quang (DSSC), cũng như anode pin lithium và cảm biến khí. Tuy nhiên, việc phát triển các mô hình lý thuyết chính xác về tương tác bề mặt giữa graphene và pha tinh thể rutile $\text{TiO}_2$ – dạng thù hình bền vững nhất về mặt nhiệt động của titania – từ lâu đã gặp phải những rào cản mang tính nguyên lý.

Khoảng trống nghiên cứu (research gap) cốt lõi xuất phát từ sự bất đối xứng hình học mạng tinh thể: trong khi bề mặt bền nhất của rutile là (110) sở hữu ô mạng cơ sở hình chữ nhật (rectangular unit cell, $a = 4{,}595\text{ \AA}$, $b = 2{,}959\text{ \AA}$, $\gamma = 90^\circ$), phần lớn các nghiên cứu lý thuyết quốc tế trước đây (như Long et al., 2012; Du et al., 2011) lại ghép nối với ô đơn vị lục giác truyền thống 2 nguyên tử của graphene ($a = b = 2{,}463\text{ \AA}$, $\gamma = 120^\circ$). Sự không tương thích góc mạng ($90^\circ$ so với $120^\circ$) này gây ra biến dạng cưỡng bức (unphysical strain) nghiêm trọng, làm sai lệch cấu trúc electron bề mặt. Bên cạnh đó, thực nghiệm chế tạo composite thực tế thường thu được graphene oxide (GO) hoặc dạng khử graphene oxide (RGO) với các nhóm chức bề mặt chưa xác định được quy luật sắp xếp cụ thể, thay vì graphene nguyên bản.

Để giải quyết triệt để vấn đề trên, luận án của NCS. Trần Thị Thoa, dưới sự hướng dẫn của PGS. Nguyễn Thị Minh Huệ và PGS. Hoàng Văn Hùng tại Trường Đại học Sư phạm Hà Nội (Chuyên ngành Hóa lí thuyết và Hóa lí, Mã số: 9440119), tập trung giải quyết các câu hỏi và giả thuyết nghiên cứu cụ thể:

  1. Câu hỏi nghiên cứu 1 ($Q_1$): Quy luật phân bố không gian và cơ chế tương tác làm bền của các nhóm hydroxyl ($-\text{OH}$) và các nhóm chức đơn ($-\text{NH}_2$, $-\text{CH}_3$, $-\text{OCH}_3$, $-\text{CHO}$, $-\text{COOH}$, epoxy) trên mặt phẳng đáy graphene diễn ra như thế nào?
  2. Câu hỏi nghiên cứu 2 ($Q_2$): Làm thế nào để triệt tiêu sai số tự tương tác (Self-Interaction Error - SIE) của phiếm hàm mật độ tiêu chuẩn khi mô tả cấu trúc dải năng lượng và độ rộng vùng cấm ($E_g$) của rutile $\text{TiO}_2$?
  3. Câu hỏi nghiên cứu 3 ($Q_3$): Xây dựng mô hình siêu ô mạng hai chiều (2D supercell) tiếp xúc dị thể giữa dẫn xuất graphene và rutile $\text{TiO}_2(110)$ với độ bất đối xứng mạng (lattice mismatch) nhỏ nhất nhằm đảm bảo độ chính xác lượng tử?

Khung lý thuyết của luận án được xây dựng vững chắc trên nền tảng cơ học lượng tử thông qua Phương pháp Phiếm hàm Mật độ (Density Functional Theory - DFT), mở rộng với phiếm hàm hiệu chỉnh tương tác phân tán van der Waals (optPBE-vdW) và phương pháp Hubbard hiệu chỉnh năng lượng tương quan mạnh ($\text{DFT}+U$). Đóng góp đột phá của luận án là lần đầu tiên thiết lập quy luật tích tụ para của nhóm hydroxyl tạo vòng lục giác liên kết hydrogen trên graphene, tối ưu hóa đồng thời bộ tham số $(U_d, U_p)$ cho rutile $\text{TiO}_2$, và giải quyết triệt để xung đột hình học mạng tinh thể bằng ô đơn vị chữ nhật 4 nguyên tử của graphene.

Literature Review và Positioning

Lịch sử nghiên cứu composite graphene/$\text{TiO}_2$ ghi nhận hai dòng nghiên cứu chính: thực nghiệm tổng hợp quang xúc tác và mô phỏng hóa học lượng tử tính toán. Về mặt thực nghiệm, quá trình khử GO bằng phương pháp hóa học hay nhiệt phân (Lerf et al., 1998; Dékány et al., 2006) luôn để lại các nhóm chức oxy hóa tồn dư. Các phân tích phổ cộng hưởng từ hạt nhân chất rắn $^{13}\text{C-MAS}$ và 2D $^{13}\text{C-SSNMRS}$ của Weiwei Cai và cộng sự (2008) đã chứng minh cấu trúc GO chứa chủ yếu carbon $sp^2$ (píc $\sim 130\text{ ppm}$), nhóm epoxy (píc $\sim 60\text{ ppm}$) và hydroxyl (píc $\sim 70\text{ ppm}$), trong đó $sp^2$, $-\text{OH}$ và epoxy xen phủ rất gần nhau. Tuy nhiên, phổ NMR thực nghiệm không thể giải mã được trật tự sắp xếp vị trí nguyên tử cục bộ của các nhóm chức này khi mật độ che phủ tăng lên.

Về mặt mô phỏng lý thuyết, các công trình quốc tế chủ yếu tập trung vào composite bề mặt anatase $\text{TiO}_2(101)$ (chẳng hạn Yamashita et al., 2017) vì anatase có cấu trúc mạng dễ thích ứng hơn với graphene. Khi nghiên cứu rutile $\text{TiO}_2(110)$, một cuộc tranh luận học thuật lớn đã nảy sinh giữa hai trường phái:

  • Trường phái ghép nối ô mạng lục giác truyền thống: Long và cộng sự (2012) cùng Du và cộng sự (2011) sử dụng ô đơn vị lục giác 2 nguyên tử của graphene đặt lên bề mặt rutile (110). Hạn chế lớn nhất của cách tiếp cận này là buộc mạng tinh thể phải chịu biến dạng nhân tạo rất lớn do sai lệch góc ($120^\circ$ so với $90^\circ$), dẫn đến sự dịch chuyển sai lệch của mức Fermi và dải dẫn/hóa trị.
  • Trường phái ô mạng chữ nhật 4 nguyên tử: Martsinovich và cộng sự (2017) cùng Suzuki và cộng sự (2017) đề xuất sử dụng ô đơn vị hình chữ nhật của graphene để tương thích với bề mặt tứ giác. Dẫu vậy, các nghiên cứu này chỉ khảo sát graphene nguyên bản hoặc mô hình GO giả định ngẫu nhiên không có cơ sở nhiệt động vững chắc.

Về phía dẫn xuất hydroxyl graphene, Lahaye và cộng sự (2009) nhận định rằng các nhóm $-\text{OH}$ ở cùng một phía không thể gắn vào hai nguyên tử carbon liền kề (vị trí ortho), trong khi Yan và cộng sự (2011) chỉ ra hai nhóm $-\text{OH}$ có xu hướng định vị ở vị trí para. Tuy nhiên, bức tranh tương tác khi số lượng nhóm $-\text{OH}$ tăng từ 3 đến 6 ($G_n\text{OH}$ với $n = 3 \div 6$) vẫn hoàn toàn bị bỏ ngỏ.

Luận án này định vị chính xác tại giao điểm của các khoảng trống học thuật trên: mở rộng nghiên cứu từ ô mạng graphene chữ nhật 4 nguyên tử sang hệ dẫn xuất hydroxyl $G_n\text{OH}$ đa chức, làm sáng tỏ bản chất liên kết hydrogen chuyển dời đỏ và xanh, đồng thời áp dụng phương pháp $\text{DFT}+U$ cải tiến đa kênh ($U_d$ và $U_p$) để mô tả chính xác cả vùng cấm lẫn hằng số mạng rutile $\text{TiO}_2(110)$.

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

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

Luận án đóng góp trực tiếp vào việc mở rộng và kiểm chứng các lý thuyết lượng tử nền tảng:

  • Mở rộng Lý thuyết Phiếm hàm Mật độ cho hệ tương quan mạnh: DFT tiêu chuẩn dựa trên gần đúng mật độ cục bộ (LDA) và gần đúng gradient tổng quát (GGA-PBE) gặp phải sai số tự tương tác (SIE - Self Interaction Error) nghiêm trọng do tính toán giả định mật độ điện tích bị phi định xứ hóa quá mức. Luận án đã mở rộng mô hình Hubbard hiệu chỉnh ($DFT+U$ theo hình thức luận Dudarev) không chỉ dừng lại ở việc áp dụng $U_d$ cho orbital Ti-$3d$ mà thiết lập sự hiệu chỉnh đồng thời $U_p$ cho orbital O-$2p$.
  • Phát triển lý thuyết liên kết Hydrogen đa kênh: Mở rộng lý thuyết liên kết hydrogen cổ điển của Linus Pauling và mô hình chuyển dịch điện tích của Coulson (1957). Luận án chứng minh sự tồn tại đồng thời của liên kết hydrogen chuyển dời đỏ thông thường ($\text{O}-\text{H}\cdots\text{O}$) làm yếu và kéo dài liên kết cộng hóa trị cho proton, cùng liên kết hydrogen chuyển dời xanh không cổ điển ($\text{O}-\text{H}\cdots\pi$) với hệ mây điện tử $\pi$ của khung carbon graphene, đóng vai trò hạ thấp tổng năng lượng hệ thống.
  • Mô hình cấu trúc nhiệt động của Graphene Oxide: Đưa ra quy luật tổ chức hình học phân tử: các nhóm $-\text{OH}$ không phân bố hỗn loạn mà tự tập hợp tại các vị trí para tạo nên cấu trúc cụm vòng 6 hydroxyl hoàn hảo (hexa-hydroxyl ring), định hình lại mô hình cấu trúc bán tinh thể của GO.

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

Khung phân tích của nghiên cứu là sự tích hợp đồng bộ giữa ba trường phái phương pháp tính toán:

  1. Lý thuyết Định lý Bloch và Hàm sóng sóng phẳng: Mô tả hệ tuần hoàn vô hạn thông qua cấu trúc mạng đảo, vùng bất khả quy Brillouin, và phép chiếu hàm sóng tăng cường (PAW - Projector Augmented Wave).
  2. Hiệu chỉnh phân tán bất định xứ phi cục bộ (vdW-DF): Ứng dụng phiếm hàm trao đổi-tương quan optPBE-vdW (tổ hợp $94{,}5268%$ PBE và $5{,}4732%$ revPBE với $\mu = 0$) để nắm bắt chính xác lực phân tán London giữa các lớp vật liệu 2D.
  3. Hình thức luận Hubbard Dudarev: Xây dựng Hamiltonian hiệu dụng $E_{\text{DFT}+U} = E_{\text{DFT}} + E_{\text{Hub}} - E_{\text{dc}}$ với tham số $U_{\text{eff}} = U - J$.

Điều kiện biên và phạm vi áp dụng (boundary conditions) được xác định rõ ràng: nhiệt độ tiệm cận $0\text{ K}$, hệ tuần hoàn 2 chiều với khoảng chân không (vacuum padding) tối thiểu $15\text{ \AA}$ dọc trục $z$ để triệt tiêu tương tác giả giữa các phiến (slabs) lặp lại định kỳ.

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

Thiết kế nghiên cứu

Nghiên cứu tuân theo thế giới quan thực chứng (positivism) và nhận thức luận lượng tử tiên nghiệm (ab initio quantum epistemological stance). Thiết kế nghiên cứu đa tầng (multi-level computational design) được triển khai qua 4 bước:

  1. Cấp độ ô mạng cơ sở (Primitive/Unit cell level): Khảo sát đối chuẩn song song ô mạng lục giác 2 nguyên tử và ô mạng chữ nhật 4 nguyên tử của graphene tinh khiết; ô mạng tứ giác khối của rutile $\text{TiO}_2$.
  2. Cấp độ dẫn xuất chức hóa (Functionalized Monolayer level): Khảo sát các cấu hình $G_1\text{F}$ và chuỗi dẫn xuất $G_n\text{OH}$ ($n = 1 \div 6$) với các vị trí thế ortho, meta, para.
  3. Cấp độ bề mặt tinh thể (Surface Slab level): Cắt và tối ưu hóa bề mặt rutile $\text{TiO}_2(110)$ theo 3 mô hình thư giãn khác nhau với độ dày thay đổi từ 1 đến 7 lớp nguyên tử.
  4. Cấp độ tiếp xúc dị thể (Heterojunction Interface level): Xây dựng siêu ô mạng composite 2D bằng cách ghép các bội số nguyên $(m \times n)$ của graphene chữ nhật với $(p \times q)$ của bề mặt rutile (110).

Quy trình nghiên cứu rigorous

Toàn bộ các tính toán hóa học lượng tử được thực hiện bằng phần mềm mô phỏng vật liệu chuyên dụng Vienna Ab-initio Simulation Package (VASP), áp dụng phương pháp sóng phẳng tăng cường tuyến tính (PAW).

Quy trình kiểm soát chất lượng và độ hội quy được thiết lập nghiêm ngặt:

  • Năng lượng cắt sóng phẳng (Plane-wave cutoff energy): $E_{\text{cut}} = 500\text{ eV}$, đảm bảo mô tả chính xác cả orbital định vị cao như C-$2p$, O-$2p$ và Ti-$3d$.
  • Lấy mẫu mạng đảo (k-point sampling): Sử dụng lưới Monkhorst-Pack với mật độ dày đặc: $15 \times 15 \times 1$ cho ô đơn vị graphene; $9 \times 9 \times 13$ cho khối rutile $\text{TiO}_2$; và $5 \times 5 \times 1$ cho các siêu ô mạng slab bề mặt và composite.
  • Tiêu chuẩn hội tụ (Convergence criteria): Năng lượng điện tử tự nhất quán (SCF) đạt độ chính xác dưới $10^{-5}\text{ eV}/\text{cell}$; lực tác dụng lên từng nguyên tử trong quá trình tối ưu hóa hình học (relaxation) nhỏ hơn $0{,}01\text{ eV/\AA}$.
  • Mô hình hóa bề mặt Slab: Để đánh giá ảnh hưởng của hiệu ứng giam hãm lượng tử và tương tác bề mặt, 3 mô hình slab rutile (110) được khảo sát:
    • FR (Full Relax): Cho phép tất cả các lớp nguyên tử tự do dịch chuyển.
    • FIL (Fix Inner Layer): Cố định các lớp nguyên tử ở tâm, chỉ giải phóng các lớp bề mặt ngoài cùng.
    • F2B (Fix Two Bottom): Cố định 2 lớp đáy dưới cùng và tối ưu các lớp trên.

Data và phân tích

Phân tích dữ liệu lượng tử bao gồm:

  • Mật độ trạng thái tổng (DOS - Density of States) và mật độ trạng thái riêng phần (PDOS - Partial Density of States) được phân tích theo từng orbital ($s, p_x, p_y, p_z, d_{xy}, d_{yz}, d_{z^2}, \dots$).
  • Cấu trúc dải năng lượng (Electronic Band Structure) được chiếu dọc theo các đường đối xứng cao trong vùng Brillouin thứ nhất ($\Gamma - \text{X} - \text{M} - \text{Y} - \Gamma$).
  • Năng lượng liên kết ($E_b$) trên từng nhóm chức được định lượng qua công thức: $$E_b = \frac{E_{\text{total}}(G_n\text{OH}) - E_{\text{total}}(\text{graphene}) - n \cdot E(\text{OH})}{n}$$
  • Phân tích hiệu ứng chênh lệch mật độ electron (Electron Density Difference - EDD) tại các mặt cắt đồng mức đẳng trị $0{,}0015\text{ e/\AA}^3$ để làm rõ sự tái phân bố điện tích và chuyển dịch electron quang sinh.
  • Mô phỏng phổ nhiễu xạ tia X (XRD) dựa trên định luật Bragg ($n\lambda = 2d\sin\theta$) để đối chiếu trực tiếp hằng số mạng và góc nhiễu xạ với dữ liệu thực nghiệm tinh thể học.

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

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

Nghiên cứu đã mang lại 5 phát hiện mang tính đột phá về cấu trúc vi mô và tính chất điện tử:

  1. Quy luật hình thành cụm Hexa-hydroxyl Para và mạng liên kết Hydrogen kép: Trích dẫn nguyên văn từ luận án:

"Trong dẫn xuất GnOH, 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 ∙∙∙ O và các liên kết hydrogen chuyển dời xanh O-H ∙∙∙ π giữa các nhóm hydroxyl." Khoảng cách liên kết cộng hóa trị $\text{C}-\text{O}$ trung bình đạt $1{,}47\text{ \AA}$, khoảng cách liên kết hydrogen $\text{O}-\text{H}\cdots\text{O}$ dao động từ $1{,}78\text{ \AA}$ đến $1{,}95\text{ \AA}$ với góc liên kết gần như tuyến tính ($160^\circ - 175^\circ$), làm tăng độ bền nhiệt động của dẫn xuất $G_6\text{OH}$ vượt trội so với các đồng phân phân tán.

  1. Trật tự độ bền nhiệt động và trạng thái Fermi của dẫn xuất $G_1\text{F}$: Độ bền liên kết của các nhóm chức đơn tuân theo thứ tự thực nghiệm khắt khe: $$G_{\text{epo}} > G_{\text{CH}3} > G{\text{NH}2} > G{\text{OCH}3} > G{\text{CHO}} \gg G_{\text{COOH}}$$ Nhóm carboxyl ($-\text{COOH}$) liên kết yếu nhất trên mặt phẳng đáy graphene, giải thích vì sao trong thực nghiệm phổ NMR, nhóm $-\text{COOH}$ chỉ xuất hiện tại các mép khuyết tật. Đặc biệt, các trạng thái điện tử xung quanh mức Fermi của $G_1\text{F}$ bị chi phối chủ yếu bởi orbital $2p_z$ của nguyên tử dị tố $X$ (liên kết trực tiếp với carbon graphene).

  2. Xác lập giá trị tham số Hubbard tối ưu đa kênh $(U_d, U_p)$: Trích dẫn từ kết quả tính toán rutile $\text{TiO}_2$:

"Đối với rutile TiO2: Up =10 eV là trị số tối ưu. Sự kết hợp của Ud và Up tối ưu (Up = 10 eV) không chỉ làm giảm sai số trong hằng số mạng mà còn làm tăng Eg nhiều hơn so với sự hiệu chỉnh của Ud." Khi tính toán DFT tiêu chuẩn, độ rộng vùng cấm rutile bị đánh giá thấp nghiêm trọng ($E_g = 1{,}85\text{ eV}$ so với thực nghiệm $3{,}03\text{ eV}$ ở $10\text{ K}$). Khi áp dụng $U_d = 7\text{ eV}$ và $U_p = 10\text{ eV}$, $E_g$ đạt chuẩn $3{,}05\text{ eV}$ đồng thời sai số hằng số mạng $a = 4{,}59\text{ \AA}$ và $c = 2{,}95\text{ \AA}$ giảm xuống dưới $1{,}2%$.

  1. Khử dao động chẵn - lẻ của tính chất bề mặt Rutile (110): Khảo sát các slab rutile (110) cho thấy năng lượng bề mặt và độ rộng vùng cấm biến thiên dao động dạng sóng chẵn-lẻ (odd-even oscillation) theo số lớp nguyên tử. Phương pháp $\text{DFT}+U$ giúp ổn định hóa và triệt tiêu dao động này khi đạt từ 5 lớp nguyên tử trở lên. Mô hình F2B tuy hội tụ nhanh nhưng gây sai số $E_g$ cục bộ lớn, do đó mô hình FR 5 lớp được chứng minh là cấu hình tối ưu nhất.

  2. Mô hình tiếp xúc dị thể hai chiều không biến dạng: Bằng việc thiết lập siêu ô mạng hình chữ nhật kết hợp giữa graphene ($m \times n$) và rutile (110) ($p \times q$), luận án đã loại bỏ hoàn toàn sai lệch góc $30^\circ$, hạ độ lệch mạng dọc các trục xuống dưới $2{,}5%$, tạo tiền đề mô tả chính xác sự dịch chuyển điện tích quang sinh qua mặt phân giới.

Implications đa chiều

  • Về mặt lý thuyết: Mở ra mô hình vi mô hoàn chỉnh giải thích cơ chế hình thành các "đảo chức hóa" (functionalized domains) trên bề mặt graphene oxide. Cung cấp dữ liệu chuẩn mực về đóng góp orbital cho các nhà hóa học lượng tử nghiên cứu hệ vật liệu lai tạp $sp^2-sp^3$.
  • Về mặt phương pháp luận: Thiết lập quy trình tính toán chuẩn kết hợp optPBE-vdW và đa tham số $\text{DFT}+(U_d, U_p)$, có thể chuyển giao trực tiếp để nghiên cứu các hệ oxide kim loại chuyển tiếp khác như $\text{ZnO}, \text{Fe}_2\text{O}_3, \text{WO}_3$ tiếp xúc với vật liệu 2D (chalcogenide, MXene).
  • Về mặt thực tiễn và công nghệ: Nhóm hydroxyl đóng vai trò "cầu nối dẫn truyền điện tích" (charge transfer bridge). Các vòng hexa-hydroxyl para làm giảm rào cản chuyển dịch electron từ $\text{TiO}_2$ sang graphene, ngăn chặn hiện tượng tái tổ hợp electron - lỗ trống ($\text{e}^-/\text{h}^+$). Đây là chỉ dẫn then chốt để các kỹ sư vật liệu tối ưu hóa quy trình khử GO (kiểm soát mật độ $-\text{OH}$) nhằm nâng cao hiệu suất lượng tử trong pin mặt trời nhạy quang (DSSC) và xúc tác phân hủy ô nhiễm nước.

Limitations và Future Research

Luận án thẳng thắn thừa nhận các giới hạn mang tính phương pháp:

  1. Giới hạn nhiệt động tại $0\text{ K}$: Toàn bộ các tính toán DFT được thực hiện trong môi trường chân không ở trạng thái tĩnh tuyệt đối ($T = 0\text{ K}$), chưa tính đến các hiệu ứng động học do dao động mạng (phonons) ở nhiệt độ phòng.
  2. Môi trường chân không lý tưởng: Chưa tích hợp mô hình dung môi tường minh (explicit solvent model) hoặc liên tục (implicit solvation) mô phỏng pha lỏng của môi trường phản ứng xử lý nước thải.
  3. Quy mô kích thước siêu ô mạng: Do giới hạn tài nguyên siêu máy tính, nồng độ nhóm chức hóa mới dừng lại ở $n = 6$ trên siêu ô mạng cơ sở, chưa thể mở rộng đến mô hình bất trật tự hoàn toàn trên diện tích micromet.

Chương trình nghiên cứu tương lai (Future Research Agenda) mở ra 4 hướng đi:

  • Ứng dụng mô phỏng Động lực học phân tử Tiên nghiệm (Ab Initio Molecular Dynamics - AIMD) tại nhiệt độ $300\text{ K} - 500\text{ K}$ để theo dõi quá trình phân rã và tái liên kết của các nhóm $-\text{OH}$.
  • Tính toán đường phản ứng xúc tác chi tiết (Nudged Elastic Band - NEB) cho phản ứng quang phân tách nước ($\text{H}_2\text{O} \rightarrow \text{H}_2 + \frac{1}{2}\text{O}_2$) trực tiếp trên bề mặt composite.
  • Khảo sát sự kết hợp đa nhóm chức đồng thời (ví dụ: hydroxyl kết hợp epoxy và amine) trên cùng một phiến graphene.
  • Tính toán tính chất truyền dẫn điện tử phi cân bằng (NEGF-DFT) để dự đoán đường đặc trưng dòng - thế (I-V) của linh kiện cảm biến khí.

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

  • Tác động học thuật: Cung cấp hệ thống dữ liệu tham chiếu chuẩn xác về cấu trúc dải, mật độ trạng thái (DOS/PDOS), và thông số mạng tinh thể cho cộng đồng nghiên cứu hóa lý thuyết tính toán. Dự báo khả năng thu hút trích dẫn cao trong các bài báo quốc tế thuộc danh mục ISI/Scopus chuyên ngành vật liệu quang điện tử.
  • Chuyển đổi công nghiệp: Định hướng cho các doanh nghiệp vật liệu và công nghệ môi trường phương pháp kiểm soát hình thái bề mặt nano RGO/$\text{TiO}_2$, từ đó hạ giá thành sản xuất tấm pin mặt trời và màng lọc quang xúc tác xử lý nước công nghiệp dệt nhuộm, hóa chất.
  • Lợi ích xã hội & Môi trường: Đóng góp trực tiếp vào mục tiêu trung hòa carbon và phát triển bền vững thông qua việc thúc đẩy giải pháp công nghệ năng lượng tái tạo (hydro xanh, pin mặt trời thế hệ mới).

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

  • Nghiên cứu sinh và Giảng viên Hóa học / Vật lý chất rắn: Tiếp cận quy trình phương pháp luận mẫu mực về xử lý sai số SIE trong $\text{DFT}+U$ và kỹ thuật chuyển đổi ô mạng chữ nhật 4 nguyên tử cho vật liệu 2D.
  • Các nhà nghiên cứu R&D vật liệu nano: Nắm bắt được nguyên lý cấu trúc vi mô để kiểm soát mức độ chức hóa bề mặt graphene oxide trong thực nghiệm tổng hợp.
  • Kỹ sư công nghệ năng lượng sạch: Ứng dụng các thông số về độ rộng vùng cấm và vị trí dải năng lượng (CBM, VBM) để thiết kế các tiếp giáp dị thể quang điện đạt hiệu suất thu nhận photon tối ưu.

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à phát hiện ra sự tự tập hợp của các nhóm hydroxyl thành các vòng hexa-hydroxyl ở vị trí para trên một phía mặt phẳng graphene, được ổn định hóa bởi cơ chế liên kết hydrogen kép (chuyển dời đỏ $\text{O}-\text{H}\cdots\text{O}$ và chuyển dời xanh $\text{O}-\text{H}\cdots\pi$). Phát hiện này mở rộng trực tiếp Lý thuyết cấu trúc Graphene Oxide của Lerf-Klinowski và Dékány, biến mô hình phân bố ngẫu nhiên thành mô hình có trật tự hình học lượng tử bền vững.

2. Đột phá về phương pháp luận so với các nghiên cứu quốc tế trước đây là gì? So sánh với nghiên cứu của Long et al. (2012) và Du et al. (2011) vốn ghép nối ô mạng lục giác của graphene gây sai lệch góc $30^\circ$, luận án đã sử dụng thành công ô mạng hình chữ nhật 4 nguyên tử ($a_3 = 2{,}463\text{ \AA}$, $a_4 = 4{,}266\text{ \AA}$) tương thích hoàn hảo với mạng tứ giác rutile $\text{TiO}_2(110)$. Đồng thời, thay vì chỉ tối ưu hóa đơn lẻ tham số $U_d$ cho Ti-$3d$, luận án thiết lập phương pháp tối ưu hóa đồng thời $(U_d = 7\text{ eV}, U_p = 10\text{ eV})$, tái tạo chính xác đồng thời cả hằng số mạng và độ rộng vùng cấm thực nghiệm ($3{,}03\text{ eV}$).

3. Phát hiện nào gây bất ngờ nhất và có bằng chứng số liệu ra sao? Kết quả bất ngờ nhất là nhóm carboxyl ($-\text{COOH}$) – vốn được xem là nhóm chức quan trọng trong hóa học hữu cơ – lại có năng lượng liên kết trên mặt phẳng đáy graphene kém bền nhất trong chuỗi $G_1\text{F}$ ($E_b$ thấp hơn rất nhiều so với epoxy và hydroxyl). Điều này cung cấp bằng chứng lượng tử giải thích vì sao phổ cộng hưởng từ hạt nhân $^{13}\text{C-SSNMRS}$ thực nghiệm của Weiwei Cai et al. chỉ ghi nhận các píc $\text{C}=\text{O}$ cường độ rất yếu tại $169$ và $193\text{ ppm}$ và hoàn toàn không có sự xen phủ không gian với carbon $sp^2$ mặt đáy.

4. Quy trình tái lập (Replication Protocol) có được cung cấp chi tiết không? Có. Luận án cung cấp đầy đủ và chi tiết các tọa độ mạng đảo, ma trận biến đổi vector mạng thực sang mạng đảo, bán kính cắt thế giả PAW, năng lượng cắt $500\text{ eV}$, lưới điểm K-points, thuật toán giải lặp tự nhất quán Kohn-Sham, và các tham số hiệu chỉnh optPBE-vdW giúp mọi phòng thí nghiệm tính toán lượng tử đều có thể tái lập chính xác $100%$ kết quả.

5. Chương trình nghiên cứu 10 năm tiếp theo được vạch ra như thế nào? Chương trình nghiên cứu hướng tới: (1) Khảo sát động học phản ứng quang xúc tác phân rã chất ô nhiễm bằng AIMD ở điều kiện nhiệt động thực; (2) Mô phỏng cơ chế chuyển dịch điện tích cực nhanh (ultrafast charge transfer dynamics) tại mặt phân giới composite dưới kích thích quang học; (3) Tích hợp học máy (Machine Learning Potentials) để dự đoán cấu trúc của các dẫn xuất graphene phức tạp chứa hàng vạn nguyên tử.

Kết luận

Luận án tiến sĩ của NCS. Trần Thị Thoa đã giải quyết trọn vẹn một bài toán học thuật hóc búa trong lĩnh vực hóa lý thuyết và vật liệu nano tính toán:

  1. Thành tựu 1: Làm sáng tỏ cấu trúc hình học và tính chất điện tử của các dẫn xuất hydroxyl graphene $G_n\text{OH}$ ($n = 1 \div 6$), chỉ ra xu hướng liên kết para tạo vòng hexa-hydroxyl bền vững nhờ mạng lưới liên kết hydrogen kép.
  2. Thành tựu 2: Xác lập thứ tự độ bền tương đối của các nhóm chức đơn $G_1\text{F}$ và chỉ ra vai trò quyết định của orbital dị tố $2p_z$ đối với các trạng thái điện tử tại mức Fermi.
  3. Thành tựu 3: Giải quyết triệt để sai số tự tương tác SIE của phiếm hàm DFT trên rutile $\text{TiO}_2$ bằng việc xác lập bộ tham số hiệu chỉnh tối ưu $(U_d = 7\text{ eV}, U_p = 10\text{ eV})$.
  4. Thành tựu 4: Khảo sát toàn diện 3 mô hình slab bề mặt rutile (110) (FR, FIL, F2B), chứng minh năng lực của $\text{DFT}+U$ trong việc loại bỏ dao động chẵn-lẻ của năng lượng bề mặt và xác định mô hình FR 5 lớp là tối ưu.
  5. Thành tựu 5: Xây dựng thành công mô hình composite 2D giữa graphene/dẫn xuất graphene và rutile $\text{TiO}_2(110)$ trên cơ sở ô đơn vị hình chữ nhật 4 nguyên tử, triệt tiêu biến dạng mạng và mở đường cho các mô phỏng quang xúc tác tiếp giáp dị thể chuẩn xác cao.

Công trình không chỉ nâng cao vị thế của nghiên cứu hóa học tính toán trong nước mà còn đóng góp dữ liệu và phương pháp luận có giá trị lâu dài cho cộng đồng khoa học vật liệu quốc tế.