计算化学公社

标题: 求助用GAFF和opls力场模拟聚合物与实验性质不符怎么调参数？ [打印本页]

作者
Author: 万年流星 时间: 2020-5-27 22:08
标题: 求助用GAFF和opls力场模拟聚合物与实验性质不符怎么调参数？
分别用GAFF和opls力场模拟了聚合物的玻璃化转变温度，发现均比实测Tg大，请问要怎么修改参数呢？

作者
Author: SherryLiu 时间: 2020-5-28 03:43
不知道楼主用的Charge是什么是参数自带的还是自己计算拟合的 OPLS力场的话一般用的是Gaussian算ChelpG 然后0.8*ChelpG进行修改

作者
Author: 万年流星 时间: 2020-5-28 07:12

SherryLiu 发表于 2020-5-28 03:43
不知道楼主用的Charge是什么是参数自带的还是自己计算拟合的 OPLS力场的话一般用的是Gaussian算ChelpG 然 ...

GAFF力场用的RESP电荷，OPLS力场用的CM5*1.2

作者
Author: naoki 时间: 2020-5-28 09:43
聚合物的Tg不太好模拟准确，Tg和实际研究聚合物的分子量有很大关系，不知你模拟和实测的分子量是否一致。而且升温和降温以及升降温速率快慢对测得Tg都有不少影响。

作者
Author: 苏玖染 时间: 2020-5-28 15:49
OPLS CM5*1.2原文也是做各种溶剂下的模拟的，模拟高分子熔体是否合适呢？

作者
Author: 万年流星 时间: 2020-5-28 16:30

苏玖染发表于 2020-5-28 15:49
OPLS CM5*1.2原文也是做各种溶剂下的模拟的，模拟高分子熔体是否合适呢？

我再优化下结构，用resp电荷做下看看

作者
Author: 万年流星 时间: 2020-5-28 16:35

naoki 发表于 2020-5-28 09:43
聚合物的Tg不太好模拟准确，Tg和实际研究聚合物的分子量有很大关系，不知你模拟和实测的分子量是否一致。而 ...

分子量选的低聚物的，5个重复单元和50个重复单元的，Tg均比实验值高20-30K，升温退火按照文献中的方法（有用materials studio，compass力场做70个重复单元的，只比实验值高几K）

作者
Author: 苏玖染 时间: 2020-5-28 16:39

万年流星发表于 2020-5-28 16:35
分子量选的低聚物的，5个重复单元和50个重复单元的，Tg均比实验值高20-30K，升温退火按照文献中的方法（ ...

不做退火，直接恒温恒压平衡模拟，到平衡状态，计算平均密度

作者
Author: wwww0wwww 时间: 2020-7-15 12:03
楼主你好，我目前也遇到类似问题，不知道应该如何优化已有的力场文件，尤其是对于聚合物熔体体系，好像没有现成的针对这类体系的电荷优化方法，似乎只有COMPASSII针对聚合物体系做过一些优化，不过参数不公开是个比较头疼的事情，不知到楼主用OPLS+RESP电荷效果怎么样？

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