24小时热门版块排行榜    

查看: 3727  |  回复: 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的回帖

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的回帖
查看全部 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的回帖

些小

铜虫 (初入文坛)

引用回帖:
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的回帖

些小

铜虫 (初入文坛)

引用回帖:
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的回帖
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 一志愿吉林大学材料学硕321求调剂 +11 Ymlll 2026-03-18 15/750 2026-03-20 19:40 by 丁丁*
[考研] 290求调剂 +4 ^O^乜 2026-03-19 4/200 2026-03-20 11:51 by 学员8dgXkO
[考研] 一志愿北京化工大学0703化学318分,有科研经历,求调剂 +3 一瓶苯甲酸 2026-03-14 3/150 2026-03-19 15:17 by 尽舜尧1
[考研] 一志愿西安交通大学材料工程专业 282分求调剂 +5 枫桥ZL 2026-03-18 7/350 2026-03-19 14:52 by 功夫疯狂
[考研] 材料考研调剂 +3 xwt。 2026-03-19 3/150 2026-03-19 11:22 by w沐阳w
[考研] 0817调剂 +3 没有答案_ 2026-03-14 3/150 2026-03-19 09:51 by Xu de nuo
[考研] 一志愿华中科技大学,080502,354分求调剂 +4 守候夕阳CF 2026-03-18 4/200 2026-03-18 22:16 by li123456789.
[考研] 085601专硕,总分342求调剂,地区不限 +5 share_joy 2026-03-16 5/250 2026-03-18 14:48 by haxia
[考研] 268求调剂 +6 简单点0 2026-03-17 6/300 2026-03-18 09:04 by 无际的草原
[基金申请] 被我言中:新模板不强调格式了,假专家开始管格式了 +4 beefly 2026-03-14 4/200 2026-03-17 22:04 by 黄鸟于飞Chao
[硕博家园] 湖北工业大学 生命科学与健康学院-课题组招收2026级食品/生物方向硕士 +3 1喜春8 2026-03-17 5/250 2026-03-17 17:18 by ber川cool子
[考研] 材料专硕326求调剂 +6 墨煜姒莘 2026-03-15 7/350 2026-03-17 17:10 by ruiyingmiao
[考研] 有没有道铁/土木的想调剂南林,给自己招师弟中~ +3 TqlXswl 2026-03-16 7/350 2026-03-17 15:23 by TqlXswl
[考研] 一志愿南京大学,080500材料科学与工程,调剂 +4 Jy? 2026-03-16 4/200 2026-03-17 11:02 by gaoqiong
[考研] 304求调剂 +5 素年祭语 2026-03-15 5/250 2026-03-16 17:00 by 我的船我的海
[考研] 一志愿211 0703方向310分求调剂 +3 努力奋斗112 2026-03-15 3/150 2026-03-16 16:44 by houyaoxu
[考研] 26考研一志愿中国石油大学(华东)305分求调剂 +3 嘉年新程 2026-03-15 3/150 2026-03-15 13:58 by 哈哈哈哈嘿嘿嘿
[考研] 0856专硕279求调剂 +5 加油加油!? 2026-03-15 5/250 2026-03-15 11:58 by 2020015
[考研] 070305求调剂 +3 mlpqaz03 2026-03-14 4/200 2026-03-15 11:04 by peike
[考研] 本科南京大学一志愿川大药学327 +3 麦田耕者 2026-03-14 3/150 2026-03-14 20:04 by 外星文明
信息提示
请填处理意见