Tổng quan nghiên cứu

Vật liệu biến hóa (metamaterials) có chiết suất âm là một trong những đột phá khoa học nổi bật của vật lý hiện đại kể từ khi Victor Veselago đặt nền móng lý thuyết vào năm 1968 và David R. Smith cùng các cộng sự chế tạo thành công mẫu thực nghiệm đầu tiên vào năm 2000. Tuy nhiên, rào cản lớn nhất ngăn cản việc thương mại hóa và ứng dụng rộng rãi loại vật liệu nhân tạo này là dải tần hoạt động có chiết suất âm thường rất hẹp (thường dưới 5% tần số trung tâm) do phụ thuộc vào cơ chế cộng hưởng từ cục bộ, đồng thời tính chất điện từ bị suy giảm nghiêm trọng khi góc phân cực của sóng tới thay đổi.

Trước thách thức đó, luận văn thạc sĩ chuyên ngành Quang học tại Trường Đại học Khoa học - Đại học Thái Nguyên do học viên Nguyễn Thị Mây thực hiện dưới sự hướng dẫn khoa học của TS. Nguyễn Thị Hiền tập trung vào mục tiêu nghiên cứu, thiết kế và tối ưu hóa cấu trúc vật liệu biến hóa mới nhằm mở rộng đồng thời vùng có độ từ thẩm âm ($\mu < 0$) và chiết suất âm ($n < 0$). Nghiên cứu được triển khai trong phạm vi dải sóng vi ba và Radar băng Ku từ 12,0 GHz đến 18,0 GHz vào năm 2018.

Ý nghĩa khoa học và thực tiễn của công trình thể hiện ở việc đề xuất thành công cấu trúc cặp đĩa hai lớp và lưới đĩa đối xứng tròn, giúp mở rộng dải tần từ thẩm âm thêm hơn 1,4 GHz so với cấu trúc một lớp truyền thống. Đồng thời, cấu trúc duy trì độ truyền qua cao trên -2 dB và hoàn toàn không phụ thuộc vào góc phân cực sóng điện từ khi quay từ 0 độ đến 30 độ. Đây là tiền đề vật lý quan trọng để phát triển các thiết bị siêu thấu kính phẳng vượt giới hạn nhiễu xạ, ăng-ten định hướng độ lợi cao và lớp phủ tàng hình đa hướng trong kỹ thuật vô tuyến.

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

Khung lý thuyết áp dụng

Luận văn xây dựng trên nền tảng điện từ học cổ điển kết hợp với lý thuyết môi trường chiết suất âm kép (Left-Handed Metamaterials). Theo định luật Veselago và Pendry, khi độ điện thẩm hiệu dụng và độ từ thẩm hiệu dụng đồng thời nhận giá trị âm ($\epsilon < 0, \mu < 0$), ba véc-tơ điện trường $\mathbf{E}$, từ trường $\mathbf{H}$ và véc-tơ sóng truyền $\mathbf{k}$ tạo thành một tam diện nghịch tuân theo quy tắc bàn tay trái. Khi đó, chiều truyền pha của sóng điện từ sẽ ngược chiều với chiều truyền năng lượng của véc-tơ Poynting, dẫn tới hiện tượng khúc xạ âm, dịch chuyển Doppler nghịch đảo và bức xạ Cherenkov ngược hướng.

Để giải quyết bài toán mở rộng dải tần, tác giả áp dụng mô hình lai hóa plasmon bậc một và bậc hai. Trong mô hình lai hóa bậc một của cấu trúc cộng hưởng từ, sự tương tác điện từ giữa hai phần tử kim loại dẫn đến sự tách mức năng lượng thành mode đối xứng (cộng hưởng điện, tần số cao) và mode bất đối xứng (cộng hưởng từ, tần số thấp). Khi nâng cấp lên mô hình lai hóa bậc hai với hệ hai lớp cấu trúc đặt liền kề dọc theo phương truyền sóng, hiện tượng hỗ cảm mạnh mẽ giữa hai lớp làm suy biến và tách mode cộng hưởng từ ban đầu thành hai mode mới riêng biệt, từ đó mở rộng đáng kể vùng tần số có độ từ thẩm âm.

