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

Trong cơ học lượng tử, phương trình Schrödinger đóng vai trò là phương trình động lực học nền tảng mô tả trạng thái và hành vi của các hệ vi mô. Tuy nhiên, phương trình này là một phương trình vi phân tuyến tính chứa các đạo hàm riêng phần với hệ số biến đổi, khiến số lượng bài toán có lời giải giải tích chính xác là rất hữu hạn (chẳng hạn như bài toán nguyên tử hydro hoặc dao động tử điều hòa đơn giản). Phần lớn các bài toán thực tế đều đòi hỏi các phương pháp giải gần đúng như lý thuyết nhiễu loạn hoặc phương pháp biến phân. Dẫu vậy, lý thuyết nhiễu loạn đòi hỏi số hạng nhiễu loạn phải rất nhỏ so với phần chính, trong khi phương pháp biến phân lại phức tạp hóa bài toán do sự xuất hiện của quá nhiều tham số và khó khăn khi áp dụng cho các trạng thái kích thích bậc cao.

Sự phát triển của vật lý vật liệu chứng minh rằng khi số chiều không gian giảm từ 3D xuống 2D, 1D hoặc 0D, tính chất vật lý của vật liệu có sự thay đổi rõ rệt. Tiếp sau graphene – loại vật liệu 2D có hạn chế lớn là độ rộng vùng cấm bằng không – dòng vật liệu bán dẫn đơn lớp TMDs (Transition Metal Dichalcogenides) nổi lên như một đối tượng nghiên cứu tiềm năng nhờ sở hữu vùng cấm trực tiếp nằm trong khoảng từ 1 eV đến 3 eV. Trong cấu trúc đơn lớp TMDs, tương tác giữa electron và lỗ trống bị giới hạn trong không gian hai chiều, tạo nên giả hạt exciton có mật độ cao, thời gian sống lâu và năng lượng liên kết lớn hơn nhiều so với vật liệu khối ba chiều.

Do năng lượng của các trạng thái kích thích exciton rất khó đo đạc chính xác trong thực nghiệm ở điều kiện thông thường, người ta thường đặt hệ vào một từ trường ngoài không đổi (từ 0 đến 60 Tesla) nhằm làm tăng năng lượng liên kết và làm nổi bật các vạch phổ quang học. Vấn đề đặt ra là cần tìm một phương pháp toán học hiệu quả, giải quyết nhanh và tốn ít tài nguyên tính toán để giải phương trình Schrödinger mô tả exciton 2D trong từ trường. Phương pháp toán tử FK (Feranchuk – Komarov, đề xuất năm 1982) cho phép đại số hóa phương trình vi phân thông qua các toán tử sinh – hủy, đưa bài toán phức tạp về dạng đại số tuyến tính với độ chính xác cao.

Khóa luận xác định các mục tiêu và nhiệm vụ nghiên cứu cụ thể:

  1. Mục tiêu nghiên cứu: Áp dụng phương pháp toán tử FK để tìm nghiệm số chính xác cho bài toán exciton hai chiều trong đơn lớp TMDs đặt trong từ trường đều.
  2. Nhiệm vụ nghiên cứu:
    • Tìm hiểu tổng quan về đối tượng exciton và đặc trưng vật liệu đơn lớp TMDs.
    • Thiết lập bài toán exciton 2D trong từ trường đều với thế tương tác Keldysh.
    • Đại số hóa và giải phương trình Schrödinger cho exciton 2D bằng phương pháp toán tử FK kết hợp lý thuyết nhiễu loạn.
    • So sánh kết quả tính toán năng lượng liên kết với dữ liệu thực nghiệm đã công bố để đánh giá độ chính xác.

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

  • Đối tượng nghiên cứu: Exciton Mott-Wannier trung hòa (gồm một electron và một lỗ trống) trong vật liệu bán dẫn đơn lớp WSe₂.
  • Phạm vi nghiên cứu: Khảo sát chuyển động tương đối hai chiều trong mặt phẳng Oxy dưới tác dụng của từ trường đều vuông góc từ 0 đến 60 Tesla, tính toán năng lượng liên kết ở các mức trạng thái lượng tử 1s, 2s và 3s.

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

