| 查看: 1118 | 回复: 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 |
回复此楼» 猜你喜欢
一志愿北京化工大学材料与化工296分求调剂
已经有15人回复
-
深圳大学硕士招生（2026秋，传感器方向）
已经有6人回复
-
290求调剂
已经有5人回复
-
一志愿河海大学085900土木水利专硕279求调剂不挑专业
已经有3人回复
-
0857环境调剂
已经有5人回复
-
一志愿浙江大学0856材料与化工求调剂
已经有4人回复
-
081200计算机科学与技术数一英一306求调剂
已经有3人回复
-
材料工程330分求调剂，一志愿985
已经有4人回复
-
求调剂
已经有4人回复
-
材料371求调剂
已经有9人回复
2楼2022-07-20 13:43:42
