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

24小时热门版块排行榜    

小木虫论坛-学术科研互动平台 » 计算模拟区 » 量子化学 » Gaussian » 金属配合物DFT计算，老是出现自旋污染，怎么解决啊
91/1返回列表
查看: 1786  |  回复: 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英文润色
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[基金申请] 年年想打听，年年打听不到。。 +9 hdzw9071 2024-06-03 10/500 2024-06-03 13:35 by bnullh
[基金申请] 前些天开会有个人见到人就搞关系，一查此人全是MDPI/Hindawi论文，鄙视! +17 zju2000 2024-06-02 17/850 2024-06-03 13:22 by thatday77
[基金申请] 为什么美国要封锁华为而不封锁丘C同呢？ +13 zju2000 2024-06-02 16/800 2024-06-03 11:54 by chemhua
[基金申请] B口人才项目 +11 WOWO159357 2024-05-29 23/1150 2024-06-03 11:18 by XUAAAAA
[硕博家园] 实验室太吵闹，无法安静学习，怎么办？ +11 utahh 2024-05-31 22/1100 2024-06-03 11:18 by utahh
[基金申请] 化学B02口青基 代表作都是什么水平的？向大佬求助 +14 arthas_007 2024-06-01 19/950 2024-06-03 11:14 by prisonrabbit
[基金申请] 离职后国自然项目不能变更单位，在新单位还能申请新的国自然项目吗 +3 宜兰 2024-05-31 11/550 2024-06-03 10:26 by lilyxjw
[硕博家园] 材料博士毕业 +3 112233ssg 2024-06-02 3/150 2024-06-03 09:59 by Lwyy1106
[基金申请] 化学口B0109(高分子合成)，拿青年基金一般需要怎样的文章水平？ +14 salmon95 2024-05-30 24/1200 2024-06-03 08:31 by LNP@mRNA
[硕博家园] 每到中夜，情难自抑 +36 sioc-sunj 2024-05-28 60/3000 2024-06-03 02:18 by 梨花儿
[考博] 求生物学的博导收留 +6 申博成功功成博� 2024-06-01 6/300 2024-06-03 01:23 by ww88
[考博] 导师不让硕转博，让我去国外读博，能理解吗？ +12 萧山幽谷 2024-05-29 20/1000 2024-06-02 12:01 by yuan0806
[教师之家] 在大地上我们只过一生---看完我的阿勒泰上头了好几天，完结那天晚上几乎失眠 +11 瞬息宇宙 2024-05-27 13/650 2024-06-01 22:15 by otani
[论文投稿] 求Sci期刊推荐 10+4 甄小鱼 2024-05-30 7/350 2024-06-01 10:41 by bobvan
[基金申请] 博后特别资助状态变化 +24 随梦而飞2017 2024-05-30 35/1750 2024-06-01 10:10 by 青岛阳仔
[考博] 申博求助?本硕双非一篇三区和四区去985工科非天坑专业有没有希望？ +4 LYS1200 2024-05-29 6/300 2024-06-01 08:50 by 学术渣渣神
[高分子] MMA预聚体光固化发雾问题求助 +3 惠亚金总 2024-05-29 10/500 2024-05-31 14:59 by 惠亚金总
[论文投稿] 高手朋友推荐比较容易投稿和录用的SCI期刊，不投稿收费SC，对分区没有要求 5+3 xintangren 2024-05-28 4/200 2024-05-29 10:46 by xintangren
[论文投稿] 核心初审被拒，理由是“选题的意义不明确，文章写得不像是科技论文”，怎么改 5+3 工藤雷花樱 2024-05-27 8/400 2024-05-29 10:09 by topedit
[基金申请] 如果您是国自然评审专家 +3 丁香园账户 2024-05-28 3/150 2024-05-29 06:44 by gaohui8888
信息提示
关闭
请填处理意见
关闭