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

Vật liệu hai chiều (2D) sở hữu nhiều tính chất vật lý và hóa học đặc thù, trở thành tâm điểm nghiên cứu của khoa học vật liệu kể từ công bố năm 2004 của Geim và Novoselov về graphene. Mặc dù graphene có cấu trúc đơn lớp phẳng với độ linh động hạt tải cao ở nhiệt độ phòng, độ dẫn nhiệt vượt trội và hệ số truyền quang học lớn, nhược điểm năng lượng vùng cấm bằng không (hoạt động như kim loại) đã hạn chế đáng kể khả năng ứng dụng trong các linh kiện bán dẫn. Nhằm khắc phục hạn chế này, các chất bán dẫn hai chiều khác có cấu trúc tương tự đã được quan tâm nghiên cứu, tiêu biểu là vật liệu dichalcogenide kim loại chuyển tiếp (TMDs - Transition Metal Dichalcogenides) và boron-nitride lục giác ($h\text{-BN}$).

Vật liệu bán dẫn 2D TMDs có công thức chung là $\text{MX}_2$, trong đó $\text{M}$ là kim loại chuyển tiếp (như $\text{Ti, V, Nb, Mo, W}$) và $\text{X}$ là nguyên tử chalcogen ($\text{S, Se, Te}$). Cấu trúc đơn lớp TMDs gồm một lớp kim loại chuyển tiếp được kẹp giữa hai lớp chalcogen theo dạng lăng trụ tam giác. Các đơn lớp TMDs thuộc nhóm VI như $\text{MoS}_2, \text{MoSe}_2, \text{WS}_2, \text{WSe}_2$ sở hữu năng lượng vùng cấm trực tiếp nằm trong khoảng từ phổ khả kiến đến gần hồng ngoại, mở ra tiềm năng lớn trong chế tạo tế bào quang điện, diode phát quang, nguồn photon đơn, laser exciton và transistor quang điện tử.

Trong các chất bán dẫn 2D, chuyển dịch quang học chủ yếu gắn liền với sự hình thành exciton – một chuẩn hạt tạo thành từ tương tác Coulomb giữa electron mang điện tích âm và lỗ trống mang điện tích dương, có cấu trúc tương tự nguyên tử hydro (exciton Mott-Wannier). Khi số chiều của hệ giảm từ 3D xuống 2D, hiệu ứng giam giữ không gian làm tương tác Coulomb giữa electron và lỗ trống tăng lên đáng kể, khiến cấu trúc phổ exciton trở nên sắc nét hơn. Tuy nhiên, năng lượng của exciton ở các trạng thái kích thích cao rất khó đo đạc bằng thực nghiệm quang học thông thường. Đặt điện trường ngoài song song với bề mặt vật liệu là giải pháp quan trọng nhằm điều chỉnh tính chất quang học, hỗ trợ quá trình ion hóa exciton (khi ion hóa nhiệt không hiệu quả trong vật liệu có năng lượng liên kết lớn như TMDs) và cho phép quan sát hiệu ứng dịch chuyển, tách mức năng lượng Stark.

Việc giải phương trình Schrödinger cho exciton 2D trong điện trường đều từng được thực hiện qua các phương pháp khác nhau:

  • A. Gelten (1974) sử dụng hệ tọa độ parabol không thứ nguyên, tách biến phương trình Schrödinger thành hai phương trình một chiều và tính trị riêng trạng thái liên kết bằng phương pháp gần đúng WKB.
  • F. L. Lederman và J. D. Dow (1976) chuyển sang hệ tọa độ parabol theo định nghĩa khác, kết hợp công thức Elliott về hệ số hấp thụ để giải số chính xác.
  • Pokutnyi et al. dùng hệ tọa độ parabol, áp dụng phương pháp ma trận Numerov giải phương trình một chiều kết hợp gần đúng WKB để tính hệ số xuyên ngầm.

