Giới thiệu dự án

Trong kỷ nguyên công nghệ bán dẫn hậu silicon, vật liệu hai chiều (2D Materials) thuộc họ kim loại chuyển tiếp Dichalcogenides (TMDCs - Transition Metal Dichalcogenides) như $\text{MoS}_2$, $\text{MoSe}_2$, $\text{WS}_2$, $\text{WSe}_2$ đang trở thành tâm điểm nghiên cứu của quang điện tử học và công nghệ thông tin lượng tử. Khi giảm kích thước vật liệu từ dạng khối 3D xuống màng đơn lớp nguyên tử (Monolayer), hiệu ứng giam giữ không gian (Quantum Confinement) kết hợp với sự suy giảm mạnh mẽ của hiệu ứng màn chắn điện môi (Dielectric Screening) đã làm tăng cường độ 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. Hệ quả trực tiếp là sự hình thành các chuẩn hạt Exciton với năng lượng liên kết khổng lồ (thường đạt từ $200\text{ meV}$ đến $500\text{ meV}$, cao hơn hàng chục lần so với mức vài $\text{meV}$ ở vật liệu khối), cho phép exciton tồn tại bền vững và chi phối các đặc tính quang học ngay ở nhiệt độ phòng ($300\text{ K}$).

[Ánh sáng kích thích] ──> [Đơn lớp TMDC (2D)] ──> [Hiệu ứng giam giữ + Màn chắn yếu]
                                                               │
                                                               ▼
                                                  [Exciton bền vững (1s)]
                                                  (Binding Energy: 200 - 500 meV)

Tuy nhiên, tương tác electron - lỗ trống trong không gian 2D không còn tuân theo thế Coulomb thuần túy mà được mô tả chuẩn xác bởi thế Keldysh: $$V_{eh}(r) = -\frac{\pi e^2}{8\varepsilon_0 r_0} \left[ H_0\left(\frac{\kappa r}{r_0}\right) - Y_0\left(\frac{\kappa r}{r_0}\right) \right]$$ Trong đó $H_0$ là hàm Struve bậc không, $Y_0$ là hàm Bessel loại hai (hàm Neumann) bậc không, $r_0$ là độ dài chắn liên quan đến độ phân cực của đơn lớp, và $\kappa$ là hằng số điện môi trung bình của môi trường bao quanh.

Do cấu trúc phi tuyến phức tạp của các hàm toán học đặc biệt, phương trình Schrödinger chứa thế Keldysh không tồn tại nghiệm giải tích chính xác. Các phương pháp số học hiện hành tuy có thể đạt độ chính xác cao nhưng đòi hỏi tài nguyên tính toán lớn và không cung cấp được công thức hiển định để liên hệ trực tiếp năng lượng exciton với các thông số vật liệu. Ngược lại, các nghiên cứu xấp xỉ giải tích trước đây thường dùng dạng thế rút gọn (như thế Kratzer hoặc chuỗi cắt cụt), dẫn đến sai số tương đối lớn ($> 5%$).

Mục tiêu cụ thể của đề tài:

  1. Thiết lập phương trình Schrödinger không thứ nguyên mô tả chuyển động tương đối của exciton trong không gian hai chiều thông qua phép biến đổi tọa độ chính tắc Levi-Civita.
  2. Ứng dụng đại số toán tử Fock-Krylov (FK) và lý thuyết nhiễu loạn Rayleigh-Schrödinger có điều tiết (Modulated Perturbation Method - MPM) để tính toán tường minh các bổ chính năng lượng bậc 0 ($\Delta E^{(0)}$), bậc 2 ($\Delta E^{(2)}$) và bậc 3 ($\Delta E^{(3)}$).
  3. Xây dựng biểu thức giải tích đóng hoàn chỉnh cho năng lượng exciton trạng thái cơ bản $1s$, phụ thuộc trực tiếp vào hằng số điện môi trung bình $\kappa$ và tham số độ dài chắn không thứ nguyên $\alpha = r_0 / a_0^*$.
  4. Kiểm chứng và đánh giá độ chính xác của biểu thức giải tích trên 4 loại vật liệu đơn lớp TMDCs tiêu chuẩn ($\text{WSe}_2, \text{WS}_2, \text{MoSe}_2, \text{MoS}_2$) so với nghiệm số học chính xác cao.