Bên cạnh đó, mô hình mạch điện LC tương đương được áp dụng để giải thích định lượng cơ chế cộng hưởng. Độ tự cảm $L$ và điện dung $C$ của đĩa kim loại và khe hở điện môi được xác định chính xác để dự đoán tần số cộng hưởng từ $f_m$ và tần số cộng hưởng điện $f_e$. Năm khái niệm cốt lõi chi phối toàn bộ nghiên cứu bao gồm: vật liệu biến hóa chiết suất âm kép, độ từ thẩm âm nhân tạo, tương tác hỗ cảm điện từ, mô hình lai hóa plasmon và tần số plasma hiệu dụng của lưới dây kim loại.

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

Nghiên cứu kết hợp chặt chẽ giữa tính toán giải tích mạch tương đương, mô phỏng trường điện từ 3D và thuật toán trích xuất tham số hiệu dụng. Cỡ mẫu nghiên cứu bao gồm các ô cơ sở tuần hoàn được khảo sát toàn diện với hằng số mạng $a_x = 8,0\text{ mm}$, $a_y = 7,5\text{ mm}$, bán kính đĩa kim loại $R = 3,0\text{ mm}$, bề dày lớp đồng $t_m = 0,036\text{ mm}$ với độ dẫn điện $\sigma = 5,88 \times 10^7\text{ S/m}$, và lớp đệm điện môi FR4 có độ dày $t_d = 0,4\text{ mm}$ với hằng số điện môi $\epsilon_r = 4,3$, hệ số tổn hao $\tan\delta = 0,025$. Khoảng cách giữa hai lớp lưới $d$ được biến thiên từ 0,2 mm đến 1,6 mm để khảo sát tương tác ngoài.

Phương pháp chọn mẫu tham số được thực hiện theo nguyên tắc quét từng biến độc lập nhằm cô lập và đánh giá chính xác tác động của từng kích thước hình học đến vị trí cộng hưởng. Để mô phỏng phản xạ $S_{11}$ và truyền qua $S_{21}$, tác giả sử dụng phần mềm CST Microwave Studio phiên bản 2017 với phương pháp tích phân hữu hạn (FIT) trong miền tần số, thiết lập điều kiện biên tuần hoàn. Lý do lựa chọn FIT miền tần số là vì phương pháp này mang lại độ chính xác cao đối với các cấu trúc cộng hưởng tuần hoàn có kích thước nhỏ hơn nhiều so với bước sóng công tác. Dữ liệu tán xạ $S$ sau đó được đưa vào thuật toán trích xuất của Chen phát triển từ phương pháp Nicolson-Ross-Weir trên nền tảng MATLAB nhằm loại bỏ tính đa trị của hàm số phức và xác định chính xác phần thực và ảo của $\epsilon, \mu, n, z$ trong thời gian nghiên cứu 2017-2018.

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

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

Nghiên cứu đã đạt được 4 phát hiện khoa học quan trọng có giá trị thực tiễn cao:

Thứ nhất, cấu trúc cặp đĩa hai lớp đã chứng minh khả năng tách vạch cộng hưởng từ xuất sắc theo đúng dự đoán của mô hình lai hóa bậc hai. Trên phổ truyền qua xuất hiện hai cực tiểu suy giảm rõ rệt tại tần số 13,2 GHz và 14,5 GHz, tạo ra vùng có độ từ thẩm âm liên tục trải rộng từ 12,8 GHz đến 15,0 GHz, mở rộng hơn 1,4 GHz (tương đương mức tăng băng thông khoảng 180%) so với cấu trúc đơn lớp.

Thứ hai, khoảng cách giữa hai lớp cấu trúc $d$ đóng vai trò quyết định cường độ tương tác ngoại. Khi giảm $d$ từ 1,6 mm xuống 0,2 mm, hệ số ghép hỗ cảm giữa hai lớp tăng mạnh, làm khoảng cách giữa hai mode cộng hưởng tách rời tăng thêm khoảng 35%, đồng thời nới rộng đáng kể dải tần từ thẩm âm.

