Tổng quan nghiên cứu

Công nghệ nano và quang tử học hiện đại đang chứng kiến sự phát triển vượt bậc nhờ vào các kỹ thuật gia công siêu chính xác, trong đó phương pháp khắc laser trực tiếp (Direct Laser Writing - DLW) đóng vai trò then chốt với khả năng tạo cấu trúc ba chiều đạt độ phân giải dưới 100 nm. Tuy nhiên, các hệ thống khắc laser hai photon (Two-Photon Absorption - TPA) thương mại sử dụng nguồn xung femto giây thường có chi phí đầu tư rất cao, lên tới khoảng 400.000 EUR, cùng cấu hình quang học phức tạp. Nhằm giải quyết rào cản kinh tế và kỹ thuật này, giải pháp khắc laser hấp thụ một photon cực thấp (Low One-Photon Absorption - LOPA) với nguồn laser liên tục giá thành khoảng 100.000 EUR đã được phát triển, giúp tiết kiệm tới 75% chi phí đầu tư thiết bị.

Vấn đề cốt lõi đặt ra trong công nghệ này là sự biến dạng của trường điện từ tại tiêu điểm khi chùm tia hội tụ qua các vật kính có khẩu độ số cao (NA > 0,7) và hiện tượng quang sai cầu do sự lệch chiết suất giữa các môi trường truyền qua. Nghiên cứu của học viên Nguyễn Mạnh Hiếu, thực hiện năm 2019 tại Trường Đại học Công nghệ – Đại học Quốc gia Hà Nội dưới sự hướng dẫn của TS. Trần Quốc Tiến (Viện Khoa học Vật liệu – Viện Hàn lâm Khoa học và Công nghệ Việt Nam), tập trung vào việc tính toán, mô phỏng chính xác các thông số trường quang tại vùng hội tụ của vật kính khẩu độ số cao. Mục tiêu chính là thiết lập mô hình tính toán số hoàn chỉnh để kiểm soát kích thước điểm hội tụ dưới 1 µm3, định lượng ảnh hưởng của phân cực, chiết suất và độ mở vật kính, từ đó cung cấp cơ sở khoa học vững chắc cho việc chế tạo vật liệu nano và linh kiện quang tử 2D, 3D.

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

Khung lý thuyết áp dụng

Nghiên cứu dựa trên nền tảng lý thuyết nhiễu xạ sóng ánh sáng và lý thuyết vector Debye - Wolf về sự hội tụ của sóng điện từ qua thấu kính có khẩu độ số cao. Khác với quang học hình học hay lý thuyết nhiễu xạ vô hướng Fresnel - Kirchhoff (vốn chỉ phù hợp khi NA < 0,7), lý thuyết vector Debye - Wolf xem xét đầy đủ bản chất vector của điện trường và từ trường khi chùm tia bị bẻ cong ở góc lớn. Mô hình tính toán tích hợp phương trình khúc xạ Snell tổng quát và hệ số truyền qua Fresnel tại bề mặt phân cách giữa hai môi trường có chiết suất khác nhau.

Các khái niệm then chốt được phát triển bao gồm: Khẩu độ số (NA = n*sinα, trong đó n là chiết suất môi trường ngâm và α là nửa góc mở cực đại), Hàm lan truyền điểm (Point Spread Function - PSF), Bán độ rộng cực đại (Full Width at Half Maximum - FWHM), và Tỷ lệ co hình học (Aspect Ratio - AR = L/T, giữa kích thước dọc trục L và kích thước ngang T). Ngoài ra, nghiên cứu mở rộng khung lý thuyết Debye - Wolf cho môi trường hấp thụ có chiết suất phức ñ = n + iκ, giải thích cơ chế quang sai pha cảm ứng Φ(θ, φ, d) khi chùm tia đi sâu vào lớp vật liệu cảm quang.

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

