이 콘텐츠를 보려면 JoVE 구독이 필요합니다. 로그인하거나 오늘 무료 평가판을 시작하세요.

방법 논문

테라헤르츠 빔 스티어링을 위한 탄소 나노튜브 기반 재구성 가능한 지능형 표면 모델링을 위한 하이브리드 ANN-Z 방법

53 조회수

⸱

DOI:

10.3791/70498

⸱

2026년 8월 28일

 ,  , 
1Technical Engineering College, Al-Bayan University, 2College of Computer Science and Mathematics, Tikrit University, 3Avicenna E-Learning Center, Tikrit University, 4Department of Automation and Artificial Intelligence Engineering, College of Information Engineering at Al-Nahrain University, 5International Applied and Theoretical Research Center (IATRC)

이 논문에서

요약

본 프로토콜은 6G 무선 통신 응용을 위해 테라헤르츠 대역(0.5-30 THz)에서 작동하는 단일벽 탄소 나노튜브 기반 광학 재구성 가능 지능형 표면의 정확한 전자기 모델링을 위한 하이브리드 인공 신경망 및 Z-변환 방법을 제시하며, 180× 반사 위상 튜닝 가능 범위가 XX를 초과하는 계산 가속화 310° 및 빔 조향 범위는 ±45°.

초록

단일벽 탄소 나노튜브 기반의 광학 재구성 가능 지능형 표면은 미래 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 단위 셀 설계

  1. SWCNT 카이랄성 선택
    1. 양자 전도도 분석을 바탕으로 SWCNT 카이랄성 (10,5)를 선택한다.
    2. 다음 공식을 사용하여 나노튜브 직경을 계산한다.
      육각형 격자 공식, d=acc*√(n²+nm+m²)/π, 구조 방정식 다이어그램.
      여기서 acc=0.142 nm는 탄소-탄소 결합 길이이다. (10,5) 카이랄성은 직경 0.60 nm와 밴드갭 1.762 eV를 가지며, 이는 테라헤르츠 동작에 최적이다.
    3. 결정학 각도 공식 θ = tan⁻¹(√3 × m/(2n + m)) = 23.4°, 수학 방정식.를 사용하여 카이랄 각도를 계산한다.
  2. 단위 셀 기하구조 정의
    1. 다음 파라미터를 적용하여 쿼츠 기판 위에 교차된 SWCNT 나노-스트립 레조네이터를 설계한다: 기판 유전율 3.8, 기판 두께 500 nm, 단위 셀 크기 1.5 × 1.5 µm2, SWCNT 스트립 폭 50 nm, SWCNT 스트립 길이 700 nm, 스트립 간 간격 120 nm. 이러한 치수는 0.5-30 THz 대역에서 서브 파장 동작을 보장한다.
  3. CST Microwave Studio 시뮬레이션 설정
    1. CST Microwave Studio를 실행하고 Microwave and RF 주파수 도메인 솔버를 사용하여 새 프로젝트를 생성한다.
    2. 정의된 파라미터를 사용하여 단위 셀 기하구조를 구축한다. x 및 y 방향에 주기적 경계 조건으로 단위 셀 경계 조건을 적용한다. 평면파 여기를 위해 z-방향을 따라 Floquet 포트를 정의한다.
    3. 주파수 스윕 범위를 0.5 THz에서 30 THz까지 0.05 THz 단계 크기로 설정하여 590개의 주파수 지점을 생성한다. 적응형 메시 세분화(adaptive mesh refinement)를 적용하여 시간 도메인 솔버를 구성하고 수렴 목표를 -40 dB로 설정한다.
    4. 전파 전자기 시뮬레이션을 실행한다. 복소 S₁₁ 반사 계수, S₂₁ 투과 계수, 반사 위상 및 흡수 스펙트럼을 CSV 파일로 내보낸다.

