| 查看: 3269 | 回复: 1 | ||
[求助]
用lammps作PMF不平衡
|
|
各位大神，我尝试用LAMMPS中fix smd命令做PMF。体系是中间为石墨烯，是石墨烯左右两边是NaCl溶液。准备分别做钠离子、氯离子和水的PMF.但运行的结果，却是左右不对称。请各位大神帮我看看出现什么问题了。谢谢！！！ units real atom_style full neighbor 2.0 bin boundary p p p read_data H-PMF.data #lj potentials pair_style lj/cut/coul/long 10.0 pair_coeff 1 1 0.07 3.55 pair_coeff 2 2 0.07 3.55 pair_coeff 3 3 0.03 2.42 pair_coeff 4 4 0.1553 3.166 pair_coeff 5 5 0 0 pair_coeff 6 6 0.1 2.583 pair_coeff 7 7 0.1 4.401 pair_coeff 8 8 0.1553 3.166 pair_coeff 9 9 0 0 pair_coeff 10 10 0.1 2.583 pair_coeff 11 11 0.1 4.401 pair_coeff 12 12 0.1553 3.166 pair_coeff 13 13 0 0 pair_modify mix arithmetic #special_bonds charmm bond_style harmonic bond_coeff 1 469 1.40 bond_coeff 2 367 1.080 bond_coeff 3 450 1.0 bond_coeff 4 450 1.0 angle_style harmonic angle_coeff 1 63 120 angle_coeff 2 35 120 angle_coeff 3 55 109.47 angle_coeff 4 55 109.47 #special_bonds lj/coul 0.0 0.0 0.5 neigh_modify delay 0 every 1 check yes kspace_style pppm 1.0e-4 #kspace_modify slab 3.0 special_bonds lj 0.0 0.0 0.5 coul 0.0 0.0 0.5 #define groups group GHE type 1 2 3 group shui1 type 4 5 12 13 group shui2 type 8 9 group LJ type 6 7 10 11 group shui union shui1 shui2 group and union shui LJ group smd-w type 12 13 neigh_modify exclude group GHE GHE #neigh_modify exclude type 9 9 #neigh_modify delay 0 every 1 #minimize 1.0e-4 1.0e-6 100 1000 #velocity and create 298.0 4928459 units box velocity GHE set 0.0 0.0 0.0 units box fix 1 GHE setforce 0 0 0 fix 2 and nvt temp 298.0 298.0 100.0 fix 3 GHE rigid/nvt group 1 GHE temp 298.0 298.0 100.0 force * off off off torque * off off off fix 4 smd-w rigid group 1 smd-w force * off off on torque * off off off fix 5 shui shake 0.0001 20 0 b 3 4 a 3 4 fix 6 smd-w smd cvel 20 -0.000005 tether NULL NULL 15 0.0 |
回复此楼» 猜你喜欢
基金中了
已经有14人回复
-
关于filecode
已经有7人回复
-
娱乐
已经有6人回复
-
化学口download_prp&fileCode的固定段好像这几天一直没变，有变的大神么？
已经有4人回复
-
关于代码变化问题，想知道的进来
已经有12人回复
-
固定端突然变了，今天
已经有10人回复
-
听说今天filecode变了
已经有47人回复
-
filecode与中标关系的预测
已经有3人回复
-
filecode变化情况
已经有22人回复
-
十年后又回来了，论文投稿求助
已经有4人回复
» 本主题相关商家推荐: (我也要在这里推广)
» 本主题相关价值贴推荐，对您同样有帮助:
- Lammps对单晶铜进行弛豫，体积忽然增大，内存异常，vmd查看发现体积基本不变
已经有13人回复
- lammps 中 做拉伸时 sss和ssp边界条件不同，结果为什么不一样？？？？
已经有5人回复
- 求助LAMMPS两次运算结果不一致问题！
已经有11人回复
- MS生成Lammps的data文件，文件中原子的类型和顺序并不是自己想要的，请问如何修改？
已经有8人回复
- MS里的forcite 模块只能用自带的那三个力场吗？想用trappe力场是不是要用lammps?
已经有6人回复
- 分子模拟中体系势能和总能量为很大的正值
已经有8人回复
- 分子动力学能量守恒的相关问题
已经有18人回复
- 想做粒子和固体（金属）相互作用，不知道lammps能否实现？
已经有6人回复
- lammps实验拉伸Al纳米线的in文件的错误提示好久找不出？
已经有8人回复
- lammps 离子液体，离子中相邻隔的挺远的原子成键了，不知什么原因
已经有8人回复
- LAMMPS 的坐标文件怎么构建呢？不可能一个一个原子的输入吧？
已经有8人回复
- 用MS建模导入LAMMPS计算，是不是不能用compass力场？
已经有9人回复
- lammps求助：计算含有不同原子体系的热导率的in文件该怎么修改
已经有7人回复
- 【求助】NVT系综在MD模拟过程中的能量变化
已经有6人回复
- 【求助】LAMMPS中非周期性边界条件和Born-Mayer-Huggins势函数不会设置
已经有4人回复
2楼2017-06-09 16:00:21
100