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

Trong cơ học lượng tử, phương trình Schrödinger mô tả phần lớn các hệ lượng tử thực tế không thể giải được một cách chính xác dưới dạng giải tích, ngoại trừ một số ít hệ đơn giản như nguyên tử hydro ba chiều, dao động tử điều hòa hoặc giếng thế vuông góc. Để giải quyết các bài toán phức tạp, các phương pháp gần đúng đã được phát triển, trong đó lý thuyết nhiễu loạn đóng vai trò trung tâm trong vật lý và hóa học lý thuyết. Tuy nhiên, đối với các hệ suy biến (nhiều trạng thái lượng tử khác nhau có cùng mức năng lượng), lý thuyết nhiễu loạn chính tắc đòi hỏi các công thức tính bổ chính năng lượng và hàm sóng rất phức tạp, phải giải các định thức thế năng bậc cao, gây khó khăn lớn cho việc tính toán giải tích lẫn lập trình số trị.

Sự phát triển của vật lý chất rắn và công nghệ vật liệu từ sau mốc tách chiết thành công graphene đơn lớp (năm 2004) đã thúc đẩy mạnh mẽ các nghiên cứu về hệ vật liệu hai chiều (2D). Trong nhóm vật liệu này, các đơn lớp bán dẫn thuộc nhóm dichalcogenide kim loại chuyển tiếp (Transition-Metal Dichalcogenides - TMD) như $\text{MoS}_2$, $\text{MoSe}_2$, $\text{WS}_2$, $\text{WSe}_2$ thu hút sự chú ý đặc biệt do sở hữu vùng cấm trực tiếp và các đặc tính quang - điện tử vượt trội so với dạng khối. Trong đơn lớp TMD, sự hình thành exciton (trạng thái liên kết giữa một electron và một lỗ trống) là cơ chế chuyển dời quang học chủ đạo, quyết định trực tiếp tính chất quang học của vật liệu.

Do sự suy giảm số chiều và hiệu ứng chắn điện môi không đồng nhất từ môi trường xung quanh, tương tác tĩnh điện giữa electron và lỗ trống trong đơn lớp TMD không còn tuân theo thế Coulomb thuần túy mà được mô tả chính xác bởi thế Keldysh. Khi đặt vật liệu trong từ trường ngoài, phổ năng lượng exciton chịu sự biến đổi phức tạp. Nhiều nghiên cứu lý thuyết trước đây (như của tác giả Ngọc Trâm, Duy Anh) hoặc chỉ áp dụng thế Coulomb, hoặc chỉ khảo sát ở nhiệt độ $0\text{ K}$ và bỏ qua ảnh hưởng của nhiệt độ. Trong khi đó, các phép đo thực nghiệm quang học hiện đại thường được tiến hành ở cả nhiệt độ phòng ($300\text{ K}$). Sự chuyển động nhiệt của khối tâm exciton trong từ trường làm xuất hiện hiệu ứng Stark chuyển động nhiệt, tạo ra một thế năng cảm ứng phụ thuộc vào cả nhiệt độ $T$ và từ trường $B$. Do đó, việc xây dựng một phương pháp nhiễu loạn mới có khả năng giải quyết bài toán nguyên tử hai chiều suy biến và tính toán phổ năng lượng exciton trong TMD dưới tác động đồng thời của từ trường và nhiệt độ là yêu cầu cấp thiết.

Mục tiêu nghiên cứu:

  1. Xây dựng lý thuyết nhiễu loạn có điều tiết (Regulated Perturbation Theory) dựa trên phương pháp toán tử Feranchuk - Komarov để giải bài toán nguyên tử hai chiều có suy biến.
  2. Áp dụng lý thuyết nhiễu loạn có điều tiết để tính toán phổ năng lượng exciton trong đơn lớp TMD ($\text{WSe}_2$) đặt trong từ trường ngoài, có xét đến ảnh hưởng của yếu tố nhiệt độ.

