湖南大学：二维碳材料成键形式对热导率的影响规律！

摘要：Graphene+ 是一种具有sp²–sp³杂化的新型碳单层，具有独特的面外半负泊松比行为和类似石墨烯（graphene）的Dirac 特性 [Yu et al, Cell Reports Physical Science 3,100790 (2022)]。在本研究中，基于第一性原理计算，以graphene为基准，对比了包括graphene+在内具有不同杂化构型的二维碳材料，揭示了软化的sp²–sp³杂化键合对晶格热输运性质的显著影响。在室温下，graphene+的热导率为170 W/mK，远低于graphene的热导率 (~3170 W/mK)，原因在于graphene+中软化的sp²–sp³键显著抑制声子的振动，导致强的非谐性和弱的声子流体动力学效应。因此，在graphene+中，热导率大幅度降低源于其中软化的sp²–sp³键合网络。该研究为graphene+热输运性质提供了深入的物理洞悉，对其在热管理领域的潜在应用具有前瞻性指导意义。

引言

碳材料一直是材料科学、凝聚态物理和热管理相关领域的热点。例如，金刚石作为典型的三维碳材料，具有自然界最高的硬度和超高热导率。此外，石墨烯（graphene）作为典型的二维碳材料，具有狄拉克锥和超高热导率。随着碳材料中越来越丰富多样的键合形式，深入理解二维碳材料的热输运性质正成为一个重要问题。

对于二维碳材料，graphene中强sp²杂化的平面键合导致其超高热导率为 3000-5000 W/mK。Choudhry et al.等人报道了graphene同素异形体的热导率降低是由于声子晶体诱导的声子分支折叠，与键合强度无关[Phys. Rev. B 100, 165401 (2019)]。此外，T-graphene、D-graphene和biphenylene由于不同的 sp²键排列而表现出多样性的热输运行为行为，热导率分别为 800、600和 166 (254) W/mK。然而，在这些研究中没有考虑 sp³杂化效应。

先前报道的涉及碳材料热输运性质的研究主要集中在具有sp²杂化键合网络的碳单层上，而对于复杂的sp²–sp³（或 sp³）杂化键合如何影响热输运性质仍然知之甚少。

最近，我们在一种sp²–sp³杂化碳单层graphene+发现了由sp³杂化诱导的几何模式切换所导致的新型负泊松比行为，即面外的半负泊松比 [Yu et al, Cell Reports Physical Science 3, 100790 (2022)]，有望为graphene+带来更多丰富的应用和前景，如冲击吸收和能量吸收等。此外graphene+还表现出优异的电子狄拉克特性，典型表现为量子反常霍尔效应和高载流子迁移率。

然而，未知的热输运性质和对其机理的基本认识限制了它的进一步发展和可能的应用，并且graphene+中sp²–sp³杂化键合有望为复杂的sp²–sp³（或 sp³）杂化键合对热输运性质的影响机理提供更加深入的理解。

成果简介

近日，湖南大学的余林凤（第一作者）、湖南大学的陈艾伶（第二作者）、湘潭大学的王晓霞（第三作者）、湘潭大学的王慧敏（第四作者）、郑州大学的秦真真（第五作者）、湖南大学的秦光照（通讯作者）开展合作研究，基于第一性原理计算，系统报道了以新型二维sp²–sp³碳同素异形体（graphene+）为代表的系列二维碳材料热输运性质研究。研究表明，软化的sp²–sp³键合抑制了graphene+中主导热导率的低频声子（0-10THz），并进一步引起了强非谐性和弱声子流体动力效应。通过对模式级声子性质的分析，graphene+的强非谐性可以归因于格林奈森（Grüneisen）参数较大，克服了较小的散射相空间，从而导致弛豫时间较短。基于动量守恒过程（N和U过程）和声子泊肃叶流的指数因子， graphene+中的弱声子流体动力效应被成功捕获。由于较弱的声子流体动力学，在graphene+中存在较强的声子-声子散射。最后，基于势阱和原子均方位移量化了原子间键强度，并通过积分的晶体轨道哈密度顿布局（ICOHP）值确定了sp²–sp³和sp²键合强度，发现sp³原子会严重软化相邻的sp²–sp³键，导致热导率的降低。结合独特的力学拉胀性质和电子狄拉克特性，深入了解以graphene+为代表的二维材料热输运行为，可以为电子器件和相关的热管理应用提供指导。

