Tổng quan về luận án

Vật lý chất rắn và khoa học vật liệu nano đương đại đang chứng kiến bước chuyển mình mạnh mẽ từ các cấu trúc ba chiều (3D) sang hệ vật liệu phân lớp hai chiều (2D). Kể từ dấu mốc lịch sử năm 2004 khi graphene được phân tách thành công từ graphite bởi Novoselov và Geim, một kỷ nguyên nghiên cứu chuyên sâu về các vật liệu nguyên tử phẳng đơn lớp đã được thiết lập. Mặc dù graphene sở hữu những đặc tính cơ học phi thường (độ bền kéo nội tại đạt 130 GPa, mô-đun Young xấp xỉ 1 TPa) cùng độ linh động hạt tải vượt trội ($>15.000\text{ cm}^2\text{V}^{-1}\text{s}^{-1}$), việc thiếu hụt một vùng cấm tự nhiên (zero-bandgap semimetal) đã tạo ra rào cản kỹ thuật nghiêm trọng trong việc chế tạo các linh kiện chuyển mạch logic bán dẫn như Transistor hiệu ứng trường (FET). Tỷ lệ đóng/ngắt ($I_{on}/I_{off}$) thấp khiến graphene khó đáp ứng yêu cầu vi điện tử hiện đại, trong khi các phương pháp can thiệp ngoài như biến dạng cơ học hay pha tạp thường phá vỡ đối xứng mạng tinh thể và suy giảm độ linh động.

Trước khoảng trống nghiên cứu (research gap) cốt lõi đó, cộng đồng vật lý lượng tử đã chuyển hướng khảo sát các hợp chất chalcogenide hai chiều, đặc biệt là nhóm monochalcogenide kim loại nhóm III ($MX$, với $M = \text{Ga, In}$ và $X = \text{S, Se, Te}$). Luận án tiến sĩ chuyên ngành Vật lý lý thuyết và Vật lý toán của tác giả Vũ Thị Tuyết Vi, dưới sự hướng dẫn khoa học của PGS.TS. Nguyễn Ngọc Hiếu và PGS.TS. Bùi Đình Hợi tại Trường Đại học Sư phạm – Đại học Huế, mang tựa đề "Nghiên cứu tính chất điện tử của vật liệu hai chiều monochalcogenide" (Mã số: 9 44 01 03), đã tiên phong giải quyết các câu hỏi nền tảng về cấu trúc nguyên tử, độ bền nhiệt động, cấu trúc dải năng lượng và cơ chế phản ứng dưới tác động trường ngoài của họ vật liệu đơn lớp $MX$ đối xứng và các dị cấu trúc bất đối xứng Janus ($M_2XY$ và $\text{GaIn}XO$).

Các câu hỏi và giả thuyết nghiên cứu chính được xác lập bao gồm:

  1. Câu hỏi nghiên cứu 1 (RQ1): Cấu trúc hình học tối ưu, phổ phonon và đặc tính đàn hồi của các đơn lớp monochalcogenide $MX$ nhóm III ($M = \text{Ga, In}$; $X = \text{S, Se, Te}$) thể hiện độ bền cơ học và ổn định nhiệt động lực học như thế nào theo tiêu chuẩn Born-Huang?
  2. Câu hỏi nghiên cứu 2 (RQ2): Việc phá vỡ tính đối xứng gương theo phương thẳng đứng thông qua việc kiến tạo cấu trúc đơn lớp Janus $M_2XY$ ($X \neq Y$) và dị hợp hợp kim Janus bậc bốn $\text{GaIn}XO$ sẽ làm biến đổi dải năng lượng, độ rộng vùng cấm ($E_g$) và phân cực điện nội tại ra sao?
  3. Giả thuyết nghiên cứu 1 (H1): Biến dạng cơ học hai trục ($\varepsilon_{xy}$) và điện trường ngoài ($E$) vuông góc với mặt phẳng sẽ kích hoạt quá trình chuyển pha điện tử từ bán dẫn vùng cấm gián tiếp sang vùng cấm trực tiếp, thậm chí gây ra quá trình chuyển tiếp bán dẫn – kim loại (semiconductor-to-metal transition).
  4. Giả thuyết nghiên cứu 2 (H2): Cấu trúc dị hợp Janus $\text{GaIn}XO$ với sự bất đối xứng ở cả phân lớp kim loại và chalcogenide sẽ kiến tạo nên gradient thế tĩnh điện cục bộ mạnh, làm phân tách thế hóa học hai bề mặt và gia tăng đáng kể độ linh động hạt tải định hướng.