Khóa luận kế thừa và phát triển phương pháp toán tử FK (FK-OM - Feranchuk-Komarov Operator Method, đề xuất năm 1982 tại Đại học Belarus), vốn đã thành công trên các hệ nguyên tử hydro và exciton 2D trong từ trường đều, để giải phương trình Schrödinger cho exciton 2D đặt trong điện trường đều.

Mục tiêu nghiên cứu: Phát triển phương pháp toán tử FK cho bài toán exciton hai chiều (2D) trong điện trường đều nhằm xác định nghiệm số chính xác của phổ năng lượng.

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

  1. Tìm hiểu tổng quan về vật liệu bán dẫn hai chiều, exciton và các phương pháp giải phương trình Schrödinger cho exciton trong điện trường.
  2. Thiết lập Hamiltonian của hệ exciton hai chiều trong điện trường đều và đưa về dạng toán tử sinh - hủy.
  3. Xây dựng bộ hàm sóng cơ sở và tính toán giải tích các yếu tố ma trận của Hamiltonian.
  4. Xây dựng chương trình tính toán bằng ngôn ngữ lập trình FORTRAN kết hợp thư viện LAPACK để tìm nghiệm số chính xác (trị riêng năng lượng).
  5. Phân tích, so sánh kết quả thu được với nghiệm giải tích và khảo sát sự phụ thuộc của phổ năng lượng vào cường độ điện trường.

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

  • Đối tượng nghiên cứu: Chuẩn hạt exciton hai chiều (mô hình Mott-Wannier) trong trường thế Coulomb tương đối giữa electron và lỗ trống, chịu tác động của điện trường đều đặt song song với mặt phẳng vật liệu.
  • Phạm vi nghiên cứu:
    • Khảo sát phương trình Schrödinger dừng (dựa trên giả thiết thời gian sống của exciton tương đối dài so với thời gian quá độ).
    • Không gian hai chiều chuẩn hóa không thứ nguyên (tọa độ vị trí, hằng số Rydberg hiệu dụng, bán kính Bohr hiệu dụng).
    • Khảo sát hai trường hợp: điện trường bằng không ($F_x = 0, F_y = 0$) để kiểm tra độ chính xác, và điện trường khác không ($F_x \neq 0, F_y = 0$) ở vùng cường độ điện trường nhỏ bắt đầu từ $5 \cdot 10^{-5}$.
    • Cắt cụt không gian vector trạng thái cơ sở tại các giá trị tham số $n_{\max} = 50$ và $n_{\max} = 80$.

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

Khóa luận sử dụng hệ thống lý thuyết lượng tử của chuẩn hạt exciton và phương pháp toán tử đại số phi nhiễu loạn:

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

  • Gần đúng khối lượng hiệu dụng: Mô tả chuyển động của electron và lỗ trống trong trường tuần hoàn của mạng tinh thể thông qua các khối lượng hiệu dụng $m_e$ và $m_h$.
  • Tách chuyển động khối tâm: Chuyển hệ tọa độ của hai hạt $(\vec{r}_e, \vec{r}_h)$ về tọa độ khối tâm $\vec{R} = (m_e \vec{r}_e + m_h \vec{r}_h)/(m_e + m_h)$ và tọa độ chuyển động tương đối $\vec{r} = \vec{r}_e - \vec{r}_h$. Hamiltonian toàn phần phân rã thành $H = H_r + H_R$, trong đó chuyển động khối tâm là hạt tự do mang khối lượng $M^* = m_e + m_h$, còn chuyển động tương đối có khối lượng rút gọn $\mu = (m_e m_h)/(m_e + m_h)$ chịu thế Coulomb $-e^2 / (4\pi\varepsilon_0 r)$ và thế điện trường $e\vec{\mathcal{E}}\cdot\vec{r}$.
  • Hệ đơn vị không thứ nguyên: Chuyển đổi phương trình Schrödinger tương đối về dạng phi thứ nguyên bằng cách chuẩn hóa theo bán kính Bohr hiệu dụng $a^* = 4\pi\varepsilon_0\hbar^2 / (\mu e^2)$, hằng số Rydberg hiệu dụng $R^* = \hbar^2 / (2\mu {a^}^2)$ và véc-tơ cường độ điện trường không thứ nguyên $\vec{F} = e a^ \vec{\mathcal{E}} / R^*$.
  • Phép biến đổi Levi-Civita: Trong bài toán tương tác Coulomb 2D, số hạng thế năng chứa $1/r = 1/\sqrt{x^2+y^2}$ ở mẫu số gây cản trở việc biểu diễn trực tiếp qua các toán tử sinh - hủy điều hòa. Khóa luận sử dụng phép biến đổi Levi-Civita từ không gian tọa độ thực $(x, y)$ sang không gian tham số $(u, v)$: $$x = u^2 - v^2, \quad y = 2uv, \quad r = u^2 + v^2$$ Định thức ma trận Jacobian của phép biến đổi là $J = 4(u^2 + v^2)$. Tích vô hướng bảo toàn được chuyển đổi thông qua điều kiện Hermite, biến đổi phương trình ban đầu thành phương trình dao động tử phi điều hòa hai chiều: $$\tilde{H}\psi(u,v) = 0 \quad \text{với} \quad \tilde{H} = r(H - E)$$

