应《网络安全法》要求，自2017年10月1日起，未进行实名认证将不得使用互联网跟帖服务。为保障您的帐号能够正常使用，请尽快对帐号进行手机号验证，感谢您的理解与支持！

24小时热门版块排行榜    

小木虫论坛-学术科研互动平台 » 计算模拟区 » 量子化学 » Gaussian » 金属配合物DFT计算，老是出现自旋污染，怎么解决啊
71/1返回列表
查看: 1749  |  回复: 13
查看全部回帖@他人 存档 新回复提醒 (忽略) 收藏    在APP中查看
【悬赏金币】回答本帖问题，作者竖走的螃蟹将赠送您 5 个金币
当前只显示满足指定条件的回帖，点击这里查看本话题的所有回帖

竖走的螃蟹

新虫 (小有名气)

[求助] 金属配合物DFT计算，老是出现自旋污染，怎么解决啊已有1人参与

关键词
%chk=Ru-CF3-TD-1.chk
%mem=5GB
%NProcShared=8
# opt td=(nstates=6) ub3lyp/genecp geom=connectivity

Charge =  0 Multiplicity = 2
C                    -3.89701  -1.49571  -0.09117
C                    -3.30042  -2.7503   -0.23496
C                    -1.93259  -2.45736  -0.19552
N                    -1.79447  -1.10611  -0.03461
H                    -3.78839  -3.70652  -0.35021
C                    -0.71996  -3.23521  -0.29002
C                    -0.6809   -4.63041  -0.4446
C                     0.54402  -5.27871  -0.5095
H                    -1.61265  -5.18196  -0.50841
C                     1.62081  -3.15155  -0.26627
C                     1.72148  -4.52898  -0.41479
H                     0.58502  -6.35762  -0.6288
H                     2.5008   -2.5251   -0.17354
H                     2.69907  -4.99678  -0.45363
N                    -2.97823  -0.50898   0.02766
C                    -5.3539   -1.1785   -0.05414
F                    -5.73139  -0.61996   1.11669
F                    -5.71962  -0.31979  -1.0311
F                    -6.08587  -2.30841  -0.21551
N                     0.4365   -2.51353  -0.21047
Ru                    0.16101  -0.33832   0.01166
C                    -0.56      1.88335   1.69652
C                    -1.00284   3.56738   0.36392
H                    -1.2167    3.82463   2.57821
N                    -0.37124   1.53614   0.37565
C                    -0.97105   3.20599   1.72837
N                    -0.64059   2.5587   -0.43918




Defaulting to unpruned grid for atomic number  44.
Defaulting to unpruned grid for atomic number  44.
Defaulting to unpruned grid for atomic number  44.
Defaulting to unpruned grid for atomic number  44.
Defaulting to unpruned grid for atomic number  44.
Defaulting to unpruned grid for atomic number  44.
Defaulting to unpruned grid for atomic number  44.
Defaulting to unpruned grid for atomic number  44.
FoFCou: FMM=F IPFlag=           0 FMFlag=      100000 FMFlg1=        2001
         NFxFlg=           0 DoJE=T BraDBF=F KetDBF=T FulRan=T
         wScrn=  0.000000 ICntrl=     500 IOpCl=  0 I1Cent=   200000004 NGrid=           0
         NMat0=    1 NMatS0=      1 NMatT0=    0 NMatD0=    1 NMtDS0=    0 NMtDT0=    0
