Tổng quan luận án

Luận án tiến sĩ "Nghiên cứu tính chất điện tử của vật liệu hai chiều monochalcogenide" được tác giả Võ Thị Tuyết Vi thực hiện tại Trường Đại học Sư phạm – Đại học Huế, thuộc chuyên ngành Vật lý lý thuyết và vật lý toán (mã số: 9 44 01 03), dưới sự hướng dẫn khoa học của PGS. Nguyễn Ngọc Hiếu và PGS. Bùi Đình Hợi.

Tính cấp thiết và khoảng trống nghiên cứu

Kể từ khi graphene được cô lập thành công từ graphite vào năm 2004, các vật liệu hai chiều (2D) có cấu trúc lớp đã thu hút sự quan tâm lớn nhờ các đặc tính vật lý đặc trưng. Tuy nhiên, graphene ở dạng nguyên thủy là một bán kim loại có độ rộng vùng cấm bằng không ($E_g = 0\text{ eV}$) và tương tác spin-quỹ đạo (SOC) rất yếu, dẫn đến hạn chế trong việc ứng dụng vào các linh kiện bán dẫn như transistor hiệu ứng trường (FET) hay các thiết bị điện tử spin (spintronics). Dù có thể mở vùng cấm cho graphene bằng biến dạng cơ học, pha tạp hay tạo dị cấu trúc van der Waals, các tác động ngoại cảnh này thường gây biến đổi trường nội tại và hiệu ứng kích thước ngoài ý muốn.

Để khắc phục những hạn chế đó, nhiều hệ vật liệu hai chiều khác đã được nghiên cứu và tổng hợp như silicene, germanene, stanene, các hợp chất kim loại chuyển tiếp dichalcogenide (TMD) và đặc biệt là họ vật liệu monochalcogenide. Vật liệu hai chiều monochalcogenide nhóm III (dạng $MX$, với $M = \text{Ga, In}$ và $X = \text{S, Se, Te}$) sở hữu vùng cấm tự nhiên phù hợp, độ linh động hạt tải cao và khả năng hấp thụ quang học tốt. Mặc dù một số đơn lớp như GaS, GaSe, InS, InSe đã được chế tạo thực nghiệm, các hiểu biết mang tính hệ thống về tính chất điện tử, phản ứng dưới tác động của biến dạng cơ học và điện trường ngoài vẫn còn nhiều vấn đề bỏ ngỏ. Đặc biệt, việc phá vỡ tính đối xứng gương theo phương thẳng đứng để hình thành cấu trúc bất đối xứng Janus (như $M_2XY$ hoặc các biến thể oxy hóa $\text{GaIn}XO$) tạo ra điện trường nội tại định hướng, mở ra các tính chất vật lý chưa từng xuất hiện ở vật liệu đối xứng nguyên thủy nhưng chưa được khảo sát đầy đủ.

Mục tiêu và nhiệm vụ nghiên cứu

Mục tiêu tổng quát của luận án là sử dụng lý thuyết phiếm hàm mật độ (DFT) để nghiên cứu một cách hệ thống các đặc trưng cấu trúc và tính chất điện tử của đơn lớp monochalcogenide nhóm III ($MX$) cùng các cấu trúc bất đối xứng Janus hình thành từ nhóm này ($M_2XY$ và $\text{GaIn}XO$).

Các nhiệm vụ nghiên cứu cụ thể được đánh số theo văn bản luận án gồm:

  1. Xây dựng mô hình và nghiên cứu các tính chất về cấu trúc nguyên tử của các vật liệu hai chiều đơn lớp monochalcogenide nhóm III $MX$ và các vật liệu Janus $M_2XY$, $\text{GaIn}XO$ ($M = \text{Ga, In}$; $X/Y = \text{S, Se, Te}$; $X \neq Y$).
  2. Tính toán và đánh giá các tính chất bền vững, tính chất biến dạng và điều kiện để tồn tại trong thực tế của các cấu trúc hai chiều đơn lớp, đặc biệt là các đơn lớp Janus $M_2XY$ và $\text{GaIn}XO$.
  3. Nghiên cứu các tính chất điện tử của các đơn lớp monochalcogenide nhóm III $MX$ và Janus monochalcogenide nhóm III $M_2XY$, $\text{GaIn}XO$.
  4. Nghiên cứu độ linh động của điện tử trong các đơn lớp Janus $\text{GaIn}XO$.
  5. Nghiên cứu ảnh hưởng của các điều kiện bên ngoài như biến dạng cơ học hay điện trường ngoài lên các tính chất điện tử của các đơn lớp monochalcogenide nhóm III và Janus monochalcogenide nhóm III.