Thứ ba, chiều dày lớp điện môi $t_d$ kiểm soát tương tác nội của từng cặp đĩa. Khi $t_d$ tăng từ 0,2 mm lên 1,0 mm, tần số cộng hưởng từ dịch chuyển dần về phía vùng tần số cao theo quy luật tỷ lệ với căn bậc hai của $t_d$, phù hợp hoàn toàn với công thức tính toán từ mạch LC lý thuyết.

Thứ tư, khi tích hợp cấu trúc cặp đĩa hai lớp với lưới dây kim loại liên tục có độ rộng $w = 1,0\text{ mm}$, vật liệu đạt được vùng chiết suất âm kép ($n < 0$) đồng thời với $\mu < 0$ và $\epsilon < 0$ trong dải tần 13,0 GHz đến 15,2 GHz, duy trì độ truyền qua cao đạt mức trên -2 dB. Đặc biệt, tính chất chiết suất âm hoàn toàn ổn định khi góc phân cực sóng điện từ quay từ 0 độ đến 30 độ.

Thảo luận kết quả

Nguyên nhân vật lý sâu xa của hiện tượng mở rộng dải tần nằm ở sự phân bố dòng điện cảm ứng bề mặt trên các đĩa kim loại. Tại mode tần số thấp (13,2 GHz), dòng điện trên hai mặt đĩa trong mỗi lớp chạy ngược chiều nhau tạo thành các vòng xoáy dòng điện khép kín đối song, sinh ra mô-men từ lưỡng cực mạnh chống lại từ trường ngoài. Ở mode tần số cao (14,5 GHz), sự tương tác hỗ cảm giữa lớp trước và lớp sau tạo ra trạng thái ghép cặp năng lượng mới, duy trì phản ứng nghịch từ mạnh mẽ trên một khoảng tần số liên tục thay vì chỉ xuất hiện tại một điểm cộng hưởng đơn lẻ như các cấu trúc vòng cộng hưởng hở (SRR) hay cặp dây bị cắt (CWP) truyền thống.

Dữ liệu nghiên cứu được biểu diễn trực quan qua các đồ thị phổ truyền qua $S_{21}$ và các đường cong phần thực của chiết suất $\text{Re}(n)$, độ từ thẩm $\text{Re}(\mu)$, độ điện thẩm $\text{Re}(\epsilon)$ theo trục tần số 12–18 GHz. Bảng số liệu tổng hợp các tham số tối ưu $a_x = a_y = 8,0\text{ mm}, R = 3,0\text{ mm}, t_d = 0,4\text{ mm}, d = 0,5\text{ mm}, w = 1,0\text{ mm}$ minh chứng rõ ràng sự vượt trội về độ suy hao thấp và tính đẳng hướng quang học. So với các công bố trước đây về cấu trúc CWP vốn chỉ hoạt động ở một góc phân cực cố định, việc chuyển sang dạng hình học đĩa tròn đã loại bỏ hoàn toàn tính dị hướng, cho phép sóng điện từ chiếu tới ở nhiều góc phân cực khác nhau mà không làm suy giảm hiệu ứng chiết suất âm.

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

