含过渡金属体系，更换大基组，结构自旋出现问题

LYFE

含过渡金属体系，使用ub3lyp/gen，过渡金属用赝式基组，其他原子使用6-31g**优化出的结构自旋呈现开壳层单重态，更换大基组ub3lyp/cc-pvtz优化，读取小基组优化结果的电子信息，但是优化结果显示原子上没有自旋分布。想请教一下这是出现了什么问题，怎么解决。

sobereva

使用赝势基组应当写genecp而非gen

倘若你确保关键词没写错，大基组下做单点时带着stable关键词做波函数稳定性测试

zjxitcc

很大概率是输入文件写错了、关键词不合理。建议晒晒输入文件（坐标可以不贴，其他内容原样复制，且用论坛上发帖功能中的代码块显示）

LYFE

1. %chk=TS-BS.chk
   
2. %nprocshared=24
   
3. %mem=3GB
   
4. #p ub3lyp/cc-pvtz nosymm pseudo=read scf=xqc int=ultrafine IOp(3/76=1000001500) IOp(3/77=0720008500) IOp(3/78=0810010000)
   
5. 
   
6. Title Card Required

复制代码

LYFE

zjxitcc 发表于 2023-5-16 10:10
很大概率是输入文件写错了、关键词不合理。建议晒晒输入文件（坐标可以不贴，其他内容原样复制，且用论坛上 ...

1. %chk=TS-BS.chk
   
2. %nprocshared=24
   
3. %mem=3GB
   
4. #p ub3lyp/cc-pvtz nosymm pseudo=read scf=xqc int=ultrafine IOp(3/76=1000001500) IOp(3/77=0720008500) IOp(3/78=0810010000)
   
5. 
   
6. Title Card Required

复制代码

zjxitcc

LYFE 发表于 2023-5-18 09:51

cc-pVTZ对Cu元素有定义，且是全电子基组，所以不应该写“ub3lyp/cc-pvtz pseudo=read”，坐标底下也不应该写Cu的赝势信息；只写ub3lyp/cc-pvtz就行。

基组变换时，若一个基组有赝势、另一个基组没有赝势（比如全电子基组），此时guess=read效果很差（因为赝势的原因，导致少了一些轨道供读取），常伴随SCF剧烈震荡或能量变化巨大，此时前一个基组下的波函数稳定性 不能类推 后一个基组的波函数稳定性。在后一个基组下，也应当做波函数稳定性分析。

顺带一提，用B3LYP泛函算双（多）核过渡金属，属于挖坑等着审稿人把投稿人扔进去的行为
《简谈量子化学计算中DFT泛函的选择》http://sobereva.com/272

LYFE

zjxitcc 发表于 2023-5-18 13:42
cc-pVTZ对Cu元素有定义，且是全电子基组，所以不应该写“ub3lyp/cc-pvtz pseudo=read”，坐标底下也不应 ...

1. #p ub3lyp/cc-pvtz nosymm scf=xqc int=ultrafine IOp(3/76=1000001500) IOp(3/77=0720008500) IOp(3/78=0810010000)

复制代码

修改了关键词，且坐标底下无Cu的赝式信息，没有再使用guess=read。继续用B3LYP是因为之前都是用此泛函为了对比，接着使用了B3LYP。但是结果显示还是没有自旋分布。

zjxitcc

LYFE 发表于 2023-6-5 18:40
修改了关键词，且坐标底下无Cu的赝式信息，没有再使用guess=read。继续用B3LYP是因为之前都是用此泛函为 ...

体系自旋多重度是啥；输出文件里最后一次出现的SCF Done下方1~2行<S**2>值是多少？

LYFE

zjxitcc 发表于 2023-6-5 19:17
体系自旋多重度是啥；输出文件里最后一次出现的SCF Done下方1~2行值是多少？

1. Mulliken charges and spin densities:
   
2.                1          2
   
3.      1  Cu   0.841725   0.000000
   
4.      2  Cu   0.841727  -0.000000
   
5.      3  O   -0.480897   0.000000
   
6.      4  O   -0.588967  -0.000000
   
7.      5  C    0.145468  -0.000000
   
8.      6  C   -0.087350   0.000000
   
9.      7  C   -0.064267  -0.000000
   
10.      8  C   -0.055238  -0.000000
    
11.      9  C   -0.012288   0.000000
    
12.     10  N   -0.213742  -0.000000
    
13.     11  H    0.163517  -0.000000
    
14.     12  H    0.153673  -0.000000
    
15.     13  H    0.161739  -0.000000
    
16.     14  H    0.155878  -0.000000
    
17.     15  C   -0.125479   0.000000
    
18.     16  N    0.108237   0.000000
    
19.     17  C   -0.129761   0.000000
    
20.     18  C    0.130359  -0.000000
    
21.     19  H    0.168635   0.000000
    
22.     20  N   -0.159641  -0.000000
    

复制代码

第二列为自旋多重度。
以下是最后一次出现SCF Done

1. SCF Done:  E(UB3LYP) =  -4802.18171400     A.U. after   21 cycles
   
2.             NFock= 21  Conv=0.78D-08     -V/T= 2.0030
   
3. <Sx>= 0.0000 <Sy>= 0.0000 <Sz>= 0.0000 <S**2>=-0.0000 S=-0.0000
   
4. <L.S>=  0.00000000000   
   
5. KE= 4.787737432010D+03 PE=-2.192152860455D+04 EE= 7.032591194270D+03
   
6. Annihilation of the first spin contaminant:
   
7. S**2 before annihilation    -0.0000,   after    -0.0000

复制代码

wzkchem5

LYFE 发表于 2023-6-7 11:52
第二列为自旋多重度。
以下是最后一次出现SCF Done

人家问的是自旋多重度，不是自旋布居数

LYFE

zjxitcc 发表于 2023-6-5 19:17
体系自旋多重度是啥；输出文件里最后一次出现的SCF Done下方1~2行值是多少？

自旋多重度设置的为开壳层单重态

zjxitcc

LYFE 发表于 2023-6-7 18:52
第二列为自旋多重度。
以下是最后一次出现SCF Done

是否检验过波函数稳定性？
PS：如果你觉得cc-pVTZ基组太大，计算过于耗时，可以改成def2TZVP。