标题: 用cp2k做结构优化和频率计算的一些问题 [打印本页] 作者Author: 轻舟~ 时间: 2025-3-18 11:33 标题: 用cp2k做结构优化和频率计算的一些问题 学生在做2-吡啶甲酸的结构优化和频率计算时，遇到许多问题。比如在CCDC里下载2-吡啶甲酸的晶胞模型（.cif文件，CCDC号130513），之后导入Multiwfn中，会提示 Warning: At least one atom has occupancy smaller than 1! All atoms will be loaded regardless of their occupancies，忽略警告，输入cp2k，并命名一个.inp文件，然后选择-1， 4 Optimizing both structure and cell，（1，默认PBE），2， 11 6-311G**，3， 11 6-311G**，4，开OT，0，输出文件，之后.inp文件里&DFT里面的内容中有这么一条： MULTIPLICITY 2 #Spin multiplicity，后面的计算报错就跟此有关，我把.inp文件放到服务器中运行，开始还能计算，而且算出来了，但是把.restart文件继续投入Multiwfn中，构造一个2*2*2超胞后，使用一样的基组和方法以及色散，开启OT，注意到生成的.inp文件中&DFT中还是有MULTIPLICITY 2 #Spin multiplicity，投入服务器计算后不久直接报错，显示自选多重度有问题。于是，我采用了一个幼稚的方法，就是自行修改这个超胞文件里的自选多重度，但是依然报错。之后我认为是不是晶报文件是不是有问题，于是就从COD下载晶胞文件，COD号2007754，采用同样的做法，从2025年3月11日开始进行原胞的结构优化，直到现在依然没有算出来，学生不知道这个过程是不是有问题,现在依次把CCDC版本的.inp文件（1、js.inp），原胞输出文件（2、js.out）,超胞.inp文件（3、jscb.inp）,超胞.out文件（4、jscb.out）;COD版本的.inp文件（5、js(COD产生).inp），COD版本的.out文件（内容太多只粘贴前两轮计算，命名为6、CODjs.out）放在帖子里，求大家解惑。