| 查看: 1202 | 回复: 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 |
回复此楼» 猜你喜欢
一个有机合成实验室都需要哪些设备？
已经有6人回复
-
2026年国自然面上资助率
已经有20人回复
-
面上再次挂了，太难了，躺也躺不了，倦也卷不过，小学校之殇！
已经有22人回复
-
微信指数没变化，科研之友没阅读
已经有18人回复
-
基础研究怎么拉横向，学校到款任务越来越多，难以完成 拉横向，都有哪些途径啊
已经有8人回复
-
HXDI做水性聚氨酯乳液，是不是特别容易出渣
已经有3人回复
-
系统今天又提示维护了，估计离放榜不远了
已经有15人回复
-
你们的时间戳变了吗
已经有4人回复
-
产物和副产物价值比较
已经有5人回复
-
时间戳他又来了
已经有17人回复
2楼2022-07-20 13:43:42