计算化学公社

标题: 请问在高斯输入文件中直接写hf/dzvp2会出错，怎么改？ [打印本页]

作者
Author: NaNaN 时间: 2020-5-14 14:50
标题: 请问在高斯输入文件中直接写hf/dzvp2会出错，怎么改？
请问各位，我在文献中看到用hf/dzvp2计算核磁，输入文件中直接写hf/dzvp2会出错，在高斯说明书中找，只能找到hf/dgdzvp2，请问二者有区别吗？在输入文件中应该如何正确表示hf/dzvp2呢？

作者
Author: sobereva 时间: 2020-5-14 15:22
不要盲目效仿文献，用这基组也没什么额外好处，HF算核磁是搞笑的，碰上内行人不挨批才怪
dzvp2在Gaussian里叫dgdzvp2，dg代表此基组是DGauss程序里定义的

作者
Author: NaNaN 时间: 2020-5-14 17:01
老师，那您的意思是用dgdzvp2和dzvp2算是一样的呗？我不是专业做核磁的，不太懂，不知道什么方法比较适合算P核磁，找了很多文献，在b3lyp/6-31G（d，p）优化的结构下，用过b972/def2tzvp、hf/6-311++G（2d，2p）、b972/pcSseg-1、hf/dgdzvp2、b3lyp/6-31G（d，p）计算核磁，算出的结果和实验值都对不上，不知道问题出在哪？算P核磁老师有建议的方法吗？

作者
Author: linq5203 时间: 2020-5-14 17:09

NaNaN 发表于 2020-5-14 17:01
老师，那您的意思是用dgdzvp2和dzvp2算是一样的呗？我不是专业做核磁的，不太懂，不知道什么方法比较适合算 ...

都对不上是差得多大，你该不会直接在GuassView里看的核磁吧...

作者
Author: snljty 时间: 2020-5-14 17:47
《谈谈如何又好又快地计算NMR化学位移》
http://bbs.keinsci.com/thread-4434-1-1.html

作者
Author: NaNaN 时间: 2020-5-15 10:28

sobereva 发表于 2020-5-14 15:22
不要盲目效仿文献，用这基组也没什么额外好处，HF算核磁是搞笑的，碰上内行人不挨批才怪
dzvp2在Gaussian ...

老师您好，昨天要回复您但发错了位置了，不知道您是否看到我的问题。您的意思是用dgdzvp2和dzvp2算出的的结果是一样的吗？我不是专业做核磁的，不太懂，不知道什么方法比较适合算P核磁，找了很多文献，在b3lyp/6-31G（d，p）优化的结构下，用过b972/def2tzvp、hf/6-311++G（2d，2p）、b972/pcSseg-1、hf/dgdzvp2、b3lyp/6-31G（d，p）计算核磁，算出的结果和实验值都对不上，用不同的方法和理想值相差5-10ppm左右，不知道问题出在哪？算P核磁老师有建议的方法吗？希望老师有空解答一下（老师写的NMR计算的博文我昨天看了，其中大多涉及H、C核磁，关于31P的较少，所以还不是很懂）

作者
Author: NaNaN 时间: 2020-5-15 10:28

linq5203 发表于 2020-5-14 17:09
都对不上是差得多大，你该不会直接在GuassView里看的核磁吧...

您好，大约差5-8ppm。服务器上跑完题，在输出文件log中看的结果

作者
Author: sobereva 时间: 2020-5-15 23:57

NaNaN 发表于 2020-5-15 10:28
老师您好，昨天要回复您但发错了位置了，不知道您是否看到我的问题。您的意思是用dgdzvp2和dzvp2算出的的 ...

一样

原理上用b972/pcSseg-1并没什么问题，起码比HF/DZVP2这种诡异的东西强

只要你不用标度法，系统性误差必定很大。差个5~10 ppm也没什么意外的

作者
Author: NaNaN 时间: 2020-5-18 15:25

sobereva 发表于 2020-5-15 23:57
一样

原理上用b972/pcSseg-1并没什么问题，起码比HF/DZVP2这种诡异的东西强

