Giới thiệu dự án

Graphene – vật liệu hai chiều (2D) cấu thành từ mạng lục giác các nguyên tử carbon lai hóa $sp^2$ – sở hữu các đặc tính vật lý vượt trội như độ bền cơ học gấp 200 lần thép, độ dẫn nhiệt đạt xấp xỉ $5000\text{ W}\cdot\text{m}^{-1}\cdot\text{K}^{-1}$, và độ linh động điện tử ở nhiệt độ phòng vượt mức $250.000\text{ cm}^2\cdot\text{V}^{-1}\cdot\text{s}^{-1}$. Mặc dù có tiềm năng to lớn trong ngành công nghiệp quang điện tử (optoelectronics) trị giá hơn 45 tỷ USD, việc ứng dụng graphene đơn lớp vào các thiết bị thu nhận ánh sáng (photodetectors, pin mặt trời, cảm biến quang phổ) gặp phải rào cản vật lý cốt lõi: khả năng hấp thụ quang học tự nhiên trong vùng khả kiến của một lớp graphene độc lập chỉ đạt mức xấp xỉ $\pi\alpha \approx 2,3%$ (với $\alpha \approx 1/137$ là hằng số cấu trúc tinh tế). Mức hấp thụ $2,3%$ đồng nghĩa với hơn $97,7%$ năng lượng photon bị phản xạ hoặc truyền qua mà không chuyển đổi thành tín hiệu hữu ích.

                  Cấu trúc Vi Hốc Cộng Hưởng Quang Học
  Ánh sáng tới I(λ = 632.8 nm)
       =======>       |  Gương Au 1  |      Khoang rỗng      | Lớp Graphene |      Khoang rỗng      |  Gương Au 2  |
                      |  (t_Au1=30nm)|        (t / 2)        |  (Monolayer) |        (t / 2)        |  (t_Au2=50nm)|
  Phản xạ R <=======  |--------------|-----------------------|--------------|-----------------------|--------------|
                      | Region (1)   |      Region (2)       |   z_grap     |      Region (3)       | Region (4)   | =====> Truyền qua T ≈ 0

Vấn đề kỹ thuật đặt ra là làm thế nào để bẫy và tập trung trường điện từ quang học tại mặt phẳng nguyên tử của graphene mà không cần chế tạo các cấu trúc đa lớp điện môi cồng kềnh, phức tạp và đắt đỏ. Đề tài "Độ hấp thụ quang học của lớp graphene bên trong một vi hốc cộng hưởng" giải quyết bài toán này thông qua mô hình vi hốc quang học Fabry-Pérot bất đối xứng được tạo bởi hai màng mỏng kim loại vàng ($Au$), định vị chính xác tấm graphene đơn lớp tại tâm hốc cộng hưởng.

Mục tiêu cụ thể của dự án:

  1. Thiết lập hệ khung lý thuyết hoàn chỉnh dựa trên hệ phương trình Maxwell và điều kiện biên dòng điện mặt 2D ($\mathbf{J} = \sigma_{grap}\mathbf{E}$) để mô hình hóa tương tác photon - graphene.
  2. Xây dựng thuật toán ma trận chuyển dời (Transfer Matrix Method - TMM) mô phỏng chính xác sự lan truyền sóng điện từ qua hệ đa lớp kim loại - điện môi - graphene.
  3. Tối ưu hóa cấu hình hình học của vi hốc ($t_{Au1} = 30\text{ nm}$, $t_{Au2} = 50\text{ nm}$, bước sóng nguồn $\lambda = 632,8\text{ nm}$) nhằm cực đại hóa độ hấp thụ chọn lọc trên tấm graphene đơn lớp.
  4. Phân tích phân bố không gian của cường độ điện trường cục bộ bên trong các miền vi hốc để giải thích cơ chế vật lý của hiện tượng cộng hưởng quang học.

Kết quả kỳ vọng định lượng:

  • Tăng độ hấp thụ quang học của graphene từ mức nền $2,3%$ lên trên $25%$ (tăng trưởng hơn $11$ lần).
  • Tổng độ hấp thụ của toàn hệ vi hốc đạt mức cực đại xấp xỉ $76%$.
  • Cường độ điện trường laser tại vị trí đặt graphene được khuếch đại gấp nhiều lần so với biên độ sóng tới.

