| 查看: 388 | 回复: 2 | |||
| 当前主题已经存档。 | |||
[交流]
自旋污染求助
|
|||
|
最近做的一个体系自旋污染大的惊人，请教高手 输入：# 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. 而且我也作了四重态，发现二重态比四重态的能量低，现在不知道该怎么办了，还请高手指点， 谢谢！ |
回复此楼» 猜你喜欢
夏威夷大学生物胶体与界面实验室诚招全奖博士生和博士后
已经有1人回复
-
求助立方晶体In(OH)3的CIF卡片
已经有2人回复
-
物理化学论文润色/翻译怎么收费?
已经有132人回复
-
大湾区大学刘天辉课题组招收2026级博士生（与中山大学物理学院联合培养）
已经有22人回复
-
PbS量子点紫外吸收
已经有0人回复
-
比利时鲁汶大学与国家留学基金委共同资助博士研究生CSC-KU Leuven PhD Scholarship
已经有0人回复
-
七嗪类物质合成求助
已经有0人回复
-
[电子教材]天津大学物理化学（第七版）上、下册
已经有1人回复
-
阴离子交换膜电解二氧化碳还原
已经有0人回复
-
伦敦大学学院Benjamin Moss博士招收人工智能、原位光谱、催化相关26年博士生
已经有11人回复
beefly
专家顾问 (职业作家)
地沟油冶炼专家
-
专家经验: +458 - QC强帖: 12
- 应助: 88 (初中生)
- 金币: 12841.4
- 散金: 27307
- 红花: 300
- 帖子: 4667
- 在线: 925.7小时
- 虫号: 408372
- 注册: 2007-06-21
- 性别: GG
- 专业: 理论和计算化学
- 管辖: 计算模拟
2楼2008-01-13 01:03:03
3楼2008-01-13 20:31:57