24小时热门版块排行榜    

北京石油化工学院2026年研究生招生接收调剂公告
查看: 2025  |  回复: 11
当前主题已经存档。
当前只显示满足指定条件的回帖,点击这里查看本话题的所有回帖

sprayan

铜虫 (初入文坛)

[交流] 【求助】[求助]用高斯算分子间作用力

想用高斯计算两个分子间的弱作用力(如氢键等)。现在将两个分子各自优化了(在别人的指导下,选择的B3LYP cc-pVDZ基组),得到能量(A),能量(B)。然后想把这两个分子合起来优化,希望得到一个能量(A+B),用能量(A)+能量(B)跟能量(A+B)的差值来衡量两者间的作用。

想问的是,怎样优化A+B. 
(1)搭建分子的时候,就按照自己预想的氢键作用点把A和B摆在一起了。有什么方法让程序自动寻找可能的作用点么?
(2)写gjf输入文件的时候,如果用三维坐标作为原子参数,应该注明为两个分子,还是直接把A+B 作为一个分子优化呢?

看出来是外行了吧,轻拍
回复此楼
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

neweroica

木虫 (著名写手)

★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★
lei0736(金币+6,VIP+0):谢谢 8-3 20:13
sprayan(金币+10,VIP+0):多谢指点。待我慢慢消化消化 8-3 09:01
引用回帖:
Originally posted by sprayan at 2009-8-2 05:06:
想用高斯计算两个分子间的弱作用力(如氢键等)。现在将两个分子各自优化了(在别人的指导下,选择的B3LYP cc-pVDZ基组),得到能量(A),能量(B)。然后想把这两个分子合起来优化,希望得到一个能量(A+B),用能量(A)+能量(B)跟能量(A+B)的差值来衡量两者间的作用。

具体方法的使用跟作用力性质有关:

1. 强氢键(X--H ... Y, X,Y=O, F,Cl等),用B3LYP还是说得过去的。为保险起见,可以尝试其他泛函,如PBE,PW91等等,然后进行比较。

2. 弱氢键(X=C, Y=S等),用密度泛函就差得较远,不过基本几何构型还凑合,但作用能误差可能会很大。

3. 范德华作用,密度泛函基本不行。但你可以尝试一些新泛函M05-2X等,或者采用DFT+D(即在DFt能上加经验色散作用校正项)。不过,这些方法都很取决于具体体系,局限性较大,使用时得小心。

4. 如果你的体系不是很大,而且机器允许,可以用更可靠的post-HF方法,如MP2,MP3,甚至MP4,或者CCSD,CCSD(T)等。不过使用这些方法时,基组不能取得太小(最小是6-31+G*),否则结果还是不合理——因为这些方法中的相关能校正项目是基于HF的波函数的。

5. 如果基组不是很大的话,计算作用能时需要考虑BSSE校正。

关于分子间作用力的量化计算,PCCP有一期有个专刊系列,你可以查一下。
引用回帖:
想问的是,怎样优化A+B. 
(1)搭建分子的时候,就按照自己预想的氢键作用点把A和B摆在一起了。有什么方法让程序自动寻找可能的作用点么?

我没听说过这样的程序。不过,如果仅仅找“可能的作用点”,一般目测即可(用化学常识判断)。要稍精确些,可以多给几种可能的初始猜测构型,然后固定A和B自身的几何构型进行优化。
引用回帖:
(2)写gjf输入文件的时候,如果用三维坐标作为原子参数,应该注明为两个分子,还是直接把A+B 作为一个分子优化呢?

不用注明两个分子——除非在计算BSSE单点(甚至优化构型,但这相当耗时)时才需要指明分子1和分子2。

[ Last edited by neweroica on 2009-8-3 at 07:16 ]
Simpler, stronger, more tolerant
7楼2009-08-03 07:12:11
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
查看全部 12 个回答

bay__gulf

金虫 (著名写手)