Phương pháp toán tử FK (Feranchuk-Komarov Operator Method)

Khác với thuyết nhiễu loạn thông thường phân chia Hamiltonian dựa trên độ lớn tương tác vật lý (yêu cầu tương tác ngoài phải rất nhỏ), phương pháp toán tử FK phân chia dựa trên hình thức toán tử đại số:

  1. Đưa các biến động lực học $(u, v, \partial/\partial u, \partial/\partial v)$ về biểu diễn toán tử sinh - hủy Dirac.
  2. Tách Hamiltonian thành hai phần: thành phần chính $H_0$ chứa toán tử trung hòa đã có nghiệm giải tích chính xác, và thành phần $V$ đóng vai trò nhiễu loạn.
  3. Đưa vào một tham số biến phân tự do $\omega > 0$ trong phép định nghĩa các toán tử sinh - hủy mới $a(\omega), a^+(\omega), b(\omega), b^+(\omega)$. Tham số $\omega$ không xuất hiện trong Hamiltonian toàn phần nhưng có mặt trong từng thành phần riêng biệt, đóng vai trò điều chỉnh độ lớn tương đối giữa $H_0$ và $V$, mở rộng bán kính hội tụ của chuỗi và tối ưu hóa độ chính xác số học.
  4. Chéo hóa toán tử hình chiếu moment động lượng $L_z = -2(a^+ a - b^+ b)$ để xây dựng bộ hàm sóng cơ sở $|n_1, n_2\rangle$.

Phương pháp tính toán và nguồn dữ liệu

  • Khóa luận thiết lập hệ phương trình trị riêng ma trận tổng quát có dạng: $$A x = \lambda B x$$ Trong đó ma trận $A$ chứa các phần tử $H_{jk, n_1 n_2}$, ma trận $B$ chứa các phần tử $R_{jk, n_1 n_2}$, vector riêng $x$ là các hệ số khai triển hàm sóng, và trị riêng $\lambda$ chính là mức năng lượng $E$ của exciton.
  • Chương trình tính toán được viết bằng ngôn ngữ FORTRAN, tích hợp các chương trình con giải bài toán trị riêng tổng quát từ gói thư viện toán học LAPACK (nằm trong bộ Intel Math Kernel Library).
  • Dữ liệu đối sánh kiểm chứng độ chính xác lấy từ nghiệm giải tích của exciton 2D không điện trường trong công trình của Yang et al. (1991).

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

Chương 1: Exciton hai chiều trong điện trường đều

