计算化学公社

标题: 求助合并Reaxff力场以及如何优化 [打印本页]

作者
Author: TinnKuka 时间: 2022-3-3 18:06
标题: 求助合并Reaxff力场以及如何优化
各位老师好，我正在学习如何针对我的体系优化Reaxff力场。我有两个体系，一个是PBT在O2中，一个是PBT和阻燃剂分子在O2中。
对第一个体系，我只需要CHO力场，所以我参考https://doi.org/10.1021/jp709896w这篇文献，觉得这个力场适合我的体系。但是我不明白如何优化力场参数，所以我的第一个问题是，利用DFT的数据，我尽可能让reaxff的数据靠近它，例如在sigma bond优化的时候，是修改D_e^sigma(sigma bond dissociation energy)这个值吗？

我的第二个问题是，O2不需要优化，但PBT分子中有苯环的存在，需要将这个力场对环进行bond、angle、torsion的优化吗？

我的第三个问题是，因为在第二个体系中，我将引入其他元素，例如P。但是并没有相似的力场，所以我需要将P的力场与CHO力场合并，那么合并力场优化的话也是需要对阻燃剂分子中P的bond、angle、torsion都要优化吗？

我的表述可能不太清晰，如果有什么问题，请告诉我，我可以补充修改，谢谢大家！

作者
Author: 413 时间: 2022-3-3 19:53
原则上要重新拟合。
但是现在reaxff力场的问题是大家用的乱七八糟

作者
Author: TinnKuka 时间: 2022-3-3 21:13

413 发表于 2022-3-3 19:53
原则上要重新拟合。
但是现在reaxff力场的问题是大家用的乱七八糟

是的，很多文献大家都是直接采用力场，但我试过直接用原始力场对oligomer，得到的结果并不太合理。请问您有没有力场拟合的相关方法、手册呢？我想学习一下。谢谢您！

作者
Author: 413 时间: 2022-3-3 21:17

TinnKuka 发表于 2022-3-3 21:13
是的，很多文献大家都是直接采用力场，但我试过直接用原始力场对oligomer，得到的结果并不太合理。请问您 ...

拟合没有手册，都写在文献里了。你看下03年那篇原始文献，再看一篇caltech的有一个教授发的jctc。

作者
Author: TinnKuka 时间: 2022-3-3 21:32

413 发表于 2022-3-3 21:17
拟合没有手册，都写在文献里了。你看下03年那篇原始文献，再看一篇caltech的有一个教授发的jctc。

好的，非常感谢您:)！请问第二篇是这篇https://doi.org/10.1021/ct5001044吗？

作者
Author: 413 时间: 2022-3-3 21:51

TinnKuka 发表于 2022-3-3 21:32
好的，非常感谢您:)！请问第二篇是这篇https://doi.org/10.1021/ct5001044吗？

是的。

作者
Author: TinnKuka 时间: 2022-3-3 21:55

413 发表于 2022-3-3 21:51
是的。

谢谢您！

作者
Author: mengxiangidea 时间: 2022-3-28 20:41
楼主你好，能请教你关于反应力场的东西吗

作者
Author: liyuan（小白） 时间: 2025-8-21 20:09
楼主有没有优化出吡啶基团的反应力场；想知道优化一个力场需要多久？有什么需要注意的问题吗

作者
Author: TinnKuka 时间: 2025-9-10 15:57

liyuan（小白）发表于 2025-8-21 20:09
楼主有没有优化出吡啶基团的反应力场；想知道优化一个力场需要多久？有什么需要注意的问题吗

我试了很多种优化的方法，都不怎么行，还是需要有经验的人来做。找van duin组合作吧

欢迎光临计算化学公社 (http://bbs.keinsci.com/)