Khóa luận được xây dựng dựa trên sự kết hợp giữa mô hình exciton hai chiều, hàm thế Keldysh, lý thuyết nhiễu loạn không suy biến và phương pháp toán tử FK:

Khung lý thuyết và mô hình tương tác

  • Mô hình exciton Mott-Wannier: Cặp electron - lỗ trống liên kết yếu trong chất bán dẫn, có bán kính tương tác lớn (khoảng 100 Å) và năng lượng liên kết nhỏ (khoảng vài trăm meV).
  • Thế tương tác Keldysh ($V_{ep}(r)$): Thay vì sử dụng thế tĩnh điện Coulomb thuần túy ($1/r$), tác giả sử dụng thế Keldysh để mô tả chính xác hiệu ứng chắn điện môi suy giảm trong vật liệu 2D: $$V_{ep}(r) = -\frac{\pi e^2}{2\kappa r_0} \left[ H_0\left(\frac{r}{r_0}\right) - Y_0\left(\frac{r}{r_0}\right) \right]$$ Trong đó $H_0$ và $Y_0$ lần lượt là hàm Struve và hàm Bessel loại hai; $r_0 = 2\pi\chi_{2D}$ là độ dài chắn phụ thuộc độ phân cực 2D; $\kappa = (\epsilon_1 + \epsilon_2)/2$ là hằng số điện môi trung bình của môi trường xung quanh.
  • Biến đổi Fourier của thế Keldysh: Để thuận tiện cho việc chuyển đổi toán tử, thế Keldysh được biểu diễn qua tích phân Fourier trong không gian xung lượng: $$V_{ep}(r) = -\frac{2}{\pi} \int \frac{dq_x dq_y}{q(1 + \alpha q)} e^{i(q_x x + q_y y)}$$ với $\alpha = r_0/a_B^*$ là tham số độ dài màn chắn không thứ nguyên.

Phương pháp nghiên cứu

  • Hệ đơn vị nguyên tử Hartree: Biến đổi phương trình Schrödinger về dạng không thứ nguyên với bán kính Bohr hiệu dụng $a_B^* = \frac{\hbar^2 \kappa}{\mu e^2}$, đơn vị năng lượng $2R^* = \frac{\mu e^4}{\hbar^2 \kappa^2}$, khối lượng rút gọn $\mu = \frac{m_e^* m_h^}{m_e^ + m_h^}$, và cường độ từ trường không thứ nguyên $\gamma = \frac{e\hbar B}{2\mu R^}$.

  • Phương pháp toán tử FK (Feranchuk – Komarov):

    1. Đưa các biến tọa độ và vi phân về toán tử sinh – hủy Dirac ($\hat{a}_x, \hat{a}_x^+, \hat{a}_y, \hat{a}_y^+$) gắn với tham số biến phân $\omega$.
    2. Chéo hóa toán tử hình chiếu moment động lượng $L_z = \hbar(\hat{a}^+\hat{a} - \hat{b}^+\hat{b})$ bằng cách định nghĩa các toán tử quay $\hat{a}, \hat{a}^+, \hat{b}, \hat{b}^+$.
    3. Chuẩn hóa hàm toán tử thế năng dạng số mũ $e^{i\vec{q}\vec{r}}$ về dạng chuẩn tắc (normal order: toán tử sinh bên trái, trung hòa ở giữa, hủy bên phải) thông qua giải hệ phương trình vi phân cấp một.
    4. Xây dựng bộ hàm cơ sở trực chuẩn $|n_1, n_2\rangle = \frac{1}{\sqrt{n_1! n_2!}} (\hat{a}^+)^{n_1} (\hat{b}^+)^{n_2} |0,0\rangle$.

  • Lý thuyết nhiễu loạn không suy biến: Khai triển năng lượng và hàm sóng theo các bậc bổ chính ($E = E^{(0)} + \Delta E^{(1)} + \Delta E^{(2)}$), trong đó năng lượng bậc không được xác định qua trị trung bình $E^{(0)} = \langle n_1 n_2 | \hat{H} | n_1 n_2 \rangle$.

  • Điều kiện cực trị xác định tham số $\omega$: Do tham số $\omega$ không có trong Hamiltonian gốc, $\omega$ được xác định chính xác tại mỗi mức từ trường thông qua nghiệm của phương trình vi phân: $$\frac{\partial E^{(0)}_{n_1 n_2}(\omega)}{\partial \omega} = 0$$

  • Nguồn dữ liệu đối chiếu: Số liệu thực nghiệm quang từ (magneto-optics) đo năng lượng exciton trong đơn lớp WSe₂ dưới từ trường 0–60 Tesla trích xuất từ công trình của Crooker et al. (Phys. Rev. Lett., 2018).

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

