[论文解读] A cellular automaton model for thermal transport in low-dimensional systems

Alejandra León|arXiv (Cornell University)|Mar 20, 2026

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框架，能捕捉声子声学、光学和弯曲通道之间的能量转移。
• 在最小计算资源下实现对复杂几何、缺陷与不规则边界的建模。
• 用石墨烯纳米带对模型进行标定和验证，以再现宽度和温度对热导率的已知趋势。

提出的方法

• 定义一个CA，每个单元表示具有固定物种标签和坐标的原子位点，并进化变量为N_a、N_o、N_f、N_T、T。
• 使用更新规则 N_i(t+1)=N_i(t)+β[ ∑_{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}，以在弹道和扩散两种极限之间插值。
• 通过总占据量校准温度 T_i = α N_Ti，且用固定权重因子 w^(s) 将 N_Ti 在声学、光学和弯曲通道之间分配，α≈4.8。
• 将计算得到的能量交换与热导率 k = J / (dT/dx) 联系起来，其中 J 是热通量，dT/dx 是梯度，置于一维/二维离散框架内。
• 演示模型在实时映射局部振动激发和获得系统温度方面的能力，同时保持对大系统的线性O(N)可扩展性。

实验结果

研究问题

• RQ1 粗粒化的元胞自动机是否能够再现低维系统中声子介导热传导的弹道、扩散和过渡两种极限？
• RQ2 空位、杂质和不规则边缘在CA框架下如何影响石墨烯纳米带的热导率和局部温度分布？
• RQ3 CA模型是否能捕捉几何引起的效应，如瓶颈和拓扑变化（如S形带）对热传导的影响？
• RQ4 在真实缺陷密度下，当系统尺寸在横向或纵向增大时，CA模型的计算可扩展性如何？

主要发现

• CA reproduces已知趋势：在原始的AGNRs中，热导率随带宽增加而增加，随温度提高而降低。
• 缺陷与不规则边缘在研究温度范围内降低热导率，与无序性研究一致。
• 模型通过温度相关的β捕捉散射和Umklapp样效应，支持从弹道到扩散的平滑过渡。
• S形几何结构显示中间瓶颈行为，随中心区域长度变化，等效热导率降低，并可能呈现热隔离效应。
• 实现能够实现实时温度和热流可视化，且对大系统具有线性(O(N))可扩展性。