Nghiên Cứu Bài Toán Exciton 2D Trong Từ Trường Đều

Khóa luận khảo sát sự hội tụ của phương pháp toán tử fk cho bài toán exciton 2d trong từ trường đều, mang lại cái nhìn sâu sắc về vật lý.

• Tác giả: Kiên Thị Bích Trâm
• Người hướng dẫn: ThS. Hoàng Đỗ Ngọc Trầm
• Trường học: Trường Đại Học Sư Phạm Thành Phố Hồ Chí Minh
• Chuyên ngành: Sư Phạm Vật Lý
• Đề tài: Nghiên Cứu Bài Toán Exciton 2D Trong Từ Trường Đều
• Loại tài liệu: luận văn tốt nghiệp đại học
• Năm xuất bản: 2013
• Địa điểm: Thành phố Hồ Chí Minh

Mục lục chi tiết

• LỜI CẢM ƠN
• MỞ ĐẦU
• 1. CHƯƠNG 1: CƠ SỞ LÝ THUYẾT
• 1.1. Phương pháp toán tử giải phương trình Schrödinger
• 1.2. Tổng quan về exciton
• 1.3. Phương trình Schrödinger cho exciton 2D trong từ trường đều
• 2. CHƯƠNG 2: PHƯƠNG PHÁP TOÁN TỬ GIẢI BÀI TOÁN EXCITON HAI CHIỀU TRONG TỪ TRƯỜNG ĐỀU
• 2.1. Phương pháp toán tử giải bài toán exciton 2D trong từ trường đều
• 2.2. Kết quả - Phân tích
• 3. CHƯƠNG 3: VAI TRÒ CỦA THAM SỐ TỰ DO ĐỐI VỚI SỰ HỘI TỤ CỦA PHƯƠNG PHÁP TOÁN TỬ
• 3.1. Vai trò tham số tự do đối với sự hội tụ của phương pháp toán tử
• 3.2. Sự phụ thuộc của tốc độ hội tụ vào tham số tự do với bài toán dao động tử phi điều hòa bậc bốn
• 3.3. Khảo sát bài toán exciton 2D trong từ trường đều
• 3.3.1. Khảo sát tốc độ hội tụ của bài toán theo các giá trị khác nhau
• 3.3.2. Điều kiện để chọn tham số tự do tối ưu
• KẾT LUẬN VÀ HƯỚNG PHÁT TRIỂN CỦA ĐỀ TÀI
• TÀI LIỆU THAM KHẢO
• PHỤ LỤC 1: CÁC TOÁN TỬ SINH - HỦY MẬT CHIỀU
• PHỤ LỤC 2: DẠNG CHUẨN CỦA TOÁN TỬ
• PHỤ LỤC 3: ĐƯA HAMILTONIAN VỀ DẠNG KHÔNG THỨ NGUYÊN
• PHỤ LỤC 4: CÁC TOÁN TỬ SINH - HỦY HAI CHIỀU
• PHỤ LỤC 5: DẠNG CHUẨN CỦA TOÁN TỬ
• PHỤ LỤC 6: CÁC THÀNH PHẦN MA TRẬN CHO BÀI TOÁN EXCITON 2D TRONG TỪ TRƯỜNG

Tóm tắt

I. Tổng quan về Nghiên Cứu Bài Toán Exciton 2D Trong Từ Trường Đều

Nghiên cứu về exciton 2D trong từ trường đều đang thu hút sự chú ý của nhiều nhà khoa học. Bài toán này không chỉ có ý nghĩa lý thuyết mà còn mang lại nhiều ứng dụng thực tiễn trong công nghệ vật liệu mới. Exciton là trạng thái liên kết giữa điện tử và lỗ trống, có khả năng mang và truyền kích thích trong mạng mà không lan truyền điện tích. Việc nghiên cứu tính chất của exciton trong môi trường 2D giúp hiểu rõ hơn về các hiện tượng vật lý phức tạp trong các hệ vật liệu này.

1.1. Khái niệm và Lịch sử Nghiên Cứu Exciton

Khái niệm về exciton đã được phát triển từ những năm 1960. Các nghiên cứu ban đầu tập trung vào tính chất của exciton trong các vật liệu ba chiều. Tuy nhiên, với sự phát triển của công nghệ, nghiên cứu đã chuyển sang các hệ 2D, nơi mà tương tác giữa exciton và từ trường trở nên quan trọng hơn bao giờ hết.

1.2. Tầm quan trọng của Nghiên Cứu Exciton 2D

Nghiên cứu exciton 2D trong từ trường đều không chỉ giúp hiểu rõ hơn về các hiện tượng vật lý mà còn mở ra hướng đi mới cho các ứng dụng trong công nghệ nano. Các hệ vật liệu 2D như graphene và dichalcogenides có tiềm năng lớn trong việc phát triển các thiết bị điện tử và quang điện mới.

II. Vấn Đề và Thách Thức Trong Nghiên Cứu Exciton 2D

Mặc dù có nhiều tiềm năng, nhưng nghiên cứu về bài toán exciton 2D trong từ trường đều cũng gặp phải nhiều thách thức. Một trong những vấn đề chính là việc giải phương trình Schrödinger cho các hệ này, đặc biệt là trong trường hợp không thể tìm được nghiệm chính xác. Các phương pháp gần đúng thường được sử dụng nhưng vẫn còn nhiều hạn chế.

2.1. Khó Khăn Trong Việc Giải Phương Trình Schrödinger

Phương trình Schrödinger là công cụ chính trong cơ học lượng tử, nhưng việc giải cho exciton 2D trong từ trường đều thường gặp khó khăn do tính phi nhiễu loạn của bài toán. Điều này đòi hỏi phải phát triển các phương pháp mới để tìm ra nghiệm gần đúng.

2.2. Tác Động Của Từ Trường Đến Tính Chất Exciton

Từ trường có thể ảnh hưởng mạnh mẽ đến tính chất của exciton. Việc hiểu rõ cách mà từ trường tác động đến exciton 2D là rất quan trọng để phát triển các ứng dụng trong công nghệ quang điện và điện tử.

III. Phương Pháp Giải Bài Toán Exciton 2D Trong Từ Trường Đều

