Tổng quan luận án

Trong bối cảnh phát triển của các hệ thống xử lý thông tin, các linh kiện quang tử và quang - điện tử cấu trúc micro và nano được định hướng ứng dụng trong các mạch vi quang - điện tử tích hợp với tốc độ xử lý và khả năng lưu trữ không phụ thuộc vào giới hạn vật lý của định luật Moore. Các cấu trúc này hướng đến việc chế tạo thế hệ linh kiện có kích thước nhỏ, khối lượng nhẹ, hiệu suất hoạt động cao, tiêu hao ít năng lượng và giảm giá thành sản xuất. Để hiện thực hóa mục tiêu này, hai hướng tiếp cận chính được xác định bao gồm: (i) sử dụng cấu trúc nhân tạo mới cho các phần tử cấu thành linh kiện; (ii) ứng dụng các vật liệu chức năng với các đặc tính quang học đặc thù.

Luận án tiến sĩ ngành Khoa học vật liệu, chuyên ngành Vật liệu quang học, quang điện tử và quang tử với tiêu đề “Nghiên cứu, thiết kế cấu trúc tinh thể quang tử 1D và 2D ứng dụng cho linh kiện lưỡng trạng thái ổn định” tập trung nghiên cứu việc thiết kế các cấu trúc nhân tạo không có sẵn trong tự nhiên nhằm phục vụ truyền thông, xử lý tín hiệu và chuyển mạch quang học.

Khoảng trống nghiên cứu được xác định từ thực tiễn: tại Việt Nam, các công trình nghiên cứu về linh kiện quang tử và quang - điện tử cấu trúc micro và nano còn rời rạc, thiếu tính hệ thống và gặp rào cản về trang thiết bị chế tạo mẫu thực nghiệm. Luận án giải quyết vấn đề này thông qua việc xây dựng mô hình lý thuyết kết hợp phương pháp tính toán và mô phỏng số có độ chính xác cao để phân tích đặc tính của cấu trúc tinh thể quang tử một chiều (1D) và hai chiều (2D).

Mục tiêu nghiên cứu của luận án được xác định gồm các nội dung sau:

  1. Tổng quan cơ sở lý thuyết về vật liệu có vùng cấm quang (Photonic Band Gap - PBG), cấu trúc tinh thể quang tử (Photonic Crystals - PhCs) và các đặc trưng của linh kiện lưỡng trạng thái quang ổn định.
  2. Đề xuất mô hình cấu trúc linh kiện quang tử mới, tính toán lý thuyết các tham số đặc trưng và đối sánh với kết quả mô phỏng.
  3. Nghiên cứu tính toán và mô phỏng có hệ thống để xác định các tham số tối ưu của cấu trúc PhCs 1D và 2D ứng dụng cho linh kiện lưỡng trạng thái quang ổn định, hướng đến hệ số phẩm chất $Q$ cao, cường độ quang kích hoạt chuyển mạch thấp và thời gian chuyển mạch nhanh.
  4. Đề xuất và thiết kế các cấu trúc kết hợp nhằm nâng cao hiệu suất làm việc và đặc tính vận hành của linh kiện lưỡng trạng thái quang ổn định.

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

  • Đối tượng nghiên cứu: Các phương pháp tính toán và mô phỏng số (FDTD, PWE, CMT); hiện tượng dẫn sóng và cộng hưởng quang học trong cấu trúc PhCs 1D và 2D; linh kiện lưỡng trạng thái quang ổn định ứng dụng trong chuyển mạch và xử lý tín hiệu toàn quang.
  • Phạm vi nghiên cứu: Khảo sát các cấu trúc tuần hoàn 1D (cách tử dẫn sóng, gương phản xạ Bragg, buồng vi cộng hưởng Fabry-Perot) và 2D (mạng tinh thể vuông, lục giác, buồng cộng hưởng khuyết tật, kênh dẫn sóng có khe hẹp) kết hợp môi trường vật liệu quang phi tuyến bậc ba.

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

Ngành khoa học quang tử (Photonics) xuất hiện từ những năm 80 của thế kỷ XIX và ghi nhận bước chuyển biến quan trọng vào năm 1987 khi E. Yablonovitch và S. John công bố khái niệm cấu trúc tinh thể quang tử (PhCs). Tính chất tuần hoàn của hằng số điện môi trong không gian tạo ra vùng cấm quang (PBG), ngăn cản sự truyền qua của sóng điện từ có tần số nằm trong dải cấm và cho phép kiểm soát, giam giữ ánh sáng mà không làm suy hao năng lượng.