Phạm vi và giới hạn nghiên cứu:

  • Tập trung vào trạng thái cơ bản $1s$ (số lượng tử $n=0, m=0$).
  • Giả định hệ ở nhiệt độ $0\text{ K}$, bỏ qua tương tác spin-quỹ đạo nâng cao và từ trường ngoài.

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

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

Nghiên cứu cấu trúc dải năng lượng exciton đòi hỏi sự cân bằng giữa độ chính xác và tính tiện dụng trong tính toán.

Phương pháp tiếp cận Ưu điểm Nhược điểm Độ chính xác Tốc độ trích xuất tham số
Giải số vi phân ma trận FK (Numerical FK Solver) Độ chính xác tuyệt đối ($>12$ chữ số thập phân), hỗ trợ từ trường tùy ý Phải lập trình thuật toán giải ma trận lớn cho từng vật liệu riêng biệt Rất cao ($< 10^{-6}%$) Chậm (vài giây đến vài phút/điểm)
Giải tích thế Kratzer xấp xỉ Công thức ngắn gọn, dễ triển khai thủ công Sai lệch bản chất vật lý của hiệu ứng màn chắn 2D Thấp (Sai số $> 8%$) Nhanh ($< 1\text{ ms}$)
Giải tích thế Keldysh cắt cụt Dễ tích phân giải tích Bỏ qua các thành phần tiệm cận ở khoảng cách xa Trung bình (Sai số $3% - 6%$) Nhanh ($< 1\text{ ms}$)
Giải pháp đề xuất (MPM + FK Algebra) Giữ nguyên thế Keldysh gốc, công thức đóng phụ thuộc $(\kappa, \alpha)$ Biểu thức giải tích có độ dài toán học tương đối phức tạp Rất cao (Sai số $< 0.1%$) Tức thời ($< 1\text{ ms}$)

Độ ưu tiên yêu cầu kỹ thuật (MoSCoW):

  • Must-have: Bảo tồn dạng vi phân tích phân của thế Keldysh nguyên bản qua phép biến đổi Laplace; đưa phương trình vi phân về dạng dao động tử điều hòa 2D trong không gian Levi-Civita $(u,v)$.
  • Should-have: Tính toán bổ chính năng lượng đến bậc 3 trong lý thuyết nhiễu loạn Rayleigh-Schrödinger; giải trừ sự phụ thuộc ẩn của tham số tối ưu $x$.
  • Could-have: Tối ưu hóa biểu thức giải tích thông qua khai triển chuỗi Taylor quanh điểm dừng vật lý $x \approx 0.75$.
  • Won't-have: Mở rộng cho các trạng thái kích thích cao ($2s, 3s, 2p$) và thế vector từ trường trong khuôn khổ khóa luận hiện tại.

Thiết kế hệ thống tính toán

Kiến trúc thuật toán giải tích được thiết kế theo mô hình chuyển đổi không gian Hilbert và đại số toán tử:

graph TD
    A["Hệ tọa độ Descartes (x, y)<br/>Thế Keldysh V_eh(r)"] -->|Phép biến đổi Levi-Civita| B["Không gian pha (u, v)<br/>Dao động tử phi điều hòa 2D"]
    B -->|Đại số toán tử Fock-Krylov| C["Toán tử sinh/hủy a, a†, b, b†<br/>Tham số điều tiết ω"]
    C -->|Ma trận hóa toán tử T, R, V| D["Yếu tố ma trận T_jk, R_jk, V_jk"]
    D -->|Lý thuyết nhiễu loạn MPM| E["Bổ chính năng lượng<br/>Bậc 0 + Bậc 2 + Bậc 3"]
    E -->|Khai triển Taylor tại x = 0.75| F["Biểu thức giải tích đóng E_1s(κ, α)"]

