计算化学公社

标题: 求助：有没有可以替代NBO分析求算二阶微扰能的方法 [打印本页]

作者
Author: Q-Chembio-llg 时间: 2020-4-6 11:00
标题: 求助：有没有可以替代NBO分析求算二阶微扰能的方法
有个问题想请教大家一下。就是我需要计算氢键donner-accepter那个二阶微扰能，但是由于种种原因没法用高斯做NBO分析，所以想问问大家有没有可以替代NBO分析计算这些数据的软件或方法以及具体的操作。谢谢大家！

作者
Author: hebrewsnabla 时间: 2020-4-6 11:14
本帖最后由 hebrewsnabla 于 2020-4-6 11:22 编辑

做NBO分析又不一定要用高斯……Multiwfn可以产生.47文件，然后直接用NBO程序就行
买个NBO应该没什么困难吧

作者
Author: snljty 时间: 2020-4-6 11:19
“各种原因”是你没买NBO还是没买Gaussian还是什么...

作者
Author: zjxitcc 时间: 2020-4-6 11:28
高斯里的NBO是比较老的版本。NBO自己也可以是独立的程序，还是新版的，无需依赖高斯，有输入文件就行。就算你不买，在这里找个pirate version也行，前提是你懂得找的技巧。

作者
Author: Q-Chembio-llg 时间: 2020-4-6 11:28

snljty 发表于 2020-4-6 11:19
“各种原因”是你没买NBO还是没买Gaussian还是什么...

唉一言难尽，我就是个本科生做做毕业论文而已，因为疫情原因没法去学校算所以只能用自己的笔记本跑，但是我装的那个G16不知道是盗版的还是其它问题，做NBO的时候L607老是报错，用sob老师说的那种产生47文件的方法它也报错，所以实在想不出法子了

作者
Author: Q-Chembio-llg 时间: 2020-4-6 11:29

hebrewsnabla 发表于 2020-4-6 11:14
做NBO分析又不一定要用高斯……Multiwfn可以产生.47文件，然后直接用NBO程序就行
买个NBO应该没什么困难吧 ...

大佬，怎么用multiwfn产生47文件啊？求指点

作者
Author: Q-Chembio-llg 时间: 2020-4-6 11:31

zjxitcc 发表于 2020-4-6 11:28
高斯里的NBO是比较老的版本。NBO自己也可以是独立的程序，还是新版的，无需依赖高斯，有输入文件就行。就算 ...

谢谢您的指点，我安装了sob老师上传的那个NBO5.0，只是现在我找不到方法产生47文件，我按sob老师给的那个在高斯里产生47文件的方法试了一下结果还是L607报错

作者
Author: zjxitcc 时间: 2020-4-6 11:37
本帖最后由 zjxitcc 于 2020-4-6 11:44 编辑

Q-Chembio-llg 发表于 2020-4-6 11:29
大佬，怎么用multiwfn产生47文件啊？求指点

这不用他人指点，你自己1分钟就会了。在multiwfn的pdf手册里搜一下.47不就行了。在Multiwfn里载入.fchk文件，输入100你就能看见.47选项了。
另外，如果是半经验波函数，要慎重。就算是个人笔记本，算个DFT下6-31G(d,p)的单点根本不是难事。

作者
Author: snljty 时间: 2020-4-6 11:41

Q-Chembio-llg 发表于 2020-4-6 11:28
唉一言难尽，我就是个本科生做做毕业论文而已，因为疫情原因没法去学校算所以只能用自己的笔记本跑，但是 ...

盗版又不是技术因素是法律因素，Gaussian自带的NBO 3.1不支持半经验的波函数做分析，不知道是不是这个问题。

作者
Author: hebrewsnabla 时间: 2020-4-6 11:47

Q-Chembio-llg 发表于 2020-4-6 11:28
唉一言难尽，我就是个本科生做做毕业论文而已，因为疫情原因没法去学校算所以只能用自己的笔记本跑，但是 ...

报错可以解决的嘛，
输入文件和报错贴出来看看

作者
Author: Q-Chembio-llg 时间: 2020-4-6 11:48

zjxitcc 发表于 2020-4-6 11:37
这不用他人指点，你自己1分钟就会了。在multiwfn的pdf手册里搜一下.47不就行了。在Multiwfn里载入.fchk文 ...

好的非常感谢您！是的我之前用b3lyp/6-31g*算过单点了，我的体系大概80个原子，一个单点大概14-15个小时左右，还可以接受。

作者
Author: Q-Chembio-llg 时间: 2020-4-6 11:49