Trên thế giới, các nhóm nghiên cứu tiêu biểu đã đóng góp các nền tảng lý thuyết và công cụ tính toán:

  • Nhóm nghiên cứu của Giáo sư J. Joannopoulos tại Viện Công nghệ Massachusetts (MIT, Hoa Kỳ) phát triển các phần mềm mã nguồn mở chuyên dụng như MIT Photonic Bands (MPB) và MIT Electromagnetic Equation Propagation (MEEP), hỗ trợ phân tích cấu trúc vùng năng lượng và mô phỏng truyền sóng điện từ.
  • Meade và cộng sự (1994) đề xuất mô hình buồng cộng hưởng và kênh dẫn sóng trên nền PhCs 2D thông qua việc tạo khuyết tật điểm và khuyết tật hàng (loại bỏ các hình trụ không khí).
  • Almeida đề xuất ý tưởng kênh dẫn sóng có khe hẹp; sau đó Robinson và cộng sự chứng minh cấu trúc buồng cộng hưởng có khe hẹp không khí (độ rộng khoảng 20 nm) trong PhCs 1D có khả năng tăng cường mật độ năng lượng ánh sáng định xứ.
  • Yamamoto và đồng nghiệp thiết kế buồng cộng hưởng khe hẹp trên PhCs 2D đạt hệ số phẩm chất $Q \approx 2 \times 10^5$. Di Falco và cộng sự công bố cấu trúc buồng cộng hưởng khe hẹp dị thường đạt $Q = 50.000$ tại bước sóng $\lambda = 1559,2\text{ nm}$, đồng thời phát triển cấu trúc kênh dẫn sóng có khe hẹp.
  • Về sự tiến triển của hệ số phẩm chất buồng cộng hưởng PhCs 2D: các công trình ban đầu ghi nhận giá trị $Q$ ở mức vài trăm; đến năm 2015, Ashfaqul Anwar Siraji và cộng sự đạt $Q = 4.459$; Z. Bushell và cộng sự (2017) đạt $Q > 10^4$; đến năm 2018, các cấu trúc trên nền vật liệu silic đạt $Q$ trong khoảng $10^9 - 10^{11}$. Các cấu trúc hình học chuẩn bao gồm: buồng cộng hưởng $H_0$ ($Q \approx 1,12 \times 10^5$ tại $\lambda = 1564\text{ nm}$, có thể tối ưu lên $1,7 \times 10^6$), buồng cộng hưởng $H_1$ ($Q$ đạt $2 \times 10^6$), buồng cộng hưởng $L_3$ (loại bỏ ba hình trụ không khí trung tâm), và buồng cộng hưởng cấu trúc dị thường ($Q$ đạt $2 \times 10^7$).

Tại Việt Nam, hướng nghiên cứu PhCs được triển khai tại Viện Khoa học vật liệu, Viện Vật lý (thuộc Viện Hàn lâm Khoa học và Công nghệ Việt Nam), Trường Đại học Bách khoa Hà Nội, Học viện Công nghệ Bưu chính Viễn thông, Trường Đại học Vinh, và Trường Đại học Khoa học - Đại học Huế. Trong đó, nhóm nghiên cứu của PGS. Phạm Văn Hội và PGS. Phạm Thu Nga tại Viện Khoa học vật liệu đã chế tạo cấu trúc PhCs 1D và 3D trên nền vật liệu silic xốp và silica ứng dụng cho cảm biến đo chất lỏng.

Vị trí của luận án được xác lập nhằm giải quyết tính chất rời rạc của các nghiên cứu trong nước, tập trung vào việc mô hình hóa lý thuyết, tính toán tham số và mô phỏng tối ưu các cấu trúc PhCs 1D và 2D cho chức năng lưỡng ổn định quang học, tạo cơ sở cho việc thiết kế linh kiện chuyển mạch quang tích hợp.

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

1. Cơ sở lý thuyết vật lý

Cấu trúc tinh thể quang tử (PhCs) và vùng cấm quang (PBG)

Cấu trúc PhCs là môi trường điện môi có chiết suất biến thiên tuần hoàn trong không gian theo một, hai hoặc ba chiều (1D, 2D, 3D). Các tham số đặc trưng gồm:

  • Hằng số mạng ($a$): Chu kỳ không gian của ô mạng cơ sở.
  • Hệ số lấp đầy ($f$): Tỷ lệ thể tích của các phần tử điện môi so với toàn bộ thể tích cấu trúc.
  • Chiết suất hiệu dụng ($n_{\text{eff}}$): Được tính thông qua hằng số điện môi hiệu dụng $\varepsilon_{\text{eff}}$: $$n_{\text{eff}} = \sqrt{\varepsilon_{\text{eff}}} = \sqrt{f \varepsilon_1 + (1 - f) \varepsilon_2}$$ Trong đó $\varepsilon_1$ và $\varepsilon_2$ lần lượt là hằng số điện môi của hai vật liệu cấu thành.