Nhiệm vụ nghiên cứu cụ thể:

  • Nhiệm vụ 1: Tìm hiểu tổng quan về bài toán nguyên tử hai chiều có suy biến, exciton trong đơn lớp TMD, bài toán exciton trong từ trường và hạn chế của lý thuyết nhiễu loạn suy biến chính tắc.
  • Nhiệm vụ 2: Nghiên cứu phương pháp toán tử Feranchuk - Komarov (FK) và phát triển lý thuyết nhiễu loạn có điều tiết, thiết lập các biểu thức giải tích tính bổ chính năng lượng và hàm sóng đến bậc hai.
  • Nhiệm vụ 3: Xây dựng các biểu thức yếu tố ma trận của toán tử Hamiltonian cho bài toán exciton đơn lớp TMD trong từ trường và trường nhiệt độ.
  • Nhiệm vụ 4: Lập trình tính toán số trị năng lượng exciton ở các trạng thái lượng tử $1s, 2s, 2p_-, 2p_+, 3s$ trong đơn lớp $\text{WSe}_2$; khảo sát vai trò và tiêu chí chọn tham số điều tiết tự do $\omega$.

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

  • Đối tượng nghiên cứu: Exciton trung hòa hai chiều (hệ liên kết gồm 1 electron mang điện tích $-e$ và 1 lỗ trống mang điện tích $+e$) trong đơn lớp vật liệu bán dẫn TMD.
  • Phạm vi nghiên cứu: Khảo sát đơn lớp $\text{WSe}_2$ đặt trong từ trường đều vuông góc với mặt phẳng cấu trúc có độ lớn từ $0\text{ Tesla}$ đến $100\text{ Tesla}$; nhiệt độ khảo sát tại $T = 0\text{ K}$ và $T = 300\text{ K}$; các mức trạng thái kích thích thấp ($1s, 2s, 2p, 3s$); tính toán bổ chính lý thuyết nhiễu loạn dừng lại ở bậc hai.

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

Khóa luận kết hợp giữa mô hình hóa giải tích cơ học lượng tử và tính toán số trị bằng lập trình. Khung lý thuyết và công cụ toán học được triển khai qua các nội dung sau:

Khung lý thuyết và mô hình vật lý

  • Bài toán nguyên tử hai chiều và tương tác trường ngoài: Phương trình Schrödinger cho chuyển động tương đối của exciton trong mặt phẳng $(x, y)$ dưới từ trường đều vuông góc phương $z$ ($\mathbf{B} = (0, 0, B)$). Hamiltonian chứa thành phần Zeeman mô tả tương tác moment từ quỹ đạo $H_{\text{Zeeman}} = -\frac{e\hbar}{2\mu} l_z B$ và thành phần nghịch từ tỉ lệ với bình phương từ trường $H_{\text{diamag}} = \frac{e^2 B^2 (x^2 + y^2)}{8\mu}$.
  • Thế tương tác Keldysh: Mô tả tương tác chắn giữa electron và lỗ trống trong màng mỏng hai chiều kẹp giữa hai môi trường điện môi: $$V_{\text{Keldysh}}(r) = -\frac{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]$$ Trong đó, $H_0$ là hàm Struve bậc 0, $Y_0$ là hàm Bessel loại hai bậc 0, $r_0 = 2\pi\chi_{2D}$ là chiều dài chắn đặc trưng cho độ phân cực của đơn lớp, và $\kappa = (\varepsilon_1 + \varepsilon_2)/2$ là hằng số điện môi hiệu dụng trung bình của môi trường bao quanh. Để phục vụ tính toán ma trận, thế Keldysh được biểu diễn qua phép biến đổi tích phân Laplace.
  • Hiệu ứng Stark chuyển động nhiệt: Khi tách chuyển động khối tâm bằng toán tử giả động lượng bảo toàn $\mathbf{P}c = \mathbf{P} - \frac{e}{2}\mathbf{B}\times\mathbf{r}$, vận tốc chuyển động nhiệt của khối tâm exciton $|\mathbf{K}| = \sqrt{2M k_B T}$ trong từ trường ngoài sinh ra thế cảm ứng Stark: $$V{ms} = -\frac{\sqrt{2 k_B T}}{M} e B x$$ Thế này chỉ xuất hiện khi nhiệt độ $T > 0\text{ K}$ và phụ thuộc tuyến tính vào từ trường $B$.

