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

Sự thu nhỏ kích thước vi mạch silicon xuống dưới ngưỡng 10 nm đang đối mặt với các giới hạn vật lý nghiêm trọng, bao gồm hiện tượng rò rỉ electron lượng tử và sự gia tăng mật độ tỏa nhiệt. Trong bối cảnh tìm kiếm vật liệu thay thế cho công nghệ bán dẫn truyền thống, graphene—cấu trúc mạng tổ ong hai chiều của các nguyên tử carbon—nổi lên như một đối tượng nghiên cứu tiềm năng nhờ độ bền cơ học cao, độ linh động điện tử lớn và khả năng dẫn nhiệt vượt trội.

Về mặt lý thuyết, graphene lớp đơn (Monolayer Graphene – MLG) sở hữu hệ thức tán sắc năng lượng tuyến tính ở vùng năng lượng thấp gần các điểm Dirac, được mô tả bằng phương trình Dirac không khối lượng với vận tốc Fermi $v_F \approx 10^8\text{ cm/s}$. Ngược lại, graphene lớp kép (Bilayer Graphene – BLG) với kiểu xếp chồng AB lại thể hiện tán sắc dạng parabol gần mức năng lượng không, tương ứng với các chuẩn hạt có khối lượng hiệu dụng $m_g \approx 0{,}033m_0$. Khi kết hợp các lớp graphene thành cấu trúc dị thể nhiều lớp, sự tương tác Coulomb xuyên lớp giữa các chuẩn hạt có bản chất khác nhau tạo ra các đặc tính kích thích tập thể (plasmon) phức tạp. Hầu hết các nghiên cứu trước đây thường giả định môi trường điện môi nền đồng nhất với tensor thế Coulomb suy biến dạng hàm mũ đơn giản $e^{-qd}$. Tuy nhiên, trong các hệ thực nghiệm, các lớp graphene thường được đặt xen kẽ giữa các chất điện môi khác nhau, khiến tương tác Coulomb nội lớp và xuyên lớp phụ thuộc phi tuyến vào hằng số điện môi của từng lớp thành phần.

Khóa luận thực hiện các mục tiêu và nhiệm vụ cụ thể sau:

  • Mục tiêu nghiên cứu: Tính toán và khảo sát phổ kích thích tập thể (phổ plasmon) cùng hệ số suy giảm hấp thụ trong hệ ba lớp graphene hỗn hợp (gồm hai lớp đơn và một lớp kép: MLG–MLG–BLG) đặt trên nền điện môi không đồng nhất ở nhiệt độ không tuyệt đối ($T = 0\text{ K}$) trong gần đúng pha ngẫu nhiên (RPA).
  • Nhiệm vụ nghiên cứu:
    1. Nghiên cứu nội dung, kiến thức về đặc tính của hệ nhiều hạt có cấu trúc nhiều lớp.
    2. Nghiên cứu, tìm mối liên hệ giữa tần số góc và vectơ sóng trong hàm điện môi động (Dynamical Dielectric Function – DDF).
    3. Nghiên cứu các tài liệu hướng dẫn sử dụng phần mềm tính toán và xử lý dữ liệu Mathematica, Origin.
    4. Nghiên cứu sử dụng ngôn ngữ lập trình C++ để hỗ trợ giải phương trình điểm không của hàm điện môi động.
    5. Khai thác các tính năng giải phương trình và vẽ đồ thị trên Mathematica, Origin và Microsoft Visual Studio.
    6. Tính toán, vẽ đồ thị phổ plasmon và xác định vùng kích thích đơn hạt (Single Particle Excitation – SPE) cho hệ ba lớp graphene trên nền điện môi không đồng nhất.
  • Đối tượng nghiên cứu: Cấu trúc ba lớp graphene phẳng song song gồm hai lớp đơn (MLG1, MLG2) và một lớp kép (BLG3) phân bố tại các tọa độ $z = 0, d, 2d$ dọc trục $Oz$, cách đều nhau một khoảng $d$.
  • Phạm vi nghiên cứu: Khảo sát lý thuyết ở nhiệt độ không tuyệt đối ($T = 0\text{ K}$), trong gần đúng bước sóng dài ($q \to 0$); bỏ qua hiệu ứng xuyên ngầm (tunneling) giữa các lớp do giả định lớp cách điện đủ dày.

