24小时热门版块排行榜    

查看: 1628  |  回复: 4

刘志凌

木虫 (小有名气)

[求助] EOMCCSD算激发态,能量为负值

EOMCCSD算激发态,能量为什么是负值????
部分输出结果:
Excitation energies and oscillator strengths:

----------------------------------------------
Excited State   1:      ?-B1G    -0.8781 eV    -1411.99 nm  f=-0.0000
Right Eigenvector
Beta  Singles Amplitudes
     I    SymI    A    SymA    Value
     8      1    20      2   -0.435558
     8      1    25      2    0.470537
     8      1    42      2    0.166284
     9      1    20      2    0.486067
     9      1    25      2   -0.505780
     9      1    42      2   -0.160251
Left  Eigenvector
Beta  Singles Amplitudes
     I    SymI    A    SymA    Value
     8      1    20      2   -0.425963
     8      1    25      2    0.459632
     8      1    42      2    0.161824
     9      1    20      2    0.491642
     9      1    25      2   -0.510968
     9      1    42      2   -0.160652
Total Energy, E(EOM-CCSD) =  -198.159381604   

----------------------------------------------
Excited State   2:      ?-B3G    -0.1147 eV    -10813.36 nm  f=0.0000
Right Eigenvector
Beta  Singles Amplitudes
     I    SymI    A    SymA    Value
     7      3    20      2   -0.648717
     7      3    25      2    0.699469
     7      3    42      2    0.245475
Left  Eigenvector
Beta  Singles Amplitudes
     I    SymI    A    SymA    Value
     7      3    20      2   -0.647863
     7      3    25      2    0.698324
     7      3    42      2    0.244545

----------------------------------------------
Excited State   3:      ?-B2G    -0.1147 eV    -10813.36 nm  f=0.0000
Right Eigenvector
Beta  Singles Amplitudes
     I    SymI    A    SymA    Value
     6      4    20      2    0.648717
     6      4    25      2   -0.699469
     6      4    42      2   -0.245475
Left  Eigenvector
Beta  Singles Amplitudes
     I    SymI    A    SymA    Value
     6      4    20      2    0.647863
     6      4    25      2   -0.698324
     6      4    42      2   -0.244545

----------------------------------------------
Excited State   4:      ?-B1G    0.0087 eV 141761.73 nm  f=0.0000
Right Eigenvector
Beta  Singles Amplitudes
     I    SymI    A    SymA    Value
     5      1    20      2    0.642877
     5      1    25      2   -0.702657
     5      1    42      2   -0.257430
Left  Eigenvector
Beta  Singles Amplitudes
     I    SymI    A    SymA    Value
     5      1    20      2    0.641932
     5      1    25      2   -0.701454
     5      1    42      2   -0.256534
Ground to excited state transition densities written to RWF  633
SavETr:  write IOETrn=   770 NScale= 10 NData=  16 NLR=2 LETran=     146.

S**2, projected HF & approx projected MPn energies after annihilation of
unwanted spin states (see manual for definitions):

spins       (S**2,0) (S**2,1)    PUHF        PMP2        PMP3        PMP4
annihilated
s+1         0.75000  0.75000 -197.308031 -198.175461 -198.089337
s+1,s+2     0.75000  0.75000 -197.308028 -198.175459 -198.089335
s+1 to s+3  0.75000  0.75000 -197.308028 -198.175459 -198.089335
s+1 to s+4  0.75000  0.75000 -197.308028 -198.175459 -198.089335
s+1 to s+5  0.75000  0.75000 -197.308028
s+1 to s+6  0.75000  0.75000 -197.308028


输入:
%mem=20GB
%NProcShared=8
#P EOMCCSD(maxcyc=300,root=1,50-50,Nstates=4)/genecp scf(tight,maxcycl
e=600)

。。。。。。
回复此楼

» 猜你喜欢

» 本主题相关价值贴推荐,对您同样有帮助:

好人一生平安
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

coolrainbow

木虫 (著名写手)

未来国家冻凉

【答案】应助回帖

★ ★
gmy1990(金币+2): 非常感谢你对论坛的贡献,奖励一下! 2011-06-12 12:28:03
刘志凌(金币+1): 我是09,不是03,03没有EOMCCSD 2011-06-12 23:52:58
如果你是0 3 的话,就会出现这种情况
技术博客:http://hi.baidu.com/coolrainbow/blog
2楼2011-06-12 12:02:28
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

