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

24小时热门版块排行榜    

Znn3bq.jpeg
小木虫论坛-学术科研互动平台 » 计算模拟区 » 量子化学 » Gaussian » 请教freq=readisotopes计算同位素效应的输入文件设置？
51/1返回列表
查看: 3812  |  回复: 13
只看楼主 @他人 存档 新回复提醒 (忽略) 收藏    在APP中查看
当前只显示满足指定条件的回帖，点击这里查看本话题的所有回帖

些小

铜虫 (初入文坛)

[求助] 请教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

» 猜你喜欢

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

» 本主题相关价值贴推荐，对您同样有帮助:
1楼 2013-12-02 16:21:42
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

些小

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

谢谢你的解答，很清楚。
我按照你的方法重新建立了输入文件。果然不报错了。但是计算结果出来，和没有定义同位素是一样的。我想计算H和D的不同IR谱图，我看文献，说freq结果是有区别的，但是我算出来是 一样的。请问你知道是怎么回事吗？？？
谢谢！
赞一下 回复此楼
5楼2013-12-06 19:24:01
已阅   回复此楼   关注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结果是有区别的，但是我算出来是 一样的。请问你知道是怎么回事吗？？？
谢谢！
赞一下(1人) 回复此楼
4楼2013-12-06 19:23:45
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
查看全部 14 个回答
普通表情 龙 兔 虎 猫
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 一志愿华中农业071010，320求调剂 +17 困困困困坤坤 2026-04-14 19/950 2026-04-17 20:08 by 关一盏灯cd
[考研] 304求调剂 +7 castLight 2026-04-16 7/350 2026-04-17 20:05 by 关一盏灯cd
[考研] 320求调剂 +4 深郊akm 2026-04-17 4/200 2026-04-17 17:32 by 一切OK
[考研] 295分求调剂 +5 ?要上岸? 2026-04-17 5/250 2026-04-17 16:51 by fenglj492
[考研] 297，工科调剂?河南农业大学本科 +14 河南农业大学-能 2026-04-14 14/700 2026-04-16 14:41 by dingyanbo1
[考研] 26药学专硕105500求调剂 +6 喽哈加油 2026-04-13 7/350 2026-04-16 14:31 by zhouxiaoyu
[基金申请] RY：中国产出的科学垃圾论文，绝对数量和比例都世界第一 +7 zju2000 2026-04-14 18/900 2026-04-16 11:36 by 欢乐颂叶蓁
[考研] 322求调剂 +8 123安康 2026-04-12 15/750 2026-04-16 11:07 by Espannnnnol
[考研] 通信工程求调剂！！！ +6 zlb770521 2026-04-14 6/300 2026-04-15 20:00 by 学员JpLReM
[考研] 调剂求收留 +34 果然有我 2026-04-10 35/1750 2026-04-15 13:05 by 西北望—风沙
[考研] 复试调剂 +21 积极向上； 2026-04-10 23/1150 2026-04-15 12:50 by 西北望—风沙
[考研] 求调剂 +12 何气正 2026-04-13 13/650 2026-04-14 14:47 by zs92450
[考研] 考研英一数一338分 +9 长江大学东校区 2026-04-13 10/500 2026-04-14 00:41 by 王珺璞
[考研] 求调剂，985材料与化工348分 +9 涵竹刘 2026-04-11 14/700 2026-04-13 22:26 by 涵竹刘
[考研] 302求调剂 +10 易！? 2026-04-13 10/500 2026-04-13 19:04 by lbsjt
[考研] 344 材料专业 求调剂211 无地域要求 +8 hualkop 2026-04-11 8/400 2026-04-12 22:24 by fqwang
[考研] 291求调剂 +8 关忆北. 2026-04-11 8/400 2026-04-12 09:32 by 逆水乘风
[考研] 270求调剂 +14 杨乐369 2026-04-11 14/700 2026-04-11 20:16 by 蓝云思雨
[考研] 调剂求助 +6 果然有我 2026-04-11 7/350 2026-04-11 16:22 by 明月此时有
[考研] 求调剂，一志愿大连理工大学354分 +5 雨声余生 2026-04-11 6/300 2026-04-11 16:12 by 雨声余生
信息提示
关闭
请填处理意见
关闭