Tổng quan luận án

Luận án tiến sĩ thuộc chuyên ngành Vật lý quang học và Kỹ thuật quang tử, tập trung nghiên cứu sự truyền và phân bố trường điện từ trong vùng hội tụ của vật kính có khẩu độ số cao nhằm phục vụ công nghệ khắc laser trực tiếp (Direct Laser Writing - DLW) trong chế tạo cấu trúc nano. Về hướng tiếp cận, công trình kết hợp chặt chẽ giữa nghiên cứu lý thuyết, phát triển mô hình tính toán số và đối chứng thực nghiệm.

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

Trong lĩnh vực công nghệ nano và quang tử học, việc chế tạo các cấu trúc micro/nano 2D và 3D giữ vai trò trung tâm. Các công cụ truyền thống như kính hiển vi điện tử quét (SEM), kính hiển vi quét xuyên hầm (STM) hay kính hiển vi lực nguyên tử (AFM) dù có độ phân giải cao nhưng chi phí vận hành lớn, phức tạp và chủ yếu giới hạn ở bề mặt hoặc không gian 2D. Ngược lại, kính hiển vi quang học (OM) cho phép nghiên cứu và chế tạo linh hoạt trong không gian ba chiều với chi phí thấp hơn, nhưng lại chịu rào cản nghiêm ngặt bởi giới hạn nhiễu xạ ánh sáng (khoảng một nửa bước sóng).

Để vượt qua giới hạn này, nhiều công nghệ đã được đề xuất:

  • Phương pháp hấp thụ một photon truyền thống (OPA): Sử dụng laser liên tục giá thành thấp nhưng ánh sáng bị suy giảm cường độ nhanh theo độ sâu do vật liệu hấp thụ mạnh, dẫn đến việc không thể tạo cấu trúc 3D dày hay đồng đều (độ dày thường chỉ đạt khoảng 1 µm).
  • Phương pháp hấp thụ hai photon (TPA): Khắc phục được độ phân giải theo trục nhờ phản ứng phi tuyến bậc hai, nhưng đòi hỏi nguồn laser xung siêu ngắn (femto giây hoặc pico giây) công suất đỉnh cao, kéo theo hệ thống quang học rất phức tạp và chi phí thiết bị lên tới khoảng 400.000 euro.
  • Phương pháp làm nghèo phát xạ kích thích (STED) hoặc in quang học điều chế hấp thụ (AMOL): Cho phép thu nhỏ điểm sáng hiệu dụng nhưng đòi hỏi kết hợp đồng thời hai hoặc ba chùm laser đồng trục (chùm kích thích và chùm làm nghèo dạng vòng xoáy hoặc lệch pha 0/π), tạo ra thách thức lớn về điều chỉnh cơ quang học.
  • Phương pháp hấp thụ một photon cực thấp (LOPA): Do Phòng thí nghiệm Quang lượng tử và Phân tử (LPQM, UMR 8537 - CNRS, Đại học Sư phạm Cachan, Pháp) đề xuất, sử dụng laser liên tục công suất thấp (vài mW) tại bước sóng có độ hấp thụ cực yếu của vật liệu cảm quang kết hợp vật kính khẩu độ số cao, cho phép giảm chi phí hệ thống xuống khoảng 100.000 euro đồng thời tận dụng hiệu ứng nhiệt cảm ứng quang học.

Khoảng trống nghiên cứu được tác giả xác định là: Dưới điều kiện hội tụ rất mạnh của vật kính có khẩu độ số cao (NA > 0,7), các lý thuyết quang học vô hướng truyền thống không còn phản ánh chính xác cấu trúc trường sáng. Đặc biệt, sự phân bố trường điện từ khi ánh sáng hội tụ qua mặt phân cách vào một môi trường có tính hấp thụ (môi trường chiết suất phức) chưa được khảo sát đầy đủ bằng mô hình toán học chặt chẽ. Việc thiếu hụt các mô hình mô phỏng định lượng về ảnh hưởng của khẩu độ số, trạng thái phân cực, dạng chùm tia tới và sự không phối hợp chiết suất dẫn đến khó khăn trong việc kiểm soát kích thước và hình học vùng hội tụ trong công nghệ DLW.

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

  1. Mục tiêu lý thuyết: Mở rộng cơ sở toán học của tích phân Debye - Wolf để mô tả sự hội tụ của sóng điện từ qua mặt phân cách phẳng vào môi trường hấp thụ; xây dựng giải thuật tính toán số phân bố trường điện từ vector trong không gian hội tụ 3D.
  2. Mục tiêu mô phỏng: Khảo sát định lượng ảnh hưởng của các tham số quang học (khẩu độ số NA, chiết suất môi trường nhúng, dạng chùm tia tới, trạng thái phân cực) đến kích thước (độ bán rộng FWHM ngang và dọc), hình dạng và tỷ lệ co (Aspect Ratio - AR) của vùng hội tụ.
  3. Mục tiêu thực nghiệm: So sánh, đối chứng kết quả tính toán mô phỏng lý thuyết với các dữ liệu thực nghiệm quan sát được trên hệ khắc laser trực tiếp, từ đó định hướng kiểm soát thông số công nghệ trong chế tạo cấu trúc nano.

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

  • Đối tượng nghiên cứu: Phân bố cường độ và phân cực của trường điện từ (trường quang) trong vùng hội tụ của vật kính có khẩu độ số cao; hàm truyền điểm (Point Spread Function - PSF); các cơ chế quang khắc laser (OPA, TPA, LOPA).
  • Phạm vi không gian: Kích thước vùng tiêu điểm bước sóng phụ (sub-wavelength) trong phạm vi từ vài micromet xuống dưới 100 nm; miền mô phỏng rời rạc hóa 3D kích thước 0,5 µm × 0,5 µm × 0,5 µm.
  • Phạm vi tham số và môi trường: Bước sóng nghiên cứu λ = 532 nm (và tham chiếu 355 nm trong hệ DLW); khẩu độ số vật kính từ NA = 0,7 đến NA = 1,515; môi trường ngâm gồm không khí (n = 1), nước (n = 1,33) và dầu ngâm (n = 1,515); môi trường hấp thụ có chiết suất phức $\tilde{n}_2 = n_2 + \mathrm{i}\kappa$.

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

