24小时热门版块排行榜    

CyRhmU.jpeg
查看: 1350  |  回复: 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的回帖
相关版块跳转 我要订阅楼主 刘志凌 的主题更新
信息提示
请填处理意见