[논문 리뷰] A cellular automaton model for thermal transport in low-dimensional systems

Alejandra León|arXiv (Cornell University)|2026. 03. 20.

Thermal properties of materials인용 수 0

한 줄 요약

본 논문은 셀룰러 오토마톤 모델을 제시하여 보손-매개 열전달을 저차원 나노구조에서 발리스틱, 확산, 전이 거동에 걸쳐 시뮬레이션하고 그래핀 나노리본에서 검증하며 대규모 시스템으로 확장 가능함을 보인다.

ABSTRACT

In this work, we formulate a theoretical model based on a cellular automaton (CA) to study thermal transport in low-dimensional nanostructures across ballistic, diffusive, and transition regimes. Unlike computationally intensive methods such as the Boltzmann Transport Equation (BTE), our model stands out for its geometrical robustness, allowing the seamless integration of substitutional impurities, vacancies, and irregular edges. We validated the model using graphene nanoribbons (AGNRs), successfully replicating the dependence of thermal conductivity on ribbon width and temperature. Results demonstrate that the model captures critical scattering and confinement effects with a linear scalability O(N). Given the increasing pressure to optimize computational resources and reduce the carbon footprint associated with AI infrastructure, this CA model emerges as a highly efficient tool for the parametric exploration and design of next-generation thermal devices.

연구 동기 및 목표

• 나선형 AI 에너지 수요 증가 속에서 나노구조에서의 음향 전달에 대한 저비용이며 물리적으로 일관된 모델의 필요성을 동기화한다.
• 음향, 광학, 구형(플렉셜) 음향 채널 간의 에너지 전이를 포착하는 거시적 CA 프레임워크를 개발한다.
• 복잡한 기하학, 결함, 불규칙 경계에 대해 최소한의 계산 자원으로 모델링을 가능하게 한다.
• 그래핀 나노리본을 사용해 폭과 온도에 따른 열전도도에 관한 알려진 경향을 재현하도록 모델을 보정하고 검증한다.

제안 방법

• 각 셀마다 고정된 종(label)과 좌표를 가진 원자 사이트를 나타내고 진화 변수 N_a, N_o, N_f, N_T, T를 갖는 CA를 정의한다.
• update rule N_i(t+1)=N_i(t)+β[ sum_{j in n(i)} N_j(t) − z_i N(t) ]를 사용해 서로 이웃한 셀 간에 진동 에너지를 재분배한다.
• β(T_m)=β_0 σ(T_m)로 온도 의존적 결합을 도입하고 σ(T_m) = (1 + e^{(T_m−T_0)/ΔT})^{-1}를 사용해 보 ballistic 과 확산 거동 사이를 보간한다.
• 온도는 총 점유에서 T_i = α N_Ti로 보정하고 α ≈ 4.8이며 N_Ti를 음향, 광학, 구형 채널에 고정 가중치 계수 w^(s)로 분배한다.
• 계산된 에너지 교환은 1D/2D 이산 프레임워크 내에서 열전도도 k = J / (dT/dx)로, 여기서 J는 열 플럭스이고 dT/dx는 기울기와 연관된다.
• 모델이 로컬 진동 여기도를 실시간으로 매핑하고 시스템 온도를 얻는 능력을 시연하며 O(N) 확장성을 유지한다.

실험 결과

연구 질문

• RQ1 저차원 시스템에서 음손-매개 열전달의 발리스틱, 확산, 전이 거동을 거친 CA가 재현될 수 있는가?
• RQ2 결함, 부품, 불규칙한 가장자리가 CA 프레임워크 내에서 그래핀 나노리본의 열전도도 및 로컬 온도 분포에 어떤 영향을 미치는가?
• RQ3 S 모양 링과 같은 기하학적 요인이 열전달에 어떤 병목 현상 및 위상적 변화를 유발하는가?
• RQ4 실제 결함 밀도 하에서 시스템 크기가 횡방향 또는 종방향으로 증가할 때 CA 모델의 계산적 확장성은 어떠한가?

주요 결과

• CA는 알려진 경향을 재현한다: 순수한 AGNR에서 폭이 증가하면 열전도도가 증가하고 온도가 증가하면 감소한다.
• 결함의 공석 및 불규칙한 가장자리는 연구 온도 범위에서 열전도도를 감소시키며 무질서 연구와 일치한다.
• 온도 의존적 β를 통해 산란 및 Umklapp 비슷한 효과를 포착하고 매끄럽게 발리스틱에서 확산으로의 전이를 지원한다.
• S-모양 기하학은 중앙 병목 현상을 보여 중앙 영역 길이가 바뀜에 따라 유효 열전도도가 감소하는 경향을 나타내며 열 차단 효과를 보일 수 있다.
• 구현은 대규모 시스템에 대해 선형( O(N) ) 확장성과 함께 실시간 온도 및 열 플럭스 시각화를 시연한다.