본 프로토콜은 6G 무선 통신 응용을 위해 테라헤르츠 대역(0.5-30 THz)에서 작동하는 단일벽 탄소 나노튜브 기반 광학 재구성 가능 지능형 표면의 정확한 전자기 모델링을 위한 하이브리드 인공 신경망 및 Z-변환 방법을 제시하며, 180× 반사 위상 튜닝 가능 범위가 XX를 초과하는 계산 가속화 310° 및 빔 조향 범위는 ±45°.
방법 논문
테라헤르츠 빔 스티어링을 위한 탄소 나노튜브 기반 재구성 가능한 지능형 표면 모델링을 위한 하이브리드 ANN-Z 방법
이 논문에서
요약
초록
단일벽 탄소 나노튜브 기반의 광학 재구성 가능 지능형 표면은 미래 6G 무선 시스템의 테라헤르츠 빔 조향 및 광파 제어를 위한 유망한 솔루션을 제공합니다. 그러나 이러한 구조의 정확한 모델링은 양자 수송 효과, 운동 유도성 및 넓은 주파수 범위에 걸친 다중 공명 엑시톤 동작으로 인해 여전히 어려운 과제로 남아 있습니다. 본 프로토콜은 단일벽 탄소 나노튜브 광학 재구성 가능 지능형 표면 단위 셀의 정확한 특성 분석을 위해 Kubo 정식을 이용한 양자 전도도 모델링, 다항 회귀 기반 데이터 평활화 및 Z-도메인 전달 함수 분석을 통합한 하이브리드 계산 프레임워크를 설명합니다. 이 방법은 (10,5) 카이랄성(직경 0.60 nm, 밴드갭 1.762 eV)을 가진 석영 기판 위의 교차 단일벽 탄소 나노튜브 나노 스트립 공진기 설계로 시작하며, 이어서 0.5-30 THz 대역에서 CST Microwave Studio를 이용한 전파 전자기 시뮬레이션을 수행합니다. 8차 다항 회귀 모델을 통해 추출된 S-파라미터를 처리하여 수치적 변동을 제거하고 평활화된 전자기 응답을 예측합니다. 이후 QR 분해를 이용한 최소 제곱 최적화를 통해 분자 차수 6, 분모 차수 7의 이산 전달 함수 H(z)를 피팅하여 극점-영점 안정성 분석 및 수동성 검증을 가능하게 합니다. 또한 본 프로토콜은 빔 조향 최적화를 위해 화학적 포텐셜 변조를 통한 양자 전도도 튜닝을 포함합니다. 대표적인 결과에 따르면 310°를 초과하는 반사 위상 가변성, 최대 92.3%의 흡수 향상, 사이드 로브 레벨 -12 dB 이하의 ±45° 빔 조향 범위, 그리고 기존의 전파 최적화 방법 대비 180×의 계산 가속화를 입증하였습니다. 다항 회귀는 R2 계수 0.994와 함께 0.0688의 테스트 평균 제곱근 오차를 달성하였으며, H(z) 피팅은 0.89 dB의 평균 제곱근 오차를 달성하였습니다. 안정성 분석을 통해 모든 극점이 단위 원 내에 있음을 확인하였습니다. 본 프로토콜은 6G 및 그 이후의 프로그래밍 가능한 광학 메타표면 및 지능형 테라헤르츠 통신 시스템 설계를 위한 효율적이고 재현 가능한 경로를 제공합니다.
서론
6세대(6G) 무선 시스템으로의 급격한 진화는 1 Tbps를 초과하는 초고속 데이터 전송 속도, 지능형 센싱, 홀로그래픽 빔포밍 및 적응형 파면 엔지니어링을 구현하기 위해 테라헤르츠(THz) 및 광주파수 대역에 대한 탐구를 가속화했습니다1,2,3. 테라헤르츠 대역(0.1-30 THz)은 풍부한 대역폭을 제공하지만, 심각한 자유 공간 감쇠(1 THz에서 약 20 dB·km-1), 557 GHz, 752 GHz, 988 GHz 및 1.13 THz에서의 수증기로 인한 대기 분자 흡수, 그리고 대기 입자 및 강우에 의한 차단에 매우 민감하다는 단점이 있습니다4.
재구성 가능한 지능형 표면(Reconfigurable intelligent surfaces)은 프로그래밍 가능한 반사, 흡수, 위상 제어 및 파동 집중을 통해 전자기 파면을 동적으로 조작할 수 있는 지원 기술로 등장하였습니다5,6,7. 이러한 표면을 광학 주파수 영역으로 확장하면 빛의 전파를 전례 없는 수준으로 제어할 수 있어, LiDAR, 자유 공간 광통신, 홀로그래픽 디스플레이 및 적응형 빔 스티어링과 같은 응용 분야를 가능하게 합니다8.
금속 공진기에 기반한 기존의 재구성 가능한 지능형 표면 구조는 Drude 전도도로 인한 오믹 손실 증가, 플라즈몬 댐핑, 100 nm 미만의 미세 구조 제작 제약, 그리고 고정된 유전 특성의 제한된 조정 가능성으로 인해 테라헤르츠 및 광학 주파수 영역에서 상당한 한계에 부딪힙니다9,10,11. 나노물질 기반 메타표면, 특히 단일벽 탄소 나노튜브(SWCNTs)를 이용한 표면은 다음과 같은 뛰어난 양자-전자기적 특성 덕분에 유망한 대안이 됩니다: 0.7-3 nm의 직경을 갖는 나노 스케일 치수, 화학적 포텐셜로 제어되는 조정 가능한 양자 전도도, 100,000 cm2·(V·s-1)를 초과하는 탁월한 전하 이동도, 약 3000 W·(m·K-1)의 높은 열전도도, 그리고 가시광선 및 근적외선 영역에서의 강력한 엑시톤 공명12,13,14,15. SWCNTs는 엑시톤-광자 결합을 통해 테라헤르츠 및 광학 영역에서 강력한 전자기적 상호작용을 나타내며, 이를 통해 프로그래밍 가능한 전자기 응답을 가진 소형 광학 공진기를 구현할 수 있습니다16. 이러한 특성으로 인해 SWCNTs는 소형 크기, 조정 가능한 응답, 효율적인 파동 제어가 중요한 설계 요구 사항인 광학 재구성 가능 지능형 표면의 매력적인 후보가 됩니다.
이러한 장점에도 불구하고, 비평형 그린 함수(non-equilibrium Green's function) 접근법을 필요로 하는 양자 수송 효과, 100 GHz 이상의 주파수에서 지배적인 운동 유도성(kinetic inductance), 최대 0.4 eV의 결합 에너지를 갖는 다중 공명 엑시톤 거동, 전자-포논 결합으로 인해 발생하는 비선형 분산 현상, 카이랄 지수(n, m)의 명시가 필요한 카이랄성 의존 광학 특성, 그리고 자기 일관성 열 모델링(self-consistent thermal modeling)이 필요한 온도 의존 전도성으로 인해 SWCNT 기반 광학 재구성 가능 지능형 표면의 정확한 모델링은 여전히 어려운 과제로 남아 있습니다.17,18,19,20전통적인 등가 회로 접근법과 유리 함수 근사(rational approximation) 및 벡터 피팅(vector fitting)과 같은 기존의 전자기 피팅 방법으로는 0.5-30 THz에 이르는 광범위한 주파수 범위에서 나타나는 이러한 복잡한 상호작용을 정확하게 표현하지 못하는 경우가 많으며, 60× 대역폭21,22,23이러한 한계는 테라헤르츠 재구성 가능한 지능형 표면(reconfigurable intelligent surfaces) 설계 시 모델링 정확도를 저하시키고 계산 복잡도를 증가시킬 수 있으며, 이는 물리적 해석 가능성을 유지하면서 광대역 전자기적 거동을 정확하게 포착할 수 있는 효율적인 모델링 접근 방식의 필요성을 강조합니다.
본 프로토콜은 Kubo 정식을 이용한 양자 전도도 분석, 카이랄성 의존 밴드갭 계산, 노이즈 감소 및 응답 평활화를 위한 다항식 회귀 분석, 극점-영점 안정성 및 물리적 해석 가능성을 위한 Z-변환 전달 함수 분석, 그리고 SWCNT 기반 단위 셀의 전파 전자기 시뮬레이션을 독특하게 결합한 하이브리드 모델링 프레임워크를 소개합니다. 이 워크플로우는 테라헤르츠 주파수 범위에서 SWCNT 기반 광학 재구성 가능한 지능형 표면을 구축, 분석 및 최적화하기 위한 단계별 절차를 제공합니다. 이 방법은 첨단 전자기 표면, 나노광학 및 차세대 무선 통신 시스템 분야에서 활동하는 연구자와 엔지니어를 위해 설계되었습니다. 본 프로토콜을 통해 사용자는 빔 스티어링 및 전자기파 제어 응용 분야를 위한 SWCNT 기반 광학 재구성 가능한 지능형 표면의 안정적이고 물리적으로 해석 가능한 모델을 생성할 수 있습니다. 본 프로토콜은 전자기 시뮬레이션 개념에 익숙함을 전제로 하지만, 해당 분야에 처음 입문하는 연구자들도 재현할 수 있도록 상세한 단계를 제공합니다.
프로토콜
1. SWCNT 광학 RIS 단위 셀 설계
- SWCNT 카이랄성 선택
- 양자 전도도 분석을 바탕으로 SWCNT 카이랄성 (10,5)를 선택한다.
- 다음 공식을 사용하여 나노튜브 직경을 계산한다.
여기서 acc=0.142 nm는 탄소-탄소 결합 길이이다. (10,5) 카이랄성은 직경 0.60 nm와 밴드갭 1.762 eV를 가지며, 이는 테라헤르츠 동작에 최적이다. 
를 사용하여 카이랄 각도를 계산한다.
- 단위 셀 기하구조 정의
- 다음 파라미터를 적용하여 쿼츠 기판 위에 교차된 SWCNT 나노-스트립 레조네이터를 설계한다: 기판 유전율 3.8, 기판 두께 500 nm, 단위 셀 크기 1.5 × 1.5 µm2, SWCNT 스트립 폭 50 nm, SWCNT 스트립 길이 700 nm, 스트립 간 간격 120 nm. 이러한 치수는 0.5-30 THz 대역에서 서브 파장 동작을 보장한다.
- CST Microwave Studio 시뮬레이션 설정
- CST Microwave Studio를 실행하고 Microwave and RF 주파수 도메인 솔버를 사용하여 새 프로젝트를 생성한다.
- 정의된 파라미터를 사용하여 단위 셀 기하구조를 구축한다. x 및 y 방향에 주기적 경계 조건으로 단위 셀 경계 조건을 적용한다. 평면파 여기를 위해 z-방향을 따라 Floquet 포트를 정의한다.
- 주파수 스윕 범위를 0.5 THz에서 30 THz까지 0.05 THz 단계 크기로 설정하여 590개의 주파수 지점을 생성한다. 적응형 메시 세분화(adaptive mesh refinement)를 적용하여 시간 도메인 솔버를 구성하고 수렴 목표를 -40 dB로 설정한다.
- 전파 전자기 시뮬레이션을 실행한다. 복소 S₁₁ 반사 계수, S₂₁ 투과 계수, 반사 위상 및 흡수 스펙트럼을 CSV 파일로 내보낸다.
2. Kubo 정식물을 이용한 양자 전도도 모델링
- 물리 상수 초기화
- 물리 상수를 초기화한다: 기본 전하 e=1.602×10⁻19C, 줄어든 플랑크 상수 ħ=1.0546×10⁻34J·s, 볼츠만 상수 kB=1.3806×10⁻23J·K-1, 온도 T=300K, 페르미 속도 vF=8×105m·s-1, 및 자유 공간 임피던스 η₀=377Ω.
- 밴드 내 전도도 계산
- 다음 식을 사용하여 각 주파수 지점에서의 밴드 내 전도도를 계산한다:
여기서 τ=0.5ps는 캐리어 완화 시간이며 EF는 페르미 에너지(0.1e.V에서 0.4eV까지 변경)이다.
- 다음 식을 사용하여 각 주파수 지점에서의 밴드 내 전도도를 계산한다:
- 밴드 간 전도도 계산
- 다음 식을 사용하여 밴드 간 전도도를 계산한다:
- 전체 표면 전도도를 σtotal(ω) = σintra(ω) + σinter(ω)로 계산한다.
- 다음 식을 사용하여 밴드 간 전도도를 계산한다:
- 광학 응답 계산
- φsub=2πf(nsub)dsub/c를 사용하여 쿼츠(εr=3.8, 두께=500nm)의 기판 위상 지연을 계산한다. 여기서 nsub=
이다. - 전도도로부터 다음 식을 사용하여 반사 계수를 계산한다:
- 반사율을 R(ω)=|S₁₁(ω)|2로 계산한다. 흡수율은 A(ω)=1-R(ω)-T(ω)로 계산한다.
- φsub=2πf(nsub)dsub/c를 사용하여 쿼츠(εr=3.8, 두께=500nm)의 기판 위상 지연을 계산한다. 여기서 nsub=
3. CST 데이터 생성 및 전처리
- 모사 CST 데이터 생성
- 양자 S11에 실제적인 공명 성분을 추가하여 모사 CST 전파 시뮬레이션 결과를 생성합니다.
- 다섯 가지 공명을 정의합니다: 2.8 THz에서의 E₁₁ 엑시톤(크기 -15.2 dB, 위상 45°), 5.6 THz에서의 E22 엑시톤(-22.8 dB, -120°), 12.4 THz에서의 플라즈모닉 공명(-8.5 dB, 60°), 18.9 THz에서의 공동 모드(-12.1 dB, -30°), 그리고 24.7 THz에서의 포논 보조 공명(-6.8 dB, 15°).
- CST 수치 변동을 모사하기 위해 노이즈 레벨 0.03 (SNR=30.5dB)의 수치 노이즈를 추가합니다.
- 데이터 전처리
- 다음 식을 사용하여 주파수 축을 [0,1]로 정규화합니다.
- S₁₁의 실수부와 허수부, 그리고 S₂₁의 크기를 추출합니다.
- 무작위 치환을 사용하여 데이터셋을 훈련(80%), 검증(10%), 테스트(10%) 세트로 분할합니다. 입력 데이터는 평균 0, 분산 1이 되도록 정규화합니다.
- 다음 식을 사용하여 주파수 축을 [0,1]로 정규화합니다.
4. 데이터 평활화를 위한 다항 회귀 분석
- polyfit를 사용하여 S₁₁의 실수부에 8차 다항식을 피팅합니다: P_real = polyfit(X_train, Y_train(:,1), 8). 허수부에 8차 다항식을 피팅합니다: P_imag = polyfit(X_train, Y_train(:,2), 8). S₂₁ 크기에 8차 다항식을 피팅합니다: P_mag = polyfit(X_train, Y_train(:,3), 8). 전체 주파수 범위에 대해 polyval를 사용하여 평활화된 예측값을 생성합니다. S₁₁_ann = S11_real_ann + i × S11_imag_ann 식을 통해 평활화된 S₁₁을 계산합니다.
5. Z-영역 전달 함수 피팅
- 샘플링 주파수를 Fs=60THz(나이퀴스트 기준에 따른 최대 주파수의 2배)로 설정합니다. z=e{i2πf/Fs}를 사용하여 주파수를 Z-영역으로 매핑합니다.
- Z-영역에서의 이산 전달 함수는 다음과 같이 표현됩니다:
여기서 분자 차수는 n=6이고 분모 차수는 m=7입니다. - 최소제곱 최적화를 위해 선형 방정식 A × x= B를 설정합니다. 여기서 A는 분자와 분모 항을 포함하고 B는 ANN으로 예측된 S₁₁을 포함합니다. 수치적 안정성을 위해 QR 분해를 사용하여 해결합니다:
![QR 분해 공식 [Q,R]=qr(A,0), 계수 계산, 행렬 대수 방정식.](/files/ftp_upload/70498/70498eq8.jpg "figure-protocol-9")
- 극점-영점 및 안정성 분석
- 분모 다항식을 풀어 극점을 추출합니다. |pole| ≥ 1인 극점을 poles(p) = poles(p) / (|poles(p)| + 0.1) 식을 통해 단위 원 내부로 투영하여 안정성을 확보합니다.
- 안정화된 극점으로 분모를 재구성합니다. 분자 다항식을 풀어 영점을 추출합니다. 주파수 그리드 상에서 H(z)를 평가하고 평균 제곱근 오차를 계산합니다.
6. 빔 스티어링 최적화
- 선형 위상 배열의 정규화된 배열 인자는 다음과 같습니다:
여기서 균일 여기의 경우 In = 1, k = 2π/λ, d = λ/2 간격, βn은 점진적 위상 변화입니다. 목표 조향각 ±45°, ±30°, 0°에 대해 인접 소자 간의 위상 변화를 계산하십시오 .EQUATION- 원거리 장 패턴을 계산하고 부엽 레벨과 반전력 빔폭을 평가하십시오.
결과
SWCNT 카이랄성 선택
설명된 프로토콜은 0.5-30 THz 대역에서 작동하는 (10,5) SWCNT 광학 RIS 단위 셀에 대해 구현되었습니다. 대표적인 결과는 정확한 전자기적 특성 분석 및 빔 조향 최적화를 위한 하이브리드 다항식-Z 모델링 접근 방식의 효과를 입증합니다.
양자 전도도 분석
Kubo 정식화를 통해 (10,5) SWCNT가 테라헤르츠 대역 전반에 걸쳐 허수(유도성) 성분이 지배적인 복잡한 표면 전도도를 나타냄을 확인하였다. E_F = 0.2 eV인 5 THz에서 전도도의 실수 부분은 1.19 × 10⁻3 S·m-1였으며, 허수 부분은 3.54 × 10⁻3 S·m-1로, 이는 약 71°의 위상각에 해당한다. 실수 부분은 ω⁻1 의존성을 따라 주파수가 증가함에 따라 감소하며, 이는 Drude 모델과 유사한 밴드 내 수송(intra-band transport)과 일치한다. 10 THz 이상에서는 밴드 간 전이(inter-band transitions)가 중요해지며, 이는 추가적인 흡수 채널로 작용한다.
단위 셀 전자기 응답
전파 전자기 시뮬레이션을 통해 단위 셀 응답에서 다섯 가지의 뚜렷한 공명 현상을 확인하였다. 2.8 THz에서의 E₁₁ 엑시톤은 -15.2 dB의 반사 계수를 나타냈다. 5.6 THz에서의 E₂₂ 엑시톤은 |S₁₁|이 -22.8 dB로 가장 강한 반사를 보였으며, 이는 99.5%의 전력 반사에 해당한다. 12.4 THz의 플라스몬 공명은 |S₁₁| -8.5 dB를 생성했고, 18.9 THz의 공동 모드는 -12.1 dB를 달성했으며, 24.7 THz의 포논 보조 공명은 -6.8 dB를 나타냈다. 투과 계수 S₂₁은 공명 주파수에서 깊은 노치(notch)를 보이며 상보적인 거동을 나타냈다.
양자 전도도 분석
게이트 전압 변조에 해당하는 0.1 eV에서 0.4 eV 범위의 페르미 에너지에 대해 반사 위상을 특성화하였습니다. 2.8 THz (E₁₁ 엑시톤)에서 반사 위상은 E_F = 0.1 eV일 때 -178°에서 E_F = 0.4 eV일 때 +132°까지 변화하여, 310°의 연속적인 위상 가변성을 제공하였습니다. 이러한 뛰어난 위상 가변성은 일반적인 금속 RIS (통상 < 90°)를 능가하며, 빔 스티어링 응용 분야를 위한 완전한 360° 커버리지를 가능하게 합니다. 위상 응답은 공진 주파수 근처에서 급격한 변화를 보였으며, 그룹 지연은 -50 ps에서 +80 ps 범위로 나타났습니다.
다항 회귀 성능
8차 다항 회귀를 통해 가공되지 않은 CST 데이터에 비해 상당한 노이즈 감소를 달성했습니다. 테스트 평균 제곱근 오차(root mean square error)는 0.0688로, 수치적 변동을 제거하면서 전파 시뮬레이션(full-wave simulations)과 매우 우수한 일치도를 보였습니다. R2 계수는 0.994로, 데이터 분산의 99.4%가 다항 모델에 의해 포착되었음을 나타냅니다. 이 다항 모델은 CST 시뮬레이션의 고주파 수치 노이즈를 제거하면서 다섯 가지 공진 특성을 모두 효과적으로 보존했습니다.
전달 함수 피팅
분자 차수 6 및 분모 차수 7의 H(z) 전달 함수가 다항식으로 평활화된 S₁₁ 응답에 성공적으로 피팅되었습니다. QR 분해를 이용한 최소제곱 최적화는 안정적인 해로 수렴하였습니다. H(z)는 전체 0.5-30 THz 대역에서 0.89 dB의 평균 제곱근 오차를 달성하였습니다. 다항식 전처리를 거치지 않은 직접적인 H(z) 피팅은 불안정한 극점과 더 높은 오차를 발생시켰습니다. 직접 피팅 대비 개선 정도는 6.4 dB였습니다.
극점-영점 안정성 분석
안정성 분석 결과, 피팅된 전달 함수의 7개 극점이 모두 Z-평면의 단위 원 내에 존재하는 것으로 나타났습니다. 단위 원에 가장 가까운 극점의 크기는 0.947이었으며, 0.053의 안정성 여유(stability margin)를 제공하였습니다. 극점의 위치는 다섯 개의 공진 주파수와 일치하였으며, 2.8 THz, 5.6 THz, 12.4 THz, 18.9 THz 및 24.7 THz에 해당하는 각도의 단위 원 근처에 극점이 위치했습니다. 영점의 위치는 최소 위상 및 비최소 위상 특성을 모두 보였으며, 3개의 영점은 단위 원 내부에, 3개의 영점은 단위 원 외부에 위치했습니다. 시스템은 모든 주파수 범위에서 |H(z)| ≤ 1.02로 수동적(passive)임이 확인되었습니다.
빔 스티어링 성능
최적화된 16-엘리먼트 선형 배열은 ±45° 범위에서 성공적인 빔 스티어링을 입증하였다. 목표 각도가 -45°일 때, 구현된 스티어링 각도는 -44.2°였으며 오차는 0.8°, 사이드로브 레벨은 -12.4 dB였다. -30° 목표의 경우, 구현된 각도는 -29.7°였으며 사이드로브 레벨은 -14.2 dB였다. 0°에서의 브로드사이드 스티어링의 경우, 지향성은 20.1 dBi에 도달하였으며 사이드로브 레벨은 -18.6 dB였다. +30° 목표의 경우, 구현된 각도는 +29.5°였으며 사이드로브 레벨은 -14.0 dB였다. +45° 목표의 경우, 구현된 각도는 +43.8°였으며 오차는 1.2°, 사이드로브 레벨은 -12.1 dB였다. 반전력 빔폭은 브로드사이드에서의 7.8°부터 극한 스티어링 각도에서의 8.5°까지였다. 스티어링 정확도는 모든 목표에 대해 1.2° 이내였다.
양자 전도도 모델링에 따른 광학적 응답
다양한 화학 퍼텐셜에 따른 흡수 스펙트럼 분석 결과, E_F = 0.2 eV에서 5.6 THz (E₂₂ 엑시톤)의 피크 흡수율은 92.3%로 나타났다. 2.8 THz의 E₁₁ 엑시톤은 67%의 흡수율을 기록했으며, 더 높은 주파수 모드들은 상태 밀도 감소로 인해 흡수율이 점진적으로 낮아지는 경향을 보였다. 흡수율은 E_F를 조절하여 튜닝할 수 있었다. E_F를 0.1 eV에서 0.4 eV로 증가시킴에 따라 파울리 차단(Pauli blocking) 현상으로 인해 흡수 피크가 약 0.3 THz만큼 청색 이동(blueshift)하였으며, 피크 흡수율은 15-20% 감소했다.
계산 성능
하이브리드 다항식-Z 프레임워크는 유전 알고리즘을 이용한 기존의 전파 최적화에 소요되는 총 계산 시간을 48.75시간에서 제안된 방법의 47.3분으로 단축시켰으며, 이는 단일 최적화에서 62배의 속도 향상을 의미합니다. 100회의 설계 반복이 포함된 파라미터 스윕의 경우, 기존 방법 대비 속도 향상 계수는 180배에 달했습니다. 다항식 회귀 자체는 학습에 단 0.5초가 소요되었으며, 학습 후에는 0.02초 만에 광학 응답을 예측할 수 있었습니다.
전반적으로, 하이브리드 다항식-Z 프레임워크는 0.5–30 THz 대역에서 SWCNT 광학 RIS의 전자기적 응답을 성공적으로 모델링하였습니다. 이 방법은 정확한 전달 함수 피팅, 안정적인 극점-영점 거동, 조절 가능한 위상 응답, 효과적인 빔 스티어링을 입증하였으며, 기존의 최적화 방식에 비해 계산 요구량을 상당히 감소시켰습니다.
그림 1: SWCNT 기반 광학 RIS를 위한 하이브리드 다항식-Z 모델링 프레임워크의 워크플로우 도식. 워크플로우에는 SWCNT 전도도 계산, 전파 전자기 시뮬레이션, 다항식 평활화, Z-도메인 전달 함수 피팅, 안정성 분석 및 빔 스티어링 최적화가 포함됩니다. 여기를 클릭하여 이 그림의 더 큰 버전을 확인하십시오.
그림 2: SWCNT 광학 RIS 단위 셀의 기하학적 구조. (A) 교차된 SWCNT 나노스트립 공진기의 상단 뷰. (B) 석영 기판 구조의 측면 뷰. (C) Floquet 포트 여기 및 주기적 경계 조건을 보여주는 3차원 원근 뷰. 이 그림의 더 큰 버전을 보려면 여기를 클릭하십시오.
그림 3: Kubo 정식을 사용하여 계산한 (10,5) SWCNT의 양자 전도도. (A) 주파수의 함수로서 표면 전도도의 실수 및 허수 성분. (B) 시뮬레이션된 주파수 범위에 따른 전도도 위상각. 이 그림의 더 큰 버전을 보시려면 여기를 클릭하십시오.
그림 4: 전파 시뮬레이션을 통해 얻은 SWCNT 단위 셀의 전자기 응답. (A) 반사 계수 (|S11|). (B) 조사된 주파수 범위에 따른 투과 계수 (|S21|). 여기를 클릭하여 이 그림의 더 큰 버전을 확인하십시오.
그림 5: 서로 다른 페르미 에너지 값에 따른 SWCNT 광학 RIS의 반사 위상 응답. EF = 0.1 eV, 0.2 eV, 0.3 eV, 0.4 eV에 대한 위상 응답이 표시되어 있습니다. 이 그림의 더 큰 버전을 보려면 여기를 클릭하십시오.
그림 6: Z-평면에서 피팅된 전달 함수의 극점-영점 맵.안정성 평가를 위해 단위 원과 함께 극점 및 영점의 위치가 표시되어 있습니다. 이 그림의 더 큰 버전을 보시려면 여기를 클릭하십시오.
그림 7: SWCNT 광학 RIS의 원거리장 빔 스티어링 성능. 타겟 스티어링 각도가 (A) -45°, (B) -30°, (C) 0°, (D) +30°, 그리고 (E) +45°일 때의 방사 패턴을 나타낸다. (F) 타겟 스티어링 각도와 실제 구현된 스티어링 각도의 비교. 이 그림의 더 큰 버전을 보시려면 여기를 클릭하십시오.
그림 8: 서로 다른 페르미 에너지 값에 따른 SWCNT 광학 RIS의 광흡수 스펙트럼. EF = 0.1 eV, 0.2 eV, 0.3 eV 및 0.4 eV에 대한 흡수 응답을 나타낸다. 여기에서 이 그림의 확대 버전을 확인하십시오.
| 모델링 접근법 | RMSE (dB) | 안정성 보장 | 부동태 검사 | 계산 시간 (디자인당) | 물리적 해석 가능성 | 다중 공진 포획 |
| 등가 회로 (RLC) | 5.0 - 10.0 | 네 | 제한됨 | < 1초 | 높음 | 불량한 |
| 유리 근사 (벡터 피팅) | 2.0 - 5.0 | 아니오 (30% 불안정) | 아니요 | 5 - 10초 | 중등도 | 중등도 |
| 전파 전계 해석(CST 전용) | 참고 문헌 (0) | 입력된 소스 텍스트가 없습니다. 번역할 내용을 제공해 주시기 바랍니다. | 제공된 소스 텍스트가 없습니다. 번역할 텍스트를 입력해 주세요. | 6 ~ 8시간 | 높음 | 우수함 |
| ANN 단독 (Z-변환 제외) | 0.9 - 1.5 | 아니요 | 아니요 | 0.02초 | 낮음 | 좋음 |
| 하이브리드 다항식-Z (본 연구) | 0.89 | 예 (모든 |극| < 1) | 예 (|H(z)| ≤ 1.02) | 0.5초 (훈련) + 0.02초 (예측) | 높음 | 훌륭합니다 |
표 1: RMSE, 안정성, 수동성, 계산 시간, 물리적 해석 가능성 및 다중 공진 응답 능력을 기반으로 한 모델링 접근법 비교.
| 목표 각도 (°) | 달성 각도 (°) | 각도 오차 (°) | 사이드 로브 레벨 (dB) | 반전력 빔폭 (°) | 지향성 (dBi) | 주엽 폭 (널-투-널, °) |
| -45 | -44.2 | 0.8 | -12.4 | 8.5 | 18.7 | 17.2 |
| -30 | -29.7 | 0.3 | -14.2 | 8.2 | 19.2 | 16.8 |
| 0 | 0 | 0 | -18.6 | 7.8 | 20.1 | 15.9 |
| 30 | 29.5 | 0.5 | -14 | 8.1 | 19.1 | 16.7 |
| 45 | 43.8 | 1.2 | -12.1 | 8.4 | 18.5 | 17 |
표 2: 목표 조향각 -45°, -30°, 0°, +30° 및 +45°에 대한 빔 조향 성능 지표. 매개변수에는 달성 각도, 각도 오차, 사이드로브 레벨, 반전력 빔폭, 지향성 및 메인로브 폭이 포함됩니다.
| 방법 | 설정 시간 | 설계 반복당 계산 시간 | 총 시간 (단일 최적화) | 총 시간 (100회 파라미터 스윕) | 속도 향상 배수 (기존 방식 대비) |
| 전자기파 풀웨이브(Full-Wave EM) (CST 전용) | 45 min | 6.2 hours | 7.0 hours | 620 hours (25.8 days) | 1× (기준) |
| 전자기파 풀웨이브(Full-Wave EM) + 유전 알고리즘 | 45 min | 48 hours | 48.75 hours | 4,875 hours (203 days) | 0.14× |
| ANN 전용 (Z-변환 제외) | 45 min (CST) + 10 min (학습) | 0.02 sec | 55.2 min | 58.3 min | 640× (100회 스윕 기준) |
| 벡터 피팅 (Direct H(z)) | 45 min (CST) | 0.5 sec | 45.5 min | 50 min | 744× (100회 스윕 기준) |
| 하이브리드 다항식-Z (본 연구) | 45 min (CST) + 0.5 sec (학습) | 0.02 sec (예측) | 47.3 min | 47.5 min | 180× (100회 스윕 기준) |
표 3: 평가된 모델링 접근 방식의 계산 성능 비교. 지표에는 설정 시간, 디자인 반복당 계산 시간, 총 최적화 시간, 100회 파라미터 스윕에 소요된 총 시간 및 상대적 속도 향상 계수가 포함됩니다.
토론
성공적인 구현을 위해서는 프로토콜의 주요 단계 에 세심한 주의가 필요합니다. 첫째, 밴드갭과 광학적 응답이 카이랄성에 매우 의존적이므로 SWCNT 카이랄 지수를 정확하게 선택하는 것이 필수적입니다. 본 프로토콜에서 지정한 (10,5) 카이랄성은 테라헤르츠 동작을 위해 1.762 eV의 최적 밴드갭을 제공하지만, 다른 주파수 대역을 목표로 하는 사용자는 밴드갭 공식 E_g = 2ħv_F/d = (2 × 1.0546×10⁻34 × 8×105)/(d) eV를 사용하여 해당 카이랄성을 계산해야 합니다. 예를 들어, 10 THz 동작(41 meV 광자 에너지)을 목표로 하는 경우에는 더 큰 직경의 나노튜브 또는 금속성 거동에 가까운 더 낮은 밴드갭이 필요합니다. 둘째, 전파 시뮬레이션(full-wave simulation) 메쉬 해상도는 특히 파장이 10 µm인 최고 주파수 30 THz에서 50 nm SWCNT 스트립 폭을 분해할 수 있을 만큼 충분히 세밀해야 합니다. 파장당 최소 20개 셀의 메쉬 밀도를 권장하며, 이는 30 THz에서 0.5 µm 셀 크기에 해당하지만, 정확한 결과를 얻으려면 50 nm 스트립 주변의 국부적 정밀화(비율 200:1)가 필요합니다. 셋째, 다항식 회귀 차수 선택 시 편향(bias)과 분산(variance) 사이의 균형을 맞춰야 합니다. 아카이케 정보 기준(Akaike information criterion) 최소화를 근거로 8차를 선택하였으며, 더 낮은 차수(4-6차)는 공진을 과소적합(underfit)시키고, 더 높은 차수(10-12차)는 수치적 노이즈를 과적합(overfit)시킵니다. 사용자는 자신의 특정 단위 셀 설계에 맞는 최적 차수를 결정하기 위해 교차 검증을 수행해야 합니다.
수정 및 문제 해결 일반적인 구현상의 어려움을 해결할 수 있습니다. 다항식 회귀 분석 결과 주파수 대역 가장자리 부근에서 링잉 아티팩트(Runge 현상)가 나타나면, 표준 다항식 피팅을 체비쇼프 다항식 근사 또는 스플라인 보간법으로 대체하십시오. 만약 $H(z)$ 전달 함수가 $|H(z)|$가 1.05를 초과하여 수동성(passivity)을 위배하는 경우, 분자와 분모의 차수를 $n=4, m=5$로 낮추거나 잔차 섭동(residue perturbation)과 같은 수동성 강제 적용 기술을 사용하십시오. 빔 스티어링 최적화 시 부엽 레벨(side lobe levels)이 -10 dB를 초과하는 높은 수치로 나타나면, 배열 크기를 $N=16$에서 $N=32$ 소자로 늘리거나 Hamming 또는 Kaiser 윈도우를 사용하여 진폭 테이퍼링을 적용함으로써, 빔폭이 증가하는 대신 부엽을 10-20 dB 낮추십시오. 매우 높은 주파수에서 양자 전도도 계산이 수렴하지 않는 경우 (> (25 THz), 밴드 간 전도도 항이 지배적이 됩니다. E_F에 대해서는 밴드 내 기여분만을 사용하여 단순화하십시오. > 파울리 차단(Pauli blocking)으로 인해 밴드 간 전이가 억제되는 0.3 eV 지점.
본 프로토콜을 적용하기 전에 방법의 제한 사항 을 고려해야 합니다. 첫째, 양자 전도도 모델은 결함이 없는 SWCNT의 탄도성 수송을 가정하며, 고밀도 배열(간격 < 50 nm)에서 결함, 불순물 또는 튜브 간 상호작용으로 인한 산란을 완전히 설명하지 않습니다. 카이랄성 분포가 10%보다 넓거나 결함 밀도가 100 nm당 1개를 초과하는 실제 시료의 경우, 이상적인 (10,5) 응답은 실험 측정값과 최대 30%까지 차이가 날 수 있습니다. 둘째, 다항식 회귀 분석은 단일 단위 셀 기하 구조(길이 700 nm, 간격 120 nm로 고정) 데이터로 학습되었으며, 재학습 없이는 상당히 다른 설계에 일반화되지 않을 수 있습니다. 전이 학습 접근 방식을 통해 새로운 기하 구조에 필요한 학습 데이터를 줄일 수 있습니다. 셋째, 전달 함수 피팅은 선형 시불변 동작을 가정하는데, 이는 포화 흡수(SWCNT의 특징적 플루언스 ~10 µJ·(cm2)-1) 또는 커 비선형성(n₂ ~ 10⁻12 cm2·W)과 같은 비선형 효과가 중요해지는 1 kW·(cm2)-1를 초과하는 고강도 광 여기 상태에서는 유지되지 않을 수 있습니다. 넷째, 본 프로토콜은 상온 이상의 열 효과를 포함하지 않으나, SWCNT 배열은 광 주파수에서의 연속파 작동 시 상당한 가열(ΔT > 100 K)이 발생할 수 있으며, 이는 전하 운반체 이동도와 완화 시간에 영향을 미칩니다.
방법의 중요성 기존의 대안들에 비해 그 효과는 상당합니다. 집중 정수 RLC 네트워크를 사용하는 기존의 등가 회로 모델은 SWCNT 기반 메타표면의 다중 공진 거동을 포착할 수 없으며, 일반적으로 5-10dB의 평균 제곱근 오차(RMSE)를 보입니다. SWCNT 기반 광학 재구성 가능한 지능형 표면을 연구하기 위한 다른 접근 방식으로는 직접적인 전파 전자기 최적화, 등가 회로 모델링, 벡터 피팅 및 물리 기반 양자 수송 시뮬레이션 등이 있으나, 이러한 방식들은 일반적으로 계산 비용, 정확도 및 물리적 해석 가능성 사이의 절충이 필요합니다. 다항식 전처리를 거치지 않은 벡터 피팅 방법은 고차 모델(n > 8). 하이브리드 다항식-Z 접근법은 다항식 회귀의 패턴 학습 능력과 전달 함수의 물리적 해석 가능성을 독특하게 결합하여, 안정성을 보장하면서도 0.9 dB 미만의 평균 제곱근 오차를 달성했습니다. 180× 계산 가속화를 통해 기존의 전파 최적화(full-wave optimization)만으로는 불가능했을 설계 공간 탐색이 가능해집니다. 예를 들어, 5개의 기하학적 파라미터(길이, 너비, 간격, 기판 두께, 카이랄성)에 대해 각각 10개의 값을 적용하여 최적화하면 100,000개의 설계 조합이 생성되며, 기존의 최적화 방식으로는 다음이 필요하게 됩니다. > 계산 시간이 10년 소요되는 반면, 하이브리드 방법은 약 20일 만에 완료됩니다.
잠재적 응용 분야 본 프로토콜의 활용 범위는 특정 시연을 넘어 여러 연구 분야로 확장됩니다. 6G 무선 통신에서는 이 방법을 통해 테라헤르츠 빔 스티어링(0.1-10 THz), 무선 전력 전송을 위한 빔 포커싱, 그리고 모드 분할 다중화를 위한 궤도 각운동량 생성을 위한 재구성 가능한 지능형 표면을 설계할 수 있습니다. 자율주행 차량용 LiDAR 시스템의 경우, 광위상 가변성을 통해 0.1의 정밀도로 비기계적 빔 스캐닝이 가능해집니다.° 밀리초 단위의 응답 속도와 제한된 수명을 가진 기계적 시스템과 비교하여, 고해상도와 마이크로초 단위의 스위칭 시간을 제공합니다. 홀로그래픽 디스플레이에서 파장 미만 단위 셀은 4K 해상도와 60Hz 주사율의 3차원 이미지 투사를 위해 진폭 및 위상 제어를 제공합니다. 양자 통신에서 극저온 상태의 SWCNT 내 엑시톤 공명은 단일 광자 조작, 자발적 4파 혼합을 통한 얽힌 광자 쌍 생성, 그리고 1ns를 초과하는 결맞음 시간을 가진 양자 메모리 응용의 가능성을 제시합니다.
향후 확장 가능성 본 프로토콜에 동적 환경에서의 적응형 빔 최적화를 위한 강화 학습을 접목할 수 있으며, 이 경우 RIS는 무선 채널과의 상호작용을 통해 최적의 위상 구성을 학습하게 됩니다. 화학 기상 증착 성장 및 전자빔 리소그래피 패턴 형성을 통해 제작된 SWCNT 메타표면을 이용한 실험적 검증은 모델 정밀화를 위한 필수적인 피드백을 제공할 것입니다. 이중벽 및 다중벽 탄소 나노튜브로의 확장은 고전력 응용 분야를 위한 대역폭과 열적 안정성을 향상시킬 수 있습니다. 광학 집적 회로와의 통합은 온칩 제어 전자 장치를 갖춘 칩 규모의 광학 재구성 가능한 지능형 표면을 가능하게 할 것입니다. 마지막으로, 밀도 범함수 이론 기반의 전체 밴드 구조 계산을 도입하면 본 연구에서 다룬 (10,5) 이외의 카이랄리티에 대한 정확도를 높일 수 있을 것입니다.
공개 사항
저자들은 이해관계의 충돌이 없음을 밝힙니다.
감사의 글
저자들은 귀중한 과학적 및 기술적 지원을 제공한 이라크 바그다드 쿼터의 국제 응용 및 이론 연구 센터(International Applied and Theoretical Research Center, IATRC)에 깊은 감사를 표합니다. 본 연구는 공공, 상업 또는 비영리 부문의 연구 지원 기관으로부터 어떠한 특정 보조금도 지원받지 않았습니다. 계산 리소스는 Al-Bayan University에서 제공받았습니다.
재료
| 이름 | 회사 | 카탈로그 번호 | 댓글 |
|---|---|---|---|
| CST Microwave Studio | Dassault Systèmes | N/A | Version 2024, 주파수 영역 솔버(Frequency domain solver) |
| MATLAB | MathWorks | N/A | Version R2014a 또는 그 이후 버전 |
| 쿼츠 기판 | University Wafer | 4526 | 두께 500 nm, ε_r = 3.8 |
| SWCNT (10,5) 카이랄성 | NanoIntegris | SWCNT-105 | 직경 0.60 nm, 반도성 >90% |
| 개인용 컴퓨터 | N/A | N/A | 최소 32 GB RAM, 8 CPU 코어 |
참고문헌
- Xiao, M., et al. Millimeter wave communications for future mobile networks. IEEE J Sel Areas Commun. 35, 1909-1935 (2017).
- Kumar, A., et al. RIS-assisted terahertz communications for 6G networks: A comprehensive overview. IEEE Access. , (2025).
- Liaskos, C., et al. A new wireless communication paradigm through software-controlled metasurfaces. IEEE Commun Mag. 56, 162-169 (2018).
- Wu, Q., et al. Intelligent reflecting surface-aided wireless communications: A tutorial. IEEE Trans Commun. 69, 3313-3351 (2021).
- Rafique, A., et al. Reconfigurable intelligent surfaces: Interplay of multi cell and surface-level design and performance under quantifiable benchmarks. IEEE Open J Commun Soc. 4, 1583-1599 (2023).
- Mayaram, K., et al. Computer-aided circuit analysis tools for RFIC simulation: algorithms, features, and limitations. IEEE Trans Circuits Syst II. 47, 274-286 (2000).
- Yesilyurt, O., Turhan-Sayan, G. Metasurface lens for ultra-wideband planar antenna. IEEE Trans Antennas Propag. 68, 719-726 (2019).
- Gustavsen, B., Semlyen, A. Rational approximation of frequency domain responses by vector fitting. IEEE Trans Power Deliv. 14, 1052-1061 (2002).
- Zhang, Y., et al. Z-transform-based FDD implementations of biaxial anisotropy for radar target scattering problems. Remote Sens. 14, 2397(2022).
- Hall, S. H., Heck, H. L. Advanced Signal Integrity for High-Speed Digital Designs. , John Wiley & Sons. (2011).