| 查看: 1146 | 回复: 16 | ||
| 当前主题已经存档。 | ||
| 当前只显示满足指定条件的回帖，点击这里查看本话题的所有回帖 | ||
peptide2882木虫 (正式写手)
|
[交流]
【讨论】菜鸟学GAMESS的心得
|
|
|
我是做合成的，计算一窍不通，由于实际的需要，要接触一下计算，选了WINGAMESS软件，想记录一下自己的学习过程，既有利于得到高手的批评指正，也有利于给更新的菜鸟一些参考。 我计算的东西是一个二肽，目的是算出其肽键由顺式异构化变成反式的能垒，以及顺反之间的能量差。 首先明确要算的内容是顺式的能量、反式的能量，以及过渡态的能量。 通过近几天努力的看以及试着运行例子文件，我认为我应该做的是：先用CHEM3D画出顺式的结构和反式的结构，借助于CHEM3D，可以将这两个结构能量优化一下，然后保存为GAMESS的输入文件； 第二步，编辑这个输入文件，用尽可能的方法基组及计算命令完善这个输入文件； 第三步，运行GAMESS，优化这两个结构，分别计算出顺式文件的局域最小能和反式文件的全局最小能，这样就得到了顺式和反式的能量。 第四步，借助于顺式和反式的优化文件，选用合适的寻找过渡态的命令来计算出过渡态的能量（过渡态要经过频率计算验证）。 有了这三个能量值，我想要的结果就出来了。这只是我现在脑子里的蓝图，具体还没操作，因为还在看说明，熟悉应该选什么样的关键词和命令，所以很想请高手为我点评，我的蓝图是否合理。谢谢阿！ [ Last edited by woshilsh on 2009-6-13 at 21:56 ] |
回复此楼
peptide2882
木虫 (正式写手)
- 应助: 1 (幼儿园)
- 金币: 3030.3
- 散金: 74
- 红花: 2
- 帖子: 520
- 在线: 96.4小时
- 虫号: 709637
- 注册: 2009-02-26
- 性别: GG
- 专业: 有机分子功能材料化学
我宁愿是后者
3楼2009-03-05 17:10:46
loovfnd
至尊木虫 (著名写手)
- 应助: 35 (小学生)
- 贵宾: 1.111
- 金币: 14361.4
- 散金: 236
- 红花: 4
- 帖子: 2790
- 在线: 187小时
- 虫号: 219674
- 注册: 2006-03-14
- 专业: 物理有机化学
2楼2009-03-05 16:51:07
peptide2882
木虫 (正式写手)
- 应助: 1 (幼儿园)
- 金币: 3030.3
- 散金: 74
- 红花: 2
- 帖子: 520
- 在线: 96.4小时
- 虫号: 709637
- 注册: 2009-02-26
- 性别: GG
- 专业: 有机分子功能材料化学
★ ★ ★
yuhuobuku(金币+3,VIP+0):欢迎参加讨论 3-7 08:08
yuhuobuku(金币+3,VIP+0):欢迎参加讨论 3-7 08:08
|
今天用上面的思路壮着胆做了一下，手动画的顺式和反式经初步的EM后保存为输入文件，寻找最小能量，HF/6-31G，SCFTYP=RHF， RUNTYP=OPTIMIZE，GUESS=HUCKEL...共38个原子的体系，顺、反式各跑了一个多小时，表观上看顺利的跑完，怀着激动的心情打开out文件看了看，没有看到希望看到的顺式能量比反式能量高的结果，二者基本都是在-1270左右（单位是kcal/mol吗？），而且最致命的是给出了一句话，“寻找全局最低点失败，未收敛”。 其实自己并不知道这样设置参数是否合理，只是看文献上类似的计算用的这中设置，就照葫芦画瓢了，包括对称性的设置我也不太懂，不知道二肽属于什么对称，默认了个C1 （有知道的朋友给咱提示一下呗），今天最大的收获就是第一次算出个东西来，并且有了个计算时间的概念，还得继续努力修改输入文件啊  如果朋友们看到有不妥的地方，非常希望能不吝赐教  |
4楼2009-03-07 00:37:51
peptide2882
木虫 (正式写手)
- 应助: 1 (幼儿园)
- 金币: 3030.3
- 散金: 74
- 红花: 2
- 帖子: 520
- 在线: 96.4小时
- 虫号: 709637
- 注册: 2009-02-26
- 性别: GG
- 专业: 有机分子功能材料化学
5楼2009-03-08 23:19:19
50