🔬 materials science

Irida-Graphene Phonon Thermal Transport via Non-equilibrium Molecular Dynamics Simulations

本研究利用非平衡分子动力学模拟揭示，新提出的二维碳同素异形体 Irida-Graphene 在室温下表现出约 215 W/mK 的本征热导率——由于其多孔的 3-6-8 环结构引起的声子散射以及降低的群速度，该热导率显著低于原始石墨烯——同时保持了各向同性的热传输特性及显著的尺寸效应。

原作者： Isaac M. Felix, Raphael M. Tromer, Leonardo D. Machado, Douglas S. Galvão, Luiz A. Ribeiro, Marcelo L. Pereira

发布于 2026-09-09

📖 1 分钟阅读☕ 轻松阅读

CC BY 4.0

原作者： Isaac M. Felix, Raphael M. Tromer, Leonardo D. Machado, Douglas S. Galvão, Luiz A. Ribeiro, Marcelo L. Pereira

原始论文采用 CC BY 4.0 许可（http://creativecommons.org/licenses/by/4.0/）。 ✨ 这是对下方论文的AI生成解释。它不是由作者撰写或认可的。如需技术准确性，请参阅原始论文。阅读完整免责声明

热量是微型机器世界中面临的一个持续挑战。随着电子设备缩小到沙粒大小，它们产生的热量变得越来越难以管理，往往威胁到要熔化它们原本旨在驱动的组件。为了解决这个问题，科学家们将目光投向了二维材料领域——这些原子层薄到本质上是平面的。在这些材料中，由碳原子排列成蜂窝状构成的单层结构，即石墨烯，长期以来一直是明星。它以惊人的效率传导热量，就像热能的超级高速公路。然而，研究人员一直在不断寻找这些材料的新变体，它们可能会提供不同的特性，或许是用某些极速性能换取其他有用的特征，或者仅仅是寻找控制热量如何在固体中移动的新方法。

在最近的一项研究中，一个研究小组将注意力转向了一种新提出的碳结构，称为 Irida-石墨烯（Irida-Graphene）。与石墨烯完美的蜂窝结构不同，这种材料由三角形、六边形和八边形的混合构成，创造出一种具有多孔、开放设计的薄片。科学家们想要了解热量是如何通过这种独特的晶格传播的。通过模拟原子运动的高性能计算机模拟，他们追踪了热能如何穿过该材料。他们发现，虽然 Irida-石墨烯仍然是良好的热导体，但其传导速度明显慢于其著名的近亲。在室温下，该材料的热传导率约为每米每开尔文 215 瓦。使用相同的方法测量时，这比纯净石墨烯观察到的约每米每开尔文 1,200 瓦有了大幅下降。

这种减速的原因在于材料自身的架构。在像石墨烯这样的完美薄片中，携带热量的振动（称为声子）可以长距离无阻碍地传播。而在 Irida-石墨烯中，不同环状形状的混合创造了一个崎岖的地形。当这些振动在材料中移动时，它们会不断撞击多孔设计带来的结构不规则性。这种散射效应减慢了它们的移动速度，就像一名跑者在充满障碍物的路径而非平坦跑道上奔跑一样。研究人员还观察到，与石墨烯相比，Irida-石墨烯中这些振动的传播速度较低，这进一步导致了热量移动能力的降低。

研究还表明，材料的大小至关重要。当薄片非常小时，热量以直接、无阻碍的方式移动。随着薄片变得更长，热量开始发生更多散射，转变为另一种传输模式。这种转变发生在一段特定的距离上，研究团队计算出该距离约为 55 纳米。尽管其导热性低于石墨烯，但该材料的表现仍优于许多其他二维物质，例如某些形式的硼或二硫化钼。这表明，在需要中等、受控的热流而非石墨烯那种极端速度的热管理系统中，Irida-石墨烯仍能找到一席之地，为设计下一代紧凑型电子设备的工程师提供一种新的工具。

技术摘要：通过非平衡分子动力学模拟研究 Irida-石墨烯的声子热输运

