Tổng quan luận án

Vật liệu graphene—mạng đơn lớp nguyên tử carbon hai chiều—được phát hiện thực nghiệm vào năm 2004 và đem lại giải Nobel Vật lý năm 2010 cho Andre Geim và Konstantin Novoselov. Graphene sở hữu các đặc tính điện tử nổi bật như độ linh động hạt tải đạt tới $2 \times 10^5\text{ cm}^2/\text{Vs}$ ở nhiệt độ phòng (gấp 100 lần silicon và 20 lần GaAs) cùng khả năng tương tác mạnh với ánh sáng khi độ hấp thụ quang của một lớp nguyên tử đơn lẻ đạt xấp xỉ 2,3%. Tuy nhiên, do không có vùng cấm năng lượng ($E_g = 0$), việc ứng dụng graphene làm kênh dẫn trong các linh kiện bán dẫn trường (MOSFET) gặp trở ngại về khả năng ngắt dòng. Ngược lại, các đặc tính động lực của hệ electron trong graphene tỏ ra phù hợp cho các nhánh quang điện tử (optoelectronics), quang tử (photonics) và đặc biệt là nano-plasmonics. Khác với kim loại truyền thống có tần số plasmon nằm ở vùng tử ngoại hoặc ánh sáng nhìn thấy, tần số plasmon của graphene nằm trong dải terahertz (THz) và hồng ngoại, đồng thời có thể điều khiển liên tục thông qua điện thế phân cực tĩnh điện để thay đổi nồng độ hạt tải.

Khoảng trống nghiên cứu then chốt được tác giả luận án chỉ ra là: phần lớn các nghiên cứu lý thuyết trước đây chủ yếu tiếp cận ở hai thái cực năng lượng. Tại giới hạn năng lượng thấp và bước sóng dài, các nghiên cứu dựa trên mô hình Dirac lý tưởng xem các mặt năng lượng xung quanh điểm K là hình nón tròn xoay đẳng hướng, dẫn đến quy luật tán sắc Dirac plasmon tỷ lệ với căn bậc hai của vectơ sóng ($q^{1/2}$). Tại giới hạn năng lượng cao, các nghiên cứu thực nghiệm ghi nhận $\pi$-plasmon liên quan đến các chuyển dời sâu trong vùng hóa trị. Mô hình Dirac không thể mô tả được $\pi$-plasmon cũng như bỏ qua sự biến dạng hình học thực tế (tính bất đẳng hướng) của các mặt đẳng năng khi năng lượng kích thích và mức độ pha tạp tăng cao. Do đó, cần có một nghiên cứu khảo sát chi tiết giới hạn áp dụng của mô hình Dirac, tính đến cấu trúc vùng năng lượng thực tế và hiệu ứng trường địa phương để làm rõ các trạng thái kích thích tập thể.

  • Mục tiêu nghiên cứu: Khảo sát các tính chất động lực học của hệ electron dẫn trong màng graphene dưới các điều kiện trường ngoài khác nhau; xác định cơ chế hình thành, tính chất phân cực và điều kiện duy trì các trạng thái kích thích tập thể (plasmon) của hệ electron trong các chế độ pha tạp từ thấp đến cao, làm rõ tiềm năng ứng dụng trong lĩnh vực nano-plasmonics.
  • Câu hỏi/Nhiệm vụ nghiên cứu:
    1. Giới hạn hoạt động của mô hình Dirac tuyến tính trong việc mô tả các hàm phản ứng điện môi và phổ tán sắc plasmon năng lượng thấp là gì khi tính đến nhiệt độ và nồng độ pha tạp hữu hạn?
    2. Hình học mặt năng lượng bất đẳng hướng và tương tác với lân cận thứ hai ảnh hưởng như thế nào đến hàm phân cực và tính phân cực của plasmon ở chế độ pha tạp mạnh?
    3. Liệu có tồn tại các mode kích thích plasmon mới khi graphene ở trạng thái pha tạp cao và khi xét đến các chuyển dời giữa hai thung lũng Dirac không tương đương?
    4. Hiệu ứng trường địa phương (LFE) tác động như thế nào đến hàm điện môi vĩ mô và phổ plasmon ở giới hạn bước sóng ngắn ứng với sự chuyển trạng thái giữa các điểm K (inter-valley)?
  • Đối tượng nghiên cứu: Hệ electron hai chiều (2DEG) của các orbital liên kết $p_z$ ($\pi$ và $\pi^*$) trong mạng tinh thể graphene thuần khiết và graphene trong cấu trúc siêu mạng (GSL) ở các chế độ pha tạp khác nhau.
  • Phạm vi nghiên cứu:
    • Không gian: Hệ vật liệu hai chiều trong mặt phẳng mạng tinh thể graphene dạng tổ ong.
    • Điều kiện vật lý và tham số: Nhiệt độ từ $0\text{ K}$ đến nhiệt độ phòng; mức năng lượng Fermi $E_F$ từ mức trung hòa điện tích ($E_F = 0$) đến chế độ pha tạp cao ($E_F = 0,5\text{ eV}$ và $E_F = 1,0\text{ eV}$); dải năng lượng photon và tần số kích thích từ vùng THz, hồng ngoại đến vùng quang học năng lượng cao; miền vectơ sóng trong toàn bộ vùng Brillouin thứ nhất (BZ).

