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

Bối cảnh thực tiễn và lý do chọn đề tài

Vật liệu graphene được hai nhà vật lý Andre Geim và Konstantin Novoselov phát hiện vào năm 2004 tại Đại học Manchester (giải Nobel Vật lý năm 2010). Đây là tinh thể hai chiều (2D) gồm các nguyên tử carbon liên kết dạng mạng lục giác với khoảng cách giữa các nguyên tử carbon $C-C$ là $0{,}142\text{ nm}$. Graphene sở hữu các đặc tính vật lý vượt trội như độ cứng gấp 200 lần thép, diện tích bề mặt lý thuyết đạt $2630\text{ m}^2/\text{g}$, độ dẫn nhiệt $5000\text{ W}/(\text{m}\cdot\text{K})$, độ linh động điện tử ở nhiệt độ phòng đạt $250000\text{ cm}^2/(\text{V}\cdot\text{s})$ tại mật độ hạt tải $2\times 10^{11}\text{ cm}^{-2}$, và độ truyền qua quang học đạt xấp xỉ $98%$.

Mặc dù có tiềm năng ứng dụng lớn trong quang điện tử, pin mặt trời và cảm biến, độ hấp thụ quang học của graphene đơn lớp trong miền ánh sáng nhìn thấy rất thấp, chỉ đạt khoảng $A = \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%$ này chưa đáp ứng được yêu cầu hiệu suất của các linh kiện thực tế như cảm biến quang hay thiết bị đốt nhiệt định xứ (local heating devices).

Các giải pháp tăng cường hấp thụ trước đây bộc lộ nhiều trở ngại:

  • Trong miền tần số terahertz (THz), việc tăng hấp thụ lên gần $100%$ đòi hỏi cấu hình phản xạ toàn phần và điều khiển mức năng lượng Fermi.
  • Ở tần số cao, các cấu trúc tinh thể quang tử dạng lỗ hoặc cấu trúc đa lớp điện môi xen kẽ đối xứng (dạng $ABABABGBABABA$) có thể nâng độ hấp thụ lên $50%$, nhưng đòi hỏi tỉ số chiết suất giữa hai lớp điện môi rất lớn ($\eta = n_A / n_B = 3{,}56$ với $s = 5$ cặp lớp) hoặc phải chế tạo chính xác độ dày từng lớp bằng đúng một phần tư bước sóng ánh sáng tới, gây khó khăn lớn cho quá trình thực nghiệm.

Nhằm giải quyết bài toán trên bằng một cấu trúc đơn giản hơn, tác giả lựa chọn đề tài nghiên cứu tính toán lý thuyết khả năng tăng cường độ hấp thụ quang học của lớp graphene đơn đặt tại tâm một vi hốc cộng hưởng (optical microcavity) tạo bởi hai lớp màng kim loại vàng mỏng.

Mục tiêu và nhiệm vụ nghiên cứu

Mục tiêu tổng quát của khóa luận là khảo sát và đánh giá sự tăng cường độ hấp thụ quang học của lớp graphene đơn lớp đặt bên trong vi hốc cộng hưởng quang học bằng phương pháp giải tích - số trị.

