| 查看: 1121 | 回复: 1 | ||
| 【悬赏金币】回答本帖问题，作者mpiyq将赠送您 5 个金币 | ||
[求助]
模拟水蒸气凝结为液滴的表面张力
|
||
|
希望控制压力为0，但是fix 3 all press/berendsen iso 0.0 0.0 1000.0 这个命令一直跑不下去，希望有大神指导一下！ 下面是我的凝结成为液滴部分in文件 # in file for water molecule SPC/E model echo screen units real atom_style full boundary p p p bond_style harmonic angle_style harmonic pair_style lj/charmm/coul/long 9.0 10.0 10.0 pair_modify mix arithmetic tail yes kspace_style ewald 0.0001 #创建原子 lattice sc 4 region box block 0 10 0 10 0 10 create_box 2 box bond/types 1 angle/types 1 extra/bond/per/atom 2 extra/angle/per/atom 2 extra/special/per/atom 2 molecule SPCE mol.SPCE create_atoms 0 box mol SPCE 10 units lattice #相互作用 pair_coeff 1 1 0.1553 3.166 pair_coeff * 2 0.000 0.000 bond_coeff 1 1000 1.0 angle_coeff 1 100 109.47 #mass mass 1 15.9994 mass 2 1.008 #MD setting neighbor 2.0 bin neigh_modify delay 0 every 2 check yes timestep 1.0 thermo 100 thermo_style custom step temp pe ke press pxx pyy pzz vol lx ly lz minimize 1e-4 1e-6 10000 10000 reset_timestep 0 #预平衡 velocity all create 200.0 54654 fix 1 all shake 0.0001 50 0 b 1 a 1 mol SPCE dump Dump all custom 5000 NVT.lammpstrj id type x y z ## 稳定压力在0，z方向不控压，盒子尺寸不变，跑足够长时间，达到压力和温度平衡 fix 3 all press/berendsen iso 0.0 0.0 1000.0 run 10000 unfix 3 fix 4 all nvt temp 300 300 100.0 run 10000 reset_timestep 0 |
回复此楼» 猜你喜欢
NSFC申报书里申请人简历中代表性论著还需要在申报书最后的附件里面再上传一遍吗
已经有14人回复
-
材料与化工一志愿南昌大学327求调剂推荐
已经有6人回复
-
化学调剂0703
已经有7人回复
-
327求调剂
已经有11人回复
-
调剂
已经有8人回复
-
梁成伟老师课题组欢迎你的加入
已经有7人回复
-
伙伴们，祝我生日快乐吧
已经有24人回复
-
中科院材料273求调剂
已经有3人回复
-
材料工程专硕274一志愿211求调剂
已经有5人回复
-
本人考085602 化学工程 专硕
已经有8人回复
2楼2022-07-20 13:43:42