Đối với cấu trúc PhCs 1D tuần hoàn yếu ($\Delta\varepsilon/\varepsilon \ll 1$ hoặc $d/a$ nhỏ), tỷ số độ rộng vùng cấm quang so với tần số trung tâm ($\omega_m$) được xác định theo công thức giải tích: $$\frac{\Delta\omega}{\omega_m} \approx \frac{2\Delta\varepsilon}{\pi\varepsilon} \sin\left(\frac{\pi d}{a}\right)$$

Phân cực và mode dẫn sóng trong cấu trúc PhCs 2D

Sóng điện từ truyền trong mặt phẳng $xy$ của PhCs 2D phân tách thành hai trạng thái phân cực độc lập:

  • Mode điện trường ngang (TE): Từ trường hướng dọc theo trục trụ ($H = H(\vec{r})\hat{z}$), điện trường nằm trong mặt phẳng $xy$.
  • Mode từ trường ngang (TM): Điện trường hướng dọc theo trục trụ ($E = E(\vec{r})\hat{z}$), từ trường nằm trong mặt phẳng $xy$.

Khảo sát giản đồ năng lượng cho thấy:

  • Cấu trúc trụ điện môi ($\varepsilon = 12$, bán kính $r = 0,2a$) trong không khí mở ra vùng cấm quang hoàn toàn đối với phân cực TM giữa dải tần số thứ nhất và thứ hai với độ rộng tương đối $\Delta\omega/\omega_m = 47%$.
  • Cấu trúc trụ không khí ($r = 0,3a$) trong nền điện môi ($\varepsilon = 12$) mở ra vùng cấm quang hoàn toàn đối với phân cực TE với độ rộng tương đối $\Delta\omega/\omega_m = 28%$.

Cơ chế cộng hưởng và hệ số phẩm chất $Q$

Hệ số phẩm chất $Q$ của buồng cộng hưởng biểu thị khả năng giam giữ năng lượng và thời gian sống của photon bên trong cấu trúc: $$Q = \omega_0 \frac{U}{-\frac{dU}{dt}} = \frac{\omega_0}{\Delta\omega}$$ Trong đó $\omega_0$ là tần số cộng hưởng trung tâm, $U$ là năng lượng điện từ tích trữ, và $\Delta\omega$ là độ rộng phổ tại nửa độ cao cực đại (FWHM).

Nguyên lý lưỡng trạng thái quang ổn định (Optical Bistability)

Hệ quang học biểu hiện đặc tính lưỡng ổn định khi thỏa mãn hai điều kiện: môi trường có tính phi tuyến quang học và tồn tại cơ chế phản hồi ngược. Mối quan hệ giữa cường độ ánh sáng đầu vào ($I_{\text{vào}}$) và cường độ đầu ra ($I_{\text{ra}}$) tuân theo biểu thức: $$I_{\text{ra}} = f(I_{\text{ra}}) \cdot I_{\text{vào}}$$ Khi hàm truyền qua $f(I_{\text{ra}})$ là một hàm phi tuyến không đơn điệu (có dạng hình chuông), đường đặc tuyến $I_{\text{ra}} - I_{\text{vào}}$ xuất hiện vùng đa trị:

  • Tại dải cường độ thấp ($I_{\text{vào}} < I_{\text{vào}1}$) và dải cường độ cao ($I_{\text{vào}} > I_{\text{vào}2}$), hệ có một trạng thái truyền qua đơn nhất.
  • Trong vùng trung gian ($I_{\text{vào}1} < I_{\text{vào}} < I_{\text{vào}2}$), một giá trị $I_{\text{vào}}$ tương ứng với hai trạng thái đầu ra ổn định (trạng thái dẫn truyền cao và trạng thái dẫn truyền thấp) cùng một trạng thái trung gian không ổn định.
  • Quá trình chuyển trạng thái diễn ra thông qua hiệu ứng trễ khi tăng $I_{\text{vào}}$ vượt ngưỡng chuyển mạch trên hoặc giảm $I_{\text{vào}}$ xuống dưới ngưỡng chuyển mạch dưới.