Khung lý thuyết của công trình được định vị vững chắc trên nền tảng cơ học lượng tử nhiều hạt, cụ thể là Lý thuyết Phiếm hàm Mật độ (Density Functional Theory – DFT) với các phép gần đúng gradient suy rộng (GGA-PBE), phiếm hàm lai Heyd-Scuseria-Ernzerhof (HSE06), kết hợp hiệu chỉnh tương tác phân tán van der Waals DFT-D2 của Grimme và tương tác spin-quỹ đạo (Spin-Orbit Coupling – SOC). Phạm vi khảo sát bao quát toàn diện 6 đơn lớp $MX$, 6 cấu trúc Janus tam nguyên $M_2XY$, và 6 cấu hình sắp xếp nguyên tử của Janus tứ nguyên $\text{GaIn}XO$. Công trình đã cung cấp dữ liệu định lượng chính xác về cấu trúc dải năng lượng, khối lượng hiệu dụng hạt tải ($m^*$), mô-đun đàn hồi 2D ($C_{2D}$), hằng số thế biến dạng ($E_d$) và độ linh động ($\mu$), tạo tiền đề lý thuyết chuẩn xác cho kỹ nghệ linh kiện quang điện tử, cảm biến khí và xúc tác quang phân tách nước.

Literature Review và Positioning

Vật lý vật liệu 2D đã trải qua các bước phát triển quan trọng từ mô hình đơn nguyên tử phẳng graphene (Wallace, 1947; Novoselov et al., 2004) đến silicene (Vogt et al., 2012), germanene (Bianco et al., 2013) và stanene (Zhu et al., 2015). Tiếp đó, nhóm vật liệu Transition Metal Dichalcogenides – TMDs dạng $MX_2$ (như $\text{MoS}_2$, $\text{WS}_2$, $\text{MoSe}_2$, $\text{WSe}_2$) đã thu hút hàng chục nghìn nghiên cứu nhờ sự xuất hiện của vùng cấm trực tiếp ở giới hạn đơn lớp (Mak et al., 2010; Splendiani et al., 2010). Tuy nhiên, TMDs thường gặp giới hạn về độ linh động hạt tải thực tế ở nhiệt độ phòng và hiệu ứng phân cực nội tại yếu do cấu trúc đối xứng phẳng hoàn hảo.

Nhánh nghiên cứu về các monochalcogenide kim loại (MCs) nhóm IV (như $\text{SnS}$, $\text{SnSe}$, $\text{GeS}$, $\text{GeSe}$) và nhóm III ($MX$) đã mở ra hướng đi đột phá nhờ tính chất dị hướng quang học và hệ số nhiệt điện ấn tượng (Fei et al., 2015). Các thực nghiệm tổng hợp vật liệu gần đây đã chứng minh khả năng chế tạo thành công các phiến nano và màng đơn lớp $\text{GaS}$ (Late et al., 2012), $\text{GaSe}$ (Zhou et al., 2014), $\text{GaTe}$ (Liu et al., 2014), $\text{InS}$ và $\text{InSe}$ (Bandurin et al., 2017). Mặc dù các công trình thực nghiệm bước đầu ghi nhận độ phát quang và khả năng ứng dụng trong dò quang phổ, bức tranh vi mô hoàn chỉnh về sự tiến hóa cấu trúc điện tử dưới các trường lực cơ – điện hỗn hợp vẫn tồn tại nhiều tranh luận khoa học chưa được giải quyết dứt điểm.

Dòng nghiên cứu / Tác giả Đối tượng & Phương pháp Phát hiện chính Khoảng trống & Hạn chế
Dòng TMDs kinh điển (Radisavljevic et al., 2011; Chhowalla et al., 2013) Đơn lớp $\text{MoS}_2$, $\text{WS}_2$; Thực nghiệm FET & DFT-PBE Độ hở vùng cấm trực tiếp 1.8–2.1 eV; Tỷ lệ $I_{on}/I_{off} \sim 10^8$. Thiếu phân cực bề mặt lưỡng cực nội tại; suy giảm tính chất quang khi chịu biến dạng kéo lớn.
Dòng Monochalcogenide nhóm III đối xứng (Zólyomi et al., 2013; Late et al., 2012) Đơn lớp $\text{GaS}$, $\text{GaSe}$, $\text{InSe}$; DFT bán cổ điển Dải hóa trị dạng "sombrero" (mũ Mexico); vùng cấm gián tiếp rộng. Chưa khảo sát có hệ thống ảnh hưởng của hiệu ứng spin-quỹ đạo (SOC) phối hợp cùng biến dạng cơ học đa trục.
Dòng Janus TMDs đột phá (Lu et al., 2017; Zhang et al., 2017) Janus $\text{MoSSe}$; Bóc tách CVD & Tổng hợp hóa học Phá vỡ đối xứng gương $z \to -z$; xuất hiện momen lưỡng cực vuông góc lớn. Chỉ tập trung vào kim loại chuyển tiếp nhóm VIB; chưa khai thác thế tĩnh điện bất đối xứng trên nền monochalcogenide nhóm III.
Nghiên cứu của Luận án (Vũ Thị Tuyết Vi et al., 2022) Đơn lớp $MX$, Janus $M_2XY$, Janus $\text{GaIn}XO$; DFT-PBE, HSE06, SOC, AIMD Xác lập quy luật biến điệu dải năng lượng toàn diện; phát hiện độ linh động hạt tải siêu cao trên Janus $\text{GaIn}XO$; kiểm soát chuyển pha bán dẫn – kim loại bằng biến dạng và điện trường. Giải quyết triệt để khoảng trống lý thuyết về họ monochalcogenide nhóm III và vật liệu Janus dị nhân tố.

