Bối cảnh và vấn đề nghiên cứu

Sự phát triển của công nghệ bán dẫn đóng vai trò nền tảng trong cuộc cách mạng vi điện tử và công nghệ thông tin hiện đại, đặt ra yêu cầu liên tục thu nhỏ kích thước linh kiện mạch tích hợp xuống quy mô nanomet nhằm gia tăng tốc độ truyền tải thông tin và tối ưu hóa năng lượng. Trong bối cảnh đó, nghiên cứu khoa học vật liệu đã có bước ngoặt quan trọng vào năm 2004 khi Novoselov và các cộng sự khám phá ra graphene – vật liệu hai chiều (2D) đầu tiên được tách từ than chì với các đặc tính dẫn điện, dẫn nhiệt và độ bền cơ học vượt trội so với vật liệu khối ba chiều (3D). Tuy nhiên, graphene lại thiếu vùng cấm năng lượng tự nhiên (zero-bandgap), gây khó khăn trong việc đóng ngắt dòng điện trong các transistor hiệu ứng trường.

Năm 2010, các nhóm nghiên cứu độc lập của Heinz và Wang đã phát hiện ra các hợp chất kim loại chuyển tiếp dichalcogenide (Transition Metal Dichalcogenides - TMDC) dạng đơn lớp sở hữu vùng cấm năng lượng trực tiếp. Đặc tính này cho phép các nhà khoa học điều chỉnh dải năng lượng thông qua việc kiểm soát số lượng lớp màng mỏng thay vì phải can thiệp vào thành phần hóa học hay pha tạp nhân tạo. Đồng thời, cấu trúc đơn lớp TMDC ghi nhận hiện tượng phát quang (photoluminescence) mạnh mẽ liên quan trực tiếp đến exciton – một chuẩn hạt lượng tử hình thành từ trạng thái liên kết tĩnh điện giữa electron tự do trên vùng dẫn và lỗ trống ở vùng hóa trị.

Để khảo sát và kiểm soát tính chất điện tử - quang học của TMDC, việc giải phương trình Schrödinger cho exciton hai chiều là phương pháp tiếp cận cốt lõi. Trong các đơn lớp 2D, tương tác Coulomb giữa electron và lỗ trống bị biến đổi do hiệu ứng suy giảm màn chắn điện môi ra môi trường xung quanh. Mô hình thế tương tác thường dùng là thế Rytova-Keldysh; tuy nhiên, hàm thế này chứa các hàm đặc biệt phức tạp (hàm Struve và Bessel), khiến việc tìm nghiệm giải tích gặp nhiều trở ngại và bắt buộc phải giải số. Năm 2022, nhóm tác giả Lý Duy Nhất đã chứng minh khả năng trích xuất các thông số cấu trúc như chiều dài màn chắn, hằng số điện môi trung bình và khối lượng rút gọn từ phổ năng lượng exciton của đơn lớp WSe2 dưới tác động của nhiệt từ. Kế thừa hướng nghiên cứu này, khóa luận "Trích xuất thông tin cấu trúc của đơn lớp TMDC từ phổ năng lượng exciton hai chiều đặt trong từ trường đều" của tác giả Vòng Nguyên Nhật Mẫn được thực hiện nhằm đề xuất một giải pháp tính toán giải tích chính xác bằng cách sử dụng thế Kratzer.

Mục tiêu nghiên cứu của khóa luận được xác định cụ thể gồm ba điểm:

  1. Sử dụng hàm thế tương tác Kratzer trong exciton 2D có khả năng tính toán giải tích để giải quyết bài toán lượng tử.
  2. Sử dụng phổ năng lượng exciton 2D thực nghiệm để trích xuất các thông tin cấu trúc vật lý của đơn lớp TMDC (tập trung vào chất vonfram đisunfua - WS2).
  3. Sử dụng các thông tin vừa trích xuất để dự đoán phổ năng lượng của exciton đặt trong từ trường đều với cường độ từ trường trung bình.

