计算化学公社

 找回密码 Forget password
 注册 Register
Views: 589|回复 Reply: 1
打印 Print 上一主题 Last thread 下一主题 Next thread

[综合交流] Gaussian中的ONIOM(Gaussian:XTB)金属体系报错

[复制链接 Copy URL]

59

帖子

0

威望

962

eV
积分
1021

Level 4 (黑子)

各位老师好,我用Gaussian 中的ONIOM(Gaussian:XTB)计算金属体系,采用的接口是SOB老师的脚本xtb.sh ,但是金属的赝势基组读取出错,请问各位老师如何解决。
我的输入文件如下

%nosave
%nproc=20
%chk=$GAUSS_SCRDIR/$JOB_NAME.chk
%mem=100GB
#p  ONIOM(PBE1PBE/genecp: External="sh /scr1/xtb_6_6_1/xtb_interface/xtb.sh" )=Embedcharge  Opt(nomicro)

Title Card Required

0 1
C                0   14.47300000   14.48100000    8.83400000 H
C                0   12.97800000   14.63400000    8.97400000 H
.
.
Pt               0   14.31600000   14.56800000   12.76200000 H

Pt 0
LANL2DZ
****
H C O N 0
6-31G(d)
****

Pt 0
LANL2DZ


报错信息如下



5万

帖子

99

威望

5万

eV
积分
112462

管理员

公社社长

2#
发表于 Post on 2023-10-18 04:10:28 | 只看该作者 Only view this author
下文提到
将Gaussian与Grimme的xtb程序联用搜索过渡态、产生IRC、做振动分析
http://sobereva.com/421http://bbs.keinsci.com/thread-10106-1-1.html
如果你要用本文的做法计算很重的元素,比如Pt,有一点主要注意,即虽然Gaussian与xtb联用时Gaussian自己并不会去做任何量化计算,但是由于你没写基组,Gaussian程序内部默认当前设的基组是STO-3G,程序会进而判断你当前体系里的所有元素是否对于STO-3G有定义,如果没有就报错。对很重的元素STO-3G都是没有定义的。为了避免这个报错,解决办法就是关键词里写上UGBS,这是一个几乎涵盖整个周期表的基组,因此就不会因为当前的基组对元素没有定义而报错了

北京科音自然科学研究中心http://www.keinsci.com)致力于计算化学的发展和传播,长期开办高质量的各种计算化学类培训:初级量子化学培训班中级量子化学培训班高级量子化学培训班量子化学波函数分析与Multiwfn程序培训班分子动力学与GROMACS培训班CP2K第一性原理计算培训班,内容介绍以及往届资料购买请点击相应链接查看。这些培训是计算化学从零快速入门以及进一步全面系统性提升研究水平的高速路!培训各种常见问题见《北京科音办的培训班FAQ》
欢迎加入“北京科音”微信公众号获取北京科音培训的最新消息、避免错过网上有价值的计算化学文章!
欢迎加入人气非常高、专业性特别强的综合性理论与计算化学交流QQ群“思想家公社QQ群”:1号:18616395,2号:466017436,3号:764390338,搜索群号能搜到哪个说明目前哪个能加,合计9000人。北京科音培训班的学员在群中可申请VIP头衔,提问将得到群主Sobereva的最优先解答。
思想家公社的门口Blog:http://sobereva.com(发布大量原创计算化学相关博文)
Multiwfn主页:http://sobereva.com/multiwfn(十分强大的量子化学波函数分析程序)
ResearchGate:https://www.researchgate.net/profile/Tian_Lu
Money and papers are rubbish, get a real life!

本版积分规则 Credits rule

手机版 Mobile version|北京科音自然科学研究中心 Beijing Kein Research Center for Natural Sciences|京公网安备 11010502035419号|计算化学公社 — 北京科音旗下高水平计算化学交流论坛 ( 京ICP备14038949号-1 )|网站地图

GMT+8, 2024-11-26 11:54 , Processed in 0.173352 second(s), 23 queries , Gzip On.

快速回复 返回顶部 返回列表 Return to list