Chương 1. Tổng quan và cơ sở lý thuyết

Chương 1 hệ thống hóa lịch sử phát hiện và bản chất vật lý của exciton từ năm 1907 (Becquerel), qua các mốc nghiên cứu của Frenkel (1931), Mott - Wannier (1937), Gross (1951) trên tinh thể Cu₂O, đến phát hiện exciton tích điện (trion) của Lampert (1958). Tác giả phân loại exciton theo hai tiêu chí: điện tích (trung hòa, trion âm, trion dương) và vật liệu (Frenkel với bán kính tương tác nhỏ cỡ 10 Å, năng lượng liên kết cỡ eV; Mott-Wannier với bán kính lớn cỡ 100 Å, năng lượng liên kết cỡ 100 meV).

Tiếp đó, chương luận giải cấu trúc của họ vật liệu TMDs đơn lớp dạng $MX_2$ ($M = \text{W, Mo, Ti...}$; $X = \text{S, Se, Te}$), liên kết cộng hóa trị nội lớp và Van der Waals giữa các lớp, nhấn mạnh sự chuyển đổi sang vùng cấm trực tiếp $1-3\text{ eV}$ khi giảm chiều dày. Các tiềm năng ứng dụng của TMDs được liệt kê chi tiết: transistor và linh kiện FET dưới 2 nm, cảm biến quang (photodetector), thiết bị ngắt quang (photoswitching), điện xúc tác sản xuất hydro sạch bằng MoS₂ thay thế bạch kim (Pt), và cảm biến sinh học định vị DNA. Cuối chương, tác giả trình bày cơ sở toán học 4 bước của phương pháp toán tử FK và biểu thức bổ chính năng lượng bậc không, bậc một, bậc hai của lý thuyết nhiễu loạn không suy biến.

Chương 2. Bài toán exciton trung hòa trong đơn lớp TMDs

Chương 2 tập trung vào các bước biến đổi giải tích và đại số hóa phương trình Schrödinger hai chiều:

  • Thiết lập toán tử Hamiltonian toàn phần của exciton dưới từ trường đều $B$: $$\hat{H} = -\left(\frac{\partial^2}{\partial x^2} + \frac{\partial^2}{\partial y^2}\right) + V_{ep}(r) + \frac{\gamma^2}{8}(x^2 + y^2) + \frac{\gamma}{2}\hat{L}_z$$
  • Chuyển đổi các biến động lực học sang hệ toán tử sinh – hủy $\hat{a}, \hat{b}$ và toán tử phụ $\hat{A}, \hat{B}$: $$\hat{A} = \frac{i}{\sqrt{2\omega}}(q_x - i q_y)\hat{a} + \frac{i}{\sqrt{2\omega}}(q_x + i q_y)\hat{b}^+$$
  • Đưa toán tử thế năng Keldysh về dạng chuẩn tắc bằng cách giải hệ phương trình vi phân cho hàm e mũ: $$V_{ep}(r) = -\frac{2}{\pi} \int \frac{dq_x dq_y}{q(1 + \alpha q)} e^{-\eta^2} e^{\eta \hat{A}^+} e^{\eta \hat{A}} e^{\eta \hat{B}^+} e^{\eta \hat{B}}$$ với $\eta = q/\sqrt{2\omega}$.
  • Xây dựng biểu thức tường minh cho các yếu tố ma trận của toán tử động năng - từ trường $\hat{T}{k_1 k_2; n_1 n_2}$ và ma trận thế năng Keldysh $\langle k_1 k_2 | V{ep} | n_1 n_2 \rangle$ qua chuỗi hàm giải tích và tích phân truy hồi $J(p, \alpha, c)$.
  • Thiết lập công thức năng lượng bậc không $E^{(0)}_{n_1 n_2} = \langle n_1 n_2 | \hat{H} | n_1 n_2 \rangle$.