Tổng quan tài liệu và vị trí của luận án

Luận án đã hệ thống hóa các hướng nghiên cứu lý thuyết và thực nghiệm liên quan đến động lực học electron và plasmonics:

  • Nền tảng lý thuyết plasmon và thực nghiệm EELS: Khái niệm dao động plasma và sự lượng tử hóa thành plasmon bắt nguồn từ các công trình của David Pines (1956) [121–123], Heinz Raether (1956–1968) về phổ tổn hao năng lượng electron (EELS), plasmon bề mặt (SPs) và sóng phân cực plasmon bề mặt (SPPs). Các kỹ thuật kích thích quang học như phản xạ toàn phần suy giảm (ATR) của Kretschmann và Otto [112], hệ cách tử phản xạ [23, 131] và hiện tượng tăng cường tín hiệu Raman bề mặt (SERS) [16, 23] đã đặt nền móng cho plasmonics kim loại.
  • Lý thuyết phản ứng điện môi và gần đúng nhiều hạt: Ehrenreich và Cohen (1959) [40] xây dựng phương pháp trường tự hợp (SCF) và gần đúng pha ngẫu nhiên (RPA) để tính hàm điện môi phụ thuộc thời gian và không gian. Lý thuyết Lindhard mở rộng kết hợp hiệu ứng trường địa phương (LFE) cho hệ không đồng nhất được phát triển bởi Adler (1962) và Wiser (1963) [4, 158, 159].
  • Cấu trúc điện tử và độ dẫn quang của graphene: Mô hình liên kết chặt (Tight-Binding - TB) cho graphene được Reich và cộng sự [133] phát triển giải tích đến lân cận thứ ba kết hợp dữ liệu tính toán từ nguyên lý đầu (DFT). Bena và cộng sự [20] phân tích sự phụ thuộc của hàm sóng TB vào hệ cơ sở nguyên tử. Nghiên cứu về độ dẫn quang phổ quát $\sigma_0 = e^2/(4\hbar)$ và độ dẫn nội dải/ngoại dải được xác lập qua các công trình của Gusynin, Sharapov, Carbotte [57], Falkovsky [42] và Nair cùng cộng sự [104]. Barbier và cộng sự [12, 13] khởi xướng mô hình siêu mạng graphene (GSL) dưới thế tĩnh điện tuần hoàn.
  • Các nghiên cứu về Dirac plasmon và $\pi$-plasmon: Das Sarma, Hwang [65, 71], Wunsch [160] và Polini [126] đã thiết lập phổ tán sắc căn bậc hai $\omega_p \propto \sqrt{q}$ cho Dirac plasmon trong gần đúng bước sóng dài. Nhóm nghiên cứu của Eberlein [40] và Kramberger [75, 90] khảo sát $\pi$-plasmon ở năng lượng cao bằng thực nghiệm EELS.

Vị trí và khoảng trống luận án giải quyết: Luận án kết nối miền chuyển tiếp giữa mô hình Dirac năng lượng thấp và mô hình toàn vùng BZ. Bằng việc áp dụng mô hình liên kết chặt có xét đến lân cận thứ hai (NNN) và phương pháp tính số hàm điện môi RPA toàn phần, luận án đã làm rõ sự phá vỡ tính đẳng hướng của phổ tán sắc plasmon khi mức Fermi tăng cao, giải thích sự xuất hiện của mode plasmon mới ở chế độ pha tạp mạnh, và bước đầu tính toán định lượng kích thích plasmon chuyển trạng thái inter-valley có xét đến hiệu ứng trường địa phương.

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

Luận án kết hợp phương pháp giải tích vi mô, lý thuyết điện từ học vĩ mô và tính toán mô phỏng số:

Khung lý thuyết và mô hình giải tích

  1. Lý thuyết điện từ vĩ mô: Hệ phương trình Maxwell cho môi trường không từ tính ($\mu = 1$), thiết lập liên hệ giữa vectơ điện dịch $\vec{D}$, điện trường $\vec{E}$, độ phân cực $\vec{P}$ và mật độ dòng cảm ứng $\vec{J}_{\text{ind}}$. Mối liên hệ một-một giữa hàm điện môi $\epsilon(q,\omega)$ và độ dẫn quang $\sigma(q,\omega)$ qua biến đổi Fourier: $$\epsilon(q,\omega) = 1 + \frac{i\sigma(q,\omega)}{\epsilon_0 \omega}$$ Điều kiện tồn tại dao động dọc (plasmon) tương ứng với các không điểm của hàm điện môi: $\text{Re}[\epsilon(q,\omega_p + i\gamma)] = 0$, trong đó $\omega_p$ là tần số dao động tập thể và $\gamma$ là tốc độ phân rã plasmon.
  2. Lý thuyết phản ứng tuyến tính và Gần đúng pha ngẫu nhiên (RPA): Áp dụng phương pháp trường tự hợp (SCF) của Ehrenreich-Cohen. Hàm điện môi được xác định qua hàm phân cực tái chuẩn hóa: $$\epsilon_{\text{RPA}}(q,\omega) = 1 - V_{2D}(q) P(q,\omega)$$ với $V_{2D}(q) = \frac{e^2}{2\epsilon_0 \epsilon_r q}$ là thế Coulomb trong không gian hai chiều.
  3. Mô hình liên kết chặt (Tight-Binding - TB):
    • Gần đúng lân cận gần nhất (NN): Xét tương tác giữa các orbital $p_z$ của carbon với năng lượng nhảy $t \approx 2,7\text{ - }3,1\text{ eV}$ ($t_{\text{cc}} = 2,67\text{ eV}$), tích phân chồng chập $s \approx 0,129$. Hàm tán sắc đối xứng electron - lỗ trống: $$E_k^l = \frac{l t |f(k)|}{1 - l s |f(k)|} \quad (l = \pm 1)$$
    • Gần đúng lân cận thứ hai (NNN): Bổ sung biên độ nhảy $t' \approx 0,1\text{ - }0,2t$ và tích phân chồng chập $s'$, phá vỡ đối xứng giữa dải dẫn và dải hóa trị: $$E_k^l = l t |f(k)| + t' |f(k)|^2$$
    • Mô hình Dirac năng lượng thấp: Tại lân cận các điểm $K, K'$, Hamiltonian hiệu dụng có dạng $H = \hbar v_F \vec{\sigma} \cdot \vec{q}$, với vận tốc Fermi $v_F = \frac{3 t a_{\text{cc}}}{2\hbar} \approx 10^6\text{ m/s}$ (khoảng cách liên kết carbon-carbon $a_{\text{cc}} \approx 0,142\text{ nm}$).
  4. Lý thuyết độ dẫn Kubo cho siêu mạng (GSL): Tính toán tensơ độ dẫn quang $\sigma_{\alpha\beta}(\omega)$ qua các phần tử ma trận toán tử vận tốc và hàm phân bố Fermi-Dirac $f(E)$.

Phương pháp tính toán và nguồn dữ liệu

  • Phương pháp giải tích: Tích phân giải tích hàm phân cực $P(q,\omega)$ cho graphene pha tạp yếu ở $T = 0\text{ K}$, phân tách phần thực $\text{Re}[P(q,\omega)]$ và phần ảo $\text{Im}[P(q,\omega)]$ trên các miền kích thích đơn hạt (SPE) nội dải (intra-band) và ngoại dải (inter-band).
  • Phương pháp số: Xây dựng thuật toán rời rạc hóa toàn bộ vùng Brillouin thứ nhất (BZ) để lấy tích phân số hàm phân cực; tìm nghiệm phương trình $\text{Re}[\epsilon_{\text{RPA}}(q,\omega)] = 0$; tính phổ tổn hao năng lượng electron: $$S(q,\omega) = -\text{Im}\left[\frac{1}{\epsilon(q,\omega)}\right]$$
  • Hệ đơn vị: Sử dụng hệ đo lường quốc tế SI trong toàn bộ các công thức tính toán.