Cơ sở lý thuyết và phương pháp nghiên cứu

Khóa luận kết hợp mô hình liên kết chặt (tight-binding model), điện động lực học vi mô và gần đúng pha ngẫu nhiên (RPA) để mô tả phản hồi điện tử của hệ nhiều hạt.

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

  • Graphene lớp đơn (MLG): Mạng tinh thể lục giác với ô nguyên tố chứa hai nguyên tử carbon $A$ và $B$, hằng số mạng $a = 2{,}46\text{ \AA}$ (độ dài liên kết C–C là $1{,}42\text{ \AA}$). Hàm sóng Bloch xây dựng từ orbital $p_z$. Phương trình ma trận Hamiltonian $2 \times 2$ dẫn đến phổ tán sắc năng lượng $E_\pm(\mathbf{k}) \approx \pm \hbar v_F |\mathbf{k}|$ quanh hai điểm Dirac $K$ và $K'$, với độ suy biến spin $g_s = 2$ và suy biến thung lũng $g_v = 2$. Năng lượng Fermi tỉ lệ với căn bậc hai của mật độ hạt tải: $E_F = \hbar v_F \sqrt{4\pi n / (g_s g_v)}$.
  • Graphene lớp kép (BLG): Cấu trúc xếp chồng AB (Bernal stacking) với ô nguyên tố gồm bốn nguyên tử $A_1, B_1, A_2, B_2$. Mô hình Slonczewski–Weiss–McClure (SWM) sử dụng tham số tương tác nội lớp $\gamma_0 = 3{,}033\text{ eV}$, tương tác liên lớp dimer $\gamma_1 = 0{,}39\text{ eV}$ giữa $B_1$ và $A_2$, tích phân xen phủ $s_0 = 0{,}129$. Bằng biến đổi Schrieffer–Wolff, Hamiltonian bốn thành phần được rút gọn thành Hamiltonian hai thành phần năng lượng thấp: $$H_{2} = -\frac{1}{2m_g} \begin{pmatrix} 0 & (p_x - i p_y)^2 \ (p_x + i p_y)^2 & 0 \end{pmatrix}$$ với khối lượng hiệu dụng $m_g = \gamma_1 / (2 v_F^2) \approx 0{,}033 m_0$, cho phổ tán sắc năng lượng dạng parabol $E(p) = \pm p^2 / (2m_g)$.
  • Hàm phân cực động $\Pi(q, \omega)$: Trong gần đúng RPA ở $T = 0\text{ K}$, hàm phân cực của MLG được trích xuất theo mô hình Hwang–Das Sarma. Hàm phân cực của BLG được tách thành thành phần chuyển mức nội vùng (intraband) và liên vùng (interband) theo tính toán của SenSarma và cộng sự.
  • Hàm điện môi động và thế tương tác Coulomb: Bỏ qua hiệu ứng xuyên ngầm, tensor hàm phản hồi có dạng đường chéo $\hat{\Pi}{ij}(q, \omega) = \delta{ij} \Pi_i^{(0)}(q, \omega)$. Hàm điện môi động của hệ được xác định qua định thức: $$\epsilon(q, \omega) = \det\left[ \hat{I} - \hat{v}(q)\hat{\Pi}(q, \omega) \right]$$ Tensor thế Coulomb $\hat{v}(q)$ được xác định chính xác từ nghiệm của phương trình Poisson với các điều kiện biên liên tục điện thế và bước nhảy dịch chuyển điện trường tại các mặt phân cách $z = 0, d, 2d$. Biểu thức phần tử thế $v_{ij}(q) = \frac{2\pi e^2}{q} f_{ij}(q)$ phụ thuộc vào các hàm hyperbolic $\sinh(qd), \cosh(qd)$ và các hằng số điện môi nền $\epsilon_1, \epsilon_2, \epsilon_3, \epsilon_4$.

