请选择 进入手机版 | 继续访问电脑版

计算化学公社

 找回密码
 现在注册!
查看: 233|回复: 0

[Gaussian/gview] 发现g16的一个bug

[复制链接]

203

帖子

0

威望

422

eV
积分
625

Level 4 (黑子)

发表于 2019-1-10 16:25:02 | 显示全部楼层 |阅读模式
对于自旋非限制性计算,当让g16(测试版本是RevB.01)同时输出fchk文件并调用NBO6计算时,生成的fchk文件中没有Beta电子的能量和分子轨道系数部分。并且fchk文件开头居然显示的是ROHF计算(虽然输入文件已经指明是UHF或UDFT)!!!

一个简单的算例是:
  1. #P UHF/STO-3G fchk=All pop=NBO6

  2. t

  3. 0 2
  4. Li

复制代码
其它一些测试:

1. 不使用fchk=All而是用%chk=***得到chk文件,再用formchk工具转换得到的fchk文件也是缺少beta电子的部分。
2. 采用pop=NBO6Read(末尾添上不添上$NBO ... $END指令部分都无所谓),仍然不能解决问题。
3. 去掉pop=NBO6Read,或者采用g16内置的pop=NBORead,则会得到正确的fchk文件——fchk文件显示是UHF计算并且包含beta电子的信息。
4. 换成g09(版本D.01),调用相同版本的NBO6,则得到完全正常的fchk文件。 另外,g16生成的fchk文件中,“Alpha MO coefficients”数据之后会出现一块“Orthonormal basis”的部分,不知那些数据是什么意思(g09没有)。

不知道大家有没有什么高见?


您需要登录后才可以回帖 登录 | 现在注册!

本版积分规则

手机版|北京科音自然科学研究中心|京公网安备 11010502035419号|计算化学公社 — 北京科音旗下高水平计算化学交流论坛 ( 京ICP备14038949-1号 )

GMT+8, 2019-1-20 18:03 , Processed in 0.101340 second(s), 23 queries .

快速回复 返回顶部 返回列表