24小时热门版块排行榜    

查看: 1316  |  回复: 11

雾溪之魅

银虫 (小有名气)

[求助] 求助CCSD(T)问题 已有2人参与

大家好,我想请求一个问题,我用CCSD(T)计算单点能正常结束后,能量反而增加,更不准确了,是什么原因呢?
因为计算中冻芯出了问题,所有在关键词中加入了iop(8/11=1),请我是由于这个关键词造成的错误吗?
回复此楼

» 猜你喜欢

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

已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

枪下游魂

木虫 (著名写手)

【答案】应助回帖

★ ★
感谢参与,应助指数 +1
gkf高: 金币+2 2014-06-26 22:53:07
IOp(8/11)
MO coefficient, orbital energy, and number of electrons test.
0 Default, same as 2 except for during BD iterations or BD=Read.
1 Just print a warning message.
2 Kill the job if any MO coefficients are greater than 1000.0 or the smallest difference between occupied and virtual orbital energies is less than 0.001. Also, kill a frozen-core job if there is significant core-valence mixing in the canonical orbitals.
00 Default, same as 10.
10 Suppress such a test (CPHF may still be done for such a case).
20 Kill the job if there is no corr. energy; e.g., if there is only 1 electron or 1 virtual spin-orbital.

根据iop(8/11=1)的描述,我觉得有一定关系,可能你的冻芯出现了选项2后面所说的那种情况。

» 本帖已获得的红花(最新10朵)

2楼2014-06-26 08:28:28
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

雾溪之魅

银虫 (小有名气)

送红花一朵
引用回帖:
2楼: Originally posted by 枪下游魂 at 2014-06-26 08:28:28
IOp(8/11)
MO coefficient, orbital energy, and number of electrons test.
0 Default, same as 2 except for during BD iterations or BD=Read.
1 Just print a warning message.
2 Kill the job if any MO c ...

嗯,这样啊,可是我要是不加上iop(8/11=1)的话就会出现l801 错误,请问这个错误还有别的解决方法吗?非常感谢
3楼2014-06-26 08:58:19
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

雾溪之魅

银虫 (小有名气)

引用回帖:
2楼: Originally posted by 枪下游魂 at 2014-06-26 08:28:28
IOp(8/11)
MO coefficient, orbital energy, and number of electrons test.
0 Default, same as 2 except for during BD iterations or BD=Read.
1 Just print a warning message.
2 Kill the job if any MO c ...

你好,我刚刚还发现某些不加IOp(8/11)的CCSD(T)单点能结果也明显增大不对,这又是什么原因呢
4楼2014-06-26 09:02:56
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

枪下游魂

木虫 (著名写手)

引用回帖:
4楼: Originally posted by 雾溪之魅 at 2014-06-26 09:02:56
你好,我刚刚还发现某些不加IOp(8/11)的CCSD(T)单点能结果也明显增大不对,这又是什么原因呢...

小弟知识有限,还需要学习。。。CCSD太高大上,没仔细研究过。。。
我想问下你这个能量的参照物是什么?是实验值?还是低理论水平下的能量?能量校正的话是不是还应该加零点能?
5楼2014-06-26 09:25:28
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

雾溪之魅

银虫 (小有名气)

引用回帖:
5楼: Originally posted by 枪下游魂 at 2014-06-26 09:25:28
小弟知识有限,还需要学习。。。CCSD太高大上,没仔细研究过。。。
我想问下你这个能量的参照物是什么?是实验值?还是低理论水平下的能量?能量校正的话是不是还应该加零点能?...

加零点能的,是优化结构中的零点能
6楼2014-06-26 09:35:27
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

zhou2009

版主 (著名写手)

【答案】应助回帖

感谢参与,应助指数 +1
你用什么样的基组?
7楼2014-06-26 14:35:49
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

雾溪之魅

银虫 (小有名气)

引用回帖:
7楼: Originally posted by zhou2009 at 2014-06-26 14:35:49
你用什么样的基组?

Def2-TZVPPD
8楼2014-06-26 14:42:58
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

zhou2009

版主 (著名写手)

【答案】应助回帖

★ ★
gkf高: 金币+2 2014-06-26 22:53:15
Def2-tzvppd是赝势基组,将内层电子(core electrons)的基函数用赝势来代替以减少基函数的数目。你的问题可能与此相关。
9楼2014-06-26 17:25:40
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

雾溪之魅

银虫 (小有名气)

引用回帖:
9楼: Originally posted by zhou2009 at 2014-06-26 17:25:40
Def2-tzvppd是赝势基组,将内层电子(core electrons)的基函数用赝势来代替以减少基函数的数目。你的问题可能与此相关。

可是为什么用CCSD(T)单点会出现这种冻芯造成的能量增大,而用DFT来算的话不会呢?
那现在这个CCSD(T)J结果不能用了,应该怎么办呢?我还是需要CCSD(T)计算单点
10楼2014-06-26 17:31:17
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 雾溪之魅 的主题更新
信息提示
请填处理意见