Vật liệu phi tuyến bậc ba được tích hợp trong cấu trúc gồm chalcogenide ($As_2S_3$) và polymer hữu cơ (DDMEBT). Các vật liệu này có khả năng chịu được mật độ công suất laser xung ở mức $\sim\text{GW/cm}^2$ do đặc tính không dẫn nhiệt và độ bền nhiệt (ngưỡng bền nhiệt của $As_2S_3$ đạt khoảng $300^\circ\text{C}$).

2. Phương pháp tính toán và mô phỏng

Luận án sử dụng ba phương pháp toán học và mô phỏng số bổ trợ lẫn nhau:

  - Tương đương mạch LC          - Miền tần số                  - Miền thời gian (Yee)
  - Phân tích ghép kênh          - Toán tử Hermit               - Bước nhảy cóc thời gian
    dẫn sóng & buồng             - Giải bài toán trị riêng      - Điều kiện biên PML
    cộng hưởng                   - Giản đồ PBG & tán sắc        - Khảo sát đáp ứng phi
                                                                  tuyến, phổ truyền/phản xạ

Lý thuyết ghép cặp mode theo thời gian (Coupled Mode Theory in Time - CMT)

Phương pháp CMT mô tả sự tương tác giữa sóng truyền trong kênh dẫn và mode dao động trong buồng cộng hưởng dựa trên mô hình tương đương mạch dao động $LC$. Phương trình vi phân mô tả biên độ trường $a(t)$ của buồng cộng hưởng có tổn hao khi ghép cặp với sóng tới $s(t)$ có dạng: $$\frac{da}{dt} = \left(j\omega_0 - \frac{1}{\tau_0} - \frac{1}{\tau_e}\right)a + \kappa s$$ Trong đó:

  • $\omega_0 = 1/\sqrt{LC}$ là tần số cộng hưởng riêng.
  • $1/\tau_0$ là hệ số suy hao nội tại của buồng cộng hưởng.
  • $1/\tau_e$ là tốc độ thoát năng lượng vào kênh dẫn sóng.
  • $\kappa$ là hệ số ghép cặp giữa kênh dẫn và buồng cộng hưởng; theo nguyên lý đảo ngược thời gian và bảo toàn năng lượng khi không có tổn hao nội ($1/\tau_0 = 0$), hệ số ghép cặp được xác định bằng: $$\kappa = \sqrt{\frac{2}{\tau_e}}$$

Phương pháp khai triển sóng phẳng (Plane Wave Expansion - PWE)

Phương pháp PWE giải hệ phương trình Maxwell trong miền tần số đối với môi trường điện môi tuần hoàn không có dòng tự do ($\vec{J} = 0$) và không có điện tích tự do ($\rho = 0$). Phương trình sóng đối với từ trường được viết dưới dạng bài toán trị riêng với toán tử Hermit: $$\vec{\nabla} \times \left( \frac{1}{\varepsilon_r(\vec{r})} \vec{\nabla} \times \vec{H}(\vec{r}) \right) = \left(\frac{\omega}{c}\right)^2 \vec{H}(\vec{r})$$ Do tính chất tuần hoàn của mạng tinh thể, điện từ trường và hằng số điện môi được khai triển thành chuỗi Fourier theo các vector mạng đảo $\vec{G}$. Việc rút gọn chuỗi Fourier thành số lượng số hạng hữu hạn chuyển phương trình vi phân thành hệ phương trình đại số ma trận, cho phép xác định chính xác cấu trúc vùng cấm quang (PBG) và đường cong tán sắc của cấu trúc PhCs.

Phương pháp đạo hàm hữu hạn trong miền thời gian (Finite-Difference Time-Domain - FDTD)

Thuật toán FDTD do Kane Yee đề xuất (1966) sử dụng lưới sai phân không gian (mạng Yee) và sơ đồ sai phân trung tâm nhảy cóc theo thời gian. Các thành phần điện trường $\vec{E}$ và từ trường $\vec{H}$ được bố trí xen kẽ nhau trong không gian và tính toán lệch pha nhau nửa bước thời gian:

  • Thành phần điện trường $\vec{E}$ được tính tại các thời điểm nguyên: $t = n\Delta t$.
  • Thành phần từ trường $\vec{H}$ được tính tại các thời điểm bán nguyên: $t = (n + 1/2)\Delta t$.