Phương pháp toán học và lý thuyết nhiễu loạn

  • Phép biến đổi tọa độ Levi-Civita: Chuyển đổi không gian tọa độ hai chiều thực $(x, y)$ sang không gian biến số $(u, v)$ theo hệ thức $x = u^2 - v^2$, $y = 2uv$, giúp chuyển phương trình vi phân dạng nguyên tử về dạng phương trình dao động tử điều hòa hai chiều.
  • Phương pháp toán tử Feranchuk - Komarov (FK): Định nghĩa các toán tử sinh ($\hat{a}^\dagger, \hat{b}^\dagger$) và toán tử hủy ($\hat{a}, \hat{b}$) phụ thuộc tham số tự do $\omega$. Xây dựng đại số kín gồm bộ ba toán tử rút gọn: $$\hat{M} = \hat{a}\hat{b}, \quad \hat{M}^\dagger = \hat{a}^\dagger\hat{b}^\dagger, \quad \hat{N} = \hat{a}^\dagger\hat{a} + \hat{b}^\dagger\hat{b} + 1$$ Thỏa mãn các hệ thức giao hoán: $[\hat{M}, \hat{N}] = 2\hat{M}$, $[\hat{N}, \hat{M}^\dagger] = 2\hat{M}^\dagger$, $[\hat{M}, \hat{M}^\dagger] = \hat{N}$.
  • Lý thuyết nhiễu loạn có điều tiết: Phân tách triệt để toán tử Hamiltonian $\hat{H}$ và toán tử chuẩn $\hat{R}$ thành phần trung hòa ($\hat{H}_0, \hat{R}_0$) và phần không trung hòa ($\hat{V}, \hat{R}_V$). Việc phân tách này giúp loại bỏ sự suy biến ngay từ gần đúng bậc không, cho phép áp dụng sơ đồ nhiễu loạn không suy biến để tính toán bổ chính mà không gặp hiện tượng kỳ dị ở mẫu số.

Cơ sở dữ liệu và phương pháp tính toán

  • Hệ hàm cơ sở: Sử dụng hệ hàm trực chuẩn của dao động tử điều hòa hai chiều $|n, m\rangle$ để tính toán các yếu tố ma trận đại số.
  • Phương pháp số: Xây dựng mã nguồn tính toán bằng ngôn ngữ Python để giải phương trình xác định tham số $\omega$, tính tích phân yếu tố ma trận thế Keldysh và lấy tổng chuỗi bổ chính bậc hai với ngưỡng cắt cụt $j_{\text{max}} = 100$.
  • Dữ liệu kiểm chứng: Sử dụng các thông số thực nghiệm và số liệu đối sánh từ công trình nghiên cứu của D. Le (Phys. Rev. B, 2019) cho đơn lớp $\text{WSe}_2$: chiều dài chắn $r_0 = 4.51\text{ nm}$, hằng số điện môi $\kappa = 4.5$, khối lượng hiệu dụng của electron $m_e = 0.36 m_0$, khối lượng hiệu dụng của lỗ trống $m_h = 0.45 m_0$.

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

Chương 1: Bài toán nguyên tử hai chiều suy biến. Exciton đơn lớp TMD

Chương 1 thiết lập mô hình vật lý và toán học của bài toán exciton hai chiều trong vật liệu TMD đặt trong trường ngoài:

  • Nguyên tử hydro 2D và trường ngoài: Trình bày phương trình Schrödinger của bài toán nguyên tử hai chiều không trường ngoài và khi có mặt từ trường $B$ cùng điện trường $E$. Năng lượng của nguyên tử hydro hai chiều bị suy biến theo số lượng tử từ $m$.
  • Tổng quan về exciton và đơn lớp TMD: Phân biệt hai mô hình exciton chính: Frenkel exciton (bán kính nhỏ, năng lượng liên kết lớn $\sim 1.5\text{ eV}$, phổ biến trong chất cách điện) và Mott - Wannier exciton (bán kính tương tác lớn, năng lượng liên kết nhỏ, phổ biến trong chất bán dẫn). Khóa luận tập trung vào mô hình exciton trung hòa (1 electron và 1 lỗ trống). Cấu trúc mạng tinh thể đơn lớp TMD có dạng lục giác $MX_2$ gồm một lớp kim loại chuyển tiếp kẹp giữa hai lớp chalcogen với độ dày cỡ $7\text{ \AA}$ (đối với $\text{MoS}_2$, khoảng cách giữa hai lớp là $6.5\text{ \AA}$).
  • Phương trình vi phân cho exciton TMD trong từ trường và nhiệt độ: Thiết lập phương trình chuyển động tương đối trong hệ tọa độ khối tâm. Sau khi áp dụng hệ đơn vị Hartree hiệu dụng ($2Ry^* = \mu e^4 / 16\pi^2\varepsilon_0^2\hbar^2$, bán kính Bohr hiệu dụng $a_B^* = 4\pi\varepsilon_0\hbar^2 / \mu e^2$) và phép biến đổi Levi-Civita $(u, v)$, phương trình Schrödinger được chuyển về dạng phương trình vi phân không thứ nguyên chứa tham số nhiệt độ $\Theta = k_B T / 2Ry^*$ và cường độ từ trường không thứ nguyên $\gamma$.
  • Biểu diễn toán tử: Sử dụng phương pháp toán tử FK để biểu diễn các toán tử tọa độ, vi phân, thế Keldysh và thế cảm ứng Stark chuyển động qua các toán tử đại số $\hat{M}, \hat{M}^\dagger, \hat{N}$, tạo tiền đề trực tiếp cho việc xây dựng lý thuyết nhiễu loạn ở chương tiếp theo.