snljty 发表于 2020-4-6 11:41
盗版又不是技术因素是法律因素，Gaussian自带的NBO 3.1不支持半经验的波函数做分析，不知道是不是这个问 ...

好的我换个方法试一下，谢谢您！

作者
Author: Q-Chembio-llg 时间: 2020-4-6 12:21

zjxitcc 发表于 2020-4-6 11:37
这不用他人指点，你自己1分钟就会了。在multiwfn的pdf手册里搜一下.47不就行了。在Multiwfn里载入.fchk文 ...

咳咳非常抱歉又来问您问题了，我产生了47文件然后用NBO5.0跑了一下，但是得到的文件中没有donner-accepter二阶微扰能的数据，不知道是哪里没有搞好，还烦请您帮我看看

作者
Author: hebrewsnabla 时间: 2020-4-6 15:05

Q-Chembio-llg 发表于 2020-4-6 12:21
咳咳非常抱歉又来问您问题了，我产生了47文件然后用NBO5.0跑了一下，但是得到的文件中没有donner-accepte ...

我看了一下，问题在于Multiwfn产生的.47文件没有fock矩阵，不能做E2PERT。用高斯产生.47就可以。

作者
Author: sobereva 时间: 2020-4-6 19:17
来来回回也不说用Gaussian自带的NBO到底出了什么错，能解决的问题想办法解决就完了，何必折腾
目前能给出E2的只有NBO这一个程序

作者
Author: Q-Chembio-llg 时间: 2020-4-6 19:32

sobereva 发表于 2020-4-6 19:17
来来回回也不说用Gaussian自带的NBO到底出了什么错，能解决的问题想办法解决就完了，何必折腾
目前能给出E ...

哦哦谢谢老师，其实就一直只是L607出错而已，只不过后来大家给的方法我试了之后都会遇到同样的报错所以一直在追问。抱歉打扰大家了，以后我一定控制发帖频率，还请大家原谅我刚入门无脑问题多哈哈。我的问题已经解决了，把sob老师给的程序删掉最后的$end之后就能正常运行输出正常的结果了。

作者
Author: Q-Chembio-llg 时间: 2020-4-6 19:32

hebrewsnabla 发表于 2020-4-6 15:05
我看了一下，问题在于Multiwfn产生的.47文件没有fock矩阵，不能做E2PERT。用高斯产生.47就可以。

好的谢谢您！我已经解决了，抱歉打扰您了

作者
Author: zjxitcc 时间: 2020-4-6 19:56

Q-Chembio-llg 发表于 2020-4-6 19:32
哦哦谢谢老师，其实就一直只是L607出错而已，只不过后来大家给的方法我试了之后都会遇到同样的报错所以一 ...

说报错的时候应该给出有效的报错信息，说Lxxx报错是不良习惯，应该改正之。

作者
Author: Q-Chembio-llg 时间: 2020-4-6 20:19

zjxitcc 发表于 2020-4-6 19:56
说报错的时候应该给出有效的报错信息，说Lxxx报错是不良习惯，应该改正之。

好的谢谢您纠错！

作者
Author: sobereva 时间: 2020-4-6 20:20

Q-Chembio-llg 发表于 2020-4-6 19:32
哦哦谢谢老师，其实就一直只是L607出错而已，只不过后来大家给的方法我试了之后都会遇到同样的报错所以一 ...

光说个L几几几出错对于解决问题通常没有丝毫帮助，最起码把最后20行内容贴出来
初学者解释不清楚就直接上传压缩后的输出文件，别人一看就知道了，省得还得猜

作者
Author: Q-Chembio-llg 时间: 2020-4-6 20:21

sobereva 发表于 2020-4-6 20:20
光说个L几几几出错对于解决问题通常没有丝毫帮助，最起码把最后20行内容贴出来
初学者解释不清楚就直接 ...

好的我知道了，谢谢sob老师！

作者
Author: Q-Chembio-llg 时间: 2020-4-10 11:36

hebrewsnabla 发表于 2020-4-6 15:05
我看了一下，问题在于Multiwfn产生的.47文件没有fock矩阵，不能做E2PERT。用高斯产生.47就可以。

好的非常感谢您！我已经用高斯解决了！

作者
Author: sobereva 时间: 2021-9-11 06:31

hebrewsnabla 发表于 2020-4-6 15:05
我看了一下，问题在于Multiwfn产生的.47文件没有fock矩阵，不能做E2PERT。用高斯产生.47就可以。

我今天更新了官网上的Multiwfn，输出的47文件可以直接包含Fock矩阵了

欢迎光临计算化学公社 (http://bbs.keinsci.com/)