Chương 1 tập trung xây dựng toàn bộ cơ sở giải tích và mô hình đại số toán tử cho bài toán exciton hai chiều:

  • Nguyên lý phương pháp toán tử FK: Trình bày 4 bước chuẩn hóa gồm: đại số hóa Hamiltonian qua toán tử sinh - hủy; phân tách Hamiltonian thành phần chính và phần nhiễu loạn thông qua tham số tự do $\omega$; xác định nghiệm bậc không từ phương trình trị riêng của toán tử trung hòa; thiết lập hệ phương trình tìm nghiệm số chính xác.
  • Phương trình Schrödinger của exciton 2D trong điện trường đều: Thiết lập Hamiltonian đầy đủ gồm động năng hai hạt, tương tác Coulomb và thế điện trường ngoài $\vec{\mathcal{E}} = (\mathcal{E}_x, 0)$. Sử dụng phép biến đổi tọa độ khối tâm và tương đối để tách thành hai phương trình độc lập. Chuyển phương trình chuyển động tương đối về dạng không thứ nguyên với hằng số điện tích hiệu dụng $Z = 1$.
  • Áp dụng phép biến đổi Levi-Civita: Áp dụng phép chuyển đổi $(x, y) \to (u, v)$ để khử tọa độ $r$ ở mẫu số của thế Coulomb. Dưới tác động của định thức Jacobian $J = 4(u^2 + v^2)$, phương trình Schrödinger được đưa về dạng Hamiltonian của hệ dao động tử phi điều hòa: $$\tilde{H} = -\frac{1}{8}\left(\frac{\partial^2}{\partial u^2} + \frac{\partial^2}{\partial v^2}\right) + F_x (u^4 - v^4) + 2F_y uv(u^2 + v^2) - E(u^2 + v^2) - 1$$
  • Đại số toán tử và tính toán các yếu tố ma trận:
    • Định nghĩa các cặp toán tử sinh - hủy mới $a, a^+, b, b^+$ phụ thuộc tham số $\omega$. Khóa luận chứng minh chi tiết trong Phụ lục 4 rằng khi có điện trường ngoài, giao hoán tử $[H, L_z] \neq 0$ (moment động lượng quỹ đạo không bảo toàn). Dù vậy, bộ hàm riêng của $L_z$ vẫn được sử dụng làm hệ cơ sở trực chuẩn $|n_1, n_2\rangle$.
    • Biểu diễn Hamiltonian thành dạng toán tử: $\tilde{H} = H^* - E R$.
    • Xác định biểu thức giải tích tường minh cho các yếu tố ma trận của $R$ và $H$. Yếu tố ma trận đường chéo của $R$ có dạng $R_{n_1 n_2}^{n_1 n_2} = -(n_1 + n_2 + 1)/\omega$.
    • Yếu tố ma trận của $H$ chứa cả phần thực và phần ảo khi có mặt điện trường ngoài $F_x, F_y$. Về mặt vật lý lượng tử, các yếu tố ma trận phức dẫn đến mức năng lượng phức $E = E_r - i\Gamma/2$, trong đó phần thực $E_r$ xác định vị trí mức năng lượng dịch chuyển, còn phần ảo $\Gamma$ mô tả tốc độ và xác suất ion hóa xuyên ngầm của exciton dưới tác dụng của điện trường.

Chương 2: Kết quả và phân tích

Chương 2 trình bày các kết quả số học thu được từ chương trình FORTRAN/LAPACK qua hai bài toán cụ thể:

1. Trường hợp điện trường bằng không ($F_x = 0, F_y = 0$)

