ACES
(含ACES II,ACES III,ACES II-MAB,CFOUR)


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。不要写为负数!