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

24小时热门版块排行榜    

小木虫论坛-学术科研互动平台 » 计算模拟区 » 量子化学 » Gaussian » 金属配合物DFT计算，老是出现自旋污染，怎么解决啊
91/1返回列表
查看: 1750  |  回复: 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的回帖

竖走的螃蟹

新虫 (小有名气)
引用回帖:
2楼: Originally posted by paramecium86 at 2023-06-24 09:36:37
这并非什么自旋污染 而是因为过渡金属体系SCF经常不容易收敛。加一些关键词就可以解决。可以看我 http://muchong.com/t-14242735-1 这个帖子 最新教程1.6版的2.2节
大体内容如下

***************************** ...

太感谢您了，我这就去试试

发自小木虫Android客户端
赞一下 回复此楼
3楼2023-06-26 15:07:41
 已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

竖走的螃蟹

新虫 (小有名气)
引用回帖:
2楼: Originally posted by paramecium86 at 2023-06-24 09:36:37
这并非什么自旋污染 而是因为过渡金属体系SCF经常不容易收敛。加一些关键词就可以解决。可以看我 http://muchong.com/t-14242735-1 这个帖子 最新教程1.6版的2.2节
大体内容如下

***************************** ...

在优化后的坐标再优化一次不知道能不能降低这种情况出现
赞一下 回复此楼
4楼2023-06-26 15:19:34
 已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

竖走的螃蟹

新虫 (小有名气)
引用回帖:
5楼: Originally posted by paramecium86 at 2023-06-26 16:15:25
几何优化和SCF是独立过程。前者已经优化好的结构并不一定可以保证让SCF一定可以收敛。...

老师再请教你个问题，我是三价Ru离子，拿羟基配位，钌本身在强配体场下多重度会变成2，而我去配位得到的是电中性的东西，不知道这个电荷数和多重度该是多少，这里我发现比较难明白，因为中性分子会理论上会占用钌离子价态，不知道这个电荷数会不会改变

发自小木虫Android客户端
赞一下 回复此楼
6楼2023-06-27 08:50:08
 已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

竖走的螃蟹

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

理论上每个羟基都会占用一个键，形成-O-Ru键

发自小木虫Android客户端
赞一下 回复此楼
7楼2023-06-27 09:05:17
 已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

竖走的螃蟹

新虫 (小有名气)
引用回帖:
8楼: Originally posted by paramecium86 at 2023-06-27 17:16:11
电中性的话净电荷就是0 ，这个和中心原子价态没啥直接关系。 而自旋多重度 最简单有效的办法首先是要数一下整个体系的电子数目。可以看我2楼提到的教程的第一节内容 就是判断自旋多重度的方法。...

老师，电中性的我知道净电荷数是0，多重度是2，但是假如是离子型，外层带一个负电荷，它的多重度没有变化也是2，这里带金属那块的净电荷是理论上是1，我设置1 2，后又出现不收敛问题了。之前中性不收敛我优化好结构后就没有出现了，但是离子型的配合物又出现问题了
赞一下 回复此楼
9楼2023-06-30 09:47:43
 已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

竖走的螃蟹

新虫 (小有名气)
送红花一朵
引用回帖:
10楼: Originally posted by paramecium86 at 2023-06-30 12:21:20
不收敛就还是按照SCF不收敛的办法去应付。解决办法都是通用的。...

老师，我发现我按照B3LYP泛函 以6-31G计算轻原子，以LANL2DZ方法计算重原子Ru计算出来的振子强度和文献的差别有十倍之大，这有什么解决方法吗
赞一下 回复此楼
11楼2023-07-06 15:46:46
 已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

竖走的螃蟹

新虫 (小有名气)
引用回帖:
10楼: Originally posted by paramecium86 at 2023-06-30 12:21:20
不收敛就还是按照SCF不收敛的办法去应付。解决办法都是通用的。...

算TD振子强度每次都特别低是怎么回事，比实际低数百倍。
%chk=Ru-ql-TD.chk
%mem=5GB
%NProcShared=8#p td=(nstates=20) ub3lyp/genecp scrf=(solvent=dichloromethane)这是我的关键词，不知道是哪里出了问题

发自小木虫Android客户端
赞一下 回复此楼
13楼2023-11-08 15:45:30
 已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 竖走的螃蟹 的主题更新
91/1返回列表
不应助 确定回帖应助 (注意：应助才可能被奖励，但不允许灌水，必须填写15个字符以上)
普通表情 龙 兔 虎 猫 论文翻译-EnPapers美助理教授团队-SCI英文润色
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[基金申请] 国社科申请书上传有误，学校已提交到省里，省里还未审核，还能退回修改嘛？ 100+3 远山晴岚 2024-05-19 6/300 2024-05-21 14:14 by 漠上藜梭
[考博] 换导师 +14 是柠檬呀！ 2024-05-18 27/1350 2024-05-21 13:27 by 是柠檬呀！
[考博] 本科出身不好是不是会被直接刷呀 +4 未来富婆蛙 2024-05-20 4/200 2024-05-21 13:20 by 慎独的小花卷
[基金申请] 太诡异了，五月底还有没有送审的。。 +12 hdzw9071 2024-05-21 12/600 2024-05-21 12:43 by dxcharlary
[基金申请] 听说面青地E09已经送了么？ +5 叉烧吃叉烧 2024-05-21 7/350 2024-05-21 11:35 by lxq99005
[找工作] 浙江的高校现在门槛都这么高吗 +11 dadqweq_qw 2024-05-16 11/550 2024-05-21 10:30 by jurkat.1640
[基金申请] 这个模块怎么成了烧香拜佛的地方了 +7 shrz98 2024-05-18 7/350 2024-05-21 10:26 by lancet0903
[硕博家园] 又见奇葩学生，可咋办 +25 appleapple2 2024-05-17 30/1500 2024-05-21 09:13 by qyws
[考博] 双非博士还是985科研助理，然后再读博 +6 lxdatj123 2024-05-18 13/650 2024-05-21 08:19 by lxdatj123
[教师之家] 青椒的未来在哪里？ +17 sxhuabin 2024-05-15 18/900 2024-05-20 22:55 by 瞬息宇宙
[基金申请] 青年拔尖出炉 +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
[硕博家园] 导师这样说到底有没有机会啊……… +9 理科生yes 2024-05-14 16/800 2024-05-20 15:28 by holypower
[基金申请] 基金评审 +4 阿呆不呆 2024-05-20 4/200 2024-05-20 15:10 by 一路向东
[论文投稿] Nano-Micro Letters投稿可以用PDF上传吗 5+4 962511369 2024-05-18 4/200 2024-05-19 15:47 by ca0yan9
[无机/物化] 请问什么溶剂能溶解二氧化锰 (金币+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
[留学DIY] 求助：医生想润 +7 safrina 2024-05-15 7/350 2024-05-17 18:50 by 卡皮吧啦123
[基金申请] 粗大，科研之友提示有一篇10多年前的文章被阅读 +9 babu2015 2024-05-16 11/550 2024-05-16 14:53 by TsolmonZ
[教师之家] 问题已解觉，谢谢大家关注！ +7 lzgj258 2024-05-15 11/550 2024-05-15 19:15 by 环境检测2024
信息提示
关闭
请填处理意见
关闭