결합 공명 그래핀 메타물질 흡수체: 튜너블 중적외선 변조

초록

중적외선 영역에서 그래핀이 적재된 메타물질 흡수체를 연구했습니다. 십자형 슬롯을 통한 결합 공명에 의해 빛-그래핀 상호작용이 크게 향상되었습니다. 흡수 피크는 페르미 준위가 증가함에 따라 현저한 청색 편이를 보이며, 흡수체의 광범위한 튜너빌리티를 가능하게 합니다. 단순한 회로 모델이 이러한 변조 거동을 잘 설명하고 예측합니다. 본 제안은 스위칭, 센싱, 변조, 생화학 검출 등 다양한 분야에 응용될 수 있습니다.

배경

플라즈모닉 메타물질(PM) 흡수체는 깊은 서브파장 스케일의 금속 나노구조로 작동합니다. 완벽한 흡수가 특정 파장에서 달성되고 조절될 수 있어, 발광체/검출기, 센서, 광열 치료, 광역학적 상호작용, 초분광 이미징 등 다양한 응용으로 이어집니다[1-7]. PM 흡수체는 또한 조절 가능한 특성을 가진 새로운 기능성 소자 설계에 유망한 플랫폼을 제공합니다. 액정, 반도체, 상변화 물질 등의 구성 요소를 도입함으로써 광학 응답을 전기적, 광학적, 열적으로 변조할 수 있습니다[8-13], 이는 새로운 유형의 변조기, 스위치, 다중 스펙트럼 검출기를 가능하게 합니다.

최근 그래핀은 플라즈모닉 물질로서의 고속 변조 능력과 튜너빌리티로 상당한 주목을 받고 있습니다[14-20]. 구체적으로 그래핀 전도도는 페르미 준위(E_F)에 의존하며, 이는 바이어스 전압을 통해 수 나노초 내에 연속적으로 조절 가능하여 근적외선 및 중적외선 영역에서 높은 변조율을 구현합니다[17, 19-24]. 그러나 단일 그래핀 층이 원자 수준의 두께에 불과하기 때문에 입사광과 플라즈모닉 공명 간의 상호작용이 매우 약합니다. 또한 파울리 차단으로 인해 중적외선 영역에서 이 상호작용은 더욱 약해집니다[22]. 결과적으로 파장 튜닝 범위와 변조 깊이가 매우 제한적입니다. 파장 이동은 일반적으로 공명 파장의 10% 미만으로[21, 22, 25-28], 광통신 및 광대역 스펙트럼 검출의 실용적 응용에는 여전히 과제로 남아 있습니다. 따라서 효율적인 전기광학 변조를 위해서는 그래핀-빛 상호작용을 크게 강화할 필요가 있습니다. 이전 연구들에서 나노 안테나 및 분할 고리 공진기 등 복잡한 나노구조 설계를 기반으로 그래핀-빛 상호작용 향상이 이론적, 실험적으로 입증되었습니다[19, 21, 22, 25, 27, 28]. 그러나 이러한 설계는 일반적으로 복잡하거나 편광 의존적이며, 작동 주파수 범위가 상대적으로 작고 튜너빌리티가 여전히 제한적입니다.

본 연구에서는 9~14 μm 범위의 변조가 가능한 그래핀 적재 흡수체를 제안했습니다. 이는 생화학 센싱 및 열화상 등 응용 분야에서 큰 관심을 끌고 있습니다[5, 29-31]. 십자형 슬롯 내부의 결합 공명은 전기장을 4개 주문만큼 향상시켜 그래핀-빛 상호작용을 강력하게 강화하며, 중심 파장에서 최대 25%의 이동을 초래합니다. 또한 전압과 기하학적 파라미터에 의해 제어되는 그래핀 유도 변조를 잘 설명하고 예측하는 단순한 LC 회로 모델을 제안합니다. 이러한 광범위한 튜너빌리티는 많은 응용 분야에서 유망할 것입니다.

방법

