- UID
- 1
- 在线时间
- 3727 小时
- 注册时间
- 2014-10-7
- 最后登录
- 2019-10-21
- 好友
- 57
- 积分
- 46931
管理员
公社社长
|
详见：http://gaussian.com/relnotes/
G16 C.01相对于G16 B.01，我感觉没有什么太大改进，值得一提的也就是以下几点
Utilities can now take the -m command-line argument to specify the amount of memory available to the utility. For example: formchk -m=1gb myfile
即使用formchk、cubegen的时候，当处理大fch文件时，不是非得用GAUSS_MEMDEF环境变量才能指定可用内存量上限，现在用-m也可以指定了
V100 support is new with [REV C]
支持了更新的GPU加速计算，但是根据之前的官方测试数据看，G16的GPU加速基本没性价比可言（相对于Terachem来说）
[REV C] NBO version 7 is supported. There are new options to the Population keyword: Pop=NPA7, Pop=NBO7, Pop=NBO7Read and Pop=NBO7Delete request Natural Population Analysis, full Natural Bond Orbital Analysis, full NBO with NBO input read from the input stream and NBO analysis of the effects of deletion of some interactions (respectively), using NBO7 via the external interface. In addition, Pop=NEDA is used to perform Natural Energy Decomposition Analysis. The analysis uses the same input information about fragments as counterpoise calculations. Deletions and optimizations with deletions now work with either NBO6 or NBO7.
全面支持了NBO7，可以通过相应关键词直接挂接NBO7使用。但说实在的，NBO7相对于NBO6，实用性上改进非常有限，所以意义其实也不大，也就是NRT的算法得到了改进。这回引入了pop=NEDA关键词，可以直接调用NBO7做NEDA能量分解，看似弥补了Gaussian能量分解上的空白，但不过我也不觉得这有什么意义，因为NEDA本身是个非主流的能量分解方法，能量项定义的方式比较古怪，而且免费的PSI4的SAPT做得已经很成熟了（参看http://sobereva.com/469，以及http://sobereva.com/252的能量分解部分），GAMESS-US的LMO-EDA能量分解又是免费的，大体系能量分解还有Multiwfn的基于力场的能量分解可用（参看http://sobereva.com/442），所以在可遇见的未来，NEDA并不会流行。Gaussian开发者真正应该做的是直接原生支持一个类Morokuma的能量分解方法，这对开发者来说理应不是特别困难，但又是不少用户急需的功能。
Information on multidimensional relaxed scans is now stored on the formatted checkpoint file with details about the axes, rather than flattened, so these can be displayed in GaussView and other programs.
原先能面扫描任务只能用gview打开输出文件才能看，可能这回打开chk文件也能看了。
[REV C] The RESP (restrained electrostatic potential) constraint can be included in computing potential-derived charges. For example, Pop=(MK,Resp=N) applies a weight of N x 10-6 Hartrees to the squared charges. Other electrostatic potential-derived charge schemes also accept this option (e.g., CHelp, HLY). N defaults to 2.
外行人看了这个，可能觉得Gaussian支持了适合分子动力学模拟的RESP电荷计算了，觉得欣喜，其实完全不是这样。一般意义上的、适合分子动力学的真正的RESP电荷怎么算在《RESP拟合静电势电荷的原理以及在Multiwfn中的计算》（http://sobereva.com/441）里做了详细介绍，读者会看到计算RESP电荷相对于MK等普通拟合静电势电荷的过程绝不仅仅是引入一个惩罚函数这么简单，而C.01只不过引入了惩罚函数，这无异于NWChem、CP2K程序实现的那种所谓的RESP电荷，其实根本不适合柔性分子的分子动力学研究目的。
PS：给我的感觉，如今ORCA版本增加0.1带来的新功能，就能顶上Gaussian发布一个大版本（如G09->G16）的程度。但对于主流用户群体，ORCA要想超过Gaussian，还很遥远，至少得十年。
|
评分
-
查看全部评分
|
发表于 2019-7-25 19:46:37
发表于 2019-7-26 09:34:47
楼主