Để hoàn thành các mục tiêu trên, khóa luận triển khai năm nhiệm vụ nghiên cứu:

  • Nhiệm vụ 1: Khảo sát tổng quan về vật liệu TMDC, tính chất exciton và các mô hình thế tương tác màn chắn.
  • Nhiệm vụ 2: Lựa chọn và chuẩn hóa hàm thế Kratzer phục vụ tính toán giải tích cho bài toán exciton hai chiều.
  • Nhiệm vụ 3: Giải phương trình Schrödinger cho exciton 2D trong hai trường hợp: không có từ trường ngoài và có từ trường đều tác dụng.
  • Nhiệm vụ 4: Sử dụng dữ liệu phổ năng lượng thực nghiệm để trích xuất các đại lượng đặc trưng của đơn lớp WS2 (năng lượng vùng cấm, năng lượng liên kết, khối lượng hiệu dụng, hằng số điện môi, bán kính kỳ vọng).
  • Nhiệm vụ 5: Đưa ra mô hình dự đoán phổ năng lượng của exciton trong đơn lớp WS2 dưới từ trường ngoài và so sánh độ chính xác với thực nghiệm.

Đối tượng nghiên cứu của đề tài là chuẩn hạt exciton 2D (trạng thái liên kết electron - lỗ trống) trong vật liệu bán dẫn đơn lớp TMDC, cụ thể là WS2. Phạm vi nghiên cứu về mặt không gian tập trung vào cấu trúc đơn lớp hai chiều đặt trên đế điện môi; phạm vi điều kiện vật lý là trạng thái dừng không từ trường và trạng thái đặt trong từ trường đều cường độ trung bình (từ 0 đến 60 T) tác dụng theo phương vuông góc với mặt phẳng đơn lớp; phạm vi dữ liệu dựa trên phổ năng lượng thực nghiệm quang từ (magneto-photoluminescence) của các trạng thái Rydberg exciton từ 1s đến 5s đã được công bố quốc tế.

Cơ sở lý thuyết và phương pháp

Khung lý thuyết của khóa luận xây dựng trên nền tảng cơ học lượng tử, vật lý chất rắn và lý thuyết chuẩn hạt exciton. Exciton được định nghĩa lần đầu vào năm 1931 bởi Yakov Frenkel như một trạng thái kích thích trung hòa điện di chuyển trong tinh thể dưới dạng sóng. Trong vật lý bán dẫn, mô hình exciton Wannier-Mott miêu tả cặp electron - lỗ trống có bán kính tương tác lớn hơn nhiều so với hằng số mạng tinh thể và năng lượng liên kết yếu (khoảng 10-30 meV trong vật liệu 3D). Tuy nhiên, khi bóc tách về hệ đơn lớp 2D TMDC (gồm một lớp nguyên tử kim loại chuyển tiếp M kẹp giữa hai lớp chalcogen X), hiệu ứng giam giữ lượng tử cùng sự thay đổi phân bố đường sức điện trường ra môi trường điện môi xung quanh khiến hiệu ứng màn chắn suy giảm mạnh. Điều này làm tăng năng lượng liên kết của exciton lên tới khoảng 500 meV.

Về mặt toán học, bài toán hai hạt electron - lỗ trống được đưa về bài toán một hạt chuyển động tương đối thông qua phép biến đổi tách khối tâm:

  • Tọa độ khối tâm: $\vec{R} = \frac{m_e \vec{r}_e + m_h \vec{r}_h}{m_e + m_h}$ với khối lượng toàn phần $M = m_e + m_h$.
  • Tọa độ tương đối: $\vec{r} = \vec{r}_e - \vec{r}_h$ với khối lượng rút gọn $\mu = \frac{m_e m_h}{m_e + m_h}$.

Khi đặt trong từ trường đều $\vec{B} = (0, 0, B)$ theo phương trục Oz, thế vectơ chọn ở chuẩn đối xứng $\vec{A}_i = \frac{1}{2} (\vec{B} \times \vec{r}_i)$. Phương trình Schrödinger tương đối xuất hiện thêm toán tử mômen động lượng $L_z = -i\hbar (x \frac{\partial}{\partial y} - y \frac{\partial}{\partial x})$ và thành phần năng lượng nghịch từ tỷ lệ với $B^2 r^2$.