Luận án đã hệ thống hóa các mốc phát triển lý thuyết và thực nghiệm trong quang học trường hội tụ mạnh và công nghệ quang khắc:

Các hướng nghiên cứu lý thuyết trường hội tụ

  • Khởi xướng từ Peter Debye [89] vào đầu thế kỷ 20 về việc biểu diễn sóng hội tụ qua sự chồng chất của các sóng phẳng đồng nhất trong một góc mở khẩu độ.
  • Wolf (1959) [90] và công trình kinh điển Richards - Wolf (những năm 1950) [84, 85]: Đặt nền móng cho việc mở rộng tích phân Debye sang mô tả vector đầy đủ của trường điện từ trong vùng hội tụ của thấu kính có NA cao, tính đến thuộc tính phân cực và sự biến đổi mặt sóng cầu Gaussian.
  • Gasper và cộng sự (1976) [92]: Đưa ra nghiệm hoàn chỉnh và các xấp xỉ tiệm cận cho trường điện và từ của sóng điện từ tùy ý truyền qua bề mặt phẳng trong môi trường đồng nhất đẳng hướng.
  • Ling và Lee [95]: Nghiên cứu sự hội tụ của sóng điện từ qua mặt giao diện dựa trên điều kiện biên phân bố dòng và phổ góc sóng phẳng theo cách tiếp cận bán hình học.
  • Ji và Hongo [96]: Khảo sát bài toán nguồn điểm qua giao diện điện môi hình cầu bằng phương pháp Maslov.
  • Stamnes (1986) [99]: Tổng hợp các hướng tiếp cận lý thuyết về nhiễu xạ và hội tụ sóng quang học.
  • Hell và cộng sự (1993) [100]: Nghiên cứu sự hội tụ trong vật liệu có chiết suất không đồng nhất bằng tích phân Fresnel - Kirchhoff, tách thành phần vector điện và phân cực để đánh giá quang sai cầu trong kính hiển vi đồng tiêu (tuy nhiên cách tiếp cận này tồn tại hạn chế do chỉ xét định lý Green mà chưa bảo đảm tính liên tục của thành phần tiếp tuyến từ trường).
  • Török và cộng sự (1995) [101]: Thiết lập nghiệm giải tích chặt chẽ cho sự truyền sóng điện từ qua giao diện phẳng giữa các môi trường có chiết suất khác nhau bằng hình thức ma trận, mở rộng lý thuyết Richards - Wolf khi có quang sai lớn.
  • Wiersma và Visser (1996) [102]: Phát triển độc lập phương pháp tính toán trường trong môi trường thứ hai bằng lý thuyết vector kéo theo (lý thuyết m).
  • Leutenegger và cộng sự (2006) [103], Lin và cộng sự [104]: Tối ưu hóa tốc độ tính toán số tích phân Debye bằng thuật toán biến đổi Fourier nhanh, giúp giảm thời gian tính toán tới 40 lần so với phương pháp tích phân trực tiếp.
  • Fedorov và Nakajima [118]: Mở rộng định luật khúc xạ Snell tổng quát cho sóng phẳng tại mặt phân cách môi trường.

Khoảng trống luận án lựa chọn giải quyết