Tham số / Đại lượng Ký hiệu Giá trị / Biểu thức trong luận án
Khoảng cách liên kết C-C $a_{\text{cc}}$ $0,142\text{ nm}$ ($1,42\text{ \AA}$)
Năng lượng nhảy lân cận gần nhất $t$ ($t_{\text{cc}}$) $2,67\text{ eV} \div 3,0\text{ eV}$
Năng lượng nhảy lân cận thứ hai $t'$ $0,1t \div 0,2t$ ($\sim 0,1\text{ eV}$)
Tích phân chồng chập orbital $s$ $\sim 0,129$ ($s < 0,2$)
Vận tốc Fermi $v_F$ $3 t a_{\text{cc}} / (2\hbar) \approx 10^6\text{ m/s}$
Độ dẫn quang phổ quát $\sigma_0$ $e^2 / (4\hbar)$
Hằng số điện môi chân không $\epsilon_0$ $8,854 \times 10^{-12}\text{ F/m}$
Chu kỳ siêu mạng A-GSL $L$ $3 N a_{\text{cc}}$ ($N = 15 \div 30$, $L \approx 12\text{ nm}$)
Thế năng điều chế siêu mạng $U_0$ $2\pi\hbar v_F / L$

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

Chương 1: Cơ sở lý thuyết nghiên cứu các tính chất động lực của hệ điện tử và các tính chất vật lý cơ bản của hệ điện tử hai chiều trong mạng graphene

Chương 1 thiết lập nền tảng lý thuyết từ phương trình Maxwell vi mô đến vĩ mô, phân tích phản ứng của môi trường vật chất với trường điện từ sóng ngang ($q \cdot E = 0$, tán sắc $q^2 = \frac{\omega^2}{c^2}\epsilon$) và sóng dọc ($q \parallel E$, điều kiện lan truyền $\epsilon(q,\omega) = 0$). Tác giả trình bày mô hình dao động tử Lorentz, độ phân cực nguyên tử $\alpha_0(\omega) = \frac{e^2/m}{\omega_0^2 - \omega^2 - i\gamma\omega}$, và làm rõ bản chất của hiệu ứng trường địa phương (LFE) khi điện trường cục bộ $E_{\text{loc}}$ khác biệt so với điện trường trung bình vĩ mô $E$.

         LỤC GIÁC (Ô THÔI 2 NGUYÊN TỬ)         CHỮ NHẬT (Ô 4 NGUYÊN TỬ A-GSL)
               Ky                                     Ky

Trong phần cấu trúc mạng graphene, tác giả phân tích sự lai hóa $sp^2$ tạo 3 liên kết $\sigma$ đồng phẳng tạo góc $120^\circ$ và orbital $p_z$ tạo liên kết $\pi$ chịu trách nhiệm cho các tính chất điện tử năng lượng thấp. Luận án xây dựng chi tiết hai cách chọn ô cơ sở:

  1. Ô nguyên thủy hình thoi chứa 2 nguyên tử (A, B) với vectơ mạng $\vec{a}1, \vec{a}2$, dẫn đến vùng BZ hình lục giác có 6 đỉnh tiếp xúc Dirac chia thành hai điểm độc lập $K, K'$: $$\vec{K} = \left(\frac{2\pi}{3 a{\text{cc}}}, \frac{2\pi}{3\sqrt{3}a{\text{cc}}}\right), \quad \vec{K}' = \left(\frac{2\pi}{3 a_{\text{cc}}}, -\frac{2\pi}{3\sqrt{3}a_{\text{cc}}}\right)$$
  2. Ô cơ sở hình chữ nhật chứa 4 nguyên tử ($a, b, c, d$) với các tọa độ xác định trong mạng thực, thu gọn vùng BZ lục giác thành BZ hình chữ nhật thông qua phép "gấp vùng" (folding), đưa 6 điểm Dirac về 2 điểm không tương đương nằm trên trục $k_y$ tại tọa độ $(0, \pm k_{yK})$ với $k_{yK} = \frac{2\pi}{3\sqrt{3}a_{\text{cc}}}$. Hamiltonian ma trận kích thước $4 \times 4$ trong không gian động lượng cho 4 nhánh năng lượng: $$E = \pm t_{\text{cc}} \sqrt{1 + 4\cos^2\left(\frac{\sqrt{3}}{2}k_y a_{\text{cc}}\right) \pm 4\cos\left(\frac{\sqrt{3}}{2}k_y a_{\text{cc}}\right)\cos\left(\frac{3}{2}k_x a_{\text{cc}}\right)}$$

Về tính chất quang học, độ dẫn quang của graphene tại $0\text{ K}$ bao gồm thành phần Drude nội dải $\sigma_D(\omega) = \frac{4 E_F}{\pi\hbar}\frac{i}{\omega + i/\tau}$ và thành phần chuyển dời ngoại dải $\sigma_I(\omega)$. Trong miền tần số $\hbar\omega < 2E_F$, độ dẫn ngoại dải bị triệt tiêu do hiệu ứng khóa Pauli (Pauli blocking).

