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

Trong vật lý chất rắn hiện đại, việc nghiên cứu các hệ vật lý thấp chiều và vật liệu cấu trúc nano đạt được nhiều bước tiến nhờ các kỹ thuật chế tạo tiên tiến như nuôi cấy chùm phân tử (MBE) và lắng đọng hơi hóa học kim loại - hữu cơ (MOCVD). Trong các mô hình cấu trúc thấp chiều thực nghiệm, dị thể bán dẫn nhiều lớp $\text{GaAs}/\text{Al}x\text{Ga}{1-x}\text{As}$ được sử dụng phổ biến do cho phép kiểm soát chặt chẽ nồng độ và đặc tính hạt tải khi điều chỉnh thành phần $x$. Trong cấu trúc này, vùng $\text{GaAs}$ đóng vai trò hố thế, còn $\text{Al}x\text{Ga}{1-x}\text{As}$ đóng vai trò rào thế. Khi lớp vật liệu đủ mỏng, chuyển động của điện tử bị giam cầm và được mô hình hóa như chuyển động trong không gian hai chiều (2D). Sự xuất hiện của các đỉnh nhọn (peak) trong phổ hấp thụ quang học của chất bán dẫn ở nhiệt độ thấp là bằng chứng thực nghiệm khẳng định sự tồn tại của exciton – một trạng thái liên kết giả hạt giữa electron và lỗ trống thông qua lực hút Coulomb.

Khi đặt hệ bán dẫn vào từ trường ngoài, exciton thể hiện nhiều hiệu ứng vật lý quan trọng như hiệu ứng Stark, hiệu ứng tách vạch Zeeman, biến đổi độ dẫn điện và sự tái cấu trúc các mức năng lượng Landau do tương tác Coulomb. Do đó, việc giải phương trình Schrödinger để xác định chính xác phổ năng lượng và hàm sóng của exciton hai chiều trong từ trường có ý nghĩa lý thuyết và ứng dụng thiết thực trong thiết kế linh kiện quang điện tử nano.

Trước đây, bài toán exciton 2D trong từ trường đã được giải quyết bằng phương pháp biến phân kết hợp khai triển chuỗi $1/N$, hoặc phương pháp lý thuyết nhiễu loạn có tính đến tiệm cận của hàm sóng, cho độ chính xác đến 7 chữ số thập phân đối với trạng thái cơ bản $1s$ và các trạng thái kích thích thấp ($2p^-, 3d^-$). Phương pháp toán tử do Feranchuk và Komarov (gọi tắt là phương pháp FK) cũng đã được áp dụng kết hợp với phép biến đổi tọa độ Levi-Civita, đạt độ chính xác lên đến 20 chữ số thập phân. Tuy nhiên, phép biến đổi Levi-Civita biến đổi năng lượng $E$ thành một tham số nằm trong toán tử Hamilton, khiến việc tìm năng lượng phải thông qua một trị riêng hình thức $Z$ bằng phương trình gián tiếp $Z(E) = \text{const}$. Cách tiếp cận gián tiếp này gặp nhiều bất lợi khi xây dựng thuật toán số và khó mở rộng cho các hệ phức tạp hơn như exciton tích điện âm ($X^-$). Ngược lại, phép biến đổi tích phân Laplace cho phép đưa biến tọa độ ở mẫu số của thế Coulomb lên tử số mà không cần đưa vào biến hình thức trung gian, giúp giải trực tiếp trị riêng năng lượng. Dù khối lượng giải tích ban đầu tăng lên, phép biến đổi Laplace trước đó mới chỉ được thử nghiệm trên exciton tự do hoặc exciton âm mà chưa từng được áp dụng cho exciton trung hòa trong từ trường đều.