Chương 2: Lý thuyết nhiễu loạn giải bài toán hai chiều suy biến

Chương 2 tập trung phát triển bộ máy toán học của lý thuyết nhiễu loạn:

  • Lý thuyết nhiễu loạn có suy biến chính tắc (Rayleigh - Schrödinger): Phân tích quy trình giải phương trình bậc $g_N$ từ định thức thế năng ma trận để tìm bổ chính bậc một $E^{(1)}$ và chỉ ra sự phức tạp khi mở rộng lên bậc hai $E^{(2)}$. Với thế Keldysh phi Coulomb, phương trình phi nhiễu loạn không thể giải chính xác giải tích, dẫn đến việc áp dụng lý thuyết chính tắc gặp trở ngại lớn về mặt tính toán.
  • Xây dựng lý thuyết nhiễu loạn có điều tiết:
    • Phương trình trị riêng bậc không có dạng: $\hat{H}0 \Phi{nm} = E_{nm}^{(0)} \hat{R}0 \Phi{nm}$, cho năng lượng bậc không: $$E_{nm}^{(0)} = \frac{h_{nm}}{r_{nm}}$$
    • Do tính chất phân tách triệt để của các toán tử không trung hòa $\langle \Phi_{nm} | \hat{V} | \Phi_{nm} \rangle = 0$ và $\langle \Phi_{nm} | \hat{R}V | \Phi{nm} \rangle = 0$, bổ chính năng lượng bậc một triệt tiêu hoàn toàn: $$\Delta E_{nm}^{(1)} = 0$$
    • Thiết lập biểu thức giải tích tường minh cho bổ chính hàm sóng bậc một $\Delta \Psi_{nm}^{(1)}$ và bổ chính năng lượng bậc hai $\Delta E_{nm}^{(2)}$: $$\Delta E_{nm}^{(2)} = \sum_{j=0}^{+\infty} \sum_{\substack{m'=-j \ (j,m') \neq (n,m)}}^{+j} \frac{\left( V_{jnmm'} r_{nm} - h_{nm} T_{jnmm'} \right) \left( V_{njmm'} r_{nm} - h_{nm} T_{njmm'} \right)}{\left( h_{jm'} r_{nm} - h_{nm} r_{jm'} \right) r_{nm}^2}$$
    • Nhờ tính không suy biến của phổ trị riêng bậc không trong không gian toán tử trung hòa, mẫu số của biểu thức bổ chính bậc hai luôn khác không mà không cần giải phương trình ma trận suy biến.
  • Tính toán yếu tố ma trận: Xác định các yếu tố ma trận $h_{nm}, r_{nm}, V_{jk}^{nm}, T_{jk}^{nm}$ trên hệ hàm cơ sở dao động tử điều hòa hai chiều $|n, m\rangle$. Trong đó, ma trận thế Keldysh $V_{jk}$ được tính toán giải tích kết hợp tích phân số qua hàm trung gian $U_{s, m}$.

Chương 3: Áp dụng lý thuyết nhiễu loạn có điều tiết tính năng lượng exciton trong đơn lớp TMD

