Tổng quan luận án

Tính cấp thiết và khoảng trống nghiên cứu

Graphene là vật liệu hai chiều đầu tiên được cô lập và chế tạo thành công trong phòng thí nghiệm (năm 2004 bởi nhóm nghiên cứu của Andre Geim và Konstantin Novoselov). Vật liệu này sở hữu cấu trúc mạng lục giác dạng tổ ong của các nguyên tử carbon với các liên kết cộng hóa trị $\sigma$ hình thành từ sự lai hóa $sp^2$ và các liên kết $\pi$ hình thành từ các orbital $p_z$. Cấu trúc vùng năng lượng của graphene có tính chất tuyến tính ở lân cận các điểm Dirac ($K$ và $K'$) và không có khe năng lượng, trong đó các hạt tải điện hành xử như các fermion Dirac hai chiều không khối lượng với vận tốc Fermi $v_F \approx 10^6\text{ m/s}$. Bên cạnh đó, quãng đường tự do trung bình của hạt tải trong graphene có thể đạt tới $l \sim 1\ \mu\text{m}$, trong khi kích thước các transistor nano thường chỉ khoảng $0{,}1\ \mu\text{m}$, biến graphene thành môi trường thuận lợi để thực hiện các quá trình truyền dẫn điện tử đạn đạo.

Tuy nhiên, việc ứng dụng graphene vào các cấu trúc và linh kiện nano gặp phải những rào cản vật lý đặc thù:

  1. Đối với chuyển tiếp lưỡng cực graphene (GBJ - Graphene Bipolar Junction): Các chuyển tiếp $p$-$n$, $n$-$p$-$n$ hay $p$-$n$-$p$ có thể được tạo ra bằng phương pháp trường tĩnh điện thông qua hệ thống các điện cực cổng đáy (back gate) và cổng trên (top gate). Trong các nghiên cứu lý thuyết trước đây, các rào thế chuyển tiếp thường được giả thiết có dạng chữ nhật hoặc hình thang. Mô hình thế chữ nhật bỏ qua độ rộng hữu hạn của vùng chuyển tiếp, trong khi các biểu thức giải tích của độ rộng vùng chuyển tiếp $D$ trong mô hình thế hình thang (như mô hình của Zhang và Fogler) bị phân kỳ khi mật độ hạt tải $n_1$ hoặc $n_2$ tiến về $0$. Do đó, việc xây dựng một mô hình rào thế trơn tru dạng Gauss mô tả sát thực tế hình học của điện trường cổng và mô tả liên tục các chế độ mật độ hạt tải mà không gặp hiện tượng kỳ dị là một khoảng trống nghiên cứu cần giải quyết.
  2. Đối với chấm lượng tử graphene (GQD - Graphene Quantum Dot): Do hiện tượng chui ngầm Klein xuất phát từ tính chirality của hàm sóng Dirac, các điện tử tới bờ thế tĩnh điện luôn có thể kết nối với các trạng thái lỗ trống mang năng lượng âm, dẫn đến việc giam cầm hạt tải hoàn toàn bằng thế tĩnh điện trong graphene là không thể thực hiện được (ngoại trừ một số điều kiện biên hoặc cấu hình đặc biệt). Hạt tải trong các chấm lượng tử graphene hình tròn (CGQD) tạo bởi thế tĩnh điện thường chỉ tồn tại ở các trạng thái giả liên kết (Quasi-Bound States - QBS) với thời gian sống hữu hạn. Việc xác định chính xác phổ năng lượng phức, thời gian giam cầm, mật độ trạng thái địa phương (LDOS) và sự phụ thuộc của chúng vào dạng thế xuyên tâm (thế hình thang, thế Lorentz, thế dạng hàm mũ) cũng như khi có từ trường ngoài đặt vào là vấn đề lý thuyết then chốt cần được làm rõ.

Mục tiêu nghiên cứu

Mục đích chính của luận án là nghiên cứu lý thuyết các tính chất điện tử và đặc trưng truyền dẫn điện của hai loại cấu trúc nano cơ bản trên nền graphene đơn lớp vô hạn:

  1. Nghiên cứu truyền dẫn điện tử đạn đạo qua chuyển tiếp lưỡng cực $n$-$p$-$n$ graphene với mô hình rào thế dạng Gauss; tính toán xác suất truyền qua, điện trở, dòng điện và shot noise.
  2. Nghiên cứu cấu trúc năng lượng, các trạng thái liên kết và giả liên kết, thời gian giam cầm hạt tải, mật độ trạng thái địa phương (LDOS) và tiết diện tán xạ trong các chấm lượng tử graphene hình tròn tạo bởi thế tĩnh điện đối xứng trục.
  3. Phát triển phương pháp ma trận truyền qua ($T$-ma trận) cho chấm lượng tử graphene hình tròn đặt trong từ trường đều vuông góc với mặt phẳng tấm graphene.
  4. Đề xuất phương pháp giải tích - số để tính trực tiếp mật độ trạng thái địa phương (LDOS) từ các hàm sóng chuẩn hóa cho thế xuyên tâm bất kỳ và đối chiếu với dữ liệu thực nghiệm hiển vi quét chui ngầm (STM).

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

  • Đối tượng nghiên cứu:
    • Chuyển tiếp lưỡng cực graphene (GBJ) cấu hình $n$-$p$-$n$ điều khiển bằng hai cổng tĩnh điện ($V_b, V_t$) với rào thế dạng Gauss.
    • Chấm lượng tử graphene hình tròn (CGQD) tạo bởi thế tĩnh điện đối xứng trục với các biên dạng thế khác nhau (thế hình thang, thế Lorentz, thế tổng quát) và CGQD trong từ trường đều vuông góc.
  • Phạm vi nghiên cứu: Luận án giới hạn trong phạm vi nghiên cứu lý thuyết và mô phỏng tính số trên cơ sở phương trình Dirac-Weyl hai chiều, gần đúng liên kết chặt bậc thấp cho electron năng lượng thấp quanh điểm Dirac, trong chế độ truyền dẫn đạn đạo ($L < l$) tại nhiệt độ thấp hoặc nhiệt độ phòng theo các điều kiện thực nghiệm được trích dẫn.

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

Các hướng nghiên cứu trước được tác giả điểm lại

Tác giả đã hệ thống hóa các công trình nghiên cứu lý thuyết và thực nghiệm trong và ngoài nước liên quan đến vật liệu carbon, graphene và các cấu trúc nano graphene:

  • Lịch sử cấu hình carbon và graphene: Radushkevich và Lukyanovich (1952) phát hiện ống carbon rỗng; Kroto, Smalley và Curl (1985) tìm ra fullerene $\text{C}_{60}$ (Giải Nobel Hóa học 1992); Iijima (1991) chế tạo ống nano carbon; Wallace (1947) tính toán cấu trúc vùng năng lượng của đơn lớp carbon bằng mô hình liên kết chặt; Boehm, Setton và Stumpp (1994) đề xuất thuật ngữ "graphene"; Geim và Novoselov (2004) tách thành công graphene đơn lớp từ graphite (Giải Nobel Vật lý 2010).
  • Hiện tượng chui ngầm Klein và tính chất Dirac: Oskar Klein (1929) phát hiện nghịch lý Klein trong phương trình Dirac tương đối tính; Katsnelson và cộng sự (2006) khảo sát chui ngầm qua rào thế vuông trong graphene và tìm ra điều kiện truyền qua hoàn toàn $T=1$ khi tới vuông góc hoặc khi thỏa mãn điều kiện cộng hưởng $q_x D = N\pi$; Cheianov và Fal'ko (2006) thiết lập công thức xác suất truyền qua cho bờ thế bước $T_{\text{step}}(\varphi) = \cos^2(\varphi)$ và bờ thế nghiêng $T(\varphi) = \exp(-\pi k_F d \sin^2\varphi)$.
  • Hiệu ứng Hall lượng tử dị thường: K. von Klitzing (1980) với hiệu ứng Hall lượng tử trong 2DEG bán dẫn; Zhang và cộng sự (2005), Novoselov và cộng sự (2005, 2007) ghi nhận hiệu ứng Hall lượng tử dị thường trong graphene ở nhiệt độ phòng với các bình nguyên bán nguyên $\sigma_{xy} = \pm 4(e^2/h)(n + 1/2)$ và mức Landau năng lượng $E_0 = 0$.
  • Nghiên cứu chuyển tiếp lưỡng cực graphene (GBJ): Özyilmaz và cộng sự (2007) đo truyền dẫn Hall lượng tử trong chuyển tiếp $n$-$p$-$n$; Huard và cộng sự (2007) khảo sát điện trở dọc $R(V_t)$ và ghi nhận sự bất đối xứng qua điểm Dirac; Gorbachev và cộng sự (2008) chế tạo chuyển tiếp $p$-$n$-$p$ đạn đạo sử dụng cổng treo "air-bridge"; Liu và Richter (2012) mô phỏng độ dẫn $G(V_t)$ bằng phương pháp liên kết chặt; Sonin (2009) tính toán hệ số Fano và độ dẫn cho rào thế bậc thang; Low và cộng sự (2009), Zhang và Fogler (2008) phân tích độ rộng chuyển tiếp hiệu dụng $D(n_1, n_2)$ có tính đến hiệu ứng chắn phi tuyến.
  • Nghiên cứu chấm lượng tử graphene (GQD): Zhao và cộng sự (2015) dùng mũi dò STM làm cổng trên để tạo chuyển tiếp $p$-$n$ hình tròn; Lee và cộng sự (2016) tạo GQD bằng điện tích sai hỏng trong lớp đế cách điện; Gutierrez và cộng sự (2016) chế tạo CGQD trên đế đồng; Hewageegana và Apalkov (2013) nghiên cứu trạng thái liên kết bền trong CGQD thế chữ nhật liên quan đến nghiệm hàm Bessel; Chen và cộng sự (2014) khảo sát thời gian giam cầm trong thế hàm mũ; De Martino và cộng sự (2007) đề xuất giam cầm bằng từ trường không đồng nhất; Giavaras và cộng sự (2009) khảo sát thế tĩnh điện kết hợp thế vector $V = V_0 r^s, A_\theta = A_0 r^t$; Maksym và cộng sự (2010) mô hình hóa hố thế tĩnh điện trong từ trường đều vuông góc.

Vị trí và khoảng trống luận án giải quyết

Mô hình/Vấn đề trước đây Hạn chế tồn tại Giải pháp và vị trí của luận án
Rào thế chữ nhật / hình thang cho GBJ Bỏ qua độ trơn tru thực tế của điện trường cổng; công thức độ rộng $D(n_1, n_2)$ bị phân kỳ ở mật độ điện tích thấp Đề xuất rào thế dạng Gauss mô tả trơn tru chuyển tiếp $n$-$p$-$n$, loại bỏ phân kỳ, tính toán đầy đủ phản xạ và giao thoa trong chế độ đạn đạo
Trạng thái giả liên kết (QBS) trong GQD Thiếu phương pháp chuẩn hóa nhất quán để tính trực tiếp LDOS cho thế xuyên tâm có hình dạng bất kỳ Phát triển phương pháp $T$-ma trận và kỹ thuật chuẩn hóa hàm sóng, tính chính xác LDOS, độ rộng mức và thời gian giam cầm $\tau$
So sánh lý thuyết - thực nghiệm GQD Các mô hình lý thuyết thường dùng thế bậc nhảy thô, khó so khớp định lượng với phổ STM Áp dụng mô hình thế thực tế (thế hình thang, thế Lorentz) và so sánh trực tiếp, khớp nối với dữ liệu STM của Gutierrez (2016)
CGQD trong từ trường đều Phương pháp $T$-ma trận chuẩn chưa được mở rộng đầy đủ cho CGQD đối xứng trục có từ trường vuông góc Xây dựng tường minh dạng hàm sóng và biểu thức $T$-ma trận cho bài toán CGQD trong từ trường đều vuông góc

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

Cơ sở lý thuyết và các phương trình chủ đạo

  • Mô hình liên kết chặt (Tight-Binding) và Phương trình Dirac-Weyl: Hệ thức tán sắc ở gần đúng bậc thấp của điện tử trong mạng tổ ong graphene: $$E(\mathbf{k}) = \pm t \sqrt{4\cos\left(\pi k_x a \sqrt{3}\right)\cos(\pi k_y a) + 4\cos^2(\pi k_y a) + 1}$$ với $t \approx 2{,}8\text{ eV}$ là tích phân nhảy, $a_0 \approx 1{,}42\text{ \AA}$, $a = a_0\sqrt{3} \approx 2{,}46\text{ \AA}$. Khai triển quanh các điểm Dirac ($\mathbf{k} = \mathbf{K} + \mathbf{q}$ với $|\mathbf{q}| \ll |\mathbf{K}|$), quan hệ tán sắc có dạng tuyến tính: $$E(\mathbf{q}) \approx \pm \hbar v_F |\mathbf{q}|$$ Phương trình Dirac hai chiều cho các giả spinor hai thành phần $\Psi(\mathbf{r}) = (\psi_A(\mathbf{r}), \psi_B(\mathbf{r}))^T$: $$\hbar v_F \boldsymbol{\sigma} \cdot \nabla \Psi(\mathbf{r}) = E\Psi(\mathbf{r})$$ với $\boldsymbol{\sigma} = (\sigma_x, \sigma_y)$ là các ma trận Pauli.
  • Tính chirality và toán tử Helicity: Toán tử helicity được định nghĩa bởi $\hat{h} = \boldsymbol{\sigma} \cdot \frac{\mathbf{q}}{|\mathbf{q}|}$, có các trị riêng $h = +1$ (electron) và $h = -1$ (lỗ trống).
  • Mức Landau trong graphene: Dưới từ trường ngoài vuông góc $B$, các mức Landau có năng lượng: $$E_n = \text{sign}(n) v_F \sqrt{2e\hbar B |n|}, \quad n = 0, \pm 1, \pm 2, \dots$$ và độ dẫn Hall lượng tử: $$\sigma_{xy} = g \frac{e^2}{h}\left(n + \frac{1}{2}\right) \quad (g = 4)$$

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

  1. Phương pháp ma trận truyền qua ($T$-ma trận):
    • Áp dụng giải phương trình Dirac-Weyl cho các cấu trúc 1 chiều (chuyển tiếp lưỡng cực) và 2 chiều đối xứng trục (chấm lượng tử hình tròn).
    • Chia miền thế biến thiên liên tục thành các lớp vi phân mỏng có thế hằng số từng đoạn, liên kết các hàm sóng tại các mặt phân cách thông qua điều kiện liên tục của spinor, từ đó xây dựng ma trận truyền qua toàn phần $T$.
  2. Phương pháp tính các đặc trưng truyền dẫn:
    • Xác suất truyền qua $T(E, \theta)$ hoặc $T(E, \varphi)$.
    • Dòng điện $I$, độ dẫn $G$, mật độ tạp âm hạt (shot noise) $S$, và hệ số Fano $F = S / (2eI)$.
    • Điện trở tổng cộng $R$, điện trở lẻ $R_{\text{odd}}$ theo điện thế cổng $V_t$ và $V_b$.
  3. Phương pháp tính trạng thái giả liên kết (QBS) và Mật độ trạng thái địa phương (LDOS):
    • Xác định năng lượng phức $E_i = \text{Re}(E_i) + i\text{Im}(E_i)$ từ các cực của ma trận $S$ hoặc các điều kiện triệt tiêu của thành phần ma trận $W$.
    • Thời gian giam cầm đặc trưng: $\tau = \hbar / |\text{Im}(E_i)|$.
    • Hàm mật độ trạng thái tổng cộng (TLDOS) và mật độ trạng thái địa phương (LDOS): $$g(E) = \frac{1}{\pi} \sum_i \frac{-\text{Im}(E_i)}{[E - \text{Re}(E_i)]^2 + [\text{Im}(E_i)]^2}$$
    • Phát triển thuật toán chuẩn hóa hàm sóng liên tục trong không gian thực để tính trực tiếp $\text{LDOS}(r, E)$.

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

Chương 1: Các tính chất điện tử của graphene

Chương 1 trình bày tổng quan lý thuyết cơ sở về vật liệu graphene đơn lớp và các cấu trúc nano dẫn xuất:

  • Cấu trúc tinh thể và vùng năng lượng: Mô tả mạng tổ ong hai chiều cấu thành từ hai mạng Bravais lục giác lồng nhau (mạng con A và B). Ô cơ sở gồm hai nguyên tử carbon với các vector cơ sở $\mathbf{a}_1, \mathbf{a}_2$ và vector mạng đảo $\mathbf{b}_1, \mathbf{b}_2$. Phân tích các quỹ đạo lai hóa $sp^2$ tạo liên kết $\sigma$ trong mặt phẳng và quỹ đạo $p_z$ tạo liên kết $\pi$ vuông góc mặt phẳng. Thiết lập hệ thức tán sắc tight-binding, chỉ ra sự suy biến tuyến tính tại các điểm Dirac $K$ và $K'$, xác định vận tốc Fermi $v_F \approx 10^6\text{ m/s}$.
  • Tính Chirality và Giả spin: Xây dựng hàm sóng dạng giả spinor hai thành phần tương ứng với xác suất có mặt của hạt trên hai mạng con. Định nghĩa toán tử helicity $\hat{h}$, chứng minh tính đối xứng electron - lỗ trống và sự liên kết chặt chẽ giữa giả spin và hướng xung lượng.
  • Hiện tượng chui ngầm Klein: Trình bày nghịch lý Klein trong điện động lực học lượng tử và cơ chế chui ngầm Klein trong graphene. Trích dẫn công thức Katsnelson cho rào thế vuông: $$T = \frac{\cos^2\varphi}{1 - \cos^2(q_x D)\sin^2\varphi}$$ với $q_x = \sqrt{(E - V_0)^2/(\hbar^2 v_F^2) - k_y^2}$. Giải thích nguyên nhân truyền qua hoàn toàn ($T=1$) khi góc tới $\varphi = 0$ và tại các đỉnh cộng hưởng $q_x D = N\pi$.
  • Hiệu ứng Hall lượng tử dị thường: Phân tích sự khác biệt giữa khí điện tử hai chiều (2DEG) trong bán dẫn thông thường (chu kỳ mức Landau $\hbar\omega_c$, các bình nguyên nguyên) và graphene (mức Landau phi tuyến $E_n \propto \text{sign}(n)\sqrt{|n|B}$, mức $E_0=0$ có độ suy biến một nửa, các bình nguyên bán nguyên $4(n+1/2)e^2/h$). Giải thích điều kiện quan sát hiệu ứng ở nhiệt độ phòng do khoảng cách mức Landau năng lượng thấp $\Delta_n \gg k_B T$.
  • Chuyển tiếp lưỡng cực graphene (GBJ): Khái quát cấu trúc GBJ tạo bởi hệ hai cổng tĩnh điện ($V_b, V_t$), các chế độ nồng độ hạt tải ($n$-$p$-$n$, $p$-$n$-$p$, $n$-$n'$-$n$, $p$-$p'$-$p$). Điểm lại các mô hình thế chữ nhật, mô hình thế hình thang (công thức Low và Zhang-Fogler) và các thí nghiệm đo điện trở dọc $R(V_t)$ của Huard, Gorbachev.
  • Chấm lượng tử graphene (GQD và CGQD): Trình bày các phương pháp thực nghiệm tạo chấm lượng tử bằng STM (Zhao), điện tích sai hỏng (Lee) và đế kim loại (Gutierrez). Phân tích bản chất trạng thái giả liên kết (QBS), công thức tính thời gian giam cầm $\tau = \hbar / |\text{Im} E_i|$ và mật độ trạng thái $g(E)$.
  • Giam cầm hạt tải bằng từ trường: Phân tích mô hình từ trường không đồng nhất (De Martino), mô hình kết hợp điện trường và từ trường $V = V_0 r^s, A_\theta = A_0 r^t$ (Giavaras), và mô hình hố thế tĩnh điện trong từ trường đều với năng lượng rơi vào khe giữa mức Landau 0 và $-1$ (Maksym).

Chương 2: Các phương pháp tính toán

Chương 2 thiết lập chi tiết các công cụ giải tích và tính toán số:

  • Phương pháp ma trận truyền qua ($T$-ma trận):
    • Khảo sát bài toán tán xạ một chiều qua rào thế $V(x)$. Chia rào thế thành $N$ đoạn nhỏ, mỗi đoạn có thế không đổi $V_j$.
    • Biểu diễn hàm sóng trong từng vùng dưới dạng tổ hợp tuyến tính của sóng truyền tới và sóng phản xạ.
    • Thiết lập ma trận liên hệ giữa hệ số sóng ở hai đầu rào thế, từ đó thu được ma trận $T$ toàn phần: $$\begin{pmatrix} A_{\text{out}} \ B_{\text{out}} \end{pmatrix} = T \begin{pmatrix} A_{\text{in}} \ B_{\text{in}} \end{pmatrix}$$
    • Trích xuất hệ số truyền qua $t$ và xác suất truyền qua $T(E, \theta) = |t|^2$.
  • Phương pháp tính các đặc trưng truyền dẫn điện tử:
    • Thiết lập công thức Landauer-Büttiker tính độ dẫn $G$: $$G = G_0 \int_{-\pi/2}^{\pi/2} T(E_F, \theta) \cos\theta, d\theta$$
    • Thiết lập công thức tính dòng điện tổng cộng $I(V)$ dưới điện áp đặt $V_{sd}$.
    • Thiết lập biểu thức tính toán mật độ tạp âm hạt (shot noise) $S$ và hệ số Fano $F$: $$F = \frac{S}{2eI} = \frac{\int T(1-T)\cos\theta, d\theta}{\int T\cos\theta, d\theta}$$
    • Phương pháp trích xuất điện trở dọc $R$ và thành phần điện trở lẻ $R_{\text{odd}}$ để so sánh trực tiếp với các phép đo thực nghiệm cổng tĩnh điện.

Chương 3: Chuyển tiếp n-p-n graphene với rào thế dạng Gauss

Chương 3 tập trung nghiên cứu cấu trúc chuyển tiếp lưỡng cực graphene $n$-$p$-$n$ sử dụng mô hình rào thế dạng Gauss:

  • Mô hình rào thế dạng Gauss:
    • Xây dựng dạng thế tĩnh điện trơn tru theo phương truyền dẫn $x$ dựa trên hình học thực tế của cổng trên có chiều dài $L$, điện thế cổng trên $V_t$, cổng đáy $V_b$, và hằng số điện môi lớp oxit cách điện.
    • Thế Gauss khắc phục triệt để sự gián đoạn của thế chữ nhật và sự phân kỳ tại điểm chuyển tiếp chế độ hạt tải của mô hình hình thang.
  • Xác suất truyền qua $T(\theta, E, L)$ và cộng hưởng Klein:
    • Tính toán và vẽ biểu đồ cực $T(\theta)$ trong khoảng góc $[-\pi/2, \pi/2]$ cho các bộ tham số $[L\text{ (nm)}, E\text{ (meV)}, V_b\text{ (V)}, V_t\text{ (V)}]$: $[25, 0, 60, -12]$, $[25, 50, 60, -12]$, $[50, 50, 60, -12]$, $[25, 0, 40, -6]$, $[25, 50, 40, -6]$, và $[50, 0, 40, -6]$.
    • So sánh trực tiếp giữa mô hình thế chữ nhật (đường nét đứt) và mô hình thế Gauss (đường nét liền): Thế Gauss cho thấy sự suy giảm mạnh của xác suất truyền qua ở các góc tới lớn, trong khi vẫn bảo toàn $T=1$ tại góc tới vuông góc $\theta = 0$ (chui ngầm Klein).
    • Phổ $T(E)$ theo năng lượng tới $E$ tại các góc tới khác nhau thể hiện các đỉnh cộng hưởng Fabry-Pérot hẹp do giao thoa lượng tử giữa hai bờ chuyển tiếp.
    • Phổ $T(L)$ theo chiều dài cổng trên $L$ cho thấy tính chất dao động tuần hoàn điều hòa của xác suất truyền qua khi hạt tải chuyển động đạn đạo trong vùng rào thế.
  • Đặc trưng điện trở, dòng điện và hệ số Fano:
    • Khảo sát sự phụ thuộc của điện trở $R$, điện trở lẻ $2R_{\text{odd}}$ và hệ số Fano $F$ vào điện thế cổng trên $V_t$ với ba giá trị thế cổng đáy $V_b = 40\text{ V}$ (đường đỏ gạch chấm), $60\text{ V}$ (đường xanh dương liền nét), và $80\text{ V}$ (đường xanh lá đứt nét).
    • Xác định chính xác vị trí điện thế cổng chuyển tiếp $V_t^{(C)}$ giữa chế độ $n$-$p$-$n$ và $n$-$n'$-$n$ ($V_t^{(C)} = -2{,}19\text{ V}$ tương ứng với các cấu hình $V_b$).
    • Điện trở $R$ đạt cực đại lân cận điểm Dirac cục bộ dưới cổng trên và thể hiện tính chất bất đối xứng rõ rệt đối với điểm cực đại khi $V_t$ quét qua vùng chuyển tiếp.
    • Hệ số Fano $F$ biến thiên nhạy cảm quanh điểm chuyển tiếp chế độ, phản ánh sự thay đổi cơ chế truyền dẫn từ chui ngầm lượng tử sang truyền dẫn đạn đạo thuần túy.

Chương 4: Chấm lượng tử graphene hình tròn tạo bởi thế tĩnh điện

Chương 4 khảo sát chi tiết cấu trúc điện tử và động học hạt tải trong chấm lượng tử graphene hình tròn (CGQD):

  • Xây dựng phương pháp $T$-ma trận cho CGQD:
    • Do tính đối xứng trục của thế giam cầm $V(r)$, phương trình Dirac-Weyl được phân rã trong hệ tọa độ cực $(r, \theta)$ với số lượng tử mômen toàn phần $j = \pm 1/2, \pm 3/2, \dots$ và số lượng tử chirality $\nu = \pm 1$.
    • Thiết lập ma trận chuyển tiếp $W$ và ma trận $T$ liên hệ giữa các nghiệm hàm sóng tại gốc $r \to 0$ và vùng tiệm cận $r \to \infty$.
    • Khảo sát ma trận $W$ trong trường hợp năng lượng $E = \bar{U} \pm \nu\Delta$.
  • Khảo sát trạng thái QBS và Tiết diện tán xạ trong thế hình thang:
    • Tính toán phổ các trạng thái liên kết và giả liên kết (QBS) cho CGQD tạo bởi thế xuyên tâm hình thang với độ cao thế $U_0 = 15$, độ dốc $\alpha = 0$, bán kính hiệu dụng $L$, số lượng tử $j = 3/2$, $\nu = +1$, $\Delta = 2$.
    • Phân tích sự phụ thuộc của năng lượng phức (gồm phần thực $\text{Re}(E)$ xác định vị trí mức và phần ảo $-\text{Im}(E)$ xác định độ rộng mức) khi độ dốc biên $\alpha$ biến thiên liên tục.
    • Tính toán tiết diện tán xạ vi phân năng lượng thấp theo góc tán xạ $\varphi$ và tiết diện tán xạ ngang $\eta(E)$ theo năng lượng tới $E$ với các giá trị $U_0 = 10, 20, 30$. Đồ thị làm rõ sự tương ứng hoàn hảo giữa các cực đại của tiết diện tán xạ và các đỉnh cộng hưởng trong mật độ trạng thái tổng cộng (TLDOS).
  • Phương pháp tính LDOS từ hàm sóng chuẩn hóa và thế Lorentz:
    • Xây dựng thuật toán tính toán trực tiếp mật độ trạng thái địa phương $\text{LDOS}(r, E)$ từ các hàm sóng liên tục được chuẩn hóa trên miền không gian thực, cho phép áp dụng cho thế xuyên tâm bất kỳ mà không bị giới hạn bởi dạng hàm giải tích.
    • Áp dụng cho CGQD tạo bởi thế dạng Lorentz: $V(r) = \frac{U_0}{1 + (r/R_0)^2}$ với $U_0 = 0{,}15\text{ eV}$ và bán kính $R_0 = 50\text{ nm}$.
  • Đối chiếu với thực nghiệm STM của Gutierrez và cộng sự (2016):
    • So sánh trực tiếp phổ $\text{LDOS}(E)$ tính toán được cho CGQD có bán kính $R = 5{,}43\text{ nm}$ (với mức năng lượng hiệu chỉnh $E_D = -0{,}347\text{ eV}$) với dữ liệu thực nghiệm STM trong tài liệu [10].
    • Đối chiếu các trạng thái ứng với $j = 1/2, 3/2, 5/2, 7/2$. Kết quả tính toán từ công thức chuẩn hóa của luận án khớp định lượng với dữ liệu đo STM về vị trí đỉnh năng lượng và độ rộng vạch phổ.
    • Phân tích độ rộng đỉnh cộng hưởng trích xuất từ thực nghiệm so với độ rộng lý thuyết $-\text{Im}(E)$, khẳng định độ chính xác của phương pháp xác định thời gian giam cầm hạt tải $\tau$.
  • Mở rộng phương pháp $T$-ma trận cho CGQD trong từ trường đều:
    • Thiết lập phương trình Dirac-Weyl có mặt từ trường vuông góc đồng nhất $B$ thông qua thế vector $\mathbf{A} = \frac{1}{2} B r \mathbf{e}_\theta$.
    • Xác định dạng giải tích của hàm sóng trong từ trường biểu diễn qua các hàm suy biến siêu bội (confluent hypergeometric functions).
    • Xây dựng biểu thức tường minh cho các phần tử của $T$-ma trận trong từ trường, mở đường cho việc tính toán phổ năng lượng và thời gian sống khi hệ chịu tác động đồng thời của thế tĩnh điện và từ trường ngoài.

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

Luận án đã đạt được các kết quả nghiên cứu và đóng góp khoa học cụ thể sau:

  1. Về mô hình chuyển tiếp lưỡng cực graphene (GBJ):

    • Đề xuất và ứng dụng thành công mô hình rào thế dạng Gauss để mô tả chuyển tiếp lưỡng cực $n$-$p$-$n$ trong chế độ truyền dẫn đạn đạo. Mô hình phản ánh sát cấu hình điện trường thực tế của cổng tĩnh điện, khắc phục triệt để nhược điểm phân kỳ toán học của mô hình thế hình thang tại các điểm chuyển tiếp chế độ nồng độ hạt tải.
    • Làm sáng tỏ cơ chế truyền dẫn qua rào thế Gauss: xác suất truyền qua đạt cực đại $T=1$ tại góc tới vuông góc do chui ngầm Klein và xuất hiện các đỉnh cộng hưởng hẹp do giao thoa sóng lượng tử.
    • Tính toán định lượng sự phụ thuộc của điện trở dọc $R$, điện trở lẻ $R_{\text{odd}}$ và hệ số Fano $F$ vào các điện thế cổng $V_b$ và $V_t$, giải thích rõ tính bất đối xứng của đường đặc trưng điện trở ghi nhận trong các thực nghiệm.

  2. Về lý thuyết chấm lượng tử graphene hình tròn (CGQD):

    • Hoàn thiện phương pháp $T$-ma trận cho các hệ hai chiều đối xứng trục, cho phép xác định chính xác phổ năng lượng phức, năng lượng vị trí mức $\text{Re}(E_i)$, độ rộng mức $-\text{Im}(E_i)$, và thời gian giam cầm hạt tải $\tau = \hbar / |\text{Im} E_i|$ trong các CGQD tạo bởi thế tĩnh điện.
    • Khảo sát chi tiết ảnh hưởng của hình dạng biên rào thế (độ dốc $\alpha$, độ sâu $U_0$, bán kính $L$) đến thời gian sống của các trạng thái giả liên kết (QBS), chỉ ra rằng độ trơn của bờ thế làm thay đổi đáng kể độ rộng mức và thời gian giam giữ electron.
    • Thiết lập mối liên hệ trực tiếp giữa tiết diện tán xạ ngang $\eta(E)$, tiết diện tán xạ vi phân và mật độ trạng thái tổng cộng TLDOS.

  3. Về phương pháp tính LDOS và đối chiếu thực nghiệm:

    • Đề xuất phương pháp mới tính trực tiếp mật độ trạng thái địa phương (LDOS) từ các hàm sóng chuẩn hóa trong không gian thực cho chấm lượng tử graphene hình tròn với thế xuyên tâm bất kỳ (đã kiểm chứng trên thế hình thang và thế Lorentz).
    • Đem lại sự phù hợp định lượng giữa kết quả tính toán lý thuyết và số liệu đo phổ hiển vi quét chui ngầm (STM) của nhóm Gutierrez và cộng sự (2016) trên chấm lượng tử graphene bán kính $5{,}43\text{ nm}$.

  4. Về chấm lượng tử graphene trong từ trường:

    • Phát triển thành công biểu thức giải tích và thuật toán $T$-ma trận cho CGQD đặt trong từ trường đều có phương vuông góc với mặt phẳng graphene, cung cấp công cụ tính toán cấu trúc năng lượng của hệ khi có sự kết hợp giữa giam cầm tĩnh điện và giam cầm từ.

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

Hạn chế của luận án

  • Các tính toán truyền dẫn và cấu trúc năng lượng được thực hiện trên mô hình graphene đơn lớp vô hạn lý tưởng, chưa xét đến các hiệu ứng tán xạ biên trong các dải nano graphene (graphene nanoribbons) có bề rộng hẹp.
  • Nghiên cứu tập trung chủ yếu vào chế độ truyền dẫn đạn đạo thuần túy ($L < l$) tại nhiệt độ không tuyệt đối hoặc nhiệt độ thấp, chưa tính đến sự tán xạ do mất trật tự mạnh, tán xạ phonon ở nhiệt độ phòng cao hoặc các hiệu ứng tương tác đa hạt Coulomb phức tạp giữa các electron trong chấm lượng tử.

Hướng nghiên cứu tiếp

  • Mở rộng phương pháp $T$-ma trận và mô hình thế trơn sang các cấu trúc graphene hai lớp (bilayer graphene) và các dị cấu trúc van der Waals hai chiều khác.
  • Nghiên cứu ảnh hưởng của thế biến dạng cơ học (strain engineering) kết hợp với thế tĩnh điện và từ trường lên động học hạt tải trong chấm lượng tử graphene.
  • Khảo sát chi tiết các hiệu ứng tương tác Coulomb (hiện tượng khóa Coulomb - Coulomb blockade) trong các chấm lượng tử graphene kích thước nhỏ hoạt động ở chế độ transistor đơn điện tử.

Giá trị tham khảo

Luận án là tài liệu chuyên khảo có giá trị học thuật cao trong lĩnh vực vật lý lý thuyết và vật lý chất rắn tính toán:

  • Đối với nghiên cứu sinh và học viên cao học chuyên ngành Vật lý lý thuyết, Vật lý chất rắn, Vật lý tính toán: Luận án cung cấp hệ thống phương pháp luận hoàn chỉnh và chi tiết về phương pháp ma trận truyền qua ($T$-ma trận) áp dụng cho phương trình Dirac-Weyl hai chiều, kỹ thuật giải bài toán tán xạ lượng tử và thuật toán tính toán các trạng thái giả liên kết (QBS) có năng lượng phức.
  • Đối với các nhà nghiên cứu trong lĩnh vực công nghệ nano và linh kiện bán dẫn hai chiều: Chương 3 và Chương 4 cung cấp các mô hình rào thế thực tế (thế Gauss, thế Lorentz), các công thức giải tích và kết quả tính số về điện trở, hệ số Fano, mật độ trạng thái địa phương (LDOS) và thời gian giam cầm hạt tải. Đây là cơ sở lý thuyết hữu ích để thiết kế, tối ưu hóa và phân tích kết quả đo đạc từ các linh kiện chuyển tiếp $n$-$p$-$n$ và chấm lượng tử graphene tạo bằng kỹ thuật cổng tĩnh điện hoặc đầu dò STM.

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

1. Tại sao hiệu ứng Hall lượng tử trong graphene lại có tính dị thường và có thể quan sát được ở nhiệt độ phòng?

Hiệu ứng Hall lượng tử trong graphene là dị thường vì các hạt tải điện hành xử như các fermion Dirac không khối lượng có hệ thức tán sắc tuyến tính. Khi đặt trong từ trường vuông góc $B$, các mức Landau có năng lượng phụ thuộc phi tuyến vào chỉ số mức: $E_n = \text{sign}(n) v_F \sqrt{2e\hbar B |n|}$. Mức Landau tại $n=0$ ($E_0 = 0$) có độ suy biến bằng một nửa các mức khác và được chia sẻ đồng thời bởi electron và lỗ trống, dẫn đến các bình nguyên độ dẫn Hall nhận các giá trị bán nguyên $\sigma_{xy} = \pm 4(e^2/h)(n + 1/2)$. Hiệu ứng này quan sát được ở nhiệt độ phòng vì khoảng cách giữa các mức Landau ở năng lượng thấp $\Delta_n = E_1 - E_0$ trong graphene rất lớn, thỏa mãn điều kiện $\Delta_n \gg k_B T$ ngay cả khi $T = 300\text{ K}$ dưới từ trường mạnh.

2. Ưu điểm nổi bật của mô hình rào thế dạng Gauss so với mô hình rào thế chữ nhật và hình thang trong chuyển tiếp n-p-n graphene là gì?

Mô hình rào thế chữ nhật bỏ qua độ rộng thực của miền chuyển tiếp và độ dốc điện trường, trong khi mô hình rào thế hình thang sử dụng biểu thức độ rộng vùng chuyển tiếp $D(n_1, n_2)$ bị phân kỳ toán học khi mật độ điện tích $n_1$ hoặc $n_2$ tiến về $0$. Mô hình rào thế dạng Gauss do luận án đề xuất khắc phục được cả hai nhược điểm này: mô tả đúng bản chất trơn tru của điện trường tạo bởi các điện cực cổng thực tế, loại bỏ hoàn toàn sự phân kỳ toán học tại điểm chuyển tiếp chế độ mật độ hạt tải, đồng thời cho phép tính toán chính xác các hiệu ứng giao thoa sóng lượng tử và chui ngầm Klein trong chế độ truyền dẫn đạn đạo.

3. Tại sao trong chấm lượng tử graphene tạo bởi thế tĩnh điện thông thường chỉ tồn tại các trạng thái giả liên kết (QBS) thay vì trạng thái liên kết bền vững?

Do hiện tượng chui ngầm Klein bắt nguồn từ tính chirality của phương trình Dirac hai chiều, khi electron bên trong chấm lượng tử va chạm với bờ rào thế tĩnh điện, luôn tồn tại các trạng thái lỗ trống mang năng lượng âm ở bên ngoài rào thế tạo thành kênh thoát liên tục cho hạt. Do đó, hạt tải không thể bị giam cầm vĩnh viễn mà chỉ tồn tại trong chấm lượng tử trong một khoảng thời gian hữu hạn ở các trạng thái giả liên kết (QBS) với mức năng lượng có phần ảo khác không ($E_i = \text{Re} E_i + i\text{Im} E_i$).

4. Thời gian giam cầm của hạt tải trong chấm lượng tử graphene được xác định như thế nào qua phổ năng lượng phức?

Thời gian giam cầm đặc trưng $\tau$ của một trạng thái giả liên kết (QBS) được xác định trực tiếp từ phần ảo của năng lượng phức theo hệ thức bất định lượng tử: $$\tau = \frac{\hbar}{|\text{Im}(E_i)|}$$ Trong đó, $\text{Re}(E_i)$ xác định vị trí đỉnh cộng hưởng trên thang năng lượng, còn $|\text{Im}(E_i)|$ tỷ lệ thuận với độ rộng nửa cực đại (FWHM) của đỉnh cộng hưởng trong hàm mật độ trạng thái địa phương (LDOS). Phần ảo càng nhỏ thì độ rộng đỉnh phổ càng hẹp và thời gian giam cầm hạt tải trong chấm lượng tử càng dài.

5. Kết quả tính toán mật độ trạng thái địa phương (LDOS) của luận án đã được kiểm chứng với dữ liệu thực nghiệm nào?

Kết quả tính toán lý thuyết về mật độ trạng thái địa phương (LDOS) từ các hàm sóng chuẩn hóa của luận án đã được đối chiếu trực tiếp với các số liệu thực nghiệm đo bằng kính hiển vi quét chui ngầm (STM) của nhóm tác giả Gutierrez và cộng sự (năm 2016, trích dẫn tài liệu [10]) trên chấm lượng tử graphene hình tròn bán kính $R = 5{,}43\text{ nm}$ với thế hiệu chỉnh Dirac $E_D = -0{,}347\text{ eV}$. Kết quả tính toán của luận án cho thấy sự tương thích chính xác về cả vị trí các mức năng lượng và độ rộng vạch phổ cho các trạng thái có số lượng tử $j = 1/2, 3/2, 5/2, 7/2$.

Kết luận

Luận án "Các tính chất truyền dẫn điện của một số cấu trúc nano graphene" của tác giả Nguyễn Thị Thuỳ Nhung đã giải quyết có hệ thống các bài toán lý thuyết về truyền dẫn điện tử đạn đạo trong chuyển tiếp lưỡng cực $n$-$p$-$n$ và cấu trúc năng lượng của chấm lượng tử graphene hình tròn. Bằng việc phát triển phương pháp ma trận truyền qua ($T$-ma trận), đề xuất mô hình rào thế dạng Gauss và phương pháp chuẩn hóa hàm sóng tính mật độ trạng thái địa phương (LDOS), luận án đã làm rõ cơ chế chui ngầm Klein, giao thoa lượng tử, các đặc trưng điện trở, hệ số Fano cũng như động học của các trạng thái giả liên kết (QBS). Các kết quả tính toán của luận án không chỉ giải quyết các hạn chế lý thuyết của các mô hình trước đó mà còn đạt được sự phù hợp định lượng với các dữ liệu thực nghiệm hiển vi quét chui ngầm (STM), đóng góp cơ sở khoa học quan trọng cho việc nghiên cứu và phát triển các linh kiện điện tử nano trên nền vật liệu graphene.