Phương pháp nghiên cứu và công cụ thực hiện

  • Phương pháp giải tích: Thiết lập biểu thức giải tích của tensor tương tác Coulomb và hàm điện môi động $\epsilon(q, \omega)$; khai triển tiệm cận trong gần đúng bước sóng dài ($q \to 0$).
  • Phương pháp giải số: Xây dựng thuật toán trên ngôn ngữ C++ (biên dịch bằng Microsoft Visual Studio) để quét lưới không gian $(q, \omega)$, tìm nghiệm thực của phương trình triệt tiêu phần thực hàm điện môi $\text{Re},\epsilon(q, \omega_p) = 0$ khi hấp thụ yếu, sau đó tính tốc độ phân rã plasmon $\gamma(q) = -\text{Im},\epsilon(q, \omega_p) / [\partial \text{Re},\epsilon(q, \omega_p)/\partial \omega]$.
  • Phần mềm bổ trợ: Wolfram Mathematica (hỗ trợ kiểm chứng biến đổi giải tích) và Origin (vẽ đồ thị phổ tán sắc plasmon, hệ số suy giảm và ranh giới vùng SPE).

Tham số mô phỏng

Hệ số điện môi của các lớp nền không đồng nhất gồm: $\text{SiO}_2$ ($\epsilon_1 = 3{,}8$), boron nitride $\text{BN}$ ($\epsilon_2 = 5{,}0$), hexagonal boron nitride $\text{hBN}$ ($\epsilon_3 = 3{,}0$) và không khí ($\epsilon_4 = 1{,}0$). Các thông số hình học và mật độ hạt tải chuẩn: khoảng cách lớp $d = 20\text{ nm}$ và $50\text{ nm}$; mật độ hạt tải $n_1, n_2, n_3$ biến thiên quanh giá trị $10^{12}\text{ cm}^{-2}$.

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

Chương 1: Lý do chọn đề tài

Chương này phân tích tổng quan bối cảnh khoa học về giới hạn của vật liệu bán dẫn silicon và các ưu điểm nổi bật của graphene. Tác giả tổng thuật lịch sử nghiên cứu phổ plasmon từ các hệ khí điện tử hai chiều (2DEG, Q2DEG) đến các cấu trúc đơn lớp (MLG), lớp kép (BLG), hệ dị thể MLG–2DEG, BLG–2DEG và các cấu trúc nhiều lớp graphene. Chương 1 nêu rõ lỗ hổng nghiên cứu: chưa có công bố nào khảo sát hệ ba lớp hỗn hợp MLG–MLG–BLG có xét đến đầy đủ tính không đồng nhất của điện môi nền ở nhiệt độ không tuyệt đối ($T = 0\text{ K}$). Trên cơ sở đó, tác giả xác định đối tượng, mục tiêu, 6 nhiệm vụ nghiên cứu cụ thể, giả thuyết khoa học, phương pháp tiếp cận và đóng góp mới của đề tài.

Chương 2: Tổng quan về graphene và phương pháp xác định phổ plasmon

Chương 2 trình bày chi tiết cơ sở lý thuyết cấu trúc điện tử của hai dạng vật liệu thành phần và phương pháp hình thức luận xác định plasmon:

  • Graphene lớp đơn (MLG): Trình bày cấu trúc mạng không gian thực, mạng đảo, vùng Brillouin thứ nhất với các điểm đối xứng cao $\Gamma, K, M$. Dẫn xuất ma trận Hamiltonian từ phương pháp liên kết chặt, giải phương trình thế kỷ tìm năng lượng $E_\pm(\mathbf{k})$, khai triển gần điểm Dirac dẫn đến phương trình Dirac không khối lượng với giả spin (pseudospin/isospin) và chỉ số chirality. Định nghĩa các đại lượng: vectơ sóng Fermi $k_F = \sqrt{4\pi n / (g_s g_v)}$, năng lượng Fermi $E_F = \hbar v_F k_F$, mật độ trạng thái $g(E) = g_s g_v |E| / [2\pi (\hbar v_F)^2]$, và tham số tương tác Wigner–Seitz $r_s = e^2 / (\epsilon \hbar v_F) \in [0; 2{,}2]$ (đặc trưng cho hệ tương tác yếu).
  • Graphene lớp kép (BLG): So sánh kiểu xếp chồng AA và AB. Trong kiểu xếp chồng AB, ma trận Hamiltonian liên kết chặt $4 \times 4$ trong gần đúng tối thiểu chứa các tham số $\gamma_0, \gamma_1$. Sử dụng phép biến đổi Schrieffer–Wolff để chiếu không gian Hilbert về hai trạng thái năng lượng thấp không dimer $A_1, B_2$, thiết lập Hamiltonian tựa Dirac chứa toán tử động lượng bậc hai với khối lượng hiệu dụng $m_g = 0{,}033m_0$.
  • Phương pháp xác định phổ plasmon: Từ phương trình sóng điện từ Maxwell, dao động plasma phân cực dọc tương ứng với điều kiện triệt tiêu của hàm điện môi động $\epsilon(q, \omega_p - i\gamma) = 0$. Khi suy giảm yếu ($\gamma \ll \omega_p$), phổ plasmon là nghiệm của $\text{Re},\epsilon(q, \omega_p) = 0$, và tốc độ suy giảm được tính qua đạo hàm riêng của phần thực và giá trị phần ảo của $\epsilon(q, \omega)$.

