24小时热门版块排行榜    

查看: 923  |  回复: 8

t13340033021

金虫 (正式写手)

[求助] Stable 关键字计算

金属体系。stable 测试波函数稳定性。


***********************************************************************
Stability analysis using singles matrix:
***********************************************************************

Eigenvectors of the stability matrix:

Eigenvector   1:  1.000-AG   Eigenvalue=-0.0122008  =0.000
     91A -> 92A        0.69385
     91B -> 92B       -0.69385

Eigenvector   2:  1.000-AG   Eigenvalue=-0.0122008  =0.000
     91A -> 93A        0.69385
     91B -> 93B       -0.69385

Eigenvector   3:  1.000-AU   Eigenvalue=-0.0115697  =0.000
     89A -> 92A       -0.42362
     89A -> 93A       -0.20992
     90A -> 92A       -0.20993
     90A -> 93A        0.42362
     89B -> 92B        0.42364
     89B -> 93B        0.20989
     90B -> 92B        0.20990
     90B -> 93B       -0.42363
SavETr:  write IOETrn=   770 NScale= 10 NData=  16 NLR=1 LETran=      62.
The wavefunction has an internal instability.



请问需要用 alter关键字来修正波函数吗?
回复此楼
甲烷无敌
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

Jasminer

铁杆木虫 (著名写手)

stable=opt
2楼2012-11-16 19:45:30
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

t13340033021

金虫 (正式写手)

测试了,无效
甲烷无敌
3楼2012-11-17 14:36:40
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

t13340033021

金虫 (正式写手)

引用回帖:
2楼: Originally posted by Jasminer at 2012-11-16 19:45:30
stable=opt

测试了,无效
甲烷无敌
4楼2012-11-17 14:36:46
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

Jasminer

铁杆木虫 (著名写手)

【答案】应助回帖

★ ★ ★ ★ ★ ★ ★ ★
感谢参与,应助指数 +1
cenwanglai: 金币+3, 谢谢回复~ 2012-11-18 16:46:27
t13340033021: 金币+5 2012-11-18 17:36:52
我不知道你说的无效到底是什么情况。
如果自动优化波函数不能使之稳定的话,我不认为简单的alter能解决问题。
如果是体系中存在多个金属原子的磁性耦合,G09的分片段指定初始波函数应可以解决这一问题。
5楼2012-11-17 23:06:48
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

t13340033021

金虫 (正式写手)

引用回帖:
5楼: Originally posted by Jasminer at 2012-11-17 23:06:48
我不知道你说的无效到底是什么情况。
如果自动优化波函数不能使之稳定的话,我不认为简单的alter能解决问题。
如果是体系中存在多个金属原子的磁性耦合,G09的分片段指定初始波函数应可以解决这一问题。

opt的确可以拿到稳定的波函数,但计算量相当于收敛了几个电子组态。我想请教您opt是否相当于程序自动寻找了正确的自旋多重度
甲烷无敌
6楼2012-11-17 23:31:27
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

Jasminer

铁杆木虫 (著名写手)

引用回帖:
6楼: Originally posted by t13340033021 at 2012-11-17 23:31:27
opt的确可以拿到稳定的波函数,但计算量相当于收敛了几个电子组态。我想请教您opt是否相当于程序自动寻找了正确的自旋多重度...

应该是相当于必要而非充分的寻找了正确的基态多重度
7楼2012-11-18 16:25:59
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

t13340033021

金虫 (正式写手)

引用回帖:
7楼: Originally posted by Jasminer at 2012-11-18 16:25:59
应该是相当于必要而非充分的寻找了正确的基态多重度...

那么是否要通过收敛每个SCF根据最终能量判断多重度呢?
甲烷无敌
8楼2012-11-18 17:37:33
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

Jasminer

铁杆木虫 (著名写手)

引用回帖:
8楼: Originally posted by t13340033021 at 2012-11-18 17:37:33
那么是否要通过收敛每个SCF根据最终能量判断多重度呢?...

对于从实际图景考虑可能的几种多重度,如果能量计算准确的话,是的。
9楼2012-11-18 22:21:28
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 t13340033021 的主题更新
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[公派出国] 售SCI一区文章,我:8.O.551.O.5.4,科目全,可伽急 +3 wsYGJnURYvm6 2026-08-31 4/200 2026-09-01 06:38 by o3EMxHUAQDtS
[硕博家园] 售SCI一区文章,我:8O5.5.1.O5.4,科目全,可伽急 +3 wsYGJnURYvm6 2026-08-31 4/200 2026-09-01 06:07 by o3EMxHUAQDtS
[找工作] 售SCI一区文章,我:8.O.55.1.O.54,科目齐全,可伽急 +3 EZ67DG6CvcKL 2026-08-31 4/200 2026-09-01 06:06 by o3EMxHUAQDtS
[考博] 售SCI一区T0P文章,我:8O.55.1.O.54,科目全,可伽急 +3 EZ67DG6CvcKL 2026-08-31 5/250 2026-09-01 05:45 by o3EMxHUAQDtS
[基金申请] 投票:  有多少人是今天查系统知道结果的? +17 爱看书的可乐 2026-08-26 19/950 2026-08-31 21:30 by xiangy672
[基金申请] 哪位高人中了,把查询到的截图贴出来让我看看,让我长长见识 +6 yuleib84 2026-08-26 7/350 2026-08-31 19:46 by 鱼翔浅底1
[基金申请] 面上意见出来了 +12 黄鸟于飞Chao 2026-08-29 23/1150 2026-08-31 18:57 by 黄鸟于飞Chao
[基金申请] 能否申诉? +7 echo8914667 2026-08-30 8/400 2026-08-31 17:00 by yihongxu
[基金申请] 基金不中,共勉 +12 eulota 2026-08-26 12/600 2026-08-31 08:42 by ZJTJZ
[基金申请] 麻烦专家们看看评委们的意见(F口面上) +6 gdd2018 2026-08-28 11/550 2026-08-30 08:58 by 超级无敌华子
[基金申请] 国自然评审意见 +13 wangmingqi 2026-08-28 19/950 2026-08-29 10:22 by Poppy1104
[基金申请] 2026年叶企孙基金 +4 bud_bud 2026-08-27 7/350 2026-08-29 07:23 by foolishmani
[基金申请] 系统查不到 +11 董八千 2026-08-26 11/550 2026-08-28 18:06 by Leogzhya
[基金申请] 基金系统什么内容也没有 30+4 winsaint 2026-08-27 9/450 2026-08-28 11:06 by maolC
[基金申请] 看板上这么多中的,有点像50人群里49个人都是骗子的那种感觉…… +5 a089 2026-08-26 6/300 2026-08-27 14:05 by jonewore
[基金申请] 为什么 国际(地区)合作与交流项目 没有放榜? 10+3 majunge000 2026-08-26 11/550 2026-08-27 08:42 by 北京莱茵编辑
[基金申请] 我不理解! +15 Edward_pc 2026-08-26 23/1150 2026-08-26 20:34 by zzuzxg
[基金申请] 2026年8月25日国自然放榜前突然收到列入评审专家邮件,有关系吗? +25 木水思豆 2026-08-25 28/1400 2026-08-26 14:53 by draco1987
[基金申请] 国合里面能看到了 +7 一怀馨秋 2026-08-26 7/350 2026-08-26 11:23 by zhaosm1982
[基金申请] 某些机构,以效率低为荣,以效率低作为存在感 +9 yuleib84 2026-08-25 10/500 2026-08-25 17:14 by alexon
信息提示
请填处理意见