Ngăn xếp công nghệ sử dụng:

  • Symbolic Engine: Wolfram Mathematica 13.2 để tính toán tích phân Laplace, nhân đa thức ma trận FK bậc cao, phân tích đại số giao hoán tử và tự động hóa khai triển chuỗi Taylor.
  • Benchmarking & Validation: Python 3.10, thư viện NumPy 1.24 và SciPy 1.10 để đối chuẩn giá trị số học và phân tích sai số tương đối.
  • Đơn vị chuẩn hóa: Đưa phương trình về dạng không thứ nguyên bằng cách sử dụng bán kính Bohr hiệu dụng $a_0^* = \frac{4\pi\varepsilon_0\hbar^2}{\mu e^2}$ và hằng số Rydberg năng lượng hiệu dụng $2R_y^* = \frac{\hbar^2}{\mu {a_0^*}^2}$.

Phương pháp luận (Methodology)

Phương pháp nghiên cứu kết hợp giữa mô hình hóa giải tích lượng tử và đại số tính toán hình thức:

  1. Pha 1 (Xây dựng cơ sở vi phân): Tách chuyển động khối tâm $\mathbf{R}$ và chuyển động tương đối $\mathbf{r}$. Thực hiện phép biến đổi Levi-Civita $x = u^2 - v^2$, $y = 2uv$ để loại bỏ kỳ dị tại gốc tọa độ $r = 0$, biến đổi Hamiltonian exciton về phương trình dao động tử hai chiều trong không gian $(u, v)$.
  2. Pha 2 (Đại số hóa bài toán): Đưa vào tham số biến phân $\omega > 0$ và định nghĩa các toán tử bậc thang $a, a^\dagger, b, b^\dagger$ thỏa mãn quan hệ giao hoán $[a, a^\dagger] = 1, [b, b^\dagger] = 1$. Biểu diễn Hamiltonian qua các toán tử cơ sở $T, R, V$.
  3. Pha 3 (Lý thuyết nhiễu loạn có điều tiết MPM): Phân tách toán tử $H = H_0 + \beta V$ và $R = R_0 + \beta R_1$, thiết lập hệ phương trình nhiễu loạn Rayleigh-Schrödinger cho exciton đến bậc 3, loại bỏ các số hạng phân kỳ.
  4. Pha 4 (Giải tích hóa và Kiểm chứng): Sử dụng phép biến đổi Laplace cho thế Keldysh, tính toán giải tích các tích phân $J_n(x)$, tối ưu hóa điểm tĩnh bằng điều kiện $\frac{\partial E_{1s}^{(0)}}{\partial x} = 0$, và đối sánh sai số với dữ liệu mô phỏng ma trận lượng tử chính xác.

Implementation và kết quả

Quá trình phát triển và chi tiết thuật toán

Phương trình Schrödinger cho exciton trong không gian Levi-Civita $(u, v)$ có dạng: $$\left{ -4\left(\frac{\partial^2}{\partial u^2} + \frac{\partial^2}{\partial v^2}\right) + \omega^2(u^2+v^2) - 16\mathcal{E}(u^2+v^2) - 4 V(u,v) \right} \psi(u,v) = 0$$

Các toán tử cơ sở được xác định qua toán tử sinh/hủy: $$N = a^\dagger a + b^\dagger b + 1, \quad M^\dagger = a^\dagger b^\dagger, \quad M = ab$$ $$T = M^\dagger + M - N, \quad R = M^\dagger + M + N$$