该工作于近日在线发表于物理学国际权威期刊《PHYSICAL REVIEW B》:

Linfeng Yu, Ailing Chen, Xiaoxia Wang, Huimin Wang, Zhenzhen Qin, and Guangzhao Qin^*, Softened sp²–sp³ bonding network leads to strong anharmonicity and weak hydrodynamics in graphene+, Physical Review B106, 125410 (2022);

图1. 不同杂化形式的二维碳结构及热导率比较。（a）graphene、（b）graphene+、（c）penta-graphene、（d）biphenylene、（e）T-graphene的俯视图和侧视图；（f）不同杂化键合形式对热导率的影响规律。

Graphene+与graphene的区别在于graphene+中的五元环沿正方形的对角线分布，形成了sp²–sp³杂化的二维碳网络。Graphene+的热导率（κ）为170 W/mK，而graphene、T-graphene、biphenylene和Penta-graphene的 κ 分别为3170、645、421(267)和437 W/mK。低 κ 意味着Graphene+的非谐性更强。进一步通过对比来自迭代(ITE)和弛豫时间近似(RTA)方法的κ，我们发现sp³杂化的引入会导致值降低，表明ITE方法与RTA方法的差异变小。这意味着声子流体动力学的削弱。因此sp³杂化的引入会导致强的非谐性和弱的流体动力效应。

图2. （a）Graphene+的归一化力常数与截止半径的依赖关系，图中插图为graphene+在300 K时的热导率随截止半径的变化。（b）Graphene+和graphene的声子色散的对比。（c）Graphene+的声子的投影态密度，（d）Graphene+中的电子在不同电子轨道的分布。

归一化力常数变化不明显，意味着graphene+中存在弱的长程相互作用。整体上，graphene+和graphene的声子振动频率范围几乎一样(0-50 THz)，这是因为具有相同的碳原子。然而，由于独特的sp²–sp³网络，graphene+表现出更复杂的声子色散。其中对热导率占据主导地位的低频声子模式被软化(晶格振动较弱)，表明sp²–sp³键合减弱，反映在力学性能中graphene+的杨氏模量为120 N/m，远低于graphene (320-350 N/m)。与graphene中的纯的sp²杂化相比，sp²–sp³杂化的碳单层中的电子在轨道中的占据分布更加不均匀，因此可能提供更大的非简谐性。

图3. 温度效应.（a）Graphene+和graphene中热导率与温度的关系。（b）Graphene+和（c）graphene中不同声子分支对热导率的百分比贡献。

基于迭代（ITE）方法，graphene+和graphene的κ 分别为~ 170和~ 3170 W/mK。Graphene+的κ比graphene低一个数量级。弛豫时间近似（RTA）方法对graphene的 κ 值为440 W/mK，比ITE方法几乎低估了一个数量级，揭示了graphene中强烈的声子流体动力效应。相反，RTA方法计算的graphene+的 κ 为~110 W/mK，比迭代的 κ 低35%，这意味着graphene+中声子流体动力学较弱。graphene的 κ 主要来自于FA支声子，在200-800 K的温度范围内，FA支声子的 κ 几乎超过50%。与graphene类似，graphene+中的FA分支所占的比重也略大于其他分支。

图4. 室温下（300 K）模式级声子性质对比分析.（a）群速度，（b）弛豫时间，（c）散射相空间，（d）格林奈森参数，（e）累积热导率和（f）四声子散射相空间的对比.