Từ bối cảnh học thuật trên, đề tài khóa luận xác định mục tiêu và các nhiệm vụ nghiên cứu cụ thể:

  • Mục tiêu nghiên cứu: Áp dụng phương pháp toán tử FK kết hợp với phép biến đổi Laplace để tìm nghiệm số chính xác cho bài toán exciton hai chiều trong từ trường đều.
  • Nhiệm vụ nghiên cứu:
    1. Tìm hiểu tổng quan lý thuyết về exciton trong vật lý bán dẫn.
    2. Nghiên cứu cơ sở phương pháp toán tử FK và các giải pháp xử lý toán học khi áp dụng cho các hệ nguyên tử, phân tử.
    3. Xây dựng chương trình tính số để tìm nghiệm chính xác năng lượng và hàm sóng cho exciton trung hòa 2D trong từ trường ở trạng thái cơ bản và một số trạng thái kích thích.
  • Đối tượng và phạm vi nghiên cứu: Đối tượng nghiên cứu là exciton trung hòa chuyển động hai chiều dưới tương tác Coulomb và từ trường đều vuông góc với mặt phẳng chuyển động; phạm vi nghiên cứu mở rộng trên toàn bộ dải cường độ từ trường từ vùng từ trường yếu đến vùng từ trường rất mạnh, không giới hạn ở điều kiện gần đúng nhiễu loạn nhỏ.

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

Khung lý thuyết của đề tài dựa trên cơ học lượng tử phi tương đối tính áp dụng cho exciton loại Wannier-Mott trong chất bán dẫn. Khác với exciton Frenkel trong chất cách điện có bán kính liên kết nhỏ cỡ ô mạng cơ sở và năng lượng liên kết lớn (~1 eV), exciton Wannier-Mott có bán kính lớn hơn nhiều lần hằng số mạng tinh thể và năng lượng liên kết nhỏ (cỡ vài meV đến vài chục meV) do hằng số điện môi $\varepsilon$ lớn và khối lượng hiệu dụng của electron ($m_e$) cũng như lỗ trống ($m_h$) nhỏ.

Toán tử Hamilton của hệ hai hạt tương tác trong từ trường ngoài đều $\mathbf{B} = (0, 0, B)$ được thiết lập qua thế vectơ đối xứng $\mathbf{A} = \frac{1}{2}[\mathbf{B} \times \mathbf{r}]$. Tác giả thực hiện phép chuyển đổi hệ tọa độ sang hệ quy chiếu khối tâm:

  • Tọa độ tương đối: $\mathbf{r} = \mathbf{r}_e - \mathbf{r}_h$ với khối lượng rút gọn $\mu = m_e m_h / (m_e + m_h)$.
  • Tọa độ khối tâm: $\mathbf{R} = (m_e \mathbf{r}_e + m_h \mathbf{r}_h) / (m_e + m_h)$ với tổng khối lượng $M = m_e + m_h$.

Phần toán tử Hamilton mô tả chuyển động khối tâm $H_{c.m.}$ có dạng dao động tử điều hòa đã biết nghiệm giải tích chính xác. Phần toán tử Hamilton mô tả chuyển động tương đối $H_r$ của cặp electron - lỗ trống được đưa về hệ đơn vị không thứ nguyên, trong đó đơn vị năng lượng là hằng số Rydberg hiệu dụng $R^* = \mu e^4 / (2\varepsilon^2 \hbar^2)$, đơn vị chiều dài là bán kính Bohr hiệu dụng $a^* = \varepsilon \hbar^2 / (\mu e^2)$, và cường độ từ trường không thứ nguyên được xác định bởi $\gamma = \hbar \omega_c / (2 R^*)$ (với $\omega_c = eB/(\mu c)$ là tần số cyclotron). Để thuận tiện cho việc khảo sát toàn dải từ trường từ 0 đến vô hạn, tác giả sử dụng biến số chuẩn hóa $\gamma' = \gamma / (\gamma + 1) \in [0, 1)$.

Phương pháp nghiên cứu cốt lõi là phương pháp toán tử FK (xây dựng bởi nhóm tác giả Feranchuk và Komarov từ năm 1979). Ý tưởng căn bản của phương pháp FK là đưa toán tử Hamilton về biểu diễn đại số của các toán tử sinh - hủy, sau đó phân rã $H = H_0 + V$. Khác với lý thuyết nhiễu loạn truyền thống (vốn phân chia dựa trên độ lớn của tương tác vật lý và đòi hỏi tham số nhiễu loạn phải nhỏ), phương pháp FK phân chia thuần túy dựa trên cấu trúc toán học của các toán tử:

  • Phần chính $H_0$ được tạo thành từ các toán tử "trung hòa" (có số toán tử sinh và hủy bằng nhau), do đó có nghiệm giải tích chính xác đóng vai trò nghiệm gần đúng bậc không.
  • Phần còn lại $V$ đóng vai trò nhiễu loạn, trong đó các tham số tương tác (như cường độ trường) xuất hiện đồng thời ở cả $H_0$ và $V$.
  • Một tham số điều chỉnh tự do $\omega$ được đưa vào định nghĩa toán tử sinh - hủy. Do toán tử Hamilton thực tế không phụ thuộc vào $\omega$, giá trị của $\omega$ được tối ưu hóa từ điều kiện cực trị năng lượng bậc không $\partial E^{(0)} / \partial \omega = 0$ để triệt tiêu hoặc cực tiểu hóa đóng góp của phần nhiễu loạn $V$, đảm bảo tốc độ hội tụ nhanh.

