24小时热门版块排行榜    

查看: 565  |  回复: 4
当前主题已经存档。

pleasant1789

木虫 (初入文坛)

[交流] 【求助】Gaussian09溶液中优化的问题 已有4人参与

Gaussian09一个很大的改进说是溶液中计算方法的改进

尝试了一下发现确实是溶液中优化速度快多了,而且不像03中很容易出错
但是同样的文件在03和09中溶液相中优化得到的结果不同,相差还挺大的

想问的问题是:
1。在09中溶液项中优化的头文件和尾文件的撰写是否跟03一样?
比如下面这个
头文件
#p b3lyp/6-31+g* opt freq=noraman scrf=(pcm,read,solvent=water)
尾文件
radii=uaks
是否正确?当然算是能算的……就是不知道对不对

2. 优化得到的结构如何处理能量问题?
Zero-point correction=                           0.128180 (Hartree/Particle)
Thermal correction to Energy=                    0.138888
Thermal correction to Enthalpy=                  0.139832
Thermal correction to Gibbs Free Energy=         0.091443         <-------  TCG
Sum of electronic and zero-point Energies=           -646.989573
Sum of electronic and thermal Energies=              -646.978865
Sum of electronic and thermal Enthalpies=            -646.977921
Sum of electronic and thermal Free Energies=         -647.026311  <------自由能

如果不考虑更高水平方法的单点能计算, 最后一项就是我们要的自由能项了。
但是如果想要得到更高方法级别下的能量,是否就读出这个结构,在溶液中再进行一次单点能计算?
比如
头文件
#p b3lyp/6-311++g** scf=tight scrf=(pcm,read,solvent=water)
尾文件
radii=uaks
TCG依然取优化文件中的数值,加上新算的单点能就是自由能。

溶剂化校正项是包含在TCG项中还是电子能项(就是SCF Done)中呢?

如果用溶液中优化的结构,在气相中进行单点能计算的话是什么情况?还有再进行一次溶剂化单点能计算吗?(这种情况当然就是假设溶剂化校正项是不包含在TCG项中)
回复此楼
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

wh-32290

铁虫 (小有名气)

我也想知道!
2楼2010-04-03 22:12:04
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

pwzhou

铁杆木虫 (正式写手)


小木虫(金币+0.5):给个红包,谢谢回帖交流
由于算法改变,所以03和09的结果不同是正常的,写法和03基本没啥差别。
3楼2010-04-04 10:30:40
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

gyli

木虫 (小有名气)

~~~~


小木虫(金币+0.5):给个红包,谢谢回帖交流
引用回帖:
Originally posted by pleasant1789 at 2010-03-02 22:47:39:
Gaussian09一个很大的改进说是溶液中计算方法的改进

尝试了一下发现确实是溶液中优化速度快多了,而且不像03中很容易出错
但是同样的文件在03和09中溶液相中优化得到的结果不同,相差还挺大的

用03的输入文件在09上提交过
运行不成问题
但未比较过二者差别
Gravitation cannot be held responsible for people falling in love.
4楼2010-04-04 14:43:44
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

fineren

木虫 (正式写手)

搬个板凳来学习
5楼2010-04-04 15:40:13
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 pleasant1789 的主题更新
普通表情 高级回复 (可上传附件)
信息提示
请填处理意见