[交流] cp2k对NaCl的结构优化问题 已有1人参与

新手使用，求大神们指教，再用cp2k进行NaCl的结构优化过程中出现的问题： 1，Nonbonded section中的 RCUT 好像对结果没什么影响，因为我设定从2.0到20.0之间的数值，结果都相同。 2，输出结果中有个xyz文件，打开后发现面上的键长和键角发生了变化（尤其是键角发生了明显变化180->178度），虽然code里面设定了keep symmetry 和keep angle，但貌似没有什么作用。 3，怎样在输出的out 文件中直接读取键长键角，是否需要加入新的print函数，应该在什么地方加入呢？ 下面是code和xyz输出图片的具体： @SET RCUT_MY 2.0 @SET SPLINE_MY 0.50 @SET ALPHA_MY 0.35 @SET OPTIMIZER LBFGS # BFGS &FORCE_EVAL   METHOD Fist   &MM     &FORCEFIELD       &CHARGE         ATOM Na         CHARGE +1.000       &END CHARGE       &CHARGE         ATOM Cl         CHARGE -1.000       &END CHARGE       &NONBONDED         &BMHFT           map_atoms NA NA           atoms NA NA           RCUT ${RCUT_MY}         &END BMHFT         &BMHFT           map_atoms NA CL           atoms NA CL           RCUT ${RCUT_MY}         &END BMHFT         &BMHFT           map_atoms CL CL           atoms CL CL           RCUT ${RCUT_MY}         &END BMHFT       &END NONBONDED       &SPLINE         EMAX_SPLINE ${SPLINE_MY}       &END SPLINE     &END FORCEFIELD     &POISSON       &EWALD         EWALD_TYPE spme         ALPHA ${ALPHA_MY}         GMAX 12         O_SPLINE 6       &END EWALD       PERIODIC XYZ       POISSON_SOLVER PERIODIC     &END POISSON   &END MM   &SUBSYS     &CELL       ABC 12*5.620 12*5.620 12*5.620       MULTIPLE_UNIT_CELL 1 1 1     &END CELL     &TOPOLOGY       #&CNETER_COORDINATES .TRUE.       #&END       COORD_FILE_NAME NaCl27.pdb       COORDINATE PDB       CONN_FILE_FORMAT OFF       MULTIPLE_UNIT_CELL 1 1 1     &END TOPOLOGY     &PRINT       &CELL       &END CELL       &TOPOLOGY_INFO         XYZ_INFO .TRUE.       &END TOPOLOGY_INFO     &END PRINT   &END SUBSYS   &PROPERTIES     &RESP     &END RESP   &END PROPERTIES &END FORCE_EVAL &GLOBAL   PROJECT NaCl27_1   RUN_TYPE GEO_OPT   PRINT_LEVEL MEDIUM &END GLOBAL &MOTION   &GEO_OPT     OPTIMIZER ${OPTIMIZER}   &END   &CELL_OPT     KEEP_ANGLES .TRUE.     KEEP_SYMMETRY .TRUE.   &END &END MOTION Screenshot_20160810_103519.png

回复此楼