计算化学公社

 找回密码 Forget password
 注册 Register
Views: 777|回复 Reply: 11
打印 Print 上一主题 Last thread 下一主题 Next thread

[Molpro] ccsd(t)计算时出现UNREASONABLE NORM错误

[复制链接 Copy URL]

108

帖子

0

威望

363

eV
积分
471

Level 3 能力者

跳转到指定楼层 Go to specific reply
楼主
本帖最后由 BinWang 于 2023-11-15 17:53 编辑

各位老师:
       大家好!我用Molpro中的ccsd(t)-F12A/AVDZ,VDZ-F12(用ccsd(t)/AVTZ,VTZ计算时出现相同错误)计算一个体系时出现如下错误:
UNREASONABLE NORM:  4.226133 CALCULATION WILL BE STOPPED

The threshold can be changed using option CC_NORM_MAX (see manual). The current value is CC_NORM_MAX= 0.20

?ERROR: NO CONVERGENCE. TRIPLES WILL NOT BE DONE. THIS CAN BE AVOIDED USING THE NOCHECK OPTION.


我把方法换成ccsd(t)-F12A/AVDZ或者ccsd(t)-F12A/VDZ-F12时出现同样错误。
请问CC_NORM_MAX= value的值从0.2改到4.3会有什么后果,结果还可信吗?有哪位老师知道引起这个问题的根本原因啊?怎么彻底的解决这个错误?
我已经贴上了输出文件,哪位老师有空帮我看一眼?提一些可靠的建议。


1.out

23.31 KB, 下载次数 Times of downloads: 5

输出文件

3811

帖子

4

威望

8000

eV
积分
11891

Level 6 (一方通行)

MOKIT开发者

2#
发表于 Post on 2023-11-15 10:09:34 | 只看该作者 Only view this author
看看是不是CCSD振幅迭代不收敛,T1诊断值是多少。
自动做多参考态计算的程序MOKIT

62

帖子

0

威望

1221

eV
积分
1283

Level 4 (黑子)

3#
发表于 Post on 2023-11-15 10:09:36 | 只看该作者 Only view this author
manual里写过cc_norm_max=50,可以直接使用

108

帖子

0

威望

363

eV
积分
471

Level 3 能力者

4#
 楼主 Author| 发表于 Post on 2023-11-15 10:35:57 | 只看该作者 Only view this author
zjxitcc 发表于 2023-11-15 10:09
看看是不是CCSD振幅迭代不收敛,T1诊断值是多少。

感谢老师的回复!计算任务没有正常结束,没有报出T1的值

108

帖子

0

威望

363

eV
积分
471

Level 3 能力者

5#
 楼主 Author| 发表于 Post on 2023-11-15 10:36:23 | 只看该作者 Only view this author
dlldwtz 发表于 2023-11-15 10:09
manual里写过cc_norm_max=50,可以直接使用

非常感谢老师的回复!

108

帖子

0

威望

363

eV
积分
471

Level 3 能力者

6#
 楼主 Author| 发表于 Post on 2023-11-15 17:44:54 | 只看该作者 Only view this author
zjxitcc 发表于 2023-11-15 10:09
看看是不是CCSD振幅迭代不收敛,T1诊断值是多少。

感谢老师的回复,确实是CCSD不收敛,输出文件没找出T1诊断信息。我贴上输出文件,麻烦您有时间帮忙看看什么原因,有什么别的办法? 1.out (23.31 KB, 下载次数 Times of downloads: 4)

108

帖子

0

威望

363

eV
积分
471

Level 3 能力者

7#
 楼主 Author| 发表于 Post on 2023-11-16 09:19:00 | 只看该作者 Only view this author
zjxitcc 发表于 2023-11-15 10:09
看看是不是CCSD振幅迭代不收敛,T1诊断值是多少。

老师您好!我把输出文件贴上去了,您有空能不能帮我看一下有什么办法啊,这个是什么原因造成我CCSD(T)结果不收敛的?http://bbs.keinsci.com/thread-41231-1-1.html

3811

帖子

4

威望

8000

eV
积分
11891

Level 6 (一方通行)

MOKIT开发者

8#
发表于 Post on 2023-11-16 17:50:57 | 只看该作者 Only view this author
本帖最后由 zjxitcc 于 2023-11-16 18:01 编辑
BinWang 发表于 2023-11-16 09:19
老师您好!我把输出文件贴上去了,您有空能不能帮我看一下有什么办法啊,这个是什么原因造成我CCSD(T)结 ...

不要改CC_NORM_MAX,不要写nocheck。你知不知道你写的wf,54,1,0是什么意思?
自动做多参考态计算的程序MOKIT

108

帖子

0

威望

363

eV
积分
471

Level 3 能力者

9#
 楼主 Author| 发表于 Post on 2023-11-22 16:51:25 | 只看该作者 Only view this author
zjxitcc 发表于 2023-11-16 17:50
不要改CC_NORM_MAX,不要写nocheck。你知不知道你写的wf,54,1,0是什么意思?

非常感谢老师的回复!WF指的是波函数的定义,54指的是54个电子,1是体系的不可约表示,0指的是自旋多重度2*S,与高斯中的定义2S+1有所不同

3811

帖子

4

威望

8000

eV
积分
11891

Level 6 (一方通行)

MOKIT开发者

10#
发表于 Post on 2023-11-22 17:18:10 | 只看该作者 Only view this author
BinWang 发表于 2023-11-22 16:51
非常感谢老师的回复!WF指的是波函数的定义,54指的是54个电子,1是体系的不可约表示,0指的是自旋多重度 ...

这个体系如果是中性的话,是98个电子。你说只有54个电子,就是在算正44+的阳离子,你说离谱不离谱。。
自动做多参考态计算的程序MOKIT

108

帖子

0

威望

363

eV
积分
471

Level 3 能力者

11#
 楼主 Author| 发表于 Post on 2023-11-22 21:51:54 | 只看该作者 Only view this author
zjxitcc 发表于 2023-11-22 17:18
这个体系如果是中性的话,是98个电子。你说只有54个电子,就是在算正44+的阳离子,你说离谱不离谱。。

,太惭愧了。
非常感谢老师的回复!

108

帖子

0

威望

363

eV
积分
471

Level 3 能力者

12#
 楼主 Author| 发表于 Post on 2023-11-22 21:53:05 | 只看该作者 Only view this author
BinWang 发表于 2023-11-22 21:51
,太惭愧了。
非常感谢老师的回复!


我纳闷了好久为啥算这么快

本版积分规则 Credits rule

手机版 Mobile version|北京科音自然科学研究中心 Beijing Kein Research Center for Natural Sciences|京公网安备 11010502035419号|计算化学公社 — 北京科音旗下高水平计算化学交流论坛 ( 京ICP备14038949号-1 )|网站地图

GMT+8, 2024-11-25 23:03 , Processed in 0.187956 second(s), 23 queries , Gzip On.

快速回复 返回顶部 返回列表 Return to list