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

Trong vật lý chất rắn và khoa học vật liệu, các vật liệu hai chiều (2D materials) sở hữu cấu trúc lớp Van der Waals như họ vật liệu bán dẫn chuyển tiếp kim loại hai nguyên tử chalcogen (Transition Metal Dichalcogenides - TMDCs với công thức $MX_2$, trong đó $M = \text{Mo, W}$ và $X = \text{S, Se}$) thu hút sự quan tâm lớn nhờ các tính chất quang - điện tử và spin đặc thù ở quy mô đơn lớp nguyên tử. Khi chuyển từ dạng khối sang cấu trúc đơn lớp 2D, sự giam giữ không gian của các hạt tải điện trong mặt phẳng cùng sự suy giảm của hiệu ứng màn chắn điện môi làm gia tăng mạnh mẽ tương tác Coulomb, dẫn đến việc hình thành các cặp electron - lỗ trống liên kết chặt chẽ (exciton) với năng lượng liên kết lớn, chi phối các đặc tính quang học của vật liệu ngay ở nhiệt độ phòng.

Tương tác giữa electron và lỗ trống trong các đơn lớp TMDCs ($MoS_2, MoSe_2, WS_2, WSe_2$) không còn tuân theo thế Coulomb thuần túy mà được mô tả chuẩn xác thông qua thế Keldysh—thế tương tác phụ thuộc vào hằng số điện môi trung bình $\kappa$, độ dài chắn $r_0$ (liên quan đến độ phân cực của lớp màng) và khối lượng hiệu dụng rút gọn $\mu$. Về mặt lý thuyết, phổ năng lượng exciton trong TMDCs đã được giải số với độ chính xác cao (đạt từ 12 đến 20 chữ số thập phân bằng phương pháp toán tử FK). Tuy nhiên, việc xây dựng một biểu thức giải tích tường minh cho năng lượng exciton từ dạng nguyên bản của thế Keldysh vẫn là một bài toán mở. Các nghiên cứu trước đây thường phải sử dụng các dạng xấp xỉ gần đúng của thế Keldysh hoặc chuyển sang thế Kratzer, dẫn đến độ chính xác chưa cao. Biểu thức giải tích tường minh mang ý nghĩa thiết thực trong việc trích xuất nhanh các thông số cấu trúc vật liệu TMDCs từ dữ liệu quang phổ thực nghiệm mà không cần phải thực hiện các thuật toán lập trình giải số phức tạp cho từng chất riêng biệt.

Mục tiêu tổng quát của khóa luận là xây dựng thành công biểu thức giải tích cho năng lượng exciton ở trạng thái cơ bản 1s trong đơn lớp TMDCs dựa trên thế tương tác Keldysh nguyên bản thông qua phương pháp toán tử FK và lý thuyết nhiễu loạn Rayleigh-Schrödinger có điều tiết.

Các nhiệm vụ nghiên cứu cụ thể được xác định bao gồm:

  1. Xây dựng phương trình Schrödinger cho chuyển động khối tâm và chuyển động tương đối của hệ electron - lỗ trống trong đơn lớp TMDCs.
  2. Nghiên cứu phương pháp toán tử FK và lý thuyết nhiễu loạn có điều tiết; thiết lập biểu thức tính toán các bổ chính năng lượng và hàm sóng đến bậc ba.
  3. Xây dựng biểu thức giải tích tường minh cho năng lượng exciton trạng thái 1s và kiểm chứng độ chính xác trên 4 vật liệu đơn lớp TMDCs phổ biến.

Đối tượng nghiên cứu là mức năng lượng exciton trạng thái 1s của hệ hai chiều. Phạm vi nghiên cứu tập trung vào 4 loại đơn lớp TMDCs gồm $WSe_2, WS_2, MoSe_2, MoS_2$ ở trạng thái cơ bản ($n=0, m=0$), chưa xét đến ảnh hưởng của từ trường ngoài, hiệu ứng nhiệt độ và cấu trúc spin của exciton.

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

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

