Tổng quan về luận án

Nghiên cứu các đặc trưng lượng tử và động học điện tử trong vật liệu hai chiều (2D) đang định hình lại ranh giới của vật lý chất cô đặc hiện đại và công nghệ nano bán dẫn. Kể từ khi Andre Geim và Konstantin Novoselov cô lập thành công graphene đơn lớp vào năm 2004 từ graphite, cấu trúc mạng tổ ong hai chiều cấu thành từ các liên kết lai hóa $sp^2$ với độ dài liên kết $a_0 \approx 1.42\text{ \AA}$ ($a = a_0\sqrt{3} \approx 2.46\text{ \AA}$) đã trở thành tâm điểm của vật lý thực nghiệm và lý thuyết. Bản chất lượng tử độc nhất của graphene bắt nguồn từ cấu trúc vùng năng lượng đặc biệt được mô tả lần đầu bởi P. R. Wallace (1947) qua mô hình liên kết chặt (tight-binding). Gần các điểm suy biến đối xứng $K$ và $K'$ tại biên vùng Brillouin thứ nhất (điểm Dirac), hệ thức tán sắc năng lượng không tuân theo phương trình Schrödinger phi tương đối tính thông thường mà thể hiện mối quan hệ tuyến tính hoàn hảo:

$$E(\mathbf{q}) \approx \pm \hbar v_F |\mathbf{q}|$$

Trong đó vận tốc Fermi $v_F = \frac{3 t a_0}{2\hbar} \approx 10^6\text{ m/s}$ (xấp xỉ $c/300$, với $t \approx 2.7\text{ eV}$ là năng lượng bước nhảy giữa các nguyên tử lân cận gần nhất). Các kích thích năng lượng thấp này hành xử như các fermion Dirac không khối lượng với tính chirality xác định và đối xứng giả spin (pseudospin).

Mạng tổ ong (Honeycomb) 
  ├── 2 mạng con A & B  ──> Giả spin (Pseudospin σ) ──> Chirality & Helicity
  └── Tán sắc tuyến tính: E = ±ħ v_F |q|
        ├── Điểm Dirac K, K' (EF = 0, không khe năng lượng)
        ├── Quãng đường tự do trung bình l ~ 1 µm (Truyền dẫn đạn đạo)
        └── Nghịch lý chui ngầm Klein (T = 1 tại góc tới trực diện φ = 0)

Khoảng trống nghiên cứu (research gap) cốt lõi xuất phát từ nghịch lý chui ngầm Klein (Klein tunneling) — một hệ quả trực tiếp của phương trình Dirac hai chiều được Oskar Klein dự đoán năm 1929 và được M. I. Katsnelson cùng V. V. Cheianov & V. I. Fal'ko xác lập năm 2006 cho graphene. Khi hạt tải Dirac đối mặt với một rào thế tĩnh điện dưới góc tới vuông góc ($\phi = 0$), xác suất truyền qua luôn đạt giá trị cực đại tuyệt đối ($T = 1$), bất kể chiều cao và độ rộng của rào thế. Tính chất này tạo ra một rào cản vật lý nghiêm ngặt: không thể giam cầm điện tử hoàn toàn bằng các rào thế tĩnh điện thông thường để tạo ra các trạng thái liên kết bền vững (bound states) như trong các chấm lượng tử bán dẫn cổ điển. Hơn nữa, trong các nghiên cứu lý thuyết về chuyển tiếp lưỡng cực graphene (Graphene Bipolar Junctions - GBJ), các mô hình truyền thống chủ yếu giả định rào thế dạng bậc thang vuông (rectangular) hoặc hình thang (trapezoidal). Các mô hình rào thế giải tích đơn giản này dẫn đến sự phân kỳ phi vật lý của điện trở tại các điểm chuyển tiếp chế độ mật độ hạt tải ($n_1, n_2 \to 0$), không phản ánh đúng cấu hình thế trơn thực tế đo được bằng thực nghiệm kính hiển vi quét chui ngầm (STM).

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 trực tiếp các điểm nghẽn lý thuyết này thông qua ba câu hỏi nghiên cứu (Research Questions - RQ) và giả thuyết khoa học (Hypotheses - H):

  • RQ1: Dạng phân bố không gian của rào thế tĩnh điện trơn (đặc biệt là thế dạng Gauss) ảnh hưởng như thế nào đến xác suất truyền qua, điện trở tiếp xúc, đặc trưng dòng - điện áp ($I\text{-}V$), shot noise và hệ số Fano trong chuyển tiếp $n\text{-}p\text{-}n$ graphene?
    • H1: Rào thế Gauss mô tả chính xác sự chuyển tiếp mật độ điện tích liên tục, loại bỏ hoàn toàn kỳ dị phân kỳ điện trở của mô hình rào góc cạnh và làm thay đổi căn bản phổ dao động cộng hưởng Fabry-Pérot của sóng điện tử Dirac.
  • RQ2: Cơ chế vi mô nào chi phối các trạng thái giả liên kết (Quasi-Bound States - QBS) và thời gian sống (confinement lifetime) $\tau$ của hạt tải trong chấm lượng tử graphene hình tròn (Circular Graphene Quantum Dots - CGQD) chịu tác động của các thế xuyên tâm trơn (hình thang, hàm mũ, Lorentz)?
    • H2: Độ dốc biên của thế giam cầm đóng vai trò quyết định: thế càng trơn thì thời gian giam cầm $\tau = \hbar/|\text{Im}(E_i)|$ càng được kéo dài do giảm thiểu tán xạ rò rỉ Klein ở các góc tới xiên.
  • RQ3: Sự kết hợp đồng thời giữa thế tĩnh điện đối xứng trục và từ trường đều vuông góc ($B \perp$ mặt phẳng graphene) tạo ra quy luật giam cầm năng lượng nào để cô lập hạt tải hoàn toàn khỏi các kênh rò rỉ liên tục?
    • H3: Khi độ sâu hố thế tĩnh điện $\Delta V$ nhỏ hơn khoảng cách giữa các mức Landau, năng lượng các trạng thái giam cầm sẽ rơi trọn vào dải khe năng lượng giữa mức Landau $0$ ($E_0 = 0$) và mức $-1$, thiết lập trạng thái giam cầm năng lượng bền vững.

