| 查看: 2631 | 回复: 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中的基组是怎样编写到输入程序中的，请指点。 |
回复此楼» 猜你喜欢
大湾区大学（筹）刘天辉课题组招收博士生（与中山大学物理学院联合培养）
已经有16人回复
-
华南师范大学广州市光谱分析与功能探针重点实验室诚招青年英才和博士后岗位研究人员
已经有11人回复
-
物理化学论文润色/翻译怎么收费?
已经有165人回复
-
华南师范大学广州市光谱分析与功能探针重点实验室诚招青年英才和博士后岗位研究人员
已经有10人回复
-
PbS量子点如何确定浓度
已经有11人回复
-
大湾区大学刘天辉课题组招收2026级博士生（与中山大学物理学院联合培养）
已经有22人回复
-
PbS量子点紫外吸收
已经有0人回复
-
求标准卡片Cu0.4In0.4Zn1.2S2 PDF#32-0340
已经有2人回复
-
比利时鲁汶大学与国家留学基金委共同资助博士研究生CSC-KU Leuven PhD Scholarship
已经有0人回复
-
量子点电致发光器件（QLED）：蓝光量子点（QDs）和纳米氧化锌电子传输材料
已经有0人回复
-
需要一个圆偏振光诱导化学反应，请问使用圆偏振光作为光源，应该需要哪些设备搭建
已经有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