Tác giả tiến hành kiểm tra tính chính xác của chương trình tính toán bằng cách so sánh mức năng lượng exciton 2D với nghiệm giải tích chính xác từ công trình [7] của Yang et al.: $E_n = -1/[2(n - 1/2)^2]$.

  • Khảo sát tham số tự do $\omega$ với $n_{\max} = 50$:
    • Năng lượng trạng thái cơ bản ($n=1$) đạt độ chính xác cao trong dải $\omega \in [1.3; 12]$.
    • Trạng thái kích thích thứ nhất ($n=2$) đạt độ chính xác trong dải $\omega \in [0.4; 4.0]$.
    • Trạng thái kích thích thứ hai ($n=3$) đạt độ chính xác trong dải $\omega \in [0.3; 2.0]$.
    • Kết quả cho thấy khi số lượng tử chính $n$ tăng, khoảng giá trị $\omega$ tối ưu bị thu hẹp dần và dịch chuyển về phía các giá trị nhỏ ($\omega < 1$). Với $n_{\max} = 50$, chương trình tính chính xác đến trạng thái $n=20$ và đạt giới hạn ở $n=25$.
  • Khảo sát ảnh hưởng của việc tăng tham số cắt cụt lên $n_{\max} = 80$:
    • Dải hội tụ $\omega$ của trạng thái $n=2$ mở rộng thành $\omega \in [0.35; 5.0]$; trạng thái $n=3$ mở rộng thành $\omega \in [0.3; 2.4]$.
    • Khả năng tính chính xác được nâng lên đến trạng thái $n=24$ và đạt tối đa ở $n=28$.
    • Rút ra quy luật: Tăng $n_{\max}$ từ 50 lên 80 không làm tăng thêm số chữ số có nghĩa sau dấu phẩy (vốn đã đạt giới hạn chính xác tuyệt đối 13 chữ số thập phân), nhưng có tác dụng mở rộng khoảng tham số $\omega$ cho nghiệm hội tụ và cho phép xác định được các mức kích thích cao hơn.

2. Trường hợp điện trường khác không ($F_x \neq 0, F_y = 0$)

Chương trình khảo sát điện trường hướng dọc trục $Ox$ bắt đầu từ cường độ nhỏ $F_x = 5 \cdot 10^{-5}$:

  • Hiện tượng tách mức năng lượng (Hiệu ứng Stark): Dưới tác dụng của điện trường ngoài, tính suy biến của các trạng thái kích thích bị phá vỡ. Mức năng lượng của trạng thái kích thích thứ $(n-1)$ bị tách thành $(2n-1)$ vạch:
    • Trạng thái cơ bản $n=1$: không bị tách vạch, mức năng lượng dịch chuyển tịnh tiến theo điện trường.
    • Trạng thái kích thích thứ nhất ($n=2$): tách thành 3 vạch năng lượng riêng biệt.
    • Trạng thái kích thích thứ hai ($n=3$): tách thành 5 vạch năng lượng.
    • Trạng thái kích thích thứ ba ($n=4$): tách thành 7 vạch năng lượng.
  • Sự xen phủ mức năng lượng: Khi tiếp tục tăng cường độ điện trường, các nhánh tách mức của các trạng thái kích thích khác nhau có xu hướng dịch chuyển và xen kẽ vào nhau, gây khó khăn cho việc phân định và gán nhãn trạng thái lượng tử nếu chỉ dựa vào giá trị trị riêng số học.

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

Kết quả tính toán định lượng

Trong trường hợp không có điện trường ngoài ($F_x = 0, F_y = 0$), phương pháp toán tử FK cho kết quả trùng khớp hoàn toàn với nghiệm giải tích chính xác của công trình [7] đến 13 chữ số thập phân sau dấu phẩy (khi đưa về dạng chuẩn khoa học):

Trạng thái ($n$) Năng lượng tính bằng FK-OM ($R^*$) Nghiệm giải tích công trình [7] ($R^*$)
$n = 1$ -2,00000000000000 -2,000000000000000000
$n = 2$ -0,22222222222222 -0,222222222222222222
$n = 3$ -0,08000000000000 -0,080000000000000000
$n = 4$ -0,04081632653061 -0,040816326530612245
$n = 5$ -0,02469135802469 -0,024691358024691358
$n = 6$ -0,01652892561983 -0,016528925619834711
$n = 7$ -0,01183431952663 -0,011834319526627219
$n = 8$ -0,00888888888889 -0,008888888888888889
$n = 9$ -0,00692041522491 -0,006920415224913495
$n = 10$ -0,00554016620499 -0,005540166204986150

