计算化学公社

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

作者
Author: sobereva    时间: 2023-5-15 19:55
使用赝势基组应当写genecp而非gen

倘若你确保关键词没写错，大基组下做单点时带着stable关键词做波函数稳定性测试
作者
Author: zjxitcc    时间: 2023-5-16 10:10
很大概率是输入文件写错了、关键词不合理。建议晒晒输入文件（坐标可以不贴，其他内容原样复制，且用论坛上发帖功能中的代码块显示）
作者
Author: LYFE    时间: 2023-5-16 11:08
  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)

  6. Title Card Required
复制代码
作者
Author: LYFE    时间: 2023-5-18 09:51
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)

  6. Title Card Required
复制代码
作者
Author: zjxitcc    时间: 2023-5-18 13:42
本帖最后由 zjxitcc 于 2023-5-18 13:46 编辑
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

作者
Author: LYFE    时间: 2023-6-5 18:40
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。但是结果显示还是没有自旋分布。
作者
Author: zjxitcc    时间: 2023-6-5 19:17
LYFE 发表于 2023-6-5 18:40
修改了关键词，且坐标底下无Cu的赝式信息，没有再使用guess=read。继续用B3LYP是因为之前都是用此泛函为 ...

体系自旋多重度是啥；输出文件里最后一次出现的SCF Done下方1~2行<S**2>值是多少？
作者
Author: LYFE    时间: 2023-6-7 18:52
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
复制代码
作者
Author: wzkchem5    时间: 2023-6-7 21:09
LYFE 发表于 2023-6-7 11:52
第二列为自旋多重度。
以下是最后一次出现SCF Done

人家问的是自旋多重度，不是自旋布居数
作者
Author: LYFE    时间: 2023-6-7 23:12
zjxitcc 发表于 2023-6-5 19:17
体系自旋多重度是啥；输出文件里最后一次出现的SCF Done下方1~2行值是多少？

自旋多重度设置的为开壳层单重态
作者
Author: zjxitcc    时间: 2023-6-8 12:04
LYFE 发表于 2023-6-7 18:52
第二列为自旋多重度。
以下是最后一次出现SCF Done


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





欢迎光临 计算化学公社 (http://bbs.keinsci.com/) Powered by Discuz! X3.3