Chương 3. Kết quả tính toán

Chương 3 trình bày các kết quả số học thu được khi giải phương trình tìm cực trị $\frac{\partial E^{(0)}}{\partial \omega} = 0$ để xác định $\omega_0$ tối ưu cho từng trạng thái lượng tử của đơn lớp $\text{WSe}_2$ với các tham số vật liệu: khối lượng hiệu dụng $m_e^* = 0.36 m_0$, $m_h^* = 0.45 m_0$.

  • Trường hợp không có từ trường ($B = 0\text{ T}$): Tính toán năng lượng liên kết exciton $E_b$ cho 3 trạng thái $1s$ ($n_1=n_2=0$), $2s$, và $3s$. Giá trị trạng thái cơ bản $1s$ đạt xấp xỉ $1.72311\text{ eV}$, sai lệch dưới $1%$ so với dữ liệu thực nghiệm.
  • Trường hợp có từ trường đều ($B = 0 - 60\text{ T}$): Năng lượng liên kết của các mức $1s, 2s, 3s$ đều tăng dần theo cường độ từ trường, phản ánh rõ nét hiệu ứng nén hàm sóng exciton trong không gian 2D.
  • Đánh giá độ lệch phổ: Đồ thị so sánh giữa năng lượng lý thuyết bậc không và thực nghiệm từ công trình [11] cho thấy đường lý thuyết bám sát các điểm đo thực nghiệm trên toàn dải từ trường từ 0 đến 60 Tesla, với độ sai lệch dao động trong khoảng từ $0%$ đến $2%$.

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

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

Các kết quả số học của khóa luận chứng minh tính hiệu quả của việc kết hợp phương pháp FK với thế tương tác Keldysh trên đơn lớp $\text{WSe}_2$:

