24小时热门版块排行榜    

查看: 1625  |  回复: 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的回帖
相关版块跳转 我要订阅楼主 雾溪之魅 的主题更新
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[硕博家园] 售SCI一区T0P文章,我:8.O.55.1.O54,科目全,可伽急 +3 23jxep3nCNZb 2026-09-14 3/150 2026-09-15 09:38 by VMDqgnlKimtZ
[公派出国] 售SCI文章,我:8O5.5.1.O.54,科目齐全,可+急 +3 6F5UbRU2I5hL 2026-09-14 3/150 2026-09-15 08:58 by C79jjtjAKjEn
[考研] 售SCI文章,我:8O5.5.1.O.54,科目齐全,可+急 +3 6F5UbRU2I5hL 2026-09-14 3/150 2026-09-15 08:14 by C79jjtjAKjEn
[考博] 售SCI-T0P文章,我:8O.5.5.1.O.54,科目齐全,可+急 +3 6F5UbRU2I5hL 2026-09-14 3/150 2026-09-15 07:58 by C79jjtjAKjEn
[考博] 售SCI一区T0P文章,我:8O.55.1.O.5.4,科目齐全,可+急 +5 s3fFTmArrBt6 2026-09-14 6/300 2026-09-15 07:38 by C79jjtjAKjEn
[考博] 售SCI一区文章,我:8O5.5.1.O5.4,科目全,可伽急 +6 s3fFTmArrBt6 2026-09-14 6/300 2026-09-15 07:14 by C79jjtjAKjEn
[找工作] 售SCI一区T0P文章,我:8.O.55.1.O54,科目全,可伽急 +3 3n8v2C8RimXI 2026-09-13 4/200 2026-09-15 06:01 by 5BDX0d0WFp7t
[博后之家] 售SCI一区T0P文章,我:8.O55.1.O.54,科目全,可十急 +6 3n8v2C8RimXI 2026-09-13 6/300 2026-09-15 05:50 by 5BDX0d0WFp7t
[教师之家] 售SCI一区T0P文章,我:8.O.55.1.O.54,科目齐全,可+急 +6 3n8v2C8RimXI 2026-09-13 6/300 2026-09-15 05:38 by 5BDX0d0WFp7t
[考博] 售SCI文章,我:8O.5.5.1O.54,科目全,可十急 +5 3n8v2C8RimXI 2026-09-13 5/250 2026-09-15 05:34 by 5BDX0d0WFp7t
[考研] 售SCI一区文章,我:8.O.551.O.5.4,科目全,可伽急 +5 QUjhNVAcOSff 2026-09-13 6/300 2026-09-15 05:13 by 5BDX0d0WFp7t
[公派出国] 售SCI一区T0P文章,我:8.O55.1.O.54,科目全,可十急 +6 QUjhNVAcOSff 2026-09-13 8/400 2026-09-15 05:01 by mibUvS8DDCwf
[考研] 售SCI一区文章,我:8.O.551.O.5.4,科目全,可伽急 +4 QUjhNVAcOSff 2026-09-13 4/200 2026-09-15 05:00 by mibUvS8DDCwf
[公派出国] 售SCI-T0P文章,我:8O.5.5.1.O.54,科目齐全,可+急 +3 LwdutQ8HoqWP 2026-09-13 6/300 2026-09-14 22:37 by vZfe6xYu34yj
[考研] 售SCI文章,我:8O5.5.1.O.54,科目齐全,可+急 +5 LwdutQ8HoqWP 2026-09-13 5/250 2026-09-14 22:25 by vZfe6xYu34yj
[教师之家] 售SCI一区文章,我:8.O.55.1.O.54,科目齐全,可伽急 +3 LwdutQ8HoqWP 2026-09-13 4/200 2026-09-14 16:13 by Aj1rhIDL5ixY
[基金申请] 求教各位大神:2026教育部人文社科青年基金项目何时公示呀? +5 云宴山人 2026-09-10 7/350 2026-09-14 12:11 by ecnu2013
[找工作] 浙江师范大学是怎么坑我的 +6 recruit123 2026-09-09 6/300 2026-09-14 12:00 by 萨仁其其格
[考博] 西北工业大学材料学院电化学传感与催化课题组招收2027年推免研究生、博士生、博士后 +3 马方园 2026-09-10 4/200 2026-09-14 09:29 by Chem张zz
[考博] 27秋季材料博士申请求助 10+3 Cahal1 2026-09-11 4/200 2026-09-14 08:50 by 北京莱茵编辑
信息提示
请填处理意见