Đối với thế tương tác Coulomb có biến động lực nằm ở mẫu số ($1/r$), tác giả sử dụng phép biến đổi tích phân Laplace: $$\frac{1}{r} = \frac{1}{\sqrt{\pi}} \int_0^\infty t^{-1/2} e^{-t r^2} , dt$$ Phép biến đổi này chuyển phần thế năng tương tác Coulomb về dạng hàm mũ của toán tử bình phương tọa độ, cho phép khai triển đại số theo các toán tử sinh - hủy mà không phải thông qua biến số hình thức phụ.

Về phương diện tính toán số, tác giả sử dụng giải thuật sơ đồ vòng lặp ma trận (iterative scheme) bậc cao thay vì tính bổ chính giải tích từng bậc của thuyết nhiễu loạn, và lập trình giải số bằng ngôn ngữ Fortran.

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

Chương 1: Tổng quan về exciton

Chương này hệ thống hóa lịch sử phát hiện và bản chất vật lý của exciton từ các mốc thực nghiệm đầu tiên: phổ hấp thụ khí hiếm của Becquerel (1907), tinh thể phân tử của Obreimov và De Haas (1929), mô hình kích thích không dịch chuyển điện tích của Frenkel (1931), mô hình exciton bán kính lớn của Mott và Wannier (1937), thực nghiệm phát hiện vạch quang phổ giống Hydro trên tinh thể $\text{Cu}_2\text{O}$ của Gross (1951), và dự đoán exciton mang điện của Lampert (1958).

Tác giả phân tích cơ chế hình thành ba trạng thái điện tích của exciton: exciton trung hòa $X_0$, exciton âm $X^-$ (trạng thái liên kết hai electron và một lỗ trống), và exciton dương $X^+$ (hai lỗ trống liên kết với một electron). Trọng tâm chương tập trung vào hệ bán dẫn hai chiều $\text{GaAs}/\text{Al}x\text{Ga}{1-x}\text{As}$, mô tả thế giếng lượng tử giam cầm điện tử thành khí điện tử 2D.

Phần trọng tâm giải tích của Chương 1 là việc thiết lập phương trình Schrödinger cho exciton 2D trong từ trường. Tác giả chứng minh chi tiết việc chọn chuẩn thế vectơ $\text{div}\mathbf{A} = 0$, tính toán giao hoán tử động lượng và từ trường, chuyển đổi hệ tọa độ sang hệ quy chiếu khối tâm $(r, R)$, tách chuyển động khối tâm và chuyển động tương đối. Cuối cùng, tác giả dẫn ra phương trình Schrödinger không thứ nguyên mô tả chuyển động tương đối: $$\left[ -\left( \frac{\partial^2}{\partial x^2} + \frac{\partial^2}{\partial y^2} \right) + \frac{\gamma^2}{4}(x^2 + y^2) + \gamma L_z - \frac{2}{\sqrt{x^2 + y^2}} \right] \psi(x, y) = E \psi(x, y)$$ trong đó $L_z = -i(x \partial/\partial y - y \partial/\partial x)$ là hình chiếu mômen động lượng quỹ đạo lên trục từ trường.

Chương 2: Phương pháp toán tử FK