Các tham số và điều kiện vận hành thuật toán FDTD:

  • Điều kiện ổn định Courant - Friedrichs - Lewy (CFL): Bước thời gian $\Delta t$ phải thỏa mãn bất đẳng thức: $$\Delta t \le \frac{\Delta x}{\sqrt{n} c_0}$$ Với $n$ là số chiều không gian ($n = 1, 2, 3$), $\Delta x$ là kích thước ô lưới không gian, và $c_0$ là vận tốc ánh sáng trong chân không.
  • Điều kiện biên hấp thụ hoàn hảo (Perfect Matched Layer - PML): Áp dụng kỹ thuật phối hợp trở kháng sóng của J. Berenger để triệt tiêu phản xạ tại các biên của không gian tính toán thông qua việc gán hệ số suy hao điện và từ tương thích: $$\eta = \sqrt{\frac{\mu}{\varepsilon}} = \sqrt{\frac{\mu_0 \mu_r^}{\varepsilon_0 \varepsilon_r^}} = \eta_0$$
  • Kỹ thuật giảm không gian tính toán bằng đối xứng: Tận dụng đối xứng gương hoặc đối xứng quay $180^\circ$ ($C_2$) để giảm dung lượng bộ nhớ máy tính và rút ngắn thời gian tính toán xuống 2 lần.

Hệ thống tính toán và kiểm chứng mô hình

Các thuật toán PWE và FDTD được thực thi thông qua hai bộ công cụ phần mềm mã nguồn mở là MPB và MEEP, vận hành trên hệ thống tính toán song song hiệu năng cao tại phòng thí nghiệm.

Độ tin cậy của mô hình mô phỏng được kiểm chứng thông qua việc tái tạo độc lập hai công trình quốc tế đã công bố:

  1. Mô phỏng cấu trúc bộ lọc quang học gồm các trụ điện môi ($\varepsilon = 11,56$, bán kính $r = 0,2a$) kết hợp các trụ biến tính ($\varepsilon = 9,5$) và buồng cộng hưởng trung tâm ($\varepsilon = 6,6$, $r = 0,05a$). Kết quả tính toán phổ truyền qua và phổ phản xạ bằng FDTD cho kết quả trùng khớp với dữ liệu công bố của nhóm tác giả H.
  2. Mô phỏng cấu trúc bộ lọc quang học chia 5 kênh đầu ra trên nền trụ điện môi ($\varepsilon = 11,56$). Dữ liệu phổ truyền qua tại các cổng ra thu được từ mô phỏng FDTD của luận án khớp với các kết quả thực nghiệm và tính toán của nhóm nghiên cứu tại Hàn Quốc.

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

Chương Tên chương theo văn bản Nội dung khoa học và trọng tâm kỹ thuật
Chương 1 Tổng quan Tổng quan cấu trúc PhCs 1D, 2D, 3D; cơ chế hình thành vùng PBG; phân cực TE/TM; cấu trúc buồng cộng hưởng ($H_0, H_1, L_3$, dị thường, khe hẹp); kênh dẫn sóng; nguyên lý chuyển mạch và lưỡng trạng thái quang ổn định.
Chương 2 Phương pháp tính toán và mô phỏng Cơ sở toán học của CMT, PWE, FDTD; cấu trúc ô mạng Yee; điều kiện biên PML; điều kiện ổn định CFL; kỹ thuật đối xứng $C_2$; kiểm chứng độ chính xác thuật toán qua các công bố quốc tế.
Chương 3 Tối ưu hóa hệ số phẩm chất và phổ cộng hưởng của cấu trúc cách tử dẫn sóng Khảo sát cấu trúc PhCs 1D cách tử dẫn sóng; tối ưu hóa các tham số hình học (chu kỳ cách tử, độ dày lớp điện môi, chiết suất) nhằm nâng cao hệ số phẩm chất $Q$ và thu hẹp bán độ rộng phổ cộng hưởng $\Delta\omega$.
Chương 4 Lưỡng trạng thái quang ổn định trong cấu trúc cách tử dẫn sóng Tích hợp vật liệu quang phi tuyến vào cấu trúc cách tử dẫn sóng 1D; thiết lập hàm truyền qua phi tuyến; xác định ngưỡng cường độ chuyển mạch quang và khảo sát đường trễ lưỡng ổn định.
Chương 5 Lưỡng trạng thái quang ổn định dựa trên sự tương tác giữa cộng hưởng và dẫn sóng khe hẹp trong cấu trúc tinh thể quang tử hai chiều Nghiên cứu cấu trúc PhCs 2D kết hợp buồng cộng hưởng và kênh dẫn sóng có khe hẹp; phân tích sự tương tác trường định xứ trong khe hẹp; tối ưu hóa cấu trúc nhằm giảm thiểu cường độ quang kích hoạt và nâng cao tốc độ chuyển mạch.