Chương 3 áp dụng lý thuyết đã xây dựng để tính toán các mức năng lượng của exciton trong đơn lớp $\text{WSe}_2$:

  • Xác định tham số tự do $\omega$: Tham số $\omega$ đóng vai trò điều tiết quy mô của phần phi nhiễu loạn và phần nhiễu loạn. Nghiên cứu đề xuất điều kiện dừng của năng lượng bậc không: $$\frac{\partial E_{nm}^{(0)}(\omega)}{\partial \omega} = 0$$ Dẫn đến phương trình xác định tham số tối ưu $\omega_0$. Với mỗi trạng thái lượng tử và mỗi giá trị từ trường $B$, phương trình này cung cấp một giá trị $\omega_0$ riêng biệt.
  • Tính toán năng lượng exciton đến bậc hai: Phân tách tổng bổ chính bậc hai thành hai phần: thành phần không phụ thuộc nhiệt độ $\Delta E^{(2)}(B, T = 0\text{ K})$ ứng với $m' = m$ và thành phần bổ chính nhiệt độ $\Delta E^{(2)}(B, T)$ ứng với $m' \neq m$ (xuất phát từ thế $V_{ms}$).
  • Kết quả tính số và đánh giá sai số: Năng lượng toàn phần $E_{nm} = E_{nm}^{(0)} + \Delta E_{nm}^{(2)}$ được tính cho các mức $1s, 2s, 2p_-, 2p_+, 3s$ trong dải từ trường $0 - 100\text{ Tesla}$. Kết quả được đối sánh chi tiết với dữ liệu số chính xác từ tài liệu tham khảo [29].
  • Khảo sát tham số điều tiết tối ưu $\omega_2$: Tác giả chứng minh rằng nếu chọn tham số $\omega_2$ tối ưu trực tiếp cho bổ chính bậc hai (thay vì $\omega_0$ từ bậc không), độ chính xác của năng lượng exciton được cải thiện vượt bậc.

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

Kết quả tính toán chính

Kết quả tính toán phổ năng lượng exciton trong đơn lớp $\text{WSe}_2$ bằng lý thuyết nhiễu loạn có điều tiết đến bậc hai được tổng hợp trong bảng đối sánh sau:

Trạng thái lượng tử Từ trường $B$ (Tesla) Năng lượng $E(B, T=0\text{ K})$ (eV) Sai số tỉ đối $\delta_{E(B, T=0)}$ so với [29] Bổ chính nhiệt $T=300\text{ K}$ Đánh giá độ chính xác
$1s$ ($n=0, m=0$) $0\text{ T}$ $-0.1687$ $< 0.2%$ Có tính toán Độ chính xác cao
$1s$ ($n=0, m=0$) $10\text{ T}$ $-0.1685$ $\sim 0.1%$ $\Delta E(B, T)$ Chính xác đến 4 chữ số thập phân khi dùng $\omega_2$
$1s$ ($n=0, m=0$) $100\text{ T}$ $-0.1512$ $< 1.0%$ $\Delta E(B, T)$ Duy trì sai số dưới $1%$
$2p_-$ ($n=1, m=-1$) $10\text{ T}$ Phù hợp giải tích $0.00016%$ $\Delta E(B, T)$ Chính xác đến 6 chữ số thập phân
$2p_+$ ($n=1, m=1$) $10\text{ T}$ Phù hợp giải tích $-0.52382%$ Sai số nhiệt lớn ($4896.5%$) Sai số năng lượng $T=0\text{ K}$ nhỏ
$2s$ ($n=1, m=0$) $\le 10\text{ T}$ Phù hợp nghiệm số $< 5.0%$ $\Delta E(B, T)$ Sai số tăng khi $B \to 100\text{ T}$
$3s$ ($n=2, m=0$) $\le 10\text{ T}$ Phù hợp nghiệm số $< 5.0%$ $\Delta E(B, T)$ Sai số tăng khi $B \to 100\text{ T}$

Chi tiết về vai trò của việc chọn tham số $\omega$ cho trạng thái $1s$ tại $B = 10\text{ Tesla}$:

  • Nghiệm số chính xác từ công trình [29]: $E_{\text{exact}} = -0.168523698\text{ eV}$.
  • Sử dụng tham số $\omega_0$ (xác định từ cực trị năng lượng bậc không): thu được kết quả đúng 2 chữ số phần thập phân.
  • Sử dụng tham số điều tiết tối ưu $\omega_2 = 0.37581305$: năng lượng tính đến bậc hai đạt giá trị chính xác đến 4 chữ số phần thập phân.

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

  1. Phương pháp luận: Ứng dụng thành công lý thuyết nhiễu loạn có điều tiết kết hợp đại số toán tử FK để giải quyết bài toán nguyên tử hai chiều suy biến; loại bỏ hoàn toàn yêu cầu giải định thức thế năng phức tạp của lý thuyết suy biến chính tắc.
  2. Tính toán giải tích và số trị: Thiết lập biểu thức giải tích tổng quát cho bổ chính năng lượng và hàm sóng đến bậc hai đối với exciton hai chiều chịu thế Keldysh trong từ trường và hiệu ứng nhiệt độ.
  3. Ý nghĩa ứng dụng: Chứng minh khả năng xây dựng các biểu thức giải tích gần đúng có độ chính xác cao cho phổ exciton TMD ở vùng từ trường thấp ($B \le 10\text{ Tesla}$), làm cơ sở lý thuyết phục vụ việc trích xuất các thông số cấu trúc vật liệu (bán kính chắn, khối lượng hiệu dụng) từ thực nghiệm quang phổ.

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

