GMX模拟熔盐时是否有合适的势函数呢？

Streicheln · 发表于 Post on 前天 23:57

本帖最后由 Streicheln 于 2026-5-9 15:25 编辑

最近刚学完sob老师的分子动力学与GROMACS培训班，想自己计算一个简单的熔盐体系（LiCl-KCl）。使用gmx自带的LJ势势函数不支持熔盐计算，所得到的扩散系数(4.125e-9)偏小(实际值1e-5左右)，而BMH势往往用于Lammps。现在想问各位老师gmx是否有合适的势函数，还是必须要自己编写itp文件。

sobereva · 发表于 Post on yesterday 10:09

不存在什么gromacs初级班，看清楚名字和里面的内容层次介绍：
北京科音分子动力学与GROMACS培训班（http://www.keinsci.com/KGMX）
从来没说过这叫什么初级班，明明包含了大量中级程度的内容

不说具体是什么熔盐，没法回答。本来靠原子电荷和GROMACS支持的LJ势也不是不能描述熔盐，例如：

Streicheln · 发表于 Post on yesterday 15:36

sobereva 发表于 2026-5-9 10:09
不存在什么gromacs初级班，看清楚名字和里面的内容层次介绍：
北京科音分子动力学与GROMACS培训班（http:/ ...

好的sob老师，我已对帖子修改。我目前就是使用LJ势函数对KCl-LiCl熔盐进行的计算，但是计算得到的值(4.125e-9)偏小(实际值1e-5左右)，我不太懂这是什么原因造成的呢？

sobereva · 发表于 Post on yesterday 20:07

Streicheln 发表于 2026-5-9 15:36
好的sob老师，我已对帖子修改。我目前就是使用LJ势函数对KCl-LiCl熔盐进行的计算，但是计算得到的值(4.12 ...

可以搜索并尝试其它能描述这种体系的力场
大部分Li、K、Cl的力场参数是描述水中离子状态（而且往往还只是低浓度状态）的，这样的参数拿来描述熔盐肯定不准
倘若实在找不到gromacs支持的参数形式的话再考虑用其它程序，诸如CP2K支持的力场函数形式丰富得多。或者尝试把某种势函数自己拟合成LJ或者gromacs支持的其它势函数形式

liuyuje714 · 发表于 Post on yesterday 20:30

低版本gmx （gmx<=2019，用cutoff-scheme=group）是支持buckingham势函数的，可以转成gmx单位进行模拟这种熔融盐。但是不支持gpu加速计算，整体模拟速度非常慢。

Streicheln · 发表于 Post on yesterday 21:25

sobereva 发表于 2026-5-9 20:07
可以搜索并尝试其它能描述这种体系的力场
大部分Li、K、Cl的力场参数是描述水中离子状态（而且往往还只 ...

谢谢SOB老师，我也看到CP2K支持很多熔盐相关力场。我在学习动力学课程之后受益匪浅，想再学习一下CP2K培训班，但是不知道是否需要先学习初级量化，还是直接学习CP2K就行？（我没有高斯这些软件的基础）

Streicheln · 发表于 Post on yesterday 21:29

liuyuje714 发表于 2026-5-9 20:30
低版本gmx （gmx

谢谢老师，我找到了一个WBK force，他也是要求GROMACS 4.6.7 double precision这个老版本，我打算试试用lampps或者cp2k来计算。

sobereva · 发表于 Post on yesterday 23:06

Streicheln 发表于 2026-5-9 21:25
谢谢SOB老师，我也看到CP2K支持很多熔盐相关力场。我在学习动力学课程之后受益匪浅，想再学习一下CP2K培 ...

必须先有初级量化班程度的基础

		自动登录 Automatic login	找回密码 Forget password
密码 Password			注册 Register

[GROMACS] GMX模拟熔盐时是否有合适的势函数呢？