Thermal Analog Computing: Application to Matrix-vector Multiplication with Inverse-designed Metastructures

기술 요약: 역설계 메타구조를 통한 열 아날로그 컴퓨팅

문제 정의
현대 클래식 컴퓨팅은 이진 논리와 디지털 스위칭에 의존하며, 이는 저항 손실과 반복적인 메모리 접근으로 인해 근본적인 에너지 한계에 직면해 있습니다. 아날로그 컴퓨팅 패러다임(예: 멤리스터, 광학 컴퓨팅)은 머신러닝 및 신호 처리의 핵심 연산인 행렬-벡터 곱셈(MVM)을 위한 에너지 효율적인 대안을 제공하지만, 열(heat)은 여전히 미개척 영역으로 남아 있습니다. 기존의 열 메타물질은 주로 열 흐름을 유도(예: 클로킹)하는 데 집중할 뿐, 계산을 수행하는 데는 초점을 맞추지 않았습니다. 또한, 전통적인 열 논리 게이트는 디지털의 한계를 그대로 반영하여 불연속적인 "온/냉(hot/cold)" 상태로 작동합니다. 본 연구는 입력 온도를 출력 전력으로 매핑하는 목표 전도도 행렬을 생성하기 위한 기하학적 구조를 결정하는 역문제(inverse problem), 즉 열 전도를 신호 전달 매체로 사용하여 연속적인 MVM을 수행하는 과제를 해결합니다.

방법론
저자들은 열화된 금속 포트에 온도장을 인가하여 입력 벡터를 인코딩하고, 지정된 단자에서 수집된 열 전력으로 출력 벡터를 나타내는 열 아날로그 컴퓨팅(thermal analog computing) 프레임워크를 제안합니다. 이 시스템은 푸리에 법칙(Fourier's law)의 선형 영역 하에서 작동하며, $P_{out} = M T_{in}$이라는 관계식을 성립시킵니다. 여기서 $M$은 오직 기하학적 구조와 재료 분포에 의존하는 전도도 행렬입니다.

주요 방법론적 구성 요소는 다음과 같습니다:

• 위상 최적화를 통한 역설계: 저자들은 2D 메타구조를 합성하기 위해 밀도 기반 위상 최적화(density-based topology optimization)를 사용합니다. 재료 밀도 $\rho(x,y)$는 고체(전도성)와 공극(void) 사이를 전환하는 연속 변수로 취급됩니다.
• 음수 행렬 요소 처리: 열 전도도는 음수가 될 수 없으므로, 본 프레임워크는 임의의 목표 행렬 $M$을 $K$개의 물리적으로 구현 가능한 양수 성분의 가중 합인 $M = \sum z_i A(\rho_i)$로 분해합니다. 양수 및 음수 계수($z_i$)는 별도의 물리적 구조로 구현되며, 최종 결과는 각 전력 출력의 디지털 사후 감산을 통해 얻어집니다.
• 미분 가능한 솔버 파이프라인: 최적화에는 자체 개발한 JAX 기반의 미분 가능한 유한 체적 열 솔버(finite-volume thermal solver)가 사용됩니다. 이를 통해 설계 파라미터(밀도 및 두께)에 대한 비용 함수의 기울기를 자동 미분을 통해 계산할 수 있는 엔드 투 엔드 미분 가능 파이프라인을 구축할 수 있습니다.
• 최적화 제약 조건: 비용 함수는 목표 행리와 재구성된 행렬 사이의 가중 $L_2$ 노름을 최소화합니다. 추가적인 제약 조건으로는 제조 가능한 바이너리 구조를 보장하기 위한 서브픽셀 평활 투영(subpixel-smoothed projection, SSP)과 비현실적인 열 플럭스 집중을 방지하기 위한 정규화가 포함됩니다.

주요 결과
본 프레임워크는 $2\times2$ 및 $3\times3$ 차원의 행렬 풀을 통해 검증되었으며, 대부분의 경우 99% 이상의 정확도를 달-성했습니다.

• 단일 구조 설계: 단순한 행렬(Identity, Directional, Two-Branch)은 단일 구조($K=1$)로 구현되었으며, 오차 1% 미만을 기록했습니다. 그러나 요소의 크기가 균일하거나 오프-대각(off-diagonal) 항이 지배적인 행렬은 단일 구조로 구현하기에 어려움이 있었습니다.
• 다중 구조 설계: 혼합된 부호와 큰 크기 차이를 가진 복잡한 행렬들이 여러 구조($K=2$ ~ $K=4$)를 사용하여 성공적으로 구현되었습니다.
  • 하다마르 행렬(Hadamard Matrix): $K=2$ 구조를 사용하여 구현되었습니다 (오차 $\approx 1\%$).
  • 회전 행렬(Rotational Matrices): $2\times2$ 회전($\theta=\pi/3$) 및 $3\times3$ 회전 토플리츠(Rotational Toeplitz) 행렬은 $K=3$ 구조를 사용하여 구현되었습니다 (오차 $\approx 3\%$).
  • 이산 푸리에 변환(DFT): $3\times3$ DFT 행렬의 실수부와 허수부는 각각 $K=4$와 $K=2$ 구조로 분해되어 구현되었으며, 각각 3.55%와 1.13%의 오차를 보였습니다.
• 성능 지표:
  • 대역폭: 열 확산에 의해 제한되며, $100\,\mu\text{m}$의 횡방향 치수에 대해 약 $1.9\,\text{kHz}$로 추정됩니다.
  • 에너지 효율: $3\times3$ 행렬 기준 약 $9.3\,\text{nJ/MAC}$로 추정됩니다. 전자 또는 광학 가속기보다 높을 수는 있으나, 본 시스템은 근본적으로 에너지 수동형(energy-passive)이며, 행렬 연산 자체를 위해 능동적인 전력을 필요로 하지 않습니다.
  • 견고성: 시스템은 최대 100 K의 온도 차이에 대해서도 선형성(1~5% 오차 이내)을 유지합니다. 상온에서의 열적 요동에 의한 신호 대 잡음비(SNR)는 약 $10^6$으로 추정됩니다.

의의 및 주장
본 논문은 열 흐름이 고충실도 선형 변환을 수행하도록 설계될 수 있음을 입증함으로써 열 아날로그 컴퓨팅의 개념 증명(proof-of-concept)을 확립했다고 주장합니다. 본 연구의 의의는 다음과 같습니다:

1. 패러다임의 전환: 열을 단순히 부산물이나 이진 스위치로 사용하는 것을 넘어, 연속적인 아날로그 신호 전달 매체로 활용합니다.
2. 인시튜(In-situ) 컴퓨팅: 열이 이미 존재하는 환경(예: 마이크로 일렉트로닉스)에서 계산을 수행할 수 있게 하여, 계산을 위한 추가적인 에너지 오버헤드 없이 온도 구배 감지 및 열 제어를 가능하게 합니다.
3. 설계 방법론: 임의의 행렬 연산을 위해 복잡한 열 메타구조를 합성할 수 있는 강력하고 미분 가능한 역설계 파이프라인을 제공합니다.

저자들은 본 구조가 속도와 에너지 밀도 제한으로 인해 고처리량의 전자 또는 광학 가속기를 대체하지는 못할 것이라고 겸허히 언급하면서도, 계산과 소산(dissipation)이 공존하는 보완적인 영역을 열어준다고 밝혔습니다. 향후 연구는 나노 스케일(포논 볼츠만 수송 방정식 해결)로 방법론을 확장하고 더 높은 차원의 행렬로 확장하는 방향으로 진행될 예정입니다.