| 查看: 2022 | 回复: 11 | |||
| 当前主题已经存档。 | |||
[交流]
【求助】[求助]用高斯算分子间作用力
|
|||
|
想用高斯计算两个分子间的弱作用力（如氢键等）。现在将两个分子各自优化了（在别人的指导下，选择的B3LYP cc-pVDZ基组），得到能量（Ａ），能量（Ｂ）。然后想把这两个分子合起来优化，希望得到一个能量（Ａ＋Ｂ），用能量（Ａ）+能量（Ｂ）跟能量（Ａ+B)的差值来衡量两者间的作用。 想问的是，怎样优化A+B. （１）搭建分子的时候，就按照自己预想的氢键作用点把Ａ和Ｂ摆在一起了。有什么方法让程序自动寻找可能的作用点么？ （２）写ｇｊｆ输入文件的时候，如果用三维坐标作为原子参数，应该注明为两个分子，还是直接把A+B 作为一个分子优化呢？ 看出来是外行了吧，轻拍  |
回复此楼» 猜你喜欢
能源环境新材料研究团队 招收 新能源、冶金、材料、化工、环境 等工科类研究生
已经有0人回复
-
武汉纺织大学教育部重点实验室调剂
已经有0人回复
-
物理化学论文润色/翻译怎么收费?
已经有215人回复
-
理论与计算化学方向调剂——初试388分（数学144）
已经有0人回复
-
能源环境新材料研究团队 研究生调剂招生
已经有0人回复
-
齐齐哈尔大学材料科学与工程学院有调剂名额
已经有0人回复
-
苏州科技大学 物理学专业 接收调剂生
已经有0人回复
-
硕士研究生调剂
已经有0人回复
-
苏州科技大学 物理学专业 接收调剂生
已经有0人回复
-
苏州科技大学 物理学专业 接收调剂生
已经有0人回复
-
高效微纳米储能材料与产业化推进实验室 苏州 招收硕士研究生（含调剂）
已经有0人回复
» 本主题相关商家推荐: (我也要在这里推广)
bay__gulf
金虫 (著名写手)
刘苏州
- 应助: 9 (幼儿园)
- 贵宾: 4.85
- 金币: 2332.8
- 红花: 1
- 帖子: 1344
- 在线: 271小时
- 虫号: 592012
- 注册: 2008-09-03
- 专业: 理论和计算化学
- 管辖: 分子模拟
2楼2009-08-02 07:29:23
zhou2009
版主 (著名写手)
- QC强帖: 15
- 应助: 337 (大学生)
- 贵宾: 4.053
- 金币: 42654.6
- 红花: 153
- 帖子: 1729
- 在线: 1274.7小时
- 虫号: 676800
- 注册: 2008-12-18
- 性别: GG
- 专业: 理论和计算化学
- 管辖: 量子化学
3楼2009-08-02 08:55:38
lsq_phi
木虫 (著名写手)
- 应助: 0 (幼儿园)
- 金币: 7759.2
- 散金: 10
- 红花: 1
- 帖子: 1079
- 在线: 234.7小时
- 虫号: 808341
- 注册: 2009-07-13
- 专业: 色谱分析
4楼2009-08-02 09:35:08
5楼2009-08-03 00:35:15
neweroica
木虫 (著名写手)
- 应助: 53 (初中生)
- 金币: 5032.5
- 散金: 301
- 红花: 31
- 帖子: 2490
- 在线: 217.1小时
- 虫号: 36775
- 注册: 2004-02-11
- 专业: 理论和计算化学
6楼2009-08-03 06:58:55
neweroica
木虫 (著名写手)
- 应助: 53 (初中生)
- 金币: 5032.5
- 散金: 301
- 红花: 31
- 帖子: 2490
- 在线: 217.1小时
- 虫号: 36775
- 注册: 2004-02-11
- 专业: 理论和计算化学
★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★
lei0736(金币+6,VIP+0):谢谢 8-3 20:13
sprayan(金币+10,VIP+0):多谢指点。待我慢慢消化消化 8-3 09:01
lei0736(金币+6,VIP+0):谢谢 8-3 20:13
sprayan(金币+10,VIP+0):多谢指点。待我慢慢消化消化 8-3 09:01
|
具体方法的使用跟作用力性质有关： 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自身的几何构型进行优化。 不用注明两个分子——除非在计算BSSE单点（甚至优化构型，但这相当耗时）时才需要指明分子1和分子2。 [ Last edited by neweroica on 2009-8-3 at 07:16 ] |
7楼2009-08-03 07:12:11
piaoma
木虫 (小有名气)
- 应助: 6 (幼儿园)
- 金币: 2725.5
- 散金: 157
- 帖子: 162
- 在线: 590.2小时
- 虫号: 201836
- 注册: 2006-02-28
- 性别: GG
- 专业: 理论和计算化学
8楼2009-08-03 20:19:11
neweroica
木虫 (著名写手)
- 应助: 53 (初中生)
- 金币: 5032.5
- 散金: 301
- 红花: 31
- 帖子: 2490
- 在线: 217.1小时
- 虫号: 36775
- 注册: 2004-02-11
- 专业: 理论和计算化学
★ ★
lei0736(金币+2,VIP+0):谢谢 8-3 22:13
lei0736(金币+2,VIP+0):谢谢 8-3 22:13
9楼2009-08-03 21:11:12
piaoma
木虫 (小有名气)
- 应助: 6 (幼儿园)
- 金币: 2725.5
- 散金: 157
- 帖子: 162
- 在线: 590.2小时
- 虫号: 201836
- 注册: 2006-02-28
- 性别: GG
- 专业: 理论和计算化学
10楼2009-08-04 12:49:32