Khung lý thuyết của luận án xây dựng trên nền tảng phương trình Dirac-Weyl 2D cho fermion phi khối lượng, lý thuyết ma trận truyền qua (Transfer Matrix Method - $T$-matrix) mở rộng cho các hệ đối xứng trục và liên tục, cùng hình thức luận tán xạ Landauer-Büttiker.

Phạm vi nghiên cứu tập trung vào màng graphene đơn lớp đạn đạo (ballistic transport regime) với quãng đường tự do trung bình $l \sim 1\text{ }\mu\text{m}$, kích thước cấu trúc nano từ $10\text{ nm}$ đến $100\text{ nm}$, dải thế cổng $V_b \in [40, 80]\text{ V}$, $V_t \in [-15, 0]\text{ V}$, khảo sát mật độ trạng thái địa phương (LDOS), năng lượng phức $E_i \in \mathbb{C}$, tiết diện tán xạ ngang $\eta$ và phổ truyền dẫn $T(E, \theta)$. Đóng góp của luận án mang tính nền tảng: giải quyết trọn vẹn bài toán động học truyền dẫn đạn đạo qua chuyển tiếp $n\text{-}p\text{-}n$ thế Gauss và cung cấp phương pháp giải tích - số trị tính toán trực tiếp LDOS từ hàm sóng chuẩn hóa, mở đường cho việc thiết kế các linh kiện nano graphene và qubit lượng tử.

Literature Review và Positioning

Bức tranh tổng quan học thuật về động học điện tử trong graphene được cấu thành từ 4 dòng nghiên cứu chính:

[Wallace 1947, Geim & Novoselov 2004] ──> Cấu trúc dải năng lượng Dirac & Chirality
                       │
                       ▼
[Klein 1929, Katsnelson 2006, Cheianov 2006] ──> Chui ngầm Klein & Tán xạ góc tới
                       │
        ┌──────────────┴──────────────┐
        ▼                             ▼
[GBJ / Ballistic Transport]    [CGQD / Giam cầm lượng tử]
- Huard (2007): Bất đối xứng R   - De Martino (2007): B-field QD
- Gorbachev (2008): Air-bridge - Hewageegana (2008): QBS thế vuông
- Zhang & Fogler (2008): Gap    - Gutierrez & Lee (2016): STM data
        │                             │
        └──────────────┬──────────────┘
                       ▼
    [LUẬN ÁN NGUYỄN THỊ THUỲ NHUNG]
    - Rào thế Gauss cho GBJ (Khử phân kỳ R)
    - Formalism tính LDOS chuẩn hóa cho CGQD
    - T-matrix mở rộng cho hệ tĩnh điện + Từ trường

1. Cấu trúc điện tử Dirac và tính chirality: Sau công trình tiên phong của Wallace (1947) và phát hiện thực nghiệm của Geim & Novoselov (2004, 2005), A. H. Castro Neto cùng các cộng sự (2009) đã hệ thống hóa các tính chất điện tử của graphene. Điểm mấu chốt là sự phân tách không gian mạng thành hai mạng con tam giác $A$ và $B$, tạo ra bậc tự do giả spin $\boldsymbol{\sigma} = (\sigma_x, \sigma_y)$ và toán tử helicity $\hat{h} = \boldsymbol{\sigma} \cdot \mathbf{q} / |\mathbf{q}|$. Toán tử này nhận trị riêng $h = +1$ cho electron vùng dẫn và $h = -1$ cho lỗ trống vùng hóa trị, bảo toàn tính chirality và triệt tiêu tán xạ ngược $180^\circ$ (Ando et al., 1998).

2. Hiện tượng chui ngầm Klein: Năm 2006, Katsnelson, Novoselov và Geim đã công bố công thức xác suất truyền qua qua bờ thế vuông bề rộng $D$, chiều cao $V_0$:

$$T(\phi) = \frac{\cos^2(\phi)}{1 - \cos^2(q_x D)\sin^2(\phi)}$$

Cùng năm, Cheianov và Fal'ko (2006) thiết lập định luật truyền qua cho thế bậc thang dốc $T(\phi) = \cos^2(\phi)$ và thế dốc tuyến tính $T(\phi) = \exp(-\pi k_F d \sin^2\phi)$, chứng minh rằng hạt tải luôn vượt qua rào thế tại góc tới trực diện $\phi = 0$.

3. Chuyển tiếp lưỡng cực graphene (GBJ): Các nhóm thực nghiệm của B. Özyilmaz et al. (2007), V. Huard et al. (2007), và R. V. Gorbachev et al. (2008) đã chế tạo thành công chuyển tiếp $p\text{-}n\text{-}p$ và $n\text{-}p\text{-}n$ đạn đạo sử dụng hệ thống cổng kép (back-gate $V_b$ và top-gate $V_t$). Huard et al. phát hiện sự bất đối xứng rõ rệt của điện trở $R(V_t)$ quanh điểm trung hòa điện tích. Về mặt lý thuyết, T. Low et al. (2009) và L. M. Zhang & M. M. Fogler (2008) đã mô hình hóa chuyển tiếp bằng rào thế hình thang với độ rộng vùng chuyển tiếp $D(n_1, n_2)$. Tuy nhiên, mô hình của Zhang & Fogler bộc lộ điểm yếu chết người: khi mật độ $n_1$ hoặc $n_2 \to 0$, tham số $D$ làm điện trở tính toán phân kỳ tới vô hạn — một kết quả mâu thuẫn trực tiếp với các quan sát thực nghiệm có giá trị điện trở hữu hạn.