Trong bức tranh tổng thể, luận án đã định vị tại giao điểm tiên phong giữa vật lý lượng tử chất rắn và kỹ nghệ cấu trúc nano: nghiên cứu chuyển tiếp từ các đơn lớp đối xứng $D_{3h}$ sang các cấu trúc phá vỡ đối xứng gương $C_{3v}$ ($M_2XY$ và $\text{GaIn}XO$). Bằng việc so sánh đa chiều với các nghiên cứu quốc tế chuẩn mực (như các tính toán của Demchenko et al. trên $\text{InSe}$ và Dong et al. trên vật liệu phân cực Janus), luận án đã chỉ ra rằng việc thay thế nguyên tử dị chủng ở hai bề mặt tạo ra một trường điện nội tại $\Delta \Phi$, dẫn đến sự phân tách mạnh của các dải năng lượng và thay đổi hoàn toàn động học dịch chuyển của electron/lỗ trống.

Đó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à hoàn thiện hệ thống luận điểm của lý thuyết dải năng lượng lượng tử trong vật liệu hai chiều thông qua các đóng góp then chốt:

                  ┌─────────────────────────────────────────────────────────┐
                  │    ĐƠN LỚP ĐỐI XỨNG MX (Nhóm không gian D3h)             │
                  │    - Cấu trúc 4 lớp nguyên tử: X - M - M - X            │
                  │    - Momen lưỡng cực vuông góc: Pz = 0                  │
                  │    - Vùng cấm gián tiếp (VBM nằm giữa Γ-K, CBM tại Γ)   │
                  └────────────────────────────┬────────────────────────────┘
                                               │ Thay thế nguyên tử mặt ngoài (X -> Y)
                                               ▼
                  ┌─────────────────────────────────────────────────────────┐
                  │    JANUS TAM NGUYÊN M2XY (Nhóm không gian C3v)           │
                  │    - Cấu trúc: X - M - M - Y (Phá vỡ đối xứng z -> -z)   │
                  │    - Xuất hiện phân cực nội tại: Pz ≠ 0, ΔΦ > 0          │
                  │    - Tách mức spin-quỹ đạo (SOC splitting) tại VBM/CBM   │
                  └────────────────────────────┬────────────────────────────┘
                                               │ Lai hóa cation dị chủng (M -> M1, M2)
                                               ▼
                  ┌─────────────────────────────────────────────────────────┐
                  │    JANUS BẬC BỐN GaInXO (Cấu trúc bất đối xứng kép)     │
                  │    - Cấu hình: X - Ga - In - O hoặc O - Ga - In - X     │
                  │    - Độ lệch thế tĩnh điện khổng lồ giữa 2 bề mặt        │
                  │    - Độ linh động electron/lỗ trống dị hướng siêu cao    │
                  └─────────────────────────────────────────────────────────┘
  1. Thuyết minh cơ chế phá vỡ đối xứng gương ($z \to -z$): Chứng minh bằng toán lượng tử rằng việc chuyển từ cấu trúc $D_{3h}$ sang $C_{3v}$ làm triệt tiêu tính suy biến không gian của các orbital nguyên tử $p_z$, dẫn đến sự xuất hiện của độ chênh lệch thế tĩnh điện chân không ($\Delta \Phi = \Phi_X - \Phi_Y$) giữa hai bề mặt vật liệu.
  2. Quy luật biến thiên năng lượng dải: Khái quát hóa mối quan hệ tỷ lệ nghịch giữa bán kính nguyên tử chalcogenide ($r_S < r_{Se} < r_{Te}$) với độ rộng vùng cấm ($E_g$). Khi đi từ hợp chất chứa S sang Te, liên kết cộng hóa trị $M-X$ dãn dài làm giảm mức độ xen phủ orbital $\sigma(M-M)$ và $p-p$, trực tiếp làm thu hẹp $E_g$.
  3. Hiệu ứng tách mức Spin-Quỹ đạo (Rashba-like SOC splitting): Xác định sự kết hợp giữa gradien thế nội tại bất đối xứng và tương tác spin-quỹ đạo mạnh của các nguyên tử nặng (như In, Te) tạo ra sự tách mức năng lượng đáng kể tại các điểm đối xứng cao $\Gamma$ và $K$ trong vùng Brillouin, làm giàu thêm cơ sở lý thuyết cho lĩnh vực Spintronics.

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

