本文献属于原创性研究（类型a），以下是根据其内容生成的中文学术报告。

研究报告：基于分子修饰与动态模拟的氢键网络结构设计实现聚酰亚胺纤维的高抗压缩性能

一、研究概况

本研究由北京化工大学、碳纤维及功能高分子教育部重点实验室的许宏杰（Hongjie Xu）、胡鹏浩（Penghao Hu）（并列第一作者），田国峰（Guofeng Tian）、齐胜利（Shengli Qi，共同通讯作者）及武德珍（Dezhen Wu）等人完成。研究成果发表于高分子科学领域权威期刊 《Polymer》 （第353卷，文献号129923），并于2026年3月27日在线发表。

二、研究背景与目的

高性能有机纤维，如高强高模聚酰亚胺（polyimide, PI）纤维，因其低密度、优异的热稳定性和耐环境老化性，在航空航天等极端工况下有着广泛应用。然而，与大多数有机纤维一样，其化学结构所固有的弱横向分子间作用力（如范德华力）导致了其抗压缩性能较差，这极大地限制了它们在结构承载领域的应用。如何突破这一化学结构的限制，开发出兼具高强高模与优异压缩性能的有机纤维，是该领域长期面临的核心挑战。

本研究的根本目的，在于利用分子动力学（molecular dynamics, MD）模拟，从原子和分子相互作用层面，系统性地探究并阐明提升聚酰亚胺纤维压缩性能的关键因素。研究团队希望借助聚酰亚胺结构的高度可设计性，通过计算机模拟高效筛选出具有高压缩潜力的化学结构，从而规避传统实验方法周期长、测试精度受限等问题，为新一代高抗压缩有机纤维的结构设计提供理论指导。

三、研究方法与详细流程

为实现上述目标，本研究设计了一套从模型构建到性能表征与机理分析的系统性计算模拟方法，具体步骤如下：

1. 化学结构与模型构建： 研究者首先设计了18种不同的聚酰亚胺化学结构。他们选用了两种常见的二酐单体：均苯四甲酸二酐（pyromellitic dianhydride, PMDA）和3,3’,4,4’-联苯四甲酸二酐（3,3′,4,4′-biphenyl tetracarboxylic dianhydride, BPDA），以及9种包含不同特征基团（如羟基、羧基、苯并咪唑等）的二胺单体，进行组合。每种聚合物模型包含10条分子链，每条链由20个重复结构单元（structural units）构成。 模拟流程首先构建了18个初始密度为0.1 g/cm³的非晶态模型（amorphous models）。之后，为模拟纤维中分子链高度取向的实际结构，使用Materials Studio软件中的“StressStrain.pl”脚本对非晶模型进行单轴拉伸，构建了拉伸比（draw ratio）从2到6的取向模型（oriented models）。这一步骤最终共生成了18个非晶模型和90个取向模型，合计108个模型，构成了一个高信息通量的研究库。

2. 单轴压缩模拟与模量计算： 这是本研究的一项技术创新。研究者通过修改“StressStrain.pl”脚本，开发了一套模拟单轴压缩的实验方法。其策略是：在模型的X轴方向上逐步施加递增的载荷使其变形，每一步变形后都进行动力学优化以获得真实的应力-应变数据。 为了比较不同体系的抗压性能，研究者提出并定义了“压缩模量”（compressive modulus）这一关键指标。他们将压缩过程前3%应变内的应力-应变数据进行线性拟合，所得斜率即为压缩模量。通过这一标准化的高通量模拟方法，得以快速、准确地评估所有108个模型的抗压缩潜力。

3. 结构与能量的多尺度分析： 在获得宏观压缩模量后，研究团队从“结构单元-分子链-亚微观结构”跨尺度视角，深入剖析了影响性能的内在机制。具体分析内容包括：

   • 结构参数分析：利用Forcite模块的长度分布函数，计算了结构单元的长度及其分布；利用回转半径函数计算了均方回转半径（mean square radius of gyration）、均方末端距（mean square end-to-end distance），并由此计算出特征比（characteristic ratio, C∞），用以衡量分子链的刚性。
   • 能量参数分析：通过能量计算工具，将体系总势能分解为成键作用与非键作用。非键作用进一步区分为分子间和分子内相互作用，并详细考察了范德华力（van der Waals forces）和静电作用能的贡献，以揭示拉伸和结构差异对体系稳定性的影响。
   • 氢键分析：特别对含有氢键供体与受体的体系，根据几何准则（氢受体距离≤2.5 Å且供体-氢-受体角度≥90°）精确计算了分子间和分子内氢键的数量和类型，以探明其在抵抗压缩形变中的作用。