Chương 3: Phổ kích thích tập thể trong hệ ba lớp MLG–MLG–BLG

Chương 3 là phần nghiên cứu trọng tâm của khóa luận, bao gồm các tính toán giải tích và kết quả mô phỏng số:

  • Dẫn xuất hàm phân cực: Trình bày dạng giải tích phần thực và phần ảo của $\Pi_{\text{MLG}}(q, \omega)$ ở $T = 0\text{ K}$ theo biến không thứ nguyên $x = q/k_F, y = \hbar\omega/E_F$; tính toán tích phân góc cho các thành phần nội vùng $\Pi_{\text{intra}}$ và liên vùng $\Pi_{\text{inter}}$ của $\Pi_{\text{BLG}}(q, \omega)$. Trong giới hạn bước sóng dài, hàm phân cực của MLG và BLG đều tỉ lệ thuận với $q^2/\omega^2$.
  • Thiết lập tensor thế Coulomb $v_{ij}(q)$: Giải phương trình Poisson cho cấu trúc 4 môi trường điện môi ($\epsilon_1, \epsilon_2, \epsilon_3, \epsilon_4$) phân cách bởi 3 mặt phẳng graphene tại $z=0, d, 2d$. Xác định tường minh các hệ số $f_{11}, f_{22}, f_{33}, f_{12}, f_{23}, f_{13}$ thông qua định thức ma trận truyền điều kiện biên $M(qd)$.
  • Phân tích kết quả phổ plasmon và hấp thụ:
    • Đặc trưng nhánh phổ: Cấu trúc ba lớp xuất hiện đúng ba nhánh plasmon riêng biệt: 1 nhánh quang học (Optical – Op) có tần số cao nhất tỉ lệ với $\sqrt{q}$ khi $q \to 0$, và 2 nhánh âm học (Acoustic – $\text{Ac}_1, \text{Ac}_2$) có tần số thấp hơn tỉ lệ tuyến tính với $q$.
    • Đặc tính suy giảm: Nhánh Op và $\text{Ac}_1$ đi vào vùng kích thích đơn hạt (SPE) khá sớm (tại $q \approx 0{,}15k_F$ và $q \approx 0{,}3k_F$) nhưng tốc độ mất mát năng lượng diễn ra chậm. Nhánh $\text{Ac}_2$ đi vào vùng SPE muộn hơn nhưng bị hấp thụ rất mạnh và biến mất hoàn toàn tại $q \approx 0{,}6k_F$.
    • Ảnh hưởng của khoảng cách lớp $d$: Khi tăng $d$ từ $20\text{ nm}$ lên $50\text{ nm}$, nhánh Op giảm tần số rõ rệt, trong khi hai nhánh $\text{Ac}_1, \text{Ac}_2$ lại tăng tần số.
    • Ảnh hưởng của mật độ hạt tải bất đối xứng: Khi giảm mật độ hạt tải ở lớp MLG giữa ($n_1 = 2n_2 = n_3 = 10^{12}\text{ cm}^{-2}$), phổ plasmon biến đổi không đáng kể. Tuy nhiên, khi giảm mật độ hạt tải ở lớp kép BLG ($n_1 = n_2 = 2n_3 = 10^{12}\text{ cm}^{-2}$), hai nhánh có tần số cao sít lại gần nhau, và nhánh $\text{Ac}_2$ bị rút ngắn chiều dài đáng kể, kết thúc tại $q \approx 0{,}25k_F$ (so với $0{,}55k_F$ khi mật độ cân bằng).
    • Ảnh hưởng của tính không đồng nhất điện môi nền: So sánh hệ thực tế ($\epsilon_1=3{,}8; \epsilon_2=5{,}0; \epsilon_3=3{,}0; \epsilon_4=1{,}0$) với hệ đồng nhất giả định ($\bar{\epsilon} = (\epsilon_1+\epsilon_4)/2 = 2{,}4$) cho thấy: điện môi không đồng nhất làm suy giảm mạnh tần số của cả 3 nhánh plasmon, trong đó mức độ dịch chuyển tần số xảy ra mạnh nhất ở nhánh Op và nhánh $\text{Ac}_1$.