Khung phân tích của luận án tích hợp đa tầng ba lý thuyết vật lý nền tảng:

  • Lý thuyết cấu trúc dải năng lượng tinh thể (Band Theory): Phân tích hình thái vùng dẫn (CBM) và vùng hóa trị (VBM), phân bố mật độ trạng thái riêng phần (PDOS) để giải mã đóng góp của các orbital $s$, $p_x, p_y$, và $p_z$.
  • Lý thuyết nhiễu loạn phiếm hàm mật độ (DFPT): Khảo sát cấu trúc dao động mạng, ma trận lực động lực học và các nhánh phonon âm (LA, TA, ZA) cùng các nhánh phonon quang nhằm xác định giới hạn ổn định cơ học.
  • Lý thuyết thế biến dạng (Deformation Potential Theory của Bardeen-Shockley): Mô hình hóa định lượng độ linh động hạt tải 2D ($\mu_{2D}$) dựa trên sự tương tác giữa electron/lỗ trống với các phonon âm tán xạ: $$\mu_{2D} = \frac{e\hbar^3 C_{2D}}{k_B T m^* m_d (E_d)^2}$$ Trong đó $e$ là điện tích nguyên tố, $\hbar$ là hằng số Planck rút gọn, $C_{2D}$ là mô-đun đàn hồi phẳng hai chiều, $k_B$ là hằng số Boltzmann, $T$ là nhiệt độ tuyệt đối ($300\text{ K}$), $m^$ là khối lượng hiệu dụng theo phương truyền dẫn, $m_d = \sqrt{m_x^ m_y^*}$ là khối lượng hiệu dụng trung bình hình học, và $E_d$ là hằng số thế biến dạng bậc một dọc theo trục chuyển dời.

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

Thiết kế nghiên cứu

Nghiên cứu được thiết kế hoàn toàn theo triết lý định lượng chuẩn xác từ nguyên lý ban đầu (Ab-Initio First-Principles Calculations), vận hành trên hệ thống phần mềm tính toán lượng tử mã nguồn mở Quantum ESPRESSO. Thiết kế đa cấp độ (multi-level design) bao gồm:

                            QUY TRÌNH MÔ PHỎNG DFT ĐA TẦNG
                                          │
            ┌─────────────────────────────┴─────────────────────────────┐
            ▼                                                           ▼
┌───────────────────────┐                                   ┌───────────────────────┐
│ TỐI ƯU HÓA HÌNH HỌC   │                                   │ KIỂM CHỨNG ĐỘ BỀN     │
│ - Tọa độ nguyên tử    │                                   │ - Tiêu chuẩn Born     │
│ - Thông số mạng a     │                                   │ - Phổ phonon (DFPT)   │
│ - Tiêu chuẩn hội tụ:  │                                   │ - Động học AIMD (300K)│
│   ΔE < 10⁻⁶ eV        │                                   │                       │
│   F < 10⁻³ eV/Å       │                                   │                       │
└───────────┬───────────┘                                   └───────────┬───────────┘
            │                                                           │
            └─────────────────────────────┬─────────────────────────────┘
                                          ▼
┌───────────────────────────────────────────────────────────────────────────────────┐
│ TÍNH TOÁN CẤU TRÚC ĐIỆN TỬ & ĐẶC TÍNH ĐỘNG HỌC HẠT TẢI                            │
│ - Năng lượng dải (GGA-PBE, Hybrid HSE06, PBE+SOC)                                 │
│ - Mật độ trạng thái (DOS / PDOS)                                                  │
│ - Thế tĩnh điện & Công thoát (Evac, EF, Φ)                                        │
│ - Khối lượng hiệu dụng (m*), Thế biến dạng (Ed) & Độ linh động hạt tải (μ)        │
└───────────────────────────────────────────────────────────────────────────────────┘

Quy trình nghiên cứu rigorous

  1. Mô hình hóa ô mạng đơn vị và chân không: Thiết lập ô mạng cơ bản với cấu trúc 4 lớp nguyên tử dạng $X-M-M-X$. Nhằm triệt tiêu hoàn toàn tương tác tĩnh điện giả mạo giữa các lớp lân cận sinh ra do điều kiện biên tuần hoàn 3D, một khoảng chân không (vacuum space) dày tới $20\text{ \AA}$ được thiết lập dọc theo trục vuông góc với mặt phẳng ($z$).
  2. Chọn lựa hàm cơ sở và hàm thế: Sử dụng phương pháp sóng phẳng giả thế (pseudopotentials) chuẩn hóa với mức năng lượng cắt (kinetic energy cutoff) cho hàm sóng sóng phẳng là $E_{cut} = 500\text{ eV}$. Phép tích phân trên vùng Brillouin được lấy mẫu thông qua lưới $k$-point dịch chuyển Monkhorst-Pack kích thước $(15 \times 15 \times 1)$.
  3. Hiệu chỉnh tương tác van der Waals và lai hóa: Để mô tả chính xác tương tác phân tán tầm xa, phương pháp hiệu chỉnh DFT-D2 của Grimme được tích hợp trong phương trình Kohn-Sham. Để khắc phục hiện tượng đánh giá thấp độ rộng vùng cấm cố hữu của phiếm hàm GGA-PBE do sai số tự tương tác, phiếm hàm lai HSE06 với tham số trộn trao đổi Hartree-Fock $a = 0.25$ và tham số sàng lọc phạm vi ngắn $\omega = 0.15\text{ \AA}^{-1}$ được áp dụng đồng thời.
  4. Quy trình kiểm tra độ tin cậy:
    • Độ bền nhiệt: Chạy mô phỏng Động lực học phân tử Ab-initio (AIMD) sử dụng bộ điều nhiệt Nosé-Hoover ở nhiệt độ phòng ($300\text{ K}$) với bước thời gian tích phân $1.0\text{ fs}$ trong tổng thời gian khảo sát tối thiểu $5\text{ ps}$ trên các siêu ô mạng $(4 \times 4 \times 1)$ và $(6 \times 6 \times 1)$.
    • Độ bền động học: Tính toán phổ phonon trên toàn bộ vùng Brillouin thông qua siêu ô mạng $(6 \times 6 \times 1)$ gồm 144 nguyên tử (72 nguyên tử $M$ và 72 nguyên tử $X$) bằng lý thuyết DFPT. Việc không xuất hiện các mode dao động có tần số ảo (imaginary frequencies) tại điểm $\Gamma$ khẳng định cấu trúc đạt cực tiểu năng lượng tự do.
    • Độ bền cơ học: Tính toán các hệ số đàn hồi độc lập $C_{11}, C_{12}, C_{22}, C_{66}$ và kiểm chứng theo tiêu chuẩn Born-Huang cho hệ tinh thể phẳng hai chiều: $C_{11} > 0$ và $C_{11}C_{22} - C_{12}^2 > 0$.

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

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

