| 查看: 2628 | 回复: 9 | |||
| 当前只显示满足指定条件的回帖，点击这里查看本话题的所有回帖 | |||
[求助]
molpro的Basis library中的基组输入已有2人参与
|
|||
|
在molpro官网的Basis library中有的元素是有A基组的，但在程序安装后的lib文件夹中对应的A基组没有该元素，则运行程序会提示找不到基组而出错。 比如：官网Basis library中有Bi元素的aug-cc-pwCVTZ-DK基组（http://www.molpro.net/info/basis ... -DK&print=1），但在/usr/local/molpro/molprop_2015_1_linux_x86_64_i8/lib中的aug_cc_pwcvtz_dk.libmol中却找不到Bi元素的记录。（所以程序找不到该元素的基组而出错） 问： Basis library中的基组是怎样编写到输入程序中的，请指点。 |
回复此楼» 猜你喜欢
求助ORCA计算单点能的输入输出文件
已经有0人回复
-
中教金源的流动池
已经有0人回复
-
物理化学论文润色/翻译怎么收费?
已经有66人回复
-
RSC ADV状态问题
已经有4人回复
-
湖南大学招收锂离子电池、离子液体、高分子与凝胶材料、电催化方向博士生多名
已经有7人回复
-
南京大学蔡亮课题组诚招2026年申请-考核制博士生2-3名（电解水制氢，XAFS谱学等）
已经有32人回复
-
华北电力大学（北京）（第一性原理计算）博士招生——学博，专博各1人
已经有40人回复
-
求助各位大佬配位化学学习教材！！！
已经有0人回复
-
大连海事大学国家级人才团队2026年博士研究生招生启事
已经有10人回复
-
内蒙古大学青年长江学者王蕾教授课题组2026年博士招聘，2-3人
已经有11人回复
-
多相催化-2026年CSC博士及洪堡合作博士后招生！
已经有0人回复
【答案】应助回帖
★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★
kent1022: 金币+20 2017-07-27 13:52:23
kent1022: 金币+20 2017-07-27 13:52:23
|
可以到 http://tyr0.chem.wsu.edu/~kipeters/basis.html 下载需要的基组。 molpro支持外部导入基组，一般读取基组分两种格式，如果是来自llibrary，那么就是.libmol基组，如果不是放到library里的，那么就是.basis. 以下各给出两例解释基组格式： 1. libmol H s aug-cc-pVQZ-DK AVQZ-DK : 7 5 1.06 4.04 5.05 6.06 7.07 basis set from gbasis database 0.82640000E+02 0.12410000E+02 0.28240000E+01 0.79770000E+00 0.25810000E+00 0.89890000E-01 0.23630000E-01 0.20131000E-02 0.15349900E-01 0.75592000E-01 0.25686570E+00 0.49737150E+00 0.29612770E+00 0.10000000E+01 0.10000000E+01 0.10000000E+01 0.10000000E+01 H p aug-cc-pVQZ-DK AVQZ-DK : 4 4 1.01 2.02 3.03 4.04 basis set from gbasis database 0.22920000E+01 0.83800000E+00 0.29200000E+00 0.84800000E-01 0.10000000E+01 0.10000000E+01 0.10000000E+01 0.10000000E+01 H d aug-cc-pVQZ-DK AVQZ-DK : 3 3 1.01 2.02 3.03 basis set from gbasis database 0.20620000E+01 0.66200000E+00 0.19000000E+00 0.10000000E+01 0.10000000E+01 0.10000000E+01 H f aug-cc-pVQZ-DK AVQZ-DK : 2 2 1.01 2.02 basis set from gbasis database 0.13970000E+01 0.36000000E+00 0.10000000E+01 0.10000000E+01 先是元素名，再是angular momentum，然后是基组reference名字（library 搜索的名字），冒号后面分别是exponents个数，contraction个数，然后是各个contraction的方式。 煮个栗子， H s aug-cc-pVQZ-DK AVQZ-DK : 7 5 1.06 4.04 5.05 6.06 7.07 说明是H s function，然后有7个exponents，5个contracted function，1.06指第一个contraction利用第一到第六个exponent，7.07表示最后一个contraction只用了第七个exponent。 接下来的一串字符 0.82640000E+02 0.12410000E+02 0.28240000E+01 0.79770000E+00 0.25810000E+00 0.89890000E-01 0.23630000E-01 0.20131000E-02 0.15349900E-01 0.75592000E-01 0.25686570E+00 0.49737150E+00 0.29612770E+00 0.10000000E+01 0.10000000E+01 0.10000000E+01 0.10000000E+01 就是首先7个exponents，然后1.06contraction里六个contraction系数，4.04里一个contraction系数（也就是1），5.05里一个，6.06根7.07两个。 2. basis s,H,82.64,12.41,2.824,0.7977,0.2581,0.08989,0.02363 c,1.6,0.0020131,0.0153499,0.075592,0.2568657,0.4973715,0.2961277 c,4.4,1.0 c,5.5,1 c,6.6,1 c,7.7,1 p,H,2.292,0.838,0.292,0.0848 c,1.1,1.0 c,2.2,1.0 c,3.3,1 c,4.4,1 d,H,2.062,0.662,0.19 c,1.1,1.0 c,2.2,1.0 c,3.3,1 f,H,1.397,0.36 c,1.1,1.0 c,2.2,1 简明扼要地多，分别是 angular momentum，element name，然后就是exponents，接下来以c打头的就是各个contraction，相信不难理解。 ------------------------- 第二种格式.basis既可以直接粘到input文件里， basis={ s,H,82.64,12.41,2.824,0.7977,0.2581,0.08989,0.02363 c,1.6,0.0020131,0.0153499,0.075592,0.2568657,0.4973715,0.2961277 c,4.4,1.0 c,5.5,1 c,6.6,1 c,7.7,1 p,H,2.292,0.838,0.292,0.0848 c,1.1,1.0 c,2.2,1.0 c,3.3,1 c,4.4,1 d,H,2.062,0.662,0.19 c,1.1,1.0 c,2.2,1.0 c,3.3,1 f,H,1.397,0.36 c,1.1,1.0 c,2.2,1 } 也可以把这些信息放到文件里来外部引入，例如 basis= { O=avqz-dk include /home/XX/molpro_basis/H_aug_cc_pvqz_dk.basis } 至于怎么在.libmol跟.basis格式之间互相转化，可以自己写script，我写了一个简易的python script可以在molpro，gaussian，cfour，DIRAC等不同软件所需基组格式之间互相转化，有空可以发到github上给大家共享。 至于您的问题，如果是correlation consistent基组的话，一般都可以在上述网站中找到，然后可以自己粘成.basis格式喂给molpro。 |
10楼2017-07-27 08:04:56
6楼2017-05-23 23:42:02
7楼2017-05-24 10:52:05
8楼2017-07-06 00:01:05