Phạm vi và giới hạn nghiên cứu: Nghiên cứu tập trung vào mô hình bán cổ điển 1D đối với sóng phẳng phân cực đơn sắc chiếu vuông góc ($\theta = 0^\circ$) ở bước sóng laser $\text{He-Ne}$ chuẩn $\lambda = 632,8\text{ nm}$. Các hiệu ứng phi tuyến bậc cao và tán xạ nhiệt độ không nằm trong phạm vi khảo sát sơ cấp.

Phân tích và thiết kế giải pháp

Phân tích hiện trạng

Trước đề tài này, các phương pháp nâng cao độ hấp thụ quang của graphene chủ yếu dựa trên tinh thể quang tử (Photonic Crystals) hoặc cấu trúc Bragg đa lớp phản xạ phân bố (DBR).

Tiêu chí so sánh Tấm Graphene Đơn Tự Do Cấu trúc Đa Lớp Điện Môi DBR (ABABGBABABA) Tinh Thể Quang Tử Dẫn Cộng Hưởng (Guided Resonance) Mô hình Vi Hốc Vàng Bất Đối Xứng (Dự án này)
Độ hấp thụ trên Graphene $2,3%$ $50,0%$ $\sim 85 - 99%$ $27,02%$
Số lượng lớp vật liệu 1 lớp 11 - 15 lớp Cấu trúc nano khắc lỗ 2D 3 lớp (Au / Air / Au)
Độ phức tạp chế tạo Rất thấp Rất cao (yêu cầu tỷ số $\delta = n_A/n_B \ge 3,56$) Cực kỳ cao (E-beam lithography) Thấp (Bốc bay nhiệt/phún xạ kim loại)
Bề dày toàn phần $0,34\text{ nm}$ Hàng micromet ($\sim 3 - 5,\mu\text{m}$) $\sim 500\text{ nm}$ $320\text{ nm} - 950\text{ nm}$
Độ nhạy sai số chế tạo Không Rất nhạy với sai lệch bề dày $1/4\lambda$ Rất nhạy với độ nhám và kích thước lỗ Dung sai hình học cao

Phân loại yêu cầu hệ thống theo khung MoSCoW:

  • Must-have (Bắt buộc): Giải chính xác nghiệm Maxwell cho sóng phẳng qua biên dẫn 2D; tích hợp hàm điện môi phụ thuộc tần số của kim loại theo mô hình Drude-Lorentz; thuật toán TMM tính toán hệ số $R$, $T$, $A$.
  • Should-have (Nên có): Khảo sát quét biến thiên khoảng cách hốc $t$ từ $0$ đến $1200\text{ nm}$; vẽ lát cắt phân bố điện trường $E(z)$ tại các điểm cộng hưởng và phản cộng hưởng.
  • Could-have (Có thể có): Mở rộng tính toán cho góc tới xiên ($\theta > 0^\circ$) và phân tách trạng thái phân cực $s$ và $p$.
  • Won't-have (Chưa thực hiện): Mô phỏng 3D FDTD toàn phần có xét đến khuyết tật biên tinh thể của graphene.

Thiết kế hệ thống

Kiến trúc mô phỏng hệ thống gồm 5 phân miền không gian dọc theo trục truyền sóng $z$:

  1. Miền 0 ($z < 0$): Không khí bán vô hạn ($n_0 = 1$), chứa sóng tới $E_i^{(0)}$ và sóng phản xạ $E_r^{(0)}$.
  2. Miền 1 ($0 \le z \le t_{Au1}$): Lớp gương vàng thứ nhất có bề dày $t_{Au1} = 30\text{ nm}$, hằng số điện môi phức $\varepsilon_{Au}(\omega)$.
  3. Miền 2 ($t_{Au1} < z < t_{Au1} + t/2$): Khoảng rỗng không khí phía trước tấm graphene.
  4. Mặt phẳng $z = t_{Au1} + t/2$: Lớp graphene đơn nguyên tử với độ dẫn mặt $\sigma_{grap} = e^2 / (4\hbar)$.
  5. Miền 3 ($t_{Au1} + t/2 < z < t_{Au1} + t$): Khoảng rỗng không khí phía sau graphene.
  6. Miền 4 ($t_{Au1} + t \le z \le t_{Au1} + t + t_{Au2}$): Lớp gương vàng phản xạ đáy có bề dày $t_{Au2} = 50\text{ nm}$.
  7. Miền 5 ($z > t_{Au1} + t + t_{Au2}$): Không khí truyền qua ($n_5 = 1$).
       +-----------------------------------------------------------------------+
       |                          HỆ THỐNG MÔ PHỎNG TMM                        |
       +-----------------------------------------------------------------------+
                                           |
                  +------------------------+------------------------+
                  |                                                 |
       +--------------------+                            +--------------------+
       |  Mô-đun Vật Lý 2D  |                            |   Mô-đun Kim Loại  |
       |  (Graphene Sheet)  |                            |   (Drude-Lorentz)  |
       | sigma = e^2 / 4hbar|                            | eps_Au(omega)      |
       +--------------------+                            +--------------------+
                  |                                                 |
                  +------------------------+------------------------+
                                           |
                               +-----------------------+
                               | Ma Trận Chuyển Dời M  |
                               | M = D01.P1.D12...D45  |
                               +-----------------------+
                                           |
                  +------------------------+------------------------+
                  |                                                 |
       +--------------------+                            +--------------------+
       |  Phổ Quang Học     |                            |  Cường Độ Trường   |
       |  R(t), T(t), A(t)  |                            |  E(z) Wave Profile |
       +--------------------+                            +--------------------+

