计算化学公社

标题: gauss03 使用赝势计算金属，NBO计算出错，怎么解决？ [打印本页]

作者
Author: chunhuaqiushi 时间: 2014-11-13 14:28
标题: gauss03 使用赝势计算金属，NBO计算出错，怎么解决？
我在用gauss03计算一个团簇（Cs）2（NO3）3时，如果用lanl2dz赝势的话，不能使用pop=nbo，一使用就出现错误“Subroutine NAOANL could not find a d-type valence orbital on atom Cs 9。”如果不使用这个关键词，优化是可以进行的，但我想知道NBO分布。我尝试了几种元素，发现有的可以计算，有的不行，例如Fe，Co，Ru都可以，但是K，Ca，Cs都不行。我发现计算Fe等元素时，不使用l607，而不能计算的元素都是要使用l607的。我从网上查到说可能是要使用NBO5.9，但我不知道gauss03是否可以。以前使用Multiwfn绘图时，是直接读取的gauss的结果，我不知道Multiwfn能不能在gauss优化结果基础上直接计算BNO。请各位指点一下。

作者
Author: sobereva 时间: 2014-11-13 15:03
这是高斯自带的NBO3.1的bug，用NBO5及之后版本可以解决（必须是那种作为link嵌入到高斯的版本）
Multiwfn进行任何与NBO有关的分析、绘图，都需要读取NBO程序的输出，而没法独立去做。

你的出错原因如提示所示，想找Cs的d基函数但是当前的基组又没有。绕过这个bug的方法就是用对Cs有d函数的赝势基组，比如def2-TZVP，或者自行给当前的基组随便加个d函数。

作者
Author: chunhuaqiushi 时间: 2014-11-13 15:57
是这样啊！我看multiwfn的手册，也只是看到怎么利用gauss结果绘图，所以想给您确认一下。谢谢您的指导，顺便说一句，您不去小木虫，是它们的损失。真的，我原先在小木虫发的求助贴就是您回复的，但我当时根本没想到居然是您给我回复的，对我帮助很大，再次谢谢您了！希望这个论坛越办越好！

作者
Author: captain 时间: 2020-5-5 22:40

sobereva 发表于 2014-11-13 15:03
这是高斯自带的NBO3.1的bug，用NBO5及之后版本可以解决（必须是那种作为link嵌入到高斯的版本）
Multiwfn ...

请问大神
我也遇到了类似的问题
用的是def2-tzvp基组
NPA分析提示Subroutine NAOANL could not find a f-type valence orbital on atom Cs
按照您说的，自行给当前的基组随便加个f函数，比如把Cs原子def2-qzvp基组中的f函数拿过来用，这样合理么？

作者
Author: sobereva 时间: 2020-5-6 08:48

captain 发表于 2020-5-5 22:40
请问大神
我也遇到了类似的问题
用的是def2-tzvp基组

可以

作者
Author: captain 时间: 2020-5-6 09:19

sobereva 发表于 2020-5-6 08:48
可以

好的，谢谢大神！

作者
Author: LQY 时间: 2021-5-15 15:03

sobereva 发表于 2014-11-13 15:03
这是高斯自带的NBO3.1的bug，用NBO5及之后版本可以解决（必须是那种作为link嵌入到高斯的版本）
Multiwfn ...

老师，您好，这个方法对结果有什么影响吗？添加这个轨道基函数有没有限制，还是说随便找一个对同一种元素有定义的就行？

作者
Author: sobereva 时间: 2021-5-15 17:18

LQY 发表于 2021-5-15 15:03
老师，您好，这个方法对结果有什么影响吗？添加这个轨道基函数有没有限制，还是说随便找一个对同一种元素 ...

基函数变了自然有影响，但对于NBO分析影响甚微
没有什么所谓的限制

作者
Author: LQY 时间: 2021-5-17 11:44

sobereva 发表于 2021-5-15 17:18
基函数变了自然有影响，但对于NBO分析影响甚微
没有什么所谓的限制

谢谢老师指点

作者
Author: zhjie221dft 时间: 2021-10-7 15:14
老师，我最近做NBO计算，分子中有Cs和Pd, 本来使用的赝势基组是LANL2DZ，如果给Cs自定义加了d,f函数后，那么对Cs定义的赝势也需要改吗？还是Cs的赝势仍然定义为LANL2DZ？谢谢

作者
Author: sobereva 时间: 2021-10-17 07:04

zhjie221dft 发表于 2021-10-7 15:14
老师，我最近做NBO计算，分子中有Cs和Pd, 本来使用的赝势基组是LANL2DZ，如果给Cs自定义加了d,f函数后，那 ...

直接用def2-SVP或更好的def2-TZVP，构建得好得多，没必要自己给一个破基组手动加极化。
给价层基组加额外的极化不用改赝势

作者
Author: zhjie221dft 时间: 2021-10-26 18:06

sobereva 发表于 2021-10-17 07:04
直接用def2-SVP或更好的def2-TZVP，构建得好得多，没必要自己给一个破基组手动加极化。
给价层基组加额 ...

谢谢老师！

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