使用CP2K培训班里面的Fe.cif，用Multiwfn扩胞后，用gv删除一部分，建立仅有两层铁原子的团簇，且铁原子数量为偶数个（我知道应该至少建立三层，但这里为了图快，就只有两层），在表面放了一个甲烷，冻结底部和四周的铁原子进行优化。
第一个问题：如果不知道体系开壳层闭壳层的话，那么无脑就把金属团簇的原子总数量弄成偶数个，那么该团簇就必然是闭壳层，会省很多事？
第二个问题：直接建立的铁胞体后，有没有必要再进行切割晶胞，将甲烷放在不同的切割面的表面？？为什么要切割晶胞？如果不切，直接使用任何一个面作为研究的表面，又会怎么样？如果把甲烷放的距离铁表面比较近或者远，优化都会进行不下去，然而靠化学直觉，把甲烷放在一个合适的位置，则会比较顺利收敛。这个“化学直觉”，对量化新手比较不友好。 (, 下载次数 Times of downloads: 10)
上传 Uploaded
点击下载Click to download
使用的输入文件如下：
! PBE0 D3 def2-SV(P) def2/J RIJCOSX opt noautostart miniprint nopop
%maxcore 1000
%pal nprocs 6 end
%geom
Constraints
省略***********
end
end
* xyz 0 1 #第三个问题：在CP2K的培训班里，讲到了bcc铁自旋多重度，设了也不起作用，于是我这里也就只设为1，这样究竟合适吗？如果不合适，那么该设置为多少呢？有些文献中，认为单质铁的自旋多重度应该是5，有些文献甚至把铁团簇的自旋多重度设为几十上百都有。如果铁表面放了一个O2分子，那么自旋多重度应该设为3吗？
省略****************
完整的输入文件在此： (, 下载次数 Times of downloads: 1)
上传 Uploaded
点击下载Click to download
第二个步骤
复制代码
2. CPU是13900K，用了6核，优化了35个小时，终于结束，但是没有冻结的铁原子向内部凹陷了。。
第四个问题：按理说，冻结的铁原子当然是稳如泰山，没有冻结的铁原子也应该几乎纹丝不动，但这里明显凹陷了，是因为PBE0不合适吗？如果PBE0不合适优化铁团簇，那用什么级别会更好呢？或者因为团簇不够大吗？或者别的什么原因呢？ (, 下载次数 Times of downloads: 9)
*****************************************************
* ERROR *
* SCF NOT CONVERGED AFTER 1230 CYCLES *
*****************************************************
Last Energy change ... 0.0000e+00 Tolerance : 1.0000e-06
Last Orbital Gradient ... 1.8242e-02 Tolerance : 5.0000e-05
Last Orbital Rotation ... 1.9619e-01
--------------------------------------------------------------------
ERROR
This wavefunction IS NOT CONVERGED!
And it is 'too far' from meeting the convergence thresholds.
Aborting ...
Please restart calculation (with larger maxiter?)
--------------------------------------------------------------------
ORCA finished by error termination in SCF
Calling Command: mpirun -np 6 /home/jing/software/orca504/orca_scf_mpi Fe.gbw b Fe
[file orca_tools/qcmsg.cpp, line 465]:
.... aborting the run