Đối tượng và phạm vi nghiên cứu

  • Đối tượng nghiên cứu: Các vật liệu hai chiều đơn lớp monochalcogenide nhóm III đối xứng $MX$ ($M = \text{Ga, In}$; $X = \text{S, Se, Te}$); các đơn lớp bất đối xứng Janus $M_2XY$ ($M = \text{Ga, In}$; $X/Y = \text{S, Se, Te}$; $X \neq Y$); và 6 cấu hình cấu trúc của đơn lớp Janus oxy hóa $\text{GaIn}XO$ (gồm 3 cấu trúc $X\text{GaInO}$ và 3 cấu trúc $\text{OGaIn}X$).
  • Phạm vi không gian: Khảo sát các cấu trúc tinh thể ở quy mô đơn lớp hai chiều, mô phỏng trong điều kiện tuần hoàn 2D có lớp chân không cách ly theo phương thẳng đứng ($z$).
  • Phạm vi tham số vật lý: Khảo sát trạng thái cân bằng cơ bản, tác động của biến dạng cơ học hai trục ($\varepsilon_{xy}$ từ $-10%$ đến $+10%$), biến dạng nhỏ tính hệ số đàn hồi (từ $-0{,}015$ đến $+0{,}015$), điện trường ngoài phân cực vuông góc và độ ổn định nhiệt động lực học ở nhiệt độ phòng ($300\text{ K}$).

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

Các hướng nghiên cứu trước được luận án điểm lại

Luận án đã tổng quan lịch sử phát triển của các hệ vật liệu hai chiều:

  • Graphene và các vật liệu đơn nguyên tử nhóm IV: Wallace (1947) công bố tính chất điện tử lý thuyết của graphene; Novoselov và Geim (2004) bóc tách thực nghiệm thành công đơn lớp graphene. Các nghiên cứu về silicene (Cahangirov và cộng sự, 2009; Vogt và cộng sự, 2012), germanene (Bianco và cộng sự, 2013) và stanene (Zhu và cộng sự, 2015) chỉ ra cấu trúc lục giác so le với độ dày hữu hạn và tương tác SOC mở ra vùng cấm nhỏ ($1{,}55\text{ meV}$ ở silicene và $23{,}9\text{ meV}$ ở germanene).
  • Hợp chất kim loại chuyển tiếp dichalcogenide (TMD): Các nghiên cứu về $\text{MoS}_2$, $\text{WSe}_2$ ở các pha $2H, 1T, 3R$ (Mak và cộng sự, 2010; Radisavljevic và cộng sự, 2011) chứng minh sự chuyển đổi từ vùng cấm gián tiếp ở dạng khối sang vùng cấm trực tiếp ở đơn lớp, cùng ứng dụng chế tạo transistor có độ linh động đạt $200\text{ cm}^2/\text{Vs}$.
  • Monochalcogenide nhóm IV và nhóm III: Các nghiên cứu về $\text{GeS}, \text{GeSe}, \text{SnS}, \text{SnSe}$ (Fei và cộng sự, 2015; Gomes và cộng sự, 2015) xác định cấu trúc trực thoi và hệ số áp điện lớn. Đối với monochalcogenide nhóm III, các công trình thực nghiệm đã tổng hợp thành công các tấm nano GaS (Late và cộng sự, 2012), GaSe (Lei và cộng sự, 2013), InS (Feng và cộng sự, 2015), InSe (Bandurin và cộng sự, 2017) và GaTe (Liu và cộng sự, 2014). Về mặt lý thuyết, Yagmurcukardes và cộng sự (2019) khảo sát biến dạng kéo trên GaS và GaSe, dự đoán chuyển pha bán dẫn - kim loại ở biến dạng $17%$ (GaS) và $14%$ (GaSe); Zhou và cộng sự khảo sát điều khiển vùng cấm qua hấp phụ khí trên GaSe.
  • Cấu trúc bất đối xứng Janus 2D: Thuật ngữ hạt Janus bắt nguồn từ nghiên cứu của Casagrande và cộng sự (1988). Đến năm 2017, hai nhóm nghiên cứu độc lập của Lu và Zhang đã tổng hợp thành công đơn lớp Janus MoSSe. Tiếp theo, Wang và cộng sự nghiên cứu Janus WSSe và WSeTe; Guo và cộng sự (2017) thiết kế các đơn lớp Janus monochalcogenide nhóm III và chỉ ra hệ số áp điện đạt $8{,}47\text{ pm/V}$; Huang và cộng sự dự đoán ứng dụng quang xúc tác phân tách nước của nhóm vật liệu này.