Dựa trên các kết quả đạt được, luận văn đưa ra 4 khuyến nghị và giải pháp hành động cụ thể cho công tác nghiên cứu và triển khai ứng dụng:

  1. Tối ưu hóa đa tham số hình học: Các nhóm nghiên cứu vật liệu biến hóa cần ứng dụng thuật toán di truyền kết hợp mô phỏng CST để tự động tối ưu tỷ số khoảng cách $d/t_d$ trong ngưỡng 0,5 đến 0,75, nhằm mục tiêu gia tăng băng thông chiết suất âm lên trên 2,5 GHz trong lộ trình nghiên cứu 6 tháng tới.
  2. Thay thế vật liệu nền điện môi cao cấp: Các phòng thí nghiệm vật lý vật liệu cần thay thế chất nền FR4 có hệ số tổn hao cao ($\tan\delta = 0,025$) bằng các loại vật liệu điện môi vi sóng chuyên dụng như Rogers RT/duroid 5880 hoặc Teflon ($\tan\delta \le 0,0009$) nhằm giảm tổn hao truyền qua xuống dưới mức 0,5 dB trong thời gian 12 tháng.
  3. Chế tạo thực nghiệm và đo kiểm vi ba: Nhóm nghiên cứu tại các trường đại học kỹ thuật cần phối hợp với các xưởng vi cơ điện tử để chế tạo tấm mạch in PCB đa lớp kích thước $200 \times 200\text{ mm}$, sau đó tiến hành đo đạc thực tế thông số tán xạ trong buồng câm vi ba bằng máy phân tích mạng véc-tơ (VNA) trong vòng 9 tháng để kiểm chứng độ chính xác so với mô phỏng.
  4. Mở rộng dải tần hoạt động lên vùng Terahertz và Quang học: Viện nghiên cứu công nghệ bán dẫn và quang điện tử cần mở rộng áp dụng mô hình cấu trúc đĩa đa lớp xuống kích thước nanômét (bán kính đĩa $R < 500\text{ nm}$) bằng kỹ thuật quang khắc chùm tia điện tử, hướng tới mục tiêu hiện thực hóa siêu thấu kính quang học khả kiến trong khung thời gian 18 đến 24 tháng.

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

Công trình luận văn này là tài liệu tham khảo chuyên môn giá trị cho 4 nhóm đối tượng trọng tâm:

  • Học viên cao học và nghiên cứu sinh ngành Vật lý Quang học, Vật lý Vô tuyến: Tài liệu cung cấp hệ thống cơ sở lý thuyết chuẩn xác về mô hình lai hóa plasmon, cách thiết lập mạch LC tương đương và bộ mã nguồn MATLAB trích xuất tham số hiệu dụng theo thuật toán Chen để phục vụ trực tiếp cho các đề tài nghiên cứu chuyên sâu.
  • Kỹ sư thiết kế ăng-ten và thiết bị siêu cao tần viễn thông: Vận dụng cấu trúc lưới đĩa hai lớp để thiết kế thấu kính hội tụ phẳng, bộ phản xạ bề mặt chọn lọc tần số (FSS) và ăng-ten mảng có băng thông rộng, tăng cường độ định hướng cho các trạm thu phát sóng di động thế hệ mới 5G và 6G.
  • Chuyên gia công nghệ quốc phòng và kỹ thuật Radar: Tham khảo các đặc trưng phân bố điện từ trường và tính chất đẳng hướng không phụ thuộc phân cực trong dải 12–18 GHz để chế tạo lớp phủ hấp thụ sóng radar, cấu trúc tàng hình điện từ và nâng cao độ phân giải của hệ thống trinh sát vô tuyến.
  • Giảng viên và nhà nghiên cứu tại các viện, trường đại học kỹ thuật: Sử dụng nội dung luận văn làm tài liệu giảng dạy chuyên đề vật liệu biến hóa, hướng dẫn thực hành mô phỏng trên phần mềm CST Studio Suite và xây dựng các bài thí nghiệm quang điện tử hiện đại.

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

  1. Tại sao cấu trúc cặp đĩa lại vượt trội hơn cấu trúc cặp dây kim loại bị cắt (CWP)? Cấu trúc cặp đĩa sở hữu tính đối xứng tròn bậc cao trong mặt phẳng, giúp các đặc tính phản xạ và truyền qua hoàn toàn không thay đổi khi quay góc phân cực từ 0 độ đến 30 độ. Trong khi đó, cấu trúc CWP dạng thanh chữ nhật bị phụ thuộc nghiêm trọng vào phương phân cực của véc-tơ điện trường, gây hạn chế lớn khi ứng dụng thực tế trong không gian mở.

  2. Mô hình lai hóa bậc hai giúp mở rộng dải tần từ thẩm âm bằng cơ chế nào? Khi đặt hai lớp đĩa kim loại liền kề dọc theo phương truyền sóng, hiện tượng hỗ cảm điện từ mạnh mẽ giữa hai lớp làm suy biến mode cộng hưởng từ đơn lẻ ban đầu thành hai mode mới riêng biệt tại 13,2 GHz và 14,5 GHz. Sự kết hợp của hai mode cộng hưởng kế tiếp nhau này giúp nới rộng dải tần từ thẩm âm liên tục thêm hơn 1,4 GHz.

  3. Vì sao phải kết hợp cấu trúc đĩa với lưới dây kim loại liên tục để tạo chiết suất âm? Cấu trúc đĩa hai lớp chỉ tạo ra độ từ thẩm âm ($\mu < 0$), trong khi để đạt được chiết suất âm kép ($n < 0$) cần đồng thời cả độ điện thẩm âm ($\epsilon < 0$). Lưới dây kim loại liên tục đóng vai trò như một môi trường plasma nhân tạo có tần số plasma cao, cung cấp $\epsilon < 0$ trên toàn bộ dải tần 12–18 GHz.

  4. Thuật toán của Chen đóng vai trò gì trong việc xử lý kết quả mô phỏng? Thuật toán của Chen cho phép trích xuất chính xác các tham số vật lý hiệu dụng gồm $\epsilon, \mu, n$ và trở kháng $z$ từ dữ liệu phản xạ $S_{11}$ và truyền qua $S_{21}$. Thuật toán giải quyết triệt để bài toán chọn nhánh nghiệm phức đối với môi trường thụ động, đảm bảo phần thực của trở kháng luôn dương và phần ảo của chiết suất thỏa mãn tính tắt dần năng lượng.

  5. Vật liệu biến hóa trong luận văn có thể chế tạo bằng công nghệ hiện có tại Việt Nam không? Hoàn toàn khả thi vì luận văn sử dụng các vật liệu phổ thông và giá thành rẻ như lá đồng có độ dày chuẩn 0,036 mm và tấm mạch in FR4 dày 0,4 mm. Các kích thước hình học ở mức milimét ($R = 3,0\text{ mm}, a = 8,0\text{ mm}$) hoàn toàn phù hợp với công nghệ quang khắc và gia công mạch in PCB tiêu chuẩn hiện nay.