Petite list used in FoFCou.
Requested convergence on RMS density matrix=1.00D-08 within 128 cycles.
Requested convergence on MAX density matrix=1.00D-06.
Requested convergence on             energy=1.00D-06.
No special actions if energy rises.
Defaulting to unpruned grid for atomic number  44.
Defaulting to unpruned grid for atomic number  44.
Defaulting to unpruned grid for atomic number  44.
Defaulting to unpruned grid for atomic number  44.
Defaulting to unpruned grid for atomic number  44.
Defaulting to unpruned grid for atomic number  44.
Defaulting to unpruned grid for atomic number  44.
Defaulting to unpruned grid for atomic number  44.
Restarting incremental Fock formation.
Restarting incremental Fock formation.
Restarting incremental Fock formation.
Restarting incremental Fock formation.
Restarting incremental Fock formation.
Restarting incremental Fock formation.
>>>>>>>>>> Convergence criterion not met.
SCF Done:  E(UB3LYP) =  -2523.16279936     A.U. after  129 cycles
            NFock=128  Conv=0.59D-06     -V/T= 2.0311
<Sx>= 0.0000 <Sy>= 0.0000 <Sz>= 0.5000 <S**2>= 0.7548 S= 0.5024
<L.S>= 0.000000000000E+00
Annihilation of the first spin contaminant:
S**2 before annihilation     0.7548,   after     0.7500
Convergence failure -- run terminated.
Error termination via Lnk1e in /share/apps/g09.D01/g09/l502.exe at Tue Jun 20 18:23:52 2023.
Job cpu time:       5 days  5 hours  3 minutes 26.5 seconds.
File lengths (MBytes):  RWF=   1697 Int=      0 D2E=      0 Chk=    175 Scr=      1

有哪位虫友知道这个怎么解决吗？跪谢
回复此楼

» 猜你喜欢
1楼2023-06-23 13:14:43
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

paramecium86

版主 (著名写手)

【答案】应助回帖

感谢参与，应助指数 +1
这并非什么自旋污染 而是因为过渡金属体系SCF经常不容易收敛。加一些关键词就可以解决。可以看我 http://muchong.com/t-14242735-1 这个帖子 最新教程1.6版的2.2节
大体内容如下

********************************************************************************
“”2.2 SCF不收敛 （l502.exe error）
报错信息一般如下
>>>>>>>>>> Convergence criterion not met.
SCF Done:  E(RB3LYP) =  -1111.12727817     A.U. after  129 cycles
            NFock=128  Conv=0.76D-04     -V/T= 2.7410
KE= 6.382296294020D+02 PE=-5.731294386816D+03 EE= 2.414551554017D+03
Convergence failure -- run terminated.
Error termination via Lnk1e in F:\Gaussian\G09W\l502.exe


A. 首先检查输入结构是不是合理。绝大部分情况下都是因为结构跟 “真实” 结构差的很远，导致的SCF难以收敛。尽量给一个能给出的最好结构 （比如先用UFF粗略优化一下，再用PM6 在优化一下。然后再用DFT计算之类的）
B. 检查体系的自旋多重度是否设置的合理。尤其是过渡金属体系。如果多重度设置的偏离 “真实” 情况很远，则非常难以收敛
C. 检查基组是否合理。尤其是包含重原子，是否合理的设置了混合基组之类的都会影响。具体参考本文1.1节设置混合基组和赝势的方法。
D. 检查当前给的初始结构是否真实存在于 “气相” 中。因为有时候气相和溶液相结构很不一样。这就会导致初始结构以为很准其实是很差的情况。如果是溶液中的结构则需要加上溶剂模型SCRF，见本文1.2 .有时候加了溶剂模型就可以帮助收敛。
E. 以上确定是没有问题的，这时候就考虑加入一些关键词帮助收敛
SCF=(XQC,MAXCONVENTIONALCYCLES=64)
这个关键词可以应用在绝大部分情形下。但是绝不要一上来就用，一定要检查之前的步骤一都没问题了再去使用。G09的话就用SCF=XQC就好。 G16可以写这一串。
SCF=NoVarAcc
高斯在计算SCF过程中，最开始的几步为了加速会降低积分精度。这样有时候就会让波函数不小心进入了一个错误的区间。很难再优化下去。所以可以用这个关键词，取消这种加速降低积分精度的方法。
SCF(VSHIFT=500) 这个一般可以用在过渡金属体系。这种体系，可能HOMO-LUMO间距过小导致波函数不容易收敛。可以计算中人为太高这个值。这个不会影响结果只会影响收敛过程。一般设定在500上下比较合理。