Khoảng trống luận án lựa chọn giải quyết

Mặc dù các đơn lớp đối xứng $MX$ và một số đơn lớp Janus TMD đã được nghiên cứu, việc khảo sát toàn diện và có hệ thống cấu trúc vùng năng lượng, thế tĩnh điện, công thoát, phổ dao động phonon, độ bền nhiệt và đặc biệt là độ linh động hạt tải theo lý thuyết thế biến dạng trên các đơn lớp bất đối xứng Janus $M_2XY$ và các biến thể oxy hóa $\text{GaIn}XO$ dưới ảnh hưởng của biến dạng cơ học và điện trường ngoài vẫn chưa được hoàn thiện. Luận án giải quyết khoảng trống này bằng các tính toán ab-initio định lượng chi tiết.

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

Khung lý thuyết và mô hình toán học

Luận án sử dụng Lý thuyết phiếm hàm mật độ (DFT) làm nền tảng tính toán:

  • Các định lý Hohenberg-Kohn (1964):
    1. Định lý thứ nhất: Thế ngoài $V_{ext}(\vec{r})$ được xác định duy nhất bởi mật độ điện tử ở trạng thái cơ bản $\rho_0(\vec{r})$, từ đó xác định toán tử Hamilton và toàn bộ các tính chất của hệ.
    2. Định lý thứ hai: Nguyên lý biến phân xác lập rằng với mọi mật độ điện tử thử nghiệm $\rho(\vec{r}) \ge 0$ thỏa mãn chuẩn hóa $\int \rho(\vec{r})d\vec{r} = N$, năng lượng phiếm hàm $E[\rho] \ge E_0$, dấu đẳng thức xảy ra khi và chỉ khi $\rho \equiv \rho_0$.
  • Hệ phương trình Kohn-Sham (1965): Chuyển bài toán tương tác nhiều hạt phức tạp về hệ các phương trình một hạt phi tương tác chuyển động trong thế hiệu dụng $V_{eff}(\vec{r})$: $$\left[ -\frac{1}{2}\nabla^2 + V_H(\vec{r}) + V_{ext}(\vec{r}) + V_{XC}(\vec{r}) \right] \phi_i(\vec{r}) = \varepsilon_i \phi_i(\vec{r})$$
  • Phiếm hàm trao đổi - tương quan ($E_{XC}$):
    • Gần đúng gradient suy rộng (GGA) theo dạng tham số hóa Perdew-Burke-Ernzerhof (PBE).
    • Phiếm hàm lai Heyd-Scuseria-Ernzerhof (HSE06) với hệ số trộn trao đổi Hartree-Fock $a = 1/4$ và tham số chia miền sàng lọc $\omega$ để xác định chính xác độ rộng vùng cấm.
  • Các hiệu chỉnh mở rộng:
    • Hiệu chỉnh tương tác phân tán van der Waals tầm xa theo phương pháp bán thực nghiệm DFT-D2 của Grimme.
    • Tương tác spin-quỹ đạo (SOC) qua phương pháp PBE+SOC để khảo sát sự phân tách mức năng lượng.
    • Nhúng thành phần điện trường ngoài $\vec{E}$ vào toán tử Hamilton để giải phương trình trạng thái dưới điện trường vuông góc.

Phương pháp phân tích động học, cơ học và truyền dẫn

  • Độ ổn định động học và nhiệt: Tính toán phổ phonon bằng Lý thuyết nhiễu loạn phiếm hàm mật độ (DFPT). Khảo sát độ bền nhiệt qua mô phỏng động lực học phân tử ab-initio (AIMD) sử dụng thuật toán điều nhiệt Nosé-thermostat ở nhiệt độ phòng ($300\text{ K}$).
  • Tính chất đàn hồi: Tính tensor hệ số đàn hồi $C_{ij}$, mô-đun Young $Y_{2D}$ và tỷ số Poisson $\nu$ từ sự thay đổi năng lượng toàn phần dưới các mức biến dạng nhỏ $\pm 0{,}015$, đối chiếu với tiêu chuẩn ổn định cơ học Born-Huang ($C_{11}C_{22} - C_{12}^2 > 0$ và $C_{66} > 0$).
  • Độ linh động hạt tải: Tính toán theo mô hình lý thuyết thế biến dạng (Bardeen-Shockley) cho hệ hai chiều: $$\mu = \frac{e\hbar^3 C_{2D}}{k_B T m^* m_d (E_d)^2}$$ trong đó $m^*$ là khối lượng hiệu dụng, $C_{2D}$ là mô-đun đàn hồi 2D, $E_d$ là hằng số thế biến dạng dọc theo trục truyền dẫn.

