24小时热门版块排行榜    

查看: 4610  |  回复: 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结果有区别。
所以请教如何定义同位素。
回复此楼
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

游子8921

铁杆木虫 (知名作家)

游子8921

引用回帖:
3楼: Originally posted by zhangmt at 2013-12-03 09:17:28
1、readisotopes不是不可以用,但是属于已被废弃,不建议使用。如果你依然喜欢用,还是可以用的,包括你算同位素效应,用法不变。
2、新的使用方法(其实已经不新了),除了温度、压强新增了关键字,对于同位素则更 ...

向高手求助,文章审稿人让我计算Equilibrium isotope effects,在小木虫里搜到了kinetic isotope effect (KIE)的帖子,但是两者的计算公式不一样。也有用readisotopes计算同位素效应的帖子,也不确定是不是适应于Equilibrium isotope effects。
请问Equilibrium isotope effects该如何计算?用高斯可以计算吗?谢谢
请教freq=readisotopes计算同位素效应的输入文件设置?
Equilibrium isotope effects(EIE).jpg


请教freq=readisotopes计算同位素效应的输入文件设置?-1
kinetic isotope effect (KIE).png

为了自己的目标而奋斗,一切皆有可能
13楼2014-02-17 17:44:52
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
查看全部 14 个回答

marson

银虫 (正式写手)

【答案】应助回帖

★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★
感谢参与,应助指数 +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
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

zhangmt

至尊木虫 (著名写手)

我叫MT

【答案】应助回帖

★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★
感谢参与,应助指数 +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
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

些小

铜虫 (初入文坛)

引用回帖:
2楼: Originally posted by marson at 2013-12-02 18:14:18
--------------------------------------------------------------------------------------
# freq=(noraman,readisotopes) ub3lyp/6-31++G(d,p) opt

Title Card Required

1 2
O                  4.5343 ...

谢谢你的解答,很清楚。
我按照你的方法重新建立了输入文件。果然不报错了。但是计算结果出来,和没有定义同位素是一样的。我想计算H和D的不同IR谱图,我看文献,说freq结果是有区别的,但是我算出来是 一样的。请问你知道是怎么回事吗???
谢谢!
4楼2013-12-06 19:23:45
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[教师之家] 售SCI一区T0P文章,我:8.O55.1.O.54,科目全,可十急 +3 cqv1zKnzd7hr 2026-11-09 4/200 2026-11-10 12:52 by C6CRn66RYYXS
[论文投稿] 售SCI一区T0P文章,我:8.O.55.1.O.5.4,科目全,可+急 +3 cqv1zKnzd7hr 2026-11-09 4/200 2026-11-10 12:30 by C6CRn66RYYXS
[考研] 售SCI一区文章,我:8O5.5.1.O5.4,科目全,可伽急 +3 cqv1zKnzd7hr 2026-11-09 4/200 2026-11-10 11:46 by YHiw1gou1Eda
[论文投稿] 售SCI一区T0P文章,我:8.O.55.1.O.5.4,科目全,可+急 +3 cqv1zKnzd7hr 2026-11-09 4/200 2026-11-10 11:35 by YHiw1gou1Eda
[教师之家] 延安大学回应“副教授留绝笔信”:确已离世,正核实举报内容 +4 yexuqing 2026-11-04 6/300 2026-11-10 11:13 by susizheng
[教师之家] 售SCI一区文章,我:8.O.55.1.O.54,科目齐全,可伽急 +3 0TEMgxdT1eIw 2026-11-09 4/200 2026-11-10 11:02 by YHiw1gou1Eda
[教师之家] 售SCI一区T0P文章,我:8.O55.1.O.54,科目全,可十急 +3 zSbvVX9JFl9J 2026-11-09 5/250 2026-11-10 09:23 by 1yv6QACPoY1E
[考博] 售SCI一区文章,我:8.O.55.1.O.54,科目齐全,可伽急 +3 zSbvVX9JFl9J 2026-11-09 4/200 2026-11-10 01:29 by FChKX8meu1ds
[文学芳草园] 厉娜 +3 张丹枫1991 2026-11-09 5/250 2026-11-09 16:53 by 张丹枫1991
[考研] 售SCI-T0P文章,我:8O.5.5.1.O.54,科目齐全,可+急 +5 ttH3f1oA6ace 2026-11-04 15/750 2026-11-09 13:19 by Equinoxhua
[硕博家园] 售SCI一区T0P文章,我:8O.55.1.O.54,科目全,可伽急 +5 5sJFDKp2kwuI 2026-11-06 16/800 2026-11-09 11:47 by GwrQVVR6Oenx
[考博] 售SCI文章,我:8O.5.5.1O.54,科目全,可十急 +5 5sJFDKp2kwuI 2026-11-06 15/750 2026-11-09 11:41 by GwrQVVR6Oenx
[硕博家园] 售SCI一区文章,我:8O5.5.1.O5.4,科目全,可伽急 +4 5sJFDKp2kwuI 2026-11-06 10/500 2026-11-09 08:00 by IhnEQauGmT3K
[公派出国] 售SCI一区文章,我:8.O.55.1.O.54,科目齐全,可伽急 +4 5sJFDKp2kwuI 2026-11-06 8/400 2026-11-09 00:14 by l48hFd8Nwae2
[教师之家] 售SCI一区T0P文章,我:8.O.55.1.O.54,科目齐全,可+急 +4 5sJFDKp2kwuI 2026-11-06 8/400 2026-11-08 23:01 by l48hFd8Nwae2
[教师之家] 售SCI一区T0P文章,我:8.O.55.1.O.5.4,科目全,可+急 +5 5sJFDKp2kwuI 2026-11-06 9/450 2026-11-08 22:46 by l48hFd8Nwae2
[硕博家园] 售SCI一区文章,我:8O5.5.1.O5.4,科目全,可伽急 +4 ttH3f1oA6ace 2026-11-04 11/550 2026-11-06 16:10 by 5sJFDKp2kwuI
[硕博家园] 售SCI一区T0P文章,我:8O.55.1.O.5.4,科目齐全,可+急 +6 kvewwGtKktqB 2026-11-04 16/800 2026-11-06 05:00 by qL1cU46YIdK0
[考研] 售SCI一区T0P文章,我:8.O55.1.O.54,科目全,可十急 +4 ttH3f1oA6ace 2026-11-04 10/500 2026-11-06 02:04 by qL1cU46YIdK0
[硕博家园] 售SCI-T0P文章,我:8O.5.5.1.O.54,科目齐全,可+急 +3 ttH3f1oA6ace 2026-11-04 9/450 2026-11-06 02:02 by qL1cU46YIdK0
信息提示
请填处理意见