Ngăn xếp công nghệ sử dụng:

  • Ngôn ngữ cốt lõi: Python 3.9.7
  • Thư viện tính toán đại số ma trận: NumPy 1.21.5
  • Thư viện xử lý tích phân hàm trường: SciPy 1.7.3 (scipy.integrate.quad)
  • Trực quan hóa dữ liệu quang học: Matplotlib 3.4.3
  • Môi trường đối chuẩn giải tích: Wolfram Mathematica 12.3

Hàm điện môi của vàng được chuẩn hóa qua mô hình mở rộng Drude: $$\varepsilon_{Au}(\omega_{opt}) = \varepsilon_b - \frac{\omega_{pl}^2}{\omega_{opt}^2 + i\omega_{opt}\gamma}$$ Trong đó hằng số điện môi nền $\varepsilon_b = 12$, tần số plasma $\hbar\omega_{pl} = 8,95\text{ eV}$, và hệ số cản dịu (damping constant) $\hbar\gamma = 72,3\text{ meV}$.

Methodology

Dự án áp dụng phương pháp nghiên cứu tính toán giải tích kết hợp mô phỏng số thông qua ma trận truyền qua (Transmission matrix $D$) và ma trận truyền pha (Propagation matrix $P$):

       +----------------------------------------------------------------------+
       | Phase 1: Xây dựng nghiệm giải tích Maxwell & điều kiện biên 2D (W1-2)|
       +----------------------------------------------------------------------+
                                          |
       +----------------------------------------------------------------------+
       | Phase 2: Lập trình hóa TMM & kiểm chuẩn độ hấp thụ nền 2.3% (W3-4)   |
       +----------------------------------------------------------------------+
                                          |
       +----------------------------------------------------------------------+
       | Phase 3: Quét tham số hốc cộng hưởng t in [0, 1200 nm] (W5-6)        |
       +----------------------------------------------------------------------+
                                          |
       +----------------------------------------------------------------------+
       | Phase 4: Tích phân tổn hao Joule & ánh xạ trường điện từ E(z) (W7-8) |
       +----------------------------------------------------------------------+
Rủi ro kỹ thuật Mức độ Chiến lược giảm thiểu
Tràn số hoặc mất ổn định số học trong ma trận truyền pha $P(z)$ với môi trường suy hao mạnh ($Au$) Trung bình Sử dụng ma trận tán xạ (Scattering Matrix $S$) để đối chứng nếu thành phần suy hao cục bộ phân kỳ.
Bỏ sót phần hấp thụ nội tại của 2 màng kim loại vàng dẫn đến tính sai phần hấp thụ của graphene Cao Tách biệt độ hấp thụ qua tích phân mật độ tổn hao điện trường: $A_{Au} = \frac{\omega}{c}\text{Im}(\varepsilon_{Au})\int
Lệch pha sóng đứng làm cực đại điện trường không rơi vào vị trí tấm graphene Cao Cố định tọa độ graphene luôn nằm tại vị trí đối xứng hình học $z = t_{Au1} + t/2$.

Implementation và kết quả

Development process