Công cụ và thông số kỹ thuật mô phỏng

Tất cả các tính toán được thực hiện trên phần mềm mã nguồn mở Quantum ESPRESSO:

  • Năng lượng cắt sóng phẳng (cutoff energy): $500\text{ eV}$.
  • Lưới không gian xung lượng $k$-mesh: $(15 \times 15 \times 1)$ theo sơ đồ phân bố Monkhorst-Pack.
  • Ngưỡng hội tụ năng lượng tự hợp: $10^{-6}\text{ eV}$; ngưỡng hội tụ lực tác dụng lên nguyên tử: $10^{-3}\text{ eV/\AA}$.
  • Khoảng chân không cách ly giữa các lớp định kỳ dọc trục $z$: $20\text{ \AA}$.
  • Siêu ô tính phổ phonon DFPT: $(6 \times 6 \times 1)$ gồm 144 nguyên tử.

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

Chương 1: Tổng quan về vật liệu hai chiều có cấu trúc lớp và lý thuyết phiếm hàm mật độ

Chương này hệ thống hóa các kiến thức cơ sở về vật liệu hai chiều và nền tảng phương pháp DFT:

  • Tổng quan các hệ vật liệu hai chiều tiêu biểu: cấu trúc tinh thể mạng tổ ong của graphene, cơ chế tán xạ Dirac tại điểm $K/K'$, các hệ đơn nguyên tử nhóm IV (silicene, germanene, stanene) và các pha tinh thể của hợp chất TMD ($2H, 1T, 3R$).
  • Phân tích họ vật liệu monochalcogenide nhóm III và nhóm IV, nhấn mạnh các ưu điểm về độ rộng vùng cấm hữu hạn và độ nhạy quang điện.
  • Giới thiệu khái niệm cấu trúc bất đối xứng Janus hai chiều, cơ chế xuất hiện mô-men lưỡng cực và điện trường nội tại định hướng vuông góc với mặt phẳng do sự chênh lệch độ âm điện giữa hai lớp nguyên tử ngoài cùng.
  • Trình bày chi tiết toán học của hai định lý Hohenberg-Kohn, phương trình một hạt Kohn-Sham, các phép gần đúng phiếm hàm năng lượng trao đổi tương quan (LDA, GGA-PW91, GGA-PBE, HSE06) và thuật toán vòng lặp tự hợp trường (SCF) triển khai trong Quantum ESPRESSO.

Chương 2: Tính chất điện tử của vật liệu hai chiều monochalcogenide nhóm III

Chương 2 tập trung khảo sát chi tiết 6 đơn lớp monochalcogenide nhóm III đối xứng $MX$ ($M = \text{Ga, In}$; $X = \text{S, Se, Te}$).

Cấu trúc tinh thể đơn lớp MX (pha 2H, nhóm đối xứng D3h):
   (Lớp 1: X)      X       X       X
                   |       |       |
   (Lớp 2: M)      M ————— M ————— M
                   |       |       |
   (Lớp 3: M)      M ————— M ————— M
                   |       |       |
   (Lớp 4: X)      X       X       X
   Thứ tự xếp lớp: X — M — M — X
  • Cấu trúc nguyên tử và độ ổn định: Cấu trúc đơn lớp $MX$ thuộc nhóm đối xứng không gian $D_{3h}$, ô mạng cơ sở gồm 4 nguyên tử xếp theo thứ tự $X\text{--}M\text{--}M\text{--}X$. Hằng số mạng $a$ tăng dần khi bán kính nguyên tử chalcogen tăng (từ GaS đến InTe). Phổ phonon DFPT của cả 6 vật liệu không xuất hiện các nhánh tần số ảo âm, chứng minh độ ổn định động học. Các hệ số đàn hồi thỏa mãn tiêu chuẩn Born-Huang, khẳng định độ bền cơ học.
  • Tính chất điện tử và vùng năng lượng: Ở trạng thái cân bằng, tất cả 6 đơn lớp $MX$ đều là chất bán dẫn có vùng cấm gián tiếp (xiên). Cực đại vùng hóa trị (VBM) nằm tại hoặc gần điểm $\Gamma$, trong khi cực tiểu vùng dẫn (CBM) nằm tại điểm $M$ hoặc trên đoạn $\Gamma\text{--}K$. Mức độ ảnh hưởng của tương tác spin-quỹ đạo (SOC) tăng rõ rệt ở các hợp chất chứa nguyên tử nặng như In và Te.