Trạng thái lượng tử Năng lượng liên kết $E_b$ ($B = 0\text{ T}$) Dải từ trường khảo sát Độ sai lệch so với thực nghiệm [11]
Trạng thái 1s ($n_1=0, n_2=0$) $1.72311\text{ eV}$ $0 - 60\text{ T}$ $< 1%$ (tại $B=0$) / $0 - 2%$ (toàn dải)
Trạng thái 2s Xác định theo $\omega_{0(2s)}$ $0 - 60\text{ T}$ $0 - 2%$
Trạng thái 3s Xác định theo $\omega_{0(3s)}$ $0 - 60\text{ T}$ $0 - 2%$

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

  • Hoàn thiện quy trình đại số hóa thế Keldysh: Chuyển đổi thành công biểu thức thế phi Coulomb phức tạp chứa hàm Bessel/Struve về dạng tích các toán tử sinh – hủy chuẩn tắc thông qua phép biến đổi Fourier và chuỗi Lie đại số.
  • Tối ưu hóa bài toán giải vi phân lượng tử: Giảm thiểu độ phức tạp tính toán so với việc giải trực tiếp phương trình vi phân đạo hàm riêng trên lưới tọa độ hoặc phương pháp biến phân nhiều tham số, chỉ cần giải phương trình cực trị một biến $\omega$ để đạt nghiệm số chính xác cao.
  • Kiểm chứng độ tin cậy của phương pháp FK: Minh chứng rằng ngay ở mức gần đúng bậc không $E^{(0)}$, phương pháp toán tử FK đã đem lại kết quả có độ chính xác cao (sai số $0-2%$) khi so sánh với phổ quang học thực nghiệm dưới từ trường mạnh lên đến 60 T.

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

  • Hạn chế:

    • Nghiên cứu mới chỉ dừng lại ở việc tính toán năng lượng ở phép gần đúng bậc không ($E^{(0)}$), chưa tính toán bổ chính năng lượng bậc một ($\Delta E^{(1)}$) và bậc hai ($\Delta E^{(2)}$) theo lý thuyết nhiễu loạn.
    • Chưa mở rộng khảo sát cho các trạng thái kích thích có moment động lượng quỹ đạo khác không ($L_z \neq 0$, ví dụ các trạng thái $2p, 3d, 3p$) và các hệ exciton tích điện (trion).
    • Khóa luận chưa tiến hành đánh giá chi tiết tốc độ hội tụ số khi tăng số lượng hàm sóng cơ sở trong ma trận biểu diễn.

  • Hướng phát triển:

    • Tiếp tục tính toán các bổ chính bậc cao hơn ($\Delta E^{(1)}, \Delta E^{(2)}$) để triệt tiêu sai số $1-2%$ hiện tại, nâng cao độ chính xác tuyệt đối cho các mức $1s, 2s, 3s$.
    • Khảo sát sự hội tụ của năng lượng khi mở rộng không gian Hilbert của bộ hàm cơ sở.
    • Xây dựng công thức giải tích gần đúng cho năng lượng exciton dựa trên các bậc bổ chính thấp nhằm cung cấp biểu thức tường minh phục vụ các nghiên cứu lý thuyết và thực nghiệm tiếp theo.

Giá trị tham khảo

  • Đối tượng sử dụng:
    • Sinh viên, học viên cao học, nghiên cứu sinh chuyên ngành Vật lý lý thuyết, Vật lý chất rắn, Sư phạm Vật lý và Khoa học vật liệu.
    • Các nhà nghiên cứu quan tâm đến phương pháp giải gần đúng phương trình vi phân lượng tử và mô phỏng tính chất quang học của vật liệu bán dẫn hai chiều (2D materials).
  • Phần nội dung giá trị nhất:
    • Chương 2 và hệ thống 7 Phụ lục giải tích: Cung cấp chi tiết từng bước biến đổi đại số đưa hàm e mũ toán tử thế Keldysh về dạng chuẩn tắc ($S_1, S_2$), cách tính tác dụng của toán tử lên hàm sóng cơ sở 2D $|n_1, n_2\rangle$, và các công thức truy hồi giải tích tích phân $J(p, \alpha, c)$.
    • Quy trình xác định tham số tối ưu $\omega$: Là tài liệu tham khảo mẫu mực về kỹ thuật tìm điểm dừng năng lượng $\partial E^{(0)}/\partial \omega = 0$ để tối ưu hóa hàm sóng dao động tử điều hòa trong bài toán phi điều hòa.

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

1. Tại sao thế tương tác Keldysh lại được sử dụng thay thế cho thế Coulomb thông thường trong bài toán exciton đơn lớp TMDs?

Trong vật liệu khối 3D, tương tác electron - lỗ trống tuân theo thế tĩnh điện Coulomb thuần túy ($1/r$). Tuy nhiên, trong đơn lớp TMDs hai chiều, sự suy giảm số chiều không gian và sự chênh lệch hằng số điện môi giữa lớp vật liệu với môi trường xung quanh làm giảm hiệu ứng chắn điện môi ở khoảng cách ngắn. Thế Keldysh (kết hợp giữa hàm Struve $H_0$ và hàm Bessel loại hai $Y_0$) phản ánh chính xác tương tác này, giúp giải thích khớp các kết quả đo quang phổ thực nghiệm mà thế Coulomb không làm được.

