24小时热门版块排行榜

返回列表

当前只显示满足指定条件的回帖，点击这里查看本话题的所有回帖

ww1987

铁杆木虫 (小有名气)

应助: 0 (幼儿园)
金币: 6781.7
散金: 4
帖子: 168
在线: 310.1小时
虫号: 658113
注册: 2008-11-19
性别: GG
专业: 理论和计算化学

[交流] 【求助】基组问题已有12人参与

请问计算单个正辛醇和对甲基苯酚这类分子用6-31G*基组好还是6-311G**基组好？？？？？
（因为我想计算精度可以的话能否用小点基组可以省时间但有人提出用大点基组）所以非常困扰！如果选取小基组的能否说明原因？

回复此楼

» 猜你喜欢

求助Cu2+1O的CIF文件(PDF: 05-0667) 已经有1人回复
KAUST（阿卜杜拉国王科技大学）MXene 器件方向博士后招聘已经有0人回复
物理化学论文润色/翻译怎么收费? 已经有157人回复
沙特阿拉伯阿卜杜拉国王科技大学（KAUST）电池方向博士后招聘已经有0人回复
福州大学新能源材料与工程研究院招收2026年入学博士已经有0人回复
海南师范大学2026年博士研究生招收 (在职想提升学历人员可报考) 申请考核制已经有0人回复
锂离子电池循环寿命衰减过快，求机理分析与改进思路已经有2人回复
求博导收留已经有5人回复
CSC & MSCA 博洛尼亚大学能源材料课题组博士/博士后招生｜MSCA经费充足、排名优已经有0人回复
德国Karlsruhe Institute of Technology招收电化学储能及联合培养CSC博士已经有1人回复
德国Karlsruhe Institute of Technology招收电化学储能及联合培养CSC博士已经有0人回复

» 本主题相关商家推荐: (我也要在这里推广)

1楼 2010-04-22 01:16:58

已阅回复此楼关注TA 给TA发消息送TA红花 TA的回帖

qq562676669

金虫 (小有名气)

应助: 0 (幼儿园)
金币: 1444.8
散金: 11
红花: 1
帖子: 77
在线: 66.4小时
虫号: 948094
注册: 2010-01-25
专业: 分析仪器与试剂

基组

★ ★ ★ ★ ★
小木虫(金币+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

赞一下(3人)

回复此楼

11楼2010-04-22 11:20:11

已阅回复此楼关注TA 给TA发消息送TA红花 TA的回帖

查看全部 14 个回答

fooo

金虫 (正式写手)

应助: 1 (幼儿园)
金币: 6828
散金: 15
红花: 1
帖子: 318
在线: 872.6小时
虫号: 710228
注册: 2009-02-27
专业: 理论和计算化学
管辖: 量子化学

★ ★ ★
小木虫(金币+0.5):给个红包，谢谢回帖交流
zhou2009(金币+2): 2010-04-22 08:25

一般而言，基组越大越好，但也要根据实际情况而定，既不能一味追求省时间，但如果速度太慢似乎也不可取。这都取决于你自己，

引用回帖:

Originally posted by ww1987 at 2010-04-22 01:16:58:
请问计算单个正辛醇和对甲基苯酚这类分子用6-31G*基组好还是6-311G**基组好？？？？？
（因为我想计算精度可以的话能否用小点基组可以省时间但有人提出用大点基组）所以非常困扰！如果选取小基组的能否 ...

赞一下(2人)

回复此楼

2楼2010-04-22 01:59:32

已阅回复此楼关注TA 给TA发消息送TA红花 TA的回帖