Cốt lõi của thuật toán dựa trên việc biểu diễn liên tục của các thành phần tiếp tuyến của điện từ trường qua mặt phân cách chứa graphene. Điều kiện biên cho trường từ tại mặt phẳng graphene ($z = 0$): $$\mathbf{e}_z \times (\mathbf{E}_2 - \mathbf{E}_1) = 0$$ $$\mathbf{e}_z \times (\mathbf{H}_2 - \mathbf{H}1) = \mathbf{J}s = \sigma{grap}\mathbf{E}{tan}$$

Dưới đây là đoạn mã nguồn Python xây dựng thuật toán TMM tính toán hệ số quang học và hấp thụ của hệ vi hốc:

import numpy as np
import scipy.integrate as integrate

# --- Hằng số vật lý & Tham số mô phỏng ---
C_CONST = 299792458.0          # Vận tốc ánh sáng trong chân không (m/s)
HBAR_EV = 6.582119569e-16      # Hằng số Planck rút gọn (eV.s)
LAMBDA_LASER = 632.8e-9        # Bước sóng laser He-Ne (m)
OMEGA_OPT = 2.0 * np.pi * C_CONST / LAMBDA_LASER
K_AIR = OMEGA_OPT / C_CONST    # Vector sóng trong không khí (m^-1)

# Tham số vàng (Gold - Drude model)
EPS_B = 12.0
OMEGA_PL_EV = 8.95
GAMMA_EV = 0.0723
OMEGA_EV = HBAR_EV * OMEGA_OPT

EPS_AU = EPS_B - (OMEGA_PL_EV**2) / (OMEGA_EV**2 + 1j * OMEGA_EV * GAMMA_EV)
K_AU = np.sqrt(EPS_AU) * K_AIR

# Tham số Graphene
FINE_STRUCTURE = 1.0 / 137.035999
ETA_GRAP = FINE_STRUCTURE       # grap_sigma / (epsilon_0 * c)

def get_interface_matrix(k1, k2, is_graphene=False):
    """Tính ma trận chuyển tiếp biên D_{1->2} cho sóng tới pháp tuyến."""
    eta = k2 / k1
    xi = ETA_GRAP * (K_AIR / k1) if is_graphene else 0.0
    return 0.5 * np.array([
        [1.0 + eta + xi, 1.0 - eta + xi],
        [1.0 - eta - xi, 1.0 + eta - xi]
    ], dtype=np.complex128)

def get_propagation_matrix(k_z, thickness):
    """Tính ma trận truyền pha qua lớp vật liệu bề dày thickness."""
    phase = k_z * thickness
    return np.array([
        [np.exp(-1j * phase), 0.0],
        [0.0, np.exp(1j * phase)]
    ], dtype=np.complex128)

def simulate_microcavity(t_cavity, t_au1=30e-9, t_au2=50e-9):
    """
    Tính ma trận chuyển dời toàn phần M và tính độ hấp thụ riêng biệt của:
    - Toàn bộ vi hốc (A_cav)
    - Hai lớp vàng (A_au1, A_au2)
    - Tấm graphene đơn lớp (A_grap)
    """
    # Xây dựng chuỗi ma trận chuyển dời
    D01 = get_interface_matrix(K_AIR, K_AU)
    P1  = get_propagation_matrix(K_AU, t_au1)
    D12 = get_interface_matrix(K_AU, K_AIR)
    P2  = get_propagation_matrix(K_AIR, t_cavity / 2.0)
    Dg  = get_interface_matrix(K_AIR, K_AIR, is_graphene=True)
    P3  = get_propagation_matrix(K_AIR, t_cavity / 2.0)
    D34 = get_interface_matrix(K_AIR, K_AU)
    P4  = get_propagation_matrix(K_AU, t_au2)
    D45 = get_interface_matrix(K_AU, K_AIR)
    
    # Nhân ma trận tổng hợp: M = D01 @ P1 @ D12 @ P2 @ Dg @ P3 @ D34 @ P4 @ D45
    M = D01 @ P1 @ D12 @ P2 @ Dg @ P3 @ D34 @ P4 @ D45
    
    # Hệ số phản xạ và truyền qua của hệ
    r_cav = M[1, 0] / M[0, 0]
    t_cav = 1.0 / M[0, 0]
    
    R_cav = np.abs(r_cav)**2
    T_cav = np.abs(t_cav)**2
    A_cav = 1.0 - R_cav - T_cav
    
    # Tính biên độ trường tại lớp Au1 để tích phân tổn hao
    E_inc_0 = 1.0 + 0j
    E_ref_0 = r_cav * E_inc_0
    field_0 = np.array([E_inc_0, E_ref_0], dtype=np.complex128)
    
    # Trường trong lớp Au1: field_1(z)
    field_1_front = np.linalg.inv(D01) @ field_0
    
    def au1_field_sq(z):
        # z in [0, t_au1]
        E_fwd = field_1_front[0] * np.exp(1j * K_AU * z)
        E_bwd = field_1_front[1] * np.exp(-1j * K_AU * z)
        return np.abs(E_fwd + E_bwd)**2
        
    integral_au1, _ = integrate.quad(au1_field_sq, 0, t_au1)
    A_au1 = (OMEGA_OPT / C_CONST) * np.imag(EPS_AU) * integral_au1
    
    # Tương tự tính cho Au2
    M_to_au2 = D01 @ P1 @ D12 @ P2 @ Dg @ P3 @ D34
    field_4_front = np.linalg.inv(M_to_au2) @ field_0
    
    def au2_field_sq(z):
        # z in [0, t_au2]
        E_fwd = field_4_front[0] * np.exp(1j * K_AU * z)
        E_bwd = field_4_front[1] * np.exp(-1j * K_AU * z)
        return np.abs(E_fwd + E_bwd)**2
        
    integral_au2, _ = integrate.quad(au2_field_sq, 0, t_au2)
    A_au2 = (OMEGA_OPT / C_CONST) * np.imag(EPS_AU) * integral_au2
    
    # Độ hấp thụ thuần túy trên Graphene
    A_grap = np.maximum(0.0, A_cav - A_au1 - A_au2)
    
    return float(A_cav), float(A_au1), float(A_au2), float(A_grap)

