美国Gaussian公司于2019年7月24日发布了Gaussian16 Revision C.01; GaussView6.1.1; TCP-Linda9.2; GMMX3.1版本。

Gaussian16 C.01版本的Linda版需要将Linda 升级至 TCP-Linda 9.2。

根据美国Gaussian公司的规定：2019年4月1日之后购买 Gaussian 16 Revision B.01, TCP-Linda 9.1,GMMX 3.0, GaussView 6.0.16 软件的用户，将在几周内收到美国Gaussian公司寄出的免费升级软件。

以下是Gaussian 16 Revision C.01 Release Notes，方便大家了解。

新模拟方法

·        [REV C] 支持NBO 7软件接口。此版本有许多新的选项：Keyword: Pop=NPA7, Pop=NBO7, Pop=NBO7Read 及 Pop=NBO7Delete。分别代表需要计算自然布居分析(Natural Population Analysis)，全自然键轨道分析(full Natural Bond Orbital Analysis)，从输入文件读取信息的全NBO分析以及考虑剔除某些作用影响之后的NBO分析。此外，Pop=NEDA 可用于分析自然能量分解分析(Natural Energy Decomposition Analysis)。此功能所使用的分子碎片输入文件信息与计算BSSE校正时一样.

·        [REV C] The RESP (restrained electrostatic potential) 电荷可以在计算potential-derived charges时使用。例如, Pop=(MK,Resp=N) 将 N x 10^-6 Hartrees电荷施加于平方电荷. 其它静电势导出的电荷方法都可以使用此选项。(如 CHelp, HLY). N 默认值为2.

·        [REV C] Pop=SaveHirshfeld 及 Pop=SaveCM5 可以保存特定形式的电荷到MM电荷中，以便于接下来的计算。

性能提升

·        Hartree-Fock及DFT计算在Linux系统下支持V100(Volta, [REV C]), P100 (Pascal, [REV B]), NVIDIA K40及K80 GPUs，上述GPUs显卡计算性能都得到了提升。具体请参考 Using GPUs

用法提升

·        [REV C] 每个振动模式的ROA不变量现在只由G16或freqchk输出(如果指定了正常模式），而不是默认输出。

·        [REV C] Utilities 模块中现在能使用 -m 参数指定所使用内存大小。例如:

formchk -m=1gb myfile -m 选项必须写在文件名或者其它参数前。

·        [REV C]  Link 0命令中%SSH 及其等效命令可用于启动Linda workers, 而不是用rsh 或 ssh。

·        [REV C] 当指定 Geom=AllCheck时一些默认项可以重新设置:

·        Field=NoChk 可用于禁止从chk文件读取外场系数。

·        Geom=GenConnectivity 强制重新计算原子连接关系，而不是使用chk文件中的信息。

·        Geom=UseStandardOrientation 使用chk文件中的标准坐标作为新作业的坐标

·        [REV C] 稳态结构优化中的一些默认项可以重新设置:

·        Opt=NGoUp=N 在只进行线性搜索之前，允许能量增加N倍。默认值为1（在第二次能量增加后只执行线性搜索）；n=-1表示在能量增加时强制执行线性搜索。

·        当接近鞍点时，Opt=NGoDown=N使程序将Hessian的最多N个特征向量与负特征值混合，形成远离鞍点的一步。默认值为3；N=-1表示关闭此功能，并且算法只执行常规的RFO步骤。

·        Opt=MaxEStep=N 表示从鞍点离开时，下一步要移动的长度为N/1000 (Bohr or radians)。对于常规优化，默认值为N=600 (0.6)，对于ONIOMOpt=Quadmac 计算，默认值为N=100 (0.1)。

·        [REV C] 多维势能面扫描的信息现在存储在fchk文件中，其中包含有关轴的详细信息，因此这些信息可以显示在GaussView和其他程序中。

·        [REV C] 本版本将在chk文件中存储和检查版本号。这可以避免由于不同版本chk文件导致的各种复杂问题。c8616程序可用于更新chk文件。并且-fixver选项可以 unfchk，即使输入格式的chk文件中没有版本，它也会自动创建一个当前版本。

·        [REV C] 支持使用4字节或8字节整数的原始二进制文件。默认为4字节(NEC系统除外)。此功能支持包括Output 关键词和formchk程序的新选项、新的Link 0命令以及新的命令行选项和环境变量。

·        [REV C] 本版本将向矩阵元素文件中添加有关ONIOM层、优化和轨迹结果的信息。它还为Output 关键词增加了新的选项，包括AO双电子积分、重叠导数、核心哈密顿量和其他矩阵和/或AO 双电子积分导数。

·        [REV C] 脚本功能增强：

·        ONIOM方法中当使用外部程序/脚本时，外部关键字(External keyword)的AllAtoms 及 ActiveAtoms用于提供所有原子的信息，或仅提供模型（高层）中的原子的信息。

·        $g16root/g16/bsd/inp2mat脚本接受Gaussian输入文件并生成一个矩阵元素文件，其中包含输入文件（坐标、基组等）所隐含的信息，而不进行完整的计算。这是由GauOpen中的python接口用来将这些信息导入到矩阵元文件中，但也可以在其他脚本中使用，以避免需要解析Gaussian输入文件。

·        testrt 工具现在输出G16使用的整数大小，以便脚本可以检查在矩阵元文件中默认使用的整数大小。

Bug修复

以下内容在C.01版本中得到修复:

·        在计算多个入射频率的Raman 或 ROA时，使用Freq=Anharmonic产生的问题已修复。

·        高角动量和纯DFT泛函计算时的内存调用问题得到修复。

·        合理地跳过了DFTB参数文件中的内容。

·        早期中止的任务chk文件使用chkchk遇到的问题已解决。

·        大分子体系中使用PM7R6方法的杂化项问题已经得到解决。

·        将GVB代码中被占轨道数的限制提高到1000，解决了大分子FMM和GVB的一些问题。

·        Grimme (D2或D3)色散校正与鬼原子(ghost atoms)的问题得到解决。

·        Punch=MO 及chkchk -p输出的轨道能量问题已经得到解决。

·        扩展名 -fck= (/fck= on Windows systems) 命令行参数问题已经修复，默认为.fck ，但 .fchk 仍然有用。

·        在一般基组输入中指定命名基组时，以前未检测到的一些错误现在可以被识别出来。

·        解决了在ONIOM模型计算时在文件上运行 formchk的问题，或在计算期间停止的ONIOM作业的chk文件上运行formchk的问题。

·        修复了.fchk (或.fch)文件中MM参数值不完整或错误的问题。

·        Field=Read在使用密度拟合时，不再尝试读取两次外场值。

·        CIS方法中的 SCF=Conventional 选项产生的错误默认值已经得到修复。

·        鬼原子(ghost atoms)存在时使用Guess=Read产生的问题已经得到修复。

·        Punch=GAMESS 现在支持 H 及更高的函数。

·        在输入文件中的原子说明部分，又可以使用-2代表Tv 平移矢量。

·        解决了原子长链分子生成内坐标的一些问题。

·        当双电子积分阈值很小时，计算单电子导数的问题已修复。

·        当使用非默认后自洽场方法(如 MP2=FreezeG2)进行Opt+Freq计算时，不能进行频率计算的问题已经得到解决。

·        对incore内存需求的低估问题（这可能导致作业默认为incore，然后导致内存耗尽）被修复。