Đối với siêu mạng graphene cấu hình armchair (A-GSL) dưới thế tĩnh điện tuần hoàn chu kỳ $L = 3 N a_{\text{cc}}$ và biên độ $U_b$, luận án xây dựng Hamiltonian khối kích thước $4N \times 4N$. Khi tăng thế rào $U_b$ lên các giá trị $U_0, 2U_0, 3U_0$ (với $U_0 = \frac{2\pi\hbar v_F}{L}$), cấu trúc vùng năng lượng xuất hiện sự dịch chuyển tương đối của các nón Dirac và hình thành các điểm cố định (pinning points) giao nhau cách đều nhau nguyên lần $U_0$. Độ dẫn quang $\sigma_{xx}, \sigma_{yy}$ của GSL suy giảm trong khoảng năng lượng $\hbar\omega \in (0, U_b)$, đồng thời xuất hiện tính dị hướng nhẹ theo phương phân cực của photon.

Chương 2: Tính toán hàm điện môi trong gần đúng RPA và khảo sát các đặc trưng plasmon của graphene trong mô hình điện tử liên kết chặt với lân cận gần nhất

Chương 2 tập trung tính toán tường minh hàm phân cực $P(q,\omega)$ và hàm điện môi $\epsilon_{\text{RPA}}(q,\omega)$ ở giới hạn gần đúng lân cận gần nhất. Bằng phương pháp giải tích cho trường hợp pha tạp yếu tại $T = 0\text{ K}$, hàm phân cực được chia thành hai phần: phần đóng góp không pha tạp $P^{(0)}(q,\omega)$ và phần biến thiên do pha tạp $P^{(1)}(q,\omega)$. Tác giả thiết lập các tích phân giải tích xác định chính xác phần thực và phần ảo:

  • Miền $\text{Im}[P(q,\omega)] = 0$ được xác định là vùng không có kích thích đơn hạt (hạt - lỗ trống), nơi plasmon có thể tồn tại mà không bị phân rã Landau không đàn hồi.
  • Bằng việc giải phương trình $\text{Re}[\epsilon_{\text{RPA}}(q,\omega_p)] = 0$, luận án chứng minh rằng trong giới hạn bước sóng dài ($q \ll k_F$), mode Dirac plasmon tuân theo quy luật tán sắc căn bậc hai: $$\omega_p(q) = \sqrt{\frac{2 e^2 E_F}{\epsilon_0 \epsilon_r \hbar^2} q}$$
    NĂNG LƯỢNG (\hbar\omega)
             0                   2kF

Khi mở rộng tính toán số cho toàn bộ vùng BZ ở nhiệt độ $T > 0\text{ K}$, luận án chỉ ra:

  1. Nhiệt độ hữu hạn làm mờ ranh giới của phân bố Fermi-Dirac, dẫn đến sự mở rộng của vùng kích thích đơn hạt và làm tăng tốc độ phân rã plasmon $\gamma$.
  2. Khi nâng cao nồng độ hạt tải (mức Fermi $E_F$ dịch chuyển sâu vào dải dẫn), mặt Fermi đẳng năng mất tính tròn xoay và bị biến dạng tam giác (trigonal warping). Do đó, phổ tán sắc plasmon tính số trong toàn vùng BZ thể hiện tính bất đẳng hướng rõ rệt, phụ thuộc vào phương của vectơ truyền sóng $\vec{q}$ (so sánh giữa phương $\Gamma\text{-}K$ và $\Gamma\text{-}M$).

Chương 3: Các đặc trưng plasmon của graphene trong chế độ pha tạp cao

Chương 3 mở rộng mô hình TB sang gần đúng lân cận thứ hai (NNN) với Hamiltonian có các yếu tố đường chéo $t_k^{AA} = 2 t' [2\cos(k_x a_{\text{cc}}) \cos(\sqrt{3}k_y a_{\text{cc}}) + \cos(2k_x a_{\text{cc}})]$, tạo ra sự bất đối xứng electron - lỗ trống và làm biến dạng mạnh các nón Dirac quanh hai điểm $K, K'$. Luận án tính toán vận tốc nhóm $\vec{v}_g(k) = \frac{1}{\hbar}\nabla_k E(k)$ và mật độ trạng thái (DOS) ở các mức pha tạp cao ($E_F = 0,5\text{ eV}$ và $E_F = 1,0\text{ eV}$).