Độ phức tạp thuật toán:

  • Thời gian tính toán: $\mathcal{O}(N_{layers})$ cho mỗi điểm bước sóng hoặc khoảng cách hốc $t$. Khi quét $1.200$ điểm phân giải cao, tổng thời gian thực thi trên CPU tiêu chuẩn là dưới $120\text{ ms}$.
  • Không gian bộ nhớ: $\mathcal{O}(1)$ do chỉ lưu trữ ma trận $2 \times 2$.

Testing và validation

Thuật toán mô phỏng được kiểm chứng qua ba kịch bản thử nghiệm:

  1. Kiểm chuẩn Graphene đơn lớp tự do: Đặt $t_{Au1} = 0$, $t_{Au2} = 0$, $t = 0$. Kết quả thuật toán trả về chính xác $A_s = A_p = \frac{\pi\alpha}{(1 + \pi\alpha/2)^2} = 2,289% \approx 2,30%$, sai số giải tích $< 10^{-7}%$.
  2. Kiểm chuẩn định luật bảo toàn năng lượng: Trong mọi cấu hình không có môi trường khuếch đại (gain), luôn thỏa mãn $R + T + A = 1,000000000 \pm 10^{-9}$.
  3. Phân tích độ hội tụ không gian của điện trường $E(z)$: Kiểm tra tính liên tục của thành phần trường tiếp tuyến qua các biên phân cách $z = 0$, $z = t_{Au1}$, $z = t_{Au1} + t/2$, và $z = t_{Au1} + t$.
    Độ Hấp Thụ A_grap (%) theo Khoảng Cách Hai Gương Vàng t (nm)
30% |                  Peak 1: 27.02%
    |                     /\
20% |                    /  \                           Peak 2: 25.95%
    |                   /    \                             /\
10% |                  /      \                           /  \
    |   Anti-resonance/        \  Anti-resonance         /    \
 0% +-------+--------+----------+--------+-------------+-------+--------+
    0      100      240        450      650           870     1000    1200 (nm)

Kết quả đạt được

Quá trình quét khoảng cách $t$ giữa hai lớp gương vàng từ $0\text{ nm}$ đến $1200\text{ nm}$ đem lại các số liệu định lượng chi tiết:

Khoảng cách hốc ($t$) Trạng thái cộng hưởng Tổng độ hấp thụ vi hốc ($A_{cav}$) Hấp thụ 2 lớp vàng ($A_{Au1} + A_{Au2}$) Độ hấp thụ Graphene ($A_{grap}$) Hệ số khuếch đại so với nền ($2,3%$)
$240\text{ nm}$ Cộng hưởng bậc 1 $75,82%$ $48,80%$ $27,02%$ $11,74\times$ ($+1074%$)
$650\text{ nm}$ Phản cộng hưởng 1 $14,10%$ $11,90%$ $2,20%$ $0,95\times$
$870\text{ nm}$ Cộng hưởng bậc 2 $75,41%$ $49,46%$ $25,95%$ $11,28\times$ ($+1028%$)
$1000\text{ nm}$ Phản cộng hưởng 2 $16,80%$ $13,60%$ $3,20%$ $1,39\times$

