24小时热门版块排行榜    

Znn3bq.jpeg
查看: 418  |  回复: 2
当前主题已经存档。

jghe

金虫 (著名写手)

[交流] 自旋污染求助

最近做的一个体系自旋污染大的惊人,请教高手
输入:# b3lyp/6-311+g(2d,p) pop=NBO scf=(maxcycle=500, tight)
Charge =  0 Multiplicity = 2

输出:
SCF Done:  E(UB+HF-LYP) =  -5204.25052453     A.U. after   32 cycles
             Convg  =    0.6735D-08             -V/T =  2.0016
             S**2   =   1.9323
Annihilation of the first spin contaminant:
S**2 before annihilation     1.9323,   after     2.5168

高斯公司回信:I agree that this is a large value for S**2 but this really is a
statement about system rather than something you can control.

   1) have you check the stability of the wavefunction?  You have a number
      of negative eigenvalues, run STABLE=Opt Guess=Read

   2) are you sure that a doublet is the lowest energy state, a large
      value for S**2 says the alpha and beta orbitals are different,
      this can occur if you have the highest beta occupied different
      from the alpha orbitals where it may be unoccupied.  Again stability
      may switch this or higher spin states may be more stable.

Try both of these.

可我做了STABLE=Opt 后发现波函数是稳定的
Eigenvectors of the stability matrix:

Eigenvector   1:   ?Spin  -A'   Eigenvalue= 0.0295043
    111A ->112A        0.93202
    110B ->111B       -0.18620
    110B ->113B       -0.10726
The wavefunction is stable under the perturbations considered.
The wavefunction is already stable.
The wavefunction is already stable.

而且我也作了四重态,发现二重态比四重态的能量低,现在不知道该怎么办了,还请高手指点, 谢谢!
回复此楼
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

beefly

专家顾问 (职业作家)

地沟油冶炼专家

★ ★
zzgyb(金币+2,VIP+0):谢谢您的参与,欢迎您常来计算模拟版。
用RO(DFT),比如ROB3LYP
beefly《西太平洋大学现代英汉词典》[bi:fli]牛肉一般地
2楼2008-01-13 01:03:03
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

jghe

金虫 (著名写手)

RoB3LYP 优化总是有问题,一圈都跑不下去, L502断掉。 有其它更好的方法没?
3楼2008-01-13 20:31:57
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 jghe 的主题更新
普通表情 高级回复 (可上传附件)
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 293调剂 +25 yj1221 2026-04-08 26/1300 2026-04-10 15:02 by 柴小白
[考研] 材料工程085601,270求调剂 +27 @ASDF1234 2026-04-08 29/1450 2026-04-10 10:48 by 紫翼精灵
[考研] 材料与化工371求调剂 +17 陪琳看海 2026-04-04 18/900 2026-04-10 07:51 by 314126402
[考研] 289求调剂 +4 L1ttleTiger 2026-04-04 4/200 2026-04-09 19:53 by xiayizhi
[考研] 一志愿中科院105500专业总分315求调剂 +6 lallalh 2026-04-09 7/350 2026-04-09 17:51 by lallalh
[考研] 296求调剂 +5 汪!?! 2026-04-09 5/250 2026-04-09 17:47 by 柠檬不酸zy
[考研] 1U盾记得记得就 +9 sanjin020722 2026-04-08 10/500 2026-04-09 14:11 by 诗与自由
[考研] 求调剂 +8 吃口冰激凌 2026-04-07 8/400 2026-04-09 08:03 by 5268321
[考研] 化学0703-一志愿211-338分求调剂 +10 vants 2026-04-05 11/550 2026-04-08 16:02 by screening
[考研] 求调剂,现在还能填的 +3 上岸小莹加油 2026-04-08 3/150 2026-04-08 14:30 by zhq0425
[考研] 计算机11408 287 求调剂 +3 LiLe5 2026-04-07 3/150 2026-04-07 23:15 by shanqishi
[考研] 0854电子信息319求调剂(接受跨专业调剂) +5 星星不眨眼喽 2026-04-05 6/300 2026-04-07 22:16 by hemengdong
[考研] 085602调剂 初试总分335 +10 19123253302 2026-04-05 10/500 2026-04-07 15:23 by 小乔同学ya
[考研] 312求调剂 +4 LR6 2026-04-06 4/200 2026-04-07 08:42 by jp9609
[考研] 331求调剂 +8 于征yz 2026-04-05 8/400 2026-04-06 00:54 by fmesaito
[考研] 工科求调剂 +15 11ggg 2026-04-03 15/750 2026-04-05 16:24 by zzx2138
[考研] 081200-11408-276学硕求调剂 +4 崔wj 2026-04-04 5/250 2026-04-05 14:06 by imissbao
[考研] 调剂 +8 熊二想上岸 2026-04-04 8/400 2026-04-05 05:27 by houyaoxu
[考研] 323求调剂 +8 李佳乐1 2026-04-04 8/400 2026-04-04 22:26 by hemengdong
[考研] 22408,264求调剂 +3 ywh729 2026-04-03 4/200 2026-04-04 11:04 by ywh729
信息提示
请填处理意见