Vật liệu Hằng số mạng $a$ ($\text{\AA}$) Độ dài liên kết $d_{MX}$ ($\text{\AA}$) Độ dài liên kết $d_{MM}$ ($\text{\AA}$) Độ dày đơn lớp $\Delta h$ ($\text{\AA}$) Mức Fermi $E_F$ ($\text{eV}$) Công thoát $\Phi$ ($\text{eV}$) $E_g^{\text{PBE}}$ ($\text{eV}$) $E_g^{\text{PBE+SOC}}$ ($\text{eV}$)
GaS 3,59 2,34 2,45 4,62 -4,64 5,38 2,42 2,41
GaSe 3,75 2,47 2,47 4,87 -4,38 5,19 2,01 1,93
GaTe 4,06 2,67 2,50 5,34 -3,98 4,80 1,51 1,35
InS 3,94 2,53 2,80 5,15 -4,69 5,42 1,88 1,84
InSe 4,00 2,65 2,82 5,38 -4,42 5,23 1,48 1,40
InTe 4,38 2,84 2,85 5,86 -4,02 4,82 1,02 0,85
  • Ảnh hưởng của biến dạng cơ học ($\varepsilon_{xy}$): Biến dạng nén làm tăng xen phủ obitan, đẩy năng lượng vùng cấm lên cao; ngược lại biến dạng kéo làm giảm khoảng cách dải năng lượng. Dưới biến dạng kéo thích hợp, xảy ra sự dịch chuyển tương đối giữa các cực trị vùng năng lượng, dẫn đến sự chuyển đổi từ vùng cấm gián tiếp sang vùng cấm trực tiếp (thẳng), và tiến tới chuyển pha bán dẫn - kim loại ở các mức biến dạng kéo lớn.

Chương 3: Vật liệu hai chiều bất đối xứng Janus monochalcogenide nhóm III

Chương 3 trình bày các kết quả nghiên cứu về 6 cấu trúc Janus đơn lớp $M_2XY$ ($M = \text{Ga, In}$; $X/Y = \text{S, Se, Te}$; $X \neq Y$).

  • Cấu trúc nguyên tử và độ bền: Cấu trúc Janus $M_2XY$ được hình thành bằng cách thay thế một lớp nguyên tử chalcogen $X$ ở một phía mặt phẳng bằng nguyên tử chalcogen $Y$. Sự bất đối xứng này phá vỡ đối xứng gương $D_{3h}$, chuyển cấu trúc về nhóm đối xứng $C_{3v}$. Phổ phonon DFPT xác nhận không có nhánh âm; mô phỏng AIMD ở $300\text{ K}$ trong suốt thời gian khảo sát cho thấy cấu trúc nguyên tử dao động quanh vị trí cân bằng mà không bị phá hủy liên kết, chứng minh tính ổn định nhiệt động lực học.
  • Tính chất điện tử và phân tách spin Rashba: Sự mất đối xứng gương kết hợp với tương tác spin-quỹ đạo (SOC) làm triệt tiêu sự suy biến spin tại các dải năng lượng gần điểm đối xứng cao, dẫn đến sự phân tách spin kiểu Rashba ($\Delta E_{split}$). Độ rộng vùng cấm tính bằng phiếm hàm PBE luôn thấp hơn so với phiếm hàm HSE06 do hạn chế tự nhiên của GGA.
Vật liệu Janus Hằng số mạng $a$ ($\text{\AA}$) Độ dày $\Delta h$ ($\text{\AA}$) Góc $\angle XMM$ ($^\circ$) Góc $\angle MMY$ ($^\circ$) $E_g^{\text{PBE}}$ ($\text{eV}$) $E_g^{\text{HSE06}}$ ($\text{eV}$) $E_g^{\text{PBE+SOC}}$ ($\text{eV}$) Phân tách $\Delta E_{split}$ ($\text{eV}$)
$\text{Ga}_2\text{SSe}$ 3,67 4,74 117,1 118,5 2,21 3,02 2,17 0,022
$\text{Ga}_2\text{STe}$ 3,82 4,98 115,2 119,8 1,85 2,54 1,73 0,084
$\text{Ga}_2\text{SeTe}$ 3,90 5,10 116,4 119,1 1,72 2,36 1,58 0,091
$\text{In}_2\text{SSe}$ 3,97 5,26 117,0 118,6 1,67 2,38 1,61 0,026
$\text{In}_2\text{STe}$ 4,15 5,50 115,0 120,2 1,32 1,93 1,18 0,098
$\text{In}_2\text{SeTe}$ 4,19 5,62 116,3 119,4 1,22 1,79 1,06 0,105
  • Điều khiển bằng biến dạng cơ học và điện trường ngoài:
    • Biến dạng hai trục ($\varepsilon_{xy}$): Làm thay đổi vị trí tương đối giữa các nhánh năng lượng, tạo ra chuyển tiếp vùng cấm gián tiếp - trực tiếp trong dải biến dạng khảo sát.
    • Điện trường ngoài ($\vec{E}$): Điện trường phân cực cùng chiều hoặc ngược chiều với điện trường nội tại làm tăng hoặc giảm rào thế tĩnh điện, làm dịch chuyển vị trí VBM và CBM, cho phép điều khiển độ rộng vùng cấm một cách liên tục.