Nghiên cứu sử dụng phương pháp mô phỏng số rời rạc hóa tích phân vector Debye trên nền tảng phần mềm MATLAB kết hợp các module tối ưu hóa bằng ngôn ngữ C.

  • Cỡ mẫu và miền không gian tính toán: Không gian quan sát 3D được thiết lập dạng lưới đồng nhất gồm 200 x 200 x 200 điểm (tương ứng Nx = Ny = Nz = 200) trong thể tích thực nghiệm 0,5 µm x 0,5 µm x 0,5 µm. Bước lưới không gian được chọn nhỏ hơn λ/10 (khoảng 53,2 nm đối với bước sóng λ = 532 nm), đảm bảo tiêu chuẩn hội tụ và độ mịn cấu trúc tương đương phương pháp sai phân hữu hạn miền thời gian (FDTD).
  • Phương pháp chia góc và lấy mẫu tích phân: Miền góc mở cực đại α (lên tới 67,5 độ) và góc phương vị φ (từ 0 đến 2π) được rời rạc hóa thành Nθ = 200 bước và Nφ = 200 bước tuyến tính.
  • Lý do lựa chọn phương pháp: Tích phân số Debye - Wolf cho phép tính toán trực tiếp sự chồng chập của hàng vạn sóng phẳng thành phần mà không làm mất đi các thành phần vector phân cực dọc (Ez). Nguồn dữ liệu mô phỏng được đối chiếu trực tiếp với kết quả thực nghiệm quang khắc LOPA tại bước sóng laser 532 nm và 355 nm để kiểm chứng độ chính xác.

Kết quả nghiên cứu và thảo luận

Những phát hiện chính

Mô phỏng 06 loại vật kính với các giá trị khẩu độ số khác nhau trong các môi trường không khí (n = 1,00), nước (n = 1,33) và dầu ngâm (n = 1,51) đã làm sáng tỏ 4 phát hiện quan trọng:

  • Tối ưu hóa kích thước voxel theo khẩu độ số: Khi tăng NA từ 0,65 lên 1,40, kích thước ngang FWHM của vùng hội tụ giảm hơn 50%, giúp thể tích hội tụ đạt mức giới hạn nhiễu xạ khoảng 0,4λ theo phương ngang và 1,33λ theo phương dọc trục z.
  • Tác động định hình của trạng thái phân cực: Phân cực tròn tạo ra đĩa sáng hội tụ đối xứng hoàn hảo quanh trục quang. Ngược lại, phân cực thẳng theo trục x khiến điểm hội tụ bị kéo dãn thành hình elip bất đối xứng với sai lệch kích thước giữa trục x và trục y từ 15% đến 20%. Phân cực hướng tâm tạo ra thành phần điện trường dọc Ez cực mạnh tại trung tâm tiêu điểm.
  • Định lượng quang sai do bất tương thích chiết suất: Khi ánh sáng truyền từ môi trường ngâm dầu (n1 = 1,51) vào lớp cảm quang có chiết suất khác biệt (n2 = 1,33 hoặc n2 = 1,60), tiêu điểm thực tế bị dịch chuyển dọc trục z một khoảng tỷ lệ thuận với độ sâu d, đồng thời cường độ đỉnh suy giảm từ 25% đến 40% kèm theo sự xuất hiện của các vân nhiễu xạ phụ.
  • Đột phá độ phân giải của kỹ thuật LOPA: Kết hợp hiệu ứng nhiệt cảm ứng quang với hấp thụ một photon cực thấp cho phép tạo ra cấu trúc polymer có kích thước thực tế đạt 57 nm và khoảng cách giữa hai vệt khắc đạt 300 nm mà không cần công đoạn ủ nhiệt 100°C sau chiếu xạ.

Thảo luận kết quả

Nguyên nhân chính dẫn đến sự biến dạng của vùng hội tụ ở vật kính NA cao là do sự đóng góp vượt trội của các tia nghiêng góc lớn (α > 60 độ), làm xuất hiện thành phần điện trường dọc trục z vốn bị triệt tiêu trong các hệ quang học NA thấp. Khi có sự không khớp chiết suất giữa môi trường dầu và đế mẫu, hàm sai pha phi tuyến làm méo mó mặt sóng cầu Gaussian, dẫn đến việc giãn rộng kích thước vệt khắc theo chiều sâu.