Chương 1. Cấu trúc tinh thể quang tử và linh kiện lưỡng trạng thái quang ổn định

Chương này tổng hợp các khái niệm cơ bản về vật liệu có cấu trúc nhân tạo, bao gồm PhCs, plasmonics và siêu vật liệu (metamaterials). Nội dung đi sâu phân loại ba dạng cấu trúc PhCs:

  • PhCs 1D: Phân tích màng đa lớp tuần hoàn, gương phản xạ Bragg (DBR) với số chu kỳ $N$, buồng vi cộng hưởng Fabry-Perot ngăn cách bởi lớp đệm có chiết suất $n_s$ và bề dày $d_s$, cấu trúc cách tử dẫn sóng hoạt động dựa trên điều kiện giao thoa phản xạ Bragg: $$a(\sin\theta_1 - \sin\theta_2) = m \lambda_n$$
  • PhCs 2D: Phân tích mạng không gian thực, mạng đảo và vùng Brillouin thứ nhất của mạng ô vuông và mạng lục giác; phân tách hai mode phân cực TE và TM; nguyên lý giam giữ và dẫn truyền ánh sáng qua việc tạo khuyết tật điểm (buồng cộng hưởng) và khuyết tật hàng (kênh dẫn sóng). Khảo sát các dạng buồng cộng hưởng có khe hẹp không khí kích thước nano (khoảng 20 nm) nhằm tăng cường tương tác ánh sáng - vật chất.
  • Linh kiện lưỡng trạng thái quang ổn định: Trình bày các cấu hình chuyển mạch quang ($1 \times 1, 1 \times 2, 2 \times 2, N \times N$), các tham số hoạt động (kích thước, thời gian chuyển mạch, lưu lượng, năng lượng chuyển mạch, công suất hao phí), và cơ chế phản hồi ngược phi tuyến bậc ba trong môi trường cộng hưởng.

Chương 2. Phương pháp tính toán và mô phỏng

Chương này trình bày chi tiết cơ sở toán học và thuật toán số của ba công cụ phân tích:

  • Thiết lập phương trình vi phân biến thiên biên độ và trao đổi năng lượng của lý thuyết ghép cặp mode theo thời gian (CMT), dẫn xuất hệ số ghép cặp giữa buồng cộng hưởng và kênh dẫn sóng.
  • Rời rạc hóa hệ phương trình vi phân Maxwell theo phương pháp khai triển sóng phẳng (PWE) để giải bài toán trị riêng xác định giản đồ vùng năng lượng PBG cho các phân cực TE và TM.
  • Triển khai thuật toán FDTD của Kane Yee trên lưới Đề-các, áp dụng điều kiện biên hấp thụ PML để hạn chế phản xạ giả tạo từ biên tính toán, xác lập bước thời gian theo tiêu chuẩn hội tụ CFL, áp dụng kỹ thuật giảm bậc đối xứng không gian ($C_2$). Hoàn thiện kiểm chứng mô hình mô phỏng thông qua việc tái tạo kết quả phổ của các linh kiện lọc quang học đã công bố trên các tạp chí quốc tế.

Chương 3. Tối ưu hóa hệ số phẩm chất và phổ cộng hưởng của cấu trúc cách tử dẫn sóng

Chương này tập trung vào việc mô hình hóa và tính toán tối ưu các tham số cấu trúc của cách tử dẫn sóng (cấu trúc PhCs 1D). Trọng tâm khảo sát gồm ảnh hưởng của chu kỳ cách tử, tỷ lệ độ dày các lớp điện môi và sự chênh lệch chiết suất đến đặc tính cộng hưởng Bragg. Các phép tính định lượng hướng tới việc thu hẹp độ rộng phổ FWHM ($\Delta\omega$) và nâng cao hệ số phẩm chất $Q$, tạo tiền đề cho việc giam giữ trường điện từ phục vụ tương tác phi tuyến.

Chương 4. Lưỡng trạng thái quang ổn định trong cấu trúc cách tử dẫn sóng

Nội dung chương 4 khảo sát sự xuất hiện của hiệu ứng lưỡng ổn định quang học khi đưa môi trường có độ cảm phi tuyến bậc ba vào cấu trúc cách tử dẫn sóng 1D. Dựa trên các phương trình CMT và mô phỏng FDTD, chương phân tích sự biến đổi của chiết suất phi tuyến dưới tác dụng của cường độ ánh sáng kích thích, từ đó xác lập các giá trị ngưỡng cường độ quang chuyển mạch ($I_{\text{vào}1}, I_{\text{vào}2}$) và khảo sát đường đặc tuyến trễ quang học của linh kiện.