INT(Acc2E=12, UltraFine)
如果用G09进行计算，也可以考虑这个关键词提高了DFT的积分精度和格点精度。这个仅仅适用在某些DFT计算中。比如M062X 这类对于格点要求比较高的泛函中。G16中这个已经是默认值不需要写。

实际计算可以用以上提到的一个或者几个关键词一起帮助收敛。
如果体系实在是太过巨大，这些都做了还是很难收敛。可以使用

SCF(CONVER=6)
这个关键词降低SCF收敛限到10-6， 而默认是10-8。但是慎用。这样的结构可能得到的根本不准。最好是优化完成后，读取波函数和结构恢复到默认的收敛限再进行一次优化。

*****************************************************************************************************
赞一下 回复此楼
2楼2023-06-24 09:36:37
 已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

paramecium86

版主 (著名写手)
引用回帖:
4楼: Originally posted by 竖走的螃蟹 at 2023-06-26 02:19:34
在优化后的坐标再优化一次不知道能不能降低这种情况出现...

几何优化和SCF是独立过程。前者已经优化好的结构并不一定可以保证让SCF一定可以收敛。
赞一下 回复此楼
5楼2023-06-26 16:15:25
 已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

paramecium86

版主 (著名写手)
引用回帖:
6楼: Originally posted by 竖走的螃蟹 at 2023-06-26 19:50:08
老师再请教你个问题，我是三价Ru离子，拿羟基配位，钌本身在强配体场下多重度会变成2，而我去配位得到的是电中性的东西，不知道这个电荷数和多重度该是多少，这里我发现比较难明白，因为中性分子会理论上会占用钌离 ...

电中性的话净电荷就是0 ，这个和中心原子价态没啥直接关系。 而自旋多重度 最简单有效的办法首先是要数一下整个体系的电子数目。可以看我2楼提到的教程的第一节内容 就是判断自旋多重度的方法。
赞一下(1人) 回复此楼
8楼2023-06-27 17:16:11
 已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

paramecium86

版主 (著名写手)
引用回帖:
9楼: Originally posted by 竖走的螃蟹 at 2023-06-29 20:47:43
老师，电中性的我知道净电荷数是0，多重度是2，但是假如是离子型，外层带一个负电荷，它的多重度没有变化也是2，这里带金属那块的净电荷是理论上是1，我设置1 2，后又出现不收敛问题了。之前中性不收敛我优化好结构 ...

不收敛就还是按照SCF不收敛的办法去应付。解决办法都是通用的。
赞一下 回复此楼

» 本帖已获得的红花（最新10朵）

竖走的螃蟹
10楼2023-06-30 12:21:20
 已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

paramecium86

版主 (著名写手)
引用回帖:
11楼: Originally posted by 竖走的螃蟹 at 2023-07-06 02:46:46
老师，我发现我按照B3LYP泛函 以6-31G计算轻原子，以LANL2DZ方法计算重原子Ru计算出来的振子强度和文献的差别有十倍之大，这有什么解决方法吗...

IR谱一般看计算出来的频率也就是峰的位置准不准。而强度的绝对值也没啥意义 还是看相对高度就可以。而且强度也没有那么准。 其次基组用的太差了。这个年代也没办法用6-31G和LanL2DZ了 太小了 至少得  6-31G(D) 和 SDD
赞一下 回复此楼
12楼2023-07-06 18:19:43
 已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

paramecium86

版主 (著名写手)
引用回帖:
13楼: Originally posted by 竖走的螃蟹 at 2023-11-08 02:45:30
算TD振子强度每次都特别低是怎么回事，比实际低数百倍。
%chk=Ru-ql-TD.chk
%mem=5GB
%NProcShared=8#p td=(nstates=20) ub3lyp/genecp scrf=(solvent=dichloromethane)这是我的关键词，不知道是哪里出了问题
...