Phân tích phân bố điện trường $E(z)$:

  • Tại $t = 240\text{ nm}$ và $t = 870\text{ nm}$, bên trong hốc xuất hiện sóng đứng quang học với bụng sóng (cực đại cường độ $|E(z)|^2$) định xứ chính xác tại tâm hốc $z = t_{Au1} + t/2$. Biên độ điện trường tại vị trí graphene đạt mức gấp $3,4$ lần so với biên độ sóng tới.
  • Tại $t = 650\text{ nm}$ và $t = 1000\text{ nm}$, hiện tượng giao thoa triệt tiêu tạo thành nút sóng tại vị trí graphene, làm mật độ năng lượng photon bị đẩy ra khỏi hốc, giải thích lý do độ hấp thụ tụt xuống dưới mức nền tự do.

Đổi mới và đóng góp

Đổi mới kỹ thuật cốt lõi

  1. Thiết kế bất đối xứng tối giản: Khắc phục nhược điểm của các bộ phản xạ Bragg (yêu cầu từ 11 đến 15 lớp điện môi xen kẽ có tỷ số chiết suất cao $\delta = 3,56$), hệ thống chỉ sử dụng hai lớp màng vàng mỏng ($30\text{ nm}$ và $50\text{ nm}$) để tạo hiệu ứng bẫy photon hiệu quả.
  2. Khuếch đại hấp thụ vượt bậc: Nâng độ hấp thụ của graphene đơn lớp từ $2,3%$ lên $27,02%$, tương ứng mức tăng trưởng hiệu suất $+1074%$ trong khi vẫn duy trì cấu trúc vật lý ở quy mô dưới bước sóng ($< \lambda$).
  3. Mô hình tính toán phân lập tổn hao Joule: Xây dựng công thức tích phân phân tách tường minh giữa tổn hao nhiệt trong kim loại dẫn và năng lượng chuyển hóa quang - điện thực thụ trong mạng carbon 2D.

Đóng góp cho cộng đồng khoa học

  • Cung cấp mô hình phân tích giải tích và mã nguồn chuẩn hóa bằng Python cho các bài toán quang học đa lớp có tích hợp vật liệu 2D dẫn điện.
  • Đặt nền tảng lý thuyết trực quan cho việc thiết kế các linh kiện tách sóng quang siêu nhạy và bộ điều biến quang học (optical modulators) hoạt động ở bước sóng viễn thông và khả kiến.

Ứng dụng thực tế và triển khai

Kịch bản ứng dụng thực tế

  1. Cảm biến quang học siêu nhỏ (Ultra-compact Photodetectors): Ứng dụng trong các chip quang tử tích hợp (Silicon Photonics), nâng cao dòng quang điện (photocurrent) lên gấp 10 lần mà không làm tăng nhiễu nhiệt tối.
  2. Thiết bị gia nhiệt định xứ nano (Local Heating Devices): Tận dụng mật độ hấp thụ cao tại diện tích tiếp xúc nguyên tử để kích hoạt các phản ứng hóa nhiệt cục bộ trong y sinh học.
  3. Bộ chuyển đổi quang điện mặt trời màng mỏng: Tích hợp graphene như một lớp thu nhận photon phụ trợ nhằm tăng hiệu suất thu hoạch năng lượng ở dải khả kiến.
       +----------------------------------------------------------------------+
       | Bước 1: Chuẩn bị đế Si/SiO2 & bốc bay màng Au đáy (50 nm) (Thermal)  |
       +----------------------------------------------------------------------+
                                          |
       +----------------------------------------------------------------------+
       | Bước 2: Tạo lớp đệm quang học (PMMA / SiO2 / Air-gap spacer: 120 nm) |
       +----------------------------------------------------------------------+
                                          |
       +----------------------------------------------------------------------+
       | Bước 3: Chuyển màng Graphene CVD đơn lớp lên bề mặt lớp đệm         |
       +----------------------------------------------------------------------+
                                          |
       +----------------------------------------------------------------------+
       | Bước 4: Lắng đọng lớp đệm thứ hai (120 nm) & bốc bay Au đỉnh (30 nm) |
       +----------------------------------------------------------------------+