在声子水平上，graphene的群速度远高于graphene+。graphene+中的FA支强的软化导致频率为0-10 THz的声子具有更低的群速度。进一步，graphene的弛豫时间远高于graphene+，即graphene+在低频区域存在较强的声子-声子散射。尽管graphene+具有比graphene更小的散射相空间，但大的Grüneisen参数克服了小的散射相空间导致强的非简谐性，由于软化的声子模式。当声子频率在0 – 10THz范围内时，graphene的 κ 增加最快，表明频率为0 – 10THz的声子主导晶格振动。因此，graphene+的 κ降低主要在于软化的低频声学声子，这与声子分支的贡献是一致的。最后我们评估了graphene+中的四声子的效应。在graphene+中，三声子散射主导声子-声子散射过程，较弱的四声子效应是由于相同原子质量的声子能带宽度较窄。因此，三声子散射主导了graphene+中的热输运性质，并且基于三声子散射得到的 κ 是可以接受的。

图5. （1）室温下归一化累积热导率作为平均自由程(MFP)的函数。Graphene+和graphene中N和U过程的对比分析。该插图是N与U过程的散射率之比(N/U)。橙色和绿色分别代表Graphene+和graphene。比较基于迭代(ITE)和弛豫时间近似(RTA)方法的石墨烯+和石墨烯的声子散射率。

Graphene+的平均自由程（MFP）略低于graphene。这意味着，软化的声子模式导致更小的MFP，这导致在热输运过程中声子的运动模式的改变。在graphene和graphene+中，N过程比U过程具有更强的散射率，graphene+中N过程对 κ 的贡献要弱于graphene。这表明与纯的sp²杂化相比，sp³杂化的引入显著抑制了声子流体动力学效应。强声子流体动力学导致graphene的MFP显著提高,这是由于RTA方法将N过程视为热阻，因此graphene中基于RTA的声子-声子散射被严重高估，这导致了对 κ 的严重低估。相反，graphene+的声子流体动力学较弱，IRE与RTA方法计算的κ与散射率改变相对较小。

图6. （a）势阱分析（b）基于分子动力学模拟计算的均方根位移(MSD)，（c）不同杂化键的积分的晶体轨道哈密顿布局(ICOHP)值的比较。

在相同的势阱深度下，graphene+的势阱越宽，表明其恢复到平衡位置所需的能量就越多，这说明graphene的恢复力小，原子键合较弱。此外，在模拟中，在graphene+中发现了比graphene更大的MSD，这意味着它在晶格振动中，偏离平衡位置更远，因此需要更多的能量来恢复偏差，即更弱的原子间键合。进一步，在graphene中发现了较低的ICOHP值，表明其轨道电子占据更大的键合态，即强的原子间键合。值得注意的是，在graphene+中，sp²–sp³键的ICOHP明显高于sp²键，而sp^2-1和sp^2-2键的ICOHP几乎没有变化，这意味着sp3杂化的引入主要软化了graphene+中的sp²–sp³键。因此，graphene+中比graphene低一个数量级的 κ 的起源可以从根本上解释为sp²–sp³键的软化。在graphene +中，sp³杂化使键合软化，导致强声子-声子散射，尤其是那些与sp³碳原子与相邻的原子的键合(sp²–sp³键)。

参考文献

[1] L. Yu, A. Chen, X. Wang, H. Wang, Z. Qin, and G. Qin, Softened sp²–sp³bonding network leads to strong anharmonicity and weak hydrodynamics in graphene+. Phys. Rev. B 106, 125410 (2022).

[2]L. Yu, Z. Qin, H. Wang, X. Zheng, and G. Qin, Half-negative poisson’s ratio in graphene+ with intrinsic dirac nodal loop, Cell Rep. Phys. Sci. 3, 100790 (2022).

湖南大学秦光照教授课题组诚聘博士后

平台情况：湖南大学机械与运载工程学院，和湖南大学汽车车身先进设计制造国家重点实验室，同时有国家超级计算长沙中心（湖大负责全面运营）的支持；湖大机械工程为国家“一流学科”，工程科学入选ESI前千分之一。实验室依山傍水，位于国家5A级风景区岳麓山下、湘江之滨（5分钟直达江边散步绿道）。长沙物价低，生活压力小，是中国最具幸福感城市和全国房价最低的省会城市之一。课题组与国内外多个一流的科研单位有长期稳定的合作关系，为团队成员的长远职业发展，包括出国交流等各方面，提供最大的指导和帮助。我们的研究方向为国际前沿，国家战略发展规划中的重点领域，湖南大学大力支持发展的新兴交叉前沿学科方向，易出重大成果。

