应《网络安全法》要求，自2017年10月1日起，未进行实名认证将不得使用互联网跟帖服务。为保障您的帐号能够正常使用，请尽快对帐号进行手机号验证，感谢您的理解与支持！

24小时热门版块排行榜    

小木虫论坛-学术科研互动平台 » 计算模拟区 » 量子化学 » Gaussian » 高斯可以计算升温下的反应吗？
101/1返回列表
查看: 2423  |  回复: 9
只看楼主 @他人 存档 新回复提醒 (忽略) 收藏    在APP中查看

xiongb123456

铜虫 (小有名气)

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

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

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

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

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

» 猜你喜欢

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

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

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

查看全部散金贴
1楼 2012-10-16 16:52:08
已阅   回复此楼   关注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 的主题更新
101/1返回列表
普通表情 龙 兔 虎 猫 高级回复 (可上传附件)
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[基金申请] 能否退出参与的面上项目解除限项 +19 koalala 2026-08-24 21/1050 2026-08-24 16:48 by Wei_ren
[基金申请] 没有任何消息-是不是就凉了 +3 图啦图啦 2026-08-24 3/150 2026-08-24 16:43 by GZUZM
[基金申请] 建议基金发布提前给出明确的时间点 +13 kulium 2026-08-21 16/800 2026-08-24 16:27 by superceng
[基金申请] 如果此刻你正在为国基感到焦虑，不妨来听听这首《基金之外》 +4 scalable 2026-08-24 4/200 2026-08-24 16:22 by truewz
[基金申请] 今日不放榜？网传国自然预计 8 月 27 日可查结果 +17 医学老男孩 2026-08-20 21/1050 2026-08-24 14:21 by refreshing11
[基金申请] 我面上完蛋了 +8 且听虎啸 2026-08-20 9/450 2026-08-24 13:52 by zzuzxg
[基金申请] 2026国自然函评费到账 +17 羊腰板 2026-08-21 18/900 2026-08-24 12:03 by iaeyuan
[基金申请] 范进中举一文的中心思想 +6 炎黄贵胄 2026-08-22 7/350 2026-08-24 11:58 by 6543yes
[基金申请] 朋友圈看到的 +7 wangzilk 2026-08-18 9/450 2026-08-24 10:50 by cmrandy
[基金申请] 让我中一个面上吧！ +13 大萍1987 2026-08-20 16/800 2026-08-24 10:23 by 太傻了
[教师之家] 跳槽后在研项目怎么办？ +5 简单化xn 2026-08-22 10/500 2026-08-23 12:38 by 简单化xn
[基金申请] 只有每年这种时候来逛逛小木虫 +24 yaoyewhu2008 2026-08-20 26/1300 2026-08-22 17:43 by kammury
[基金申请] filecode，4个jtjc了 +13 ziyangfang 2026-08-19 16/800 2026-08-22 17:08 by WH3796
[基金申请] 人气不行了 +8 fansofjerry 2026-08-21 8/400 2026-08-22 16:30 by zyqchem
[基金申请] 今天基金会出结果吗？20260819 +16 kkkl_v 2026-08-19 17/850 2026-08-22 16:12 by 阿布Abu
[基金申请] 时间戳又变了 +13 wuchongjun 2026-08-20 19/950 2026-08-21 17:21 by 紫杉醇
[基金申请] 应该是下周三26日公布了吧？ +4 哈哈蛤？ 2026-08-21 4/200 2026-08-21 10:58 by Vivilian
[基金申请] 今天放榜没戏了吧 +9 yuleib84 2026-08-19 11/550 2026-08-21 10:06 by gltch
[基金申请] 基金啊基金 +4 longfie172 2026-08-20 4/200 2026-08-21 08:58 by mark mao
[基金申请] 明天放榜？ +5 Shxjjxjkx 2026-08-18 5/250 2026-08-18 18:14 by -大大大大大-
信息提示
关闭
请填处理意见
关闭