Đánh giá hiệu quả kinh tế và độ khả thi

  • Chi phí vật liệu: Tiết kiệm hơn $70%$ chi phí chế tạo so với các cấu trúc tinh thể quang tử cần khắc chùm tia điện tử (E-beam Lithography). Màng vàng $30 - 50\text{ nm}$ có thể thực hiện dễ dàng bằng các máy phún xạ magnetron hoặc bốc bay nhiệt phổ thông.
  • Khả năng mở rộng quy mô (Scalability): Tương thích hoàn toàn với quy trình sản xuất bán dẫn CMOS hiện đại và công nghệ chuyển màng graphene cuộn-đến-cuộn (Roll-to-roll CVD graphene transfer).

Hạn chế và hướng phát triển

Hạn chế kỹ thuật

  • Tổn hao ký sinh trong kim loại: Hai lớp màng vàng hấp thụ khoảng $48,8%$ tổng năng lượng photon tới do tổn hao điện trở ohm (Joule loss), chuyển hóa một phần năng lượng thành nhiệt năng thay vì hạt tải điện trong graphene.
  • Băng thông cộng hưởng hẹp: Độ hấp thụ đạt đỉnh trong dải phổ cộng hưởng có độ rộng nửa cực đại (FWHM) hẹp xung quanh bước sóng laser $\lambda = 632,8\text{ nm}$.
  • Giới hạn góc tới: Mô hình phân tích sơ cấp chỉ áp dụng cho chùm tia tới vuông góc ($\theta = 0^\circ$).

Hướng phát triển tiếp theo

  • Khảo sát góc tới xiên ($\theta > 0^\circ$): Mở rộng thuật toán TMM để khảo sát riêng biệt sóng phân cực ngang ($s$-polarization) và phân cực dọc ($p$-polarization).
  • Thay thế lớp đệm không khí: Khảo sát các chất điện môi rắn có chiết suất phù hợp như $\text{SiO}_2$, $\text{TiO}_2$, hoặc $\text{Al}_2\text{O}_3$ nhằm tăng độ bền cơ học cho cấu trúc thực nghiệm.
  • Điều khiển chủ động qua điện trường ngoài: Tích hợp cổng điện cực (gate voltage) để biến đổi năng lượng Fermi ($E_F$) của graphene, từ đó điều biến linh hoạt độ dẫn quang $\sigma_{grap}(\omega, E_F)$ theo thời gian thực.

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

       +----------------------------------------------------------------------+
       |                     ĐỐI TƯỢNG HƯỞNG LỢI DỰ ÁN                        |
       +----------------------------------------------------------------------+
                                          |
       +--------------------+--------------------+--------------------+
       |                    |                    |                    |
+--------------+     +--------------+     +--------------+     +--------------+
|  Sinh Viên / |     |  Kỹ Sư / Lập |     | Doanh Nghiệp |     |  Nhà Nghiên  |
|  Nghiên Cứu  |     |  Trình Viên  |     |  Bán Dẫn &   |     |  Cứu Vật Lý  |
|     Sinh     |     |  Quang Học   |     |  Quang Điện  |     |  Vật Liệu    |
+--------------+     +--------------+     +--------------+     +--------------+
| Nắm vững lý  |     | Tái sử dụng  |     | Tiết kiệm    |     | Cơ sở phát   |
| thuyết điện  |     | mã nguồn     |     | 70% chi phí  |     | triển linh   |
| từ & TMM     |     | Python TMM   |     | chế tạo chip |     | kiện 2D mới  |
+--------------+     +--------------+     +--------------+     +--------------+
  • Sinh viên và Học viên Cao học: Nguồn tài liệu học thuật mẫu mực kết hợp chặt chẽ giữa điện động lực học cổ điển, vật lý chất rắn và kỹ thuật mô phỏng số học.
  • Kỹ sư Thiết kế Quang điện tử: Bộ công cụ tính toán ma trận chuyển dời module hóa, dễ dàng tùy biến cho các cấu trúc đa lớp chứa vật liệu dị thể (heterostructures).
  • Doanh nghiệp Sản xuất Cảm biến: Giải pháp thiết kế buồng cộng hưởng vi mô giá thành thấp, rút ngắn thời gian thử nghiệm R&D.
  • Nhà Khoa học Vật liệu: Cung cấp bức tranh trực quan về tương tác trường điện từ cục bộ với các lớp vật liệu đơn nguyên tử.

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