Mặc dù bài toán sóng truyền qua giao diện phẳng giữa hai môi trường điện môi trong suốt đã được nghiên cứu sâu, song hiện tượng sóng điện từ hội tụ qua thấu kính NA cao vào một môi trường hấp thụ (nơi chiết suất mang giá trị phức và vector sóng trở thành vector phức) chưa được xử lý toán học thấu đáo. Trong công nghệ khắc laser trực tiếp, vật liệu cảm quang luôn có tính hấp thụ nhất định làm suy giảm biên độ sóng theo độ sâu, ảnh hưởng trực tiếp đến ngưỡng năng lượng polyme hóa. Luận án đã lấp đầy khoảng trống này bằng cách mở rộng biểu diễn Debye - Wolf cho môi trường hấp thụ và phân tích chi tiết ảnh hưởng của phân cực vector lên cấu trúc tiêu điểm.

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

Khung lý thuyết và khái niệm cốt lõi

  1. Nhiễu xạ và tiêu chuẩn phân giải:
    • Nguyên lý Huygens - Fresnel và phương pháp đới cầu Fresnel: Giải thích định tính và định lượng sự sai lệch phương truyền thẳng của sóng ánh sáng khi qua các lỗ tròn, đĩa tròn, khe hẹp và cách tử phẳng.
    • Nhiễu xạ trên tinh thể: Định luật Vulf - Bragg với hiệu quang lộ cực đại $\Delta L = 2d \sin\varphi = k\lambda$.
    • Tiêu chuẩn phân giải Rayleigh: Giới hạn phân giải ngang tối thiểu $T/2 = 0,61\lambda/\mathrm{NA}$ và dọc $L/2 = 2\lambda/\mathrm{NA}^2$, với khẩu độ số $\mathrm{NA} = n \sin\alpha$ ($n$ là chiết suất, $\alpha$ là nửa góc mở cực đại).
  2. Lý thuyết vector Debye - Richards - Wolf:
    • Ánh sáng ở lối ra của vật kính khẩu độ số cao biến đổi thành mặt sóng cầu Gaussian bán kính $f$ (tiêu cự vật kính).
    • Biểu thức tích phân trường điện từ tại điểm quan sát bất kỳ $P(x_2, y_2, z_2)$ trong hệ tọa độ Descartes: $$E(x_2, y_2, z_2) = -\frac{\mathrm{i}C}{\lambda} \int_{0}^{\alpha} \int_{0}^{2\pi} \sin\theta A(\theta, \varphi) B(\theta, \varphi) \mathbf{P}(\theta, \varphi) \exp\left[\mathrm{i}kn(z_2 \cos\theta + x_2 \sin\theta \cos\varphi + y_2 \sin\theta \sin\varphi)\right] \mathrm{d}\theta \mathrm{d}\varphi$$ Trong đó: $A(\theta, \varphi)$ là hàm biên độ chùm tia tới; $B(\theta, \varphi) = \sqrt{\cos\theta}$ là hệ số lọc quang bảo toàn năng lượng qua thấu kính; $\mathbf{P}(\theta, \varphi) = \mathbf{T}(\theta, \varphi) \mathbf{P}_0(\theta, \varphi)$ là trạng thái phân cực trong vùng hội tụ với $\mathbf{T}(\theta, \varphi)$ là toán tử ma trận biến đổi không gian $3 \times 3$.
  3. Khai triển Debye - Wolf mở rộng cho môi trường hấp thụ:
    • Không gian ảnh gồm môi trường 1 không suy hao (chiết suất thực $n_1$) và môi trường 2 có suy hao (chiết suất phức $\tilde{n}_2 = n_2 + \mathrm{i}\kappa$).
    • Vector sóng trong môi trường 2 là vector phức $\mathbf{k}_2 = \mathbf{k}'_2 + \mathrm{i}\mathbf{k}''_2$, trong đó phần ảo $\mathbf{k}''_2$ đặc trưng cho sự suy giảm biên độ luôn vuông góc với mặt phân cách phẳng.
    • Trường điện từ được biểu diễn qua ma trận Jacobi $\mathbf{J}0 = (k_1 / k'2)^2$ và thỏa mãn phương trình sóng không thuần nhất: $$E_2(x, y, z) = -\frac{\mathrm{i}}{2\pi k_1} \iint{\Omega_1} \frac{k'2}{s{1z}} \mathbf{T}^{(2)} a(s{1x}, s_{1y}) \exp[-\mathrm{i}d(k_1 s_{1z} - k'2 s'{2z})] \exp(\mathrm{i}k'2 s'{2z} z) \exp[\mathrm{i}k_1(s_{1x}x + s_{1y}y)] \exp[-k''2(z + d)] \mathrm{d}s{1x} \mathrm{d}s_{1y}$$

