求助使用反应力场模拟两种加热方式下的键能

ocean6

各位老师好，最近我想用ReaxFF来模拟微波加热和传统NVT加热沥青过程中的产物各种化学键之间的能量，加热过程中为了尽可能的减小温度带来的影响（控制变量），NVT加热中采用线性拟合微波的温度变化，用来说明两种加热方式对于各种键之间的抑制和促进作用。我的问题有两个：
(1) 上述的这种思路是否正确，或者说可行
(2) 如果可行的话，我是不是要对反应后的产物片段进行提取，然后进行DFT的计算呢，DFT计算的时候是不是要考虑分子片段不要太大呢

ocean6

谢谢老师，关于我说的NVT加热线性拟合就是简介的您发的文献链接上的，reax_tools确实好用，您提的建议非常好，如果我想知道反应物的主要的化学键在哪个热解反应阶段的断裂情况，该怎样进行呢？比如说像下面这种

Graphite

什么是“NVT加热中采用线性拟合微波的温度变化”？
倒是有人用电场+NVE来表示微波加热（https://doi.org/10.1016/j.fuel.2023.127550），也有比较正弦电场、恒电场的。
如果要比较“各种化学键之间的能量”，lammps+reaxff本身不好直接提取，要么你写个脚本提取bonds文件的键级。要么考虑研究分子结构的变化、键数量的变化之类，用reaxtools。
有关反应后的产物片段，两种工况下得到的产物大概率不一样，比如微波可能导致结构更碎片化，或更多的5/7元环而不是6元环之类，关键产物的分子式都不一样，算DFT没法比。
不过说实话，如果只是为了定性或半定量说明两种工况的结果不一样，reaxtools基础功能的四张图就够了。http://bbs.keinsci.com/thread-54458-1-1.html