| 查看: 3782 | 回复: 13 | |||
| 当前只显示满足指定条件的回帖，点击这里查看本话题的所有回帖 | |||
[求助]
请教freq=readisotopes计算同位素效应的输入文件设置？
|
|||
|
用高斯03计算同位素效应。 计算CH2=CH-CH2-CD3 输入关键词freq=readisotopes. 下面分子坐标后加温度压强，再加原子量。 总是出错，显示failed to get reasonable mass for atom 1 atomic number 6。 我再三核对了输入文件的元素顺序应该没有错。 输入文件的顺序为 c c 1 c 2 c 3 H 1 H 1 H 2 H 3 H 3 H 4 H 4 H 4 我后面定义原子数的时候写的如下 12 12 12 12 1 1 1 1 1 2 2 2 计算就总是出错。翻看旧帖。有人说readisotopes的关键词如果只是定义压强和温度已经废弃不用。但是我要算同位素效应。请各位指教。 另外如果我把readisotopes的关键词去掉，只是在后面写上 12 12 12 12 1 1 1 1 1 2 2 2 计算可以通过。但是freq结果和没有定义D用H时是一样的。翻开旧帖应该是freq结果有区别。 所以请教如何定义同位素。 |
回复此楼» 收录本帖的淘帖专辑推荐
 Materials Modeling |
» 猜你喜欢
【211博士招生】环境化学、地学、毒理方向，擅长污染物迁移转化降解、计算模拟者优先
已经有17人回复
-
求助如何利用vasp计算磁光克尔效应
已经有1人回复
-
物理化学论文润色/翻译怎么收费?
已经有291人回复
-
070300总分298求调剂
已经有2人回复
-
一志愿厦门大学化学学硕307求调剂
已经有19人回复
-
理论与计算化学方向调剂——初试388分（数学144）
已经有2人回复
-
调剂招生
已经有7人回复
-
河南理工大学化工学院招收调剂生啦！
已经有1人回复
-
硕士研究生调剂
已经有0人回复
-
投稿期刊求助
已经有1人回复
-
江西师大化学工程学院调剂今晚0点开放！ 缺额100个，速冲！
已经有0人回复
» 本主题相关商家推荐: (我也要在这里推广)
» 本主题相关价值贴推荐，对您同样有帮助:
- 高斯计算时假设分子在什么温度下？
已经有9人回复
- Gaussian频率计算，如何加温度？
已经有13人回复
- 【原创】GAUSSION计算常见错误及解决方案
已经有32人回复
zhangmt
至尊木虫 (著名写手)
我叫MT
- QC强帖: 5
- 应助: 99 (初中生)
- 金币: 6961.8
- 散金: 10406
- 红花: 49
- 帖子: 1761
- 在线: 763.2小时
- 虫号: 880392
- 注册: 2009-10-22
- 性别: GG
- 专业: 理论和计算化学
【答案】应助回帖
★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★
感谢参与，应助指数 +1
gkf高: 金币+5 2013-12-03 10:25:40
些小: 金币+10, ★★★很有帮助, 谢谢 2013-12-03 15:18:37
感谢参与，应助指数 +1
gkf高: 金币+5 2013-12-03 10:25:40
些小: 金币+10, ★★★很有帮助, 谢谢 2013-12-03 15:18:37
|
1、readisotopes不是不可以用，但是属于已被废弃，不建议使用。如果你依然喜欢用，还是可以用的，包括你算同位素效应，用法不变。 2、新的使用方法（其实已经不新了），除了温度、压强新增了关键字，对于同位素则更简单，就是在你写坐标的时候，在相应的元素名称后面写上(iso=n)即可，举例而言： #p hf/3-21g freq title 0 1 O(iso=17) 0.0 0.0 0.0 H 1.0 0.0 0.0 H(iso=2) 0.0 0.0 1.0 ========== iso后面可以接两种格式的数据，如果是实数，则严格按照你输入的数据执行，比如，你写(iso=17.000)，那么计算的时候O同位素就采用17.000，如果你写的是整数，比如上面的例子，那么O同位素采用的是实际相应的精确同位素质量17O(16.99913)，2H(2.01410), 同时你没有特别说明的那个H，则自动采用1H的精确质量(1.00783)，我这么解释，你清楚了没？ |
3楼2013-12-03 09:17:28
【答案】应助回帖
★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★
感谢参与，应助指数 +1
gmy1990: 金币+2 2013-12-02 19:38:09
些小: 金币+10, ★★★很有帮助, 谢谢 2013-12-03 15:17:32
感谢参与，应助指数 +1
gmy1990: 金币+2 2013-12-02 19:38:09
些小: 金币+10, ★★★很有帮助, 谢谢 2013-12-03 15:17:32
|
-------------------------------------------------------------------------------------- # freq=(noraman,readisotopes) ub3lyp/6-31++G(d,p) opt Title Card Required 1 2 O 4.53438900 -2.22112600 0.11517600 H 0.00000000 0.00000000 0.00000000 298.15,1.0 16 2 ---------------------------------- 温度和压强应该按照上面的那种格式来输入，如果还是不正确的话你换做笛卡尔坐标试试~ |
2楼2013-12-02 18:14:18
4楼2013-12-06 19:23:45
5楼2013-12-06 19:24:01
