| 查看: 1146 | 回复: 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 |
回复此楼» 猜你喜欢
一志愿华中科技大学071000，求调剂
已经有4人回复
-
生物学一志愿985，分数349求调剂
已经有4人回复
-
生物学调剂
已经有4人回复
-
求调剂院校信息
已经有4人回复
-
085600材料与化工306
已经有4人回复
-
286求调剂
已经有10人回复
-
328求调剂，英语六级551，有科研经历
已经有9人回复
-
一志愿北京化工大学070300 学硕336求调剂
已经有4人回复
-
286分人工智能专业请求调剂愿意跨考！
已经有8人回复
-
资源与环境 调剂申请(333分)
已经有5人回复
» 本主题相关商家推荐: (我也要在这里推广)