Phương pháp tính toán số và công cụ mô phỏng

  • Thuật toán số: Rời rạc hóa tích phân Debye liên tục thành tổng hữu hạn các sóng phẳng. Góc cực $\theta$ được chia thành $N_\theta$ bước đều từ 0 đến $\alpha$ với bước nhảy $\Delta\theta = \alpha / N_\theta$; góc phương vị $\varphi$ chia thành $N_\varphi$ bước đều từ 0 đến $2\pi$ với bước nhảy $\Delta\varphi = 2\pi / N_\varphi$.
  • Lưới không gian: Vùng quan sát được mô hình hóa theo dạng lưới 3D $(N_x, N_y, N_z)$ tương tự phương pháp sai phân hữu hạn trong miền thời gian (FDTD) với khoảng cách mắt lưới nhỏ hơn $\lambda / 10$.
  • Tham số cài đặt mô phỏng: $N_\theta = N_\varphi = 200$; $N_x = N_y = N_z = 200$; kích thước thể tích khảo sát là $0,5,\mu\text{m} \times 0,5,\mu\text{m} \times 0,5,\mu\text{m}$; góc mở $\alpha \approx 67,5^\circ$.
  • Môi trường phần mềm: Các đoạn mã thuật toán được lập trình bằng ngôn ngữ C tích hợp chạy trên nền phần mềm MATLAB (chi tiết mã nguồn trong Phụ lục A của văn bản).

Nguồn dữ liệu và thiết bị thực nghiệm

  • Dữ liệu thực nghiệm được thu thập từ hệ thống khắc laser trực tiếp (DLW) gồm: Nguồn laser bước sóng $355\text{ nm}$ hoặc $532\text{ nm}$; hệ dẫn quang có thấu kính chuẩn trực, bản chậm $\lambda/2$, bộ tách chùm phân cực (PBS) và bản chậm $\lambda/4$; vật kính hội tụ khẩu độ số cao; hệ dịch chuyển mẫu 3D bằng gốm áp điện (PZT); đầu thu huỳnh quang dùng photodiode tuyết lở (APD); hệ thống máy tính điều khiển qua phần mềm Igor và LabVIEW.

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

CHƯƠNG I: CƠ SỞ LÝ THUYẾT

Chương I tập trung thiết lập nền tảng quang sóng và công nghệ chế tạo cấu trúc vi mô, nano bằng kỹ thuật quang khắc laser:

  • Khảo sát bản chất sóng của ánh sáng qua các hiện tượng nhiễu xạ sóng cầu và sóng phẳng. Phân tích định lượng cường độ sáng qua lỗ tròn, đĩa tròn, khe hẹp đơn và cách tử phẳng gồm $N$ khe hẹp. Tại cách tử, giữa hai cực đại chính xuất hiện $N - 1$ cực tiểu phụ và $N - 2$ cực đại phụ.
  • Giới thiệu cấu hình nguyên lý của hệ khắc laser trực tiếp (DLW) với khả năng dịch chuyển quét không gian 3 chiều qua bộ áp điện PZT, cho phép tạo cấu trúc có độ phân giải dưới $100\text{ nm}$.
  • Phân loại và so sánh cơ chế của hai nhóm vật liệu cảm quang polymer:
    • Cảm quang âm: Trùng hợp và đóng rắn tại vùng chiếu sáng, bền vững sau khi rửa dung môi, chuyên dùng cho cấu trúc quang tử 3D.
    • Cảm quang dương: Bị cắt mạch và rửa trôi tại vùng chiếu sáng, thường chỉ áp dụng cho cấu trúc 1D, 2D hoặc làm khuôn mẫu lift-off.
  • Phân tích cơ chế kích thích:
    • OPA: Dùng laser liên tục giá rẻ nhưng giới hạn bề dày khoảng $1,\mu\text{m}$.
    • TPA/MPA: Dùng xung femto giây tại bước sóng khoảng $800\text{ nm}$ (như laser Ti-Sapphire), kích thích phi tuyến bậc ba tại tiêu điểm kích thước $< 1,\mu\text{m}^3$, chi phí hệ thống khoảng 400.000 euro.
    • LOPA: Kỹ thuật do phòng thí nghiệm LPQM đề xuất, sử dụng laser liên tục bước sóng khả kiến tại rìa hấp thụ yếu của polymer, giảm chi phí còn 100.000 euro. Quá trình quang khắc LOPA tạo ra hiệu ứng nhiệt cảm ứng tại chỗ giúp hoàn thành liên kết mà không cần ủ nhiệt ngoài ở $100^\circ\text{C}$, đồng thời hạn chế tích lũy năng lượng lan truyền sang vùng lân cận.
    • Ứng dụng nhiệt cảm ứng trong LOPA trên chất cảm quang dương đã tạo ra cấu trúc vượt qua giới hạn nhiễu xạ với kích thước đạt tới $57\text{ nm}$ và khoảng cách hai cấu trúc là $300\text{ nm}$. Ngoài ra, phương pháp còn cho phép in màu plasmonic trên màng vàng mỏng thông qua việc laser quét gây co màng tạo hạt nano vàng cộng hưởng ở vùng $500\text{ nm}$.
  • Tổng hợp các ứng dụng của cấu trúc quang tử chế tạo bằng DLW: Lò xo vi cơ, buồng cộng hưởng vi mô, tinh thể quang tử điều chỉnh dải cấm, giá đỡ tế bào sinh học và vật liệu nhân tạo metamaterials (MMs) 3D có khả năng hấp thụ hoàn toàn sóng điện từ.

