davi 发表于 2024-10-31 22:57
老师，您意思是直接用文献中阴阳离子的所有参数就行了吗，我看的文献说的要貌似考虑极化再×0.8，那对于 ...

算力允许，同时跑一把，测试一下g(r)稳定性就有数了，CMD跑起来也快

davi 发表于 2024-10-31 20:21
老师，我还想请教个问题。我想要计算的是LiFSI盐在醚类混合溶剂(1:1 THF:2MeTHF)中Li+ g(r)的情况。是否 ...

阴阳离子没有必要，因为力场开发的时候甚至还考虑了阴离子的极化，比RESP应该更精确。中性的有机分子可以全部用RESP，不讲究的话用真空环境就行，讲究的话用0RESP2(0.5)。

slxc920113 发表于 2024-10-31 09:56
AuToFF也可以计算CM1A电荷，LJ对于拟合参数的溶液环境没有那么大的影响。

老师，我还想请教个问题。我想要计算的是LiFSI盐在醚类混合溶剂(1:1 THF:2MeTHF)中Li+ g(r)的情况。是否应该是先去结构优化Li+与FSI-离子(再上溶剂环境)，优化完用sob老师程序计算resp电荷就行，对于溶剂，也是一样的道理(不需要设置溶剂环境了吧)，然后算resp电荷，最后将文献的力场参数(除了电荷)补充上就行

slxc920113 发表于 2024-10-31 09:56
AuToFF也可以计算CM1A电荷，LJ对于拟合参数的溶液环境没有那么大的影响。

老师，如果我将离子电荷按0.7缩放了，还有必要进行溶剂分子的resp电荷计算吗

davi 发表于 2024-10-31 14:53
明白了，谢谢老师，我先直接试试CM1A看看效果。老师，还有个问题，Li+采用Autoff产生的LJ参数，我去看了 ...

AuToFF也可以计算CM1A电荷，LJ对于拟合参数的溶液环境没有那么大的影响。

slxc920113 发表于 2024-10-31 09:50
CM1A电荷对于常规的有机分子还能接受，但是我一般保险起见会替换成RESP电荷。

明白了，谢谢老师，我先直接试试CM1A看看效果。老师，还有个问题，Li+采用Autoff产生的LJ参数，我去看了看引用文献，貌似是做的是溶剂水，那我是四氢呋喃这种溶剂，LJ参数是不是不太可靠。

davi 发表于 2024-10-31 14:47
对于四氢呋喃和二甲基四氢呋喃这种，直接采用LigParGen生成的参数可靠吗，老师。我用Autoff产生的唯一区 ...

CM1A电荷对于常规的有机分子还能接受，但是我一般保险起见会替换成RESP电荷。

slxc920113 发表于 2024-10-31 09:33
OPLS系列的力场自带电荷能不用就不要用（无法保证总的电荷正确），LigPargen是用CM1A方法计算的电荷覆盖 ...

对于四氢呋喃和二甲基四氢呋喃这种，直接采用LigParGen生成的参数可靠吗，老师。我用Autoff产生的唯一区别仅有电荷。感觉直接用Autoff的OPLS/M应该是可行的吧

davi 发表于 2024-10-31 14:27
老师，您好。我看了看Autoff选OPLS/M产生的参数是一样的，但是电荷貌似不太一样，这是为什么呢？第二个问 ...

OPLS系列的力场自带电荷能不用就不要用（无法保证总的电荷正确），LigPargen是用CM1A方法计算的电荷覆盖原来自带的电荷，但是这个方法对于特殊的成键电荷误差也很大，如果条件允许就自己计算RESP电荷。
OPLS-AA/M比OPLS-AA/L新，优先选择。但是支持的原子类型比OPLS-AA/L少不少，可以将OPLS-AA/L作为缺失参数时候的备选力场。
离子体系建议用CL&P等专用的力场。

slxc920113 发表于 2024-10-31 09:21
保持OPLS-AA的版本相同，Ligpargen的参数来自OPLS-AA/M，在AuToFF中你需要选择OPLS-AA/M，而不是OPLS-AA/L

老师，您好。我看了看Autoff选OPLS/M产生的参数是一样的，但是电荷貌似不太一样，这是为什么呢？第二个问题是OPLS/L或者OPLS/M我到底该用哪一个呢，不太理解选取的依据是啥。我目前是想要做LiFSI溶质在THF/2MeTHF里面计算配位数。

保持OPLS-AA的版本相同，Ligpargen的参数来自OPLS-AA/M，在AuToFF中你需要选择OPLS-AA/M，而不是OPLS-AA/L