Trong trường hợp có điện trường ngoài ($F_x \neq 0$), kết quả tính toán số minh chứng cụ thể cho hiện tượng tách mức Stark ở trạng thái kích thích thứ nhất ($n=2$) thành 3 nhánh vạch phổ phụ thuộc vào cường độ điện trường:

Cường độ điện trường ($F_x$) Vạch mức 1 ($R^*$) Vạch mức 2 ($R^*$) Vạch mức 3 ($R^*$)
$5,0 \cdot 10^{-5}$ -0,237226004943158 -0,222222221297009 -0,203101165733148
$1,0 \cdot 10^{-4}$ -0,237226487639043 -0,222222129707035 -0,203100900341492
$1,5 \cdot 10^{-4}$ -0,237227950258994 -0,222221852237120 -0,203101165733148
$2,0 \cdot 10^{-4}$ -0,237274679057189 -0,222213032852702 -0,203074378636976
$2,5 \cdot 10^{-4}$ -0,237419783176928 -0,222186203810867 -0,203048693601639
$3,0 \cdot 10^{-4}$ -0,237656201450600 -0,222143849914059 -0,202861323800138
$3,5 \cdot 10^{-4}$ -0,237987029540199 -0,222089182432556 -0,202676725833476
$4,0 \cdot 10^{-4}$ -0,239454736091213 -0,221891696685213 -0,201828727048276
$4,5 \cdot 10^{-4}$ -0,240084857165791 -0,221828368090637 -0,201454956899363
$5,0 \cdot 10^{-4}$ -0,241522127130521 -0,221726320749949 -0,200585919085719

Đóng góp học thuật

  1. Phát triển phương pháp: Khóa luận đã mở rộng thành công phạm vi ứng dụng của phương pháp toán tử FK sang bài toán exciton hai chiều chịu tác dụng của điện trường đều. Đây là bước phát triển tiếp nối các công trình trước đó vốn chỉ áp dụng FK-OM cho trường ngoài là từ trường.
  2. Khử kỳ dị Coulomb hiệu quả: Việc kết hợp thành công phép biến đổi tọa độ Levi-Civita đã chuyển đổi hiệu quả bài toán thế Coulomb kỳ dị $1/r$ về hệ dao động tử phi điều hòa, cho phép đại số hóa hoàn toàn bài toán qua hệ toán tử sinh - hủy mà không cần sử dụng hệ tọa độ parabol phức tạp hay các xấp xỉ bán cổ điển WKB như các nghiên cứu trước.
  3. Công bố khoa học: Toàn bộ kết quả tính toán biểu thức giải tích của các yếu tố ma trận đã được nghiệm thu và công bố trong bài báo khoa học: Phạm Thị Mỹ Hảo, Nguyễn Thị Thùy Trang, Hoàng Đỗ Ngọc Trâm, "Yếu tố ma trận cho exciton hai chiều trong điện trường", đăng trên Tạp chí Khoa học Đại học Sư phạm TP. Hồ Chí Minh.

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

Hạn chế được ghi nhận trong văn bản

  • Tối ưu hóa tham số $\omega$ khi có điện trường: Tác giả chưa tìm được quy tắc lựa chọn tham số $\omega$ tối ưu trong trường hợp điện trường khác không để đạt độ chính xác cao tương đương trường hợp $F=0$. Nguyên nhân được xác định là do ma trận Hamiltonian khi có điện trường xuất hiện các yếu tố phức, trong khi tham số $\omega$ đưa vào tính toán hiện tại chỉ nhận giá trị thực.
  • Chưa trích xuất xác suất ion hóa: Khóa luận chưa tính toán cụ thể giá trị số của độ rộng mức $\Gamma$ (tỷ lệ ion hóa xuyên ngầm của exciton) từ phần ảo của các yếu tố ma trận phức.
  • Giới hạn phạm vi điện trường: Mới chỉ khảo sát thành công ở dải điện trường có giá trị nhỏ ($F_x \le 5 \cdot 10^{-4}$). Khi điện trường tăng cao, hiện tượng tách mức Stark mạnh làm các vạch phổ của các trạng thái kích thích khác nhau đan xen phức tạp, gây khó khăn cho việc nhận diện và lựa chọn nghiệm số.
  • Giới hạn số trạng thái kích thích: Do giới hạn thời gian tính toán và cấu hình cắt cụt cơ sở ($n_{\max} = 80$), chương trình mới khảo sát được trạng thái kích thích tối đa đến $n=28$.