CHƯƠNG II: NGHIÊN CỨU CƠ BẢN VỀ PHÂN BỐ ÁNH SÁNG TRONG VÙNG HỘI TỤ CỦA MỘT VẬT KÍNH CÓ KHẨU ĐỘ SỐ CAO

Chương II phát triển khung giải tích và thuật toán số mô tả cấu trúc trường quang trong không gian hội tụ:

  • Phân tích sự hình thành trường giao thoa từ các tia sáng nhiễu xạ qua khẩu độ tròn. Vùng hội tụ thực tế có dạng hình elip đối xứng trục với kích thước ngang $T$ và dọc $L$, tạo nên tỷ lệ co $\mathrm{AR} = L/T$.
  • Đánh giá giới hạn của tiêu chuẩn Rayleigh: Công thức vô hướng xem vùng hội tụ là đĩa Airy bỏ qua tính chất vector và chiều truyền năng lượng (vector Poynting), không thể giải thích hiện tượng vùng hội tụ bị rỗng hoặc bất đối xứng dưới khẩu độ số cao ($\mathrm{NA} > 0,7$).
  • Thiết lập hệ phương trình vector Richards - Wolf với hệ tọa độ cầu liên kết với tọa độ Descartes qua vector đơn vị $\mathbf{s} = (\sin\theta \cos\varphi, \sin\theta \sin\varphi, \cos\theta)$.
  • Phân tích ma trận biến đổi phân cực $\mathbf{T}(\theta, \varphi)$ qua hai phép quay tọa độ $R$ và $C$: $$\mathbf{P}(\theta, \varphi) = \begin{bmatrix} 1 + (\cos\theta - 1)\cos^2\varphi & (\cos\theta - 1)\cos\varphi \sin\varphi & -\sin\theta \cos\varphi \ (\cos\theta - 1)\cos\varphi \sin\varphi & 1 + (\cos\theta - 1)\sin^2\varphi & -\sin\theta \sin\varphi \ \sin\theta \cos\varphi & -\sin\theta \sin\varphi & \cos\theta \end{bmatrix} \begin{bmatrix} p_x(\theta, \varphi) \ p_y(\theta, \varphi) \ p_z(\theta, \varphi) \end{bmatrix}$$
  • Xây dựng mô hình Debye - Wolf mở rộng cho trường hợp hội tụ qua giao diện phẳng vào môi trường hấp thụ: Sử dụng định luật Snell phức, thiết lập các điều kiện biên cho thành phần tiếp tuyến của điện trường và từ trường, đưa thừa số suy giảm $\exp[-k''_2(z + d)]$ trực tiếp vào nhân tích phân.
  • Đề xuất thuật toán rời rạc hóa tổng sóng phẳng hai lớp trên không gian lưới 3D, lập trình hoàn chỉnh bằng ngôn ngữ C chạy trên MATLAB.

CHƯƠNG III: KẾT QUẢ MÔ PHỎNG VỀ PHÂN BỐ TRƯỜNG QUANG TRONG VÙNG HỘI TỤ CỦA VẬT KÍNH CÓ KHẨU ĐỘ SỐ CAO

Chương III trình bày chi tiết các dữ liệu mô phỏng định lượng thu được từ thuật toán đã xây dựng:

  • Khảo sát 06 loại vật kính khẩu độ số cao trong các môi trường khác nhau với chùm tia tới phân cực tròn đồng nhất, bước sóng $\lambda = 532\text{ nm}$:
    • Khi tăng NA, kích thước ngang ($T$) và kích thước dọc ($L$) của vùng hội tụ (tính theo độ bán rộng FWHM) đều giảm mạnh. Kích thước dọc $L$ giảm từ $1,8,\mu\text{m}$ xuống còn khoảng $0,4,\mu\text{m}$; kích thước ngang $T$ giảm từ $0,6,\mu\text{m}$ xuống dưới $0,3,\mu\text{m}$.
    • Tỷ lệ co $\mathrm{AR} = L/T$ giảm rõ rệt từ mức xấp xỉ 3,5 - 4,0 xuống khoảng 2,0 tại $\mathrm{NA} = 1,4$.
    • Mặc dù góc mở $\alpha$ của vật kính trong không khí nhỏ hơn trong nước, nhưng kích thước vùng hội tụ với $\mathrm{NA} = 0,9$ ($n = 1$) lại nhỏ hơn so với vật kính có góc $\alpha$ tương đương trong môi trường nước do bước sóng hiệu dụng bị co ngắn trong môi trường có chiết suất cao.
    • So sánh giá trị FWHM tính từ lý thuyết vector Debye luôn nhỏ hơn giá trị dự đoán từ công thức vô hướng Rayleigh.

Bảng tổng hợp tham số các chế độ chùm tia đầu vào được khảo sát:

Chế độ chùm tia Biểu thức giải tích biên độ $A(\theta, \varphi)$ Đặc điểm phân bố tại vùng hội tụ
Chùm đồng nhất $A(\theta, \varphi) = 1$ Kích thước điểm lấy nét nhỏ nhất, đối xứng
Chùm Gaussian $A(\theta, \varphi) = \exp\left(-\frac{x_1^2 + y_1^2}{w^2}\right)$ Điểm lấy nét có kích thước lớn hơn chùm đồng nhất, tỷ lệ co nhỏ hơn, có dạng gần đẳng hướng
Chùm Laguerre - Gaussian (bậc 1, $p=1$) $A(\theta, \varphi) = \beta \frac{\sin\theta}{\sin\alpha} \exp\left[-\left(\beta \frac{\sin\theta}{\sin\alpha}\right)^2\right] L_1^1\left[2\left(\beta \frac{\sin\theta}{\sin\alpha}\right)^2\right]$ Cường độ đỉnh trên trục quang bị triệt tiêu/chia làm hai đỉnh dọc trục $z_2$, xuất hiện cấu trúc vòng bên

Ghi chú: $w$ là eo chùm Gaussian, $\beta$ là tỷ lệ giữa bán kính vật kính và eo chùm (thiết lập $\beta = 2,0$ với đa thức $L_1^1(x) = 2 - x$).

  • Nghiên cứu trạng thái phân cực của chùm tia tới (tại $\mathrm{NA} = 1,515$, $\lambda = 532\text{ nm}$):

Bảng trạng thái vector phân cực đầu vào $\mathbf{P}_0 = [p_x, p_y, p_z]^T$:

Trạng thái phân cực Thành phần $p_x(\theta, \varphi)$ Thành phần $p_y(\theta, \varphi)$ Thành phần $p_z(\theta, \varphi)$
Tuyến tính theo trục x $1$ $0$ $0$
Tuyến tính theo trục y $0$ $1$ $0$
Phân cực tròn phải (RC) $1/\sqrt{2}$ $\mathrm{i}/\sqrt{2}$ $0$
Phân cực tròn trái (LC) $1/\sqrt{2}$ $-\mathrm{i}/\sqrt{2}$ $0$
Phân cực hướng tâm (Radial) $\cos\varphi$ $\sin\varphi$ $0$
Phân cực phương vị (Azimuthal) $-\sin\varphi$ $\cos\varphi$ $0$
  • Hiện tượng xuất hiện thành phần điện trường dọc $E_{z2}$ khi chùm phân cực thẳng theo trục x:
    • Khi chùm phân cực tuyến tính $x$ đi qua vật kính có NA cao, sự bẻ cong mạnh của các tia sáng làm xoay vector phân cực, tạo ra thành phần trường điện từ dọc trục quang $E_{z2}$ đáng kể tại vùng hội tụ bên cạnh thành phần ngang $E_{x2}$.
    • Thành phần $E_{x2}$ giữ nguyên pha qua mọi hướng nên đóng góp chính vào cường độ tâm.
    • Ngược lại, thành phần $E_{z2}$ tạo bởi các tia từ nửa trên trục x (hướng dương) ngược pha $\pi$ so với các tia từ nửa dưới trục x, dẫn đến cấu trúc phân bố cường độ tổng hợp bị kéo dài bất đối xứng dọc theo phương phân cực của ánh sáng tới trên mặt phẳng $(x_2 y_2)$.
  • Đối chiếu mô phỏng và thực nghiệm: Các kết quả tính toán số về dạng hình học và kích thước FWHM của điểm hội tụ phù hợp với quan sát thực nghiệm khi thay đổi khẩu độ số và phân cực chùm sáng.

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

Đóng góp về mặt lý luận khoa học

  1. Luận án đã mở rộng thành công lý thuyết nhiễu xạ vector Debye - Richards - Wolf cho trường hợp chùm ánh sáng hội tụ qua mặt phân cách phẳng vào môi trường hấp thụ có chiết suất phức ($\tilde{n} = n + \mathrm{i}\kappa$). Mô hình giải quyết đồng thời sự chồng chất sóng phẳng không thuần nhất và định luật khúc xạ Snell phức.
  2. Chứng minh một cách định lượng sự hình thành và vai trò của thành phần điện trường dọc $E_z$ trong điều kiện hội tụ khẩu độ số cao ($\mathrm{NA} > 0,7$), chỉ rõ nguyên nhân gây ra sự biến dạng hình học bất đối xứng của hàm lan truyền điểm (PSF) khi sử dụng chùm phân cực tuyến tính.
  3. Làm sáng tỏ sự sai lệch giữa tiêu chuẩn Rayleigh vô hướng và lý thuyết vector Debye, khẳng định tính bắt buộc của việc sử dụng mô hình vector Debye trong các bài toán quang học nano và quang khắc độ phân giải cao.

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

  1. Xây dựng và hoàn thiện chương trình tính toán số tối ưu bằng ngôn ngữ C chạy trên môi trường MATLAB, cho phép mô phỏng chính xác phân bố cường độ 3D của trường quang trong vùng tiêu điểm với các tham số linh hoạt về khẩu độ số, môi trường chiết suất, mode chùm tia (Gaussian, Laguerre - Gaussian) và các trạng thái phân cực khác nhau.
  2. Cung cấp cơ sở dữ liệu tính toán định lượng về kích thước FWHM ($T, L$) và tỷ lệ co ($\mathrm{AR} = L/T$) phục vụ việc tối ưu hóa quy trình khắc laser trực tiếp (DLW). Kết quả chỉ ra rằng chùm Gaussian thích hợp để tạo điểm lấy nét đẳng hướng, trong khi chùm Laguerre - Gaussian bậc 1 có thể tạo điểm hội tụ có cường độ triệt tiêu ở tâm phục vụ kỹ thuật siêu phân giải.
  3. Hỗ trợ cơ sở lý thuyết và công nghệ cho phương pháp hấp thụ một photon cực thấp (LOPA), giúp chế tạo các cấu trúc quang tử 2D, 3D với chi phí thiết bị thấp (khoảng 100.000 euro) thay thế cho hệ TPA xung femto giây đắt đỏ (400.000 euro), đồng thời tận dụng hiệu ứng nhiệt cảm ứng quang học để chế tạo cấu trúc nano kích thước tới $57\text{ nm}$ trên chất cảm quang dương.

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