4. Giam cầm trong chấm lượng tử graphene (GQD): P. Hewageegana và V. Apalkov (2008) chỉ ra rằng trong CGQD thế vuông, trạng thái liên kết thực sự chỉ xuất hiện tại các điều kiện cộng hưởng giao thoa cực kỳ hạn chế (khi chiều cao thế là nghiệm của hàm Bessel), còn lại hầu hết là trạng thái giả liên kết (QBS) có thời gian sống hữu hạn. A. De Martino et al. (2007) và C. Giavaras et al. (2009) đề xuất giam cầm bằng từ trường không đồng nhất và trường hỗn hợp. Gần đây hơn, J. Gutierrez et al. (2016, Nature Physics) và C. Lee et al. (2016) đã dùng đầu dò STM tạo ra các CGQD hình tròn bán kính $R \approx 50\text{ nm}$ trên chất nền và đo trực tiếp phổ LDOS, ghi nhận biên thế thực nghiệm luôn biến thiên trơn tru cỡ $\sim 1\text{ nm}$.

Vị thế và đóng góp vượt trội của luận án: Luận án định vị chính xác tại điểm giao thoa giữa hạn chế lý thuyết của mô hình Zhang-Fogler và nhu cầu giải thích dữ liệu STM chính xác cao của Gutierrez et al. (2016). Bằng cách đưa vào mô hình thế Gauss liên tục cho GBJ và phát triển công cụ ma trận $T$ vi phân giải tích cho thế xuyên tâm bất kỳ, luận án đã loại bỏ triệt để điểm kỳ dị phân kỳ điện trở của Zhang & Fogler, đồng thời đạt được sự trùng khớp định lượng hoàn hảo với dữ liệu phổ LDOS thực nghiệm của Gutierrez et al. mà không cần dựa vào các tham số điều chỉnh tùy tiện.

Đóng góp lý thuyết và khung phân tích

Đóng góp cho lý thuyết

Luận án mở rộng sâu sắc khuôn khổ điện động lực học tương đối tính ứng dụng trong vật lý chất rắn, mang lại 3 đóng góp lý thuyết căn bản:

                  ĐÓNG GÓP LÝ THUYẾT NỀN TẢNG
                               │
       ┌───────────────────────┼───────────────────────┐
       ▼                       ▼                       ▼
[Khử phân kỳ phi vật lý]  [Quy luật QBS & Sống sót] [Giam cầm mức Landau]
Rào thế Gauss cho GBJ:    Xác lập tau = hbar/|Im(E)|   Kết hợp hố thế tĩnh điện
Khử điểm kỳ dị n1,n2->0   LDOS tính từ hàm sóng        Delta V < Delta E_n
trong Zhang-Fogler (2008) chuẩn hóa cho thế trơn       triệt tiêu kênh rò rỉ Klein
  1. Hoàn thiện lý thuyết truyền dẫn qua chuyển tiếp lưỡng cực thực tế: Chứng minh rằng rào thế Gauss phản ánh chính xác bản chất sàng lọc tĩnh điện phi cục bộ của các điện cực nano. Luận án bác bỏ giả định thế góc cạnh thô sơ, chỉ ra rằng sự làm mịn rào thế dẫn đến việc chọn lọc góc truyền qua hẹp hơn, giải thích cơ chế suy giảm các đỉnh cộng hưởng Fabry-Pérot bậc cao.
  2. Chuẩn hóa khái niệm và định lượng trạng thái giả liên kết (QBS): Thiết lập mô hình toán học giải tích xác định chính xác năng lượng phức $E_i = \text{Re}(E_i) - i\text{Im}(E_i)$ của electron trong hố thế đối xứng trục. Luận án làm sáng tỏ rằng phần ảo $\text{Im}(E_i)$ không đơn thuần là tham số giải tích mà phản ánh trực tiếp xác suất phân rã lượng tử và bề rộng vạch phổ Lorentz $\Gamma = 2|\text{Im}(E_i)|$ trong hàm mật độ trạng thái: $$g(E) = \frac{1}{\pi} \sum_{i} \frac{|\text{Im}(E_i)|}{[E - \text{Re}(E_i)]^2 + [\text{Im}(E_i)]^2}$$
  3. Xác lập điều kiện giam cầm năng lượng tuyệt đối trong từ trường: Mở rộng lý thuyết Maksym et al. (2012), luận án chứng minh chặt chẽ rằng sự kết hợp giữa hố thế tĩnh điện đối xứng trục có độ sâu $\Delta V$ và từ trường ngoài vuông góc $B$ sẽ bẻ cong các mức năng lượng cục bộ. Khi $\Delta V < \Delta E_n$ (khoảng cách giữa các mức Landau), electron bị khóa hoàn toàn trong dải cấm giữa mức $n=0$ và $n=-1$, vô hiệu hóa hiện tượng rò rỉ Klein.

Khung phân tích độc đáo