Hướng phát triển đề tài

  • Nâng cấp chương trình tính toán nhanh và tối ưu hơn thông qua việc khai thác sâu các gói con chuyên dụng khác của thư viện LAPACK.
  • Nghiên cứu thuật toán nhận dạng và tự động phân loại nhánh nghiệm của phương trình đặc trưng nhằm khắc phục khó khăn khi các mức kích thích xen phủ.
  • Mở rộng khảo sát phổ năng lượng và xác suất ion hóa exciton ở các vùng cường độ điện trường lớn hơn.
  • Mở rộng phương pháp toán tử FK cho các hệ hạt phức tạp hơn trong cấu trúc bán dẫn 2D, tiêu biểu là bài toán exciton mang điện tích (trion âm và trion dương hai chiều) đặt trong điện trường ngoài.

Giá trị tham khảo

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

  • Sinh viên, học viên ngành Sư phạm Vật lý, Vật lý lý thuyết và Vật lý tính toán: Cung cấp tài liệu tham khảo chi tiết về kỹ thuật lượng tử hóa đại số, quy trình khử kỳ dị thế Coulomb bằng biến đổi Levi-Civita và phương pháp chuyển bài toán phương trình vi phân đạo hàm riêng sang bài toán ma trận trị riêng tổng quát.
  • Người nghiên cứu vật liệu bán dẫn thấp chiều (2D TMDs): Cung cấp các số liệu định lượng chuẩn xác về phổ năng lượng exciton không trường ngoài và các bảng dữ liệu mô tả sự tách mức Stark trong điện trường đều, làm cơ sở lý thuyết đối sánh với các phép đo phổ quang phát quang (PL) hoặc phổ hấp thụ thực nghiệm.
  • Phần giá trị tham khảo chuyên sâu: Toàn bộ 7 phụ lục tính toán giải tích ở cuối khóa luận (từ Phụ lục 1 đến Phụ lục 7) cung cấp chi tiết từng bước đại số giao hoán tử, công thức tác dụng của toán tử sinh - hủy và biểu thức giải tích của các yếu tố ma trận $H$ và $R$, rất hữu ích cho việc tự xây dựng mã nguồn mô phỏng cơ học lượng tử.

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

1. Phép biến đổi Levi-Civita giải quyết khó khăn toán học nào trong bài toán này?

Trong bài toán exciton 2D, thế tương tác Coulomb chứa số hạng $1/r = 1/\sqrt{x^2+y^2}$ ở mẫu số, gây kỳ dị tại gốc tọa độ và không thể biểu diễn trực tiếp dưới dạng đa thức của các toán tử sinh - hủy điều hòa. Phép biến đổi Levi-Civita chuyển tọa độ $(x, y)$ sang $(u, v)$ theo hệ thức $x = u^2 - v^2$ và $y = 2uv$, dẫn đến khoảng cách $r = u^2 + v^2$. Kết hợp với định thức Jacobian $J = 4(u^2+v^2)$, phương trình Schrödinger được nhân thêm thừa số $r$, triệt tiêu hoàn toàn mẫu số kỳ dị và đưa phương trình về dạng dao động tử phi điều hòa hai chiều.

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