刘苏州

★ ★ ★ ★ ★ ★
lei0736(金币+1,VIP+0):谢谢 8-2 13:27
sprayan(金币+5,VIP+0):谢谢回复。那范德华力该用什么算呢? 8-2 10:32
量化方法不能计算vdW 力
2楼2009-08-02 07:29:23
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

zhou2009

版主 (著名写手)

★ ★ ★ ★ ★ ★ ★
lei0736(金币+2,VIP+0):谢谢 8-2 13:27
sprayan(金币+5,VIP+0):谢谢回复 8-2 10:33
如果算氢键,优化A+B时:
(1)搭建分子的时候,就按照自己预想的氢键作用点把A和B摆在一起优化。
(2)写gjf输入文件的时候,如果用三维坐标作为原子参数,直接把A+B 作为一个分子优化。
3楼2009-08-02 08:55:38
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

lsq_phi

木虫 (著名写手)


sprayan(金币+1,VIP+0):谢谢参与 8-2 10:33
不知道,但关注中...:dnd:
4楼2009-08-02 09:35:08
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
普通表情 高级回复 (可上传附件)
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 323求调剂(计算机视觉和大模型项目经历) +3 chaoxiicy 2026-03-31 3/150 2026-04-05 10:33 by zhq0425
[考研] 材料调剂 +11 一样YWY 2026-04-02 13/650 2026-04-04 23:10 by 无际的草原
[考研] 286求调剂 +3 草木不言 2026-04-04 3/150 2026-04-04 22:40 by lbsjt
[考研] 323求调剂 +8 李佳乐1 2026-04-04 8/400 2026-04-04 22:26 by hemengdong
[考研] 材料科学与工程调剂 +19 深V宿舍吧 2026-03-30 20/1000 2026-04-04 22:13 by hemengdong
[考研] 282电子信息0854专硕调剂 +4 202451007219 2026-04-02 6/300 2026-04-04 21:55 by laoshidan
[考研] 278求调剂 +3 依旧! 2026-04-02 4/200 2026-04-04 20:27 by 蓝云思雨
[考研] 285求调剂 +4 AZMK 2026-04-04 5/250 2026-04-04 16:45 by cql1109
[考研] 一志愿武理材料工程302调剂环化或化工 +19 Doleres 2026-03-31 20/1000 2026-04-04 16:44 by 啊俊!
[考研] 一志愿085404,总分291,四级已过,求调剂 +5 阿俊阿俊阿俊 2026-04-04 7/350 2026-04-04 13:23 by 莲菜就是藕吧
[考研] 265求调剂 +20 梁梁校校 2026-04-01 21/1050 2026-04-04 00:38 by userper
[考研] 286求调剂 +8 lim0922 2026-04-02 8/400 2026-04-03 20:19 by rzh123456
[考研] 285求调剂 +5 AZMK 2026-04-03 8/400 2026-04-03 18:17 by AZMK
[考研] 材料专硕调剂 +18 椰椰。 2026-03-29 18/900 2026-04-03 16:45 by 玲玲0606
[考研] 22408 266求调剂 +3 masss11222 2026-04-02 3/150 2026-04-02 18:11 by 笔落锦州
[考研] 270调剂 +7 maxjxbsk 2026-04-02 7/350 2026-04-02 09:50 by yulian1987
[考研] 286求调剂 +5 Sa67890. 2026-04-01 7/350 2026-04-01 19:50 by 6781022
[考研] 290求调剂 +5 dfffsar 2026-03-29 5/250 2026-04-01 19:45 by 6781022
[考研] 材料与化工(0856)304求B区调剂 +8 邱gl 2026-03-30 16/800 2026-04-01 17:58 by 邱gl
[考研] 288资源与环境专硕求调剂,不限专业,有学上就行 +25 lllllos 2026-03-30 26/1300 2026-04-01 09:52 by 一只好果子?
信息提示
请填处理意见