Chương 2 trình bày có hệ thống cơ sở lý thuyết và quy trình 4 bước của phương pháp toán tử FK thông qua bài toán mẫu dao động tử phi điều hòa $H = -\frac{1}{2}\frac{d^2}{dx^2} + \frac{1}{2}x^2 + \lambda x^4$ ($\lambda > 0$):

  1. Bước 1: Đưa tọa độ và đạo hàm về các toán tử sinh $a^+$ và hủy $a$ chứa tham số thực $\omega > 0$, đưa toán tử Hamilton về dạng chuẩn (toán tử sinh nằm bên trái, toán tử hủy nằm bên phải).
  2. Bước 2: Phân tách Hamilton thành phần chính $H_0(a^+a, \lambda, \omega)$ chứa các toán tử trung hòa và phần nhiễu loạn $V(a^+, a, \lambda, \omega)$.
  3. Bước 3: Giải phương trình trị riêng cho $H_0$ để tìm phổ năng lượng bậc không $E_n^{(0)}$ trên hệ trạng thái cơ sở dao động tử điều hòa $|n\rangle$; thiết lập điều kiện $\partial E_n^{(0)} / \partial \omega = 0$ để xác định phương trình đại số bậc 3 tìm tham số tối ưu $\omega$.
  4. Bước 4: Xây dựng sơ đồ vòng lặp số để xác định năng lượng chính xác $E_n$ và hệ số khai triển hàm sóng mà không dùng tham số khai triển nhỏ.

Tác giả so sánh tốc độ hội tụ giữa phương pháp FK và lý thuyết nhiễu loạn tiêu chuẩn. Kết quả cho thấy phương pháp FK hội tụ chính xác đến 10 chữ số thập phân cho mọi giá trị $\lambda$, kể cả trạng thái kích thích ($n=4$). Ngược lại, lý thuyết nhiễu loạn phân kỳ khi $\lambda \ge 0.1$ ở trạng thái cơ bản và phân kỳ ngay từ $\lambda = 0.03$ ở trạng thái kích thích.

Cuối chương, tác giả khảo sát 4 vấn đề kỹ thuật khi mở rộng phương pháp FK cho hệ nguyên tử, phân tử:

  • Khử tọa độ ở mẫu số (phân tích ưu - nhược điểm giữa biến đổi Levi-Civita và biến đổi tích phân Laplace).
  • Đưa toán tử hàm mũ chứa toán tử sinh - hủy về dạng chuẩn thông qua lý thuyết đại số Lie kín.
  • Xây dựng hệ hàm sóng cơ sở thỏa mãn tính đối xứng và các định luật bảo toàn của hệ vật lý (bảo toàn hình chiếu mômen động lượng $L_z$).
  • Lựa chọn tham số tự do $\omega$ và cấu trúc sơ đồ vòng lặp để bảo đảm tốc độ hội tụ.

Chương 3: Phương pháp toán tử FK giải bài toán exciton 2D trong từ trường đều

Đây là chương trọng tâm thực hiện các tính toán mới của khóa luận. Tác giả áp dụng quy trình FK kết hợp biến đổi Laplace giải phương trình Schrödinger cho exciton 2D:

  • Biểu diễn đại số: Đưa bài toán hai chiều về hai cặp toán tử sinh - hủy $(a, a^+)$ theo phương $x$ và $(b, b^+)$ theo phương $y$ với $\omega_x = \omega_y = \omega$. Sau khi đưa thế tương tác Coulomb qua tích phân Laplace, biểu thức toán tử hàm mũ được chuẩn hóa thành dạng giải tích dựa trên đại số kín của các toán tử $M = \frac{1}{2}(a+ib)^2$, $M^+ = \frac{1}{2}(a^+-ib^+)^2$ và $N = a^+a + b^+b + 1$.
  • Phân tách toán tử: Tách toán tử Hamilton thành phần chính $H_0$ chứa các toán tử trung hòa $a^+a, b^+b$ và phần nhiễu loạn $V$.
  • Hệ hàm sóng cơ sở và nghiệm bậc không: Xây dựng hệ hàm sóng cơ sở trực chuẩn $|k, m\rangle$ là hàm riêng của $H_0$ và $L_z$, trong đó $m$ là số lượng tử hình chiếu mômen quỹ đạo ($m = 0, \pm 1, \pm 2, \dots$), $k$ là số lượng tử bán kính ($k = 0, 1, 2, \dots$). Năng lượng bậc không được xác định bởi: $$E_{k,m}^{(0)} = H_{0,kk} = \left(\frac{1}{\omega} + \frac{\gamma^2 \omega}{4}\right)(2k + |m| + 1) + \gamma m - I_0(k, m, \omega)$$ Giá trị $\omega$ được xác định ban đầu từ phương trình $\partial E_{k,m}^{(0)} / \partial \omega = 0$.
  • Sơ đồ vòng lặp ma trận và tính số: Thiết lập công thức ma trận đầy đủ cho $H_{kk'}$ và $V_{kk'}$, lập trình trên Fortran thực hiện phép lặp đến bậc $s = 300$. Tác giả tính toán phổ năng lượng cho các trạng thái $1s$ ($k=0, m=0$), $2p^-$ ($k=0, m=-1$), $3d^-$ ($k=0, m=-2$) và $5d^-$ ($k=2, m=-2$) trên dải từ trường $\gamma' \in [0, 0.9]$.
  • Ý nghĩa các số lượng tử: So sánh với kết quả phép biến đổi Levi-Civita của công trình trước, tác giả thiết lập mối liên hệ tường minh: $N = k + |m| + 1$ (với $N$ là số lượng tử chính). Do đó, $k = N - |m| - 1$ chính là số lượng tử hướng tâm (bán kính), còn $m$ là số lượng tử quỹ đạo.

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