老师，算核磁适合用多大的模型呢？
如果用大模型的话，输入文件这样写# nmr=(giao,nomicro) oniom(hf/gen:uff) nosymm 交题出错，错在哪了呢？
还有老师博文中提到的标度法列举的都是小分子，适合用于分子筛算核磁吗？

作者
Author: thanhtam 时间: 2020-5-18 16:03
本帖最后由 thanhtam 于 2020-5-18 16:04 编辑

NaNaN 发表于 2020-5-18 15:25
老师，算核磁适合用多大的模型呢？
如果用大模型的话，输入文件这样写# nmr=(giao,nomicro) oniom(hf/ge ...

你不提供报错信息，别人怎么知道错在哪儿……

上面sob老师说过了，HF算核磁是搞笑的，还在用

“大模型”是多大？说清楚。

作者
Author: NaNaN 时间: 2020-5-18 17:57

thanhtam 发表于 2020-5-18 16:03
你不提供报错信息，别人怎么知道错在哪儿……

上面sob老师说过了，HF算核磁是搞笑的，还在用

这是216T大模型报错信息QPErr --- A syntax error was detected in the input line.
# nmr=(giao,nomicro) oniom(hf/genecp:uff
         '
Last state= "NMR2"
TCursr= 167612 LCursr=       12
Error termination via Lnk1e in /home/wzc/G09_Linux/g09/l1.exe 报错信息QPErr --- A syntax error was detected in the input line.
# nmr=(giao,nomicro) oniom(hf/genecp:uff
         '
Last state= "NMR2"
TCursr= 167612 LCursr=       12
Error termination via Lnk1e in /home/wzc/G09_Linux/g09/l1.exe
看见sob老师说的了，但是还没找到合适的方法，所以先用这个算算，看见合适的再改正

作者
Author: sobereva 时间: 2020-5-18 23:43

NaNaN 发表于 2020-5-18 17:57
这是216T大模型报错信息QPErr --- A syntax error was detected in the input line.
# nmr=(giao,nomi ...

nomicro怎么能写到nmr里，明明是opt里的选项
根本没必要用ONIOM，恰当弄个簇模型直接QM计算就完了，用ONIOM自找麻烦。起码NMR的时候用ONIOM绝对多余。
显然也可以用于分子筛

作者
Author: NaNaN 时间: 2020-5-20 08:36

sobereva 发表于 2020-5-18 23:43
nomicro怎么能写到nmr里，明明是opt里的选项
根本没必要用ONIOM，恰当弄个簇模型直接QM计算就完了，用ON ...

谢谢老师回复。还有两个问题，一是老师说的QM在输入文件中怎么表示呢？可以说的具体一点吗？（高斯说明书在学校，暂时不能查找），或者老师有相关博文吗？
二是老师关于算核磁的博文中，标度法参数汇总表中都是算H，C，N，B的核磁，没有P的相关参数，怎么将标度法应用到P的核磁呢？

作者
Author: zjxitcc 时间: 2020-5-20 09:02

NaNaN 发表于 2020-5-20 08:36
谢谢老师回复。还有两个问题，一是老师说的QM在输入文件中怎么表示呢？可以说的具体一点吗？（高斯说明书 ...

QM是个统称，这里是全量子计算、无需ONIOM的意思，不是让你写个QM关键词作为方法，你正常写泛函名称就行。

“高斯说明书在学校，暂时不能查找”不成立，说明书到处都有，例如查泛函名称，上官网gaussian.com/dft就行，网页命名规律就是gaussian.com/关键词信息，再如查看scf就上gaussian.com/scf。

作者
Author: sobereva 时间: 2020-5-21 07:16

NaNaN 发表于 2020-5-20 08:36
谢谢老师回复。还有两个问题，一是老师说的QM在输入文件中怎么表示呢？可以说的具体一点吗？（高斯说明书 ...

你搜不到P的参数就没法用，除非你自己找一批可靠的实验数据自己拟合一个

Gaussian电子版pdf手册论坛里就有，这里也有http://sobereva.com/g09。现在根本没人看打印版的

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