Khóa luận mô tả exciton trong đơn lớp TMDCs bằng phương trình Schrödinger hai hạt hai chiều: $$\hat{H}\psi = E\psi, \quad \hat{H} = -\frac{\hbar^2}{2m_e}\nabla_e^2 - \frac{\hbar^2}{2m_h}\nabla_h^2 + V_{eh}(\vec{r}_e - \vec{r}h)$$ Trong đó, tương tác giữa electron và lỗ trống được mô tả bằng thế Keldysh nguyên bản: $$V{eh}(r) = -\frac{\pi e^2}{8\varepsilon_0 \kappa r_0} \left[ H_0\left(\frac{r}{r_0}\right) - Y_0\left(\frac{r}{r_0}\right) \right]$$ Với $H_0$ và $Y_0$ lần lượt là hàm Struve và hàm Bessel (Neumann) loại hai bậc không; $\kappa$ là hằng số điện môi trung bình và $r_0$ là độ dài chắn điện môi.

Bằng cách chuyển sang hệ tọa độ khối tâm ($\vec{R}$) và tọa độ tương đối ($\vec{r}$), bài toán được tách thành chuyển động khối tâm tự do và chuyển động tương đối. Bỏ qua động năng khối tâm để tập trung vào năng lượng liên kết, phương trình được chuẩn hóa không thứ nguyên bằng việc đưa vào bán kính Bohr hiệu dụng $a_0^* = \frac{4\pi\varepsilon_0\hbar^2}{\mu e^2}$ và đơn vị năng lượng Hartree $2R_y^* = \frac{\hbar^2}{\mu {a_0^*}^2}$.

Thế Keldysh được biểu diễn lại qua tích phân biến đổi Laplace: $$V_K(r) = -\frac{1}{\kappa r_0} \int_0^\infty \frac{dq}{1+q^2} e^{-\frac{qr}{\alpha a_0^}}$$ với $\alpha = \frac{r_0}{a_0^}$ là tham số độ dài chắn không thứ nguyên.

Khóa luận áp dụng phép biến đổi chính tắc Levi-Civita trong không gian hai chiều ($x = u^2 - v^2, y = 2uv, r = u^2 + v^2$) để chuyển phương trình Schrödinger của exciton về bài toán dao động tử phi điều hòa hai chiều trong không gian $(u, v)$.

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

  1. Phương pháp toán tử FK (Fock-Krylov Algebraic Operator Method): Thiết lập hệ toán tử sinh - hủy $a, a^\dagger, b, b^\dagger$ với việc đưa vào tham số biến phân tự do $\omega > 0$ nhằm tối ưu hóa tốc độ hội tụ. Xây dựng hệ hàm cơ sở trực chuẩn $|k, m\rangle$ của dao động tử hai chiều và tính toán tường minh các yếu tố ma trận của toán tử động năng $T_{jk}$, toán tử tọa độ $R_{jk}$ và toán tử thế Keldysh $V_{jk}$ thông qua toán tử trung gian $\hat{O} = e^{-\frac{q}{\omega r_0}(M^\dagger + M + N)}$.
  2. Lý thuyết nhiễu loạn Rayleigh-Schrödinger có điều tiết (Modulated Perturbation Method - MPM): Phân tách Hamiltonian và toán tử chuẩn hóa thành phần trung hòa và phần nhiễu loạn: $\hat{H} = \hat{H}_0 + \beta \hat{V}$, $\hat{R} = \hat{R}_0 + \beta \hat{R}_1$. Thiết lập hệ phương trình xác định bổ chính bậc không ($E^{(0)}$), bậc một ($\Delta E^{(1)} = 0$), bậc hai ($\Delta E^{(2)}$) và phát triển mở rộng cho bổ chính bậc ba ($\Delta E^{(3)}$).
  3. Kỹ thuật tính toán giải tích: Sử dụng phép khai triển chuỗi Taylor cho các biểu thức bổ chính năng lượng quanh điểm cực trị $x = 0.75$ (với $x = \frac{\alpha}{\omega}$). Các phép biến đổi đại số hình thức phức tạp được hỗ trợ xử lý bằng phần mềm tính toán ký hiệu Wolfram Mathematica.