Khung phân tích của luận án tích hợp liền mạch 4 trụ cột lý thuyết:

  • Phương trình Dirac-Weyl 2D không khối lượng: Mô tả động học của spinor 2 thành phần $\Psi = (\psi_A, \psi_B)^T$ dưới tác dụng của thế ngoài $V(\mathbf{r})$: $$\left[ v_F \boldsymbol{\sigma} \cdot \left( -i\hbar\boldsymbol{\nabla} + e\mathbf{A} \right) + V(\mathbf{r})\mathbb{I} \right] \Psi(\mathbf{r}) = E \Psi(\mathbf{r})$$
  • Phương pháp Ma trận truyền qua vi phân ($T$-Matrix Formalism): Đối với thế đối xứng trục $V(r)$, phương trình 2D được quy về hệ phương trình vi phân thuần nhất 1 chiều theo biến bán kính $r$ cho các thành phần spinor xuyên tâm ứng với số lượng tử mômen xung lượng toàn phần $j = \pm 1/2, \pm 3/2, \dots$: $$\frac{d}{dr} \begin{pmatrix} f_j(r) \ g_j(r) \end{pmatrix} = \mathcal{M}(r, E, j) \begin{pmatrix} f_j(r) \ g_j(r) \end{pmatrix}$$
  • Hình thức luận tán xạ sóng thành phần (Partial Wave Scattering Theory): Biểu diễn ma trận tán xạ $S$, ma trận chuyển tiếp $T$, tiết diện tán xạ vi phân $\frac{d\sigma}{d\phi}$ và tiết diện tán xạ ngang $\eta$ (transport cross-section): $$\eta(E) = \int_{-\pi}^{\pi} (1 - \cos\phi) \frac{d\sigma}{d\phi} d\phi$$
  • Lý thuyết vận tải lượng tử Landauer-Büttiker và Khung nhiễu Shot Noise: Tính toán trực tiếp độ dẫn điện $G$, cường độ nhiễu dòng $S_{noise}$ và hệ số Fano $F = S_{noise} / (2e I)$ ở chế độ đạn đạo không cân bằng.

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

┌─────────────────────────────────────────────────────────────────────────┐
│                      QUY TRÌNH TÍNH TOÁN & PHÂN TÍCH                    │
└─────────────────────────────────────────────────────────────────────────┘
                                     │
     ┌───────────────────────────────┴───────────────────────────────┐
     ▼                                                               ▼
[CHUYỂN TIẾP n-p-n THẾ GAUSS]                       [CHẤM LƯỢNG TỬ TRÒN CGQD]
  1. Phân đoạn vi phân thế Gauss                      1. Khai triển sóng theo góc: j = ±1/2, ±3/2...
  2. T-matrix từng bậc phẳng: T = ∏ T_i               2. Tích phân Dirac-Weyl 1D theo r
  3. Tính T(E, θ) qua liên tục spinor                 3. Giải năng lượng phức Det[M(E_i)] = 0
  4. Tích phân Landauer: G, R, I-V, Fano              4. Chiết xuất LDOS chuẩn hóa & so sánh STM

Thiết kế nghiên cứu

Luận án tuân thủ triết lý Chủ nghĩa duy thực phản biện (Critical Realism) kết hợp với Thực chứng luận lượng tử (Quantum Positivism): các trạng thái lượng tử vi mô (spinor, pha giao thoa, mức Landau) là thực tại vật lý khách quan, có thể truy xuất và kiểm chứng chính xác thông qua các đại lượng vĩ mô và vi mô quan sát được (phổ LDOS từ STM, độ dẫn $G$, hệ số Fano $F$).

Thiết kế nghiên cứu sử dụng phương pháp mô phỏng số lượng tử chính xác cao (high-precision quantum numerical computation) trên nền tảng giải tích giải phương trình vi phân Dirac. Hệ thống được mô hình hóa ở cấp độ đơn lớp nguyên tử, giải quyết trên không gian cấu hình liên tục với độ chia lưới không gian đạt dưới $0.05\text{ nm}$, đảm bảo hội tụ nghiệm số tuyệt đối.

Quy trình nghiên cứu rigorous

Quy trình thuật toán ma trận truyền qua ($T$-matrix) được chuẩn hóa qua 5 bước nghiêm ngặt:

  1. Rời rạc hóa rào thế liên tục: Một rào thế trơn bất kỳ $V(x)$ (dạng Gauss $V(x) = V_0 \exp(-x^2/2\sigma^2)$) hoặc thế xuyên tâm $V(r)$ (dạng Lorentz $V(r) = U_0 / [1 + (r/R_0)^2]$) được chia thành $N$ lớp mỏng vi phân ($N \ge 1000$). Trong mỗi lớp thứ $m \in [x_m, x_{m+1}]$, thế năng được coi là một hằng số $V_m = V(x_m)$.
  2. Thiết lập nghiệm cục bộ: Tại mỗi phân đoạn $m$, hàm sóng Dirac được biểu diễn bằng tổ hợp tuyến tính chính xác của các sóng phẳng spinor lan truyền tiến và lùi: $$\Psi_m(x, y) = \left[ A_m \begin{pmatrix} 1 \ s_m e^{i\theta_m} \end{pmatrix} e^{i q_{x,m} x} + B_m \begin{pmatrix} 1 \ -s_m e^{-i\theta_m} \end{pmatrix} e^{-i q_{x,m} x} \right] e^{i k_y y}$$ với $s_m = \text{sign}(E - V_m)$ và $q_{x,m} = \sqrt{[(E - V_m)/\hbar v_F]^2 - k_y^2}$.
  3. Khớp điều kiện biên liên tục: Tại mỗi ranh giới $x = x_m$, áp dụng điều kiện liên tục của spinor $\Psi_m(x_m) = \Psi_{m+1}(x_m)$, xây dựng ma trận truyền cục bộ $T_m$ kết nối $(A_{m+1}, B_{m+1})^T = T_m (A_m, B_m)^T$.
  4. Nhân ma trận tổng hợp: Ma trận truyền toàn phần của cấu trúc là tích có thứ tự của toàn bộ chuỗi ma trận thành phần: $$\mathcal{T}{total} = \prod{m=N}^{1} T_m = \begin{pmatrix} \mathcal{T}{11} & \mathcal{T}{12} \ \mathcal{T}{21} & \mathcal{T}{22} \end{pmatrix}$$
  5. Triệt tiêu sai số và kiểm chứng tính Unitarity: Xác suất truyền qua $T$ và phản xạ $R$ được trích xuất từ các phần tử của $\mathcal{T}{total}$: $$t = \frac{\mathcal{T}{11}\mathcal{T}{22} - \mathcal{T}{12}\mathcal{T}{21}}{\mathcal{T}{22}}, \quad r = -\frac{\mathcal{T}{21}}{\mathcal{T}{22}}, \quad T = 1 - |r|^2 = 1 - \left|\frac{\mathcal{T}{21}}{\mathcal{T}{22}}\right|^2$$ Tính bảo toàn xác suất $T(E, \theta) + R(E, \theta) = 1.0000000$ được kiểm tra nghiêm ngặt tại mọi điểm năng lượng và góc tới.

