|
|
基组
★ ★ ★ ★ ★
小木虫(金币+0.5):给个红包，谢谢回帖交流
zhou2009(金币+2): 2010-04-22 12:36
aylayl08(金币+2):谢谢讨论 2010-04-22 20:47
基组用的越精确，计算的结果也就越精确，相应的时间当然也长，个人建议根据自己的电脑配置，能算多精确就算多精确，误差越小当然越好。
1..基组 基组有两种，一种是全电子基组，另一种是价电子基组。价电子基组对内层的电子用包含了相对论效应的赝势（缩写PP，也称为模型势，有效核势，缩写ECP）进行近似。以铟原子为例，可以输入下面的命令来观察铟原子3-21G基组的形式：
# 3-21G GFPrint
Title
0 2
In
输出结果为：
3-21G (6D, 7F)
S 3 1.00 N=1层，s轨道
.1221454700E+05 .6124760000E-01
.1848913600E+04 .3676754000E+00
.4063683300E+03 .6901359000E+00
SP 3 1.00 N=2层，s和p轨道
.5504422550E+03 -.1127094330E+00 .1523702990E+00
.1197743540E+03 .8344349530E-01 .6096507530E+00
.3866926950E+02 .9696880360E+00 .3970249500E+00
D 3 1.00 N=3层，d轨道
.1021735600E+03 .1205559000E+00
.2839463200E+02 .4884976000E+00
.8924804500E+01 .5850190000E+00
SP 3 1.00 N=4层，s和p轨道
.6572360380E+01 .4284830560E+00 .1091304620E-01
.2502157560E+01 -.4633643610E+00 .5036758870E+00
.9420245940E+00 -.8219679320E+00 .5581808650E+00
D 3 1.00 N=4层，d轨道
.4535363700E+01 .2508574000E+00
.1537148100E+01 .5693113000E+00
.4994922600E+00 .3840635000E+00
以上都是内层电子轨道，每个轨道用三个高斯型函数拟合，也就是3-21G中的3。
SP 2 1.00 N=5层，s和p轨道
.1001221380E+01 -.4364171950E+00 -.2316333510E-01
.1659704190E+00 .1189893480E+01 -.9903308880E+00
这是价电子轨道，每个轨道用二个高斯型函数拟合，也就是3-21G中的2
SP 1 1.00 N=6层，s和p轨道
.5433974090E-01 .1000000000E+01 .1000000000E+01
这是空轨道，每个轨道用一个高斯型函数拟合，也就是3-21G中的1。可见3-21G是全电子基组.
再看看铟原子LanL2DZ基组的结果：
LANL2DZ (5D, 7F)
S 2 1.00
.4915000000E+00 -.4241868100E+01
.3404000000E+00 .4478482600E+01
S 1 1.00
.7740000000E-01 .1000000000E+01
P 2 1.00
.9755000000E+00 -.1226473000E+00
.1550000000E+00 .1041757100E+01
P 1 1.00
.4740000000E-01 .1000000000E+01
这些都是价电子轨道。在ECP上加上了双zeta（DZ）基组。所以LanL2DZ是价电子基组。
区分核势和价电子基组的方法是：
核势第一行有一个被核势代替的内壳层电子个数，而且就是固定的那么几个偶数。像上面的例子，核势代替In的46个内壳层电子。
如何得到Gaussian中ECP基组的赝势呢？只要在Route Section行添加IOp(3/18=1)就行了。例如要输出Cl原子完整的ECP基组CEP-121G，就必须用到GFPrint和IOp(3/18=1)关键字，像下面这样：
#p HF/CEP-121G GFPrint IOp(3/18=1)
Title Card Required
0 2
Cl
其中GFPrint输出CEP-121G的价电子基组部分,这是输出结果中Cl原子的CEP-121G基组部分，输出结果中第一部分是GFPrint结果，第二部分是IOp(3/18=1)的结果：
第一部分：
Standard basis: CEP-121G (6D, 10F)
AO basis set:
第二部分
Pseudopotential Parameters |
|
回复此楼