Để khắc phục nhược điểm phức tạp của hàm thế Rytova-Keldysh $V_{RK}(r) = -\frac{\pi e^2}{2 \kappa r_0} [H_0(r/r_0) - Y_0(r/r_0)]$, tác giả áp dụng hàm thế Kratzer - dạng thế ban đầu dùng mô tả phổ dao động quay của phân tử hai nguyên tử: $$V(r) = -\frac{A}{r} + \frac{B}{r^2} + C$$ Đối với bài toán exciton Wannier-Mott, hàm thế được điều chỉnh thành dạng thế Coulomb màn chắn: $$V(r) = -\frac{e^2}{4\pi\varepsilon_0 \kappa r} + \frac{\delta e^2 r_0}{4\pi\varepsilon_0 \kappa r^2}$$ Chuyển sang hệ đơn vị không thứ nguyên (độ dài đo theo bán kính Bohr hiệu dụng $a_0 = \frac{4\pi\varepsilon_0 \hbar^2}{\mu e^2}$, năng lượng đo theo Rydberg hiệu dụng $R_y = \frac{\mu e^4}{16\pi^2\varepsilon_0^2 \hbar^2}$), hàm thế Kratzer có dạng ngắn gọn: $$V(r) = -\frac{2}{\kappa r} + \frac{\xi_0}{2 r^2}$$ với $\xi_0 = \frac{2\delta r_0}{a_0}$ là tham số hiệu chỉnh không thứ nguyên đặc trưng cho mức độ chắn điện môi.

Phương pháp nghiên cứu kết hợp giữa giải tích thuần túy, phương pháp số học và lý thuyết nhiễu loạn:

  • Phương pháp giải tích chính xác: Khi không có từ trường ngoài, phương trình vi phân được chuyển sang tọa độ cực $(r, \phi)$, tách biến hàm sóng $\Psi(r, \phi) = R(r)\Phi(\phi)$ với $\Phi(\phi) = \frac{1}{\sqrt{2\pi}} e^{im\phi}$. Khử các điểm kỳ dị tại gốc tọa độ và vô cùng, đưa phương trình bán kính về phương trình vi phân hypergeometric suy biến để tìm nghiệm dưới dạng đa thức Laguerre liên kết $L_{n-|m|-1}^{2\xi}(r)$.
  • Phương pháp lý thuyết nhiễu loạn dừng (Stationary Perturbation Theory): Tính toán bổ chính năng lượng đến bậc II đối với thành phần từ trường.
  • Phương pháp nhiễu loạn Hypervirial (Hypervirial Perturbation Method - HPM): Kết hợp mối quan hệ siêu vi giữa các toán tử và định lý Hellmann-Feynman để xây dựng hệ thức truy hồi cho các mômen tọa độ $\langle q^j \rangle$, giải trực tiếp các hệ số nhiễu loạn bậc cao mà không cần tích phân hàm sóng trung gian.
  • Nguồn dữ liệu: Tác giả sử dụng số liệu thực nghiệm về phổ quang phát quang phân cực của đơn lớp WS2 trong từ trường từ 0 đến khoảng 65 T của nhóm Goryca, Crooker và Chernikov (công bố trên các tạp chí Nature Communications và Physical Review Letters). Tọa độ đỉnh phổ năng lượng của các trạng thái Rydberg $ns$ ($n = 1, 2, 3, 4, 5$) tại từ trường $B = 0$ được số hóa bằng phần mềm Origin với độ chính xác 5 chữ số thập phân.

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

CHƯƠNG 1: CƠ SỞ LÝ THUYẾT