Luận án đã công bố hệ thống kết quả định lượng chuẩn xác, làm sáng tỏ các tính chất vật lý chưa từng được báo cáo đồng bộ trước đây:

                            PHỔ BIẾN ĐIỆU VÙNG CẤM DƯỚI TRƯỜNG LỰC NGOÀI
   Eg (eV)
     ▲
     │       /────── Bán dẫn vùng cấm gián tiếp (Indirect)
     │      /
     │     /─────── Bán dẫn vùng cấm trực tiếp (Direct)
     │    /
     │   /───────── Kim loại hóa (Semiconductor-to-Metal Transition, Eg = 0)
     │  /
     └────────────────────────────────────────────────────────► Biến dạng εxy / Điện trường E
  1. Tiến hóa tham số mạng và năng lượng liên kết: Hằng số mạng cân bằng ($a$) của các đơn lớp $MX$ đối xứng tăng đơn điệu theo kích thước bán kính nguyên tử chalcogenide: từ $3.59\text{ \AA}$ ($\text{GaS}$) lên $4.38\text{ \AA}$ ($\text{InTe}$). Năng lượng liên kết kết tụ ($E_{coh}$) luôn mang giá trị âm lớn, khẳng định quá trình hình thành đơn lớp tỏa nhiệt mạnh và có tính khả thi cao trong thực nghiệm tổng hợp hóa học.
  2. Sự dịch chuyển loại vùng cấm: Hầu hết các đơn lớp $MX$ nguyên bản là chất bán dẫn vùng cấm gián tiếp với cực đại vùng hóa trị (VBM) nằm trên đoạn $\Gamma-K$ và cực tiểu vùng dẫn (CBM) định vị tại điểm $\Gamma$. Tuy nhiên, khi chuyển sang cấu trúc Janus tam nguyên $M_2XY$ (như $\text{Ga}_2\text{SSe}, \text{Ga}_2\text{STe}, \text{In}_2\text{SSe}$), sự tái phân bố mật độ điện tích đã chuyển đổi một số cấu trúc sang bán dẫn vùng cấm trực tiếp tại điểm $\Gamma$, với độ rộng vùng cấm HSE06 trải dài từ $1.15\text{ eV}$ đến $2.85\text{ eV}$.
  3. Hiện tượng phân cực bề mặt và thế tĩnh điện: Trong các cấu trúc Janus $\text{GaIn}XO$, sự chênh lệch độ âm điện giữa Oxi ($3.44$) và các nguyên tố chalcogenide ($S: 2.58, Se: 2.55, Te: 2.10$) đã tạo ra một gradient thế tĩnh điện khổng lồ. Công thoát điện tử tại mặt O ($\Phi_O$) lớn hơn đáng kể so với mặt $X$ ($\Phi_X$), tạo nên bước nhảy thế $\Delta \Phi = \Phi_X - \Phi_O$ vượt mức $1.5\text{ eV}$.
  4. Độ linh động hạt tải siêu cao và tính dị hướng: Đối với các đơn lớp Janus $\text{GaIn}XO$ (như $\text{SGaInO}$ và $\text{SeGaInO}$), khối lượng hiệu dụng của electron ($m_e^*$) ghi nhận các giá trị rất nhỏ ($\sim 0.12 - 0.25\text{ }m_0$), kết hợp cùng hằng số thế biến dạng thấp ($E_d < 3.0\text{ eV}$) dọc theo trục chuyển vị, đã đem lại độ linh động electron lý thuyết đạt mức $2.5 \times 10^3\text{ cm}^2\text{V}^{-1}\text{s}^{-1}$ đến $8.2 \times 10^3\text{ cm}^2\text{V}^{-1}\text{s}^{-1}$ ở nhiệt độ phòng.
  5. Hiệu ứng chuyển pha bán dẫn - kim loại dưới biến dạng cơ học: Khi áp dụng biến dạng kéo/nén hai trục đối xứng ($\varepsilon_{xy}$) trong phạm vi $\pm 10%$, cấu trúc vùng năng lượng của $MX$ và $M_2XY$ trải qua sự biến điệu mạnh mẽ. Biến dạng kéo làm dịch chuyển mức VBM tiến về mức Fermi ($E_F$), kích hoạt sự đóng hoàn toàn vùng cấm ($E_g = 0\text{ eV}$) ở mức biến dạng tới hạn $\varepsilon_{cr} \approx +8%$, đánh dấu bước chuyển pha bán dẫn – kim loại có thể đảo ngược.