Kết quả và đóng góp của nghiên cứu

Khóa luận đã giải quyết bài toán kích thích tập thể cho hệ dị thể ba lớp graphene hỗn hợp trên nền điện môi thực tế.

Đại lượng / Hiện tượng Đặc tính trong hệ ba lớp MMBLG So sánh / Ý nghĩa vật lý
Số lượng nhánh plasmon 3 nhánh phân biệt: 1 nhánh quang học (Op) và 2 nhánh âm học ($\text{Ac}_1, \text{Ac}_2$). Phù hợp với bậc của định thức hàm điện môi $3 \times 3$ trong gần đúng RPA.
Quy luật tán sắc bước sóng dài Nhánh Op: $\omega \sim \sqrt{q}$; Nhánh Ac: $\omega \sim q$. Nhánh Op đặc trưng cho dao động đồng pha của hạt tải; các nhánh Ac đặc trưng cho dao động ngược pha.
Đặc tính hấp thụ trong vùng SPE Nhánh Op và $\text{Ac}_1$ suy giảm chậm sau khi vào SPE; nhánh $\text{Ac}_2$ suy giảm nhanh và tắt hẳn ở $q \approx 0{,}6k_F$. Sự hấp thụ năng lượng phụ thuộc vào vị trí tương đối giữa nhánh phổ và biên SPE nội vùng/liên vùng.
Tác động của khoảng cách lớp ($d$) $d$ tăng từ $20\text{ nm} \to 50\text{ nm}$: tần số Op giảm, tần số $\text{Ac}_1, \text{Ac}_2$ tăng. Tương tác Coulomb xuyên lớp giảm khi tăng khoảng cách, làm yếu liên kết pha giữa các lớp.
Tác động của mật độ hạt tải BLG ($n_3$) $n_3$ giảm xuống $0{,}5 \times 10^{12}\text{ cm}^{-2}$: hai nhánh cao sít lại, nhánh $\text{Ac}_2$ kết thúc sớm ở $q \approx 0{,}25k_F$. Chuẩn hạt có khối lượng trong BLG quyết định độ dài tồn tại của nhánh plasmon âm học thấp nhất.
Ảnh hưởng của điện môi không đồng nhất Tần số của cả 3 nhánh phổ đều giảm so với nền đồng nhất $\bar{\epsilon}$. Hiệu ứng chắn Coulomb thay đổi do sự phân cực bất đối xứng của các lớp điện môi $\text{SiO}_2/\text{BN}/\text{hBN}/\text{kh^ong kh'i}$.

Đóng góp mới của tác giả

  • Thiết lập hoàn chỉnh mô hình lý thuyết tính toán phổ plasmon và tốc độ phân rã cho hệ dị thể ba lớp hỗn hợp MLG–MLG–BLG có xét đến cấu trúc điện môi nền nhiều lớp thực tế.
  • Chỉ ra sự bất đối xứng trong phản hồi plasmon khi thay đổi mật độ hạt tải: lớp BLG chi phối mạnh cấu trúc phân nhánh và vùng tồn tại không suy giảm của plasmon âm học, trong khi lớp MLG thứ hai ít ảnh hưởng hơn.
  • Cung cấp dữ liệu định lượng phục vụ công bố khoa học (bài báo "Kích thích tập thể trong hệ 3 lớp graphene hỗn hợp ở nhiệt độ không tuyệt đối" gửi đăng tại Tạp chí Khoa học Trường Đại học Tân Trào, đề tài C2021-16-06 do ĐHQG-HCM tài trợ).

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

  • Hạn chế của khóa luận:
    • Mô hình tính toán chỉ thực hiện ở nhiệt độ không tuyệt đối ($T = 0\text{ K}$), chưa khảo sát ảnh hưởng của nhiệt độ hữu hạn đến sự mở rộng nhiệt của phân bố Fermi–Dirac.
    • Bỏ qua tán xạ hạt tải với dao động mạng tinh thể (phonon) và tán xạ do tạp chất pha tạp.
    • Chưa xét đến sự xuất hiện của độ rộng vùng cấm năng lượng (bandgap opening) trong graphene dưới tác dụng của điện trường ngoài hoặc chất nền bất đối xứng.
    • Mô hình giới hạn ở cấu trúc ba lớp, chưa tổng quát hóa cho hệ đa lớp tùy ý ($N > 3$).
  • Hướng phát triển tiếp theo:
    • Mở rộng bài toán sang thang nhiệt độ hữu hạn ($T > 0$) để đánh giá sự dịch chuyển tần số plasmon và độ mở rộng vạch phổ do hiệu ứng nhiệt.
    • Bổ sung tương tác electron–phonon và đánh giá cơ chế suy giảm plasmonic Landau.
    • Khảo sát cấu trúc graphene lớp kép có vùng cấm điều khiển được (gapped BLG) và nghiên cứu vận tốc sóng âm trong hệ.
    • Tính toán bổ sung các hiệu ứng phân cực spin, hiệu ứng tương quan–trao đổi cục bộ và ảnh hưởng tái cấu trúc vùng năng lượng do tương tác với chất nền.

Giá trị tham khảo

Tài liệu này hữu ích cho các đối tượng nghiên cứu và học tập sau:

  • Sinh viên, học viên cao học và nghiên cứu sinh: Thuộc các chuyên ngành Vật lý lý thuyết và Vật lý toán, Vật lý chất rắn, Vật lý kỹ thuật, Khoa học và Kỹ thuật vật liệu đang thực hiện các đề tài về vật liệu hai chiều (2D materials), dị thể van der Waals và quang tử plasmon (plasmonics).
  • Nội dung tham khảo trọng tâm:
    • Phương pháp giải tích thiết lập tensor thế tương tác Coulomb $v_{ij}(q)$ cho hệ phân lớp trên nền điện môi bất đối xứng qua phương trình vi phân Poisson.
    • Phương pháp thu gọn Hamiltonian bằng biến đổi Schrieffer–Wolff cho cấu trúc năng lượng của BLG.
    • Quy trình lập trình số trên C++ kết hợp Mathematica để giải hệ phương trình phi tuyến xác định điểm không của hàm điện môi động và vẽ giản đồ tán sắc.

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

1. Cấu trúc hệ ba lớp trong khóa luận được bố trí như thế nào?

Hệ gồm ba lớp graphene đặt song song dọc theo trục $Oz$ tại các vị trí $z = 0$ (lớp đơn MLG1), $z = d$ (lớp đơn MLG2) và $z = 2d$ (lớp kép BLG3), cách đều nhau khoảng cách $d$. Hệ được nhúng trong môi trường điện môi gồm 4 lớp phân tách: lớp dưới cùng là $\text{SiO}_2$ ($\epsilon_1 = 3{,}8$), lớp giữa MLG1 và MLG2 là $\text{BN}$ ($\epsilon_2 = 5{,}0$), lớp giữa MLG2 và BLG3 là $\text{hBN}$ ($\epsilon_3 = 3{,}0$) và môi trường trên cùng là không khí ($\epsilon_4 = 1{,}0$).

2. Sự khác biệt cơ bản về cấu trúc dải năng lượng giữa MLG và BLG trong mô hình là gì?

Ở vùng năng lượng thấp gần các điểm Dirac, MLG thể hiện hệ thức tán sắc tuyến tính $E(k) = \pm \hbar v_F k$, mô tả hạt Dirac không khối lượng với vận tốc $v_F \approx 10^8\text{ cm/s}$. Trong khi đó, BLG (kiểu xếp chồng AB) thể hiện hệ thức tán sắc dạng parabol $E(p) = \pm p^2/(2m_g)$, mô tả các chuẩn hạt có khối lượng hiệu dụng $m_g \approx 0{,}033m_0$, được xác định từ tham số tương tác liên lớp dimer $\gamma_1 = 0{,}39\text{ eV}$ thông qua biến đổi Schrieffer–Wolff.

3. Việc tăng khoảng cách giữa các lớp từ 20 nm lên 50 nm làm thay đổi các nhánh phổ plasmon ra sao?

Khi tăng khoảng cách $d$ từ $20\text{ nm}$ lên $50\text{ nm}$ (với mật độ hạt tải cố định $n_1 = n_2 = n_3 = 10^{12}\text{ cm}^{-2}$), tần số của nhánh quang học (Op) giảm rõ rệt, trong khi tần số của hai nhánh âm học ($\text{Ac}_1, \text{Ac}_2$) lại tăng lên. Điều này chứng minh tương tác Coulomb xuyên lớp giảm khi tăng khoảng cách, dẫn đến sự dịch chuyển ngược chiều giữa dao động đồng pha và ngược pha của các lớp hạt tải.

4. Mật độ hạt tải của lớp nào có ảnh hưởng quyết định đến hình dạng phổ plasmon?

Sự suy giảm mật độ hạt tải ở lớp kép BLG ($n_3$) gây ảnh hưởng mạnh nhất: khi giảm $n_3$ xuống còn $0{,}5 \times 10^{12}\text{ cm}^{-2}$, hai nhánh tần số cao sít lại gần nhau và nhánh âm học thấp nhất ($\text{Ac}_2$) bị suy giảm chiều dài nghiêm trọng, kết thúc tại $q \approx 0{,}25k_F$ thay vì $0{,}55k_F$. Ngược lại, việc giảm mật độ hạt tải ở lớp đơn ở giữa ($n_2$) hầu như không làm biến đổi đáng kể hình dạng các đường phổ.

5. Vì sao điện môi nền không đồng nhất làm giảm tần số phổ plasmon so với nền đồng nhất?

Sự xuất hiện của các lớp điện môi có hằng số điện môi lớn xen kẽ ($\epsilon_2 = 5{,}0, \epsilon_1 = 3{,}8, \epsilon_3 = 3{,}0$) làm tăng hiệu ứng chắn tĩnh điện (dielectric screening), làm suy giảm cường độ thế tương tác Coulomb nội lớp và xuyên lớp giữa các hạt tải so với mô hình nền đồng nhất lấy trung bình $\bar{\epsilon} = 2{,}4$. Do đó, năng lượng kích thích tập thể của toàn hệ thống bị kéo xuống mức tần số thấp hơn.

Kết luận

Khóa luận đã hoàn thành việc xây dựng mô hình lý thuyết và tính toán số phổ kích thích tập thể cùng hệ số hấp thụ cho hệ ba lớp graphene hỗn hợp MLG–MLG–BLG trên nền điện môi không đồng nhất ở nhiệt độ $T = 0\text{ K}$ trong gần đúng RPA. Kết quả nghiên cứu làm sáng tỏ sự tồn tại của ba nhánh plasmon (một nhánh quang học và hai nhánh âm học) với các quy luật suy giảm năng lượng và phụ thuộc khoảng cách lớp khác biệt. Nghiên cứu đã chứng minh vai trò then chốt của lớp graphene kép và tính không đồng nhất của chất nền điện môi trong việc điều khiển tần số cũng như vùng tồn tại của dao động plasma, đóng góp cơ sở dữ liệu quan trọng cho các nghiên cứu lý thuyết và ứng dụng plasmonic của vật liệu hai chiều.