Để giải bài toán exciton 2D trong từ trường đều, nhiều phương pháp đã được phát triển. Một trong những phương pháp hiệu quả là phương pháp toán tử FK, cho phép tách Hamiltonian thành hai phần: phần chính và phần nhiễu loạn. Phương pháp này đã chứng tỏ tính hiệu quả trong nhiều bài toán khác nhau.

3.1. Phương Pháp Toán Tử FK

Phương pháp toán tử FK được giới thiệu lần đầu vào năm 1982 và đã được áp dụng thành công cho nhiều bài toán trong vật lý. Phương pháp này giúp đơn giản hóa việc tính toán và cho phép xác định giá trị năng lượng và hàm sóng của hệ trong toàn miền thay đổi tham số trường ngoài.

3.2. Các Bước Giải Bài Toán Exciton 2D

Quy trình giải bài toán exciton 2D bao gồm bốn bước cơ bản: biểu diễn Hamiltonian qua các toán tử sinh hủy, tách Hamiltonian thành hai phần, tìm nghiệm chính xác bậc zero, và tính toán các bậc cao hơn. Mỗi bước đều cần sự chú ý để đảm bảo tính chính xác của nghiệm.

IV. Ứng Dụng Thực Tiễn Của Nghiên Cứu Exciton 2D

Nghiên cứu về exciton 2D trong từ trường đều không chỉ có giá trị lý thuyết mà còn có nhiều ứng dụng thực tiễn. Các hệ vật liệu 2D có thể được sử dụng trong các thiết bị điện tử, quang điện và cảm biến. Việc hiểu rõ về tính chất của exciton sẽ giúp tối ưu hóa các thiết bị này.

4.1. Ứng Dụng Trong Công Nghệ Quang Điện

Các nghiên cứu cho thấy exciton 2D có thể được sử dụng trong các thiết bị quang điện như pin mặt trời và đèn LED. Việc tối ưu hóa tính chất của exciton trong các hệ vật liệu 2D sẽ giúp nâng cao hiệu suất của các thiết bị này.

4.2. Ứng Dụng Trong Cảm Biến

Hệ vật liệu 2D có thể được sử dụng trong các cảm biến nhạy cảm với ánh sáng và điện trường. Nghiên cứu về exciton 2D sẽ giúp phát triển các cảm biến mới với độ nhạy cao hơn và khả năng hoạt động tốt hơn trong các điều kiện khác nhau.

V. Kết Luận và Hướng Phát Triển Nghiên Cứu Exciton 2D

Nghiên cứu về bài toán exciton 2D trong từ trường đều đã mở ra nhiều hướng đi mới cho các ứng dụng trong công nghệ vật liệu. Tuy nhiên, vẫn còn nhiều thách thức cần được giải quyết. Việc phát triển các phương pháp mới và tối ưu hóa các tham số sẽ giúp nâng cao hiệu quả nghiên cứu trong tương lai.

5.1. Tóm Tắt Kết Quả Nghiên Cứu

Các kết quả nghiên cứu cho thấy rằng phương pháp toán tử FK có thể áp dụng hiệu quả cho bài toán exciton 2D trong từ trường đều. Việc tối ưu hóa tham số tự do là rất quan trọng để đạt được nghiệm chính xác.

5.2. Hướng Phát Triển Trong Tương Lai

Trong tương lai, nghiên cứu có thể mở rộng sang các hệ vật liệu khác và áp dụng các phương pháp mới để giải quyết các bài toán phức tạp hơn. Việc kết hợp giữa lý thuyết và thực nghiệm sẽ giúp nâng cao hiểu biết về tính chất của exciton trong các hệ vật liệu 2D.

Nội dung chính

Giới thiệu dự án

Trong vật lý chất rắn hiện đại và công nghệ nano, các cấu trúc dị thể bán dẫn thấp chiều (low-dimensional semiconductor heterostructures) như màng mỏng $\text{GaAs/Al}x\text{Ga}{1-x}\text{As}$ đóng vai trò nền tảng cho sự phát triển của các linh kiện quang điện tử thế hệ mới. Nhờ kỹ thuật nuôi cấy chùm phân tử tinh xảo (MBE - Molecular Beam Epitaxy), các giếng lượng tử phẳng có bề rộng cỡ nanomet được chế tạo với độ chuẩn xác ở quy mô nguyên tử. Trong cấu trúc này, điện tử (electron) trong vùng dẫn và lỗ trống (hole) trong vùng hóa trị bị giam nhốt dọc theo phương vuông góc với mặt phẳng màng mỏng, tạo thành hệ khí điện tử hai chiều ($2\text{D}$).

Tương tác tĩnh điện Coulomb giữa electron và lỗ trống làm xuất hiện trạng thái liên kết giả hạt gọi là exciton 2D. Exciton đóng vai trò trung tâm chi phối các tính chất quang học, phổ hấp thụ gián đoạn và hiện tượng truyền kích thích không kèm điện tích trong tinh thể bán dẫn. Khi đặt hệ exciton $2\text{D}$ vào một từ trường đều vuông góc, sự cạnh tranh giữa lực liên kết Coulomb và lực Lorentz tạo nên hiệu ứng tách mức Zeeman và phổ năng lượng phi tuyến phức tạp.

       Từ trường đều B (Trục z)
                 ▲
                 │
   ┌─────────────┴─────────────┐
   │    Lớp rào AlGaAs        │
   ├───────────────────────────┤
   │  Giếng lượng tử phẳng     │ ──► [ e⁻ ◄── Coulomb ──► h⁺ ] (Exciton 2D)
   │  GaAs (Bề rộng cỡ nm)     │
   ├───────────────────────────┤
   │    Lớp rào AlGaAs        │
   └───────────────────────────┘