Implications đa chiều

  • Về mặt lý thuyết: Cung cấp mô hình toán – lý chuẩn mực mô tả sự biến điệu cấu trúc dải năng lượng lượng tử của các mạng tinh thể phân lớp chịu ứng suất cơ nhiệt và điện trường ngoài, bổ sung mắt xích quan trọng vào lý thuyết chất rắn phi đối xứng.
  • Về mặt phương pháp: Chuẩn hóa quy trình tính toán phối hợp giữa DFPT, AIMD và lý thuyết thế biến dạng 2D, có thể chuyển giao áp dụng trực tiếp cho các họ vật liệu 2D mới nổi khác như MXenes, Pnictogens hay Perovskites đơn lớp.
  • Ứng dụng thực tiễn trong công nghệ vi điện tử: Các thông số độ linh động cao và công thoát biến tính là dữ liệu gốc quý giá để thiết kế các transistor hiệu ứng trường kênh dẫn 2D thế hệ mới, linh kiện logic tiêu thụ năng lượng cực thấp (sub-10 nm gate FETs).
  • Ứng dụng trong quang điện tử và năng lượng sạch: Độ rộng vùng cấm nằm trong vùng quang phổ khả kiến – tử ngoại cùng với sự phân tách thế bề mặt nội tại biến các đơn lớp Janus $M_2XY$ và $\text{GaIn}XO$ thành vật liệu xúc tác quang lý tưởng cho quá trình tách hydro từ nước ($\text{H}_2\text{O} \to \text{H}_2 + \frac{1}{2}\text{O}_2$) mà không cần chất trợ xúc tác ngoài.

Limitations và Future Research

Mặc dù đạt được những bước tiến lý thuyết đột phá, công trình vẫn ghi nhận một số giới hạn học thuật cần được mở rộng trong tương lai:

  1. Giới hạn xấp xỉ nhiệt độ không tuyệt đối ($0\text{ K}$): Hầu hết các tính toán cấu trúc điện tử chuẩn của DFT được thực thi ở điều kiện tĩnh $T = 0\text{ K}$. Mặc dù độ bền nhiệt đã được kiểm chứng qua AIMD ở $300\text{ K}$, sự dịch chuyển dải năng lượng do giãn nở nhiệt và tán xạ electron-phonon phụ thuộc nhiệt độ chưa được lượng hóa đầy đủ qua phương pháp phiếm hàm Allen-Heine-Cardona.
  2. Bỏ qua hiệu ứng kích thích đa hạt (Many-Body Excitonic Effects): Các tính toán dải năng lượng dừng lại ở mức phiếm hàm lai HSE06. Để thu được phổ hấp thụ quang học có độ chính xác tuyệt đối, cần triển khai phương pháp Many-Body Perturbation Theory kết hợp gần đúng $GW$ và phương trình Bethe-Salpeter (GW-BSE) nhằm tính đến năng lượng liên kết exciton.
  3. Mô hình hóa khuyết tật và ảnh hưởng môi trường: Nghiên cứu tập trung vào các phiến tinh thể lý tưởng tuyệt đối. Trong thực tế chế tạo, sự xuất hiện của các sai hỏng mạng (vacancy defects), ranh giới hạt (grain boundaries) và tương tác điện môi từ chất nền ($\text{SiO}_2$, $\text{h-BN}$) sẽ tạo ra các mức bẫy điện tử cục bộ.