Chương 1 hệ thống hóa cơ sở vật lý của exciton và vật liệu TMDC:

  • Khái niệm và tính chất exciton: Làm rõ cơ chế hình thành exciton trực tiếp qua chuyển dời quang học từ đỉnh dải hóa trị lên đáy dải dẫn. Tác giả phân loại exciton theo hai tiêu chí: (1) Mức độ hoạt động quang học: Exciton sáng (Bright Exciton - BE) có thời gian sống ngắn, nhạy cảm với ánh sáng do electron và lỗ trống cùng spin, cùng động lượng; Exciton tối (Dark Exciton - DE) có thời gian sống dài do bị cấm chuyển dời bởi spin (spin-forbidden) hoặc mômen động lượng (momentum-forbidden). (2) Bản chất liên kết vật liệu: Phân biệt Frenkel exciton trong điện môi và Wannier-Mott exciton trong bán dẫn.
  • Đơn lớp TMDC: Mô tả cấu trúc mạng lục giác gồm một lớp kim loại chuyển tiếp (W, Mo, Ti...) liên kết cộng hóa trị với hai lớp nguyên tử chalcogen (S, Se, Te), các lớp xếp chồng liên kết yếu qua lực Van der Waals. Cấu trúc này cho phép bóc tách thành đơn lớp 2D, ống nano (1D), dải nano (1D) hoặc chấm lượng tử (0D).
  • Phép tách khối tâm: Trình bày chi tiết toán tử Hamilton cho hệ exciton 2D trong hai trường hợp: không có từ trường và đặt trong từ trường đều theo phương thẳng đứng. Thiết lập phương trình chuyển động tương đối độc lập với tọa độ khối tâm.
  • Khảo sát các mô hình thế: Phân tích sự biến đổi của đường sức điện trường khi chiều không gian giảm từ 3D xuống 2D, dẫn đến sự khác biệt giữa hằng số điện môi trong lớp $\varepsilon_{2D}$ và ngoài lớp $\varepsilon_{env}$. Tác giả so sánh định lượng thế Coulomb thuần túy, thế Rytova-Keldysh và thế Kratzer. Đồ thị so sánh cho thấy thế Kratzer với tham số $\xi_0 \approx 0.21$ tiệm cận rất sát với thế Keldysh ở khoảng cách lớn hơn hằng số mạng $a = 3.25\ \text{Å}$ của đơn lớp TMDC trong chất nền boron nitride lục giác (hBN, $r_0 = 10\ \text{Å}$), khẳng định tính khả thi của việc thay thế thế Keldysh bằng thế Kratzer để giải giải tích.

CHƯƠNG 2: EXCITON 2D TRONG ĐƠN LỚP TMDC DƯỚI TÁC DỤNG THẾ KRATZER

Chương 2 tập trung giải giải tích bài toán exciton 2D khi không có từ trường ngoài tác dụng và trích xuất các thông số cơ bản:

  • Giải phương trình Schrödinger: Trong tọa độ cực, đặt $\xi = \sqrt{m^2 + \xi_0}$, phương trình bán kính được đưa về dạng vi phân chính tắc: $$r g''(r) + (2\xi + 1 - r) g'(r) + \left[ \frac{2\lambda}{\kappa} - \xi \right] g(r) = 0$$ Nghiệm được biểu diễn qua chuỗi lũy thừa $g(r) = \sum_{j=0}^{\infty} a_j r^j$. Hệ thức truy hồi cho các hệ số khai triển được xác lập: $$a_{j+1} = \frac{j + \xi - \frac{2\lambda}{\kappa}}{(j + 1)(j + 2\xi + 1)} a_j$$
  • Điều kiện ngắt chuỗi: Để hàm sóng hội tụ tại vô cùng, chuỗi lũy thừa phải ngắt thành đa thức bậc $n - |m| - 1$ (với $n$ là số lượng tử chính, $m$ là số lượng tử từ). Từ đó, tác giả tìm ra biểu thức năng lượng dạng không thứ nguyên: $$E = -\frac{1}{2\kappa^2 \left( n - |m| - \frac{1}{2} + \sqrt{m^2 + \xi_0} \right)^2}$$
  • Hàm sóng và các kỳ vọng bán kính: Hàm sóng chuẩn hóa $\Psi_{nm}(r, \phi)$ được giải tích thông qua đa thức Laguerre liên kết. Từ hệ số chuẩn hóa $N_{nm}$, tác giả dẫn xuất công thức giải tích cho bán kính trung bình $\langle r \rangle$ (Bảng 2) và bán kính trung bình bình phương $\langle r^2 \rangle$ (Bảng 3) cho các trạng thái $n = 1, 2, 3$.
  • Trích xuất thông số phổ thực nghiệm: Năng lượng đỉnh quang học tuân theo biểu thức $E_{opt} = E_{gap} + E_{bind}$. Tác giả lấy giá trị năng lượng đỉnh thực nghiệm của WS2 tại các trạng thái 1s, 2s, 3s từ phổ của Chernikov et al. Lập tỉ số hiệu năng lượng giữa các trạng thái để triệt tiêu năng lượng vùng cấm $E_{gap}$, phương trình chỉ còn phụ thuộc vào tham số $\xi_0$, cho kết quả: $$\xi_0 = 0.36$$ Từ $\xi_0$, tác giả thiết lập mối quan hệ phụ thuộc giữa hằng số điện môi trung bình $\kappa$ và khối lượng rút gọn $\mu$: $$\frac{\kappa^2}{\mu} \approx \frac{7.60681 e^4}{32\pi^2\varepsilon_0^2\hbar^2}$$ Kết quả này cho phép tính toán trực tiếp năng lượng vùng cấm $E_{gap}$ và năng lượng liên kết exciton $E_{bind}$ của WS2 ở các trạng thái từ 1s đến 5s, đạt độ tương thích cao với mô hình số học dùng thế Keldysh trong tài liệu tham khảo [28] (Bảng 4).