Data và phân tích

Toàn bộ các biểu thức tích phân động học được lập trình và xử lý số trị trên môi trường phần mềm tính toán chuyên dụng (Fortran/MATLAB) với độ chính xác kép (double precision):

  • Dòng điện tổng cộng qua chuyển tiếp: $$I = \frac{g e W}{h^2 v_F} \int_{-\infty}^{\infty} dE \left[ f(E, \mu_L) - f(E, \mu_R) \right] (E - U_L) \int_{-\pi/2}^{\pi/2} T(E, \theta) \cos\theta d\theta$$ với thừa số suy biến $g = 4$ ($g_s = 2$ cho spin, $g_v = 2$ cho thung lũng $K, K'$).
  • Độ dẫn tuyến tính $G$ và Điện trở $R$: $$G = \frac{g e^2 W}{h^2 v_F} |\mu_0 - U_L| \int_{-\pi/2}^{\pi/2} T(\theta) \cos\theta d\theta, \quad R = \frac{1}{G}$$
  • Điện trở lẻ đối xứng $2R_{odd}$: Trích xuất phần bất đối xứng qua công thức $2R_{odd} = R(V_t) - R(-V_t)$ để phân tích sự chênh lệch mật độ hạt tải hai phía tiếp xúc.
  • Hệ số Fano (Fano Factor): $$F = \frac{S_{noise}}{2 e I} = \frac{\int_{-\pi/2}^{\pi/2} T(\theta)[1 - T(\theta)] \cos\theta d\theta}{\int_{-\pi/2}^{\pi/2} T(\theta) \cos\theta d\theta}$$ Đặc trưng cho sự chuyển đổi giữa cơ chế dẫn Poisson ($F = 1$) và dẫn đạn đạo thuần nhất ($F \to 0$, hoặc $F = 1/3$ tại điểm Dirac cho cấu trúc ngắn rộng).

Phát hiện đột phá và implications

                   5 PHÁT HIỆN ĐỘT PHÁ CỦA LUẬN ÁN
                                 │
  ┌──────────────────────────────┼──────────────────────────────┐
  ▼                              ▼                              ▼
[1. Khử phân kỳ R]       [2. Lọc góc tới & Pha]        [3. Thời gian sống QBS]
Rào Gauss khử triệt để   Gauss thu hẹp phổ góc T(θ);   Thế trơn Lorentz/thang
điểm kỳ dị n1,n2 -> 0    triệt tiêu dao động giả       kéo dài tau = hbar/|Im(E)|
của Zhang-Fogler (2008)  của thế góc cạnh              so với thế bước vuông
  │                                                             │
  └──────────────────────────────┬──────────────────────────────┘
                                 ▼
                 ┌───────────────────────────────┐
                 ▼                               ▼
       [4. Khớp dữ liệu STM]          [5. B-field Confinement]
       Trùng khớp phổ LDOS của        Hố thế tĩnh điện + B vuông góc
       Gutierrez et al. (2016)        khóa e- vào khe Landau (0, -1),
       (R=50nm, U0=0.43eV)            triệt tiêu hoàn toàn Klein leak

Những phát hiện then chốt

  1. Giải quyết dứt điểm nghịch lý phân kỳ điện trở trong chuyển tiếp lưỡng cực: So sánh đối chiếu trực tiếp giữa mô hình rào thế vuông/hình thang và rào thế Gauss với các tham số điện áp thực nghiệm $[L(\text{nm}), E(\text{meV}), V_b(\text{V}), V_t(\text{V})] = [25, 0, 60, -12]$ và $[50, 50, 40, -6]$ chứng minh: rào thế Gauss mô tả quá trình chuyển tiếp liên tục giữa chế độ $n\text{-}p\text{-}n$ và $n\text{-}n_0\text{-}n$ tại các giá trị điện thế cổng tới hạn $V_t^{(C)} = -2.19\text{ V}$ (tương ứng $V_b = 40, 60, 80\text{ V}$). Đồ thị điện trở $R(V_t)$ và điện trở lẻ $2R_{odd}(V_t)$ trong mô hình Gauss hoàn toàn trơn tru, không xuất hiện các đỉnh phân kỳ vô hạn phi vật lý như trong lý thuyết của Zhang & Fogler (2008).
  2. Cơ chế chọn lọc góc truyền qua và triệt tiêu giao thoa Fabry-Pérot giả tạo: Khảo sát biểu đồ phân bố góc $T(\theta)$ chỉ ra rằng rào thế Gauss đóng vai trò như một bộ lọc góc (angular collimator) tự nhiên. Trong khi thế chữ nhật tạo ra vô số đỉnh cộng hưởng nhân tạo mỏng đạt $T = 1$ ở các góc tới xiên lớn do phản xạ sắc nhọn tại các cạnh đứng, thế Gauss làm suy giảm nhanh xác suất truyền qua khi $|\theta| > \pi/6$, tập trung truyền dẫn đạn đạo trong nón góc hẹp quanh trục vuông góc $\theta = 0$.
  3. Quy luật kéo dài thời gian sống $\tau$ của hạt tải theo độ trơn rào thế trong CGQD: Bằng việc giải phương trình trị riêng phi Hermite để tìm phổ năng lượng phức của các trạng thái giả liên kết (QBS), luận án phát hiện: khi chuyển từ thế giam cầm biên sắc nét (thế bậc thang $\alpha = 0$) sang thế trơn (thế hình thang nghiêng $\alpha > 0$ và thế Lorentz), phần ảo $|\text{Im}(E_i)|$ giảm rõ rệt. Điều này đồng nghĩa với việc độ rộng vạch cộng hưởng trong phổ mật độ trạng thái địa phương (LDOS) co hẹp lại, và thời gian giam cầm lượng tử $\tau = \hbar / |\text{Im}(E_i)|$ tăng lên đáng kể, chứng minh rằng rào thế trơn trong thực nghiệm có khả năng "bẫy" electron Dirac lâu hơn nhiều so với dự đoán lý thuyết cổ điển.
  4. Tái lập chính xác tuyệt đối phổ LDOS thực nghiệm STM của nhóm Gutierrez (2016): Luận án đã áp dụng phương pháp giải tích mới tính LDOS từ các hàm sóng chuẩn hóa cho chấm lượng tử graphene hình tròn có bán kính $R = 50\text{ nm}$, độ sâu thế $U_0 = 0.43\text{ eV}$ và mức hiệu chỉnh phông Dirac $E_D = -0.347\text{ eV}$. Kết quả tính toán của luận án cho các trạng thái với số lượng tử mômen toàn phần $j = 1/2, 3/2, 5/2, 7/2$ trùng khớp hoàn hảo với dữ liệu thực nghiệm đo bằng STM trong công trình của J. Gutierrez et al. (Nature Physics 2016), vượt trội hơn hẳn các công thức không chuẩn hóa trước đó.
  5. Thiết lập trạng thái giam cầm tuyệt đối bằng từ trường đều: Trong hệ CGQD đặt dưới từ trường vuông góc $B$, các mức năng lượng lượng tử hóa phân tách thành các mức Landau tương đối tính $E_n = \text{sign}(n) v_F \sqrt{2e\hbar B |n|}$. Luận án chứng minh rằng khi đặt một hố thế tĩnh điện nông có độ sâu $\Delta V < E_1 - E_0$, các mức QBS biến đổi hoàn toàn thành các trạng thái liên kết thực sự nằm cô lập trong khe năng lượng giữa mức $n=0$ và $n=-1$, không có bất kỳ trạng thái liên tục nào bên ngoài để chui ngầm rò rỉ.
Thông số / Hiện tượng Mô hình rào thế Vuông / Góc cạnh (Prior Art) Mô hình rào thế Gauss / Trơn (Luận án) Ý nghĩa vật lý & Thực nghiệm
Điện trở tại điểm tới hạn $V_t^{(C)}$ Phân kỳ phi vật lý ($R \to \infty$) khi $n_1, n_2 \to 0$ Biến thiên trơn tru, liên tục, hữu hạn tuyệt đối Phù hợp hoàn hảo với dữ liệu đo của Huard et al. (2007)
Phổ góc truyền qua $T(\theta)$ Đỉnh cộng hưởng giả tạo lan rộng ở góc lớn Bộ lọc nón góc hẹp ($T \to 0$ khi $|\theta| > \pi/6$) Loại bỏ phản xạ biên nhân tạo
Thời gian giam cầm hạt tải $\tau$ Cực ngắn do rò rỉ Klein ở mọi góc Tăng mạnh theo độ trơn thế ($\alpha > 0$, thế Lorentz) Cho phép bẫy hạt tải Dirac trong thời gian đủ dài
Mô phỏng phổ LDOS STM Lệch pha năng lượng và sai khác độ rộng vạch Khớp chính xác phổ thực nghiệm của Gutierrez (2016) Xác nhận tính đúng đắn của chuẩn hóa hàm sóng
Giam cầm lượng tử dưới từ trường Rò rỉ liên tục nếu chỉ dùng điện trường Giam cầm tuyệt đối trong khe Landau ($n=0 \to -1$) Nền tảng cho kiến trúc Qubit Graphene

Implications đa chiều

  • Về mặt lý thuyết vật lý chất cô đặc: Cung cấp chuẩn mực mới về việc mô hình hóa các rào thế nano trong vật liệu 2D, chứng minh rằng không thể bỏ qua hiệu ứng trơn hóa thế năng do sàng lọc tĩnh điện khi phân tích động học fermion Dirac.
  • Về mặt phương pháp luận tính toán: Thiết lập bộ công cụ thuật toán $T$-matrix vi phân chuẩn hóa, có thể chuyển giao trực tiếp để nghiên cứu các vật liệu Dirac và Weyl 2D/3D mới nổi như Silicene, Germanene, Phosphorene, và các chất cách điện Topo (Topological Insulators).
  • Về mặt công nghệ và thiết bị bán dẫn nano: Định hướng quy trình thiết kế transistor lưỡng cực graphene (Graphene Bipolar Transistors) vận hành ở chế độ đạn đạo siêu cao tần (THz) với mức tiêu tán năng lượng và hệ số nhiễu shot noise cực thấp.
  • Về mặt điện toán lượng tử (Quantum Computing): Đặt nền móng cho việc chế tạo các chấm lượng tử graphene giam cầm electron bằng từ trường để hiện thực hóa các qubit spin đơn electron có thời gian kết hợp pha (coherence time) dài nhờ tương tác spin-quỹ đạo yếu và không có đồng vị từ tính của hạt nhân Carbon-12.

Limitations và Future Research

  1. Mô hình đơn hạt độc lập (Single-particle approximation): Luận án giả định các electron chuyển động độc lập trong trường thế tĩnh điện hiệu dụng, chưa tính toán tường minh tương tác đa hạt electron-electron (Coulomb correlation) và hiệu ứng trao đổi (exchange interactions) vốn có thể làm tái chuẩn hóa vận tốc Fermi $v_F$ ở nồng độ hạt tải cực thấp.
  2. Giả thiết màng graphene đơn lớp hoàn hảo: Nghiên cứu tập trung vào graphene lý tưởng không biên (infinite sheet), chưa xét đến ảnh hưởng tán xạ vùng biên (armchair/zigzag edge states) và các khuyết tật mạng nguyên tử (vacancies, Stone-Wales defects) hay pha tạp hóa học ngẫu nhiên.
  3. Giới hạn nhiệt độ không tuyệt đối ($T = 0\text{ K}$): Phần lớn các dẫn xuất dòng và shot noise sử dụng phân bố hàm bước Fermi; ở nhiệt độ phòng ($T = 300\text{ K}$), tán xạ tán xạ nhiệt với phonon âm và phonon quang (electron-phonon scattering) sẽ làm giảm quãng đường tự do trung bình đạn đạo.
  4. Cấu hình từ trường đồng nhất: Luận án khảo sát từ trường ngoài $B$ đều vuông góc; chưa mở rộng cho các cấu hình vi từ trường không đồng nhất (inhomogeneous magnetic barriers) sinh ra từ các vi cấu trúc sắt từ bề mặt.

Chương trình nghiên cứu 10 năm tiếp theo (Future Research Agenda):

  • Hướng 1: Tích hợp phương pháp phiếm hàm mật độ kết hợp hàm Green không cân bằng (DFT-NEGF) để mô hình hóa tương tác đa hạt trong GBJ thế Gauss.
  • Hướng 2: Mở rộng mô hình sang hệ graphene hai lớp xoắn góc ma thuật (Magic-Angle Twisted Bilayer Graphene - MATBG) để khảo sát vận tải lượng tử trong điều kiện tương quan điện tử mạnh và siêu dẫn topo.
  • Hướng 3: Nghiên cứu động học spin trong chấm lượng tử graphene đa cổng kết hợp vi từ trường để thiết kế cổng logic lượng tử 2-qubit (CNOT gate).
  • Hướng 4: Khảo sát truyền dẫn đạn đạo ở nhiệt độ phòng có tính đến tương tác electron-phonon bằng phương pháp phương trình động học Boltzmann lượng tử.

Tác động và ảnh hưởng

  • Tác động học thuật sâu rộng: Các công trình công bố từ luận án cung cấp lời giải tham chiếu chuẩn mực cho cộng đồng vật lý lượng tử 2D, giải quyết các tranh cãi học thuật kéo dài về tính phân kỳ của điện trở chuyển tiếp và mở ra hướng tiếp cận giải tích chính xác cho phổ LDOS. Dự kiến các bài báo liên quan tạo lập chỉ số trích dẫn cao trong các tạp chí chuyên ngành hàng đầu (như Physical Review B, 2D Materials, Carbon).
  • Chuyển đổi công nghiệp bán dẫn và R&D nano: Cung cấp thông số thiết kế chính xác (chiều dài cổng $L$, độ dốc thế $\sigma$, điện áp cổng tối ưu) cho các phòng thí nghiệm R&D công nghệ cao (Intel, TSMC, Samsung) trong cuộc đua thu nhỏ linh kiện bán dẫn xuống dưới mốc $2\text{ nm}$, vượt qua giới hạn vật lý của silicon truyền thống.
  • Định hình chính sách công nghệ cao: Cung cấp luận cứ khoa học thực chứng cho các chiến lược phát triển vật liệu tiên tiến và công nghệ lượng tử quốc gia, khẳng định vai trò then chốt của vật liệu carbon cấu trúc nano trong kỷ nguyên hậu bán dẫn silicon.
  • Lợi ích kinh tế - xã hội: Góp phần hiện thực hóa các mạch tích hợp nano siêu nhanh, siêu tiết kiệm năng lượng, cắt giảm công suất tiêu thụ của các trung tâm dữ liệu toàn cầu và giảm thiểu phát thải nhiệt ra môi trường.

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

  • Nghiên cứu sinh & Nhà nghiên cứu trẻ (Doctoral & Postdoc Researchers): Tiếp cận một khung phương pháp luận ma trận $T$ vi phân hoàn chỉnh, minh bạch từ nền tảng giải tích đến mã nguồn tính số, giúp rút ngắn thời gian phát triển mô hình nghiên cứu các vật liệu 2D mới.
  • Các giáo sư & Chuyên gia vật lý lý thuyết (Senior Academics): Sở hữu một giải pháp lý thuyết triệt để giải quyết dứt điểm nghịch lý phân kỳ điện trở và công cụ tính toán LDOS chuẩn hóa để so sánh với các dữ liệu thực nghiệm tiên tiến nhất.
  • Kỹ sư R&D linh kiện nano bán dẫn (Semiconductor Industry R&D): Nắm vững các công thức định lượng về đặc tính $I\text{-}V$, hệ số Fano và góc truyền qua để tối ưu hóa thiết kế hình học điện cực cổng cho transistor đạn đạo.
  • Các nhà phát triển công nghệ điện toán lượng tử (Quantum Tech Developers): Ứng dụng cơ chế giam cầm hạt tải Dirac trong khe Landau để chế tạo các chấm lượng tử graphene có thời gian sống dài dùng làm qubit trạng thái rắn tin cậy.

Câu hỏi chuyên sâu

1. Đóng góp lý thuyết độc đáo nhất của luận án là gì và đã mở rộng lý thuyết nào?

Đóng góp lý thuyết độc đáo nhất là việc thiết lập mô hình rào thế Gauss liên tục cho chuyển tiếp lưỡng cực graphene (GBJ) và xây dựng phương pháp giải tích tính LDOS trực tiếp từ hàm sóng chuẩn hóa cho chấm lượng tử graphene hình tròn (CGQD). Luận án đã mở rộng trực tiếp lý thuyết truyền dẫn qua chuyển tiếp $p\text{-}n$ graphene của Cheianov & Fal'ko (2006) và khắc phục triệt để sai lầm phân kỳ phi vật lý của mô hình hình thang do Zhang & Fogler (2008) và Low et al. (2009) đề xuất, mang lại bức tranh vật lý nhất quán và chính xác tuyệt đối.

2. Phương pháp luận của luận án có điểm gì đột phá so với các nghiên cứu quốc tế trước đây?

So với phương pháp thế bước vuông của Hewageegana & Apalkov (2008) vốn chỉ tính được các trạng thái giả liên kết trên các bậc thế gián đoạn, và phương pháp gần đúng WKB của Katsnelson et al. (2006) vốn thất bại tại các điểm quay cổ điển, phương pháp Ma trận truyền qua vi phân ($T$-matrix) kết hợp chuẩn hóa hàm sóng toàn phần của luận án cho phép xử lý mọi dạng thế xuyên tâm trơn $V(r)$ bất kỳ. Phương pháp này giải trực tiếp hệ phương trình Dirac 1D mà không cần xấp xỉ thô, bảo toàn tính Unitarity của dòng hạt tải ở mức tuyệt đối ($T + R = 1$).

3. Phát hiện bất ngờ nhất trong quá trình nghiên cứu là gì?

Phát hiện bất ngờ nhất là rào thế Gauss đóng vai trò như một bộ lọc nón góc cực kỳ chọn lọc, triệt tiêu gần như hoàn toàn các đỉnh cộng hưởng Fabry-Pérot giả tạo ở góc tới xiên lớn. Trước đây, các mô hình thế vuông dự đoán vô số đỉnh cộng hưởng sắc nét ($T = 1$) ở các góc xiên. Luận án chứng minh rằng tính chất làm trơn biên thế Gauss đã "dập tắt" các phản xạ góc nhọn, phản ánh đúng dữ liệu đo thực tế không thấy các dao động điện trở nhiễu loạn này.

4. Luận án có cung cấp quy trình tái lập (Replication Protocol) minh bạch không?

Có. Toàn bộ hệ thống phương trình giải tích vi phân Dirac, biểu thức giải ma trận phân đoạn $T_m$, điều kiện biên liên tục của spinor, giải thuật tìm nghiệm phức $\det[\mathcal{M}(E_i)] = 0$ và công thức tích phân số Simpson cho dòng $I$, độ dẫn $G$, hệ số Fano $F$ đều được trình bày chi tiết từng bước trong Chương 2, Chương 3 và Chương 4, cho phép bất kỳ nhóm nghiên cứu độc lập nào cũng có thể lập trình tái lập chính xác 100% các kết quả số trị.

5. Lộ trình phát triển khoa học 10 năm tới được phác thảo như thế nào?

Lộ trình 10 năm tập trung vào 3 mũi nhọn: (1) Mở rộng formalism $T$-matrix sang các dị cấu trúc Van der Waals đa lớp xoắn góc ma thuật (Moiré superlattices); (2) Kết hợp trường biến dạng cơ học (strain engineering / pseudo-magnetic fields) với thế tĩnh điện để giam cầm hạt tải mà không cần từ trường ngoài; (3) Thiết kế và mô phỏng hoàn chỉnh kiến trúc chip xử lý lượng tử đơn electron dựa trên ma trận chấm lượng tử graphene.

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 đại diện cho một công trình nghiên cứu xuất sắc, chuẩn mực và đột phá trong lĩnh vực vật lý lý thuyết và vật lý tính toán chất cô đặc, với 5 đóng góp cốt lõi được khẳng định:

  1. Xác lập mô hình rào thế Gauss cho chuyển tiếp lưỡng cực graphene (GBJ): Khắc phục dứt điểm điểm kỳ dị phân kỳ điện trở của mô hình Zhang-Fogler, mô tả chính xác sự chuyển tiếp liên tục giữa các chế độ nồng độ hạt tải $n\text{-}p\text{-}n$, $p\text{-}n\text{-}p$, $n\text{-}n_0\text{-}n$.
  2. Khám phá quy luật lọc góc và dập tắt cộng hưởng giả tạo: Chứng minh rào thế Gauss giới hạn nón truyền dẫn đạn đạo trong phạm vi góc hẹp $|\theta| \le \pi/6$, phản ánh đúng cơ chế động học lượng tử thực nghiệm.
  3. Phát triển công cụ tính toán LDOS chuẩn hóa cho chấm lượng tử graphene hình tròn (CGQD): Tái lập chính xác tuyệt đối phổ đo quét chui ngầm STM của nhóm Gutierrez et al. (Nature Physics 2016) trên chấm lượng tử bán kính $R = 50\text{ nm}$.
  4. Định lượng thời gian giam cầm lượng tử của trạng thái giả liên kết (QBS): Chứng minh mối liên hệ tỷ lệ thuận giữa độ trơn của bờ thế giam cầm (thế hình thang dốc, thế Lorentz) và thời gian sống $\tau = \hbar / |\text{Im}(E_i)|$ của electron Dirac.
  5. Hoàn thiện cơ chế giam cầm hạt tải tuyệt đối bằng từ trường vuông góc: Chứng minh sự khóa chặt các mức năng lượng lượng tử vào dải khe Landau ($n = 0 \to -1$), mở đường cho công nghệ chế tạo linh kiện nano và qubit graphene hiệu năng cao.

Công trình không chỉ nâng tầm hiểu biết về điện động lực học Dirac hai chiều mà còn đặt nền móng lý thuyết vững chắc cho các ứng dụng công nghệ lượng tử và bán dẫn nano trong tương lai.