24小时热门版块排行榜    

CyRhmU.jpeg
查看: 1774  |  回复: 9

xiongb123456

铜虫 (小有名气)


[交流] 高斯可以计算升温下的反应吗?

虫友们 我最近在计算一个C-H键活化的反应
其中有分别有两种金属参与反应 据文献报道
其中一种金属可以引发C-H键活化的发生
而另一种在室温下不可以发生C-H键活化(只有在高温的情况下才可以)

正巧我的计算结果与这个实验相符
室温下 第一种金属可以发生C-H活化 反应前后 ΔG<0
而另一种金属 反应前后 ΔG>0
分子的吉布斯自由能是用(溶剂化能+吉布斯自由能的校正项)校正出来的

现在问题来了!我希望能把那种室温下不能发生反应的金属 进行高温计算!
从而分析是否升温后 会ΔG<0;
那么我应该如何在高斯的输入文件中 体现出升温 温度的改变呐?
因为我还算了溶剂化能 那么在溶剂化情况下 也可以升温吗?
高斯说明书上就简单地说了 在命令行添加temperature=?
我觉得这种方法好像不是很合理

最后对于这种计算不同温度的反应! 虫友们有什么好办法呐?
PS:我用的是高斯03
回复此楼

» 猜你喜欢

» 本主题相关商家推荐: (我也要在这里推广)

» 本主题相关价值贴推荐,对您同样有帮助:

» 抢金币啦!回帖就可以得到:

查看全部散金贴

已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

gongxd325

至尊木虫 (文坛精英)



xiongb123456(金币+1): 谢谢参与
Gaussian只能基于0K的结构获得不同温度的热力学函数
2楼2012-10-16 18:06:56
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

狼之一族

新虫 (初入文坛)



xiongb123456(金币+1): 谢谢参与
引用回帖:
2楼: Originally posted by gongxd325 at 2012-10-16 18:06:56
Gaussian只能基于0K的结构获得不同温度的热力学函数

那么溶剂效应 只算了一个单点 可以获得不同温度下的溶剂化能吗?
3楼2012-10-16 23:18:39
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

caby张

新虫 (初入文坛)



xiongb123456(金币+1): 谢谢参与
我觉得有一个笨办法,如果你的分子小,你可以像做实验那样,在每升高10度(看你升高的梯度决定)就算一次溶剂化能,分析能量变化
4楼2012-10-17 10:01:47
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

xiongb123456

铜虫 (小有名气)


引用回帖:
4楼: Originally posted by caby张 at 2012-10-17 10:01:47
我觉得有一个笨办法,如果你的分子小,你可以像做实验那样,在每升高10度(看你升高的梯度决定)就算一次溶剂化能,分析能量变化

嗯 那我可以拿水分子试试先哈
不过我记得我在计算溶剂化能的输入文件中 加入temperature=?时
貌似 溶剂化能 没有改变的说
5楼2012-10-17 10:25:04
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

caby张

新虫 (初入文坛)


溶剂化能是不变,你最后的G不是加上了气相的校正值了的嘛,你在气相下算不同温度的单点能就OK了
6楼2012-10-17 10:43:11
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

xiongb123456

铜虫 (小有名气)


引用回帖:
6楼: Originally posted by caby张 at 2012-10-17 10:43:11
溶剂化能是不变,你最后的G不是加上了气相的校正值了的嘛,你在气相下算不同温度的单点能就OK了

嗯 溶剂化能是不变的 我刚刚用了一个小分子甲烷 算他在甲醇中不同温度的溶剂化能 发现溶剂化能没有改变!
在气相下算不同温度的单点能? 那输入文件是不是应该写成
#p 基组 freq temperature=?  不需要opt了!
7楼2012-10-17 11:01:27
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

caby张

新虫 (初入文坛)


你这个结构是优化好的就直接算不同温度下的频率就行了,不用OPT
8楼2012-10-17 11:46:21
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

xiongb123456

铜虫 (小有名气)


引用回帖:
8楼: Originally posted by caby张 at 2012-10-17 11:46:21
你这个结构是优化好的就直接算不同温度下的频率就行了,不用OPT

如果按照你说的 那么只要将已经优化好的分子的 chk文件
进行freqchk 就可以了!
高斯说明书210页
freqchk 工具用于从 checkpoint 文件中获取频率和热化学数据,以及可选的温度,气压,换算因子,和/或同位素取代等信息。
直接freqchk solvent.chk   Checkpoint 文件名可位于命令行
Write Hyperchem files?  n  
Temperature (K)? [0=>298.15]  300   改变温度
Pressure (Atm)? [0=>1 atm]  1.5    改变气压
Scale factor for freqs during thermochem? [0=>1/1.12]  1    不换算
Do you want to use the principal isotope masses? [Y]: n  
For each atom, give the integer mass number.
In each case, the default is the principal isotope.
Atom number 1, atomic number 8: [16]    回车,接受默认值
Atom number 2, atomic number 1: [1] 2    以整数形式指定同位素质量
...
以下的频率输出内容同上。注意如果定义了换算因子,只会反映在热化学数据中,频率不做换算。
9楼2012-10-17 13:02:45
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

luolun2008

至尊木虫 (著名写手)



xiongb123456(金币+1): 谢谢参与
引用回帖:
9楼: Originally posted by xiongb123456 at 2012-10-17 13:02:45
如果按照你说的 那么只要将已经优化好的分子的 chk文件
进行freqchk 就可以了!
高斯说明书210页
freqchk 工具用于从 checkpoint 文件中获取频率和热化学数据,以及可选的温度,气压,换算因子,和/或同位素取代 ...

楼主回答的很详细。
10楼2012-10-18 14:19:34
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 xiongb123456 的主题更新
普通表情 高级回复(可上传附件)
信息提示
请填处理意见