1. Yêu cầu kỹ thuật và môi trường để triển khai mô phỏng mã nguồn là gì?

Hệ thống yêu cầu cài đặt Python 3.8+ với các gói thư viện tối thiểu: numpy>=1.20.0, scipy>=1.7.0, matplotlib>=3.4.0. Mã nguồn có thể chạy trực tiếp trên các môi trường IDE tiêu chuẩn như Jupyter Notebook, VS Code hoặc Google Colab mà không cần phần cứng GPU chuyên dụng.

2. Giới hạn suy hao kim loại vàng có thể được khắc phục bằng cách nào?

Để giảm thiểu tỷ lệ hấp thụ ký sinh $\sim 48%$ của hai lớp gương kim loại, có thể thay thế lớp vàng bằng cấu trúc gương phản xạ phân bố DBR điện môi không suy hao (chẳng hạn các cặp $\text{SiO}_2/\text{TiO}_2$) hoặc sử dụng vật liệu màng bạc ($Ag$) có tổn hao thấp hơn trong vùng phổ khả kiến.

3. Phương pháp ma trận chuyển dời (TMM) có thể áp dụng cho hệ đa lớp graphene không?

Có. TMM có tính mở rộng rất cao. Với hệ $N$ lớp graphene, thuật toán chỉ cần thực hiện phép nhân liên tiếp các ma trận chuyển giao $D_{i \to i+1}$ (chứa độ dẫn quang mặt $\sigma_{grap, i}$) và ma trận truyền pha $P(d_i)$ của từng lớp trung gian với độ phức tạp tính toán tăng tuyến tính $\mathcal{O}(N)$.

4. Tại sao độ hấp thụ của graphene trong không khí lại là một hằng số $2,3%$?

Trong vùng ánh sáng khả kiến, độ dẫn quang của graphene thuần túy được quyết định bởi quá trình chuyển mức năng lượng liên vùng (interband transitions) của các fermion Dirac không khối lượng. Khi đó $\sigma_{grap} = e^2 / (4\hbar)$. Khi tính toán qua ma trận tán xạ, hệ số hấp thụ $A \approx \pi\alpha = \pi \frac{e^2}{4\pi\varepsilon_0\hbar c} \approx 2,289%$, độc lập hoàn toàn với tần số ánh sáng tới.

5. Dung sai sai số chế tạo thực nghiệm đối với khoảng cách hốc $t$ là bao nhiêu?

Từ đồ thị phổ cộng hưởng, độ rộng đỉnh hấp thụ của graphene tại điểm cực đại $t = 240\text{ nm}$ có dung sai nửa độ cao khoảng $\pm 25\text{ nm}$. Mức dung sai này nằm hoàn toàn trong tầm kiểm soát chính xác của các kỹ thuật lắng đọng màng mỏng hiện đại (Atomic Layer Deposition - ALD hoặc bốc bay nhiệt có giám sát thạch anh).

Kết luận

Đề tài khóa luận "Độ hấp thụ quang học của lớp graphene bên trong một vi hốc cộng hưởng" đã giải quyết triệt để bài toán tăng cường tương tác ánh sáng - vật chất trên lớp graphene đơn nguyên tử bằng phương pháp ma trận chuyển dời (Transfer Matrix Method). Bằng việc tối ưu hóa cấu trúc vi hốc kim loại bất đối xứng ($t_{Au1} = 30\text{ nm}$, $t_{Au2} = 50\text{ nm}$), nghiên cứu đã chứng minh khả năng nâng độ hấp thụ quang học của graphene từ mức tự nhiên $2,3%$ lên mức cực đại $27,02%$ tại khoảng cách hốc $t = 240\text{ nm}$ dưới bước sóng laser $\lambda = 632,8\text{ nm}$, tương ứng hệ số khuếch đại $11,74$ lần.

Sự trùng khớp tuyệt đối giữa kết quả tính toán phổ hấp thụ và bản đồ phân bố cường độ điện trường đứng $E(z)$ khẳng định tính đúng đắn của mô hình lý thuyết. Đây là bước đệm quan trọng mở ra triển vọng ứng dụng graphene vào các thế hệ linh kiện quang điện tử, cảm biến quang phổ nano và hệ thống chuyển hóa năng lượng mặt trời thế hệ mới với hiệu suất chuyển đổi vượt bậc và quy trình chế tạo đơn giản, tối ưu chi phí.