1. | 问: | 在Post-HF计算中,为什么CPU的使用率很低? |
| 答: | |
| | 加上选项ABCDTYPE=AOBASIS。对于有些耦合簇计算方法(见手册),定义CC_PROGRAM=1还可以进一步减少硬盘读写,提高CPU使用率。加上这些选项后,通常能把Post-HF计算速度提高几倍。 |
2. | 问: | 程序说输入文件有错,但却总是找不到输入错误。 |
| 答: | |
| | 各个版本的ACES程序都有一些不好的输入约定。如果不遵守这些约定,程序在读输入文件的时候会导致错误。例如: 输入文件结尾必须有多个空行(大概三四个) 元素符号、坐标变量一律用大写字母(程序对此进行检查) 直角坐标中,元素符号前不能有空格 Z矩阵输入中,不能连用多个空格(程序对此进行检查,因此不要试图通过插入多个空格来使Z矩阵变得整齐),也不接受带负号的Z矩阵变量 当自定义基组时,原子在基组、ECP中的顺序必须和坐标中一致(程序对此不进行检查!!!) ECPDATA中的赝势函数类型(如p,s-p,p-d,等)必须用小写字母,否则读IIII积分时会出现莫名其妙的错误 CFour 2.0beta版bug:CH4这类Td点群的分子,HCH键角一定要用109.47122度,有效数字多了会导致对称性错误。另外,有些高对称分子会出现莫名其妙的对称性错误,即便关闭对称性也无法解决,必须改用1.0版。在以后版本的ACES III/Cfour中,也许能避免这些问题。 |
3. | 问: | 计算结果明显不合理(或不收敛),但却找不到输入错误。 |
| 答: | |
| | 1. 默认的占据不合理。需要指定占据数。 2. 由程序缺陷导致。各个版本ACES程序的HF模块都有个缺陷,在计算过渡金属体系,或者使用了3-zeta以上的大基组时,经常会得到错误的HF结果。当然此后的post-HF计算当然也就不对了。错误的征兆可能有: 1) HF不收敛; 2) HF收敛但能量极其不合理; 3) HF打印的HOMO、LUMO、芯轨道等轨道能量不合理; 4) 在打印的HF轨道中,FULLSYM和COMPSYM不可约表示不符合对称关系(如线形分子Delta轨道,在FULLSYM可能会错误地显示为两个SG+)。 解决方法是为HF加入阻尼选项: ACES II Florida版用:DAMP_TYP=DAVIDSON ACES II MAB版和CFour用:SCF_EXPSTART=6,SCF_EXPORDER=8,SCF_DAMP=500 个别情况下,SCF_DAMP还需要改为600、700等。但即便如此,也不能100%地保证HF就是对的。如果怀疑有问题,最好用其他从头算程序,用同样的基组对同样的结构,在HF级别做个单点计算进行验证。 3. 当自定义基组时,还有另一种可能: 程序要求原子在基组、ECP部分出现的顺序(扣除虚原子后),必须和坐标中一致,程序对此不进行检查(见上)。如果不一致,计算结果当然是错的。 |
4. | 问: | (Cfour)在使用ECP的情况下,解析频率计算得到奇怪的结果。 |
| 答: | |
| | ECP二阶导数计算目前还不支持对称性,因此需要关闭对称性(SYM=OFF),或者用解析梯度做数值差分(VIB=FINDIF)。 |
5. | 问: | (Cfour) 如何通过交换/转动轨道来调整占据? |
| 答: | |
| | 截止到2.0beta版还没有这个功能,CFour猜出来是啥就是啥。 |
6. | 问: | (CFour)用DROPMO指定post-HF冻结轨道的时候,如何知道HF的轨道编号? |
| 答: | |
| | HF轨道有三种编号系统: (1) 轨道能量从低到高编号; (2) 按照不可约表示分别编号; (3) 同(2),但占据轨道和非占据轨道分别编号。 在输入文件中冻结轨道的时候使用(1),但输出文件可能会打印(2)或(3),不要混淆。 |
7. | 问: | (Cfour)非谐振计算如何指定特殊基组? |
| 答: | |
| | 非谐振计算必须用特殊的关键词%basis让非谐模块定位基组信息的位置。例如,HF: CCSD ANHAM F H 1 R1
R1 = 1.02
*ACES2(CALC=CCSD,REFERENCE=RHF VIB=ANALYTIC,ANHARM=VPT2,ANH_ALGORITH=PARALLEL BASIS=SPECIAL FROZEN_CORE=OFF,ABCDTYPE=AOBASIS,CC_PROGRAM=ECC MEM_UNIT=GB,MEMORY=8)
F:3-21G H:3-21G
%anh_keywords CALC=CCSD,REFERENCE=RHF BASIS=SPECIAL FROZEN_CORE=OFF,ABCDTYPE=AOBASIS,CC_PROGRAM=ECC MEM_UNIT=GB,MEMORY=8
%basis F:3-21G H:3-21G |
类似的特殊关键词还有%ecp,%masses,%isotopes,参见源代码cubic/dumpgeoms2.f。
注意: (1) 写ECP信息的 输出格式有BUG(cubic/wrbasis.f),需要在do循环之前打印一空行。 (2) 非谐性质计算在使用数值Hessian时存在问题,必须使用解析Hessian。若同时使用了ECP,由于程序的BUG还必须关闭对称性。 |
8. | 问: | (CFour)在解析频率计算中,耦合微扰方程不收敛。 |
| 答: | |
| | 增加求解CPHF方程的迭代次数,并适当降低收敛阈值,如: CPHF_MAXCYC=200,CPHF_CONVER=8 注意:8表示1.0e-8。不要写为负数! |