CHƯƠNG 3: EXCITON 2D TRONG ĐƠN LỚP TMDC ĐẶT TRONG TỪ TRƯỜNG ĐỀU DƯỚI TÁC DỤNG THẾ KRATZER

Chương 3 giải quyết bài toán exciton 2D trong từ trường đều và so sánh các phương pháp tính toán:

  • Xây dựng bài toán nhiễu loạn: Đưa phương trình Schrödinger trong từ trường về dạng không thứ nguyên, chứa toán tử mômen động lượng $L_z$ và thành phần từ trường bậc hai $\lambda r^2$.
  • Phương pháp lý thuyết nhiễu loạn dừng: Tác giả tính toán bổ chính năng lượng đến bậc II cho trạng thái 1s. Bổ chính bậc I được tính qua tích phân kỳ vọng $\langle \Psi | r^2 | \Psi \rangle$, trong khi bổ chính bậc II yêu cầu lấy tổng trên các trạng thái kích thích $n' \neq n$ (tính đến $n' = 4$).
  • Phương pháp nhiễu loạn Hypervirial (HPM): Sử dụng toán tử chuyển đổi $r = \alpha q$ để đưa Hamiltonian về dạng $H = p^2 + V(q) + \lambda q^2$. Bằng cách lựa chọn toán tử trung gian $O = q^{j+1} \frac{d}{dq} - q^j \frac{d^2}{dq^2}$, áp dụng định lý Hellmann-Feynman và hệ thức siêu vi $\langle \Psi | [H, O] | \Psi \rangle = 0$, tác giả xây dựng công thức truy hồi tuyến tính trực tiếp cho các hệ số năng lượng nhiễu loạn bậc I, II, III mà không cần tính hàm sóng bổ chính.
  • Đánh giá và so sánh hai phương pháp: Dữ liệu Bảng 5 chứng minh phương pháp Hypervirial perturbation cho kết quả năng lượng trạng thái 1s có độ hội tụ và độ chính xác vượt trội so với lý thuyết nhiễu loạn dừng bậc II khi so sánh với phổ thực nghiệm từ 0 đến 60 T.
  • Trích xuất thông tin cấu trúc và dự đoán phổ: Sử dụng giá trị khối lượng hiệu dụng exciton $\mu = 0.2 m_0$ từ tài liệu thực nghiệm, tác giả xác định được hằng số điện môi $\kappa$ và bán kính trung bình bình phương $\sqrt{\langle r^2 \rangle}$ của các trạng thái. Cuối cùng, tác giả ngoại suy đường cong phổ năng lượng exciton WS2 dưới từ trường đến 60 T cho các trạng thái Rydberg (Hình 8 và Bảng 7).