Nghiên cứu đạt được các kết quả định lượng cụ thể sau 300 vòng lặp tính toán trên chương trình Fortran:

  • Đối với trạng thái cơ bản $1s$ ($k=0, m=0$), giá trị năng lượng hội tụ ổn định và khớp với công trình giải bằng biến đổi Levi-Civita ở 1 đến 2 chữ số thập phân đầu tiên.
  • Đối với các trạng thái kích thích ($2p^-, 3d^-, 5d^-$), tốc độ hội tụ nhanh hơn rõ rệt so với trạng thái cơ bản, đạt độ chính xác từ 5 đến 7 chữ số thập phân (một số điểm đạt tới 9 chữ số thập phân khi tinh chỉnh tham số $\omega$).

Bảng dưới đây tổng hợp so sánh mức năng lượng tính bằng phương pháp FK kết hợp biến đổi Laplace với công trình đối chiếu sử dụng biến đổi Levi-Civita [10] tại vòng lặp thứ 300:

Trạng thái Cấu hình $(k, m)$ Tham số từ trường $\gamma'$ Năng lượng FK (Laplace) Năng lượng đối chiếu [10] Độ hội tụ / Sai số tương đối
$1s$ $k=0, m=0$ $0.1$ Số liệu hội tụ vòng 300 Trùng khớp 2 chữ số thập phân $\sim 10^{-2}$
$1s$ $k=0, m=0$ $0.5$ Số liệu hội tụ vòng 300 Trùng khớp 2 chữ số thập phân $\sim 10^{-2}$
$2p^-$ $k=0, m=-1$ $0.1$ Số liệu hội tụ vòng 300 Trùng khớp 5-6 chữ số thập phân $\sim 10^{-5} - 10^{-6}$
$2p^-$ $k=0, m=-1$ $0.7$ Số liệu hội tụ vòng 300 Trùng khớp 5-6 chữ số thập phân $\sim 10^{-5} - 10^{-6}$
$3d^-$ $k=0, m=-2$ Các giá trị $\gamma'$ Số liệu hội tụ vòng 300 Trùng khớp 5-7 chữ số thập phân $\sim 10^{-6} - 10^{-7}$
$5d^-$ $k=2, m=-2$ Các giá trị $\gamma'$ Số liệu hội tụ vòng 300 Trùng khớp 5-7 chữ số thập phân Sai số $\sim 0.1 \times 10^{-6}$