合作导师：秦光照，男，任职于湖南大学机械与运载工程学院和汽车车身先进设计制造国家重点实验室，教授，博士生导师，岳麓学者，Vebleo Fellow，湖南省海外高层次人才，国家留学基金委建设高水平大学公派研究生项目及中外合作项目评审专家，担任多个知名国际期刊编辑和审稿人。获德国亚琛工业大学最佳博士论文奖（博士学位最高荣誉SUMMA CUM LAUDE），国家优秀自费（非公派）留学生奖学金，以及2021中国新锐科技人物卓越影响奖。曾先后在郑州大学、中国科学院、德国亚琛工业大学和美国南卡罗莱纳大学从事科研工作。主持中国国家自然科学基金项目（NSFC）、美国博士后创新项目、德国jara高性能计算项目等。主要从事微纳尺度传热及智能热管理方面的研究工作。共计合作出版发表高水平SCI论文80余篇，论文总被引2000+次，H因子25。另有学术专著3部，软件著作权2个。更多详细信息请见个人主页：http://qgz.sxl.cn/；http://grjl.hnu.edu.cn/p/2019216

研究领域及兴趣：使用超级计算机集群服务器，以高性能模拟仿真为基本手段，结合机器学习、多尺度模拟、3D打印、和热导率精密测量等关键技术，在微纳尺度及原子尺度下，对热智能系统进行设计与制造，以获得具有先进传热性能的材料和结构，并采用多种手段对其进行高效调控，广泛应用于国家战略发展规划中的多个重要领域。主要研究领域为微纳尺度传热及智能热管理技术，致力于从事微纳尺度下的热能传输过程、传热机制、和相关热功能材料及器件的研究与开发。在半导体微电子器件热管理、高效能量输运与转换、以及新型热能节能环保技术等领域具有广泛应用，例如新能源（如热电转换）、制冷、电子器件散热等。所开展的新材料热物性测量，能量传输机制建模和性质计算等，属于多学科交叉领域。

课题组在微纳尺度传热及智能热管理领域拟开展研究方向包括但不限于以下几方面 ：

1.先进热功能材料设计的高通量机器学习及数据驱动智能设计研究

2.多尺度微结构设计与微系统制造及其能量输运过程研究

3.探索超越传统手段的新型调控方法并探索传热性能调控极限

博士后职位要求：申请者需已经获得博士学位或即将获得博士学位，具有物理、材料、化学、机械、能动、电子、计算机或相关专业背景，在国际主流学术刊物上以第一作者或通讯作者发表过研究论文，具有独立思考和创造性思维能力，良好的中英文阅读和写作能力。

具有传热，固体物理，光学或力学等知识背景；有理论建模或计算、编程经验，例如第一性原理计算、分子动力学模拟等；或者有材料制备及相关热物性测量仪器操作经验；

其他相关情况：

1）按政策享受长沙市人才补贴（9-16万）。

2）博士后在站期间计算工作年限从进站之日起计算，博士后期间经评审后可认定为副研究员。

3）在站期间成果丰富，入选湖南大学“杰出博士后奖”(每年评选一次)者，可直接申请副教授岗位。

4）入选“博士后创新人才支持计划”、“博士后国际交流计划”（引进、派出项目）或同等竞争程度的博士后人才项目且研究成果突出者，出站后直接入选“杰出博士后”，申请留校不受本校来源比例限制。

5）合作导师将大力协助支持申请博士后科学基金、国家自然科学基金及省市各级课题。

6）在站期间，学校提供博士后公寓或租房补贴，其子女享受我校教职工子女入托、入学待遇。

7）博士后出站后，业绩优秀者可以续聘，或根据工作意愿推荐到合作单位任职。

8）可根据具体情况推荐到美国和欧洲知名高校课题组进行访问研究，海外留学经历受国家认可。

9）更多信息请见湖南大学官方公告

应聘及申请方式

招聘长期有效，请将个人简历以及相关材料（如代表作，研究简述等）发送至邮箱gzqin@hnu.edu.cn