Hạn chế

  • Chuỗi tổng vô hạn trong biểu thức bổ chính năng lượng bậc hai phải thực hiện cắt cụt ở $j_{\text{max}} = 100$ trong quá trình tính toán số trị.
  • Nghiên cứu mới dừng lại ở việc thiết lập và tính toán bổ chính đến bậc hai ($\Delta E^{(2)}$), chưa tính toán các bậc bổ chính cao hơn ($E^{(3)}, E^{(4)}$).
  • Bổ chính nhiệt độ xuất phát từ hiệu ứng Stark chuyển động ($V_{ms}$) có sai số tỉ đối rất lớn khi so sánh với công trình đối chứng (sai số thấp nhất là $34.3%$ ở mức $1s$, $B = 0\text{ T}$ và tăng rất cao ở các trạng thái kích thích khác).
  • Ở vùng từ trường rất cao ($100\text{ Tesla}$) đối với các trạng thái kích thích cao ($2s, 3s$), sai số tính toán còn tương đối lớn so với vùng từ trường thấp.

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

  • Tiếp tục mở rộng xây dựng lý thuyết nhiễu loạn có điều tiết cho các bậc bổ chính cao hơn và khảo sát sự hội tụ của chuỗi năng lượng.
  • Khảo sát chuyên sâu cơ chế ảnh hưởng của tham số tự do $\omega$ lên tốc độ hội tụ của các bậc bổ chính.
  • Kiểm tra, hiệu chỉnh lại phương pháp tính và các yếu tố ma trận của thành phần bổ chính phụ thuộc nhiệt độ ($V_{ms}$) nhằm giảm thiểu sai số.
  • Hoàn thiện việc xây dựng công thức giải tích gần đúng hoàn chỉnh cho năng lượng exciton TMD phụ thuộc từ trường, hỗ trợ nghiên cứu lý thuyết và phân tích thực nghiệm quang phổ.

Giá trị tham khảo

  • Đối tượng tham khảo: Sinh viên, học viên cao học, nghiên cứu sinh chuyên ngành Vật lý lý thuyết, Vật lý chất rắn, Sư phạm Vật lý và Khoa học vật liệu tính toán đang nghiên cứu về vật liệu hai chiều hoặc các phương pháp gần đúng trong cơ học lượng tử.
  • Nội dung giá trị:
    • Chương 1 và Phụ lục A, C, D: Cung cấp chi tiết các bước biến đổi vi phân và toán tử (tách chuyển động khối tâm, phép biến đổi tọa độ Levi-Civita, thiết lập đại số kín của toán tử sinh hủy FK hai chiều).
    • Chương 2: Cung cấp khung mẫu toán học hoàn chỉnh về lý thuyết nhiễu loạn có điều tiết, có thể tái sử dụng cho các bài toán lượng tử 2D khác có thế tương tác phi Coulomb hoặc suy biến.
    • Chương 3: Phương pháp lựa chọn và tối ưu hóa tham số điều tiết tự do $\omega$ để đạt độ chính xác số trị cao mà không cần tăng bậc nhiễu loạn.

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

1. Tại sao tương tác giữa electron và lỗ trống trong đơn lớp TMD không thể dùng thế Coulomb thuần túy mà phải sử dụng thế Keldysh?
Do cấu trúc đơn lớp TMD có độ dày rất nhỏ (cỡ $7\text{ \AA}$), khi số chiều bị suy giảm xuống hai chiều, các đường sức điện trường giữa electron và lỗ trống bị trải rộng ra môi trường điện môi bao quanh. Hiệu ứng chắn điện môi không đồng nhất này khiến tương tác tĩnh điện bị biến đổi, không còn tỷ lệ nghịch thuần túy với khoảng cách $1/r$ như thế Coulomb mà phải mô tả bằng thế Keldysh chứa hàm Struve $H_0$ và hàm Bessel $Y_0$.