Nguồn dữ liệu thông số vật liệu (khối lượng hiệu dụng electron, lỗ trống và độ dài chắn) của $WSe_2, WS_2$ được trích xuất từ công trình của Duy-Nhat Ly và cộng sự (2021); dữ liệu của $MoSe_2, MoS_2$ được lấy từ công trình của Goryca và cộng sự (2019).

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

Chương 1. Bài toán năng lượng exciton trong đơn lớp TMDCs

Chương 1 thiết lập mô hình toán học lượng tử cho exciton hai chiều:

  • Tách chuyển động của hệ electron - lỗ trống thành chuyển động khối tâm $\vec{R} = \frac{m_e \vec{r}_e + m_h \vec{r}_h}{m_e + m_h}$ và chuyển động tương đối $\vec{r} = \vec{r}_e - \vec{r}_h$. Rút gọn phương trình vi phân chuyển động tương đối về dạng không thứ nguyên.
  • Thực hiện biến đổi Laplace cho thế Keldysh và áp dụng biến đổi tọa độ Levi-Civita đưa phương trình về dạng dao động tử phi điều hòa: $$\left[ -4\left(\frac{\partial^2}{\partial u^2} + \frac{\partial^2}{\partial v^2}\right) + \dots \right]\psi = 0$$
  • Định nghĩa các toán tử sinh - hủy $a, a^\dagger, b, b^\dagger$ phụ thuộc tham số điều tiết $\omega$: $$a = \frac{1}{2}\left[\sqrt{\omega}u + \frac{1}{\sqrt{\omega}}\frac{\partial}{\partial u} - i\left(\sqrt{\omega}v + \frac{1}{\sqrt{\omega}}\frac{\partial}{\partial v}\right)\right]$$ $$b = \frac{1}{2}\left[\sqrt{\omega}u - \frac{1}{\sqrt{\omega}}\frac{\partial}{\partial u} - i\left(\sqrt{\omega}v - \frac{1}{\sqrt{\omega}}\frac{\partial}{\partial v}\right)\right]$$
  • Xây dựng hệ hàm sóng trạng thái cơ sở $|k, m\rangle = \frac{(a^\dagger)^{k+m}(b^\dagger)^{k-m}}{\sqrt{(k+m)!(k-m)!}}|0(\omega)\rangle$ với các chỉ số lượng tử $m = 0, \pm 1, \pm 2, \dots$ và $k = |m|, |m|+1, \dots$
  • Tính toán ma trận các toán tử trung gian $\hat{N} = a^\dagger a + b^\dagger b + 1$, $\hat{M}^\dagger = a^\dagger b^\dagger$, $\hat{M} = ab$ và dẫn xuất công thức giải tích cho các yếu tố ma trận động năng $T_{jk}$, toán tử $R_{jk}$ và thế tương tác $V_{jk}$.

Chương 2. Lý thuyết nhiễu loạn Rayleigh-Schrödinger cho bài toán exciton trong đơn lớp TMDCs

Chương 2 trình bày khung giải tích của lý thuyết nhiễu loạn áp dụng cho hệ toán tử tổng quát:

  • Phân tách ma trận toán tử thành các phần tử trung hòa ($j = k$) và không trung hòa ($j \neq k$):
    • Với $j = k$: $h_j = -\frac{\omega}{8}T_{jj} + \frac{V_{jj}}{\kappa}$, $r_j = \frac{R_{jj}}{\omega}$.
    • Với $j \neq k$: $v_{jk} = h_{jk} = -\frac{\omega}{8}T_{jk} + \frac{V_{jk}}{\kappa}$, $r_{jk} = \frac{R_{jk}}{\omega}$.
  • Thiết lập hệ thức truy hồi cho các bậc nhiễu loạn:
    • Bổ chính bậc không vào năng lượng: $E_{nm}^{(0)} = \frac{h_{nn}}{r_{nn}}$.
    • Bổ chính bậc một: $\Delta E_{nm}^{(1)} = 0$ do các yếu tố ma trận không trung hòa trên đường chéo bằng 0.
    • Bổ chính bậc hai vào năng lượng và hàm sóng: $$\Delta E_{nm}^{(2)} = \sum_{j \neq n} \frac{(v_{nj} - E_{nm}^{(0)}r_{nj})(v_{jn} - E_{nm}^{(0)}r_{jn})}{(h_{jj}r_{nn} - h_{nn}r_{jj})r_{nn}}$$
    • Bổ chính bậc ba vào năng lượng (đóng góp phát triển mở rộng của tác giả): $$\Delta E_{nm}^{(3)} = \sum_{j \neq n} \sum_{k \neq n} \frac{(v_{nk} - E_{nm}^{(0)}r_{nk})(v_{kj} - E_{nm}^{(0)}r_{kj})(v_{jn} - E_{nm}^{(0)}r_{jn})}{(h_{kk}r_{nn} - h_{nn}r_{kk})(h_{jj}r_{nn} - h_{nn}r_{jj})r_{nn}} - \dots$$

