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

24小时热门版块排行榜    

小木虫论坛-学术科研互动平台 » 计算模拟区 » 量子化学 » Gaussian » 【求助】基组问题
142/2返回列表上一页12
查看: 2321  |  回复: 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
普通表情 龙 兔 虎 猫 高级回复 (可上传附件)
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 085601材料工程专硕求调剂 +5 慕寒mio 2026-03-16 5/250 2026-03-17 21:31 by hmn_wj
[考研] 材料专硕306英一数二 +7 z1z2z3879 2026-03-16 9/450 2026-03-17 17:31 by ccjequ
[考研] 材料工程专硕调剂 +5 204818@lcx 2026-03-17 5/250 2026-03-17 17:27 by Little-xue
[考研] 344求调剂 +4 knight344 2026-03-16 4/200 2026-03-17 17:27 by ruiyingmiao
[考研] 308求调剂 +4 是Lupa啊 2026-03-16 4/200 2026-03-17 17:12 by ruiyingmiao
[考研] 材料专硕326求调剂 +6 墨煜姒莘 2026-03-15 7/350 2026-03-17 17:10 by ruiyingmiao
[考研] 285化工学硕求调剂（081700） +9 柴郡猫_ 2026-03-12 9/450 2026-03-17 10:18 by Sammy2
[考研] 药学383 求调剂 +3 药学chy 2026-03-15 4/200 2026-03-16 20:51 by 元子^0^
[考研] 304求调剂 +4 ahbd 2026-03-14 4/200 2026-03-16 16:48 by 我的船我的海
[考研] 277材料科学与工程080500求调剂 +3 自由煎饼果子 2026-03-16 3/150 2026-03-16 14:10 by 运气yunqi
[考研] 中科大材料专硕319求调剂 +3 孟鑫材料 2026-03-13 3/150 2026-03-14 18:10 by houyaoxu
[基金申请] 有必要更换申报口吗 20+3 fannyamoy 2026-03-11 3/150 2026-03-14 00:52 by zhanghaozhu
[考研] 290求调剂 +9 ADT 2026-03-11 9/450 2026-03-13 21:55 by JourneyLucky
[考研] 0703化学一志愿211 总分320求调剂 +5 玛卡巴卡啊哈 2026-03-11 5/250 2026-03-13 21:40 by JourneyLucky
[考研] 四川大学085601材料工程专硕 初试294求调剂 +4 祝我们好在冬天 2026-03-11 4/200 2026-03-13 21:39 by peike
[考研] 274求调剂 +3 S.H1 2026-03-12 3/150 2026-03-13 15:15 by JourneyLucky
[考研] 277求调剂 +4 anchor17 2026-03-12 4/200 2026-03-13 11:15 by 白夜悠长
[考研] 求调剂 资源与环境 285 +3 未名考生 2026-03-10 3/150 2026-03-13 10:31 by houyaoxu
[考博] 26读博 +4 Rui135246 2026-03-12 10/500 2026-03-13 07:15 by gaobiao
[考研] 081200-11408-276学硕求调剂 +3 崔wj 2026-03-12 4/200 2026-03-12 19:33 by 求调剂zz
信息提示
关闭
请填处理意见
关闭