| 查看: 1992 | 回复: 7 | ||||
幸福男人金虫 (小有名气)
|
[求助]
g09中优化构型的能量如何表示
|
|
使用G09程序包，先用小基组b3lyp/6-31g(d,p)对构型进行优化和频率分析，得到的能量如下： Zero-point vibrational energy 349660.1 (Joules/Mol) 83.57077 (Kcal/Mol) Zero-point correction= 0.133179 (Hartree/Particle) Thermal correction to Energy= 0.147309 Thermal correction to Enthalpy= 0.148677 Thermal correction to Gibbs Free Energy= 0.081790 (用A表示) Sum of electronic and zero-point Energies= -346.650693 Sum of electronic and thermal Energies= -346.636563 Sum of electronic and thermal Enthalpies= -346.635194 Sum of electronic and thermal Free Energies= -346.702081 (用B表示) 再用大基组b3lyp/6-311+g(2df,2p)计算单点能： SCF Done: E(RB3LYP) = -346.895863675 （用C表示） 最后考虑到溶剂效应： SCF Done: E(RB3LYP) = -346.904863644 （用D表示） A.U. after 12 cycles Convg = 0.5522D-08 -V/T = 2.0047 KE= 3.452718482840D+02 PE=-1.488925327817D+03 EE= 4.562252983460D+02 SMD-CDS (non-electrostatic) energy (kcal/mol) = -0.66 （用E表示） (included in total energy above) 那么优化构型的最终能量该如何表示呢？（由于反应温度为135摄氏度，想用自由能表示）我初学高斯，还请各位赐教。 |
回复此楼» 收录本帖的淘帖专辑推荐
 高斯学习中的基础问题 |
» 猜你喜欢
招贤纳博（已结束）
已经有28人回复
-
表征测试机构后台私聊我
已经有0人回复
-
物理化学论文润色/翻译怎么收费?
已经有112人回复
-
Coordination Chemistry Reviews
已经有0人回复
» 本主题相关价值贴推荐，对您同样有帮助:
Total energy 和 Gibbs free energy的关系
已经有11人回复- 基组or结构问题，B3lyp优化成功，B3PW91/WB97XD优化全出错
已经有16人回复
- 结构优化及频率计算问题
已经有4人回复
- 请教用gaussian优化构型后得到的是log文件，如何把结构导出？
已经有7人回复
- Material studio Dmol3 优化结构和 Gaussian09 NMR 计算
已经有9人回复
- 求助：优化弱相互作用的团簇构型的S1态应该使用什么方法？
已经有12人回复
- gaussian不同半经验优化得到的能量相差较大
已经有4人回复
- 关于H2和CO生成甲醛的反应络合物优化构型的计算
已经有4人回复
- g09优化时出现L1错误
已经有21人回复
- 在gaussian计算中，先用小基组计算，在改成大基组计算出现的问题
已经有23人回复
- G03中的(Polarized solute)-Solvent 的值在G09中怎么找到
已经有10人回复
- 请教高斯09算溶剂效应及能量校正问题，多送金币+小红花
已经有19人回复
- G09 固定坐标优化 出错
已经有6人回复
- 溶剂化构型优化问题
已经有5人回复
- 关于TD-DFT优化构型并计算频率的问题
已经有12人回复
- td做激发态优化出错
已经有19人回复
- 计算时候总是不能收敛如何解决？
已经有5人回复
- G09 TDDFT优化第一激发三重态的几何 多重度如何写
已经有11人回复
- 高斯对多分子体系的构型优化
已经有4人回复
- 关于结构优化中的对称性问题
已经有9人回复
- 【求助】关于激发态优化及发射光谱的计算
已经有36人回复
- 【讨论】结构用opt优化不了，而用opt=z-matrix却能优化
已经有7人回复
- 【求助】如何计算有机小分子的三重态能级（T1-S0）？
已经有12人回复
幸福男人
金虫 (小有名气)
- 应助: 2 (幼儿园)
- 金币: 1202.7
- 帖子: 172
- 在线: 119.3小时
- 虫号: 1418413
- 注册: 2011-09-26
- 性别: GG
- 专业: 理论和计算化学
2楼2012-04-30 22:46:33
幸福男人
金虫 (小有名气)
- 应助: 2 (幼儿园)
- 金币: 1202.7
- 帖子: 172
- 在线: 119.3小时
- 虫号: 1418413
- 注册: 2011-09-26
- 性别: GG
- 专业: 理论和计算化学
3楼2012-05-01 06:49:06
幸福男人
金虫 (小有名气)
- 应助: 2 (幼儿园)
- 金币: 1202.7
- 帖子: 172
- 在线: 119.3小时
- 虫号: 1418413
- 注册: 2011-09-26
- 性别: GG
- 专业: 理论和计算化学
|
在小木虫上我看到：溶液下溶质的总自由能= Esoln（溶液下溶质的总能量）+Gnes（溶解增加的自由能）+deltaGcorr_gas（气相下溶质的热自由能校正）。Esoln+Gnes的值在PCM计算后在log文件中作为Total free energy in solution: with all non electrostatic terms 项自动给出，而deltaGcorr_gas是在气相下作freq计算后得到的值。注意，在这里溶液和气相下溶质的结构需要分别在液相和气相重作结构优化，很有可能是不一样的。 这种方法意味着：溶液下溶质的总自由能需要在液相重作结构优化，要用大基组吗？ 如果只用气相几何，考虑到溶剂对吉布斯自由能的校正，该如何计算溶液下溶质的总自由能呢？ |
4楼2012-05-01 08:32:43
richen0564
金虫 (著名写手)
- 应助: 2 (幼儿园)
- 金币: 1272.7
- 散金: 324
- 红花: 13
- 帖子: 1690
- 在线: 364小时
- 虫号: 351842
- 注册: 2007-04-22
- 性别: GG
- 专业: 有机分子功能材料化学
5楼2012-05-03 12:26:12
楼主现在应该明白了吧、我现在写文章也在这个问题上很困扰、请问最后的吉布斯自由能该如何表示?是不是大基组计算单点能+Thermal correction to Gibbs Free Energy?
6楼2012-09-04 09:47:28
幸福男人
金虫 (小有名气)
- 应助: 2 (幼儿园)
- 金币: 1202.7
- 帖子: 172
- 在线: 119.3小时
- 虫号: 1418413
- 注册: 2011-09-26
- 性别: GG
- 专业: 理论和计算化学
7楼2012-09-05 05:40:19
8楼2012-09-10 08:35:44