2. Vai trò của tham số tự do $\omega$ trong phương pháp toán tử FK là gì và nó được xác định như thế nào?

Tham số $\omega$ được đưa vào khi định nghĩa các toán tử sinh – hủy Dirac đóng vai trò điều tiết độ rộng hàm sóng cơ sở, qua đó phân phối lại độ lớn giữa thành phần chính và thành phần nhiễu loạn của Hamiltonian. Do $\omega$ không xuất hiện trong Hamiltonian vật lý ban đầu, giá trị tối ưu của nó ($\omega_0$) được xác định bằng cách giải phương trình vi phân điều kiện cực trị năng lượng bậc không: $\frac{\partial E^{(0)}_{n_1 n_2}(\omega)}{\partial \omega} = 0$.

3. Mục đích của việc chuyển từ các toán tử sinh – hủy Dirac sang toán tử $\hat{a}, \hat{b}$ là gì?

Toán tử hình chiếu moment động lượng quỹ đạo lên trục $Oz$ ($\hat{L}_z$) giao hoán với Hamiltonian của bài toán exciton 2D trong từ trường vuông góc. Việc định nghĩa lại hệ toán tử $\hat{a}, \hat{b}$ giúp chéo hóa hoàn toàn toán tử moment động lượng về dạng trung hòa $\hat{L}_z = \hbar(\hat{a}^+\hat{a} - \hat{b}^+\hat{b})$, cho phép sử dụng trực tiếp các trạng thái lượng tử $|n_1, n_2\rangle$ làm bộ hàm sóng cơ sở chung cho cả $\hat{H}$ và $\hat{L}_z$.

4. Kết quả tính toán năng lượng liên kết trong khóa luận khớp với thực nghiệm ở mức độ nào?

Kết quả tính toán năng lượng liên kết trạng thái cơ bản $1s$ của exciton trong đơn lớp $\text{WSe}_2$ khi không có từ trường đạt $1.72311\text{ eV}$, sai lệch dưới $1%$ so với thực nghiệm của Crooker et al. (2018). Trên toàn dải từ trường ngoài từ 0 đến 60 Tesla, mức năng lượng bậc không của các trạng thái $1s, 2s, 3s$ bám sát đồ thị thực nghiệm với độ sai lệch chỉ từ $0%$ đến $2%$.

5. Nguyên nhân chính gây ra độ lệch 0 - 2% giữa kết quả lý thuyết trong khóa luận và dữ liệu thực nghiệm là gì?

Tác giả mới chỉ tính toán năng lượng ở mức gần đúng bậc không $E^{(0)}$ mà chưa cộng dồn các số hạng bổ chính bậc một $\Delta E^{(1)}$ và bậc hai $\Delta E^{(2)}$ theo lý thuyết nhiễu loạn không suy biến. Việc tính thêm các bổ chính bậc cao được xác định là hướng đi tiếp theo để đưa sai số về mức tiệm cận tuyệt đối.

Kết luận

Khóa luận tốt nghiệp của sinh viên Trần Vĩnh Phát đã ứng dụng thành công phương pháp đại số toán tử FK kết hợp biến đổi Fourier thế Keldysh để giải phương trình Schrödinger cho exciton 2D trong đơn lớp bán dẫn $\text{WSe}_2$ dưới tác dụng của từ trường đều $0 - 60\text{ T}$. Bằng các phép biến đổi thuần đại số và tối ưu hóa tham số $\omega$, nghiên cứu đã xác định chính xác phổ năng lượng liên kết của các trạng thái $1s, 2s, 3s$ với độ sai lệch so với dữ liệu thực nghiệm chỉ từ $0 - 2%$. Đây là một công trình nghiên cứu nghiêm túc, có giá trị học thuật cao trong việc phát triển các công cụ toán tử giải tích phục vụ tính toán cấu trúc điện tử và tính chất quang học của vật liệu hai chiều.