2. Kubo 정식물을 이용한 양자 전도도 모델링

  1. 물리 상수 초기화
    1. 물리 상수를 초기화한다: 기본 전하 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Ω.
  2. 밴드 내 전도도 계산
    1. 다음 식을 사용하여 각 주파수 지점에서의 밴드 내 전도도를 계산한다:
      정적 평형 공식, 복소 전도도 방정식, 물리 연구, 수학적 표기법.
      여기서 τ=0.5ps는 캐리어 완화 시간이며 EF는 페르미 에너지(0.1e.V에서 0.4eV까지 변경)이다.
  3. 밴드 간 전도도 계산
    1. 다음 식을 사용하여 밴드 간 전도도를 계산한다:
      복소 공식이 포함된 정적 평형 개념; 물리 분석을 위한 수학 방정식 포함.
    2. 전체 표면 전도도를 σtotal(ω) = σintra(ω) + σinter(ω)로 계산한다.
  4. 광학 응답 계산
    1. φsub=2πf(nsub)dsub/c를 사용하여 쿼츠(εr=3.8, 두께=500nm)의 기판 위상 지연을 계산한다. 여기서 nsub=비유전율의 제곱근, √εr, 공식; 전자기 이론, 재료 특성.이다.
    2. 전도도로부터 다음 식을 사용하여 반사 계수를 계산한다:
      분광 결과, S11(ω) 공식, 전자기파 거동 특성 분석.
    3. 반사율을 R(ω)=|S₁₁(ω)|2로 계산한다. 흡수율은 A(ω)=1-R(ω)-T(ω)로 계산한다.

3. CST 데이터 생성 및 전처리

  1. 모사 CST 데이터 생성
    1. 양자 S11에 실제적인 공명 성분을 추가하여 모사 CST 전파 시뮬레이션 결과를 생성합니다.
    2. 다섯 가지 공명을 정의합니다: 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°).
    3. CST 수치 변동을 모사하기 위해 노이즈 레벨 0.03 (SNR=30.5dB)의 수치 노이즈를 추가합니다.
  2. 데이터 전처리
    1. 다음 식을 사용하여 주파수 축을 [0,1]로 정규화합니다.
      정규화된 주파수 공식 \(z_{f_{norm}}=\frac{f-f_{min}}{f_{max}-f_{min}}\)을 나타내는 수식.
    2. S₁₁의 실수부와 허수부, 그리고 S₂₁의 크기를 추출합니다.
    3. 무작위 치환을 사용하여 데이터셋을 훈련(80%), 검증(10%), 테스트(10%) 세트로 분할합니다. 입력 데이터는 평균 0, 분산 1이 되도록 정규화합니다.

4. 데이터 평활화를 위한 다항 회귀 분석

  1. 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-영역 전달 함수 피팅

  1. 샘플링 주파수를 Fs=60THz(나이퀴스트 기준에 따른 최대 주파수의 2배)로 설정합니다. z=e{i2πf/Fs}를 사용하여 주파수를 Z-영역으로 매핑합니다.
  2. Z-영역에서의 이산 전달 함수는 다음과 같이 표현됩니다:
    디지털 신호 처리에서 사용되는 다항식이 포함된 이산 시간 전달 함수 공식 H(z).
    여기서 분자 차수는 n=6이고 분모 차수는 m=7입니다.
  3. 최소제곱 최적화를 위해 선형 방정식 A × x= B를 설정합니다. 여기서 A는 분자와 분모 항을 포함하고 B는 ANN으로 예측된 S₁₁을 포함합니다. 수치적 안정성을 위해 QR 분해를 사용하여 해결합니다:
    QR 분해 공식 [Q,R]=qr(A,0), 계수 계산, 행렬 대수 방정식.
  4. 극점-영점 및 안정성 분석
    1. 분모 다항식을 풀어 극점을 추출합니다. |pole| ≥ 1인 극점을 poles(p) = poles(p) / (|poles(p)| + 0.1) 식을 통해 단위 원 내부로 투영하여 안정성을 확보합니다.
    2. 안정화된 극점으로 분모를 재구성합니다. 분자 다항식을 풀어 영점을 추출합니다. 주파수 그리드 상에서 H(z)를 평가하고 평균 제곱근 오차를 계산합니다.

