24小时热门版块排行榜    

查看: 2273  |  回复: 7
当前只显示满足指定条件的回帖,点击这里查看本话题的所有回帖

tahaomei

铜虫 (小有名气)

[求助] 柔性扫描两个H2分子之间的势能曲线,该怎么写gauss配置文件? 已有3人参与

考虑bsse矫正,请问该怎么写配置文件呢?gaussian教程看了不少,还是不清楚。谢谢啦
回复此楼
从是西方,过十万亿佛土,有世界名曰极乐。其土有佛,号阿弥陀,今现在说法。设我得佛,十方众生,至心信乐,欲生我国,乃至十念,若不生者,不取正觉。唯除五逆
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

tahaomei

铜虫 (小有名气)

我现在不算H2,算C2 分子。这么写对不?

%chk=c2.chk
%mem=8GB
# opt=modredundant mp2/aug-cc-pvdz counterpoise=2

c2 c2 interaction

0 1
x 0.0 0.0 0.0
c 0.628765 0.0 0.0 1
c -0.628765 0.0 0.0 1
x 0.0 0.0 rcc
c 0.0 0.628765 rcc 2
c 0.0 -0.628765 rcc 2

rcc 2.0 S 10 0.3
从是西方,过十万亿佛土,有世界名曰极乐。其土有佛,号阿弥陀,今现在说法。设我得佛,十方众生,至心信乐,欲生我国,乃至十念,若不生者,不取正觉。唯除五逆
7楼2014-07-10 11:42:37
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
查看全部 8 个回答

oinkmasta

金虫 (小有名气)

【答案】应助回帖

感谢参与,应助指数 +1
BSSE矫正在命令行输入counterpoise=n就可以,n指的是输入的分子(或聚合物的单体等)的数量。矫正值在输出文件(以文件形式打开,或者gview里面stream出来)可以读出。
柔性扫描参照 http://www.gaussian.com/g_tech/g_ur/k_opt.htm,在页面上搜索scan就能找到你要的说明了。主要就是要用到opt=addredundant,以及在分子结构那一块后面要 空行,加一行命令说明addredundant要做什么,空行 。
有问题欢迎继续提出~
16届本科在北美
2楼2014-07-09 23:03:30
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

tahaomei

铜虫 (小有名气)

引用回帖:
2楼: Originally posted by oinkmasta at 2014-07-09 23:03:30
BSSE矫正在命令行输入counterpoise=n就可以,n指的是输入的分子(或聚合物的单体等)的数量。矫正值在输出文件(以文件形式打开,或者gview里面stream出来)可以读出。
柔性扫描参照 http://www.gaussian.com/g_t ...

问题是,两个H2分子之间不同的定位,是否会影响扫描结果?

比如说,两个H2分子最开始的时候,可以是相互垂直的,也可以使相互平行的。

还有,是否要以H2分子的中心为X原子来写配置文件?
从是西方,过十万亿佛土,有世界名曰极乐。其土有佛,号阿弥陀,今现在说法。设我得佛,十方众生,至心信乐,欲生我国,乃至十念,若不生者,不取正觉。唯除五逆
3楼2014-07-10 07:32:50
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

oinkmasta

金虫 (小有名气)

【答案】应助回帖

分子不太能被看成单个原子吧..扫描方面,一般只能扫描 坐标/键长/键角/二面角,因为要控制那一个自由度,优化其他的自由度。
你扫描想要得到的是什么信息?  一般来说我会先优化这两个H2分子,得到一个geometry(翻成“几何”太别扭了),然后比如说我想知道氢原子A和氢原子B之间间距对能量的影响,我就对它们的间距(即键长,不用担心有没有键)进行扫描。   每个扫描得到的都是在控制 间距=xxx 的情况下的优化结果。也就是说,你不需要考虑它垂直或者平行有什么影响,Gaussian都是直接优化把能量降下去的,很有可能两种情况优化出来完全是同一个geometry。
16届本科在北美
4楼2014-07-10 08:29:09
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
信息提示
请填处理意见