Chương 4: Oxy hóa vật liệu hai chiều Janus monochalcogenide nhóm III

Chương 4 khảo sát quá trình oxy hóa hình thành các đơn lớp Janus $\text{GaIn}XO$ ($X = \text{S, Se, Te}$) với 6 cấu hình hình học: 3 cấu trúc $X\text{GaInO}$ và 3 cấu trúc $\text{OGaIn}X$.

  • Đặc trưng cấu trúc, cơ học và năng lượng: Năng lượng gắn kết ($E_{coh}$) của tất cả các cấu hình đều mang giá trị âm lớn (từ $-4{,}8\text{ eV/nguyên tử}$ đến $-5{,}6\text{ eV/nguyên tử}$), khẳng định khả năng tự hình thành thuận lợi về mặt nhiệt động học. Các hằng số đàn hồi $C_{11}, C_{12}, C_{22}, C_{66}$ và mô-đun Young $Y_{2D}$ đều thỏa mãn tiêu chuẩn bền vững cơ học Born-Huang.
  • Tính chất điện tử và sự bất đối xứng công thoát: Sự xuất hiện của lớp nguyên tử oxy có độ âm điện rất lớn tạo ra một điện trường nội tại mạnh vuông góc với mặt phẳng. Điều này dẫn đến sự phân tách mức chân không giữa hai bề mặt ($\Delta \Phi = |\Phi_O - \Phi_X|$), làm công thoát ở mặt chứa oxy ($\Phi_O$) khác biệt rõ rệt so với mặt chứa chalcogen ($\Phi_X$). Một số cấu hình duy trì trạng thái bán dẫn với độ rộng vùng cấm xác định, trong khi một số cấu hình chuyển sang trạng thái kim loại ($M$).
  • Đặc trưng truyền dẫn hạt tải: Luận án tính toán chi tiết khối lượng hiệu dụng $m^*$, mô-đun đàn hồi $C_{2D}$, hằng số thế biến dạng $E_d$ và độ linh động hạt tải $\mu$ dọc theo hai trục trực giao $x$ và $y$. Kết quả cho thấy độ linh động hạt tải có tính dị hướng cao giữa hai phương, và độ linh động của điện tử trong một số cấu hình đạt giá trị vượt trội so với lỗ trống.
Cấu hình Janus Hướng Khối lượng hiệu dụng $m^*/m_0$ Mô-đun đàn hồi $C_{2D}$ ($\text{N/m}$) Thế biến dạng $E_d$ ($\text{eV}$) Độ linh động hạt tải $\mu$ ($\text{cm}^2/\text{Vs}$)
$\text{SGaInO}$ (e) $x$ 0,21 74,52 3,12 864,2
$y$ 0,23 74,52 2,85 920,5
$\text{SeGaInO}$ (e) $x$ 0,19 68,14 2,94 1120,4
$y$ 0,20 68,14 2,71 1250,8
$\text{OGaInS}$ (e) $x$ 0,28 81,20 4,20 345,6
$y$ 0,31 81,20 3,95 362,1

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

Đóng góp mới về mặt khoa học và lý luận

Các kết quả mới được luận án tự nêu bao gồm:

  1. Hệ thống hóa dữ liệu vi mô của đơn lớp $MX$: Xác lập đầy đủ bộ thông số cấu trúc mạng, hằng số đàn hồi, thế tĩnh điện, công thoát và cấu trúc vùng năng lượng của 6 đơn lớp monochalcogenide nhóm III ($MX$) ở trạng thái cân bằng và dưới tác động của biến dạng cơ học hai trục ($\varepsilon_{xy}$).
  2. Làm rõ cơ chế vật lý trong đơn lớp Janus $M_2XY$: Chứng minh rằng việc phá vỡ đối xứng gương theo phương thẳng đứng tạo ra điện trường nội tại định hướng và làm xuất hiện hiện tượng phân tách spin Rashba rõ nét khi xét tương tác spin-quỹ đạo (SOC). Xác lập quy luật điều khiển độ rộng vùng cấm và chuyển tiếp pha vùng cấm gián tiếp - trực tiếp thông qua biến dạng cơ học và điện trường ngoài.
  3. Khám phá tính chất của hệ Janus oxy hóa $\text{GaIn}XO$: Xác định tính ổn định nhiệt động và cơ học của 6 cấu hình $\text{GaIn}XO$, làm sáng tỏ sự phân hóa tính chất bán dẫn - kim loại, sự chênh lệch thế chân không ($\Delta \Phi$) giữa hai bề mặt bất đối xứng, và định lượng độ linh động hạt tải dị hướng theo lý thuyết thế biến dạng.

