求gamess能量分解输入文件的格式

小菜鸟求关照

brothers 发表于 2015-2-2 16:52
我可不是老师，嘿嘿~gamess才学了不到半个月

您说的那个是哪一部分呢？

brothers

小菜鸟求关照 发表于 2015-2-2 16:56
您说的那个是哪一部分呢？

搜索basnam关键词，一定有这个例子

$contrlscftyp=rohf icharg=+3 mult=4 runtyp=gradient
         pp=read ispher=1 $end
$system mwords=1 $end
$guess guess=huckel $end
$basis basnam(1)=metal, ligO,ligO,ligO,ligO,ligO,ligO,
                  ligH,ligH,ligH,ligH,ligH,ligH,
                  ligH,ligH,ligH,ligH,ligH,ligH $end
$data
Cr+3(H2O)6 complex...SBKJC & 6-31G(d) geometry
Th

CHROMIUM   24.0    .0000000000 .0   .0000000000
OXYGEN      8.0    .0000000000 .0  2.0398916104
HYDROGEN    1.0    .7757887450 .0  2.6122732372
$end
!       core potential basis forChromium
$metal
sbkjc

$end
!       normal 6-31G(d) for oxygenligands
$ligO
n31 6
d 1 ; 1 0.8 1.0

$end
!       unpolarized basis forhydrogens
$ligH
n31 6

$end
$ecp
Cr-ecp SBKJC
O-ecp none
   ...snipped...   there must be 6 O's given here
O-ecp none
H-ecp none
   ...snipped...   there must be 12 H's given here
H-ecp none
$end

小菜鸟求关照

brothers 发表于 2015-2-2 16:57
搜索basnam关键词，一定有这个例子 $contrlscftyp=rohf icharg=+3 mult=4 runtyp=gradient
         pp= ...

嗯呢 谢谢

小菜鸟求关照

我看了，还是不理解赝势基组在gamess的定义，比如我从emsl上获得的定义基组写到输入文件里，上面给出的是混合基组的格式，和这个不一样吧，basis basnam(1)=H,Heavy,H,Heavy,H,H,H 这部分的分组名称我也不太明白？
brothers   
看我给你的那个gamess自带的例子当中，
$basis basnam(1)=metal, ligO,ligO,ligO,ligO,ligO,ligO,
                  ligH,ligH,ligH,ligH,ligH,ligH,
                  ligH,ligH,ligH,ligH,ligH,ligH $end
意思很明显啊，对每个原子都需要定义一个代号，然后在$data $end后面，注意是后面，利用basnam定义的这个代号来识别基组。比如用$ligH 基组  空行 $end 来为H原子指明基组内容。这个基组内容可以是gamess自定义的形式，比如n 31G；也可以是从网站上copy下来的一系列本征值那种形式。  在$data $end那部分是不用输入基组内容的，只需要输入对称性和原子坐标。


现在我去试一试。

小菜鸟求关照

$contrl scftyp=rhf runtyp=eda cctyp=CCSD(T) ispher=1 pp=read $end
$system mwords=800 $end
$guess guess=huckel $end
$basis basnam(1)=gas,ligF,ligCl,ligH,ligH $end
$lmoeda matom(1)=3,2 mcharg(1)=0,0 mmult(1)=1,1 $end
$DATA
请问一下我这里哪部分错了 总是出现
ERROR READING INPUT GROUP $BASIS    *****
THE PROBLEM IS WITH THIS INPUT LINE, NEAR THE X MARKER
$BASIS BASNAM(1)=GAS,LIGF,LIGCL,LIGH,LIGH $END                                 
               X
谢谢各位老师啦！

youyno

sobereva 发表于 2015-1-20 19:18
以前随手写的：

!-----LMOEDA能量分解例子，用6-31+G*。必须写nosym=1。用了MP、CC、DFT才有色散能。如 ...

这样的输入直接放在window版本下的gamess-us里面好像不行？是不是和输入格式有关系？

brothers

是没问题的。
你自己好好看看格式要求

zidu113

sobereva 发表于 2015-1-20 19:18
以前随手写的：

!-----LMOEDA能量分解例子，用6-31+G*。必须写nosym=1。用了MP、CC、DFT才有色散能。如 ...

老师您好，请问自定义基组时，使用$basis basnam(1)，关键词basnam后面的括号中的（1）是什么意思？

sobereva

zidu113 发表于 2020-4-29 16:02
老师您好，请问自定义基组时，使用$basis basnam(1)，关键词basnam后面的括号中的（1）是什么意思？

这是GAMESS-US输入文件里数组的定义规则，代表定义第一个元素。但如果定义时后面跟着多个元素值，将被依次认为是从1、2、3...号元素的值。手册里都有说明

zidu113

sobereva 发表于 2020-4-30 07:17
这是GAMESS-US输入文件里数组的定义规则，代表定义第一个元素。但如果定义时后面跟着多个元素值，将被依 ...

老师，谢谢您的及时回复，因为没有找到手册，所以对数组的定义规则还是有些疑问。根据您的提示，我理解basnam(1)中的括号里的值如果是2，比如basnam(2)=metal, ligO, 。。。，那么就是从2个元素开始依次用metal, ligO, 。。。进行赋值，对吗？

sobereva

zidu113 发表于 2020-5-1 13:49
老师，谢谢您的及时回复，因为没有找到手册，所以对数组的定义规则还是有些疑问。根据您的提示，我理解ba ...

手册就在官网上啊，程序包里也都自带了，怎么会找不到
是

zidu113

sobereva 发表于 2020-5-1 19:40
手册就在官网上啊，程序包里也都自带了，怎么会找不到
是

，谢谢您的解答。