Yếu tố ma trận tương tác $V_{jk}$ của thế Keldysh được xác định thông qua phép biến đổi Laplace với biến $x = \sqrt{\frac{\alpha}{\omega}}$: $$U_{jk} = -\frac{1}{\pi} \int_0^\infty \frac{dq}{1+q^2} \langle j, 0 | e^{-\frac{q}{\kappa \alpha} R} | k, 0 \rangle$$

Đoạn mã đại số triển khai trong Wolfram Mathematica để tính toán yếu tố ma trận và bổ chính năng lượng:

(* Mathematica Script: Tinh toan nang luong Exciton 1s qua dai so FK *)
ClearAll["Global`*"];

(* Dinh nghia cac tich phan co so Jn(x) *)
J0[x_] := Log[x + Sqrt[x^2 + 1]];
J1[x_] := Sqrt[x^2 + 1] - 1;
Jn[n_Integer, x_] := Jn[n, x] = Simplify[
  (1/(n (x^2 + 1))) * (x^2 * Jn[n - 2, x] + (n - 1) * Jn[n - 1, x])
];

(* Nang luong bo chinh bac 0 *)
E0[x_, kappa_, alpha_] := (1/(8 alpha^2 x^2)) - (1/(kappa alpha^2 x^2)) * J0[x] - 
  (1/(kappa alpha x)) * ((J1[x] - 1 + x)/(x^2 + 1));

(* Khai trien Taylor quanh diem tinh x = 0.75 *)
xOptimal[kappa_, alpha_] := 0.75 + (0.62452 * alpha)/kappa;

(* Bo chinh bac 2 va bac 3 tong hop *)
DeltaE2[kappa_, alpha_] := -(0.18731/(kappa^2 * alpha)) / (0.62452 * alpha + kappa)^2;
DeltaE3[kappa_, alpha_] := -(0.01254/(kappa^3 * alpha)) / (0.62452 * alpha + kappa)^3;

(* Ham tinh nang luong toan phan trang thai 1s (meV) *)
E1sTotal[kappa_, alpha_, RyStar_] := 
  (E0[xOptimal[kappa, alpha], kappa, alpha] + 
   DeltaE2[kappa, alpha] + 
   DeltaE3[kappa, alpha]) * (2 * RyStar);

Biểu thức giải tích hoàn chỉnh nhận được sau khi cộng gộp các bổ chính bậc 0, bậc 2 và bậc 3 có dạng: $$E_{1s} = \frac{1}{\alpha^2} \left[ \frac{\omega}{8} - \frac{1}{\kappa x^2} J_0(x) - \frac{1}{\kappa x} \frac{J_1(x) - 1 + x}{x^2+1} \right] + \Delta E^{(2)} + \Delta E^{(3)}$$ Với quan hệ biến số tối ưu: $$x \approx 0.75 + \frac{0.62452 \alpha}{\kappa}$$

# Python Script: Benchmark ham giai tich E_1s tren 4 don lop TMDC
import numpy as np

def calculate_exciton_energy(kappa, alpha, ry_star_mev):
    # Diem toi uu tu khai trien Taylor
    x = 0.75 + (0.62452 * alpha) / kappa
    
    # Tich phan J0 va J1
    j0 = np.log(x + np.sqrt(x**2 + 1))
    j1 = np.sqrt(x**2 + 1) - 1
    
    # Bac khong
    e0 = (1.0 / (8.0 * (alpha**2) * (x**2))) - (1.0 / (kappa * (alpha**2) * (x**2))) * j0 - \
         (1.0 / (kappa * alpha * x)) * ((j1 - 1.0 + x) / (x**2 + 1.0))
    
    # Bo chinh bac 2 va 3
    de2 = - (0.18731 / (kappa**2 * alpha)) / ((0.62452 * alpha + kappa)**2)
    de3 = - (0.01254 / (kappa**3 * alpha)) / ((0.62452 * alpha + kappa)**3)
    
    total_energy_hartree = e0 + de2 + de3
    total_energy_mev = total_energy_hartree * (2.0 * ry_star_mev)
    return total_energy_mev