Kiến nghị và giải pháp đề xuất

  • Đề xuất sử dụng biến dạng cơ học hai trục như một giải pháp kỹ thuật hiệu quả (strain engineering) để chuyển đổi đặc tính vùng cấm từ gián tiếp sang trực tiếp, giúp tăng cường hiệu suất phát xạ quang trong các linh kiện quang điện tử.
  • Đề xuất ứng dụng cấu trúc Janus bất đối xứng ($M_2XY$ và $\text{GaIn}XO$) vào các thiết bị phân tách nước quang xúc tác (photocatalytic water splitting) và cảm biến điện trường, nhờ cơ chế điện trường nội tại tự phát giúp phân tách hiệu quả các cặp điện tử - lỗ trống quang sinh và sự chênh lệch thế hóa khử ở hai bề mặt.
  • Khuyến nghị sử dụng các đơn lớp $\text{SeGaInO}$ với độ linh động điện tử cao ($> 1000\text{ cm}^2/\text{Vs}$) trong thiết kế kênh dẫn cho transistor hiệu ứng trường thế hệ mới.

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

Hạn chế

  • Các mô phỏng trong luận án chủ yếu thực hiện trên các mô hình tinh thể hai chiều lý tưởng, chưa khảo sát các yếu tố thực tế như khuyết tật mạng điểm (vacancies), sai hỏng biên hạt, hoặc ảnh hưởng của các loại đế tựa (substrates) thường dùng trong chế tạo linh kiện.
  • Nghiên cứu tập trung vào trạng thái đơn lớp độc lập, chưa xét đến tương tác ghép tầng trong các cấu trúc nhiều lớp (multilayer) phức tạp.

Hướng nghiên cứu tiếp

  • Mở rộng nghiên cứu các dị cấu trúc van der Waals (vdW heterostructures) ghép nối giữa các đơn lớp Janus monochalcogenide với graphene, TMDs hoặc các vật liệu 2D khác.
  • Khảo sát sâu hơn các tính chất quang học phi tuyến, hiệu ứng áp điện (piezoelectric effect) và hoạt tính xúc tác điện hóa trong môi trường dung dịch cụ thể.
  • Tính toán ảnh hưởng của pha tạp nguyên tử và khuyết tật cấu trúc lên độ bền và cơ chế truyền dẫn hạt tải.

Giá trị tham khảo

  • Đối tượng thụ hưởng:
    • Nghiên cứu sinh và học viên cao học: Thuộc các chuyên ngành Vật lý lý thuyết, Vật lý chất rắn, Khoa học vật liệu và Hóa học tính toán cần tài liệu tham khảo về mô phỏng ab-initio hệ thấp chiều.
    • Giảng viên và nhà nghiên cứu: Làm việc trong lĩnh vực vật liệu nano 2D, mô phỏng cấu trúc dải năng lượng và cơ chế truyền dẫn lượng tử.
    • Các kỹ sư và nhà thực nghiệm bán dẫn: Tìm kiếm các thông số vật lý chuẩn xác (hằng số mạng, công thoát, độ linh động hạt tải, ngưỡng biến dạng) để định hướng thiết kế linh kiện quang điện tử, cảm biến áp điện và spintronics.
  • Phần nội dung đáng tham khảo nhất:
    • Quy trình chi tiết tính toán độ bền cơ học qua hệ số đàn hồi Born-Huang và độ ổn định nhiệt động qua phổ phonon DFPT / mô phỏng AIMD.
    • Phương pháp luận và công thức tính toán độ linh động hạt tải hai chiều dựa trên lý thuyết thế biến dạng Bardeen-Shockley kết hợp DFT.
    • Bảng dữ liệu định lượng về các thông số cấu trúc, công thoát và độ rộng vùng cấm của các hệ đơn lớp $MX$, $M_2XY$ và $\text{GaIn}XO$.

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

1. Tại sao các đơn lớp monochalcogenide nhóm III $MX$ nguyên thủy không xuất hiện hiện tượng phân tách spin Rashba, trong khi cấu trúc Janus $M_2XY$ lại có?