Kết quả và đóng góp

Khóa luận đã định lượng chính xác các thông số vật lý của đơn lớp WS2 dựa trên mô hình giải tích thế Kratzer kết hợp phổ thực nghiệm. Kết quả so sánh giữa hai phương pháp nhiễu loạn trong từ trường đều tại trạng thái cơ bản 1s được trích xuất từ văn bản:

Từ trường $B$ (T) Năng lượng thực nghiệm $E_{exp}$ (eV) Kết quả tính bằng Nhiễu loạn dừng (eV) Sai số Nhiễu loạn dừng (eV) Kết quả tính bằng Hypervirial (eV) Sai số Hypervirial (eV)
0 2.06450 2.06450 0.00000 2.06450 0.00000
10 2.06465 2.06478 +0.00013 2.06466 +0.00001
20 2.06510 2.06562 +0.00052 2.06514 +0.00004
30 2.06585 2.06702 +0.00117 2.06593 +0.00008
40 2.06690 2.06898 +0.00208 2.06704 +0.00014
50 2.06825 2.07150 +0.00325 2.06846 +0.00021
60 2.06990 2.07458 +0.00468 2.07019 +0.00029

Bảng số liệu trên chỉ rõ phương pháp Hypervirial perturbation duy trì sai số ở mức dưới $0.3\ \text{meV}$ tại từ trường cực đại 60 T, trong khi phương pháp nhiễu loạn dừng bậc II có sai số tích lũy lên tới gần $4.7\ \text{meV}$.

Các kết quả định lượng chính mà tác giả đã đạt được:

  • Xác định tham số hiệu chỉnh không thứ nguyên của thế Kratzer cho WS2: $\xi_0 \approx 0.36$.
  • Thiết lập biểu thức giải tích tường minh cho phổ năng lượng và bán kính kỳ vọng $\langle r \rangle$, $\langle r^2 \rangle$ của exciton 2D theo đa thức Laguerre liên kết.
  • Tính toán năng lượng liên kết exciton $E_{bind}$ và độ rộng vùng cấm $E_{gap}$ của đơn lớp WS2 khi không có từ trường, cho kết quả tương đương với mô hình thế Keldysh giải số trong tài liệu [28].
  • Trích xuất thành công khối lượng hiệu dụng exciton $\mu \approx 0.2 m_0$, tham số hằng số điện môi trung bình $\kappa$ và kích thước không gian $\sqrt{\langle r^2 \rangle}$ của các trạng thái 1s, 2s, 3s.
  • Dự đoán chính xác phổ năng lượng exciton của WS2 dưới từ trường đều từ 0 đến 60 T cho các trạng thái 1s và 2s.

Đóng góp khoa học tự nêu của tác giả:

  • Chứng minh tính ứng dụng hiệu quả của thế Kratzer trong bài toán exciton 2D của vật liệu TMDC. Thay vì phải giải số phức tạp như thế Rytova-Keldysh, thế Kratzer cho phép thu được nghiệm giải tích hoàn toàn mà vẫn đảm bảo độ chính xác vật lý cao.
  • Tiên phong đưa phương pháp nhiễu loạn Hypervirial kết hợp định lý Hellmann-Feynman vào giải bài toán exciton trong đơn lớp TMDC đặt dưới từ trường ngoài, tối ưu hóa tốc độ và độ chính xác tính toán so với phương pháp nhiễu loạn dừng truyền thống.

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