그림 1a와 같이, 주기 Λ = 8 μm로 패턴화된 금속 패치가 유전체 스페이서로 분리된 금속 기판 위에 배열되어 있습니다. 단일 그래핀 층이 패치와 스페이서 사이에 샌드위치 구조로 위치합니다. 기판은 매우 두꺼워 반사 거울 역할을 합니다. 스페이서 층의 두께는 t_d = 520 nm, 금속 패치의 두께는 t_m = 100 nm입니다. 그림 1b는 단위 셀의 평면도입니다. 두 개의 서브유닛이 편광 독립성을 지원하기 위해 대각선 대칭으로 배열되어 있습니다. 십자형 슬롯이 각 정사각형 패치에 식각되어 네 개의 동일한 작은 조각으로 나눕니다. S_1과 S_2의 작은 조각 크기는 각각 l_1 = 1.5 μm, l_2 = 1.7 μm입니다. 두 서브유닛의 슬롯 폭은 a = 20 nm입니다. 금속 재료는 금(Au)으로 선택되었으며, 그 광학적 특성은 ω_p = 1.369 × 10^16 Hz, τ = 1.224 × 10^14 Hz인 드루드 모델 ε(ω) = 1 - ω_p^2/(ω(ω + τ))로 기술됩니다[32]. 유전체 스페이서는 황화아연(ZnS)으로 구성되며, 중적외선 영역에서 무시할 수 있는 손실로 광학 지수 n = 2.2를 가집니다[33].

유한 차분 시간 영역(FDTD; Lumerical FDTD Solutions) 방법을 사용하여 반사율 스펙트럼과 전자기장 분포를 계산했습니다. 시뮬레이션은 x, y 방향으로 주기적 경계 조건, z 방향으로 완벽 매칭 레이어 조건으로 수행되었습니다. 단일 그래핀 층은 표면 전도도 접근법으로 2차원 구조로 모델링되었습니다[34]. 그래핀 층의 표면 전도도 σ_g는 쿠보 공식으로 계산되며, 이는 밴드간 항 σ_inter와 밴드내 항 σ_intra를 포함합니다[35].

여기서 e와 ξ는 전자의 전하와 에너지, ℏ는 디랙 상수, ω는 각진동수, f_d ≡ 1/(e^{(ξ-E_F)/k_BT}+1)는 페르미-디랙 분포, T는 절대 온도, Γ은 산란율, k_B는 볼츠만 상수, E_F는 페르미 준위입니다. 계산에서 T = 300 K, Γ = 10 meV로 설정했습니다[28]. 그래핀 층 근처의 메시 크기는 0.25 nm, 슬롯에서는 2.5 nm입니다. 그래핀의 유효 유전율은 다음과 같이 표현됩니다.

ε_g = 1 + iσ_g/(ε_0 ω t_g)

여기서 ε_0는 진공 유전율, t_g는 그래핀 층 두께입니다. 식 (1)과 (2)는 그래핀의 광학 상수가 E_F에 따라 변함을 보여줍니다. 이 변화는 흡수 주파수의 튜너빌리티를 가져오며, 나노구조의 결합 공명에 의해 그 범위가 크게 확대되어 소자에서 요구되는 전압을 실질적으로 낮춥니다.

결과 및 논의

그림 2a는 정상 입사에서 x-편광파(φ = 0)에 대한 흡수 스펙트럼을 보여줍니다. 페르미 준위 E_F = 0 eV일 때, 파장 λ = 12.4 μm와 13.3 μm에서 두 개의 흡수 피크가 관측됩니다. 12.1~13.5 μm 범위의 입사광은 나노구조에 의해 거의 흡수됩니다. E_F가 증가함에 따라 공명은 짧은 파장 쪽으로 이동합니다. E_F = 0.2 eV에서 흡수 피크는 11.8 μm와 12.46 μm로 이동하며, 각각 4.8%와 6%의 상대적 이동을 보입니다. 한편 피크 2의 흡수율은 감소하는데, 이는 높은 E_F에서 메타물질과 공기 사이의 임피던스 불일치에 기인합니다[28]. 흥미롭게도 페르미 준위가 계속 증가함에 따라 피크 2가 피크 1보다 더 빠르게 청색 편이합니다. 이 관찰된 거동은 나중에 회로 모델로 설명하겠습니다.

변조는 매개변수 M = Δλ/λ_0로 정량화할 수 있습니다. 여기서 λ_0는 E_F = 0 eV에서의 공명 파장, Δλ는 E_F 변화에 의한 파장 이동입니다. 그림 2a에서 E_F가 0.6 eV에 도달할 때 피크 1과 피크 2에 대해 각각 M_1 = 20.1%, M_2 = 25.5%를 보입니다. 공명의 변조 범위는 이전 연구들[19, 21, 22, 25-28]보다 훨씬 넓습니다. 낮은 E_F에서 이러한 큰 변조는 많은 응용 분야에서 매우 바람직합니다. 별도 계산에 따르면 흡수 피크는 스페이서 두께가 감소함에 따라 청색 편이합니다(추가 파일 1). 따라서 변조의 적절한 시작점을 설정하기 위해 두께를 최적화할 수 있습니다. 또한 제안된 메타물질의 광학 응답은 그림 2b와 같이 편광 독립적입니다. 설계의 대칭성 덕분에 편광 각도 φ가 0°에서 90°로 변해도 흡수 스펙트럼은 변하지 않습니다.

