阿夸夸 发表于 2026-2-15 21:25
方便的话 可以把输入文件上传一下，根据你的描述“加UKS且注释掉Multiplicity”感觉输入文件写的可能有问题

您好老师
比如我试过这样写：

1. &DFT
   
2.     BASIS_SET_FILE_NAME  BASIS_MOLOPT
   
3.     POTENTIAL_FILE_NAME  POTENTIAL
   
4. #   WFN_RESTART_FILE_NAME cp2k ads slab-RESTART.wfn
   
5.     CHARGE    0 #Net charge
   
6.     MULTIPLICITY    1 #Spin multiplicity

复制代码

也试过：

1. &DFT
   
2.     BASIS_SET_FILE_NAME  BASIS_MOLOPT
   
3.     POTENTIAL_FILE_NAME  POTENTIAL
   
4. #   WFN_RESTART_FILE_NAME cp2k ads slab-RESTART.wfn
   
5.     CHARGE    0 #Net charge
   
6.     UKS
   
7.     MULTIPLICITY    5 #Spin multiplicity

复制代码

也试过注释掉或者直接删掉自旋多重度：

1. &DFT
   
2.     BASIS_SET_FILE_NAME  BASIS_MOLOPT
   
3.     POTENTIAL_FILE_NAME  POTENTIAL
   
4. #   WFN_RESTART_FILE_NAME cp2k ads slab-RESTART.wfn
   
5.     CHARGE    0 #Net charge
   
6.     UKS
   
7.     #MULTIPLICITY    1 #Spin multiplicity

复制代码

northern_mar 发表于 2026-2-15 20:02
您好！有个问题想请教一下您。
您说的第三点“3.Orca计算时使用PBE-D3(BJ)做单点能任务进行自旋多重度扫描 ...

自旋多重度=2S+1=alpha电子-beta电子+在inp输入文件里面设置的，比如分子体系可能有一行是* xyz   0   1，0代表体系电荷，1代表体系的自旋多重度。对于很多分子体系而言一般是闭壳层的，没有未成对电子，所以自旋多重度=2*0+1=1，也叫单重态，开壳层体系的自旋多重度得进行测试（我是这样认为的，文献里面也是这么做的，如果有更好的做法请指教），比如奇数个电子就测试二重态、四重态等等，偶数个电子如果是开壳层就测试三重态、五重态、七重态等等，对于含有过渡金属的开壳层体系，最高的可能的自旋多重度应当为过渡金属所有电子都不成对的情况（我记得社长在一个帖子里面是这样回复的，具体的表述记不太清楚了你可以搜一下），然后从众多测试中挑选出能量最低的那一个自旋多重度，只不过换了一个词，把逐一测试并挑选能量最低的自旋多重度的过程称为了“扫描”，本质还是量子化学计算，没有专门的软件，写个脚本在传算上一次性上传多个任务就行

方便的话 可以把输入文件上传一下，根据你的描述“加UKS且注释掉Multiplicity”感觉输入文件写的可能有问题

您好！有个问题想请教一下您。
您说的第三点“3.Orca计算时使用PBE-D3(BJ)做单点能任务进行自旋多重度扫描，发现Spin 5时能量最低”是如何实现的？我也有个体系，不知道自旋多重度该设怎么设，如果有办法直接用程序进行扫描那可太好了。现在只能一个一个尝试比较能量大小。