Định hướng nghiên cứu tiếp nối:

  • Mở rộng tính toán hiệu ứng kích thích đa hạt $GW_0\text{-BSE}$ để xác định chính xác năng lượng liên kết exciton và phổ điện quang phân cực.
  • Nghiên cứu mô hình hóa các dị cấu trúc van der Waals đa tầng ghép nối giữa Janus $M_2XY$ và các vật liệu 2D khác (như Graphene, $\text{MoS}_2$, $\text{WS}_2$) nhằm tối ưu hóa sự phân tách hạt tải quang điện.
  • Khảo sát tương tác giao diện và ảnh hưởng của hiệu ứng chất nền điện môi cao ($\text{HfO}_2$, $\text{Al}_2\text{O}_3$) lên độ linh động thực tế của linh kiện FET.

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

  • Tác động học thuật sâu rộng: Kết quả của luận án đã được công bố trên 08 công trình khoa học uy tín, bao gồm 06 bài báo chuyên ngành trên các tạp chí quốc tế thuộc danh mục ISI/SCIE danh giá và 02 bài báo trên tạp chí chuyên ngành trong nước. Tỷ lệ công bố quốc tế dày đặc khẳng định chất lượng học thuật đạt chuẩn mực toàn cầu, đóng góp nguồn trích dẫn chuyên môn cao cho cộng đồng vật lý tính toán.
  • Đột phá công nghệ vật liệu: Mở ra hướng đi mới cho các phòng thí nghiệm thực nghiệm trong và ngoài nước trong việc định hướng tổng hợp hóa học màng mỏng bán dẫn đơn lớp bằng phương pháp bóc tách pha lỏng, lắng đọng pha hơi hóa học (CVD) hoặc epitaxy chùm phân tử (MBE).
  • Đóng góp phát triển nguồn nhân lực chất lượng cao: Đánh dấu sự làm chủ toàn diện các kỹ thuật mô phỏng cơ học lượng tử phức tạp cao trên hệ thống máy tính hiệu năng cao (HPC) tại Việt Nam, khẳng định năng lực hội nhập quốc tế của nền vật lý lý thuyết nước nhà.

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

                                  MA TRẬN HƯỞNG LỢI HỌC THUẬT & KỸ NGUYỆN
┌──────────────────────────────────────────────┬──────────────────────────────────────────────┐
│ CỘNG ĐỒNG HỌC THUẬT & NGHIÊN CỨU SINH        │ KỸ SƯ R&D CÔNG NGHỆ BÁN DẪN & VI ĐIỆN TỬ     │
│ - Tiếp cận bộ thông số hằng số mạng, phonon, │ - Ứng dụng số liệu khối lượng hiệu dụng m*,  │
│   cấu trúc dải năng lượng chuẩn xác tuyệt đối│   độ linh động μ để thiết kế kênh dẫn FET.   │
│ - Tái sử dụng quy trình tối ưu hóa DFT-D2,   │ - Tận dụng gradient thế tĩnh điện Janus      │
│   HSE06 và DFPT trên Quantum ESPRESSO.       │   trong thiết kế cảm biến và tế bào quang điện│
├──────────────────────────────────────────────┼──────────────────────────────────────────────┤
│ NHÀ HOẠCH ĐỊNH CHÍNH SÁCH KHOA HỌC           │ CHUYÊN GIA NĂNG LƯỢNG TÁI TẠO & QUANG HÓA    │
│ - Căn cứ khoa học định hướng đầu tư trọng    │ - Dữ liệu vị trí dải năng lượng VBM/CBM      │
│   điểm vào công nghệ bán dẫn và vật liệu mới │   so với mức thế oxy hóa khử của nước phục   │
│   theo các chương trình công nghệ quốc gia.  │   vụ chế tạo vật liệu xúc tác tách Hydro xanh│
└──────────────────────────────────────────────┴──────────────────────────────────────────────┘

Câu hỏi chuyên sâu

1. Đóng góp lý thuyết độc đáo và đột phá nhất của luận án là gì?

Đóng góp mang tính cách mạng nhất là việc thiết lập mô hình vi mô toàn diện về sự phá vỡ đối xứng gương hai chiều trong các đơn lớp Janus tam nguyên $M_2XY$ và tứ nguyên $\text{GaIn}XO$. Bằng cách kết hợp phiếm hàm lai HSE06 với hiệu ứng spin-quỹ đạo (SOC), luận án đã chứng minh rằng sự bất đối xứng ở hai bề mặt không chỉ tạo ra một mômen lưỡng cực tĩnh điện nội tại không đổi $\Delta \Phi$, mà còn giải phóng sự suy biến dải năng lượng, trực tiếp biến đổi vật liệu từ bán dẫn vùng cấm gián tiếp thành bán dẫn vùng cấm trực tiếp với độ linh động hạt tải vượt trội.

2. Sự đổi mới về mặt phương pháp luận tính toán so với các công trình quốc tế cùng thời kỳ?

So với các nghiên cứu trước đây thường chỉ dựa trên phép gần đúng chuẩn GGA-PBE (vốn đánh giá thấp độ rộng vùng cấm từ $30 - 50%$) hoặc bỏ qua tương tác van der Waals tầm xa, luận án đã xây dựng một giao thức mô phỏng tích hợp 5 tầng:

  1. Tối ưu hóa cấu trúc có gắn hiệu chỉnh tán xạ Grimme DFT-D2;
  2. Sàng lọc dao động mạng qua siêu ô mạng $6 \times 6 \times 1$ với lý thuyết DFPT;
  3. Khảo sát động học phân tử AIMD $300\text{ K}$ kiểm chứng độ bền nhiệt;
  4. Tính toán cấu trúc dải chính xác cao bằng phiếm hàm lai HSE06 kết hợp ma trận SOC;
  5. Trích xuất hằng số thế biến dạng $E_d$ và độ linh động $\mu$ theo lý thuyết Bardeen-Shockley.

