多晶B6N6、BN、BC3和C3N4纳米片的力学性能、断裂韧性和裂纹扩展行为的分子动力学研究

《Results in Surfaces and Interfaces》：Molecular dynamics study on the mechanical properties, fracture toughness and crack propagation behavior of polycrystalline B6N6, BN, BC3, and C3N4 nanosheets

　　

本研究采用分子动力学（MD）模拟来研究多晶（polycrystalline）B₆N₆、BN、BC₃和C₃N₄纳米片和纳米管的力学性能、断裂韧性（fracture toughness）和裂纹扩展行为。系统分析了晶粒数量、温度、应变率、边缘预裂纹和圆形切口缺陷对杨氏模量（Young's modulus）、抗拉强度（tensile strength）、失效应变（failure strain）和临界断裂韧性（critical fracture toughness）的影响。结果表明，增加晶粒数量会降低抗拉强度和杨氏模量，其中BN纳米片表现出最高的整体力学性能。多晶B₆N₆表现出最大的断裂韧性，其次是BN、BC₃和C₃N₄，反映了原子结构和键合在失效前能量吸收中的作用。热分析表明，BN由于原子振动减少和更强的B-N键，在高温下保持优异的机械稳定性。对预裂纹和切口的分析表明，BN对裂纹扩展最敏感，而BC₃和C₃N₄表现出裂纹不敏感行为，失效通常从晶界（grain boundary）而不是裂纹尖端开始。应变率效应表明，较高的应变率通过限制原子重排和裂纹生长来提高断裂韧性。对于纳米管，增加直径会提高杨氏模量，但会降低抗拉强度、失效应变和断裂韧性，其中C₃N₄纳米管对温度最敏感。这些发现为多晶纳米片和纳米管的力学行为提供了详细的见解，指导设计具有优化强度、韧性和热稳定性的纳米材料，用于先进应用。

**论文解读：多晶B₆N₆、BN、BC₃和C₃N₄纳米片及纳米管的力学性能与断裂行为分子动力学研究**

**研究背景与意义**

二维（2D）材料在实际合成中（如化学气相沉积、剥离或溶液处理）往往形成多晶结构，包含大量晶界（grain boundary，GB）和原子尺度缺陷，这些缺陷显著改变其物理、化学及电子性质。尽管对单晶石墨烯及六方氮化硼（hexagonal boron nitride，h-BN）等材料已有大量研究，但关于多晶硼碳氮（B-C-N）化合物（如B₆N₆、BN、BC₃、C₃N₄）的系统比较仍十分匮乏。这些材料在高温电子器件、热管理系统、柔性光电器件及光催化领域有广泛应用，其力学性能（强度、韧性）和热稳定性是实际可靠性的关键。然而，晶粒尺寸、温度、缺陷类型及应变率如何影响多晶B-C-N纳米片和纳米管的断裂韧性及裂纹扩展行为，尚缺乏统一理解。为此，研究人员（Morteza Ghorbanzadeh Ahangari 等）通过分子动力学（Molecular Dynamics，MD）模拟，系统对比四种多晶纳米片（B₆N₆、BN、BC₃、C₃N₄）及相应纳米管的力学性能、断裂韧性和裂纹扩展机制，旨在为设计具有优化强度与热稳定性的纳米材料提供理论指导。该研究成果发表在《Results in Surfaces and Interfaces》。

**主要技术方法**

研究人员采用大型原子/分子并行模拟器（LAMMPS）进行MD模拟，原子间相互作用使用优化的Tersoff势函数。多晶纳米片通过AtomSK软件利用Voronoi分割方法生成，晶粒数量范围4–28，纳米片尺寸约为150×150??2。多晶纳米管由平面纳米片沿特定方向卷曲而成，直径范围12.7–31.8??，长度50??。模拟施加重力场下的单轴拉伸，应变率范围2.5×10⁶–2.5×10⁹?s^-1，温度控制范围100–1500?K（采用Nose-Hoover热浴）。为研究缺陷影响，引入边缘预裂纹（长度12–40??）和圆形切口（直径9–23??）缺陷。使用Open Visualization Tool（OVITO）可视化原子构型和变形机制。

**研究结果**

