Gromacs 模拟硝酸盐的黏度偏差太大是什么原因?

peixing

老师您好：        我模拟NaNO3-KNO3-Ca(NO3)2体系黏度，采用0.8的RESP因子，150摄氏度与实验值偏差50%,550℃偏差-20%，为什么会出现低温正偏差高温负偏差的情况呢？密度偏差均在10%以内。调整过RESP因子出现类似的偏差规律，该怎么调整参数来减小偏差呢？

sobereva

本来粘度就是不容易算准的量，受非键作用（包括原子电荷和范德华作用）以及模拟粘度本身的手段很敏感。当前力场只能让密度准确到10%以内，想粘度能算得不错，基本不太可能。看看那些模拟熔盐粘度不错的文章里的力场参数怎么搞出来的。

peixing

  老师好   我用同样的力场模拟二元硝酸盐的密度黏度偏差都在20%以内，现在三元盐150度-6%，250度-18%，350度-1%，450度29%。我现在掌握到的经验是增大RESP电荷黏度增大，但是减小电荷会减小450度的偏差却增大250度偏差，这似乎形成了前后冲突。我想问一下老师，每个温度点的黏度偏差是否存在联系？模拟结束后怎么看黏度是否稳定了呢？谢谢老师！

sobereva

peixing 发表于 2022-5-13 15:55
老师好   我用同样的力场模拟二元硝酸盐的密度黏度偏差都在20%以内，现在三元盐150度-6%，250度-18%，350 ...

简单形式的势函数往往注定难以把不同温度区间的性质同时描述好。比如水模型，有的模型就专门是用来模拟冰点附近的性质的。
可以考虑针对不同温度区间给出不同的参数
另外，影响结果的非键相互作用的因素绝对不只是RESP电荷的整体的scale factor，LJ参数、不同原子的电荷相对比例，都有可优化的油水。可以以一系列温度下模拟的粘度和实验粘度的整体差值作为目标函数，通过非线性优化的方式对参数进行优化（但这就牵扯到自己写脚本了）。