Dữ liệu tính toán được minh họa trực quan thông qua các bản đồ phân bố cường độ 2D trên các mặt phẳng cắt ngang (x2y2), mặt phẳng đứng (x2z2, y2z2) và đồ thị phân bố cường độ 1D dọc theo các trục đối xứng. Kết quả đối chứng cho thấy sự tương thích rất cao giữa mô hình vector lý thuyết và hình ảnh thực nghiệm chụp từ kính hiển vi đồng tiêu (CLSM), vượt trội hoàn toàn so với mô hình vô hướng (vốn gây sai số từ 30% đến 35% về kích thước tiêu điểm).

Đề xuất và khuyến nghị

  • Ứng dụng vật kính nhúng dầu khẩu độ số cao: Khuyến nghị các phòng thí nghiệm chế tạo nano sử dụng vật kính có NA từ 1,40 trở lên kết hợp môi trường dầu ngâm chiết suất n = 1,51 nhằm thu hẹp kích thước vệt khắc xuống dưới 200 nm, nâng cao độ phân giải quang học thêm 35% trong giai đoạn 2024-2025 (Chủ thể: Kỹ sư thiết kế hệ quang học và cán bộ nghiên cứu).
  • Chuẩn hóa hệ thống điều khiển phân cực chùm tia: Thiết lập tấm trễ pha λ/4 và λ/2 trên đường truyền quang trước khi chùm tia đi vào vật kính để duy trì trạng thái phân cực tròn hoặc phân cực hướng tâm đồng nhất, triệt tiêu hiện tượng méo elip của vệt khắc trong vòng 6 tháng tới (Chủ thể: Kỹ sư vận hành hệ thống khắc laser DLW).
  • Tích hợp module bù quang sai cầu chủ động: Phát triển thuật toán bù pha số học dựa trên bộ điều chế ánh sáng không gian (SLM) nhằm triệt tiêu quang sai pha Φ(θ, φ, d) khi khắc cấu trúc 3D ở độ sâu từ 10 µm đến 100 µm trong vật liệu cảm quang (Chủ thể: Nhóm nghiên cứu mô phỏng và kỹ thuật phần mềm).
  • Khai thác vật liệu lai vô cơ - hữu cơ sol-gel cho kỹ thuật LOPA: Ứng dụng vật liệu lai có khả năng điều chỉnh chiết suất linh hoạt và độ hấp thụ cực thấp tại bước sóng 532 nm nhằm chế tạo tinh thể quang tử 3D và linh kiện metamaterial với hiệu suất tăng 40% trong lộ trình 1-2 năm tới (Chủ thể: Các viện nghiên cứu vật liệu và doanh nghiệp chế tạo vi mạch).

Đối tượng nên tham khảo luận văn

  • Nghiên cứu sinh và học viên cao học ngành Vật lý quang học, Quang tử học, Vật liệu và Linh kiện Nano: Nắm vững kỹ thuật lập trình mô phỏng tích phân vector Debye - Wolf và phương pháp xử lý trường điện từ 3D trong MATLAB để ứng dụng vào các đề tài nghiên cứu chuyên sâu.
  • Kỹ sư thiết kế hệ thống thiết bị laser và quang khắc vi điện tử: Sử dụng các công thức tính toán khẩu độ số, hệ số truyền qua Fresnel và tối ưu hóa hệ thấu kính hội tụ mạnh trong việc chế tạo các máy khắc laser trực tiếp thế hệ mới.
  • Nhà nghiên cứu phát triển vật liệu Polymer và Hóa học quang học: Hiểu rõ cơ chế hấp thụ một photon cực thấp (LOPA) và hiệu ứng nhiệt cảm ứng quang để tối ưu hóa thành phần chất cảm quang âm và dương ở kích thước dưới 100 nm.
  • Chuyên gia chế tạo linh kiện cảm biến và siêu vật liệu Metamaterial: Ứng dụng quy trình chế tạo cấu trúc quang tử 2D, 3D và cấu trúc plasmonic màng vàng để nâng cao độ nhạy cho các cảm biến sinh học và hóa môi trường.

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

Tại sao cần sử dụng lý thuyết vector Debye thay vì lý thuyết vô hướng Fresnel?