振子强度低数百倍是啥意思。只看峰相对的峰高即可。绝对值没有什么实际的意义。还是说你要用振子强度去算寿命什么的
赞一下 回复此楼
14楼2023-11-08 17:19:38
 已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 竖走的螃蟹 的主题更新
71/1返回列表
不应助 确定回帖应助 (注意：应助才可能被奖励，但不允许灌水，必须填写15个字符以上)
普通表情 龙 兔 虎 猫 论文翻译-EnPapers美助理教授团队-SCI英文润色
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考博] 急求博导收留 +15 8139173 2024-05-15 17/850 2024-05-21 11:43 by 82年拉菲
[基金申请] 听说面青地E09已经送了么？ +5 叉烧吃叉烧 2024-05-21 7/350 2024-05-21 11:35 by lxq99005
[论文投稿] Angew拒稿重投经历大修返回又送审了 48+4 YOONAZHA 2024-05-18 6/300 2024-05-21 10:35 by jurkat.1640
[硕博家园] 蹲一个男朋友 +23 伊伊莎贝拉 2024-05-17 38/1900 2024-05-21 10:32 by 未来可期简直不�
[硕博家园] 民办高校入职后稳定吗？ +15 905452934 2024-05-16 23/1150 2024-05-21 10:20 by balabiubiu
[基金申请] 审不上青基又非升即走的青椒 和 牢里踩缝纫机的犯人哪个活的更舒服一点？ +10 非非飞远了 2024-05-20 11/550 2024-05-21 09:47 by 芙兰朵露
[教师之家] 另一半是老师会加分哦 +14 acogw9004 2024-05-14 14/700 2024-05-21 08:54 by bxwisdom_886
[考博] 双非博士还是985科研助理，然后再读博 +6 lxdatj123 2024-05-18 13/650 2024-05-21 08:19 by lxdatj123
[基金申请] 申请基金代表性成果 +14 lancet0903 2024-05-17 20/1000 2024-05-21 00:23 by dxcharlary
[硕博家园] 被导师抓着一个褥羊毛该咋办 +13 未来可期简直不� 2024-05-15 25/1250 2024-05-20 23:03 by ca0yan9
[基金申请] 青年拔尖出炉 +6 babu2015 2024-05-14 7/350 2024-05-20 20:20 by hema6g13
[考博] 海南大学国家级人才田新龙教授团队氢燃料电池与器件方向招收科研助理 (重点支持方向) +6 Dreammzp 2024-05-17 9/450 2024-05-20 18:40 by Dreammzp
[有机交流] 除DMSO +5 Spiralup 2024-05-15 6/300 2024-05-20 17:49 by 刘洪振
[教师之家] 普通高校有编制讲师与211高校无编制专职科研岗，应该怎么选 +19 1137025804 2024-05-15 25/1250 2024-05-20 13:32 by Zhangsanxi
[考博] 【2025 申博】材料或者冶金工程 +4 枫落孤城 2024-05-19 5/250 2024-05-20 10:52 by 枫落孤城
[基金申请] 有知道工材e01情况的吗？ +13 xiaopang8958 2024-05-17 19/950 2024-05-20 09:52 by 蜕变123~
[无机/物化] 请问什么溶剂能溶解二氧化锰 (金币+5) +4 这是春天 2024-05-15 4/200 2024-05-19 11:24 by wanghuawei
[考博] 双非院士组还是211普通老师 +9 timergoi 2024-05-14 10/500 2024-05-18 08:31 by 3001160025
[基金申请] 粗大，科研之友提示有一篇10多年前的文章被阅读 +9 babu2015 2024-05-16 11/550 2024-05-16 14:53 by TsolmonZ
[基金申请] 有人看了我的科研之友 +9 jurkat.1640 2024-05-14 12/600 2024-05-15 16:04 by arthas_007
信息提示
关闭
请填处理意见
关闭