Tham số $\omega$ là một số thực dương được đưa vào khi định nghĩa các toán tử sinh - hủy mới. $\omega$ hoàn toàn không xuất hiện trong Hamiltonian toàn phần nên không làm biến đổi bản chất vật lý của hệ. Tuy nhiên, nó có mặt trong cả thành phần chính $H_0$ lẫn thành phần nhiễu loạn $V$, đóng vai trò như một tham số điều chỉnh độ lớn tương đối giữa hai thành phần này. Việc chọn $\omega$ phù hợp giúp làm nhỏ thành phần nhiễu loạn $V$, tối ưu hóa bán kính hội tụ của chuỗi và giúp thuật toán đạt độ chính xác cao với số lượng hàm cơ sở hữu hạn.

3. Việc tăng tham số cắt cụt từ $n_{\max} = 50$ lên $n_{\max} = 80$ mang lại hiệu quả gì?

Theo khảo sát ở Chương 2, việc tăng $n_{\max}$ từ 50 lên 80 không làm tăng thêm số chữ số thập phân có nghĩa của kết quả (vốn đã đạt chính xác tuyệt đối 13 chữ số thập phân sau dấu phẩy). Tuy nhiên, tăng $n_{\max}$ mang lại hai lợi ích kỹ thuật quan trọng: mở rộng khoảng giá trị $\omega$ cho phép thu được nghiệm hội tụ chính xác, và nâng cao số lượng mức kích thích có thể tính toán được (mở rộng từ trạng thái $n=20$ lên $n=24$, với khả năng khảo sát tối đa đạt $n=28$).

4. Hiệu ứng Stark biểu hiện như thế nào trên phổ năng lượng của exciton 2D?

Khi đặt trong điện trường đều song song, tính đối xứng của hệ bị suy giảm và tính suy biến năng lượng bị phá vỡ. Ngoại trừ trạng thái cơ bản ($n=1$) không bị tách mức, các trạng thái kích thích thứ $(n-1)$ đều bị tách thành $(2n-1)$ vạch phổ năng lượng riêng biệt. Cụ thể trong khóa luận: trạng thái kích thích thứ nhất ($n=2$) tách thành 3 vạch, trạng thái kích thích thứ hai ($n=3$) tách thành 5 vạch, và trạng thái kích thích thứ ba ($n=4$) tách thành 7 vạch.

5. Tại sao các yếu tố ma trận của Hamiltonian lại xuất hiện phần ảo khi có điện trường ngoài?

Khi có điện trường ngoài $F_x, F_y$, các biến đổi đại số của các toán tử sinh - hủy dẫn đến sự xuất hiện của các phần tử ma trận chứa đơn vị ảo $i$. Ý nghĩa vật lý lượng tử của hiện tượng này là hệ không còn ở trạng thái dừng thuần túy mà trở thành trạng thái giả dừng (quasi-steady state) với mức năng lượng phức $E = E_r - i\Gamma/2$. Phần thực $E_r$ thể hiện vị trí mức năng lượng bị dịch chuyển bởi điện trường, còn phần ảo $\Gamma$ tương ứng với độ rộng mức, phản ánh xác suất hạt xuyên ngầm qua hàng rào thế thế năng để ion hóa.

Kết luận

Khóa luận tốt nghiệp của tác giả Phạm Thị Mỹ Hảo đã phát triển thành công phương pháp toán tử FK kết hợp với phép biến đổi tọa độ Levi-Civita để giải phương trình Schrödinger cho exciton hai chiều đặt trong điện trường đều. Kết quả tính toán số học trên ngôn ngữ FORTRAN sử dụng thư viện LAPACK cho thấy độ chính xác đạt tới 13 chữ số thập phân ở trường hợp không điện trường, đồng thời mô tả rõ nét quy luật tách mức Stark $(2n-1)$ vạch ở các trạng thái kích thích khi có điện trường ngoài. Công trình đóng góp một hướng tiếp cận đại số phi nhiễu loạn hiệu quả cho vật lý tính toán các hệ bán dẫn thấp chiều, đồng thời mở ra hướng nghiên cứu tiếp theo về xác suất ion hóa và bài toán các phức hợp trion trong điện trường.