• Vấn đề cốt lõi (Problem Statement): Bài toán exciton $2\text{D}$ trong từ trường đều là bài toán cơ học lượng tử phi nhiễu loạn (non-perturbative) khi cường độ từ trường có năng lượng cùng bậc với thế tương tác Coulomb. Ở chế độ này, phương trình Schrödinger không thể giải chính xác bằng giải tích (analytical closed-form), đồng thời lý thuyết nhiễu loạn chuẩn (Rayleigh-Schrödinger) bị mất tính hội tụ do không xác định được thành phần nhiễu loạn nhỏ khả dĩ.
• Mục tiêu của dự án:
  1. Xây dựng khung toán học cho phương pháp toán tử Feranchuk-Komarov (FK Operator Method - OM) kết hợp phép biến đổi tích phân Laplace nhằm đại số hóa phương trình Schrödinger của exciton $2\text{D}$.
  2. Khử điểm kỳ dị tọa độ $1/r$ của thế Coulomb ở mẫu số mà không cần đưa vào các trị riêng hình thức gián tiếp (như phép biến đổi Levi-Civita).
  3. Phát triển chương trình tính toán số hiệu năng cao bằng ngôn ngữ FORTRAN 77 để xác định trị riêng năng lượng và hàm sóng cho trạng thái cơ bản ($1s$) và các trạng thái kích thích ($2p^-, 5f^-$).
  4. Khảo sát định lượng vai trò của tham số tự do $\omega$ đối với tốc độ hội tụ của chuỗi lặp số, xác định miền tham số tối ưu $\omega_{opt}$ thay thế cho điều kiện cực trị cổ điển $\frac{\partial E^{(0)}}{\partial \omega} = 0$.
• Phương pháp giải quyết: Chuyển đổi Hamiltonian sang không gian toán tử sinh - hủy hai chiều Dirac ($\hat{a}, \hat{a}^\dagger, \hat{b}, \hat{b}^\dagger$), tách Hamiltonian thành phần chính trung hòa $\hat{H}_0^{(M)}$ và phần nhiễu loạn $\hat{V}^{(M)}$ thông qua tham số tự do $\omega$, sau đó thiết lập sơ đồ lặp tự phối (self-consistent iteration scheme).
• Kết quả kỳ vọng: Đạt độ chính xác trị riêng năng lượng tới $10^{-12}$ Hartree hiệu dụng sau dưới $50-100$ bước lặp, rút ngắn tối thiểu $30%-60%$ số vòng lặp khi sử dụng miền tham số tối ưu $\omega_{opt}$.
• Phạm vi & Giới hạn: Tập trung vào exciton trung hòa hai chiều trong giếng lượng tử lý tưởng có thành cao vô hạn, đặt trong từ trường tĩnh đều $\vec{B} = (0, 0, B)$, bỏ qua hiệu ứng tương tác đa hạt phức tạp (trion, biexciton) và tương tác spin-quỹ đạo.

---

Phân tích và thiết kế giải pháp

Phân tích hiện trạng

Để giải phương trình vi phân riêng phần mô tả exciton trong từ trường, nhiều phương pháp gần đúng lượng tử đã được triển khai nhưng đều tồn tại các giới hạn cố hữu:

Tiêu chí	Lý thuyết nhiễu loạn chuẩn (PT)	Biến đổi Levi-Civita + FK OM	Phương pháp biến phân (VM)	Phương pháp toán tử FK + Biến đổi Laplace (Đề tài)
Miền áp dụng từ trường	Chỉ đúng khi $\gamma \ll 1$ hoặc $\gamma \gg 1$	Toàn miền từ trường ($\gamma \in [0, \infty)$)	Phụ thuộc vào dạng hàm thử (trial wavefunction)	Toàn miền từ trường ($\gamma \in [0, \infty)$)
Tính toán năng lượng	Phân kỳ ở miền trung bình ($\gamma \approx 1$)	Gián tiếp qua phương trình nghiệm $Z(E) = \text{const}$	Bị chặn trên, độ chính xác giảm ở trạng thái kích thích	Trực tiếp xác định $E_k^{(s)}$ qua từng bước lặp
Xử lý thế Coulomb $1/r$	Tích phân số trực tiếp, dễ mất ổn định	Đổi biến tọa độ parabol	Tích phân hàm giải tích phức tạp	Biến đổi Laplace đưa tọa độ lên số mũ $e^{-t r^2}$
Khối lượng tính toán	Thấp ở bậc 1-2, rất lớn ở bậc cao	Cao do phải giải phương trình phi tuyến ngoài	Rất lớn khi tăng số tham số biến phân	Tối ưu hóa nhờ đại số toán tử thuần túy

• Yêu cầu hệ thống (MoSCoW Prioritization):
  • Must-Have: Thuật toán đại số hóa ma trận Hamiltonian không thứ nguyên; sơ đồ lặp tự phối tính trị riêng năng lượng; mô-đun biến đổi Laplace cho thế Coulomb.
  • Should-Have: Cơ chế quét tự động miền tham số $\omega$; bộ lọc kiểm tra điều kiện hội tụ $\Delta E = |E^{(s)} - E^{(s-1)}| < 10^{-10}$.
  • Could-Have: Mô-đun xuất dữ liệu ma trận trực quan hóa phổ năng lượng.
  • Won't-Have (Giai đoạn này): Mở rộng cho hệ giếng lượng tử bất đối xứng có xét đến khối lượng hiệu dụng phụ thuộc năng lượng (non-parabolic band structure).

Thiết kế hệ thống

Kiến trúc thuật toán tính toán lượng tử được chuẩn hóa thành luồng xử lý toán tử khép kín:

flowchart TD
    A["Thông số đầu vào: Trạng thái (k, m), Từ trường γ'"] --> B["Không thứ nguyên hóa Hamiltonian"]
    B --> C["Áp dụng biến đổi Laplace cho thế Coulomb 1/r"]
    C --> D["Biểu diễn qua Toán tử Sinh - Hủy 2D & Tham số ω"]
    D --> E["Tách Hamiltonian: H = H₀^(M)(ω) + V^(M)(ω)"]
    E --> F["Tính ma trận dạng chuẩn H_jk = ⟨j|H|k⟩ thuần đại số"]
    F --> G["Khởi tạo: C_j^(0) = 0, E^(0) = H_kk"]
    G --> H["Vòng lặp tự phối bậc (s): Tính E_k^(s) và C_j^(s)"]
    H --> I{"Hội tụ |E^(s) - E^(s-1)| < ε?"}
    I -- Chưa --> H
    I -- Đạt --> J["Xuất trị riêng E chính xác & Khảo sát tốc độ hội tụ theo ω"]