Kết luận

  • Thiết kế và tối ưu thành công cấu trúc vật liệu biến hóa đĩa hai lớp và lưới đĩa hai lớp hoạt động hiệu quả trong dải tần số sóng Radar băng Ku từ 12,0 GHz đến 18,0 GHz.
  • Vận dụng thành công mô hình lai hóa plasmon bậc hai và lý thuyết mạch LC tương đương để giải thích tường minh cơ chế ghép hỗ cảm và tách vạch cộng hưởng từ.
  • Mở rộng thành công dải tần có độ từ thẩm âm ($\mu < 0$) thêm hơn 1,4 GHz và đạt được vùng chiết suất âm kép ($n < 0$) liên tục từ 13,0 GHz đến 15,2 GHz với độ truyền qua cao trên -2 dB.
  • Khắc phục triệt để nhược điểm phụ thuộc vào góc phân cực của sóng điện từ, duy trì tính đẳng hướng hoàn hảo từ 0 độ đến 30 độ nhờ hình học đối xứng tròn.
  • Xây dựng hoàn chỉnh quy trình tính toán, mô phỏng trên phần mềm CST Studio 2017 và trích xuất tham số hiệu dụng bằng thuật toán Chen trên MATLAB.

Công trình của tác giả Nguyễn Thị Mây dưới sự hướng dẫn của TS. Nguyễn Thị Hiền tại Trường Đại học Khoa học - Đại học Thái Nguyên là một đóng góp học thuật xuất sắc vào lĩnh vực nghiên cứu vật liệu biến hóa tại Việt Nam. Trong giai đoạn 6 đến 12 tháng tới, các nhóm nghiên cứu nên tiếp tục phát triển chế tạo mẫu thực nghiệm và mở rộng cấu trúc sang dải tần Terahertz để ứng dụng vào công nghệ viễn thông siêu cao tần. Độc giả quan tâm hãy kết nối với các nhóm nghiên cứu quang học và vật lý vô tuyến để đẩy mạnh hợp tác triển khai ứng dụng công nghệ vật liệu biến hóa tiên tiến vào thực tiễn.