완벽한 흡수 메커니즘은 공명에서의 장 분포로 명확히 설명됩니다. 그림 1에 표시된 잘 알려진 금속-절연체-금속(MIM) 구조[3, 32, 36-38]로 인해 국소 표면 플라즈몬 폴라리톤(SPP)이 자극되어 각 패치에 콤팩트한 자기 공명을 형성합니다. 그림 3a, b는 E_F = 0.2 eV에서 공명 파장 λ_1 = 11.8 μm와 λ_2 = 12.46 μm에서 그래핀 층의 정규화된 자기장 |H|^2를 보여줍니다. SPP가 강하게 국소화되어 두 서브유닛이 독립적으로 작동할 수 있습니다. 그러나 각 서브유닛 내부의 분할 슬롯 폭이 좁기 때문에 네 조각의 공명이 실제로 서로 결합되어 있습니다. 이 결합은 슬롯 내부의 전기장을 엄청나게 증가시킵니다(그림 3c, d). 입사광이 x 편광이므로 y 방향 슬롯의 E 필드만 여기서 명백합니다. 공명 결합에 의해 향상된 E 필드 강도는 입사광 E_inc보다 4개 주문 더 큽니다. 반면 이전 연구에서 변조에 사용된 가장 강화된 필드는 패치 가장자리에 있었습니다. 그림 3e, f는 그림 3c, d의 흰색 선을 따라 슬롯과 가장자리에서의 향상 정도를 날카롭게 비교합니다.

이러한 장 분포는 우리 제안에서 변조가 왜 그렇게 큰지 잘 설명합니다. 섭동 이론에 기반하여, 그래핀 유도 공명 이동은 Δω = -iσ_g ∫_S |E_s|^2 dS / W_0로 평가됩니다[22]. 여기서 |E_s|^2는 그래핀 층의 전기장 강도, W_0는 저장 에너지, S는 그래핀이 덮인 영역입니다. 공명의 스펙트럼 이동(Re(Δω))은 σ_g의 허수부에 의해 결정되며, 이는 중적외선 영역에서 실수부보다 훨씬 큽니다[22, 28]. 그림 3c-f에서 명확히 보듯이, 좁은 슬롯 내 전기장 향상은 가장자리에서의 향상보다 10배 이상 큽니다. 결과적으로 적분 값은 패치 슬롯에서 크게 향상된 E 필드에 의해 주로 기여되어, 금속 가장자리에서만 향상된 E 필드를 가진 이전 경우들[21, 22, 25, 27, 28]보다 훨씬 큰 피크 이동을 초래합니다.

장 분포와 위 논의를 바탕으로 튜닝 거동을 연구하기 위한 LC 회로 모델을 제안합니다. 그림 4a와 같이, L_i와 C_i(i = 1, 2)는 그림 1b의 패치 S_i에 대한 인덕턴스와 커패시턴스입니다. 슬롯 폭 a가 매우 크고 그래핀 층이 없을 때, 슬롯과 그래핀에 의한 효과는 무시할 수 있습니다. 그러면 L_i와 C_i는 흡수 스펙트럼에서 얻은 공명 파장으로 피팅하여 별도 계산으로 결정할 수 있습니다[37, 39, 40]. 결과는 서브유닛 S_1에 대해 L_1 = 0.07 pH, C_1 = 350 aF, 서브유닛 S_2에 대해 L_2 = 0.075 pH, C_2 = 380 aF입니다. 각 서브유닛 내부의 슬롯 유도 결합 효과는 션트 커패시턴스 C_c로 설명되며, 이는 슬롯 폭 a가 증가함에 따라 감소하는 것으로 밝혀졌습니다. 우리 경우 a = 20 nm에서 C_c = 290 aF이며, a가 10 nm 증가할 때마다 200 aF, 180 aF, 135 aF가 됩니다. 공명 파장은 회로 임피던스를 0으로 하여 얻습니다: λ_i^0 = 2πc_0√(L_i C_i^0). 여기서 c_0는 진공 중 빛 속도, "i"는 서브유닛 S_i를 나타내며, C_i^0 = C_i + C_c입니다.