Phát hiện quan trọng của chương là việc khảo sát hàm phổ mất mát năng lượng electron $S(q,\omega)$ (EELS) và các không điểm của hàm điện môi $\epsilon_{\text{RPA}}(q,\omega)$ theo các phương tác động của vectơ sóng $\vec{q}$ (phương $O_x$ và $O_y$):

  • Trong chế độ pha tạp mạnh, ngoài nhánh Dirac plasmon truyền thống có năng lượng cao hơn, xuất hiện một nhánh plasmon thứ hai ở vùng năng lượng thấp.
  • Luận án phân tích hàm mật độ chuyển trạng thái $S_{q,\omega}(k_x, k_y)$ trên nền các đường đẳng năng và chứng minh: sự xuất hiện của mode plasmon mới là hệ quả kết hợp của tính bất đẳng hướng hình học trong từng nón Dirac và sự không tương đương giữa các trạng thái điện tử của hai thung lũng Dirac độc lập trong vùng Brillouin khi mức Fermi bị đẩy lên cao.
Đặc tính so sánh Graphene pha tạp thấp ($E_F \to 0$) Graphene pha tạp cao ($E_F \ge 0,5\text{ eV}$)
Mô hình điện tử phù hợp Mô hình Dirac tuyến tính (NN) Mô hình liên kết chặt lân cận thứ hai (NNN)
Hình học mặt Fermi Mặt nón tròn xoay đẳng hướng Bị biến dạng tam giác (Trigonal warping)
Đối xứng electron - lỗ trống Bảo toàn ($E(k) = -E(-k)$) Bị phá vỡ do số hạng $t'$ ($t' \approx 0,1 \div 0,2t$)
Số lượng mode plasmon 01 mode Dirac plasmon duy nhất 02 mode plasmon (xuất hiện thêm mode năng lượng thấp)
Đặc tính tán sắc plasmon Đẳng hướng ($\omega_p \propto \sqrt{q}$) Bất đẳng hướng mạnh (phụ thuộc phương $\vec{q}$)
Cơ chế chuyển dời chủ đạo Nội dải lân cận điểm Dirac Chuyển dời phức hợp nội thung lũng và liên thung lũng

Chương 4: Hàm điện môi có tính đến hiệu ứng trường địa phương. Áp dụng cho kích thích plasmon ứng với sự chuyển trạng thái giữa các điểm K trong graphene