beefly

专家顾问 (职业作家)

地沟油冶炼专家

1. 说明你指定的参考态为开壳层激发态。

2. 50-50对开壳层计算没有意义。
beefly《西太平洋大学现代英汉词典》[bi:fli]牛肉一般地
3楼2011-06-12 22:43:05
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

coolrainbow

木虫 (著名写手)

未来国家冻凉

【答案】应助回帖

刘志凌(金币+1): 问过他人,用stable=opt,解决了 2011-06-15 09:26:18
引用回帖:
Originally posted by coolrainbow at 2011-06-12 12:02:28:
如果你是0 3 的话,就会出现这种情况

我说的是

分子说明0 3

就会出现这种情况
技术博客:http://hi.baidu.com/coolrainbow/blog
4楼2011-06-13 09:17:33
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

beefly

专家顾问 (职业作家)

地沟油冶炼专家

先把比eom-cc便宜的cis、tddft弄懂,得到正常的结果,再做eom-cc计算。

这些方法都是先做一个参考态的hf或dft计算,参考态通常是基态,或闭壳层激发态。激发态通过线性响应理论得到。先检查你的参考态是如何定义的。
beefly《西太平洋大学现代英汉词典》[bi:fli]牛肉一般地
5楼2011-06-13 11:51:58
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 刘志凌 的主题更新
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考博] 售SCI一区T0P文章,我:8O.55.1.O.54,科目全,可伽急 +3 EZ67DG6CvcKL 2026-08-31 3/150 2026-08-31 22:14 by o3EMxHUAQDtS
[考研] 售SCI文章,我:8O.5.5.1O.54,科目全,可十急 +3 EZ67DG6CvcKL 2026-08-31 3/150 2026-08-31 22:12 by o3EMxHUAQDtS
[基金申请] 投票:  有多少人是今天查系统知道结果的? +17 爱看书的可乐 2026-08-26 19/950 2026-08-31 21:30 by xiangy672
[基金申请] 基金未中,这种答复是模板吗? +6 zhaosm1982 2026-08-27 7/350 2026-08-31 21:18 by qdxxmc
[基金申请] 面上意见出来了 +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
[基金申请] 29号明天会评吗 +4 笨笨唐 2026-08-28 4/200 2026-08-31 09:30 by huixian257
[基金申请] 为什么到现在没收到通知? +5 tannykie 2026-08-29 5/250 2026-08-30 21:05 by purplejack
[基金申请] 有没有仍没收到信息的 +7 德尚中行 2026-08-27 8/400 2026-08-30 20:52 by purplejack
[基金申请] 麻烦专家们看看评委们的意见(F口面上) +6 gdd2018 2026-08-28 11/550 2026-08-30 08:58 by 超级无敌华子
[基金申请] 国自然面上复盘~欢迎讨论 (金币+15) +15 晴天加油 2026-08-26 16/800 2026-08-29 18:28 by symmetry
[基金申请] 国自然评审意见 +13 wangmingqi 2026-08-28 19/950 2026-08-29 10:22 by Poppy1104
[基金申请] 看板上这么多中的,有点像50人群里49个人都是骗子的那种感觉…… +5 a089 2026-08-26 6/300 2026-08-27 14:05 by jonewore
[基金申请] 我不理解! +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
[基金申请] 出来了 +9 trojank 2026-08-26 9/450 2026-08-26 14:25 by 宝贝虫子
[基金申请] 为什么国自然不能直接公布 +4 bjdxyxy 2026-08-26 4/200 2026-08-26 13:12 by qingmu1201
[基金申请] 国际合作可查了,中了面上 (EPI+1)(金币+50) +18 Ldrop2023 2026-08-26 18/900 2026-08-26 11:15 by cmrandy
[基金申请] 在坚冰还盖着北海的时候,我看到了怒放的梅花。 (金币+10) +6 ziyangfang 2026-08-25 9/450 2026-08-25 20:26 by huagongfeihu
[基金申请] 某些机构,以效率低为荣,以效率低作为存在感 +9 yuleib84 2026-08-25 10/500 2026-08-25 17:14 by alexon
信息提示
请填处理意见