6. 빔 스티어링 최적화

  1. 선형 위상 배열의 정규화된 배열 인자는 다음과 같습니다:
    ​안테나 설계를 위한 배열 인자 방정식; AF(θ) 공식; 공학 수학 방정식
    여기서 균일 여기의 경우 In = 1, k = 2π/λ, d = λ/2 간격, βn은 점진적 위상 변화입니다. 목표 조향각 ±45°, ±30°, 0°에 대해 인접 소자 간의 위상 변화를 계산하십시오 .EQUATION
    1. 원거리 장 패턴을 계산하고 부엽 레벨과 반전력 빔폭을 평가하십시오.

액세스가 제한되었습니다. 이 콘텐츠를 보려면 로그인하거나 체험판을 시작하세요.

결과

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 ...

액세스가 제한되었습니다. 이 콘텐츠를 보려면 로그인하거나 체험판을 시작하세요.

토론

성공적인 구현을 위해서는 프로토콜의 주요 단계 에 세심한 주의가 필요합니다. 첫째, 밴드갭과 광학적 응답이 카이랄성에 매우 의존적이므로 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개 ...

액세스가 제한되었습니다. 이 콘텐츠를 보려면 로그인하거나 체험판을 시작하세요.

공개 사항

저자들은 이해관계의 충돌이 없음을 밝힙니다.

감사의 글

저자들은 귀중한 과학적 및 기술적 지원을 제공한 이라크 바그다드 쿼터의 국제 응용 및 이론 연구 센터(International Applied and Theoretical Research Center, IATRC)에 깊은 감사를 표합니다. 본 연구는 공공, 상업 또는 비영리 부문의 연구 지원 기관으로부터 어떠한 특정 보조금도 지원받지 않았습니다. 계산 리소스는 Al-Bayan University에서 제공받았습니다.

액세스가 제한되었습니다. 이 콘텐츠를 보려면 로그인하거나 체험판을 시작하세요.

재료

이 논문에 사용된 재료 목록
이름회사카탈로그 번호댓글
CST Microwave StudioDassault SystèmesN/AVersion 2024, 주파수 영역 솔버(Frequency domain solver)
MATLABMathWorksN/AVersion R2014a 또는 그 이후 버전
쿼츠 기판University Wafer4526두께 500 nm, ε_r = 3.8
SWCNT (10,5) 카이랄성NanoIntegrisSWCNT-105직경 0.60 nm, 반도성 >90%
개인용 컴퓨터N/AN/A최소 32 GB RAM, 8 CPU 코어

참고문헌

  1. Xiao, M., et al. Millimeter wave communications for future mobile networks. IEEE J Sel Areas Commun. 35, 1909-1935 (2017).
  2. Kumar, A., et al. RIS-assisted terahertz communications for 6G networks: A comprehensive overview. IEEE Access. , (2025).
  3. Liaskos, C., et al. A new wireless communication paradigm through software-controlled metasurfaces. IEEE Commun Mag. 56, 162-169 (2018).
  4. Wu, Q., et al. Intelligent reflecting surface-aided wireless communications: A tutorial. IEEE Trans Commun. 69, 3313-3351 (2021).
  5. 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).
  6. 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).
  7. Yesilyurt, O., Turhan-Sayan, G. Metasurface lens for ultra-wideband planar antenna. IEEE Trans Antennas Propag. 68, 719-726 (2019).
  8. Gustavsen, B., Semlyen, A. Rational approximation of frequency domain responses by vector fitting. IEEE Trans Power Deliv. 14, 1052-1061 (2002).
  9. Zhang, Y., et al. Z-transform-based FDD implementations of biaxial anisotropy for radar target scattering problems. Remote Sens. 14, 2397(2022).
  10. Hall, S. H., Heck, H. L. Advanced Signal Integrity for High-Speed Digital Designs. , John Wiley & Sons. (2011).

액세스가 제한되었습니다. 이 콘텐츠를 보려면 로그인하거나 체험판을 시작하세요.

재인쇄 및 허가

태그

탄소 나노튜브 표면양자 전도도 모델링쿠보 공식다항 회귀Z-영역 분석전자기 시뮬레이션빔 스티어링 최적화광자 메타표면
문의하기
도서관에 추천하기

개인정보 처리방침

이용 약관

정책

저작권 © 2026 MyJoVE Corporation. 모든 권리 보유.