2차원 그래핀 층은 기본적으로 인덕터 역할을 합니다. 그림 3에서 보듯 그래핀 층의 주요 기여는 전기장이 강화된 슬롯 위치에서 옵니다. 슬롯 폭이 동작 파장과 그래핀 플라즈몬 파장보다 훨씬 작으므로 준정적 근사가 유효합니다. 슬롯 양단의 전압 V와 전류 I는 V = aE, I = 2l_i t_g(σ_g - iωε_0)E로 평가할 수 있습니다. 여기서 E는 그래핀 층의 전기장입니다. 따라서 그래핀 층의 기여를 설명하는 인덕턴스 L_g = -1/ω Im(V/I)[41]를 도입할 수 있으며, 이는 다음과 같이 표현됩니다.

L_g = a / (2l_i ω^2 ε_0 |Re(ε_g)| t_g) (i = 1, 2)

이 인덕터는 그림 4a와 같이 병렬 요소로 작용합니다. 결과적으로 한 패치의 총 인덕턴스는 1/L_i' = 1/L_i + 1/L_g로 얻습니다. 그래핀 층이 있는 각 서브유닛의 최종 공명 파장은 다음과 같습니다.

λ_i' = 2πc_0√(L_i' C_i^0) (i = 1, 2)

각 서브유닛이 독립적으로 작동하므로 메타물질의 총 임피던스는 두 서브유닛 임피던스의 병렬 연결로 얻을 수 있습니다.

이 LC 모델은 E_F 증가에 따른 공명의 청색 편이를 예측합니다. 식 (1)과 (2)에서 유도하면, 더 높은 E_F에서 그래핀의 |Re(ε_g)| 값이 더 커지며, 이는 식 (3)에서 더 작은 L_g를 제공합니다. 인덕터의 병렬 연결로 인해 최종 인덕턴스 L_i'가 더 작아지며, 식 (4)에서 더 짧은 공명 파장을 유도합니다. 계산 결과는 그림 4b에 요약되어 있으며, FDTD 시뮬레이션으로 얻은 공명 파장과 잘 일치합니다. 작은 편차는 우리 LC 모델이 각 패치 가장자리의 약한 필드 기여를 무시하기 때문입니다(그림 3c-f). LC 모델은 또한 기하학적 파라미터가 공명의 청색 편이에 어떻게 영향을 미치는지 보여줍니다. 식 (4)를 미분하면 ∂λ_i'/∂L_i' ∝ 1/√L_i'입니다. 작은 √L_i' 값이 이 청색 편이의 민감도를 높이는 데 유리함을 알 수 있습니다. 인덕터가 병렬로 연결되고 L_i가 고정되어 있으므로, 작은 총 인덕턴스 L_i'는 작은 그래핀 인덕턴스 L_g를 의미합니다. 식 (3)에 따르면 튜닝 범위를 늘리기 위해 슬롯 폭 a는 작아야 하고 패치 크기 l은 커야 합니다. 그림 4c는 S_1 내부 슬롯 폭이 50에서 20 nm로 감소할 때 E_F = 0.4 eV에서 공명의 청색 편이가 약 6%에서 15%로 증가함을 보여줍니다. 반면 슬롯 폭을 a = 20 nm로 고정하면, 패치 크기가 1.5에서 1.8 μm로 변할 때 공명이 15%에서 22%로 증가합니다(그림 4d). FDTD 시뮬레이션과의 좋은 일치는 이러한 단순한 회로 모델이 관련 메타물질 소자 연구에 효율적인 방법임을 입증합니다.

결론

결론적으로, 우리는 큰 변조 범위를 가진 편광 독립적 광대역 메타물질 흡수체를 설계했습니다. 두 공명 모두 E_F가 0에서 0.6 eV로 증가할 때 중심 파장의 20.1%와 25.5%까지 튜닝 범위에 도달합니다. 이러한 큰 변조는 각 금속 패치의 십자형 슬롯 내 결합 공명에 의해 엄청나게 향상된 그래핀-빛 상호작용에서 비롯됩니다. 이 효과는 LC 모델에서 그래핀 도입 인덕터로 잘 설명됩니다. 이러한 단순한 모델은 다양한 기하학적 파라미터 하에서 변조 거동을 예측하며, 결과가 FDTD 시뮬레이션과 잘 일치합니다. 본 제안은 광통신, 센싱, 열화상 등 잠재적 응용에 유익합니다.

약어

E_F:

페르미 준위

FDTD:

유한 차분 시간 영역

MIM:

금속-절연체-금속

PM:

플라즈모닉 메타물질

ZnS:

황화아연