24小时热门版块排行榜    

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

fcy110918

铁杆木虫 (正式写手)

[求助] Gaussian计算π-π堆积弱相互作用的势能曲线 已有3人参与

各位小木虫计算大神:

小虫本来学生物的,由于实验需要,接触Guassian和DFT不久,很多地方一片迷茫啊。目前想考虑利用Guassian计算分子间的π-π堆积结构的势能曲线,比如说两个蒽分子堆积在一起时势能与距离的曲线,所想要的结果如附件所示,
本人的分子已经opt,也获得了freq,但是不会设置scan中的坐标(不知道是不是这个功能可以得到像下面一样的数据),也不会如何控制两个分子间的距离,最想得到的结果就是和下面的附件图一样(暂时只会用Gaussian)
诚心希望得到大神们的指点

Gaussian计算π-π堆积弱相互作用的势能曲线
w75h2013403_1473375893_866.jpg


Gaussian计算π-π堆积弱相互作用的势能曲线-1
Screen Shot 2016-11-02 at 20.04.29.png


Gaussian计算π-π堆积弱相互作用的势能曲线-2
Screen Shot 2016-11-02 at 20.04.43.png
回复此楼

» 收录本帖的淘帖专辑推荐

自己感觉好的东西

» 猜你喜欢

» 本主题相关商家推荐: (我也要在这里推广)

» 本主题相关价值贴推荐,对您同样有帮助:

已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

fcy110918

铁杆木虫 (正式写手)

送红花一朵
引用回帖:
2楼: Originally posted by 小范范1989 at 2016-11-03 08:38:41
说一下我的理解,不一定对哈,因为之前我了解过一点点,
1:首先,你要找到两个分子,平行排列。(当然也不一定非要平行,向你后面的分子,垂直)。
这个怎么弄呢?拿两个平行的来举例子:
通过高斯view画一个苯 ...

谢谢大侠的回复!由于我要用到的分子比单纯两个苯环要复杂的多,所以手动调节比较麻烦呢,应该是会有可以自动设置的方法!
还是很感激你的帮助哦!!!
4楼2016-11-03 19:21:36
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
查看全部 13 个回答

小范范1989

木虫 (著名写手)

说一下我的理解,不一定对哈,因为之前我了解过一点点,
1:首先,你要找到两个分子,平行排列。(当然也不一定非要平行,向你后面的分子,垂直)。
这个怎么弄呢?拿两个平行的来举例子:
通过高斯view画一个苯环,然后把他的坐标复制一下,再放到文件的后面,这样是不是就是两个重合的苯环了?然后,你在把刚刚复制的那个苯环的Z坐标,都加上1或者2或者3或者4....这不就是两者之间的距离的调节了吗?
如果是垂直,我感觉是不是先通过高斯view画出来,然后调节一个分子的Z坐标,同时加多少,这样就能调节距离了吧?
应该找个过程不需要优化吧,这个我不确定,我感觉是直接做单点计算,当然,我的理解不一定对。
2:有了坐标之后,就要计算能量,
把不同距离的分子,保存不同的文件,比如说0.5.gjf   1.gjf   1.5.gjf......然后做单点计算(我感觉单点,不需要优化,不一定对)。
3:这个地方计算的泛函比较重要,因为必须要有D3校正,如果是b3lyp之类的,结果会是一条一直下降的线,不会有凹的表现。
选泛函,这个在你那个书的169页,第21也ppt、

4:个人见解,不一定对,还望指点

» 本帖已获得的红花(最新10朵)

It doesn't matter how slow you are, as long as you're determined to get there, you'll get there.
2楼2016-11-03 08:38:41
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

ra2ghgzh

木虫 (正式写手)

学痴

【答案】应助回帖

★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ...
感谢参与,应助指数 +1
fcy110918: 金币+50, ★★★很有帮助, 辛苦了! 2016-11-03 19:22:35
fcy110918: 金币+100, ★★★很有帮助, 再次谢谢啦! 2016-11-04 00:39:52
下面是两个平行苯环的PES计算
%mem=600MW
%nprocshared=16
# opt=z-matrix b3lyp/cc-pvdz

Title Card Required

