| 查看: 1741 | 回复: 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 |
回复此楼» 猜你喜欢
售SCI-T0P文章，我:8O.5.5.1.O.54,科目齐全,可+急
已经有4人回复
-
售SCI文章，我:8O.5.5.1O.54,科目全,可十急
已经有5人回复
-
售SCI一区T0P文章，我:8.O55.1.O.54,科目全,可十急
已经有5人回复
-
售SCI文章，我:8O5.5.1.O.54,科目齐全,可+急
已经有6人回复
-
售SCI一区T0P文章，我:8.O55.1.O.54,科目全,可十急
已经有4人回复
-
售SCI文章，我:8O5.5.1.O.54,科目齐全,可+急
已经有4人回复
-
小木虫上这么多卖论文的，真有人买论文么？感觉没必要啊
已经有12人回复
-
售SCI一区文章，我:8O5.5.1.O5.4,科目全,可伽急
已经有4人回复
-
各位道友，我要去昆明玩几天，回来见。
已经有8人回复
-
售SCI一区T0P文章，我:8O.55.1.O.5.4,科目齐全,可+急
已经有4人回复
