| 查看: 1913 | 回复: 2 | ||
[求助]
高分子中二氧化硅夹杂模拟 bond atoms missing
|
|
请教一下大家 模拟高分子体系里有纳米sio2。 是在ms里建模pcff势能，然后导出在lammmps里算，class2势能。 之前只算高分子都是没问题的。 加了sio2之后，先nve没问题，但是npt时运行一段时间（几万到几十万）总是报错： bond atoms x x missing 或者out of range atoms cannot compute pppm。 看了一下一般都是sio2表面的氢原子丢失了。 看轨迹系统也比较稳定没有炸开。 改小步长只能延迟报错时间。 不知道可能是什么原因呢  哦，MS里建模SiO2是导入自带的structure的glasses的SiO2_21A_3d。不知道是不是这里有问题 下面是文件： #################################################### # Initialization variable fname index R1L4 variable tem equal 500 log ${fname}.log units real atom_style full boundary p p p #################################################### # Forcefield parameters pair_style lj/class2/coul/long 9.5 9.5 pair_modify tail yes bond_style class2 angle_style class2 dihedral_style class2 improper_style class2 kspace_style pppm 5.0e-6 neigh_modify delay 0 every 1 check yes dielectric 1.0 special_bonds lj/coul 0.0 0.0 1.0 comm_modify cutoff 20.0 #***************************试着改的 read_data ${fname}.data #################################################### # Equilibriation step 1 - Langevin dynamics at 500K thermo_style custom step temp press ke pe vol #thermo_modify lost/bond ignore #**********************************试着改的 thermo 100 timestep 0.5 #**********************************试着改的 fix 1 all nve/limit 0.4 velocity all create ${tem} 1234 fix 2 all temp/rescale 1 ${tem} ${tem} 1 1.0 fix 3 all momentum 1 linear 1 1 1 angular dump dump0 all atom 1000 all0_*.atom run 10000 unfix 1 unfix 2 write_data ${fname}_end1.txt ################################################### min_style sd #***************************试着改的 minimize 1.0e-5 1.0e-7 10000 100000 #################################################### fix 1 all npt temp ${tem} ${tem} 100 iso 0 0 1000 run 800000 unfix 1 undump dump0 write_data ${fname}_end2.txt 一般运行到这里就报错了“run 800000 ”能运行几万到十几万步 |
回复此楼» 猜你喜欢
浙大物院表面物理课题组长期招收普博及直博学生
已经有259人回复
-
【征稿】Springer Nature热电专题 Discover Materials / Discover Chemistry
已经有47人回复
-
物理学I论文润色/翻译怎么收费?
已经有121人回复
-
25届硕士申请，材料计算方向
已经有2人回复
-
课题组招1-2名硕士
已经有1人回复
-
机时越来越贵之后，我把一部分 DFT 流程搬到了 AI 计算平台上
已经有0人回复
-
斯坦福大学与Elsevier世界前2%学者新排名查询
已经有1人回复
-
27申博推荐：南京航空航天大学国际前沿院光学方向招收博士生！
已经有19人回复
-
华南师范大学+双一流学科物理学+广东省杰青团队2027年博士招生+机器学习驱动材料设计
已经有1人回复
-
apl投稿状态
已经有2人回复
-
华南师范大学+双一流学科物理学+广东省杰青团队2027年博士招生+机器学习驱动材料设计
已经有0人回复
» 本主题相关商家推荐: (我也要在这里推广)
2楼2018-12-13 16:20:54
awuli张
新虫 (小有名气)
- 应助: 0 (幼儿园)
- 金币: 253
- 散金: 30
- 帖子: 97
- 在线: 22.5小时
- 虫号: 15793324
- 注册: 2019-07-31
- 专业: 化学工程及工业化学
3楼2019-09-27 00:45:12
50