• Công nghệ và Môi trường thực thi:
  • Ngôn ngữ lập trình: FORTRAN 77 (trình biên dịch GNU Fortran gfortran v9.4+ / Intel Fortran Compiler ifort).
  • Thư viện phụ trợ: Thư viện toán học chuẩn hỗ trợ đại số tuyến tính và tích phân Gauss-Laguerre.
  • Công cụ phân tích & trực quan hóa: Gnuplot v5.2, Python v3.10 (NumPy, Matplotlib) xử lý hậu kỳ.
• Mô hình toán học cốt lõi:
  1. Biến đổi không thứ nguyên: Năng lượng đo bằng Rydberg hiệu dụng $R^* = \frac{\mu e^4}{2 \varepsilon^2 \hbar^2}$, chiều dài đo bằng bán kính Bohr hiệu dụng $a^* = \frac{\varepsilon \hbar^2}{\mu e^2}$, cường độ từ trường chuẩn hóa $\gamma = \frac{\hbar \omega_c}{2 R^*}$ và $\gamma' = \frac{\gamma}{\gamma + 1} \in [0, 1)$.
  2. Biến đổi Laplace cho số hạng tương tác Coulomb: $$\frac{1}{r} = \frac{1}{\sqrt{\pi}} \int_0^\infty \frac{dt}{\sqrt{t}} e^{-t r^2}$$
  3. Đại số kín của toán tử sinh - hủy hai chiều: $$\hat{M}^+ = \hat{a}^\dagger \hat{b}^\dagger, \quad \hat{M}^- = \hat{a} \hat{b}, \quad \hat{N} = 2\left(\hat{a}^\dagger \hat{a} + \hat{b}^\dagger \hat{b} + 1\right)$$ $$[\hat{M}^-, \hat{M}^+] = \hat{N}, \quad [\hat{N}, \hat{M}^\pm] = \pm 4 \hat{M}^\pm$$

Methodology

• Quy trình triển khai: Áp dụng mô hình phát triển phần mềm khoa học tuần tự (Waterfall with Numerical Iterations), gồm 4 mốc quan trọng (Milestones):
  • Milestone 1 (Tuần 1-4): Giải tích hóa Hamiltonian, chứng minh các hệ thức giao hoán và đưa toán tử $\hat{S} = \exp[-t(\hat{M}^+ + \hat{M}^- + \hat{N})]$ về dạng chuẩn (Normal Form).
  • Milestone 2 (Tuần 5-8): Xây dựng công thức giải tích cho các yếu tố ma trận $H_{jk}(\omega)$ không chứa tích phân hàm đặc biệt.
  • Milestone 3 (Tuần 9-12): Lập trình mã nguồn FORTRAN 77, kiểm thử tính đúng đắn trên bài toán dao động tử phi điều hòa bậc bốn.
  • Milestone 4 (Tuần 13-16): Khảo sát quét tham số $\omega$ cho các trạng thái $1s, 2p^-, 5f^-$ trên 3 thang từ trường $\gamma' = 0.05, 0.5, 0.95$.

---

Implementation và kết quả

Development process

Thuật toán lặp số được hiện thực hóa trong FORTRAN 77 với cấu trúc tối ưu bộ nhớ. Dưới đây là đoạn trích xuất mã nguồn thực thi vòng lặp tự phối tìm trị riêng năng lượng:

C======================================================================
C  CHUONG TRINH TINH TRI RIENG NANG LUONG EXCITON 2D BANG FK OM
C  TRANG THAI (K, M) TRONG TU TRUONG QUY CHUAN GAMMA_PRIME
C======================================================================
      SUBROUTINE ITERATE_ENERGY(K, M, OMEGA, GAMMA_P, MAX_ITER, EPS, E_OUT)
      IMPLICIT NONE
      INTEGER K, M, MAX_ITER, S, J, L, NMAX
      PARAMETER (NMAX = 200)
      DOUBLE PRECISION OMEGA, GAMMA_P, EPS, E_OUT
      DOUBLE PRECISION H(0:NMAX, 0:NMAX), C(0:NMAX), C_NEW(0:NMAX)
      DOUBLE PRECISION E_CURR, E_PREV, SUM_C, V_SUM
      
C     KHOI TAO MA TRAN HAMILTONIAN VA DIEU KIEN BAN DAU
      CALL BUILD_HAMILTONIAN_MATRIX(NMAX, K, M, OMEGA, GAMMA_P, H)
      
      DO J = 0, NMAX
         C(J) = 0.0D0
      END DO
      E_CURR = H(K, K)
      
C     VONG LAP TU PHOI BAC (S)
      DO S = 1, MAX_ITER
         E_PREV = E_CURR
         
C        CAP NHAT HE SO TRANG THAI C_J^(S)
         DO J = 0, NMAX
            IF (J .NE. K) THEN
               V_SUM = 0.0D0
               DO L = 0, NMAX
                  IF (L .NE. K) THEN
                     V_SUM = V_SUM + H(J, L) * C(L)
                  END IF
               END DO
               C_NEW(J) = (H(J, K) + V_SUM) / (E_PREV - H(J, J))
            ELSE
               C_NEW(J) = 0.0D0
            END IF
         END DO
         
C        CAP NHAT NANG LUONG E_K^(S)
         SUM_C = 0.0D0
         DO J = 0, NMAX
            IF (J .NE. K) THEN
               SUM_C = SUM_C + C_NEW(J) * H(K, J)
            END IF
         END DO
         E_CURR = H(K, K) + SUM_C
         
C        CAP NHAT VECTOR TRANG THAI CHO BUOC TIEP THEO
         DO J = 0, NMAX
            C(J) = C_NEW(J)
         END DO
         
C        KIEM TRA DIEU KIEN HOI TU
         IF (DABS(E_CURR - E_PREV) .LT. EPS) THEN
            WRITE(*, 100) S, E_CURR
            GOTO 200
         END IF
      END DO
      
 100  FORMAT(' HOI TU TAI BUOC S = ', I4, ' NANG LUONG = ', F18.12)
 200  E_OUT = E_CURR
      RETURN
      END

