应《网络安全法》要求，自2017年10月1日起，未进行实名认证将不得使用互联网跟帖服务。为保障您的帐号能够正常使用，请尽快对帐号进行手机号验证，感谢您的理解与支持！

24小时热门版块排行榜    

小木虫论坛-学术科研互动平台 » 计算模拟区 » 量子化学 » Gaussian » 【求助】基组问题
142/2返回列表上一页12
查看: 2644  |  回复: 13
只看楼主 @他人 存档 新回复提醒 (忽略) 收藏    在APP中查看

qq562676669

金虫 (小有名气)

基组

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

ww1987

铁杆木虫 (小有名气)
★ ★
aylayl08(金币+2):谢谢讨论 2010-04-22 20:47
本人计算的参数为静电势参数 预计6-31G 可以满足！  而对赝势我想 不需要考虑  因为我计算的是C O H三种原子 并不涉及第三周期及以上的元素同时原子个数不多。
赞一下(1人) 回复此楼
12楼2010-04-22 12:44:27
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

xxjiang

银虫 (小有名气)
★
小木虫: 金币+0.5, 给个红包，谢谢回帖
小弟正在学习gaussian，请问各位大仙谁有基组说明书呀，能否发给我一份，xingxing151908@163.com谢谢
赞一下 回复此楼
13楼2012-08-15 19:36:36
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

二沅~

新虫 (小有名气)
★
小木虫: 金币+0.5, 给个红包，谢谢回帖
我算的苯酚结构优化，用的是6-311++G(d,p)，17分钟算完的
赞一下 回复此楼
14楼2015-11-15 22:33:34
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 ww1987 的主题更新
142/2返回列表上一页12
普通表情 龙 兔 虎 猫 高级回复 (可上传附件)
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[基金申请] 能否退出参与的面上项目解除限项 +13 koalala 2026-08-24 15/750 2026-08-24 13:04 by lugia2
[基金申请] 2026国自然函评费到账 +17 羊腰板 2026-08-21 18/900 2026-08-24 12:03 by iaeyuan
[基金申请] 建议基金发布提前给出明确的时间点 +12 kulium 2026-08-21 15/750 2026-08-24 12:02 by iaeyuan
[基金申请] 朋友圈看到的 +7 wangzilk 2026-08-18 9/450 2026-08-24 10:50 by cmrandy
[基金申请] 让我中一个面上吧！ +13 大萍1987 2026-08-20 16/800 2026-08-24 10:23 by 太傻了
[基金申请] 放榜前的不淡定 20+4 snowwithsea 2026-08-19 17/850 2026-08-24 10:20 by echo8914667
[基金申请] 2026年的国家社科基金项目通讯评审的新规则与新动向、新挑战 +4 process2012 2026-08-23 5/250 2026-08-23 19:58 by jurkat.1640
[教师之家] 跳槽后在研项目怎么办？ +5 简单化xn 2026-08-22 10/500 2026-08-23 12:38 by 简单化xn
[基金申请] 只有每年这种时候来逛逛小木虫 +24 yaoyewhu2008 2026-08-20 26/1300 2026-08-22 17:43 by kammury
[基金申请] filecode，4个jtjc了 +13 ziyangfang 2026-08-19 16/800 2026-08-22 17:08 by WH3796
[基金申请] 时间戳今天，20号变了 +5 archvillain 2026-08-20 5/250 2026-08-22 06:12 by hui_daxiao
[基金申请] 看来今天不会放榜了？ +8 chengyan1220 2026-08-21 11/550 2026-08-21 17:52 by dcqxinyang
[基金申请] 时间戳又变了 +13 wuchongjun 2026-08-20 19/950 2026-08-21 17:21 by 紫杉醇
[基金申请] 我面上完蛋了 +7 且听虎啸 2026-08-20 8/400 2026-08-21 12:31 by 酷酷墨镜
[基金申请] 应该是下周三26日公布了吧？ +4 哈哈蛤？ 2026-08-21 4/200 2026-08-21 10:58 by Vivilian
[论文投稿] 投稿咨询 +5 wwm09 2026-08-17 7/350 2026-08-21 10:11 by 期刊论文帮手
[基金申请] 今天放榜没戏了吧 +9 yuleib84 2026-08-19 11/550 2026-08-21 10:06 by gltch
[基金申请] 估计是周四 +3 archvillain 2026-08-18 3/150 2026-08-21 01:48 by jnhyjjm
[基金申请] 重要消息，中午系统在维护 +11 yuleib84 2026-08-18 12/600 2026-08-20 11:09 by xskun
[基金申请] 明天放榜？ +5 Shxjjxjkx 2026-08-18 5/250 2026-08-18 18:14 by -大大大大大-
信息提示
关闭
请填处理意见
关闭