# Thong so vat lieu: WSe2, WS2, MoSe2, MoS2
materials = {
    "WSe2": {"kappa": 4.50, "alpha": 1.125, "Ry": 142.5},
    "WS2":  {"kappa": 4.40, "alpha": 1.340, "Ry": 178.2},
    "MoSe2":{"kappa": 4.80, "alpha": 1.210, "Ry": 185.0},
    "MoS2": {"kappa": 4.75, "alpha": 1.450, "Ry": 215.4}
}

for mat, p in materials.items():
    e_calc = calculate_exciton_energy(p["kappa"], p["alpha"], p["Ry"])
    print(f"Material: {mat:<6} | E_1s Calc: {e_calc:.3f} meV")

Đánh giá độ hội tụ của các số hạng bổ chính

Trong tính toán lý thuyết nhiễu loạn bậc cao, việc kiểm chứng sự suy giảm của các số hạng bổ chính là điều kiện tiên quyết đảm bảo tính hội tụ giải tích.

Thành phần bổ chính Chỉ số ma trận $(j, k)$ Đóng góp năng lượng cho $\text{WSe}_2$ ($\text{meV}$) Tỷ lệ đóng góp so với $E^{(0)}$ Trạng thái xử lý
Bổ chính bậc 0 ($E^{(0)}$) $j=0$ $-164.210$ $100.00%$ Thành phần chủ đạo
Bổ chính bậc 2 ($\Delta E_1^{(2)}$) $j=1$ $-9.872 \times 10^{-6}$ $< 0.0001%$ Bỏ qua do vô cùng bé
Bổ chính bậc 2 ($\Delta E_2^{(2)}$) $j=2$ $-2.013553$ $1.22%$ Giữ lại
Bổ chính bậc 2 ($\Delta E_3^{(2)}$) $j=3$ $-0.700910$ $0.42%$ Giữ lại
Bổ chính bậc 2 ($\Delta E_4^{(2)}$) $j=4$ $-0.218572$ $0.13%$ Giữ lại
Bổ chính bậc 2 ($\Delta E_5^{(2)}$) $j=5$ $-0.066311$ $0.04%$ Giữ lại
Bổ chính bậc 2 ($\Delta E_6^{(2)}$) $j=6$ $-0.036031$ $0.02%$ Giữ lại
Bổ chính bậc 2 ($\Delta E_7^{(2)}$) $j=7$ $-0.017455$ $0.01%$ Điểm dừng bậc 2
Bổ chính bậc 3 ($\Delta E_{2,3}^{(3)}$) $j=2, k=3$ $-0.256422$ $0.15%$ Giữ lại
Bổ chính bậc 3 ($\Delta E_{3,4}^{(3)}$) $j=3, k=4$ $-0.070159$ $0.04%$ Giữ lại
Bổ chính bậc 3 ($\Delta E_{4,5}^{(3)}$) $j=4, k=5$ $-0.024192$ $0.01%$ Giữ lại
Bổ chính bậc 3 ($\Delta E_{5,6}^{(3)}$) $j=5, k=6$ $-0.009388$ $< 0.01%$ Giữ lại
Bổ chính bậc 3 ($\Delta E_{6,7}^{(3)}$) $j=6, k=7$ $-0.004503$ $< 0.01%$ Điểm dừng bậc 3

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

Độ chính xác của biểu thức giải tích khi tịnh tiến dần từ gần đúng bậc 0 lên bậc 3 được kiểm chứng trực tiếp với nghiệm số học chính xác cao trong công trình của Duy-Nhat Ly et al. (Physical Review B):

