| 查看: 1243 | 回复: 0 | |||
| 【悬赏金币】回答本帖问题，作者懒羊羊要当大厨将赠送您 5 个金币 | |||
[求助]
LAMMPS计算外加电场作用下reaxff
|
|||
|
各位老师好，我在计算微波作用下焦炭CO2气化反应时，根据参考文献通过外加电场来模拟微波的作用。令人沮丧的是，在外加电场作用下，反应不仅没有得到促进，反而呈现抑制作用，这与文献计算结果和实验结果相反。文献使用的计算软件是ADS，在偶然的机会下，我发现在在ADS官方网站上有这样一句话“LAMMPS计算外加电场作用下反应力场计算，仅考虑了电场的力的作用，而没有考虑电场能量的作用（https://v.fermitech.com.cn/wiki/doku.php?id=adf%3Aapplywallpotential&sid_for_share=99125_4）”。这确实有可能造成我的计算结果与实验和文献相反，希望各位老师能为我指点。下面是我计算所使用的in文件 units real dimension 3 boundary p p f neighbor 2 bin neigh_modify every 10 delay 0 check yes timestep 0.1 atom_style full pair_style reaxff NULL safezone 5.0 mincap 600 minhbonds 1000 variable T equal 3000.0 variable r_p string "result/0.04" variable Ez equal 0.04 log ${r_p}/logout.lammps read_data H2O_c.data fix zwalls all wall/reflect zlo EDGE zhi EDGE pair_coeff * * ffield.reax.CHON C N O H fix 1 all qeq/reaxff 1 0.0 10.0 1.0e-6 reaxff min_style sd minimize 1e-8 1e-8 10000 10000 velocity all create 300.0 4928459 rot yes dist gaussian reset_timestep 0 thermo 200 thermo_style custom step temp pe ke v_Ez thermo_modify flush yes fix 2 all nvt temp 300.0 300.0 10 fix 3 all reaxff/species 1 1 200 ${r_p}/NVT.out element C N O H cutoff 3 3 0.5 fix 4 all reaxff/bonds 200 ${r_p}/bonds_nvt.txt dump 1 all xyz 200 ${r_p}/NVT.lammpstrj dump_modify 1 element C N O H run 200000 unfix 2 unfix 3 unfix 4 undump 1 reset_timestep 0 fix 5 all nvt temp 300 ${T} 10.0 fix 6 all reaxff/species 1 1 200 ${r_p}/heat.out element C N O H cutoff 3 3 0.5 fix 7 all reaxff/bonds 200 ${r_p}/bonds_heat.txt run 200000 unfix 5 unfix 6 unfix 7 reset_timestep 0 fix e all efield 0.0 0.0 v_Ez fix t all nve fix 8 all temp/berendsen ${T} ${T} 1000 dump 2 all xyz 200 ${r_p}/gas.lammpstrj dump_modify 2 element C N O H fix 9 all reaxff/species 1 1 200 ${r_p}/gas.out element C N O H cutoff 3 3 0.5 fix 10 all reaxff/bonds 200 ${r_p}/bonds_gas.txt run 5000000 |
回复此楼» 猜你喜欢
2025冷门绝学什么时候出结果
已经有3人回复
-
天津工业大学郑柳春团队欢迎化学化工、高分子化学或有机合成方向的博士生和硕士生加入
已经有4人回复
-
康复大学泰山学者周祺惠团队招收博士研究生
已经有6人回复
-
AI论文写作工具：是科研加速器还是学术作弊器？
已经有3人回复
-
孩子确诊有中度注意力缺陷
已经有6人回复
-
2026博士申请-功能高分子，水凝胶方向
已经有6人回复
-
论文投稿，期刊推荐
已经有4人回复
-
硕士和导师闹得不愉快
已经有13人回复
-
请问2026国家基金面上项目会启动申2停1吗
已经有5人回复
-
同一篇文章，用不同账号投稿对编辑决定是否送审有没有影响？
已经有3人回复