Nhiệm vụ nghiên cứu cụ thể bao gồm:

  1. Tổng quan cơ sở lý thuyết về quang học sóng, hệ phương trình Maxwell, tenxơ ứng suất Maxwell (Maxwell's stress tensor) và lực bức xạ ánh sáng.
  2. Xây dựng thuật toán ma trận chuyển dời (Transfer Matrix Method - TMM) cho sóng phân cực phẳng ($s$ và $p$) truyền qua hệ phân lớp chứa lớp dẫn điện bề mặt 2D (graphene).
  3. Thiết lập mô hình tính toán quang học cho vi hốc cộng hưởng tạo bởi hai lớp màng vàng có bề dày bất đối xứng ($t_{Au1} = 30\text{ nm}$ và $t_{Au2} = 50\text{ nm}$), ở giữa là môi trường không khí chứa một lớp graphene đơn tại tâm.
  4. Khảo sát độ hấp thụ quang học tổng thể của vi hốc, độ tiêu tán năng lượng trong từng lớp vàng và độ hấp thụ riêng của lớp graphene theo khoảng cách giữa hai gương vàng.
  5. Tính toán phân bố trường điện từ laser bên trong các miền của vi hốc cộng hưởng nhằm kiểm chứng cơ chế cộng hưởng quang học tăng cường hấp thụ.

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

  • Đối tượng nghiên cứu: Độ hấp thụ quang học, hệ số phản xạ, hệ số truyền qua và sự phân bố điện trường laser của lớp graphene đơn lớp và hệ vi hốc cộng hưởng kim loại.
  • Phạm vi không gian: Cấu trúc vi hốc phân lớp một chiều đặt trong môi trường không khí ($\varepsilon_{air} = 1$), gồm lớp vàng thứ nhất dày $t_{Au1} = 30\text{ nm}$, lớp không khí thứ nhất dày $t/2$, lớp graphene đơn tại tâm, lớp không khí thứ hai dày $t/2$, và lớp vàng thứ hai dày $t_{Au2} = 50\text{ nm}$. Khoảng cách $t$ giữa hai lớp vàng được khảo sát biến thiên từ $0$ đến $1000\text{ nm}$.
  • Phạm vi nguồn sáng: Chùm bức xạ laser đơn sắc có bước sóng khả kiến $\lambda = 632{,}8\text{ nm}$ chiếu theo phương vuông góc với bề mặt cấu trúc.

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

Khung lý thuyết và mô hình vật lý

Nghiên cứu sử dụng hệ thống lý thuyết điện từ học cổ điển kết hợp mô hình độ dẫn mặt lượng tử của graphene:

  • Hệ phương trình Maxwell và phương trình sóng d'Alembert: Mô tả sự lan truyền của sóng điện từ đơn sắc phẳng trong môi trường đồng nhất và đẳng hướng.
  • Tenxơ ứng suất Maxwell: Thiết lập biểu thức lực bức xạ và áp suất bức xạ tác dụng lên bản mỏng: $P_z = \varepsilon_0 (E_i^2 + E_r^2 - E_t^2)$.
  • Mô hình độ dẫn quang học bề mặt của graphene: Graphene đơn lớp được xem như một mặt phẳng dẫn điện 2D vô cùng mỏng tại $z = 0$, có độ dẫn điện mặt $\sigma_{grap} = \frac{e^2}{4\hbar}$.
  • Mô hình Drude - Lorentz cho màng kim loại vàng: Hàm điện môi phức của vàng $\varepsilon_{Au}(\omega)$ tại tần số quang học được xác định theo công thức: $$\varepsilon_{Au}(\omega_{opt}) = \varepsilon_b - \frac{\omega_{pl}^2}{\omega_{opt}^2 + i\gamma\omega_{opt}}$$ Trong đó, các tham số vật lý của vàng được sử dụng gồm: hằng số điện môi nền $\varepsilon_b = 12$, tần số plasmon $\hbar\omega_{pl} = 8{,}95\text{ eV}$, và hằng số tắt dần $\hbar\gamma = 72{,}3\text{ meV}$.

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

Khóa luận sử dụng phương pháp nghiên cứu lý thuyết kết hợp mô phỏng số:

  1. Phương pháp ma trận chuyển dời (Transfer Matrix Method - TMM): Dựa trên mô hình ma trận của Zhan và các cộng sự (2013) dành cho hệ đa lớp chứa graphene. Phương pháp liên hệ vector biên độ trường $(E_i, E_r)^T$ giữa các miền thông qua tích của ma trận chuyển tiếp qua mặt phân cách (Transmission matrix $D$) và ma trận truyền qua môi trường đồng nhất (Propagation matrix $P$).
  2. Phương pháp ma trận tán xạ (Scattering Matrix - S): Biến đổi từ ma trận chuyển dời $M$ để biểu diễn liên hệ trực tiếp giữa sóng tới và sóng phản xạ/truyền qua, hỗ trợ xác định hệ số phản xạ biên độ $r = M_{21}/M_{11}$ và hệ số truyền qua biên độ $t = 1/M_{11}$.
  3. Phương pháp tích phân tiêu tán Joule: Độ hấp thụ quang học của các lớp kim loại vàng được tính bằng tích phân phần ảo của hàm điện môi nhân với bình phương môđun điện trường dọc theo bề dày lớp kim loại: $$A_{Au} = \frac{\omega_{opt}}{c} \text{Im}(\varepsilon_{Au} - 1) \int_{0}^{t_{Au}} |E(z)|^2 dz$$ Từ đó, độ hấp thụ quang học của riêng lớp graphene được tách biệt chính xác: $$A_{grap} = A_{cav} - A_{Au1} - A_{Au2}$$
  4. Công cụ tính toán và xử lý số liệu: Thuật toán ma trận và tích phân giải tích được xây dựng và tính toán số trị bằng ngôn ngữ lập trình Python kết hợp phần mềm Wolfram Mathematica.

Nội dung chính theo từng phần của khóa luận

Phần 1: Cơ sở lý thuyết

Nội dung phần này hệ thống hóa các nền tảng vật lý phục vụ việc thiết lập mô hình tính toán:

  • Giới thiệu cấu trúc hình học mạng lục giác phẳng của graphene và các dạng thù hình carbon liên quan (fullerene 0D, ống nano carbon 1D, graphite 3D). Nêu bật các tính chất cơ, nhiệt, điện và quang học chuẩn của graphene.
  • Trình bày hệ phương trình Maxwell ở dạng vi phân và tích phân trong chân không và môi trường điện môi.
  • Khảo sát sự truyền sóng phẳng đơn sắc, thiết lập nghiệm sóng d'Alembert và chứng minh mối quan hệ trực giao giữa vector cường độ điện trường $\vec{E}$, vector cường độ từ trường $\vec{H}$ và vector sóng $\vec{k}$.
  • Thiết lập biểu thức tenxơ ứng suất Maxwell $T_{ij}$, áp dụng định luật bảo toàn động lượng của trường điện từ để tính toán mật độ lực Lorentz và áp suất bức xạ quang học tác dụng lên các lớp phẳng.

Phần 2: Sự hấp thụ quang học của lớp graphene trong không khí

Phần 2 đi sâu vào việc phát triển toán học cho phương pháp ma trận chuyển dời đối với graphene đơn lớp:

  • Khảo sát điều kiện biên điện từ học đối với sóng phân cực vuông góc (phân cực $s$) và phân cực song song (phân cực $p$) khi đi qua mặt phân cách chứa mật độ dòng điện mặt graphene $\vec{J} = \sigma_{grap} \vec{E}$.
  • Xây dựng dạng tổng quát của ma trận chuyển tiếp $D_{1\to 2}$ qua lớp graphene: $$D_{1\to 2} = \frac{1}{2} \begin{pmatrix} 1 + \eta_m + \xi_m & 1 - \eta_m - \xi_m \ 1 - \eta_m + \xi_m & 1 + \eta_m - \xi_m \end{pmatrix}$$ Trong đó các hệ số $\eta_m, \xi_m$ phụ thuộc vào góc tới, chiết suất môi trường và phân cực $m \in {s, p}$.
  • Thiết lập ma trận truyền pha $P(\Delta z)$ qua lớp điện môi có bề dày $\Delta z$.
  • Áp dụng tính toán độ phản xạ $R$, độ truyền qua $T$ và độ hấp thụ $A = 1 - R - T$ của một lớp graphene đơn tự do đặt trong không khí. Kết quả giải tích chứng minh rằng ở điều kiện chiếu vuông góc: $$R = \frac{\pi^2\alpha^2}{(2+\pi\alpha)^2}, \quad T = \frac{4}{(2+\pi\alpha)^2}, \quad A = \frac{\pi\alpha}{\left(1 + \frac{\pi\alpha}{2}\right)^2} \approx 2{,}3%$$ Kết quả này khẳng định độ hấp thụ của graphene đơn lớp trong không khí là một hằng số quang học không phụ thuộc vào tần số ánh sáng tới.

Phần 3: Sự hấp thụ quang của lớp graphene bên trong vi hốc cộng hưởng

Đây là nội dung trọng tâm của khóa luận, thực hiện tính toán số trị và phân tích mô hình vi hốc cộng hưởng:

  • Thiết lập ma trận chuyển dời toàn phần của hệ vi hốc gồm 5 lớp (gương vàng 1, khoang không khí 1, graphene, khoang không khí 2, gương vàng 2): $$M = D_{0\to 1} P(z_1) D_{1\to 2} P(z_2 - z_1) D_{2\to 3} P(z_3 - z_2) D_{3\to 4} P(z_4 - z_3) D_{4\to 5}$$
  • Khảo sát sự phụ thuộc của độ hấp thụ tổng $A_{cav}$ theo khoảng cách $t$ giữa hai gương kim loại trong dải từ $0$ đến $1000\text{ nm}$. Đường cong hấp thụ thể hiện hai đỉnh cực đại rõ rệt: đỉnh thứ nhất tại $t = 240\text{ nm}$ và đỉnh thứ hai tại $t = 870\text{ nm}$, đạt giá trị xấp xỉ $76%$.
  • Tính toán sự phân bố cường độ điện trường $|E(z)|^2$ bên trong các miền của vi hốc bằng các phương trình truyền sóng thành phần.
  • Thực hiện tích phân số xác định độ tiêu hao năng lượng trong hai lớp gương vàng $A_{Au1}$, $A_{Au2}$ và suy ra độ hấp thụ thuần của graphene $A_{grap}$.
  • So sánh phân bố điện trường laser tại các khoảng cách cộng hưởng ($t = 240\text{ nm}$ và $t = 870\text{ nm}$) với các khoảng cách không cộng hưởng ($t = 650\text{ nm}$ và $t = 1000\text{ nm}$), chứng minh hiện tượng khuếch đại điện trường tại tâm vi hốc chính là nguyên nhân trực tiếp dẫn tới sự tăng cường độ hấp thụ quang.

Kết quả và đóng góp khoa học

Bảng tổng hợp các thông số mô hình và kết quả tính toán

Dưới đây là các thông số vật lý của mô hình và kết quả tính toán quang học được trích xuất từ văn bản khóa luận:

Thông số / Đại lượng Giá trị trong mô hình Đơn vị / Ghi chú
Bước sóng nguồn laser ($\lambda$) $632{,}8$ $\text{nm}$ (ánh sáng đỏ khả kiến)
Bề dày lớp vàng thứ nhất ($t_{Au1}$) $30$ $\text{nm}$
Bề dày lớp vàng thứ hai ($t_{Au2}$) $50$ $\text{nm}$
Hằng số điện môi nền của vàng ($\varepsilon_b$) $12$ Không thứ nguyên
Năng lượng plasmon của vàng ($\hbar\omega_{pl}$) $8{,}95$ $\text{eV}$
Hằng số tắt dần của vàng ($\hbar\gamma$) $72{,}3$ $\text{meV}$
Độ dẫn mặt graphene ($\sigma_{grap}$) $e^2 / (4\hbar)$ $\Omega^{-1}$
Độ hấp thụ graphene đơn lớp tự do ($A_0$) $\approx 2{,}3$ $%$
Độ hấp thụ cực đại của toàn vi hốc ($A_{cav}^{max}$) $\approx 76$ $%$ (tại $t = 240\text{ nm}$ hoặc $870\text{ nm}$)
Độ hấp thụ cực đại của graphene tại $t = 240\text{ nm}$ $\approx 27$ $%$
Độ hấp thụ cực đại của graphene tại $t = 870\text{ nm}$ $\approx 26$ $%$

Đóng góp cụ thể của tác giả

  1. Khẳng định tính khả thi của cấu trúc vi hốc kim loại bất đối xứng: Chứng minh về mặt lý thuyết rằng không cần sử dụng hệ đa lớp điện môi phức tạp hay chiết suất dị thường, chỉ cần đặt graphene đơn lớp trong vi hốc cộng hưởng tạo bởi hai gương vàng mỏng ($30\text{ nm}$ và $50\text{ nm}$) cũng có thể nâng độ hấp thụ quang của graphene từ $2{,}3%$ lên tới $27%$ tại bước sóng $632{,}8\text{ nm}$. Mức tăng này tương đương gấp hơn $11{,}7$ lần so với graphene trong môi trường tự do.
  2. Làm rõ cơ chế vật lý thông qua phân bố trường: Bằng việc giải và vẽ phân bố điện trường $E(z)$ trong toàn bộ các miền không gian, tác giả chỉ ra rằng tại các khoảng cách cộng hưởng ($t = 240\text{ nm}$ và $t = 870\text{ nm}$), điện trường tại vị trí đặt graphene ở tâm vi hốc được tăng cường vượt bậc so với biên độ sóng tới. Ngược lại, tại các vị trí không cộng hưởng ($t = 650\text{ nm}$ và $t = 1000\text{ nm}$), điện trường bên trong vi hốc bị triệt tiêu, giải thích thỏa đáng sự sụt giảm độ hấp thụ.
  3. Ứng dụng thành công phương pháp ma trận chuyển dời mở rộng: Xây dựng mã nguồn tính toán số bằng Python mô tả chính xác tương tác quang học của vật liệu hai chiều kết hợp vật liệu kim loại có tán sắc phức theo mô hình Drude - Lorentz.

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

Hạn chế của đề tài

  • Nghiên cứu mới dừng lại ở mức độ tính toán lý thuyết và mô phỏng số một chiều, chưa tiến hành chế tạo mẫu thực nghiệm và đo đạc quang phổ thực tế.
  • Khảo sát chỉ giới hạn ở trường hợp sóng phẳng đơn sắc chiếu vuông góc ($\theta = 0^\circ$) với một bước sóng cố định $\lambda = 632{,}8\text{ nm}$.
  • Môi trường đệm bên trong vi hốc được giả định là không khí ($\varepsilon_{air} = 1$), chưa xét đến ảnh hưởng của các lớp đế điện môi thực tế dùng để nâng đỡ màng graphene (như $\text{SiO}_2$, $\text{PMMA}$, hoặc $\text{h-BN}$).

Hướng nghiên cứu mở rộng

Tác giả đề xuất hướng phát triển tiếp theo của đề tài:

  • Mở rộng mô hình tính toán lý thuyết cho trường hợp chùm bức xạ laser chiếu xiên góc bất kỳ so với pháp tuyến mặt phẳng tới (xét riêng phân cực $s$ và phân cực $p$).
  • Khảo sát vi hốc cộng hưởng được lấp đầy bởi các loại vật liệu điện môi khác nhau thay vì môi trường không khí nhằm tìm kiếm cấu hình tối ưu hơn cho các thiết bị quang điện tử tích hợp.

Giá trị tham khảo học thuật

Khóa luận là tài liệu tham khảo có giá trị học thuật đối với:

  • Sinh viên và học viên ngành Sư phạm Vật lý, Vật lý Quang học, Khoa học Vật liệu: Cung cấp tài liệu tổng quan mạch lạc về quang học của vật liệu 2D, phương pháp giải hệ phương trình Maxwell có điều kiện biên dẫn mặt, và cách vận dụng lý thuyết tenxơ ứng suất Maxwell.
  • Người nghiên cứu mô phỏng quang học nano (Nanophotonics): Phần trình bày chi tiết về phương pháp ma trận chuyển dời (TMM) cho hệ phân lớp chứa graphene và kỹ thuật tách độ hấp thụ kim loại bằng tích phân tiêu tán Joule là tài liệu mẫu rõ ràng để lập trình mô phỏng các cấu trúc cộng hưởng quang học khác.

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

1. Tại sao độ hấp thụ quang học của graphene đơn lớp trong không khí lại là hằng số 2,3%?

Theo kết quả tính toán giải tích từ điều kiện biên điện từ và phương pháp ma trận chuyển dời, độ dẫn quang học bề mặt của graphene ở miền ánh sáng nhìn thấy có giá trị lượng tử hóa $\sigma_{grap} = e^2 / (4\hbar)$. Khi chiếu vuông góc trong không khí, độ hấp thụ được tính trực tiếp từ hằng số cấu trúc tinh tế $\alpha = e^2 / (4\pi\varepsilon_0 \hbar c) \approx 1/137$ theo công thức $A = \pi\alpha / (1 + \pi\alpha/2)^2 \approx \pi\alpha \approx 2{,}3%$. Giá trị này hoàn toàn không phụ thuộc vào tần số hay bước sóng ánh sáng tới.

2. Hai lớp màng vàng tạo vi hốc có cấu trúc và kích thước như thế nào?

Vi hốc cộng hưởng trong khóa luận được thiết kế bất đối xứng gồm: lớp vàng phía trước (nơi chùm tia laser chiếu vào) có bề dày $t_{Au1} = 30\text{ nm}$ và lớp vàng phía sau có bề dày $t_{Au2} = 50\text{ nm}$. Lớp graphene đơn được đặt chính giữa, cách mỗi lớp vàng một khoảng bằng $t/2$.

3. Những khoảng cách nào giữa hai lớp vàng giúp độ hấp thụ của graphene đạt cực đại?

Độ hấp thụ của lớp graphene đạt cực đại tại hai giá trị khoảng cách giữa hai lớp vàng:

  • Đạt khoảng $27%$ khi khoảng cách $t = 240\text{ nm}$.
  • Đạt khoảng $26%$ khi khoảng cách $t = 870\text{ nm}$.

4. Tại sao độ hấp thụ tổng của vi hốc là 76% nhưng độ hấp thụ riêng của graphene chỉ đạt 27%?

Độ hấp thụ tổng của vi hốc $A_{cav} = 1 - R_{cav} - T_{cav}$ bao gồm năng lượng bị hấp thụ bởi cả lớp graphene và hai lớp màng vàng kim loại ($A_{cav} = A_{grap} + A_{Au1} + A_{Au2}$). Do vàng là kim loại có phần ảo của hằng số điện môi khác không ($\text{Im}(\varepsilon_{Au}) \neq 0$), sóng điện từ truyền qua vàng sẽ bị tiêu hao năng lượng do hiệu ứng Joule. Phần hấp thụ còn lại sau khi trừ đi tổn hao trong hai lớp vàng mới chính là độ hấp thụ thực tế trong lớp graphene.

5. Kết quả phân bố điện trường chứng minh điều gì?

Tính toán phân bố điện trường $E(z)$ cho thấy tại $t = 240\text{ nm}$ và $t = 870\text{ nm}$, biên độ điện trường laser tại tâm vi hốc (nơi đặt graphene) được tăng cường mạnh mẽ so với điện trường tới. Ngược lại, tại $t = 650\text{ nm}$ và $t = 1000\text{ nm}$, sóng phản xạ trong hốc giao thoa triệt tiêu khiến điện trường tại tâm rất nhỏ. Điều này chứng minh sự gia tăng hấp thụ quang học của graphene là kết quả trực tiếp của hiện tượng cộng hưởng khuếch đại vi hốc.

Kết luận

Khóa luận tốt nghiệp của tác giả Lương Hoàng Sang đã nghiên cứu thành công bài toán tăng cường độ hấp thụ quang học của lớp graphene đơn đặt tại tâm vi hốc cộng hưởng kim loại vàng bất đối xứng bằng phương pháp ma trận chuyển dời. Kết quả tính toán lý thuyết với nguồn laser bước sóng $632{,}8\text{ nm}$ cho thấy độ hấp thụ của graphene tăng từ mức tự nhiên $2{,}3%$ lên đạt đỉnh $27%$ tại khoảng cách hai gương $t = 240\text{ nm}$ và $26%$ tại $t = 870\text{ nm}$. Sự phù hợp giữa phân bố trường điện từ và phổ hấp thụ đã làm sáng tỏ cơ chế cộng hưởng quang học trong việc nâng cao hiệu suất tương tác ánh sáng - vật chất trên vật liệu hai chiều. Công trình là tài liệu tham khảo hoàn chỉnh, chuẩn xác về mặt giải tích và phương pháp số cho các nghiên cứu tiếp theo về linh kiện quang điện tử nano dựa trên graphene.