Khi vật kính có khẩu độ số cao (NA > 0,7), các tia sáng hội tụ ở góc nghiêng rất lớn làm xuất hiện thành phần điện trường dọc trục z. Lý thuyết vô hướng bỏ qua thuộc tính vector và hướng truyền năng lượng của sóng, dẫn đến sai số kích thước tiêu điểm lên tới 35%, trong khi tích phân vector Debye mô tả chính xác 100% phân bố trường 3D.

Kỹ thuật khắc laser LOPA mang lại lợi thế kinh tế nào so với TPA?

Phương pháp TPA đòi hỏi nguồn laser xung femto giây đắt tiền với chi phí khoảng 400.000 EUR. Kỹ thuật LOPA chỉ sử dụng laser liên tục công suất thấp (vài mW) tại bước sóng hấp thụ yếu như 532 nm, giúp hạ giá thành hệ thống xuống khoảng 100.000 EUR mà vẫn tạo ra cấu trúc 3D có độ phân giải tương đương.

Khẩu độ số NA tác động như thế nào đến độ phân giải không gian của hệ DLW?

Khẩu độ số tỷ lệ thuận với góc mở α của chùm tia. Khi tăng NA từ 0,65 lên 1,40, kích thước điểm hội tụ ngang FWHM giảm hơn 50%, cho phép kích thước vệt hội tụ đạt giới hạn nhiễu xạ khoảng 0,4λ theo phương ngang và 1,33λ theo phương dọc, tạo điều kiện chế tạo cấu trúc nhỏ tới 57 nm.

Phân cực của chùm laser ảnh hưởng ra sao đến hình dạng vệt khắc?

Chùm laser phân cực tròn tạo ra điểm hội tụ đối xứng hoàn hảo quanh trục quang. Trong khi đó, phân cực thẳng theo phương x làm điểm hội tụ bị méo dạng elip bất đối xứng với độ chênh lệch giữa hai trục từ 15% đến 20%, ảnh hưởng tiêu cực đến tính đồng nhất của cấu trúc nano.

Hiện tượng lệch chiết suất giữa các môi trường gây ra tác hại gì?

Khi ánh sáng truyền qua mặt phân cách giữa hai môi trường có chiết suất khác nhau (như từ dầu n = 1,51 vào mẫu n = 1,33), quang sai cầu xuất hiện làm tiêu điểm bị dịch chuyển dọc trục, giảm từ 25% đến 40% cường độ sáng cực đại và làm biến dạng cấu trúc 3D khi khắc sâu.

Kết luận

  • Xây dựng thành công thuật toán mô phỏng 3D trường điện từ tại vùng hội tụ của vật kính khẩu độ số cao dựa trên tích phân vector Debye - Wolf mở rộng.
  • Định lượng chi tiết tác động của khẩu độ số NA (từ 0,65 đến 1,40), trạng thái phân cực và sự bất tương thích chiết suất đến hình dạng và kích thước tiêu điểm.
  • Chứng minh tính khả thi vượt trội của công nghệ khắc laser LOPA tại bước sóng 532 nm, cho phép tạo cấu trúc nano đạt độ phân giải thực tế 57 nm.
  • Cắt giảm 75% chi phí đầu tư thiết bị quang khắc 3D (từ 400.000 EUR xuống còn 100.000 EUR) nhờ việc thay thế laser xung femto giây bằng laser liên tục.
  • Thiết lập hệ thống dữ liệu mô phỏng và thực nghiệm chuẩn xác phục vụ chế tạo linh kiện quang tử, tinh thể quang tử và vật liệu metamaterial tại Việt Nam.

Lộ trình tiếp theo hướng tới việc mở rộng mô hình tính toán cho các chùm tia đặc biệt như Bessel, chùm tia xoáy pha (vortex) và tích hợp các thuật toán bù pha chủ động trong giai đoạn tới. Quý độc giả, nhà khoa học và doanh nghiệp công nghệ quan tâm đến giải pháp quang khắc nano chi phí thấp có thể tham khảo chi tiết toàn văn luận văn để cùng hợp tác phát triển các ứng dụng thực tiễn.