2. Bổ chính năng lượng bậc một trong lý thuyết nhiễu loạn có điều tiết của đề tài có giá trị bằng bao nhiêu và vì sao?
Bổ chính năng lượng bậc một $\Delta E_{nm}^{(1)}$ luôn bằng 0. Nguyên nhân là do phương pháp toán tử FK đã phân tách Hamiltonian thành phần trung hòa $\hat{H}0, \hat{R}0$ và phần không trung hòa $\hat{V}, \hat{R}V$. Các yếu tố ma trận chéo của các toán tử không trung hòa trên hệ hàm cơ sở trực chuẩn đều đồng nhất bằng không ($\langle \Phi{nm} | \hat{V} | \Phi{nm} \rangle = 0$ và $\langle \Phi{nm} | \hat{R}V | \Phi{nm} \rangle = 0$).

3. Tham số tự do $\omega$ có vai trò gì và được xác định bằng cách nào?
Tham số tự do $\omega$ xuất hiện trong định nghĩa toán tử sinh hủy FK, có vai trò điều tiết quy mô giữa phần phi nhiễu loạn và phần nhiễu loạn mà không làm thay đổi nghiệm chính xác của phương trình Schrödinger ban đầu. Tham số này có thể được xác định bằng cách giải phương trình cực trị năng lượng bậc không $\partial E^{(0)}/\partial \omega = 0$ (cho ra $\omega_0$), hoặc lựa chọn giá trị tối ưu trực tiếp ở bậc hai ($\omega_2$) để cực tiểu hóa sai số so với nghiệm chính xác.

4. Tại sao yếu tố nhiệt độ lại làm xuất hiện thế cảm ứng Stark chuyển động trong phương trình Schrödinger?
Khi exciton chuyển động nhiệt trong mặt phẳng đơn lớp với vận tốc khối tâm $\mathbf{v}{cm}$ (liên hệ với động lượng khối tâm $|\mathbf{K}| = \sqrt{2M k_B T}$), từ trường ngoài $\mathbf{B}$ vuông góc với mặt phẳng sẽ tạo ra một điện trường cảm ứng Lorentz $\mathbf{E}{ind} = \mathbf{v}{cm} \times \mathbf{B}$ tác dụng lên exciton. Tương tác lưỡng cực điện của exciton trong điện trường cảm ứng này sinh ra thế năng Stark chuyển động nhiệt $V{ms}$, tỉ lệ thuận với $\sqrt{T}$ và $B$.

5. Độ chính xác của phổ năng lượng exciton tính toán trong khóa luận đạt mức độ nào ở nhiệt độ $0\text{ K}$?
Ở $T = 0\text{ K}$ và vùng từ trường dưới $10\text{ Tesla}$, năng lượng exciton ở các trạng thái $1s$ và $2p$ có sai số tỉ đối dưới $1%$ so với nghiệm chính xác; đặc biệt trạng thái $2p_-$ ở $10\text{ Tesla}$ đạt sai số chỉ $0.00016%$. Đối với trạng thái $1s$ tại $10\text{ Tesla}$, khi sử dụng tham số điều tiết tối ưu $\omega_2 = 0.37581305$, kết quả tính toán chính xác đến 4 chữ số thập phân so với giá trị đối chứng.

Kết luận

Khóa luận đã xây dựng thành công lý thuyết nhiễu loạn có điều tiết dựa trên phương pháp toán tử Feranchuk - Komarov để giải bài toán nguyên tử hai chiều suy biến, khắc phục được sự phức tạp của lý thuyết suy biến chính tắc. Phương pháp được áp dụng hiệu quả để tính toán phổ năng lượng của exciton trung hòa trong đơn lớp bán dẫn $\text{WSe}_2$ dưới tác động của từ trường ngoài đến $100\text{ Tesla}$ và hiệu ứng nhiệt độ. Kết quả tính toán ở nhiệt độ $0\text{ K}$ đạt độ chính xác cao so với nghiệm số chính xác, khẳng định tiềm năng ứng dụng của lý thuyết trong việc xây dựng các biểu thức giải tích phục vụ nghiên cứu quang học vật liệu hai chiều.