计算化学公社

标题: 请教几个GAMESS-US进行能量分解分析时遇到的问题 [打印本页]

作者
Author: wangpeng 时间: 2017-6-23 00:24
标题: 请教几个GAMESS-US进行能量分解分析时遇到的问题
大家好！我最近在试着用GAMESS-US进行能量分解分析来研究一些成键的细节，虽然看了一些GAMESS官方的使用手册和网上一些相关的帖子，但是在实际运算中还是经常出错，有下面几个初级问题想向大家请教一下，还请大家不吝赐教，谢谢！
1. 我想用LMO-EDA来分析过渡金属炔化物中TM-C键的细节（例子请见后面的输入文件），不知道是否合适呢？
2. 我看到网上有些人建议用MP2进行LMO-EDA分析，因为我所用的构型是用DFT优化的，不知道继续用DFT进行EDA计算的话会不会和用MP2计算EDA的结果区别比较大？
3. 我在进行构型优化时给过渡金属用的Lanl2TZ基组，给C，H，O，N等主族元素用的是TZVP基组，请问用同样的基组计算EDA会不会有问题呢？另外，写输入文件时基组的书写方式还有些不明白的地方，我在EMSL上可以找到TZVP（DFT orbital）基组，不知道这个是不是就可以直接在GAMESS代替TZVP基组呢？我看有些GAMESS输入文件的例子中的TZVP基组和EMSL上的并不一样，是不是还要在使用时进行一下修改呢？
4. 我参考着一些GAMESS输入文件的范例，自己写了一些输入文件，但是很多计算时会出现错误，比如下面这个，而且好像在输出文件中找不到明显的错误提示，所以还请大家指点一下
(, 下载次数 Times of downloads: 74) (, 下载次数 Times of downloads: 21)

作者
Author: sobereva 时间: 2017-6-23 02:53

1 当前的片段划分显得莫名其妙

2 没必要用MP2，原理上并不会更好，本身MP2算配合物也不算理想

3 意义不明，没有任何问题。
甭用TZVP（DFT orbital），弄清楚到底要用def-TZVP还是def2-TZVP。

4 有提示-SOSCF- AND -DIIS- ARE MUTUALLY EXCLUSIVE CONVERGERS
SOSCF=.t. DIRSCF=.T.别同时用

作者
Author: wangpeng 时间: 2017-6-23 15:53

sobereva 发表于 2017-6-23 02:53
1 当前的片段划分显得莫名其妙

2 没必要用MP2，原理上并不会更好，本身MP2算配合物也不算理想

谢谢sob老师指点

，对于片段划分我还是有点不明白，比如在我这个分子Au-CCH中，我想看一下Au和C之间的能量，是不是直接划分为1和3两个片段就行了，还是有其他的规则呢？希望您能再次抽空指点一下，谢谢！

作者
Author: sobereva 时间: 2017-6-23 16:03

wangpeng 发表于 2017-6-23 15:53
谢谢sob老师指点，对于片段划分我还是有点不明白，比如在我这个分子Au-CCH中，我想看一下Au和C之间的 ...

你得按照顺序来划分
Au当前体系是第四个原子，应当划分成3,1

作者
Author: wangpeng 时间: 2017-7-3 16:09

sobereva 发表于 2017-6-23 16:03
你得按照顺序来划分
Au当前体系是第四个原子，应当划分成3,1

多谢sob老师指点哈，这几天我用不同的方法算了一下EDA，结果看起来挺合理的，不过还有个小问题想向sob老师请教一下，按照ADF中的ETS-NOCV方法分析，总结合能会被分解成了静电相互作用，Pauling排斥作用，轨道相互作用和色散作用四项，请问是不是可以把LMO-EDA结果中的EX项加上REP项等同于ETS-NOCV结果里的Pauling排斥作用，另外LMO-EDA中的POL项是否可以等同于ETS-NOCV中的轨道相互作用一项呢？

作者
Author: sobereva 时间: 2017-7-4 01:58

wangpeng 发表于 2017-7-3 16:09
多谢sob老师指点哈，这几天我用不同的方法算了一下EDA，结果看起来挺合理的，不过还有个小问题想向sob老 ...

单从物理意义上看有这样的对应性，但未必实际能定性对上，我没检验过

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