Chương 5. Lưỡng trạng thái quang ổn định dựa trên sự tương tác giữa cộng hưởng và dẫn sóng khe hẹp trong cấu trúc tinh thể quang tử hai chiều

Chương 5 mở rộng bài toán sang cấu trúc PhCs 2D, tập trung vào cấu hình kết hợp giữa buồng cộng hưởng và kênh dẫn sóng có khe hẹp kích thước nano. Bằng cách tận dụng hiệu ứng ánh sáng chậm và sự tăng cường mật độ năng lượng điện trường bên trong khe hẹp có chiết suất thấp, cấu trúc cho phép tăng cường hiệu ứng quang phi tuyến. Kết quả tính toán và mô phỏng xác định các tham số hình học tối ưu nhằm giảm cường độ ngưỡng quang đầu vào cần thiết cho quá trình chuyển mạch, đồng thời đảm bảo hệ số phẩm chất $Q$ cao và thời gian đáp ứng nhanh ở thang pico giây (ps).

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

Đóng góp về mặt học thuật và lý thuyết

  • Xây dựng hệ thống tính toán và mô phỏng hoàn chỉnh kết hợp giữa lý thuyết ghép cặp mode theo thời gian (CMT), phương pháp khai triển sóng phẳng (PWE) và phương pháp đạo hàm hữu hạn trong miền thời gian (FDTD) để khảo sát các cấu trúc PhCs 1D và 2D.
  • Cung cấp các phân tích định lượng về sự phụ thuộc của vùng cấm quang (PBG), hệ số phẩm chất $Q$ và phân bố trường điện từ vào các tham số hình học (hằng số mạng $a$, bán kính lỗ/trụ $r$, độ rộng khe hẹp, hệ số lấp đầy $f$) đối với cả hai trạng thái phân cực TE và TM.

Đóng góp về mặt thiết kế và ứng dụng thực tiễn

  • Đề xuất mô hình cấu trúc cách tử dẫn sóng 1D và cấu trúc PhCs 2D kết hợp kênh dẫn sóng - buồng cộng hưởng có khe hẹp nano đạt hệ số phẩm chất $Q$ cao và giảm thiểu cường độ quang kích hoạt chuyển mạch.
  • Cung cấp các bộ thông số thiết kế tối ưu phục vụ việc chế tạo các phần tử chức năng trong mạch tích hợp toàn quang như cổng logic quang, bộ nhớ quang học, bộ lọc dải hẹp và bộ chuyển mạch quang tốc độ cao.
  • Thiết lập hệ thống kiểm chứng tính toán mô phỏng trên nền tảng phần mềm mã nguồn mở MEEP và MPB chạy trên hệ thống máy tính hiệu năng cao, đóng vai trò tài liệu tham khảo kỹ thuật cho các nghiên cứu tiếp theo trong nước.

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

  • Hạn chế của luận án: Do điều kiện trang thiết bị công nghệ chế tạo mẫu vi/nano trong nước còn hạn chế, các kết quả của luận án chủ yếu tập trung vào xây dựng mô hình lý thuyết, tính toán giải tích và mô phỏng số, chưa thực hiện chế tạo thực nghiệm toàn diện cho các linh kiện PhCs 2D đề xuất. Ngoài ra, các yếu tố về băng thông và thời gian chuyển mạch (ở thang pico giây) được khảo sát dưới điều kiện kích thích nguồn laser đơn sắc lý tưởng.
  • Hướng nghiên cứu tiếp theo: Phát triển quy trình công nghệ chế tạo thực nghiệm các mẫu cách tử dẫn sóng 1D và PhCs 2D trên nền vật liệu silic hoặc màng mỏng phi tuyến; mở rộng nghiên cứu kết hợp giữa cấu trúc PhCs với cấu trúc cộng hưởng plasmon bề mặt và buồng vi cộng hưởng lai để tối ưu hóa hơn nữa kích thước và hiệu suất linh kiện.

Giá trị tham khảo

Luận án là tài liệu chuyên khảo có giá trị đối với các nhóm đối tượng:

  • Nghiên cứu sinh, học viên cao học và giảng viên: Chuyên ngành Khoa học vật liệu, Vật lý chất rắn, Quang học, Quang điện tử và Công nghệ viễn thông có thể tham khảo phương pháp luận thiết lập mô hình tính toán, các dẫn xuất giải tích phương trình Maxwell, và kỹ thuật mô phỏng số sử dụng MEEP và MPB.
  • Kỹ sư và nhà thiết kế linh kiện quang điện tử: Tham khảo các cấu hình hình học, tỷ lệ tham số mạng và phân tích ngưỡng năng lượng chuyển mạch để định hướng thiết kế các mạch vi quang tích hợp, cảm biến quang học dải hẹp và các phần tử xử lý tín hiệu toàn quang.

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

