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

24小时热门版块排行榜    

小木虫论坛-学术科研互动平台 » 计算模拟区 » 量子化学 » Gaussian » 请问Gaussian 03可以算同位素吗？该怎么设置呢？
71/1返回列表
查看: 1959  |  回复: 6
只看楼主 @他人 存档 新回复提醒 (忽略) 收藏    在APP中查看

独孤立啊

铜虫 (正式写手)

[求助] 请问Gaussian 03可以算同位素吗？该怎么设置呢？ 已有1人参与

如题，我在Gaussian 03手册上看到了可以设置同位素的信息，我想做H2O的氘代 HOD，尝试了两种方法，都没有成功：

尝试的方法一：只在input文件中增加了同位素设置，H(Iso=2)。结果是计算结束后，查看原子质量，还是1并不是2，并没哟进行同位素取代计算。
尝试的方法二：在input文件中增加了同位素设置，H(Iso=2)，并增加了关键词freq=readisotopes。结果是报L101错误，详细报错信息是WANTED A FLOATING POINT NUMBER AS INPUT.
  FOUND A STRING AS INPUT.
，并在log文件中出现一些乱码。
Gaussian手册上也说明同位素可以设置成整数，并且我也尝试过改为Iso=2.0，依然报错。

期待有做过的小伙伴赐教

谢谢
回复此楼

» 猜你喜欢
1楼 2017-10-21 09:45:07
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

beefly

专家顾问 (职业作家)

地沟油冶炼专家

【答案】应助回帖

★ ★ ★ ★ ★ ★ ★ ★ ★ ★
感谢参与，应助指数 +1
zhou2009: 金币+5 2017-10-22 08:24:37
独孤立啊: 金币+5, ★★★很有帮助 2017-10-23 13:19:30
gaussian频率计算保存fchk文件，然后用我的UniMoVib程序，并在输入里指定同位素质量

https://github.com/zorkzou/UniMoVib
赞一下(1人) 回复此楼
beefly《西太平洋大学现代英汉词典》[bi:fli]牛肉一般地
2楼2017-10-21 16:55:22
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

独孤立啊

铜虫 (正式写手)
引用回帖:
2楼: Originally posted by beefly at 2017-10-21 16:55:22
gaussian频率计算保存fchk文件，然后用我的UniMoVib程序，并在输入里指定同位素质量

https://github.com/zorkzou/UniMoVib

感谢您的回复及提供可以计算的程序，请问直接用Gaussian 03不能计算吗？有什么限制呢？看到手册里有介绍但是就是没有计算起来，不知道是不是版本的限制。
赞一下 回复此楼
3楼2017-10-23 13:19:22
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

beefly

专家顾问 (职业作家)

地沟油冶炼专家
引用回帖:
3楼: Originally posted by 独孤立啊 at 2017-10-23 13:19:22
感谢您的回复及提供可以计算的程序，请问直接用Gaussian 03不能计算吗？有什么限制呢？看到手册里有介绍但是就是没有计算起来，不知道是不是版本的限制。...

gaussian自身也可以，但是不知道你的输入为什么不行。也许是版本太低？
赞一下(1人) 回复此楼
beefly《西太平洋大学现代英汉词典》[bi:fli]牛肉一般地
4楼2017-10-23 20:42:17
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

独孤立啊

铜虫 (正式写手)
引用回帖:
4楼: Originally posted by beefly at 2017-10-23 20:42:17
gaussian自身也可以，但是不知道你的输入为什么不行。也许是版本太低？...

这是我的输入文件，请问有什么问题吗？
%nprocshared=2
%mem=500MB
# opt=(calcfc,ts,noeigentest) freq=readisotopes rb3lyp/6-311++g(d,p) nosymm pop=npa

6h2o-GP-++-D-TS-LL-1019

0 1
O                  1.15351800   -1.30301300   -1.45714400
H                  1.61299900   -0.47964800   -1.22115000
H                  1.12191500   -1.77451500   -0.61015200
O                  0.98527100   -1.26551600    1.38147900
H                  1.34893100   -0.37695600    1.21575100
H(Iso=2)           0.64042600   -1.26025600    2.27896200
O                  1.77933400    1.08889200    0.06954300
H                  2.45900400    1.76474300    0.14504600
H                  0.90703700    1.55373500   -0.02347400
O                 -1.39518000   -0.36527700   -1.61070000
H                 -0.50103200   -0.77765500   -1.68796800
H                 -1.87176100   -0.58533000   -2.41611600
O                 -0.75017500    1.96287600   -0.17445200
H                 -1.06961500    1.34633300   -0.86093400
H                 -1.19073800    1.63050400    0.62524900
O                 -1.86075700    0.00120100    1.63145000
H(Iso=2)          -2.81066800   -0.15392400    1.61541000
H(Iso=2)          -1.49643300   -0.59223600    0.95946100