**1. 晶粒数量对力学性能的影响**
通过比较4–28个晶粒的多晶纳米片，发现增加晶粒数量导致所有材料的抗拉强度和杨氏模量下降。BN纳米片表现出最高的整体力学性能（24晶粒时抗拉强度75.70?GPa，杨氏模量356.30?GPa），而C₃N₄最低。B₆N₆尽管模量低（245.15?GPa），但断裂韧性最高（15.9?J/m³），归因于其八面体环拓扑结构在拉伸时发生非仿射重排，延迟失效并吸收更多能量。BN的强度降低幅度最小（约15%），归因于其低能晶界与部分离子键；C₃N₄降幅最大（约18%），源于大环结构导致的悬空键弱连接。

**2. 温度效应**
在100–1500?K范围内，所有材料的杨氏模量和抗拉强度均随温度升高而降低。BN的热稳定性最高（模量降幅15%，强度降幅58%），B₆N₆的强度最敏感（降幅84%）。温度敏感性顺序与声子频率和键合非调和性相关：BN因高光学声子频率（E_2g ~1370?cm^-1）和低非调和性而稳定；B₆N₆的八面体环引入低频扭转模式（~200–400?cm^-1），促进热激活提前断裂。

**3. 边缘预裂纹对断裂韧性的影响**
引入长度12–40??的边缘预裂纹后，各材料的临界断裂韧性（K_IC）随裂纹长度线性增加（不符合线弹性断裂力学），但BN的K_IC值最高，且裂纹敏感性最大。B₆N₆的K_IC略低，但其应力-应变曲线积分得到的断裂能量更高，体现能量吸收与裂纹尖端应力集中之间的互补关系。BC₃和C₃N₄表现出裂纹不敏感行为，失效常起始于晶界而非裂纹尖端。裂纹宽度增加也提高K_IC，时间敏感性顺序为C₃N₄ > BC₃ > B₆N₆ > BN。

**4. 圆形切口缺陷的影响**
增大圆形切口直径（9–23??）使杨氏模量和抗拉强度线性下降，下降斜率顺序为C₃N₄ < BC₃ < B₆N₆ < BN。BN对切口最敏感，裂纹从切口径向扩展；C₃N₄最不敏感，裂纹扩展速度在BC₃中最高。

**5. 多晶纳米管的力学性能**
对于直径12.7–31.8??的纳米管，增大直径提高杨氏模量但降低抗拉强度、失效应变和断裂韧性（韧性下降）。C₃N₄纳米管性能变化率最大（模量变化49%，强度变化41%），B₆N₆变化率最小。

**6. 温度对纳米管的影响**
在100–1500?K范围内，纳米管的力学性能下降，C₃N₄纳米管模量降幅最大（70%），B₆N₆最小（34%）。C₃N₄纳米管对温度最敏感。

**讨论与结论**

研究人员指出，多晶纳米片和纳米管的力学行为受晶界密度、温度、应变率及缺陷形态的复杂耦合控制。BN在强度、模量和热稳定性方面综合最优，适合需要刚度和强度的应用；B₆N₆因高断裂韧性适用于能量吸收结构；BC₃和C₃N₄对裂纹不敏感，但绝对性能较低。研究成果否定了缺陷总是降低性能的简单认知，揭示了不同化学配比带来的独特变形机制（如八面体铰链效应、大环呼吸模式等）。

**结论翻译**（摘自原文“4. Conclusion”节选）：
MD模拟揭示，多晶B₆N₆、BN、BC₃和C₃N₄纳米片和纳米管的力学性能、断裂韧性和裂纹扩展行为受晶粒数量、温度、应变率和结构缺陷的强烈影响。BN纳米片表现出最高的杨氏模量、抗拉强度和热稳定性，而B₆N₆表现出最大的断裂韧性。对预裂纹和圆形切口的分析表明，BN对裂纹扩展最敏感，而BC₃和C₃N₄表现出裂纹不敏感断裂行为，失效常起始于晶界而非裂纹尖端。提高应变率通过限制原子重排和裂纹生长增强断裂韧性。对于纳米管，增大直径提高刚度但降低抗拉强度、失效应变和断裂韧性，其中C₃N₄纳米管对温度最敏感。