1. Hai điều kiện cơ bản để hình thành hiện tượng lưỡng trạng thái quang ổn định là gì?

Hiện tượng lưỡng trạng thái quang ổn định đòi hỏi hệ quang học phải có: (i) môi trường vật liệu có tính chất quang phi tuyến (chiết suất phụ thuộc vào cường độ ánh sáng kích thích, ví dụ phi tuyến bậc ba Kerr) và (ii) cơ chế phản hồi ngược quang học (thực hiện qua buồng cộng hưởng hoặc cấu trúc cách tử phản xạ).

2. Sự khác biệt về vùng cấm quang (PBG) giữa hai mode phân cực TE và TM trong PhCs 2D là gì?

Đối với cấu trúc PhCs 2D dạng các trụ điện môi ($\varepsilon = 12, r = 0,2a$) đặt trong không khí, vùng cấm quang hoàn toàn chỉ xuất hiện đối với mode phân cực TM (với độ rộng tương đối $\Delta\omega/\omega_m = 47%$). Ngược lại, đối với cấu trúc các trụ không khí ($r = 0,3a$) đặt trong nền điện môi ($\varepsilon = 12$), vùng cấm quang hoàn toàn xuất hiện đối với mode phân cực TE (với độ rộng tương đối $\Delta\omega/\omega_m = 28%$).

3. Phương pháp FDTD xử lý hiện tượng phản xạ sóng không mong muốn tại biên tính toán bằng cách nào?

FDTD áp dụng điều kiện biên hấp thụ hoàn hảo (PML) do J. Berenger đề xuất. Kỹ thuật này thiết lập sự phối hợp trở kháng sóng tại mặt phân cách giữa môi trường mô phỏng và lớp biên hấp thụ, giúp sóng điện từ lan truyền đến biên bị suy hao hoàn toàn mà không phản xạ ngược trở lại vùng tính toán.

4. Hệ số phẩm chất $Q$ của buồng cộng hưởng quang học được xác định như thế nào trong mô phỏng?

Hệ số phẩm chất $Q$ được định lượng theo tỷ số giữa tần số cộng hưởng trung tâm $\omega_0$ và bán độ rộng phổ cộng hưởng tại nửa độ cao cực đại $\Delta\omega$ (FWHM) theo biểu thức $Q = \omega_0 / \Delta\omega$, hoặc tính theo tỷ lệ giữa năng lượng tích trữ trong buồng cộng hưởng và năng lượng tiêu hao sau mỗi chu kỳ dao động ($Q = \omega_0 U / (-dU/dt)$).

5. Những loại vật liệu quang phi tuyến nào được đề cập để tích hợp trong cấu trúc PhCs nhằm tránh hư hại dưới laser cường độ cao?

Luận án đề cập đến các vật liệu điện môi phi tuyến bậc ba gồm $As_2S_3$ và polymer DDMEBT. Các vật liệu này có khả năng chịu được mật độ công suất laser xung lên tới $\sim\text{GW/cm}^2$ nhờ đặc tính không dẫn nhiệt, ngưỡng phá hủy nhiệt cao (khoảng $300^\circ\text{C}$ đối với $As_2S_3$) và thời gian kích hoạt của xung laser đủ ngắn để không gây tích tụ nhiệt phá hủy cấu trúc.

Kết luận

Luận án đã thực hiện nghiên cứu cơ bản có hệ thống về mô hình hóa, thiết kế và mô phỏng các cấu trúc tinh thể quang tử 1D và 2D ứng dụng cho linh kiện lưỡng trạng thái quang ổn định. Thông qua việc kết hợp các phương pháp toán học và mô phỏng số gồm CMT, PWE và FDTD, tác giả đã làm sáng tỏ ảnh hưởng của các tham số cấu trúc đến hệ số phẩm chất $Q$, vùng cấm quang PBG và ngưỡng năng lượng chuyển mạch quang. Các kết quả đạt được cung cấp cơ sở dữ liệu và cấu hình thiết kế cho việc phát triển các thế hệ linh kiện chuyển mạch và xử lý tín hiệu toàn quang trong các mạch vi quang - điện tử tích hợp.