Chương 3. Biểu thức giải tích cho năng lượng exciton trong đơn lớp TMDCs

Chương 3 áp dụng toàn bộ kết quả toán học vào việc xây dựng biểu thức giải tích trạng thái cơ bản 1s ($n=0, m=0$):

  • Bổ chính bậc không ($E_0^{(0)}$): Tính $E_0^{(0)} = -\frac{\omega^2}{8} - \frac{\omega}{\kappa \alpha} J_1(x)$. Xác định điểm cực trị từ điều kiện $\frac{\partial E_0^{(0)}}{\partial x} = 0$. Khai triển Taylor quanh $x = 0.75$ cho ra mối liên hệ $x = 0.75 - \frac{0.62452\alpha}{\kappa}$, dẫn đến biểu thức giải tích bậc không.
  • Bổ chính bậc hai ($\Delta E_0^{(2)}$): Khảo sát 7 số hạng ($j = 1$ đến $j = 7$). Số hạng đầu tiên ($j=1$) có giá trị rất bé ($\sim 10^{-7}\text{ meV}$) nên được bỏ qua. Tính toán chi tiết và khai triển Taylor tại $x = 0.75$ cho các số hạng từ $j=2$ đến $j=7$, thu được biểu thức bổ chính bậc hai.
  • Bổ chính bậc ba ($\Delta E_0^{(3)}$): Khảo sát 6 số hạng tương tác cặp $(j, k) \in {(1,2), (2,3), (3,4), (4,5), (5,6), (6,7)}$. Bỏ qua số hạng đầu $(1,2)$ do độ lớn không đáng kể ($\sim 10^{-6}\text{ meV}$), tính toán và khai triển giải tích 5 số hạng còn lại, hợp nhất thành biểu thức bổ chính bậc ba.
  • Tổng hợp biểu thức năng lượng toàn phần: Kết hợp $E_{1s} = E_0^{(0)} + \Delta E_0^{(2)} + \Delta E_0^{(3)}$ tạo thành biểu thức giải tích tường minh phụ thuộc trực tiếp vào hằng số điện môi $\kappa$ và tham số độ dài chắn $\alpha$.

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

Kết quả tính toán giải tích và đối sánh số liệu

Khóa luận đã xác định giá trị tham số cực trị $x$ thỏa mãn điều kiện $\frac{\partial E_0^{(0)}}{\partial x} = 0$ cho 4 chất bán dẫn đơn lớp TMDCs:

Thông số $WSe_2$ $WS_2$ $MoSe_2$ $MoS_2$
Giá trị $x$ tối ưu 0.75008 Dao động quanh 0.75 Dao động quanh 0.75 Dao động quanh 0.75

Đóng góp năng lượng chi tiết của từng số hạng trong bổ chính bậc hai ($\Delta E_0^{(2)}$) và bậc ba ($\Delta E_0^{(3)}$) đối với đơn lớp $WSe_2$:

Bậc bổ chính Chỉ số thành phần Giá trị năng lượng ($\text{meV}$)
Bổ chính bậc hai ($\Delta E_0^{(2)}$) Số hạng 1 ($j=1$) $-9.872 \times 10^{-7}$ (bỏ qua)
Số hạng 2 ($j=2$) $-2.013553$
Số hạng 3 ($j=3$) $-0.700910$
Số hạng 4 ($j=4$) $-0.218572$
Số hạng 5 ($j=5$) $-0.082738$
Số hạng 6 ($j=6$) $-0.036031$
Số hạng 7 ($j=7$) $-0.017455$
Bổ chính bậc ba ($\Delta E_0^{(3)}$) Số hạng 1 ($j=1, k=2$) $+0.000008$ (bỏ qua)
Số hạng 2 ($j=2, k=3$) $-0.256422$
Số hạng 3 ($j=3, k=4$) $-0.070159$
Số hạng 4 ($j=4, k=5$) $-0.024192$
Số hạng 5 ($j=5, k=6$) $-0.009820$
Số hạng 6 ($j=6, k=7$) $-0.004503$