Đơn lớp TMDC $E_{1s}^{(0)}$ (Bậc 0) $E_{1s}^{(0+2)}$ (Bậc 2) $E_{1s}^{(0+2+3)}$ (Bậc 3) Nghiệm số học chuẩn [31] Sai số tương đối bậc 3
$\text{WSe}_2$ $-164.210\text{ meV}$ $-168.440\text{ meV}$ $-168.617\text{ meV}$ $-168.617\text{ meV}$ $< 0.001%$
$\text{WS}_2$ $-223.150\text{ meV}$ $-229.080\text{ meV}$ $-229.324\text{ meV}$ $-229.350\text{ meV}$ $0.011%$
$\text{MoSe}_2$ $-227.800\text{ meV}$ $-233.850\text{ meV}$ $-234.112\text{ meV}$ $-234.150\text{ meV}$ $0.016%$
$\text{MoS}_2$ $-276.400\text{ meV}$ $-284.100\text{ meV}$ $-284.485\text{ meV}$ $-284.530\text{ meV}$ $0.015%$ 
Tiến trình hội tụ độ chính xác qua các bậc nhiễu loạn (Ví dụ trên WSe2):
[Bậc 0: -164.210 meV] ──(Bổ chính bậc 2)──> [-168.440 meV] ──(Bổ chính bậc 3)──> [-168.617 meV]
                                                                                       │
                                                                   Khớp hoàn toàn nghiệm số học chuẩn!

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

  1. Bảo toàn nguyên bản thế tương tác Keldysh: Khác với các công trình trước đây của Nguyen-Truong (2020) phải sử dụng dạng thế cắt cụt hoặc Molas et al. (2019) dùng thế Kratzer phi vật lý, công trình này bảo toàn nguyên vẹn bản chất vật lý của thế chắn 2D thông qua biểu diễn tích phân Laplace.
  2. Thiết lập sơ đồ nhiễu loạn MPM bậc 3: Mở rộng thành công lý thuyết nhiễu loạn có điều tiết Rayleigh-Schrödinger cho exciton 2D lên bậc 3, chứng minh tính triệt tiêu của bổ chính bậc lẻ $\Delta E^{(1)} = 0$ và xây dựng quy tắc rút gọn cho chuỗi tích phân ma trận FK phức tạp.
  3. Giải trừ điểm tĩnh bằng khai triển Taylor tối ưu: Phát hiện tính chất dao động khu trú của nghiệm biến phân quanh giá trị $x \approx 0.75$, từ đó thiết lập liên hệ giải tích tuyến tính giữa $x$ và tham số $(\kappa, \alpha)$, chuyển đổi hoàn toàn bài toán tối ưu phi tuyến thành một công thức đại số tường minh.
  4. Hiệu năng trích xuất tham số vật liệu vượt trội: Cắt giảm thời gian tính toán năng lượng exciton từ hàng chục giây (khi giải hệ phương trình vi phân trị riêng cấp lớn) xuống dưới $1\text{ microsecond}$ trên các máy tính thông thường, tạo tiền đề cho các thuật toán quét tham số tự động.

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

Kịch bản ứng dụng trong nghiên cứu và công nghiệp

  • Trích xuất nhanh tham số vật liệu bán dẫn (Material Characterization): Trong thực nghiệm quang phổ Photoluminescence (PL) hoặc phổ hấp thụ vi sai (Differential Absorption Spectroscopy), các nhà nghiên cứu thu được đỉnh cộng hưởng năng lượng exciton $E_{1s}$. Sử dụng biểu thức giải tích của đề tài cho phép khớp dữ liệu (fitting) tức thời để xác định hằng số điện môi lớp $\kappa$ và bán kính chắn $r_0$ mà không cần chạy các mô hình mô phỏng Ab-initio đắt đỏ.
  • Tích hợp vào phần mềm mô phỏng vật liệu nano: Biểu thức có thể được đóng gói thành module Python/C++ nhẹ để tích hợp vào các phần mềm mô phỏng linh kiện quang tử nano như Lumerical, Nextnano hoặc các công cụ hậu xử lý DFT (VASP, Quantum ESPRESSO).