Hạn chế của nghiên cứu

  • Mô hình tích phân Debye giả định thấu kính không có quang sai ($\Phi(s_x, s_y) = 0$), ngoại trừ quang sai sinh ra do không phối hợp chiết suất tại mặt phân cách phẳng.
  • Kỹ thuật hội tụ bằng một vật kính đơn lẻ luôn tạo ra vùng lấy nét có kích thước dọc ($L$) lớn hơn kích thước ngang ($T$), chưa thể đạt được độ phân giải trục hoàn toàn tương đương độ phân giải ngang nếu không bổ sung hệ quang học đối xứng.
  • Trong mô phỏng chùm Laguerre - Gaussian, sự chênh lệch cường độ thực tế giữa vòng trong và vòng ngoài khiến cường độ tại tâm chưa đạt mức triệt tiêu hoàn toàn bằng 0 như lý thuyết lý tưởng.

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

  • Mở rộng mô phỏng cho hệ thống quang học 4Pi (sử dụng hai vật kính đặt đối diện nhau) nhằm can thiệp vào trường ánh sáng dọc trục quang, thu nhỏ kích thước dọc $L$ để triệt tiêu tỷ lệ co $\mathrm{AR}$.
  • Khảo sát sâu hơn sự tương tác giữa trường phân bố hội tụ với các màng kim loại mỏng để tối ưu hóa quá trình in màu plasmonic và chế tạo vật liệu biến tính metamaterials 3D.
  • Phát triển các cấu trúc tinh thể quang tử tích hợp plasmon bề mặt ứng dụng làm cảm biến quang học có độ nhạy cao phục vụ phát hiện hóa chất độc hại trong môi trường và cảm biến sinh học tại Việt Nam.

Giá trị tham khảo

Luận án là tài liệu tham khảo chuyên sâu cho các đối tượng và định hướng nghiên cứu sau:

  • Nghiên cứu sinh và học viên cao học ngành Quang học, Vật lý kỹ thuật, Khoa học vật liệu và Kỹ thuật nano: Tham khảo phương pháp toán học mở rộng tích phân Debye - Wolf và thuật toán rời rạc hóa mô phỏng trường quang vector.
  • Các kỹ sư và nhóm nghiên cứu phòng thí nghiệm quang khắc: Sử dụng các bảng tra cứu tham số FWHM, tỷ lệ co $\mathrm{AR}$ và ảnh hưởng của phân cực để tối ưu hóa nguồn laser, lựa chọn vật kính có NA phù hợp trong hệ DLW, LOPA, TPA.
  • Các nhà nghiên cứu phát triển linh kiện quang điện tử và cảm biến: Ứng dụng quy trình chế tạo cấu trúc nano polymer và cấu trúc hạt kim loại plasmonic trong thiết kế cảm biến quang hóa môi trường và linh kiện quang tử nano.

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

1. Tại sao tiêu chuẩn phân giải Rayleigh không còn chính xác khi áp dụng cho vật kính có khẩu độ số cao?

Tiêu chuẩn Rayleigh được xây dựng dựa trên lý thuyết nhiễu xạ vô hướng Fresnel - Kirchhoff đối với đĩa Airy của khẩu độ tròn. Lý thuyết này bỏ qua tính chất vector của ánh sáng, bỏ qua hướng lan truyền năng lượng (vector Poynting) và giả định góc mở nhỏ. Khi vật kính có khẩu độ số cao ($\mathrm{NA} > 0,7$), mặt sóng bị uốn cong thành mặt cầu sâu, các tia sáng hội tụ ở góc lớn làm quay vector phân cực và xuất hiện thành phần điện trường dọc $E_z$. Do đó, chỉ có lý thuyết vector Debye - Richards - Wolf mới phản ánh đúng hình dạng và kích thước thực của vùng hội tụ.

