| 查看: 1120 | 回复: 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 |
回复此楼» 猜你喜欢
伙伴们，祝我生日快乐吧
已经有24人回复
-
中科院材料273求调剂
已经有3人回复
-
材料工程专硕274一志愿211求调剂
已经有5人回复
-
本人考085602 化学工程 专硕
已经有8人回复
-
求老师收留调剂
已经有5人回复
-
机械专硕调剂
已经有3人回复
-
东华理工大学化材专业26届硕士博士申请
已经有6人回复
-
欢迎申博同学联系
已经有7人回复
-
东南大学364求调剂
已经有3人回复
-
326求调剂
已经有4人回复
2楼2022-07-20 13:43:42