输出报错如下：
End of file reading temperature and pressure.
  WANTED A FLOATING POINT NUMBER AS INPUT.
  FOUND A STRING AS INPUT.
Atom number 1, atomic number 8: [16] 
Error termination via Lnk1e in /home/abc/opt/g03/l101.exe at Thu Oct 19 11:40:00 2017.
Job cpu time:  0 days  0 hours  0 minutes  0.8 seconds.
File lengths (MBytes):  RWF=     13 Int=      0 D2E=      0 Chk=      1 Scr=      1
赞一下 回复此楼
5楼2017-10-24 07:37:39
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

beefly

专家顾问 (职业作家)

地沟油冶炼专家
引用回帖:
5楼: Originally posted by 独孤立啊 at 2017-10-24 07:37:39
这是我的输入文件，请问有什么问题吗？
%nprocshared=2
%mem=500MB
# opt=(calcfc,ts,noeigentest) freq=readisotopes rb3lyp/6-311++g(d,p) nosymm pop=npa

6h2o-GP-++-D-TS-LL-1019

0 1
O            ...

你要是用freq=readisotopes，就要在结构之后，按照标准格式提供温度压强和原子量，这时候iso不起作用
赞一下(1人) 回复此楼
beefly《西太平洋大学现代英汉词典》[bi:fli]牛肉一般地
6楼2017-10-24 09:30:07
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

匿名

本帖仅楼主可见
赞一下
7楼2018-04-27 11:31:12
已阅   申请QC强帖   回复此楼   编辑   查看我的主页
相关版块跳转 我要订阅楼主 独孤立啊 的主题更新
71/1返回列表
普通表情 龙 兔 虎 猫
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[基金申请] 2026国自然函评费到账 +16 羊腰板 2026-08-21 17/850 2026-08-23 23:02 by tianxiaochun
[基金申请] 让我中一个面上吧！ +12 大萍1987 2026-08-20 15/750 2026-08-23 20:44 by 大萍1987
[基金申请] 2026年的国家社科基金项目通讯评审的新规则与新动向、新挑战 +4 process2012 2026-08-23 5/250 2026-08-23 19:58 by jurkat.1640
[教师之家] 跳槽后在研项目怎么办？ +5 简单化xn 2026-08-22 10/500 2026-08-23 12:38 by 简单化xn
[基金申请] 今天放榜吗？ +15 布布和一二 2026-08-19 16/800 2026-08-23 09:55 by 张春生
[基金申请] filecode，4个jtjc了 +13 ziyangfang 2026-08-19 16/800 2026-08-22 17:08 by WH3796
[基金申请] 今天基金会出结果吗？20260819 +16 kkkl_v 2026-08-19 17/850 2026-08-22 16:12 by 阿布Abu
[基金申请] 时间戳今天，20号变了 +5 archvillain 2026-08-20 5/250 2026-08-22 06:12 by hui_daxiao
[基金申请] 建议基金发布提前给出明确的时间点 +10 kulium 2026-08-21 13/650 2026-08-21 21:58 by alongwaytogo
[基金申请] 今日不放榜？网传国自然预计 8 月 27 日可查结果 +16 医学老男孩 2026-08-20 20/1000 2026-08-21 21:13 by Ldrop2023
[基金申请] 看来今天不会放榜了？ +8 chengyan1220 2026-08-21 11/550 2026-08-21 17:52 by dcqxinyang
[基金申请] 时间戳又变了 +13 wuchongjun 2026-08-20 19/950 2026-08-21 17:21 by 紫杉醇
[基金申请] 我面上完蛋了 +7 且听虎啸 2026-08-20 8/400 2026-08-21 12:31 by 酷酷墨镜
[基金申请] 应该是下周三26日公布了吧？ +4 哈哈蛤？ 2026-08-21 4/200 2026-08-21 10:58 by Vivilian
[基金申请] 基金啊基金 +4 longfie172 2026-08-20 4/200 2026-08-21 08:58 by mark mao
[基金申请] 估计是周四 +3 archvillain 2026-08-18 3/150 2026-08-21 01:48 by jnhyjjm
[基金申请] 时间戳变了，能看出什么问题？ +18 基诺咪客 2026-08-17 23/1150 2026-08-20 17:19 by Godzela
[基金申请] 朋友圈看到的 +6 wangzilk 2026-08-18 8/400 2026-08-19 10:55 by Haru815
[基金申请] 明天放榜？ +5 Shxjjxjkx 2026-08-18 5/250 2026-08-18 18:14 by -大大大大大-
[基金申请] 今天维护系统维护 祝所有人 高中 +8 gjjjzhong 2026-08-18 9/450 2026-08-18 13:01 by 家与远方
信息提示
关闭
请填处理意见
关闭