四、主要结果与分析

1. 压缩模量的对比与筛选： 模拟结果清楚地显示，随着拉伸比的增加，所有体系的压缩模量均呈上升趋势。通过对比不同二酐体系发现，PMDA系列的整体压缩模量普遍高于BPDA系列。其中，PMDA/2,2’-二氨基-4,4’-联苯二酚（2,2′-ddb）体系表现最为突出，在拉伸比为4时，其压缩模量从5.52 GPa飙升至23.13 GPa，提升幅度高达319%。此外，另一个含羧基的体系PMDA/4,4’-二氨基联苯-2,2’-二甲酸（dabpa）也展现出超过20 GPa的理论压缩模量峰值。这两个体系被识别为具有最高抗压缩潜力的结构。

2. 分子链刚性的关键作用： 对18种体系特征比的分析表明，PMDA体系的分子链比BPDA体系更刚硬，这解释了其整体压缩模量更高的原因。进一步将体系分为柔性、半刚性和刚性三组后，发现刚性组中的PMDA/2,2’-ddb和PMDA/dabpa正是压缩模量最高的两个体系。这清晰地指向第一个关键结论：刚性的主链结构是获得高压缩模量的必要条件。刚性棒状单元的高旋转位垒，从结构单元层面限制了分子链在压缩应力下的形变。

3. 拉伸效应的结构演化机制： 以性能差异显著但化学组成完全相同的同分异构体PMDA/2,2′-ddb和PMDA/3,3′-ddb为对象，深入研究了拉伸的影响。结果发现，伴随拉伸，两个体系的结构单元都被均匀拉长，其长度分布向完全伸展的尺寸集中，同时分子链的取向度不断提高。当拉伸比达到4时，PMDA/2,2’-ddb体系的结构单元平均长度增长趋于稳定，初始的低长度侧峰消失，表明构象伸展接近极限。这恰好对应于其压缩模量达到峰值的时刻。

4. 能量与氢键网络的协同作用： 能量分析揭示了更深层的机制。在PMDA/2,2’-ddb体系中，分子间的范德华力在拉伸比为4时达到最大值（绝对值），与此同时，其分子内链段间的静电相互作用也显著强于PMDA/3,3’-ddb。这表明更强的电子共轭结构使构象调整需要更高的能量。 最关键的发现来自对氢键的分析。与直觉相反，尽管PMDA/3,3’-ddb拥有更多的总氢键和分子内氢键，但其压缩模量远低于PMDA/2,2’-ddb。深入分析发现，随着拉伸，PMDA/2,2’-ddb体系中的分子间氢键数量显著增加，而PMDA/3,3’-ddb主要形成分子内氢键。这揭示了一个重要的竞争关系：在氢键供体数量一定的前提下，分子内氢键与分子间氢键可能存在此消彼长的竞争。强化分子间氢键，能更有效地抵抗分子链间的滑移，从而显著提升体系抵抗纵向压缩变形的能力，这比增加分子内氢键更为关键。

五、研究结论与价值

本研究成功构建了一套从分子设计、高通量模拟到跨尺度机理分析的完整理论框架，得出以下核心结论：提升聚酰亚胺纤维压缩性能的关键不仅仅在于拉伸作用本身，其根本在于构建具备刚性主链、强分子内共轭结构以及丰富分子间氢键网络的化学体系。 本工作的科学价值在于，它首次通过系统的分子动力学模拟，清晰地区分并阐明了分子内与分子间氢键在高性能纤维压缩行为中的不同乃至竞争性角色，修正了“氢键数量越多，压缩性能越好”的简单认知。在应用价值上，该研究预测了PMDA/2,2’-ddb和PMDA/dabpa两种具有超过20 GPa超高理论压缩模量的聚酰亚胺结构，为实验化学家合成新型高抗压缩纤维指明了具体的分子设计方向，有望推动有机纤维在先进结构承载领域的更广泛应用。

六、研究亮点

• 创新性方法：开发了基于修改脚本的单轴压缩模拟方法，并提出了以“前3%应变拟合斜率”来定义和计算压缩模量的评估策略，为高分子材料压缩性能的模拟研究提供了新的标准化工具。
• 高通量设计：一次性构建了包含18种化学结构、6个拉伸梯度的108个模型，实现了对聚酰亚胺抗压性能的高通量筛选，极大提升了研究效率。
• 反直觉的机理发现：揭示了在抵抗压缩形变时，连接不同分子链的分子间氢键比仅局限于单链内部的分子内氢键更为有效，这一发现深化了对纤维中弱相互作用协同机制的理解，是该领域一个重要的概念性进步。