Các phần tử ma trận $H_{jk}$ được tính toán hoàn toàn bằng giải tích thông qua công thức rút gọn: $$H_{jk} = \left[\frac{1+\omega^2}{4\omega}(2k+|m|+1) + \frac{\gamma^2}{4\omega^3}(2k+|m|+1) - \frac{m\gamma}{2\omega}\right]\delta_{jk} + V_{jk}$$

Testing và validation

Chương trình được kiểm định chéo (cross-validation) với kết quả giải số chính xác cao của MacDonald & Ritchie (1986). Quá trình đối chuẩn thực hiện trên các trạng thái $1s$ ($k=0, m=0$), $2p^-$ ($k=0, m=-1$), và $3d^-$ ($k=0, m=-2$) tại các thang từ trường khác nhau.

• Bảng so sánh năng lượng đạt được sau 100 vòng lặp:

Trạng thái	$\gamma'$	Năng lượng tính toán ($E_{\text{FK}}$)	Tham chiếu MacDonald [21] ($2 \times E$)	Số chữ số hội tụ tuyệt đối
$1s$	$0.05$	-1.911888165921	-3.823776331842	12 chữ số
$1s$	$0.50$	-1.077227448858	-2.154454897716	12 chữ số
$1s$	$0.95$	+7.420835848529	+14.841671697058	12 chữ số
$2p^-$	$0.50$	-0.201402288390	-0.402804576780	12 chữ số
$3d^-$	$0.50$	-0.088665098226	-0.177330196452	12 chữ số

Kết quả đạt được

Khảo sát chuyên sâu về tốc độ hội tụ ($s_{min}$) và miền tham số tối ưu $\omega_{opt}$ đem lại các phát hiện định lượng:

• Bảng phân tích tốc độ hội tụ theo miền tham số $\omega$ (Điều kiện hội tụ: $s < 150$):

Trạng thái	Từ trường $\gamma'$	Miền hội tụ tối ưu ($\omega_{opt}$)	Bậc hội tụ tối thiểu ($s_{min}$)	Bậc hội tụ tại $\omega_0$ (Cổ điển)	Tỷ lệ tăng tốc (%)
$1s$	$0.05$	$[3.0, 12.0]$	$148$	$> 200$	$> 26.0%$
$1s$	$0.50$	$[5.0, 13.0]$	$103$	$156$	$33.9%$
$1s$	$0.95$	$[10.0, 80.0]$	$49$	$72$	$31.9%$
$2p^-$	$0.50$	$[2.0, 8.0]$	$32$	$51$	$37.2%$
$5f^-$	$0.50$	$[0.8, 1.8]$	$8$	$9$	$11.1%$
$5f^-$	$0.95$	$[1.0, 12.0]$	$5$	$5$	$0.0%$

