| 查看: 1119 | 回复: 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 |
回复此楼» 猜你喜欢
泣血叩求调剂恩，愿以丹心报师恩
已经有3人回复
-
274求调剂0856材料化工
已经有6人回复
-
0856化工原理
已经有3人回复
-
291分兽医专硕 一志愿华中农
已经有3人回复
-
工科0856专硕化学工程269能调剂吗
已经有8人回复
-
341求调剂
已经有3人回复
-
接受26届调剂生
已经有23人回复
-
311求调剂
已经有5人回复
-
308求调剂
已经有3人回复
-
271求调剂
已经有3人回复
2楼2022-07-20 13:43:42