Chương 4 nghiên cứu bài toán ở giới hạn bước sóng ngắn, nơi tính không đồng nhất của mật độ điện tử trong ô mạng không thể bị bỏ qua. Phép gần đúng RPA thông thường (xem trường tự hợp là trường trung bình vĩ mô) được mở rộng bằng cách đưa vào ma trận điện môi đầy đủ $\epsilon_{\vec{G},\vec{G}'}(\vec{q},\omega)$ với $\vec{G}, \vec{G}'$ là các vectơ mạng đảo: $$\epsilon_{\vec{G},\vec{G}'}(\vec{q},\omega) = \delta_{\vec{G},\vec{G}'} - V_{2D}(\vec{q}+\vec{G}) P_{\vec{G},\vec{G}'}(\vec{q},\omega)$$ Hằng số điện môi vĩ mô $\epsilon_M(\vec{q},\omega)$ có tính đến hiệu ứng trường địa phương (LFE) được xác định từ nghịch đảo của ma trận điện môi: $$\epsilon_M(\vec{q},\omega) = \frac{1}{\left[\epsilon^{-1}\right]_{\vec{0},\vec{0}}(\vec{q},\omega)}$$

Luận án áp dụng công thức mở rộng này cho trường hợp kích thích với vectơ sóng trao đổi lớn $\vec{q} \approx \vec{K} - \vec{K}'$, tương ứng với quá trình kích thích electron chuyển trạng thái giữa hai điểm Dirac khác nhau (inter-valley transition). Kết quả tính số hàm EELS dự đoán sự tồn tại của một mode kích thích tập thể liên thung lũng (inter-valley plasmon) có quy luật tán sắc tuyến tính theo vectơ sóng ($\omega_p \propto q$).

Kết quả và những đóng góp mới

  1. Đóng góp mới về mặt lý luận và học thuật:
    • Xác lập giới hạn vật lý của mô hình Dirac đối với hệ electron hai chiều trong graphene: chỉ ra rằng khi mức độ pha tạp tăng cao, tính bất đẳng hướng của các mặt đẳng năng $\pi$ chi phối trực tiếp đến hàm phân cực, làm mất tính đẳng hướng của hệ thức tán sắc plasmon.
    • Phát hiện và giải thích bản chất vật lý của một mode plasmon mới ở vùng năng lượng thấp trong graphene pha tạp mạnh. Luận án chứng minh rằng tính bất đẳng hướng của các nón Dirac kết hợp với sự không tương đương giữa hai thung lũng $K, K'$ trong vùng Brillouin là điều kiện cần và đủ để hình thành mode dao động tập thể này.
    • Xây dựng quy trình mở rộng lý thuyết trường tự hợp để tính toán ma trận điện môi RPA có xét đến hiệu ứng trường địa phương (LFE) cho graphene, kiểm chứng tính khả dĩ của mode inter-valley plasmon có quan hệ tán sắc tuyến tính ở giới hạn bước sóng ngắn.
  2. Đóng góp về mặt công bố khoa học:
    • 01 công trình trên tạp chí Physica E: Low-dimensional Systems and Nanostructures (2014) về ảnh hưởng của nhiệt độ, độ pha tạp và tính bất đẳng hướng của dải năng lượng lên đặc trưng plasmon.
    • 01 công trình trên tạp chí Physica Status Solidi (b) (2016) công bố phát hiện về mode plasmon mới trong chế độ pha tạp cao.
    • 01 công trình phối hợp trên Journal of Physics: Condensed Matter (2014) về độ dẫn quang của siêu mạng graphene dưới tác động của hệ điện cực tuần hoàn.
  3. Kiến nghị và giải pháp đề xuất:
    • Đề xuất khai thác các mode plasmon bất đẳng hướng và mode plasmon mới ở vùng tần số THz/hồng ngoại trong thiết kế các linh kiện nano-plasmonics có khả năng điều khiển định hướng truyền sóng bằng điện trường ngoài.
    • Kiến nghị tiếp tục phát triển thuật toán tính số LFE toàn phần cho các cấu hình pha tạp bất đối xứng và cấu trúc dị thể 2D.

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

  • Hạn chế được chỉ ra từ văn bản:
    • Lý thuyết nhiều hạt trong luận án chủ yếu dừng lại ở gần đúng pha ngẫu nhiên (RPA), vốn bỏ qua các tương tác tương quan - trao đổi bậc cao giữa các electron ngoài RPA (chưa xét hiệu ứng trường cục bộ động lực học $G(q,\omega)$).
    • Nghiên cứu tập trung vào mô hình màng graphene phẳng lý tưởng, chưa đưa vào tính toán tường minh ảnh hưởng của các sai hỏng mạng (defects), tạp chất tích điện ngẫu nhiên trong chất nền, và tương tác tán xạ electron - phonon đối với thời gian sống thực tế của plasmon.
    • Kết quả tính toán hàm điện môi có xét hiệu ứng trường địa phương (LFE) ở Chương 4 mới là các kết quả số ban đầu áp dụng cho chuyển dời giữa hai điểm $K$, cần được hoàn thiện thành chương trình tính toán tổng quát cho toàn bộ vùng BZ.
  • Hướng nghiên cứu tiếp:
    • Hoàn thiện chương trình tính số hàm điện môi có xét LFE cho toàn bộ không gian xung lượng và các trạng thái kích thích phức hợp.
    • Khảo sát ảnh hưởng của các cơ chế tán xạ thực tế (electron-phonon, độ nhám bề mặt) lên sự suy giảm và độ rộng vạch phổ của mode plasmon mới.
    • Mở rộng phương pháp sang các hệ siêu mạng graphene hai chiều phức tạp hơn và các hệ vật liệu hai chiều dị thể (van der Waals heterostructures).

Giá trị tham khảo

  • Đối tượng thụ hưởng:
    • Nghiên cứu sinh, học viên cao học và giảng viên: Thuộc các chuyên ngành Vật lý lý thuyết, Vật lý chất cô đặc, Vật lý tính toán và Khoa học vật liệu nano.
    • Kỹ sư và nhà nghiên cứu công nghệ: Trong các lĩnh vực quang điện tử (optoelectronics), linh kiện vi sóng/terahertz và cảm biến nano-plasmonics.
  • Các nội dung có giá trị tham khảo chuyên sâu:
    • Phương pháp giải tích tính hàm phân cực 2D: Hệ thống công thức chi tiết trong Chương 2 và các phần Phụ lục (biến đổi Fourier thế Coulomb 2D, tích phân hàm chồng chập trạng thái, tách phần thực/ảo của hàm phân cực RPA).
    • Kỹ thuật chọn ô cơ sở hình chữ nhật: Phương pháp biến đổi Hamiltonian và gấp vùng BZ từ ô hình thoi sang ô chữ nhật 4 nguyên tử (Chương 1 và Chương 3), rất hữu ích cho mô phỏng số dải năng lượng và siêu mạng.
    • Thuật toán tính độ dẫn Kubo và ma trận LFE: Khung tính toán độ dẫn quang tensơ cho siêu mạng GSL và ma trận hàm điện môi nghịch đảo $\epsilon_{\vec{G},\vec{G}'}^{-1}$ cho hệ không đồng nhất.

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

1. Tại sao graphene lại có ưu thế hơn kim loại truyền thống trong các ứng dụng nano-plasmonics?

Kim loại truyền thống có mật độ electron cố định, tần số plasmon nằm ở dải ánh sáng nhìn thấy hoặc tử ngoại, và có quãng đường truyền sóng ngắn do suy giảm nhiệt mạnh. Graphene có tần số dao động plasmon nằm trong vùng terahertz (THz) và hồng ngoại, có khả năng giam giữ ánh sáng ở kích thước nhỏ hơn bước sóng, và đặc biệt nồng độ hạt tải có thể điều khiển liên tục thông qua điện trường tĩnh điện bên ngoài mà không cần thay đổi cấu trúc vật liệu.

2. Giới hạn của mô hình Dirac tròn xoay trong việc mô tả plasmon của graphene là gì?

Mô hình Dirac chỉ mô tả chính xác trạng thái điện tử trong khoảng năng lượng hẹp sát điểm tiếp xúc ($E \approx 0$) với giả định mặt đẳng năng là hình nón tròn xoay lý tưởng. Khi graphene được pha tạp mạnh (mức Fermi $E_F$ lớn) hoặc năng lượng kích thích cao, cấu trúc vùng năng lượng thực tế bị biến dạng tam giác (trigonal warping) và mất tính đối xứng tròn xoay. Việc tiếp tục dùng mô hình Dirac sẽ bỏ qua tính bất đẳng hướng của phổ tán sắc plasmon và không thể mô tả được các mode kích thích mới.

3. Cơ chế vật lý nào tạo ra mode plasmon mới ở chế độ pha tạp cao trong Chương 3?

Mode plasmon mới xuất hiện ở vùng năng lượng thấp khi graphene bị pha tạp mạnh ($E_F \ge 0,5\text{ eV}$) là kết quả đồng thời của hai yếu tố: (i) tính bất đẳng hướng của các mặt đẳng năng bên trong từng nón Dirac do tương tác với các nguyên tử ở lân cận thứ hai (NNN), và (ii) sự không tương đương giữa các trạng thái điện tử tại hai thung lũng Dirac độc lập ($K$ và $K'$) trong vùng Brillouin khi mức Fermi nâng cao.

4. Hiệu ứng khóa Pauli (Pauli blocking) tác động như thế nào đến độ dẫn quang của graphene?

Hiệu ứng khóa Pauli ngăn cản các chuyển dời điện tử ngoại dải (inter-band transitions) từ dải hóa trị lên dải dẫn khi các trạng thái đích đã bị chiếm đóng bởi electron. Do đó, ở nhiệt độ $0\text{ K}$, độ dẫn quang ngoại dải của graphene pha tạp bị triệt tiêu hoàn toàn trong dải năng lượng photon $\hbar\omega < 2E_F$, và độ dẫn trong miền này chủ yếu do đóng góp Drude nội dải chi phối.

5. Vai trò của hiệu ứng trường địa phương (LFE) trong kích thích plasmon inter-valley là gì?

Khi xét kích thích plasmon có bước sóng ngắn ứng với sự chuyển dời giữa hai điểm Dirac $K$ và $K'$, vectơ sóng trao đổi $\vec{q}$ có độ lớn tương đương khoảng cách giữa hai thung lũng trên mạng đảo. Ở thang kích thước nguyên tử này, điện trường cục bộ dao động mạnh trong ô mạng cơ sở và không thể xem là điện trường trung bình vĩ mô đồng nhất. Việc tính đến ma trận điện môi đầy đủ có LFE cho phép mô tả chính xác thế chắn vi mô và dự đoán mode inter-valley plasmon có tán sắc tuyến tính.

Kết luận

Luận án đã thực hiện nghiên cứu lý thuyết và tính toán mô phỏng toàn diện về các tính chất động lực học của hệ electron hai chiều trong mạng tinh thể graphene thuần khiết và siêu mạng graphene. Công trình đã làm sáng tỏ giới hạn của mô hình Dirac, chứng minh bản chất phân cực bất đẳng hướng của phổ plasmon dưới ảnh hưởng của nhiệt độ và nồng độ pha tạp, đồng thời phát hiện sự hình thành của một mode plasmon năng lượng thấp mới trong chế độ pha tạp mạnh bắt nguồn từ tính bất đẳng hướng thung lũng. Bên cạnh đó, việc mở rộng lý thuyết hàm điện môi RPA có xét đến hiệu ứng trường địa phương (LFE) đã cung cấp cơ sở lý thuyết cho kích thích plasmon liên thung lũng bước sóng ngắn, đóng góp dữ liệu định lượng quan trọng cho định hướng ứng dụng graphene trong lĩnh vực quang điện tử và nano-plasmonics dải tần terahertz.