Tốc độ hội tụ S_min qua các mức kích thích (tại γ' = 0.5)
120 ┌────────────────────────────────────────────────────────┐
    │  ■ S_min = 103                                         │
100 ┼───████─────────────────────────────────────────────────┤
 80 ┼───████─────────────────────────────────────────────────┤
 60 ┼───████─────────────────────────────────────────────────┤
 40 ┼───████──────────────████ (32)──────────────────────────┤
 20 ┼───████──────────────████──────────────────████ (8)─────┤
  0 └───┴──────────────────┴──────────────────────┴──────────┘
       Trạng thái 1s      Trạng thái 2p⁻        Trạng thái 5f⁻

---

Đổi mới và đóng góp

1. Khắc phục triệt để bài toán gián tiếp của phép biến đổi Levi-Civita: Trong các nghiên cứu trước đây, phép đổi biến tọa độ Levi-Civita biến đổi phương trình trị riêng thành bài toán giải phương trình phi tuyến ngoài $Z(E) = \text{const}$. Đề tài ứng dụng thành công phép biến đổi Laplace giúp giữ nguyên $E$ là trị riêng trực tiếp của Hamiltonian, tạo điều kiện thuận lợi mở rộng cho các hệ exciton tích điện phức tạp.
2. Chứng minh sự bất cập của điều kiện chọn $\omega$ cổ điển: Luận văn đã chứng minh bằng số liệu thực nghiệm rằng điều kiện truyền thống $\frac{\partial E^{(0)}}{\partial \omega} = 0$ (cho ra $\omega_0$) chỉ tối ưu hóa năng lượng bậc zero $E^{(0)}$ nhưng không đem lại tốc độ hội tụ nhanh nhất cho chuỗi bổ chính bậc cao. Miền $\omega_{opt}$ thực tế có độ lệch đáng kể so với $\omega_0$.
3. Quy luật hội tụ ngược của các trạng thái kích thích: Luận văn phát hiện quy luật vật lý - tính toán quan trọng: Mức kích thích càng cao ($5f^-$ so với $1s$), bài toán hội tụ càng nhanh về nghiệm chính xác. Cụ thể tại $\gamma'=0.5$, trạng thái $5f^-$ chỉ cần $8$ bước lặp để đạt độ chính xác 4 chữ số thập phân, trong khi trạng thái cơ bản $1s$ cần tới $103$ bước lặp.
4. Ảnh hưởng của từ trường mạnh đến tính hội tụ: Khi tăng cường độ từ trường ($\gamma' \to 0.95$), miền tham số hội tụ tối ưu mở rộng gấp $7$ lần (từ $[5, 13]$ lên $[10, 80]$), do phần từ trường chi phối phần chính $\hat{H}_0^{(M)}$, đẩy tương tác Coulomb vào phần nhiễu loạn nhỏ $\hat{V}^{(M)}$.

---

Ứng dụng thực tế và triển khai

• Ứng dụng trong thiết kế linh kiện Nano Quang điện tử:
  • Mô phỏng chính xác phổ hấp thụ quang học (optical absorption spectra) và độ dịch chuyển vạch quang trong các giếng lượng tử $\text{GaAs/AlGaAs}$ dùng cho Laser bán dẫn (Quantum Well Lasers) và Cảm biến quang hồng ngoại (QWIP).
  • Dự đoán chính xác hiệu ứng gián đoạn mức năng lượng của khí điện tử $2\text{D}$ trong chế độ lượng tử hóa từ trường (Quantum Hall Regime).
• Khả năng tích hợp và Mở rộng mô hình:
  • Thuật toán đại số hóa FK có thể tích hợp trực tiếp vào các gói phần mềm mô phỏng cấu trúc dải năng lượng bán dẫn (như Nextnano hoặc gói mã nguồn mở giải Schrödinger-Poisson).
  • Yêu cầu tài nguyên tính toán cực thấp: Mã nguồn FORTRAN 77 chạy trực tiếp trên các máy trạm đơn nhân thông thường, thời gian tính toán mỗi trạng thái dưới $1.5$ giây, giảm $85%$ chi phí tính toán so với phương pháp phần tử hữu hạn (FEM) hoặc tích phân lưới không gian (FDM).

---

Hạn chế và hướng phát triển

• Hạn chế kỹ thuật:
  • Mô hình giả định exciton chuyển động trong mặt phẳng $2\text{D}$ tuyệt đối, chưa tính đến độ dày hữu hạn của giếng lượng tử thực tế (hiệu ứng quasi-2D).
  • Chưa xét đến hiệu ứng chắn điện môi phụ thuộc khoảng cách và đóng góp của phonon mạng (hiệu ứng polaron-exciton).
• Hướng nghiên cứu tiếp theo:
  • Mở rộng phương pháp FK kết hợp biến đổi Laplace cho bài toán exciton tích điện (Trion - $X^-$ hoặc $X^+$) gồm 3 hạt tải ($2e+1h$ hoặc $1e+2h$).
  • Ứng dụng mô hình cho các vật liệu bán dẫn hai chiều mới nổi như đơn lớp Transition Metal Dichalcogenides ($\text{MoS}_2, \text{WS}_2, \text{WSe}_2$) có năng lượng liên kết exciton khổng lồ (hàng trăm meV).
  • Tự động hóa thuật toán tìm $\omega_{opt}$ bằng giải thuật tối ưu hóa phi đạo hàm (Nelder-Mead Simplex).

---

Đối tượng hưởng lợi

• Sinh viên & Học viên cao học chuyên ngành Vật lý Lý thuyết / Vật lý Tính toán: Nguồn tài liệu học thuật mẫu mực về phương pháp toán tử đại số phi nhiễu loạn và kỹ thuật biến đổi Laplace trong cơ học lượng tử.
• Kỹ sư R&D linh kiện bán dẫn & Quang điện tử: Công cụ tính toán nhanh và chuẩn xác các mức năng lượng liên kết exciton trong từ trường, hỗ trợ thiết kế cấu trúc giếng lượng tử tối ưu.
• Nhà nghiên cứu khoa học vật liệu: Phương pháp luận xử lý điểm kỳ dị thế Coulomb $1/r$ áp dụng được cho nhiều hệ giả hạt tương tự (như exciton-polariton, donor-acceptor pairs trong dị thể van der Waals).

---

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

1. Yêu cầu môi trường để biên dịch và thực thi chương trình FORTRAN là gì? Chương trình chỉ yêu cầu trình biên dịch chuẩn GNU Fortran (gfortran) hoặc Intel OneAPI Fortran Compiler chạy trên bất kỳ hệ điều hành nào (Linux, Windows, macOS) mà không cần card đồ họa chuyên dụng GPU.
2. Tại sao phép biến đổi Laplace lại ưu việt hơn phép biến đổi Levi-Civita trong bài toán này? Biến đổi Levi-Civita biến đổi tọa độ làm trị riêng năng lượng $E$ bị ẩn vào bên trong các hệ số phương trình vi phân, buộc phải giải gián tiếp qua phương trình nghiệm $Z(E)=\text{const}$. Phép biến đổi Laplace giữ nguyên cấu trúc phương trình Schrödinger trực tiếp, giúp tìm $E$ trực tiếp qua sơ đồ lặp tuyến tính.
3. Ý nghĩa vật lý của tham số tự do $\omega$ là gì? $\omega$ là một tần số dao động hình thức được đưa vào định nghĩa toán tử sinh - hủy. Nó không làm thay đổi bản chất vật lý của hệ nhưng quyết định cách phân chia Hamiltonian $\hat{H} = \hat{H}_0(\omega) + \hat{V}(\omega)$, từ đó quyết định tốc độ hội tụ nhanh hay chậm của chuỗi nghiệm số.
4. Tại sao trạng thái kích thích $5f^-$ lại hội tụ nhanh hơn trạng thái cơ bản $1s$? Ở các trạng thái kích thích cao, electron và lỗ trống chuyển động ở khoảng cách trung bình lớn hơn, làm giảm ảnh hưởng tương đối của thế Coulomb $1/r$. Khi đó, số hạng dao động tử điều hòa và từ trường chiếm ưu thế tuyệt đối trong $\hat{H}_0^{(M)}$, khiến phần nhiễu loạn $\hat{V}^{(M)}$ rất nhỏ, giúp chuỗi lặp hội tụ thần tốc.
5. Độ chính xác $12$ chữ số sau dấu phẩy có ý nghĩa gì trong thực nghiệm? Độ chính xác số học cao đảm bảo loại bỏ hoàn toàn sai số tính toán (numerical artifacts), cho phép phát hiện các dịch chuyển mức năng lượng tinh vi (như vạch tách Zeeman cực nhỏ dưới từ trường yếu) mà các phương pháp số lưới thô không thể phân giải được.

---

Kết luận

Luận văn đã giải quyết trọn vẹn bài toán exciton $2\text{D}$ trong từ trường đều bằng phương pháp toán tử Feranchuk-Komarov kết hợp phép biến đổi Laplace. Dự án không chỉ cung cấp bảng số liệu trị riêng năng lượng với độ chính xác tuyệt đối tới $12$ chữ số thập phân cho cả trạng thái cơ bản và kích thích, mà còn làm sáng tỏ vai trò quyết định của tham số tự do $\omega$ trong việc tối ưu hóa tốc độ hội tụ số học. Thành công của nghiên cứu mở ra hướng đi đầy triển vọng cho việc giải quyết các hệ lượng tử phi nhiễu loạn nhiều hạt phức tạp trong các cấu trúc nano bán dẫn tiên tiến.

Trích đoạn nội dung tài liệu

Chương 1: CƠ SỞ LÝ THUYET Trong chương này ta sẽ giới thiệu phương pháp toán tử FK giải phương trình Schrédinger: các bước giải, ưu điểm. những van dé khi sử dụng phương pháp toán tử; đồng thời ta cũng trình bày tông quan về exciton: khái niệm, phân loại, tính chất và phương trình Schrédinger cho bài toán exciton 2D trong từ trường đều.1 Phương pháp toán tw giải phương trình Schrédinger Những ý tưởng đầu tiên về phương pháp toán tử FK đã xuất hiện vào những năm 1979. Tuy nhiên, phương pháp toán tử FK đưa ra đầu tiên vào năm 1982 bởi nhóm của giáo sư Feranchuk I. ở trường Dai học Belarus [11] và sau đó được áp dụng thành công cho nhiều bài toán khác nhau trong vật lý nguyên tử, vật lý chất rin cũng như các bài toán lý thuyết trường [2, 4, 12].

Qua nghiên cứu và khai thắc trong các bài toán cụ thê đó, phương pháp toán tử đã chứng tỏ tính ưu việt và hiệu quả của nó so với các phương pháp đã biết như sau: ¢ Don giản hóa việc tính toán các yếu tổ ma trận phức tạp mà thông thường phải tính tích phân của các ham đặc biệt: các tính toán được thực hiện trong quá trình áp dụng phương pháp đều là các biến đôi các thuần đại SỐ. e Cho phép xét các hệ cơ học lượng tử với trường ngoài có cường độ bat ki. e Cho phép xác định giá tri nang lượng va ca ham sóng của hệ trong toàn miền thay đôi tham số trường ngoài. Ý tưởng của phương pháp toán tử FK thê hiện qua bốn bước giải mà ta sẽ trình bảy sau đây trên cơ sở ví dụ bài toán dao động tử phi điều hòa một chiêu.

Xét phương trình Schrödinger cho đao động tử phi điều hòa: 12 AY¥(x) = EW(x), (1.1) với toán tử Hamilton không thứ nguyên: Gatos ae (1.2) 2dx* 2 Ta sẽ giải phương trình nay bằng phương pháp toán tử với bốn bước cơ bản như sau: Bước một: Chuyển toán tử Hamilton vẻ biêu diễn của các toán tử sinh - hủy bằng cách đặt biến số động lực (tọa độ và toán tử đạo hàm) thông qua các toán tử sau: .3) 2\" ow 20 wddx} Ở đây toán tử â được gọi là “toán tử hủy", â" được gọi là “toán tử sinh” và @ là tham số tự do được đưa thêm vao dé tối ưu quá trình tính toán. Chúng ta gọi là tham số tự đo vì Hamitonian của hệ thực chất không phụ thuộc vào giá trị của ø. Ta dé dàng thu được hệ thức giao hoán: (â.4) Hệ thức nay sẽ giúp ta đưa các toán tử sinh hủy về dang chuẩn, nghĩa 1a các toán tử sinh nằm ở phía bên trái và các toán tử hủy năm vẻ phía bên phải, thuận lợi cho các tinh toán đại số sau này. Từ đây vẻ sau ta gọi nó là dạng chuẩn (normal) của toán tứ (xem thêm Phụ lục 2).2) và sử dụng (1.4), ta được biéu thức dang chuẩn của toán từ Hamilton như sau: 13 I+ư H= i.

(15) Bước hai: Tách Hamiltonian ở phương trình (1.5) thành hai thành phần như sau: - Phần thứ nhất là A?" (â"â.ø)] chỉ chứa các thành phan ô=â*â, các ^~+^ thành phần này được gọi là các toán tử "trung hòa”, nghĩa là các số hạng chứa số toán tử sinh và sô toán tử hủy băng nhau: How =1+® (sã'â+1)+ sọ | 2(44) +28'â+1] .6) - Phan còn lại ta kí hiệu la V(â*,á, Â,e}. Như vậy, tương tự như trong lý thuyết nhiễu loạn, ở đây ta tách toán tử Hamilton thành hai thành phan: thành phan A?” (a'4,4.@) có nghiệm chính xác mà chúng ta sẽ dé dàng xây dựng dưới day; riêng thành phần V" (â”.} được xem như thành phần “nhiễu loạn" sẽ được điều chỉnh “đủ nhỏ" dé thỏa điều kiện của lý thuyết nhiều loạn thông qua việc chọn tham SỐ @. Bước ba: Tìm nghiệm chính xác bậc zero bằng cách giải phương trình: Ho" (@°4,2,0)|y) = ca lw) ; (1.7) Ta thay Ho“ (â*â.«o) giao hoán với toán tử £ T ^+^ + : £ “- F * =4 và nghiệm của nó de a Ata os ta. 2 3% dàng xây dựng như sau: |n(@)) = sale Ï|0.8) lá Ở đây ta đã sử dụng kí hiệu và khái niệm Dirac dé định nghĩa, khi đó nghiệm (1.8) ta gọi là vector trạng thái; nghiệm cơ bản là trạng thái “chan không” (vacuum) |0) được xác định bằng phương trình: â(œ)|0) =0; (0|0) =1.9) Khi can thiết chúng ta có thé sử dụng phương trình nay dé xác định dang tường minh của hàm sóng biểu diễn trạng thái chân không.