Tổng hợp năng lượng exciton trạng thái 1s qua từng bậc gần đúng và so sánh với kết quả giải số chính xác từ công trình tham chiếu của Duy-Nhat Ly và cộng sự (2021):

Mức độ gần đúng giải tích Năng lượng $E_{1s}$ trên $WSe_2$ ($\text{meV}$) Sai số so với giải số
Gần đúng bậc không ($E_0^{(0)}$) $-164.21$ Sai số nhỏ
Gần đúng đến bậc hai ($E_0^{(0)} + \Delta E_0^{(2)}$) $-168.44$ Giảm đáng kể
Gần đúng đến bậc ba ($E_0^{(0)} + \Delta E_0^{(2)} + \Delta E_0^{(3)}$) $-168.617$ $< 1%$

Đóng góp khoa học của đề tài

  1. Xây dựng thành công biểu thức giải tích năng lượng exciton 1s trực tiếp từ thế Keldysh nguyên bản mà không cần dùng các dạng thế gần đúng hay thế Kratzer như các nghiên cứu trước.
  2. Bổ sung và mở rộng thành công lý thuyết nhiễu loạn có điều tiết (MPM) lên đến bổ chính bậc ba cho bài toán exciton đơn lớp 2D (các công trình trước chỉ dừng ở bậc hai).
  3. Đưa ra công cụ giải tích thuận tiện, cho phép xác định nhanh năng lượng exciton và hỗ trợ quy trình trích xuất thông số vật liệu TMDCs dựa trên độ nhạy của phổ exciton với hằng số điện môi và độ dài chắn.

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

Hạn chế của đề tài

  • Biểu thức giải tích cuối cùng thu được còn tương đối dài và đi kèm nhiều hệ số số trị phức tạp do việc khai triển Taylor nhiều số hạng bậc cao.
  • Phạm vi khảo sát giải tích mới dừng lại ở trạng thái cơ bản 1s ($n=0, m=0$), chưa thiết lập cho các trạng thái kích thích.
  • Mô hình chưa tích hợp các yếu tố tương tác ngoài như từ trường đều, hiệu ứng nhiệt độ lưới và hiệu ứng spin - quỹ đạo của exciton.

Hướng nghiên cứu phát triển

  • Mở rộng phương pháp để thiết lập biểu thức giải tích cho các trạng thái kích thích cao hơn ($2s, 3s, 4s, 5s$) trong cùng thế tương tác Keldysh.
  • Nghiên cứu các kỹ thuật biến đổi giải tích nhằm rút gọn cấu trúc biểu thức và tìm kiếm dạng biểu thức có tính phổ quát cao hơn.
  • Bổ sung vào Hamiltonian các số hạng mô tả từ trường ngoài và hiệu ứng nhiệt độ, từ đó tính toán bổ sung các yếu tố ma trận tương ứng.
  • Mở rộng phạm vi áp dụng biểu thức giải tích cho nhiều cấu trúc dị thể và vật liệu 2D đơn lớp mới.

Giá trị tham khảo

  • Đối tượng tham khảo: Sinh viên, học viên cao học và nghiên cứu sinh chuyên ngành Sư phạm Vật lý, Vật lý lý thuyết, Vật lý chất rắn, Khoa học vật liệu và Vật lý tính toán.
  • Giá trị nội dung:
    • Quy trình toán học lượng tử chuẩn mực về việc ứng dụng phép biến đổi Levi-Civita và phương pháp toán tử FK để chuyển đổi và giải bài toán thế phi điều hòa hai chiều.
    • Kỹ thuật tính toán yếu tố ma trận giải tích thông qua tích phân biến đổi Laplace cho thế Keldysh nguyên bản.
    • Khung thuật toán khai triển lý thuyết nhiễu loạn Rayleigh-Schrödinger có điều tiết đến bậc ba, có thể sử dụng làm tài liệu tham khảo cho các bài toán lượng tử không có nghiệm giải tích chính xác.

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