Khóa luận tự xác định các hạn chế kỹ thuật xuất phát từ bản chất mô hình toán học và phạm vi nghiên cứu:

  • Do sử dụng hàm thế Kratzer ở dạng quy ước không thứ nguyên, mô hình không trích xuất trực tiếp được đại lượng chiều dài màn chắn $r_0$ như khi sử dụng hàm thế Keldysh. Nguyên nhân là do chưa thể xác định độc lập hệ số hiệu chỉnh $\delta$ trong biểu thức liên hệ $\xi_0 = \frac{2\delta r_0}{a_0}$.
  • Phạm vi tương thích của mô hình dự đoán phổ năng lượng trong từ trường bị thu hẹp dần ở các trạng thái kích thích cao: Trong khi trạng thái 1s và 2s giữ được độ chính xác cao trong toàn dải từ trường từ 0 đến 60 T, thì trạng thái 3s chỉ dự đoán tốt trong khoảng từ 0 đến 20 T, và trạng thái 4s chỉ bảo đảm độ chính xác trong khoảng từ 0 đến 5 T (do hiệu ứng xen phủ và các bậc nhiễu loạn cao chưa xét đến).

Hướng nghiên cứu tiếp theo được tác giả đề xuất:

  • Khảo sát và áp dụng các dạng biến thể khác của hàm thế Kratzer (như thế Kratzer sửa đổi - Modified Kratzer Potential, thế Kratzer-Fues) nhằm nâng cao độ chính xác ở các trạng thái kích thích cao.
  • Tìm kiếm và thử nghiệm các dạng thế Coulomb màn chắn khác có khả năng giải tích nhưng vẫn giữ được thông số độ dài màn chắn $r_0$.
  • Mở rộng phạm vi tính toán và kiểm chứng mô hình trên các vật liệu TMDC đơn lớp khác ngoài WS2 (chẳng hạn như MoS2, MoSe2, WSe2).

Giá trị tham khảo

Bản khóa luận tốt nghiệp là tài liệu tham khảo có giá trị chuyên môn cao đối với:

  • Sinh viên, học viên cao học và nghiên cứu sinh thuộc các chuyên ngành Vật lý lý thuyết, Vật lý tính toán, Vật lý chất rắn và Khoa học vật liệu nano đang thực hiện các đề tài liên quan đến cấu trúc dải năng lượng và tính chất quang học của vật liệu hai chiều.
  • Giảng viên giảng dạy học phần Cơ học lượng tử nâng cao hoặc Vật lý bán dẫn thấp chiều: Khóa luận cung cấp một ví dụ mẫu mực về việc áp dụng đa thức Laguerre liên kết giải phương trình vi phân và kỹ thuật sử dụng định lý Hellmann-Feynman trong bài toán nhiễu loạn.
  • Các phần đáng tham khảo trực tiếp:
    • Quy trình biến đổi toán học đưa phương trình vi phân tọa độ cực về dạng chuẩn tắc của đa thức Laguerre liên kết trong Chương 2.
    • Thuật toán truy hồi của phương pháp nhiễu loạn Hypervirial (HPM) trong Chương 3, có thể trực tiếp chuyển hóa thành mã lập trình tính toán đại số trên máy tính.
    • Phương pháp lập tỉ số triệt tiêu tham số chưa biết để trích xuất các hằng số vật lý từ dữ liệu thực nghiệm quang học.

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

1. Tại sao tác giả lựa chọn thế Kratzer thay vì thế Rytova-Keldysh vốn phổ biến hơn khi nghiên cứu exciton 2D trong TMDC? Thế Rytova-Keldysh mô tả chính xác hiệu ứng màn chắn điện môi 2D nhưng biểu thức toán học chứa hàm Struve và hàm Bessel bậc không, dẫn đến việc không thể giải giải tích chính xác phương trình Schrödinger mà phải dùng phương pháp giải số phức tạp. Ngược lại, thế Kratzer có dạng đại số đơn giản hơn (tỉ lệ nghịch với $r$ và $r^2$), cho phép tìm ra nghiệm giải tích chính xác dưới dạng đa thức Laguerre liên kết, giúp việc tính toán phổ năng lượng và trích xuất thông số vật lý trở nên thuận tiện và nhanh chóng mà sai số so với thực nghiệm là không đáng kể.