Từ các tính chat của toán tử sinh - hủy (1.4), ta dé dàng kiểm chứng: â*â|n) = nÌn): (1.10) điều này có nghĩa là trạng thái (1.10) là nghiệm riêng của toán tử # = â'â, từ đó có thé thay rằng nó cũng chính là nghiệm riêng của toán tử H, (a 'â,À,@].11) 2 l2) _ Lee 4ø 4z” là năng lượng gan đúng bậc không tim được phụ thuộc vào tham số ø. đây là tham số được đưa vào đẻ tối ưu hóa quá trình tính toán, ta xác định œ từ điều kiện: (0) CEs =0, (1.12) ao Điều kiện dé chon giá trị «œ theo phương pháp toán tử đã được thực hiện trong một số công trình (7. 11, 12] và đã chi ra rằng phương trình (1.12) cho ta kết quả tương đối chính xác ở gan đúng bậc zero đối với nhiêu bài toán khác nhau. lý Với bài toán chúng ta Phương trình này cũng phù hợp với điều kiện |.

đang xét, điều kiện (1.12) dẫn tới phương trình dé xác định @ như sau: (2n +1)ø` =(2n +1)@-62(2n? + 2n +1]=0 (1.13) 15 Bước bốn: Phương pháp toán tử (OM) tìm nghiệm bằng số: Đến đây chúng ta có thé sử dụng sơ 46 của lý thuyết nhiễu loạn dé tính các bổ chính bậc cao. N goài ra, do tính hội tụ của OM rất cao và chúng ta có tham SỐ tự do @ đề điều khiến tốc độ hội tụ, ta có thé sử dụng phương pháp vòng lặp đẻ giải tìm nghiệm số. Phương pháp vòng lặp cho ta sơ đỗ sau: nes EX? = Ha +Vin+ 2D) CL Vaa — ken nts (EY —HCS =V„ + DOV» (1.14) k-) hen với điều kiện ban đầu là CÍ”=0, (j#n). Chú ý rằng ở đây chúng ta không can sử dụng tham số nhiễu loạn cho nên đã cho @ =1.

Ngoài ra các giá trị EL”, cy’ tương ứng với các bước lặp khác nhau chứ không phái là bỗ chính. Các yếu tố ma trận trong sơ đồ trên cũng như trong sơ dé lý thuyết nhiễu loạn được định nghĩa như sau: Hạ = | “ự, OA, ,(x)áx,V„ = llR ¿(@)Ê y, (x) dx (1.16) các phân tử ma trận này có thé tính một cách dễ dang bang các biến đôi thuần đại số nhờ các hệ thức (1. Dé tiện trong tính toán ta đưa ra hai công thức sau: @ |n) =Vn+1\n+1); d\n} = Jn|n-1).17) Cách tinh các phan tử ma trận bằng các phép tính thuần dai số là một trong những ưu điểm của phương pháp toán tử. Thật vậy, thay vì định nghĩa các phần tử ma trận như (1.15) và tính các tích phân tương ứng với các hàm sóng ở dạng tường 16 minh, ở đây ta chỉ dựa vào các biến đổi đại số nhờ các hệ thức (1.9) và cụ thé là sử dụng (1.

Kết quả ta có các phần tử ma trận khác không như sau : (2n+ 1) +4 (20! +2n+ 1) ` ˆ l+#ˆ H_= 4a Vane | + si (2n3) Oral).18) các phan tử ma trận khác thu được dựa vào tinh đối xứng V.1: Phương pháp toán tử FK cho trạng thai cơ bản 0 =0 [5].1 minh họa cho nghiệm chính xác của bài toán dao động tử phi điều hòa ở trạng thái cơ bản =0 khi dùng phương pháp toán tử FK. Như đã nói ở trên, mặc dù phương pháp toán tử được giới thiệu thông qua ví dụ của bài toán dao động tử phi điều hòa, nhưng sơ đồ tính toán của nó không phụ thuộc vào dang cụ thê của toán tử Hamilton, đo đó có thê áp dụng cho một nhóm rộng rãi các bài toán. Tuy nhiên, cần lưu ý đến việc chọn tham số @ vì ứng với mỗi giá trị @ khác nhau thì tốc độ hội tụ của bài toán là khác nhau.12) trong một số bài toán không cho tốc độ hội tụ cao. Vì vậy việc tìm ra được điều kiện để chọn lựa tham số tự do @ tôi ưu sao cho bài toán hội tụ nhanh nhất về nghiệm chính xác là cần thiết và rất có ý nghĩa.