Quy luật vật lý rút ra từ đồ thị phổ năng lượng:

  1. Khi cường độ từ trường $\gamma$ tăng lên, năng lượng của tất cả các trạng thái exciton đều tăng đơn điệu.
  2. Trong vùng từ trường mạnh ($\gamma' \to 1$), tương tác Coulomb đóng vai trò phá vỡ sự suy biến của các mức năng lượng Landau, làm tách rõ các vạch phổ.
  3. Đối với các trạng thái có cùng số lượng tử chính $N$, trạng thái nào có số lượng tử từ $m$ lớn hơn thì sở hữu mức năng lượng cao hơn.

Bảng quy đổi và xác định ý nghĩa các số lượng tử giữa hai phương pháp:

| Trạng thái phổ | Số lượng tử chính $N$ | Số lượng tử từ $m$ | Số lượng tử bán kính $k$ ($k = N - |m| - 1$) | Ký hiệu quy ước | | :--- | :--- | :--- | :--- | :--- | | $1s$ | 1 | 0 | 0 | $s$ ứng với $m=0$ | | $2p^-$ | 2 | -1 | 0 | $p^-$ ứng với $m=-1$ | | $3d^-$ | 3 | -2 | 0 | $d^-$ ứng với $m=-2$ | | $4d^-$ | 4 | -2 | 1 | $d^-$ ứng với $m=-2$ | | $5d^-$ | 5 | -2 | 2 | $d^-$ ứng với $m=-2$ | | $6p^+$ | 6 | +1 | 4 | $p^+$ ứng với $m=+1$ |

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

  • Lần đầu tiên triển khai thành công việc kết hợp phép biến đổi tích phân Laplace với phương pháp toán tử FK cho bài toán exciton trung hòa 2D trong từ trường đều.
  • Khắc phục nhược điểm của phép biến đổi Levi-Civita: tìm trực tiếp trị riêng năng lượng mà không cần giải phương trình hàm ẩn qua biến hình thức $Z(E) = \text{const}$, tạo tiền đề thuận lợi cho việc lập trình mở rộng sang các bài toán phức tạp hơn như exciton tích điện âm ($X^-$).
  • Xác lập tường minh quy tắc chuyển đổi và ý nghĩa vật lý của bộ số lượng tử $(k, m)$ trong biểu diễn Laplace sang bộ số lượng tử bán kính và quỹ đạo chuẩn $(N, n_r, m)$.

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

Khóa luận chỉ ra một số hạn chế cụ thể trong quá trình tính toán:

  • Tốc độ xử lý của chương trình Fortran còn chậm, do đó số bước lặp của sơ đồ vòng lặp chỉ mới dừng lại ở $s = 300$.
  • Tốc độ hội tụ của trạng thái cơ bản $1s$ chậm hơn đáng kể so với các trạng thái kích thích, tại vòng lặp 300 mới chỉ đạt độ chính xác từ 1 đến 2 chữ số thập phân.
  • Việc lựa chọn tham số tự do $\omega$ từ điều kiện đạo hàm bậc không $\partial E^{(0)} / \partial \omega = 0$ chưa phải là tối ưu tuyệt đối cho các trạng thái kích thích bậc cao; tác giả vẫn phải sử dụng phương pháp dò thử số học để tìm giá trị $\omega$ tối ưu lân cận nhằm đảm bảo tốc độ hội tụ.

Các hướng phát triển tiếp theo được đề xuất:

  • Cải tiến cấu trúc thuật toán và tối ưu hóa mã nguồn chương trình tính số để tăng tốc độ xử lý, cho phép tăng số vòng lặp $s$ lên cao hơn trong thời gian tính toán ngắn.
  • Ứng dụng các thư viện tính toán số học đa độ chính xác (Multi Precision) để khai thác và mở rộng số chữ số thập phân hội tụ sau dấu phẩy.
  • Thực hiện nghiên cứu so sánh định lượng chi tiết giữa sơ đồ vòng lặp (iterative scheme) và sơ đồ lý thuyết nhiễu loạn truyền thống trong khuôn khổ phương pháp toán tử FK.

Giá trị tham khảo

Tài liệu là tài liệu chuyên khảo hữu ích cho sinh viên, học viên cao học và nghiên cứu sinh thuộc các chuyên ngành Sư phạm Vật lý, Vật lý lý thuyết, Vật lý chất rắn và Vật lý tính toán đang nghiên cứu về cơ học lượng tử ứng dụng và vật lý hệ thấp chiều.

Các nội dung có giá trị tham khảo chuyên sâu bao gồm:

  • Toàn bộ quy trình giải tích khử kỳ dị thế Coulomb $1/r$ bằng phép biến đổi Laplace và kỹ thuật chuẩn hóa toán tử hàm mũ trên đại số Lie kín (trình bày chi tiết tại Chương 3 và các Phụ lục 4, 5).
  • Kỹ thuật xây dựng hệ hàm sóng cơ sở bảo toàn tính đối xứng hai chiều và hình chiếu mômen động lượng $L_z$ (Phụ lục 6).
  • Thuật toán số vòng lặp tính ma trận Hamilton bậc cao không cần tính tích phân số phức tạp, có thể tái lập trên các ngôn ngữ lập trình khoa học như Fortran, C/C++ hoặc Mathematica.

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

1. Điểm khác biệt căn bản giữa phương pháp toán tử FK và lý thuyết nhiễu loạn truyền thống là gì?
Lý thuyết nhiễu loạn truyền thống đòi hỏi phải tách toán tử Hamilton dựa trên độ lớn vật lý của các số hạng ($H = H_0 + V$ với $V \ll H_0$), và chuỗi bổ chính sẽ phân kỳ khi tương tác trường ngoài lớn. Ngược lại, phương pháp FK phân chia toán tử hoàn toàn dựa trên hình thức toán học: phần $H_0$ thu gom tất cả các toán tử trung hòa (có số toán tử sinh và hủy bằng nhau) nên luôn tìm được nghiệm giải tích chính xác; tham số tương tác xuất hiện đồng thời ở cả phần chính và phần nhiễu loạn; đồng thời việc đưa vào tham số tự do $\omega$ được tối ưu hóa giúp triệt tiêu phần nhiễu loạn, cho phép tìm nghiệm hội tụ ở mọi cường độ trường bất kỳ.

2. Vì sao tác giả lựa chọn phép biến đổi Laplace thay cho phép biến đổi Levi-Civita dù khối lượng tính toán ban đầu lớn hơn?
Phép biến đổi Levi-Civita đưa bài toán về dạng dao động tử phi điều hòa nhưng làm thay đổi vai trò của năng lượng: $E$ trở thành một tham số bên trong toán tử và phải tìm gián tiếp thông qua phương trình ẩn $Z(E) = \text{const}$. Phép biến đổi Laplace biểu diễn trực tiếp thế $1/r$ qua tích phân suy rộng của toán tử hàm mũ, cho phép giữ nguyên $E$ là trị riêng trực tiếp của phương trình Schrödinger. Điều này giúp giải thuật tính số trực tiếp hơn, không phụ thuộc vào các biến hình thức phụ và thuận lợi hơn khi phát triển thuật toán cho các hệ phức hợp nhiều hạt như exciton âm ($X^-$).

3. Ý nghĩa vật lý của hai số lượng tử $k$ và $m$ trong luận văn được xác lập như thế nào?
Thông qua việc đối chiếu phổ năng lượng với kết quả giải bằng phép biến đổi Levi-Civita, tác giả xác định $m$ chính là số lượng tử từ (hình chiếu mômen động lượng quỹ đạo lên trục $z$, nhận các giá trị $0, \pm 1, \pm 2, \dots$), tương ứng với các trạng thái quang học $s, p, d$. Số lượng tử $k$ đóng vai trò là số lượng tử bán kính ($k = n_r \ge 0$). Số lượng tử chính $N$ của hệ được xác định theo hệ thức: $N = k + |m| + 1$.

4. Tại sao ở vòng lặp thứ 300, trạng thái cơ bản $1s$ chỉ hội tụ đến 2 chữ số thập phân trong khi các trạng thái kích thích đạt độ chính xác từ 5 đến 7 chữ số thập phân?
Trạng thái cơ bản $1s$ chịu ảnh hưởng mạnh nhất của tương tác Coulomb tầm ngắn tại gốc tọa độ, khiến phần đóng góp của các yếu tố ma trận không trung hòa trong khai triển Laplace giảm chậm hơn qua từng bước lặp ma trận. Các trạng thái kích thích có hàm sóng phân bố rộng hơn, bán kính trung bình lớn hơn nên tốc độ suy giảm của các yếu tố ma trận nhiễu loạn diễn ra nhanh hơn, dẫn đến tốc độ hội tụ của sơ đồ vòng lặp đối với các trạng thái kích thích nhanh hơn đáng kể so với trạng thái cơ bản.

Kết luận

Khóa luận tốt nghiệp đã vận dụng thành công phương pháp toán tử FK kết hợp với phép biến đổi tích phân Laplace để giải phương trình Schrödinger cho exciton hai chiều trong từ trường đều trên toàn bộ dải cường độ trường. Bằng các phép biến đổi đại số thuần túy và thuật toán vòng lặp ma trận lập trình trên Fortran, tác giả đã tính toán được nghiệm số năng lượng cho cả trạng thái cơ bản và các trạng thái kích thích với độ chính xác cao mà không cần giải gián tiếp qua phương trình trị riêng hình thức. Nghiên cứu đồng thời làm sáng tỏ mối liên hệ định lượng và ý nghĩa vật lý của bộ số lượng tử bán kính và quỹ đạo trong mô hình biến đổi Laplace. Kết quả này khẳng định tính khả thi của hướng tiếp cận bằng biến đổi Laplace trong việc khảo sát các cấu trúc lượng tử thấp chiều có tương tác Coulomb.