2. Giá trị tham số hiệu chỉnh $\xi_0$ của đơn lớp WS2 được xác định như thế nào trong khóa luận? Tham số $\xi_0$ được xác định thuần túy từ phổ thực nghiệm tại từ trường $B = 0$. Bằng cách lập tỉ số hiệu năng lượng giữa các trạng thái kích thích: $\frac{E_{2s} - E_{1s}}{E_{3s} - E_{2s}}$, thành phần năng lượng vùng cấm $E_{gap}$ và hằng số thứ nguyên năng lượng tự triệt tiêu. Phương trình thu được chỉ còn chứa ẩn số duy nhất là $\xi_0$. Thay các giá trị năng lượng thực nghiệm đọc được từ phổ của Chernikov et al. vào phương trình, tác giả tính toán được $\xi_0 \approx 0.36$.

3. Vì sao phương pháp nhiễu loạn Hypervirial (HPM) lại vượt trội hơn phương pháp lý thuyết nhiễu loạn dừng bậc II? Phương pháp nhiễu loạn dừng truyền thống đòi hỏi phải tính tích phân các yếu tố ma trận chuyển dời giữa trạng thái cơ bản và vô hạn các trạng thái kích thích khác nhau; việc chỉ lấy đến trạng thái $n' = 4$ khiến sai số tích lũy tăng nhanh khi từ trường tăng (sai số đạt tới 4.68 meV tại 60 T). Trong khi đó, phương pháp Hypervirial sử dụng định lý Hellmann-Feynman và hệ thức giao hoán siêu vi để xây dựng hệ thức truy hồi đại số trực tiếp cho các hệ số nhiễu loạn đến bậc cao (bậc III) mà không cần tích phân hàm sóng, giúp việc tính toán đơn giản, hội tụ nhanh và duy trì sai số cực nhỏ (chỉ 0.29 meV tại 60 T).

4. Khóa luận dự đoán phổ năng lượng exciton trong từ trường đạt độ chính xác trong những phạm vi nào? Mô hình tính toán của khóa luận dự đoán rất tốt phổ năng lượng exciton trong dải từ trường từ 0 đến 60 T đối với hai trạng thái 1s và 2s (đường lý thuyết khớp sát với điểm thực nghiệm). Đối với trạng thái kích thích cao 3s, mô hình chỉ đảm bảo độ chính xác trong khoảng từ trường từ 0 đến 20 T. Đối với trạng thái 4s, độ chính xác chỉ duy trì được trong khoảng từ trường hẹp từ 0 đến 5 T; vượt qua ngưỡng này, các hiệu ứng phi tuyến bậc cao bắt đầu làm sai lệch kết quả tính toán.

5. Vì sao mô hình trong khóa luận không trích xuất trực tiếp được chiều dài màn chắn $r_0$ của đơn lớp WS2? Trong biểu thức hàm thế Kratzer không thứ nguyên, tham số hiệu chỉnh gắn liền với công thức $\xi_0 = \frac{2\delta r_0}{a_0}$. Mặc dù đã tính ra được giá trị số của $\xi_0 \approx 0.36$ và bán kính Bohr hiệu dụng $a_0$, nhưng hệ số điều chỉnh hình học $\delta$ chưa được xác định độc lập từ mô hình lý thuyết. Do đó, tác giả không thể tách riêng giá trị của độ dài màn chắn $r_0$ nếu không giả định trước hệ số $\delta$.

Kết luận

Khóa luận tốt nghiệp của tác giả Vòng Nguyên Nhật Mẫn đã giải quyết thành công bài toán exciton hai chiều trong đơn lớp bán dẫn TMDC WS2 bằng việc áp dụng hàm thế tương tác Kratzer có nghiệm giải tích chính xác. Thông qua việc kết hợp kỹ thuật đa thức Laguerre liên kết khi không có từ trường và phương pháp nhiễu loạn Hypervirial khi đặt trong từ trường đều, công trình đã xác lập một quy trình tính toán hiệu quả, hội tụ cao để trích xuất các thông số cấu trúc vi mô quan trọng như năng lượng vùng cấm, năng lượng liên kết, khối lượng hiệu dụng và kích thước exciton. Đây là một công trình nghiên cứu vật lý lý thuyết và vật lý tính toán hoàn chỉnh, mang lại giá trị phương pháp luận rõ ràng cho việc khảo sát tính chất quang - điện tử của các vật liệu hai chiều thế hệ mới.