logo

您所在位置网站首页 > 海量文档  > 生活休闲 > 科普知识

高斯常见问题.doc 19页

本文档一共被下载: ,您可全文免费在线阅读后下载本文档。

  • 支付并下载
  • 收藏该文档
  • 百度一下本文档
  • 修改文档简介
全屏预览

下载提示

1.本站不保证该用户上传的文档完整性,不预览、不比对内容而直接下载产生的反悔问题本站不予受理。
2.该文档所得收入(下载+内容+预览三)归上传者、原创者。
3.登录后可充值,立即自动返金币,充值渠道很便利
特别说明: 下载前务必先预览,自己验证一下是不是你要下载的文档。
  • 内容提供方 tianma2015(上传创作收益人)
  • 发布时间:2015-09-13
  • 需要金币160(10金币=人民币1元)
  • 浏览人气
  • 下载次数
  • 收藏次数
  • 文件大小:317.5 KB
下载过该文档的会员
高斯常见问题,高斯问题,杜高斯贝银行出金问题,高斯奥特曼,高斯枪,高斯分布,高斯奥特曼国语版全集,高斯奥特曼剧场版,高斯模糊,高斯狙击步枪
你可能关注的文档:
一、混合基组 #p b3lyp/GenECP opt? ? Gen 自己在输入文件末尾定义基组,ECP 使用赝势 空行 g03 opt 空行 0 1 ...??... ...??... 空行 C H O 0? ? 碳 氢 氧 数字0? ? GEN 6-31G**? ?基组名 ****? ?? ???结束符号 Ru P Cl 0? ?Ru P Cl 数字0? ???GEN Lanl2dz? ???基组名 ****? ?? ???结束符号 空行 Ru P Cl 0? ? Ru P Cl 数字0? ? ECP Lanl2dz? ?? ?基组 ===================== 建议你查看 GEN 关键字。上述,如不加ECP亦可。是我感觉你要是maxcycle=200还是999死的话,就减小一下步长试试,可能会有用的!maxstep=1该先看下收敛指标是不是很接近了,如果很接近就继续做(读chk)或取出最后一个坐标继续算,如果相差很远,可以加非线形实验下,再一个就是加大循环,缩短步长。L202是重新定位坐标,计算对称性,检查变量的,它出错应该就是输入结构有问题。我算过的用HF方法虚频超过1000,而用DFT只要100-300.所以说方法很重要。但是,只要你找到的TS虚频与反应的方向一致并且IRC已经验证,这个TS就是可信的。虚频也不能太小,比如小于100.过渡态与你的产物相差甚远??说明在这条反应路径上还有过渡态你没有找到寻找过渡态不是一件容易的事(对于我和大多数刚涉及量化的人来说),因此我希望通过写这个经验小结能对大家有些帮助。 ? ?1.首先遇到的问题是,用哪种方法来寻找过渡态? ? ?GAUSSIAN提供的方法是QSTN和TSN方法。两种方法各有优点和缺点。QSTN方法特别QST3方法要求输入反应物,过渡态的猜测结构,产物这三者的结构。特别麻烦。但很管用,一般不会出现不收敛的情况。对于TSN(对应关键词为OPT=TS)方法,只要求输入过渡态的初始结构,但这个初始结构非常的关键,如果结构不好,则很容易出现不收敛的情况。所以我建议,如果是刚开始做过度态的话,用QSTN方法是好的选择,等有了“感觉”之后,再用TSN方法。 ? ?2. 怎么解决经常出现的错误? ? ?在找过度态的时候,经常碰到的一些问题就是(1)不收敛,(2)有一个错误的本征值(错误信息为:there is a wrong sign eigenvalue in hessian matrix.....),(3)和LINK9999错误导致退出。对于不收敛的情况,可以分为两类,比如提示信息里的CONVERGENCE FAILER 提醒收敛到了10(-5),而此时你设定的SCF循环次数也仅仅是64步,那么完全有希望通过加大SCF循环次数来达到收敛的目的。倘若只收敛到10(-3)或10(-2),此时加大循环次数可能就没用了。结果还是CONVERGE FAILER。此时可采用SCF=QC,来达到强制收敛的目的 。因为SCF=QC(LINK508)的计算量比默认的L502要大,所以不到万不得以就不用它了。出现第二个错误可以直接用关键词OPT=NOEIGEN 来实现。LINK9999出错是因为已经走完了默认的步数,但还未完成。系统会自动跳出。出现这种情况大多数就是因为优化步数和SCF步数超过了默认值。可用OPT(MAXCYCLE=100)和SCF(MAXCYCLE=300)来改错。 ? ?3.怎么样控制过渡态的优化,使得过渡态不至于收敛到其他的分子结构中去? ? ?我用GAUSS VIEW 可以解决这个问题 ,当刚开始运行GAUSSIAN时,你用GVIEW去打开输出文件时,你可以看到你的过渡态的初始输入结构,当一个循环过后(从上一个LINK502到下一个LINK502),你再打开输出文件,你就可以清晰地看到优化一步后分子的构型,这样就可以随时监控过度态分子的结构,倘若已经有收敛到其他分子构型的趋势时,你就可以把它给KILL了,而不至于需要等全部工作结束后,打开输出文件才知道已经不是想要的过渡态了。 如果收敛到其它的构型上去,可以考虑缩小OPT的步长.iop(1/8=2或3)即可。 ? ?4. 还需要加其它的关键词吗? ? ?建议在OPT中加入CALCFC。这样可以加大找到过渡态的几率。本人深有体会! ? ?5. 如何控制寻找过渡态的步长? ? ?这个可以用IOP来实现。具体相应的IOP是iop(1/8)=2,此语句说明在找过度态的时候,以2A(A为基本单位长度)为单位来寻找过渡态。我用高斯03做了两个稳定点A和B的过渡态TS,发现频率文件的振动方向还不足以证明TS就是A和B的过渡态,所以,我又接着做了IRC,结果发现IRC正方向走出的结构与A结构差别很小,而反方向走出的结构与

发表评论

请自觉遵守互联网相关的政策法规,严禁发布色情、暴力、反动的言论。
用户名: 验证码: 点击我更换图片

“原创力文档”前称为“文档投稿赚钱网”,本站为“文档C2C交易模式”,即用户上传的文档直接卖给(下载)用户,本站只是中间服务平台,本站所有文档下载所得的收益归上传人(含作者)所有【成交的100%(原创)】。原创力文档是网络服务平台方,若您的权利被侵害,侵权客服QQ:3005833200 电话:19940600175 欢迎举报,上传者QQ群:784321556