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

24小时热门版块排行榜    

小木虫论坛-学术科研互动平台 » 计算模拟区 » 量子化学 » Gaussian » 请问Gaussian 03可以算同位素吗？该怎么设置呢？
71/1返回列表
查看: 1503  |  回复: 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返回列表
普通表情 龙 兔 虎 猫
信息提示
关闭
请填处理意见
关闭