| 查看: 1152 | 回复: 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 |
回复此楼» 猜你喜欢
化工学硕294分，求导师收留
已经有29人回复
-
求调剂
已经有10人回复
-
考研求调剂
已经有13人回复
-
求调剂
已经有3人回复
-
人工智能320调剂08工类还有机会吗
已经有17人回复
-
考研英一数一338分
已经有10人回复
-
求助调剂，跨调
已经有15人回复
-
085600材料与化工329分求调剂
已经有20人回复
-
085600材料与化工349分求调剂
已经有15人回复
-
求调剂
已经有13人回复
» 本主题相关商家推荐: (我也要在这里推广)