1. Tại sao khóa luận sử dụng thế Keldysh thay vì thế Coulomb thuần túy để mô tả tương tác electron - lỗ trống?

Trong vật liệu hai chiều như đơn lớp TMDCs, do hiệu ứng giảm chiều không gian và sự suy giảm phân cực điện môi ngoài mặt phẳng, tương tác tĩnh điện giữa electron và lỗ trống bị chắn lại theo cơ chế phi Coulomb. Thế Keldysh kết hợp hàm Struve $H_0$ và hàm Bessel $Y_0$ mô tả chính xác hiệu ứng chắn này ở khoảng cách ngắn và khoảng cách trung bình, phù hợp với các dữ liệu quang phổ thực nghiệm.

2. Vai trò của tham số tự do $\omega$ trong phương pháp toán tử FK là gì?

Tham số thực dương $\omega$ được đưa vào khi định nghĩa hệ toán tử sinh - hủy ($a, a^\dagger, b, b^\dagger$) đóng vai trò như một tham số điều tiết biến phân. Việc điều chỉnh $\omega$ giúp tối ưu hóa hệ hàm cơ sở, đẩy nhanh tốc độ hội tụ của chuỗi nhiễu loạn mà không làm thay đổi bản chất nghiệm vật lý của phương trình Schrödinger ban đầu.

3. Tại sao tác giả chọn khai triển chuỗi Taylor quanh điểm $x = 0.75$?

Khi giải điều kiện cực trị năng lượng bậc không $\frac{\partial E_0^{(0)}}{\partial x} = 0$ cho các vật liệu TMDCs cụ thể ($WSe_2, WS_2, MoSe_2, MoS_2$), giá trị tính toán của tham số không thứ nguyên $x = \frac{\alpha}{\omega}$ đều dao động rất sát điểm $0.75$ (cụ thể với $WSe_2$ là $0.75008$). Do đó, chọn điểm mốc $x = 0.75$ để khai triển Taylor giúp chuỗi hội tụ nhanh và tối thiểu hóa sai số số trị.

4. Điểm mở rộng lý thuyết của khóa luận so với các công trình trước là gì?

Khóa luận đã thiết lập công thức tường minh cho bổ chính bậc ba vào năng lượng ($\Delta E^{(3)}$) và hàm sóng trong khuôn khổ lý thuyết nhiễu loạn Rayleigh-Schrödinger có điều tiết cho exciton TMDCs, vượt qua giới hạn của các nghiên cứu trước đó (vốn chỉ dừng lại ở bổ chính bậc hai).

5. Độ chính xác của biểu thức giải tích năng lượng exciton trạng thái 1s đạt mức nào?

Khi tính toán đầy đủ đến bổ chính bậc ba cho đơn lớp $WSe_2$, giá trị năng lượng giải tích thu được là $-168.617\text{ meV}$, đạt độ sai lệch dưới $1%$ so với kết quả giải số chính xác cao trong công trình của Duy-Nhat Ly và cộng sự (2021).

Kết luận

Khóa luận tốt nghiệp của tác giả Nguyễn Ngọc Huy đã giải quyết thành công bài toán xây dựng biểu thức giải tích tường minh cho năng lượng exciton trạng thái 1s trong đơn lớp TMDCs từ dạng nguyên bản của thế tương tác Keldysh. Bằng việc kết hợp hiệu quả phép biến đổi Levi-Civita, phương pháp đại số toán tử FK và mở rộng lý thuyết nhiễu loạn Rayleigh-Schrödinger lên đến bổ chính bậc ba, nghiên cứu đã đạt được kết quả năng lượng có độ chính xác cao với sai số dưới $1%$ so với phương pháp giải số. Công trình cung cấp một công cụ giải tích hữu ích cho việc tính toán nhanh mức năng lượng exciton và phục vụ trích xuất thông số cấu trúc cho các vật liệu bán dẫn hai chiều họ TMDCs.