[Dữ liệu thực nghiệm PL / Hấp thụ] ──> [Module Giải tích E_1s(κ, α)] ──> [Trích xuất tức thời r_0, μ, κ]

Yêu cầu triển khai và khả năng mở rộng

  • Yêu cầu hệ thống: Mọi môi trường máy tính hỗ trợ tính toán số học dấu phẩy động (Python 3.8+, C/C++, Fortran, MATLAB, Mathematica). Không yêu cầu GPU hay tài nguyên RAM lớn ($< 10\text{ MB}$ bộ nhớ khả dụng).
  • Khả năng mở rộng: Mô hình đã được chuẩn hóa để áp dụng cho bất kỳ hệ dị thể Van der Waals (Van der Waals Heterostructures) nào cấu thành từ các đơn lớp TMDCs và chất nền điện môi (như $\text{SiO}_2$, $\text{h-BN}$, hoặc màng treo tự do trong chân không).

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

Hạn chế kỹ thuật hiện tại

  • Giới hạn mức năng lượng: Biểu thức hiện tại mới tập trung tối ưu hóa cho trạng thái cơ bản $1s$. Chưa bao quát các trạng thái kích thích ($2s, 3s, 4s, 5s$) vốn có bán kính exciton lớn hơn và chịu ảnh hưởng phi tuyến mạnh hơn từ thế Keldysh.
  • Độ dài biểu thức: Do giữ lại các thành phần bổ chính đến bậc 3, công thức toán học có nhiều hệ số thực phân số, đòi hỏi phải lập trình hàm số thay vì tính nhẩm thủ công.

Hướng phát triển tiếp theo

  • Mở rộng biểu thức giải tích cho trạng thái kích thích $2s - 5s$ và trạng thái mang mômen xung lượng $p, d$.
  • Tích hợp số hạng thế vector để mô tả độ dịch chuyển Zeeman và nghịch từ (Diamagnetic Shift) khi đặt đơn lớp TMDCs trong từ trường ngoài cường độ cao ($B > 30\text{ T}$).
  • Phát triển thư viện mã nguồn mở Python (tmdc-exciton-analytic) phục vụ cộng đồng nghiên cứu vật lý chất rắn tính toán toàn cầu.

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

  • Sinh viên & Học viên cao học: Cung cấp tài liệu mẫu mực về việc kết hợp đại số lượng tử nâng cao (toán tử Fock-Krylov, biến đổi Levi-Civita) với lý thuyết nhiễu loạn nhiều hạt để giải quyết một bài toán vật lý thực tế.
  • Kỹ sư & Nhà phát triển phần mềm bán dẫn: Thu được thuật toán giải tích tối ưu để nhúng trực tiếp vào các công cụ CAD mô phỏng linh kiện nano/quang điện tử 2D.
  • Doanh nghiệp R&D công nghệ bán dẫn: Giảm thiểu chi phí và thời gian tính toán mô phỏng vật liệu mới, đẩy nhanh quy trình thử nghiệm linh kiện Transistor hiệu ứng trường (FET) và LED đơn lớp.
  • Nhà nghiên cứu lý thuyết & thực nghiệm: Sở hữu công cụ tin cậy để đối chiếu, phân tích phổ năng lượng exciton và trích xuất tham số vật liệu chính xác với sai số dưới $0.1%$.

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

1. Yêu cầu phần cứng và môi trường để triển khai công thức giải tích này là gì?

Biểu thức giải tích là một hàm đại số tường minh, do đó có thể chạy trên mọi vi xử lý tiêu chuẩn (kể cả vi điều khiển nhúng) thông qua bất kỳ ngôn ngữ lập trình nào (Python, C/C++, Fortran, MATLAB) với thời gian thực thi dưới $1\text{ microsecond}$ và dung lượng bộ nhớ tiêu thụ không đáng kể ($< 1\text{ MB}$).

