当前位置: 首页 > 量子化学 >请问Gaussian 03可以算同位素吗?该怎么设置呢?

请问Gaussian 03可以算同位素吗?该怎么设置呢?

作者 独孤立啊
来源: 小木虫 300 6 举报帖子
+关注

如题,我在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,依然报错。

期待有做过的小伙伴赐教

谢谢 返回小木虫查看更多

今日热帖
  • 精华评论
  • beefly

    gaussian频率计算保存fchk文件,然后用我的UniMoVib程序,并在输入里指定同位素质量

    https://github.com/zorkzou/UniMoVib

  • 独孤立啊

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

    https://github.com/zorkzou/UniMoVib

    感谢您的回复及提供可以计算的程序,请问直接用Gaussian 03不能计算吗?有什么限制呢?看到手册里有介绍但是就是没有计算起来,不知道是不是版本的限制。

  • beefly

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

    gaussian自身也可以,但是不知道你的输入为什么不行。也许是版本太低?

  • 独孤立啊

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

  • 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不起作用

  • 梁子XDU

    楼主你好,我非化学专业的,想问下Gaussian 这个软件能不能计算双原子分子某特定电子态时的所用振动能级和转动能级能量大小?

猜你喜欢
下载小木虫APP
与700万科研达人随时交流
  • 二维码
  • IOS
  • 安卓