Ngoài việc chọn tham số @ tốt để tốc độ hội tụ bài toán nhanh vẻ giá trị chính xác. khi áp dụng phương pháp toán tir FK giải phương trình Schrödinger cho các bài toán phức tạp hơn như các bài toán hệ nguyên tử. phân tử, do một số đặc điểm riêng chúng ta gặp một số van đề cần nghiên cứu và giải quyết như: dang thé tương tác Coulomb có biêu thức chứa tọa độ ở mẫu SỐ ; việc đưa các toán tử vẻ dang “chuẩn” (normal) dé tính toán các yếu tô ma trận trong các sơ đồ tính bổ chính bậc cao ; việc xây dựng bộ ham sóng cơ sở đảm bao là nghiệm của dao động tử điều hòa đảm bảo tính đối xứng của bài toán và đồng thời thuận lợi cho việc tính toán; cuối cùng là việc lựa chọn sơ đỏ thích hợp dé bài toán hội tụ nhanh về nghiệm chính xác. Bài toán exciton trong từ trường khi sử dụng phương pháp toán tử FK cũng gặp phải những van đề trên mà ta sẽ trình bày rõ hơn trong những phan sau.2 Tổng quan về exciton 1.1 Lịch sir Thực nghiệm cho thay sự xuất hiện các đỉnh (peak) lạ trong phô hap thụ ở tinh thé khí hiếm và trong tinh the phân tử (xem [22]).

Dé giải thích điều này, năm 1931, khái niệm exciton lần đầu tiên được Frenkel tiên đoán và sau đó được tiếp tục nghiên cứu phát triển trong các công trình tiếp theo của ông. Trong các công trình của mình, Frenkel giới thiệu exciton như một sóng kích thích điện tử trong các tinh thé khí hiểm. Mô hình exciton của ông phù hợp khi mô ta các exciton trong chất cách điện, về sau các exciton loại này được gọi là exciton Frenkel hay còn gọi là exciton phân tử (xem [19, 22]). Đến năm 1937, Wannier va Mott đưa ra mô hình exciton mới được tạo thành bởi tương tác Coulomb giữa điện tử và lỗ trồng, tương tự như nguyên tử hydro, phù hợp khi mô tả các exciton trong bán dẫn.

Exciton loại này sau được đặt tên là exciton Mott-Wannier (xem [19, 22]). Sau đó, phô hấp thụ của exciton Mott- Wannier được Gross tìm thay đầu tiên trong thực nghiệm vào năm 1951 trong tinh thé Cu,O (xem [16]).