2. Sự khác biệt căn bản về mặt kỹ thuật và chi phí giữa phương pháp OPA, TPA và LOPA trong khắc laser trực tiếp là gì?

  • OPA dùng laser liên tục giá rẻ nhưng vật liệu hấp thụ mạnh, ánh sáng suy giảm nhanh theo độ sâu nên chỉ chế tạo được cấu trúc 1D, 2D mỏng khoảng $1,\mu\text{m}$.
  • TPA dùng laser xung siêu ngắn (femto giây) công suất đỉnh lớn tại bước sóng không hấp thụ tuyến tính (~$800\text{ nm}$), kích thích phi tuyến bậc ba giúp tạo cấu trúc 3D nhưng hệ thống phức tạp và giá thành cao (~400.000 euro).
  • LOPA dùng laser liên tục công suất thấp (vài mW) tại dải bước sóng hấp thụ cực yếu của vật liệu kết hợp vật kính NA cao. Ánh sáng xuyên sâu hàng trăm micromet mà không đổi dạng tiêu điểm, chi phí hệ thống giảm xuống khoảng 100.000 euro và tận dụng được hiệu ứng nhiệt cảm ứng quang học.

3. Tại sao chùm ánh sáng phân cực tuyến tính theo trục x lại tạo ra vùng hội tụ bị kéo dài bất đối xứng?

Khi chùm phân cực thẳng theo trục x đi qua vật kính có NA cao, các tia sáng bị bẻ cong hội tụ mạnh về tiêu điểm. Sự quay của vector điện trường tạo ra thành phần phân cực dọc trục quang $E_{z2}$. Trong khi thành phần ngang $E_{x2}$ giữ cùng pha thì thành phần $E_{z2}$ từ nửa trên và nửa dưới trục x lại ngược pha nhau một góc $\pi$. Sự giao thoa và phân bố của các thành phần vector này làm cho tổng cường độ phát xạ tại mặt phẳng tiêu cự $(x_2 y_2)$ bị kéo dài dọc theo phương phân cực x của chùm tới.

4. Thuật toán mô phỏng trong luận án đã rời rạc hóa tích phân Debye như thế nào?

Tích phân Debye liên tục được chuyển đổi thành tổng hai lớp hữu hạn của các sóng phẳng: góc $\theta$ được chia thành $N_\theta = 200$ bước từ 0 đến $\alpha$ với bước $\Delta\theta = \alpha/N_\theta$; góc phương vị $\varphi$ chia thành $N_\varphi = 200$ bước từ 0 đến $2\pi$ với bước $\Delta\varphi = 2\pi/N_\varphi$. Không gian quan sát được chia lưới 3D $(N_x, N_y, N_z) = (200, 200, 200)$ với kích thước mắt lưới nhỏ hơn $\lambda/10$ trên thể tích $0,5,\mu\text{m} \times 0,5,\mu\text{m} \times 0,5,\mu\text{m}$, được lập trình bằng ngôn ngữ C chạy trên MATLAB.

5. Ý nghĩa của tỷ lệ co (Aspect Ratio - AR) trong vùng hội tụ là gì và nó thay đổi thế nào theo NA?

Tỷ lệ co $\mathrm{AR} = L/T$ là tỷ số giữa kích thước tiêu điểm theo chiều dọc ($L$) và chiều ngang ($T$) tính theo độ bán rộng FWHM. Vì một vật kính đơn chỉ hội tụ ánh sáng từ một phía dọc theo trục quang nên $L$ luôn lớn hơn $T$. Khi tăng khẩu độ số NA từ 0,7 lên 1,4, kích thước dọc $L$ giảm mạnh hơn kích thước ngang $T$, dẫn đến tỷ lệ co AR giảm từ khoảng 3,5 - 4,0 xuống còn khoảng 2,0, giúp vùng hội tụ tiệm cận hình dạng đẳng hướng cần thiết cho quang khắc 3D.

Kết luận

Luận án đã giải quyết hoàn chỉnh bài toán lý thuyết và mô phỏng số về phân bố trường quang vector trong vùng hội tụ của vật kính có khẩu độ số cao khi truyền vào môi trường hấp thụ. Bằng việc mở rộng tích phân Debye - Wolf và đối chứng với dữ liệu thực nghiệm, công trình đã làm sáng tỏ ảnh hưởng cốt lõi của khẩu độ số, trạng thái phân cực và các mode chùm tia tới cấu trúc tiêu điểm bước sóng phụ. Những kết quả đạt được cung cấp cơ sở khoa học và tham số kỹ thuật chuẩn xác, trực tiếp phục vụ việc tối ưu hóa công nghệ khắc laser trực tiếp và phát triển kỹ thuật hấp thụ một photon cực thấp (LOPA) trong chế tạo vật liệu nano và linh kiện quang tử.