| 查看: 1215 | 回复: 0 | ||
| 【悬赏金币】回答本帖问题，作者周三啊将赠送您 10 个金币 | ||
周三啊铁虫 (小有名气)
|
[求助]
关于recenter命令
|
|
请问各位大佬，我在一篇文献中看到关于CNT吸附在石墨烯上的结构，文献中提到集合优化时在NVE条件下用minimizer命令，之后在X方向移动CNT，并保持CNT的质心保持在一定的X值，测出CNT在X方向上的sliding potential energy traces 。我按照他文献的过程，先进行minimizer，之后用fix recenter，发现结果是先执行minimizer，此时的CNT的质心已经变化了，请问是我的in script哪里出错了 # Initialization units metal boundary p p p atom_style full processors * * 1 # domain decomposition over x and y # read lammps data file read_data CNT66_f_3.1_AB.data mass 1 12.0107 # carbon mass (g/mole) | membrane mass 2 12.0107 # carbon mass (g/mole) | adsorbate # Separate atom groups group membrane type 1 group adsorbate type 2 ######################## Potential defition ######################## pair_style hybrid/overlay airebo 3.0 ilp/graphene/hbn 16.0 #################################################################### pair_coeff * * airebo CH.airebo NULL C # chemical pair_coeff * * ilp/graphene/hbn BNCH.ILP C C # long range #################################################################### # Neighbor update settings neighbor 2.0 bin neigh_modify every 1 neigh_modify delay 0 neigh_modify check yes ##timestep timestep 0.001 compute cntcom adsorbate com compute gracom membrane com variable cntcomx equal c_cntcom[1] variable cntcomy equal c_cntcom[2] variable cntcomz equal c_cntcom[3] variable gracomx equal c_gracom[1] variable gracomy equal c_gracom[2] variable gracomz equal c_gracom[3] ##minimize thermo 100 velocity adsorbate create 0.1 12345 fix subf membrane setforce 0.0 0.0 0.0 fix thermostat all nve fix fi_com membrane recenter INIT NULL NULL min_style fire minimize 1.0e-6 1.0e-6 1000 1000 thermo 100 thermo_style custom step etotal pe ke v_cntcomx v_cntcomy v_cntcomz v_gracomx v_gracomy v_gracomz thermo_modify line one format float %.10f thermo_modify flush yes norm no lost warn dump 11 all atom 100 cnt_gra.lammpstrj run 1000  HIU8)FDYZPW~M%(I6TL2523.png |
回复此楼» 猜你喜欢
售SCI一区T0P文章，我:8.O.55.1.O.5.4,科目全,可+急
已经有5人回复
-
售SCI一区T0P文章，我:8O.55.1.O.5.4,科目齐全,可+急
已经有6人回复
-
售SCI一区T0P文章，我:8.O.55.1.O.54,科目齐全,可+急
已经有4人回复
-
售SCI一区T0P文章，我:8.O.55.1.O.5.4,科目全,可+急
已经有6人回复
-
售SCI-T0P文章，我:8O.5.5.1.O.54,科目齐全,可+急
已经有4人回复
-
售一区SCI文章T0P，我:8O.551.O54,科目全,可十急
已经有5人回复
-
售SCI文章，我:8O.5.5.1O.54,科目全,可十急
已经有7人回复
-
售SCI一区T0P文章，我:8.O.55.1.O.5.4,科目全,可+急
已经有5人回复
-
售SCI文章，我:8O.5.5.1O.54,科目全,可十急
已经有8人回复
-
售一区SCI文章T0P，我:8O.551.O54,科目全,可十急
已经有5人回复
» 本主题相关商家推荐: (我也要在这里推广)
40