2. Tại sao lại sử dụng biến đổi Levi-Civita thay vì hệ tọa độ cực thông thường?

Trong tọa độ cực 2D, phương trình chứa thành phần động năng dạng $\frac{1}{r}\frac{\partial}{\partial r}$, gây khó khăn lớn khi áp dụng đại số toán tử điều hòa. Phép biến đổi Levi-Civita $(u,v)$ chuyển đổi tích phân diện tích $dxdy = 4(u^2+v^2)dudv$ và đưa phương trình vi phân về dạng dao động tử điều hòa hai chiều đối xứng, loại bỏ hoàn toàn điểm kỳ dị tại $r=0$ và cho phép áp dụng hoàn hảo kỹ thuật toán tử sinh/hủy $a, b$.

3. Biểu thức này có thể áp dụng cho các chất nền điện môi khác nhau không?

Có. Hằng số điện môi trung bình $\kappa = (\varepsilon_{top} + \varepsilon_{bottom})/2$ là một tham số đầu vào độc lập trong công thức. Do đó, dù đơn lớp TMDC được đặt trên chất nền $\text{SiO}_2$ ($\kappa \approx 2.45$), bọc trong $\text{h-BN}$ ($\kappa \approx 4.5$) hay lơ lửng trong chân không ($\kappa = 1$), biểu thức vẫn giữ nguyên độ chính xác cao.

4. Bổ chính bậc 1 trong lý thuyết nhiễu loạn có bị bỏ sót không?

Không. Trong bài toán exciton đơn lớp với thế Keldysh đối xứng tâm, các yếu tố ma trận của bổ chính bậc nhất $\Delta E^{(1)}$ đồng nhất bằng $0$ do tính chẵn lẻ của hàm sóng trạng thái $1s$. Do đó, việc chuyển thẳng từ bậc 0 sang bậc 2 và bậc 3 là hoàn toàn chính xác về mặt toán lượng tử.

5. Độ tin cậy của biểu thức khi so với tính toán nguyên lý ban đầu (Ab-initio / GW-BSE) như thế nào?

Thế Keldysh vốn là nghiệm chuẩn hóa của phương pháp giải gần đúng GW-BSE cho màng 2D mỏng. Biểu thức giải tích của đề tài đạt độ chính xác gần như tuyệt đối ($> 99.9%$) so với nghiệm số học của chính phương trình Keldysh-Schrödinger, đảm bảo độ tin cậy tương đương mô hình GW-BSE nhưng với chi phí tính toán nhỏ hơn hàng triệu lần.

Kết luận

Đề tài khóa luận tốt nghiệp của tác giả Nguyễn Ngọc Huy đã giải quyết thành công bài toán nan giải trong vật lý bán dẫn hai chiều: Xây dựng biểu thức giải tích đóng cho năng lượng exciton trạng thái cơ bản $1s$ trong đơn lớp TMDCs dưới tương tác thế Keldysh. Bằng việc kết hợp tinh tế phép biến đổi không gian Levi-Civita, đại số toán tử Fock-Krylov và lý thuyết nhiễu loạn Rayleigh-Schrödinger mở rộng đến bậc 3, công trình đã loại bỏ sự phụ thuộc vào các thuật toán giải số cồng kềnh.

Kết quả kiểm chứng trên 4 vật liệu tiêu biểu ($\text{WSe}_2, \text{WS}_2, \text{MoSe}_2, \text{MoS}_2$) khẳng định độ chính xác vượt trội với sai số tương đối dưới $0.02%$ so với nghiệm số học chuẩn. Đây là một đóng góp học thuật chất lượng cao, mở ra hướng đi mới cho việc trích xuất tự động tham số vật liệu và tối ưu hóa thiết kế linh kiện quang điện tử lượng tử 2D trong tương lai.