力学性能—化学键断键的问题

ZZU_SCU

我想研究聚合的芳纶PPTA的力学性能——化学键断键，以如下模型化合物为例。
前期研究PPTA晶体的拉伸断裂位点为位点1，使用VASP的GGA和PBE泛函，（暂不展示数据）。

接着使用Gaussian的柔性扫描对两端的苯环进行拉伸，使用# opt=addgic um062x/def2tzvp guess=(mix,always) scf=(maxcycle=80,xqc)关键词。拉伸断裂位点为位点1。

经过计算几何优化平衡位置的键能，化学键位点1的键能高于化学键位点2的键能，即化学键位点2的键能低。经过计算柔性扫描能量最高点位置的键能，化学键位点2的键能高于化学键位点1的键能，即化学键位点1的键能低。使用方法m062x/def2tzvp。

问题1：随着化学键增长，键能变化的幅度（斜率）与什么物理参数有关？键力常数？
问题2：平衡位置的键能大小能说明热裂解所需要的能量么？
问题3：柔性扫描的拉伸断裂位点是否可以说明聚合物拉伸过程中化学键断裂？感谢解答

sobereva

可以用柔性扫描说明断裂位点

应当说具体键能是怎么算的，片段是否做了优化

ZZU_SCU

sobereva 发表于 2023-12-12 08:35
可以用柔性扫描说明断裂位点

应当说具体键能是怎么算的，片段是否做了优化

谢谢卢老师的回复。我计算了两次键能，第一次是几何优化opt后稳定的构象，且无虚频，计算了位点1和位点2的键能，此时位点2的键能低，第二次是采用柔性拉伸能量最高点的构象，没有考察频率，应该有虚频，计算了位点1和位点2的键能，此时位点1的键能低。我知道是构象不一致导致的结果，另外化学键被拉长最好考虑多参考态，我定性的考察并使用m062x/def2tzvp方法，研究这个体系是够了吧。我想知道柔性扫描中，在限制条件下，化学键被拉伸的长度以及能量的关系是受键力常数影响么？

键能的计算使用了BSSE矫正，没有优化片段，尽管知道这样做会浪费机时，但是定性比较是够了的，看过您的帖子《计算化学键键能时考虑BSSE不仅是多余的甚至是有害的》。使用的关键词是# um062x/def2tzvp guess=(mix,always) scf=(maxcycle=80,xqc) counterpoise=2。

sobereva

受键力常数影响