计算化学公社

 找回密码 Forget password
 注册 Register
Views: 10736|回复 Reply: 3
打印 Print 上一主题 Last thread 下一主题 Next thread

[GAMESS-US] GAMESS软件求助

[复制链接 Copy URL]

58

帖子

0

威望

612

eV
积分
670

Level 4 (黑子)

跳转到指定楼层 Go to specific reply
楼主
本帖最后由 Tinker 于 2014-12-10 23:51 编辑

求助各位大大,有两个问题,

1. 第一个问题是关于GAMESS里边怎么使用混合基组,比如高斯里边我可以写,
B3LYP/GENECP 。。。。

下边写
Ir 0
Lanl2dz
****
C H 0
6-31g*
****

Ir 0
Lanl2dz


而在GAMESS,他们给我的输入文件是这样的 在分子信息后边写价电子基组,之后再写个$ECP去写赝势
IRIDIUM    77.0      9.277565000         1.202463000         6.804552000
S   3
  1      2.3500000             -1.6784642        
  2      1.5820000              2.0952553        
  3      0.5018000              0.4162934        
S   4
  1      2.3500000              1.6464467        
  2      1.5820000             -2.2748150        
  3      0.5018000             -1.0494357        
  4      0.2500000              1.2167791        
S   1
  1      0.0598000              1.0000000        
P   3
  1      2.7920000             -0.3889212        
  2      1.5410000              0.9077516        
  3      0.5285000              0.4691443        
P   2
  1      0.5100000             -0.1170669        
  2      0.0980000              1.0489002        
P   1
  1      0.0290000              1.0000000        
D   2
  1      1.2400000              0.5087022        
  2      0.4647000              0.5862102        
D   1
  1      0.1529000              1.0000000   


$ECP
  IR-ECP GEN     60     4
  5      ----- g-ul potential     -----
     -0.1578014        1    823.5880147        
  -1517.5270446        2    364.6613336        
   -316.5306529        2     55.7082801        
    -91.8880941        2     12.0464544        
     -9.2241773        2      3.5120610        
  6      ----- s-ul potential     -----
      3.1578014        0    188.0490770        
     26.8322577        1    340.4194712        
    800.4250007        2    128.2373673        
    369.4050683        2     33.8644961        
    242.4171899        2      4.7560005        
   -118.2173282        2      3.9649974        
  5      ----- p-ul potential     -----
      2.1578014        0    289.7291139        
     61.9678610        1     87.4633789        
    269.0581986        2     30.4363766        
    231.1654793        2      4.0553412        
   -133.6952667        2      3.5525341        
  5      ----- d-ul potential     -----
      3.1578014        0    136.4017106        
     45.9349803        1     95.0776925        
    359.0344668        2     49.2258410        
    176.4740119        2     15.0874145        
     54.5155286        2      4.0405764        
  5      ----- f-ul potential     -----
      3.9546197        0    127.3507908        
     52.9773655        1     66.2364374        
    274.8643383        2     34.4299229        
    137.2047338        2     10.1995721        
     14.8633305        2      2.5409702        
O-ECP NONE
C-ECP  NONE
C-ECP  NONE
C-ECP  NONE
H NONE
H NONE
H NONE
H NONE
H NONE
H NONE
H NONE
H NONE
H NONE
O-ECP NONE
$END


一定要通过这种方法写混合基组吗?

还有一个问题,最近在练习使用GAMESS做MCSCF的计算,我师兄的做法是,是先做HF计算,得到的轨道信息作为MCSCF的初始轨道,贴在MCSCF输入文件的下边,是否一定需要这样?


44

帖子

0

威望

663

eV
积分
707

Level 4 (黑子)

2#
发表于 Post on 2014-12-11 00:15:17 | 只看该作者 Only view this author
GAMESS里已经定义好的ECP很少,大部分要像你所说的方法输入,具体有哪些你可以看手册。

是的。需要guess=moread 然后用$VEC提供初始轨道,就像Gaussian里面guess=read,不过轨道信息在check文件中,而不是写在输入文件里。做casscf也常常需要guess=(read,alter),gamess貌似没有alter这个命令,可以手动交换。

评分 Rate

参与人数
Participants 2
eV +6 收起 理由
Reason
Tinker + 3 谢谢
sobereva + 3

查看全部评分 View all ratings

58

帖子

0

威望

612

eV
积分
670

Level 4 (黑子)

3#
 楼主 Author| 发表于 Post on 2014-12-11 01:42:35 | 只看该作者 Only view this author
jjjspring 发表于 2014-12-11 00:15
GAMESS里已经定义好的ECP很少,大部分要像你所说的方法输入,具体有哪些你可以看手册。

是的。需要guess ...

谢谢!我再试试。

找了 两个说明先贴在这里
http://www.spec.org/hpc96/docs/R ... s/gamess/input.html
http://myweb.liu.edu/~nmatsuna/gamess/input/ORMAS.html

729

帖子

21

威望

5233

eV
积分
6382

Level 6 (一方通行)

4#
发表于 Post on 2015-1-2 09:35:15 | 只看该作者 Only view this author
基组部分:可以设置rungms脚本中的外部基组文件变量,然后就可以在$BASIS中直接引用自定义基组名称。不过近两年的rungms脚本把这个环境变量丢掉了,需要用两年前的rungms脚本。

ECP部分:没办法,目前还得对每个原子(不是元素!)挨个定义

评分 Rate

参与人数
Participants 1
eV +2 收起 理由
Reason
sobereva + 2

查看全部评分 View all ratings

本版积分规则 Credits rule

手机版 Mobile version|北京科音自然科学研究中心 Beijing Kein Research Center for Natural Sciences|京公网安备 11010502035419号|计算化学公社 — 北京科音旗下高水平计算化学交流论坛 ( 京ICP备14038949号-1 )|网站地图

GMT+8, 2025-8-12 21:51 , Processed in 0.203037 second(s), 21 queries , Gzip On.

快速回复 返回顶部 返回列表 Return to list