0 1
C                  0.00000000    1.39499067    0.00000000
C                 -1.20809735    0.69749533    0.00000000
C                 -1.20809735   -0.69749533    0.00000000
C                  0.00000000   -1.39499067    0.00000000
C                  1.20809735   -0.69749533    0.00000000
C                  1.20809735    0.69749533    0.00000000
H                  0.00000000    2.49460097    0.00000000
H                 -2.16038781    1.24730049    0.00000000
H                 -2.16038781   -1.24730049    0.00000000
H                  0.00000000   -2.49460097    0.00000000
H                  2.16038781   -1.24730049    0.00000000
H                  2.16038781    1.24730049    0.00000000
C                  0.00000000    1.39499067    R
C                 -1.20809735    0.69749533    R
C                 -1.20809735   -0.69749533    R
C                  0.00000000   -1.39499067    R
C                  1.20809735   -0.69749533    R
C                  1.20809735    0.69749533    R
H                  0.00000000    2.49460097    R
H                 -2.16038781    1.24730049    R
H                 -2.16038781   -1.24730049    R
H                  0.00000000   -2.49460097    R
H                  2.16038781   -1.24730049    R
H                  2.16038781    1.24730049    R
Variables:
R 2.0 S 30 0.1
得到的结果是
Gaussian计算π-π堆积弱相互作用的势能曲线-3
计算结果不吻合的原因可能是计算方法和基组没选择好

» 本帖已获得的红花(最新10朵)

我永远当不了老板
3楼2016-11-03 11:37:37
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

fcy110918

铁杆木虫 (正式写手)

送红花一朵
引用回帖:
3楼: Originally posted by ra2ghgzh at 2016-11-03 11:37:37
下面是两个平行苯环的PES计算
%mem=600MW
%nprocshared=16
# opt=z-matrix b3lyp/cc-pvdz

Title Card Required

0 1
C                  0.00000000    1.39499067    0.00000000
C                 -1 ...

大侠,你的回复很有帮助,非常感谢!
不过我有一些想继续请教的地方啊,
1.你的坐标中如何能够归零啊,比如Z坐标很多都是零,我都没有试过哈,
2. 那些R坐标都是自己手动改的吗?还是在gview里面可以直接设置?
3.最后语句Variables: R 2.0 S 30 0.1 和网上的一些例子,比如 B 1 8 S 21 -0.100000,两者有什么异同点呢?R.B.S什么的都是啥意思啊!

问题有点多哈,辛苦啦!!
5楼2016-11-03 19:29:55
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 26药学专硕105500求调剂 +7 喽哈加油 2026-04-13 8/400 2026-04-19 20:21 by Equinoxhua
[考研] 297,工科调剂? +11 河南农业大学-能 2026-04-14 11/550 2026-04-19 20:07 by Equinoxhua
[考研] 291求调剂 +11 关忆北. 2026-04-14 11/550 2026-04-19 17:16 by 中豫男
[考研] 求调剂 +10 小聂爱学习 2026-04-16 12/600 2026-04-19 16:51 by 中豫男
[考研] 307中医考研调剂 +9 于以采蘩 2026-04-14 9/450 2026-04-19 08:41 by 烟雨流涯
[考研] 生物学调剂 +10 纸扇zhishan 2026-04-13 10/500 2026-04-18 21:19 by zqndavala
[考研] 320求调剂 +5 深郊akm 2026-04-17 5/250 2026-04-18 19:52 by 王珺璞
[考研] 接受任何调剂 +6 也就是栗子 2026-04-17 7/350 2026-04-18 17:20 by 涵竹刘
[考研] 260求调剂 +4 Zyt1314520.. 2026-04-17 5/250 2026-04-18 08:28 by babysonlkd
[考研] 一志愿华中农业071010,320求调剂 +17 困困困困坤坤 2026-04-14 19/950 2026-04-17 20:08 by 关一盏灯cd
[有机交流] 二苯甲酮酸类衍生物 50+3 小白爱主人 2026-04-17 6/300 2026-04-17 18:47 by kf2781974
[考研] 一志愿中科大材料与化工,353分还有调剂学校吗 +10 否极泰来2026 2026-04-15 12/600 2026-04-17 17:54 by mapenggao
[考研] 295分求调剂 +5 ?要上岸? 2026-04-17 5/250 2026-04-17 16:51 by fenglj492
[考研] 一志愿沪9,生物学326求调剂 +9 刘墨墨 2026-04-15 9/450 2026-04-16 17:14 by 崔崔崔cccc
[考研] 297,工科调剂?河南农业大学本科 +14 河南农业大学-能 2026-04-14 14/700 2026-04-16 14:41 by dingyanbo1
[基金申请] RY:中国产出的科学垃圾论文,绝对数量和比例都世界第一 +7 zju2000 2026-04-14 18/900 2026-04-16 11:36 by 欢乐颂叶蓁
[考研] 求调剂学校 +14 不会吃肉 2026-04-13 16/800 2026-04-15 21:59 by noqvsozv
[考研] 考研调剂 +13 长弓傲 2026-04-13 14/700 2026-04-14 14:44 by zs92450
[考研] 085408光电信息工程专硕355一志愿长春光机所调剂 +6 王ymaa 2026-04-13 13/650 2026-04-14 11:33 by 王ymaa
[考研] 085600材料与化工329分求调剂 +24 叶zilin 2026-04-13 25/1250 2026-04-14 09:20 by 试管破裂
信息提示
请填处理意见