计算化学公社

标题: 使用lammps模拟TIP4P水溶液求助 [打印本页]

作者
Author: zhangsw17 时间: 2021-4-13 17:44
标题: 使用lammps模拟TIP4P水溶液求助
各位好，我最近刚开始学习分子动力学模拟，想要模拟盐溶液在固体表面的吸附情况，参考文献中用到了TIP4P刚性水分子，我用MS建立了键角104.52度，键长0.9572A的水分子，然后用AC盒子构建了包含水分子和K+，Cl-的溶液体系，用build layer将溶液放在固体表面，然后导入lammps计算，我看了很久的手册，还是有一些问题想要请教大家
1）文献中对于这个体系用到的势函数类似于u=Σqiqj/r+Aexp(-Br)-C/r^4+D/r^6，如果想在lammps中实现是不是要用hybrid/overlay来拼凑出文献中的势函数形式？lammps中有这种LJ6-4形式的势函数吗，或者是单独的C/r^4和C/r^6这样的势函数，我只找到了buck类型是Aexp(-Br)-C/r^6。
2）导入MS的模型之后是否需要能量最小化minimize之后再计算？这一过程是不是会改变水分子的构型，手册中提到TIP4P模型用fix shake来固定水分子的键长键角，这一命令不能在minimize之前使用。

作者
Author: sobereva 时间: 2021-4-14 06:25
TIP4P是很过时的水模型。同为四点水模型的OPC完爆TIP4P。如今还用TIP4P可能会遭到内行审稿人吐槽

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