问题陈述
纳米材料的热管理对于纳米电子器件和热电器件的进步至关重要。尽管原始石墨烯以其卓越的热导率而闻名，但发现并表征具有多样化拓扑结构的各种新型二维（2D）碳同素异形体，对于扩展用于特定应用的材料库至关重要。最近，一种名为 Irida-石墨烯（Irida-G）的新型二维碳同素异形体被提出，其特征是由 3、6 和 8 元碳环组成的平坦片层拓扑结构。尽管 Irida-G 具有结构稳定性和金属特性，但其热输运特性仍有待探索。理解这些性质对于确定其在热管理应用中的可行性，以及理解其独特的孔隙结构如何影响与标准石墨烯相比的声子散射是必要的。

方法论
作者采用了使用 LAMMPS 软件包的经典反应分子动力学（MD）模拟。原子间作用力通过第二代反应经验键级（REBO）势能进行建模，该势能已通过 GULP 软件计算声子色散进行了验证；不存在虚频的结果证实了该模型的结构稳定性。

为了研究热输运，研究采用了反向非平衡分子动力学（RNEMD）方法。模拟设置包括：

系统配置： 在 $xy $平面内采用周期性边界条件，在$ z$ 方向采用自由边界。系统规模从 13,824 到超过 100 万个原子不等，长度从 10 nm 到 800 nm 不等。
热流诱导： 系统被划分为若干层。通过在“热区”和“冷区”之间交换动能来建立稳态温度梯度。
分析： 利用基于时间平均热流和温度梯度的傅里叶定律计算热导率 ( $\kappa$ )。研究还通过速度自相关函数分析了振动态密度（VDOS），并通过色散关系计算了声子群速度，以阐明热输运机制。

主要贡献与结果

本征热导率： 模拟显示，在室温下直接从 REBO 势能计算，Irida-G 具有约 130 W/mK 的本征热导率 ( $\kappa_\infty$ )。然而，考虑到经典力场往往会低估实验值（如石墨烯，REBO 势能得出约为 1200 W/mK，而实验值为 1500–5000 W/mK），作者应用了一种半经验修正方法。这种基于石墨烯已知偏差进行缩放的修正法表明，Irida-G 的“真实”本征热导率约为 215 W/mK。
尺寸效应与输运机制： 研究确定了取决于系统长度 ( $L$ $L$ ) 的三个不同热输运机制：
- 弹道输运机制 ( $L < 30$ nm)： 热导率随长度线性增加。
- 弹道-扩散过渡区 ( $30 < L < 100$ nm)： 这是一个对长度依赖性降低的过渡区域。
- 扩散输运机制 ( $L > 100$ nm)： 热导率在很大程度上不再随长度变化。
  有效声子平均自由程估计约为 55 nm。
降低机制： Irida-G 的热导率显著低于原始石墨烯（在相同的模拟框架下约为 1200 W/mK）。这种降低归因于其孔隙结构（三角形和八边形的排列）引起的增强声子散射。
声子动力学：
- 群速度： 计算出 Irida-G 的平均声子群速度为 410 m/s，明显低于石墨烯的 780 m/s。由于晶格热导率与群速度的平方成正比，这一差异是热导率降低的主要因素。
- 振动谱： VDOS 分析表明，弯曲（ZA）声子，特别是在低频（低于 10 THz）下的声子，是主要的载热子。这些模式表现出更长的波长和更少的散射，与在其他二维材料中观察到的趋势相似。
各向同性： 热输运性质被发现是几乎各向同性的，在 $x$ 和 $y$ 方向之间几乎没有差异。

意义与主张
本文声称，虽然 Irida-G 的热导率低于原始石墨烯，但其数值（~215 W/mK）仍与其他已知二维材料（如六方氮化硼 h-BN、PHA-石墨烯和双苯并烯网络）具有竞争力。研究表明，Irida-G 独特的拓扑排列引入了显著的声子散射，从而有效地调节了其热学性质。作者得出结论，Irida-G 在需要中等热导率的特定热管理应用中具有潜力，并具有其他新兴二维材料的特征。这项工作建立了关于碳同素异形体中结构孔隙率如何影响热输运的基础理解，弥合了理论提议与功能材料表征之间的差距。

类似论文