3. Phát hiện thực nghiệm/tính toán nào bất ngờ và phản trực giác nhất?

Phát hiện bất ngờ nhất là hiện tượng đảo chiều đáp ứng dải năng lượng dưới tác động của biến dạng cơ học: trái ngược với các vật liệu bán dẫn truyền thống vốn có vùng cấm thu hẹp khi bị nén mạng tinh thể, ở một số đơn lớp monochalcogenide chứa Indium ($\text{InS}, \text{InSe}$), biến dạng nén hai trục lại làm gia tăng độ rộng vùng cấm cục bộ do sự tái sắp xếp định hướng của các liên kết cộng hóa trị $\text{In}-\text{In}$ vuông góc với mặt phẳng. Ngoài ra, việc áp dụng biến dạng kéo $\varepsilon_{xy} \approx +8%$ đã gây ra hiện tượng kim loại hóa hoàn toàn mà không làm sụp đổ độ bền động học của mạng tinh thể.

4. Giao thức tái lập nghiên cứu (Replication Protocol) có được đảm bảo không?

Hoàn toàn minh bạch và chuẩn hóa. Luận án công bố chi tiết toàn bộ các tham số kỹ thuật đầu vào: phần mềm tính toán (Quantum ESPRESSO), chuẩn giả thế, mức năng lượng cắt sóng phẳng ($500\text{ eV}$), lưới tích phân Monkhorst-Pack ($15 \times 15 \times 1$), chiều dày lớp chân không ($20\text{ \AA}$), ngưỡng hội tụ năng lượng ($10^{-6}\text{ eV}$) và ngưỡng lực tối ưu ($10^{-3}\text{ eV/\AA}$). Mọi nhà nghiên cứu độc lập đều có thể tái lập chính xác các kết quả cấu trúc dải và phổ phonon này.

5. Lộ trình phát triển học thuật 10 năm được phác thảo ra sao?

Luận án mở ra chương trình nghiên cứu dài hạn trong thập kỷ tới:

  • Giai đoạn 2022–2025: Phát triển các thuật toán học máy (Machine Learning Interatomic Potentials) dựa trên dữ liệu DFT của luận án để sàng lọc nhanh hàng nghìn dị cấu trúc van der Waals 2D.
  • Giai đoạn 2026–2029: Phối hợp cùng các nhóm thực nghiệm chế tạo linh kiện FET và cảm biến quang dựa trên màng mỏng $\text{GaSe}$ và Janus $\text{Ga}_2\text{SSe}$.
  • Giai đoạn 2030–2035: Tích hợp vật liệu monochalcogenide vào cấu trúc chip tích hợp quang - điện tử 3D nguyên khối (Monolithic 3D Integration), thương mại hóa các linh kiện tách hydro năng lượng mặt trời.

Kết luận

Luận án tiến sĩ của tác giả Vũ Thị Tuyết Vi là một công trình nghiên cứu khoa học công phu, mẫu mực và mang tính đột phá cao trong lĩnh vực Vật lý chất rắn tính toán. Công trình đã đúc kết các giá trị học thuật cốt lõi:

  1. Xác lập bức tranh toàn cảnh và chính xác tuyệt đối về đặc tính cấu trúc, độ bền cơ học, độ bền nhiệt động và phổ phonon của 6 đơn lớp monochalcogenide đối xứng nhóm III ($MX$).
  2. Tiên phong giải mã các tính chất điện tử mới nổi của hệ vật liệu bất đối xứng gương Janus tam nguyên $M_2XY$ và dị hợp bậc bốn $\text{GaIn}XO$, khẳng định sự tồn tại của trường điện nội tại mạnh và bước nhảy thế tĩnh điện $\Delta \Phi$.
  3. Khám phá quy luật biến điệu dải năng lượng lượng tử dưới tác động hiệp đồng của biến dạng cơ học đa trục và điện trường ngoài, làm sáng tỏ cơ chế chuyển pha bán dẫn – kim loại có thể điều khiển được.
  4. Lượng hóa độ linh động hạt tải siêu cao của các cấu trúc Janus, mở ra tiềm năng to lớn cho việc chế tạo linh kiện vi điện tử kênh dẫn nano tốc độ cao.
  5. Mở ra 3 nhánh nghiên cứu thực nghiệm và lý thuyết mới: kỹ nghệ biến dạng 2D (Strain-tronics), điện tử học spin dựa trên phân cực Janus (Janus Spintronics), và quang xúc tác tách nước thế hệ mới.
  6. Minh chứng chất lượng nghiên cứu đỉnh cao thông qua 08 công trình khoa học đã xuất bản, trong đó có 06 bài báo trên các tạp chí quốc tế thuộc danh mục SCIE uy tín, tạo dấu ấn học thuật mạnh mẽ và khẳng định vị thế khoa học của nghiên cứu vật lý Việt Nam trên trường quốc tế.