Hiện tượng phân tách spin Rashba đòi hỏi hai điều kiện đồng thời: có tương tác spin-quỹ đạo (SOC) và sự phá vỡ tính đối xứng nghịch đảo không gian hoặc đối xứng gương. Ở đơn lớp $MX$ nguyên thủy, cấu trúc tinh thể có tính đối xứng gương theo phương thẳng đứng ($D_{3h}$), do đó tính suy biến spin được bảo toàn. Ở cấu trúc Janus $M_2XY$, việc thay thế một mặt nguyên tử chalcogen bằng một nguyên tố khác đã phá vỡ đối xứng gương ($C_{3v}$), sinh ra điện trường nội tại định hướng vuông góc, dẫn đến sự phân tách mức spin $\Delta E_{split}$ khi xét hiệu ứng SOC.

2. Biến dạng cơ học hai trục ($\varepsilon_{xy}$) làm thay đổi cấu trúc vùng năng lượng của đơn lớp $MX$ và $M_2XY$ theo cơ chế nào?

Biến dạng nén ($\varepsilon_{xy} < 0$) làm giảm độ dài liên kết, tăng độ xen phủ giữa các obitan nguyên tử, đẩy dải dẫn lên cao và làm rộng vùng cấm. Ngược lại, biến dạng kéo ($\varepsilon_{xy} > 0$) làm giãn mạng tinh thể, làm suy yếu tương tác liên kết và thu hẹp vùng cấm. Đồng thời, biến dạng kéo làm thay đổi tốc độ dịch chuyển năng lượng khác nhau giữa các điểm đối xứng cao (như $\Gamma, K, M$), dẫn đến sự hoán đổi vị trí cực trị VBM/CBM, tạo ra chuyển tiếp từ vùng cấm gián tiếp sang trực tiếp và chuyển pha bán dẫn - kim loại.

3. Ưu điểm của việc sử dụng phiếm hàm lai HSE06 so với phiếm hàm chuẩn GGA-PBE trong luận án là gì?

Phiếm hàm GGA-PBE có xu hướng đánh giá thấp đáng kể độ rộng vùng cấm của chất bán dẫn do hiện tượng tự tương tác phi vật lý và thiếu thành phần trao đổi Hartree-Fock chính xác. Phiếm hàm lai HSE06 khắc phục nhược điểm này bằng cách tích hợp $25%$ thành phần trao đổi Hartree-Fock tầm ngắn kết hợp với phép sàng lọc điện tử, giúp tái hiện chính xác hơn độ rộng vùng cấm thực nghiệm của các hệ vật liệu hai chiều.

4. Sự khác biệt về công thoát ($\Delta \Phi$) giữa hai bề mặt của đơn lớp Janus oxy hóa $\text{GaIn}XO$ có ý nghĩa vật lý gì?

Sự chênh lệch công thoát $\Delta \Phi = |\Phi_O - \Phi_X|$ phản ánh sự tồn tại của một rào thế tĩnh điện và mô-men lưỡng cực nội tại mạnh hướng từ bề mặt này sang bề mặt kia do độ âm điện rất lớn của nguyên tử oxy. Về mặt ứng dụng, độ chênh lệch thế này giúp phân tách tức thời các cặp điện tử - lỗ trống quang sinh về hai phía bề mặt đối diện, ngăn chặn sự tái hợp hạt tải và hỗ trợ tối ưu các phản ứng oxy hóa khử độc lập trong công nghệ quang xúc tác tách nước.

Kết luận

Luận án của tác giả Võ Thị Tuyết Vi đã thực hiện một nghiên cứu lý thuyết hệ thống và toàn diện bằng phương pháp phiếm hàm mật độ (DFT) về cấu trúc, độ bền cơ - nhiệt, tính chất điện tử và đặc trưng truyền dẫn hạt tải của họ vật liệu hai chiều monochalcogenide nhóm III $MX$ cùng các cấu trúc bất đối xứng Janus $M_2XY$ và $\text{GaIn}XO$. Công trình đã làm sáng tỏ quy luật biến đổi cấu trúc vùng năng lượng dưới tác động của biến dạng cơ học hai trục và điện trường ngoài, đồng thời định lượng hiệu ứng phân tách spin Rashba và độ linh động hạt tải dị hướng sinh ra từ sự phá vỡ đối xứng gương. Các kết quả tính toán ab-initio trong luận án cung cấp cơ sở dữ liệu vi mô chuẩn xác và đáng tin cậy, góp phần định hướng cho các nghiên cứu thực nghiệm tiếp theo trong việc thiết kế và chế tạo các linh kiện quang điện tử, spintronics và vật liệu chuyển đổi năng lượng thế hệ mới.