收集高斯问题Word下载.docx
- 文档编号:21158489
- 上传时间:2023-01-28
- 格式:DOCX
- 页数:41
- 大小:76.56KB
收集高斯问题Word下载.docx
《收集高斯问题Word下载.docx》由会员分享,可在线阅读,更多相关《收集高斯问题Word下载.docx(41页珍藏版)》请在冰豆网上搜索。
wh1125
我也想弄明白,顺便帮你顶!
gamble22
帮顶啦
:
cool:
天空空
2009-07-30
Bauschlicher'
sANObasissetfor(ScTiVCrMnFeCoNiCu)
youcanfindithere:
EMSLBasisSetExchangehttps:
//bse.pnl.gov/bse/portal
zhouzhongjunzzj
2009-08-01
那里是有基组但是不算不过去啊
2009-08-02
QUOTE:
Originallypostedbyzhouzhongjunzzjat2009/8/116:
21:
那里是有基组但是不算不过去啊
Gaussian的积分程序处理ANO基组效率极低,可以考虑使用MOLCAS,它的SEWARD程序是专门为ANO基组设计,积分效率很高。
基组
基组是量子化学专用语。
量子化学中的基组是在量子化学中用于描述体系波函数的若干具有一定性质的函数,基组是量子化学从头计算的基础,在量子化学中有着非常重要的意义。
基组的概念最早脱胎于原子轨道,随着量子化学的发展,现在量子化学中基组的概念已经大大扩展,不局限于原子轨道的原始概念了。
在量子化学计算中,根据体系的不同,需要选择不同的基组,构成基组的函数越多,基组便越大,对计算的限制就越小,计算的精度也越高,同时计算量也会随基组的增大而剧增。
1.斯莱特型基组
斯莱特型基组就是原子轨道基组,基组由体系中各个原子中的原子轨道波函数组成。
斯莱特型基组是最原始的基组,函数形式有明确的物理意义,但是这一类型的函数,数学性质并不好,在计算多中心双电子积分时,计算量很大,因而随着量子化学理论的发展,斯莱特型基组很快就被淘汰了。
2.高斯型基组
高斯型基组用高斯函数替代了原来的斯莱特函数。
高斯型函数在计算中有较好的性质,可以将三中心和四中心的双电子积分轻易转化为二中心的双电子积分,因而可以在相当程度上简化计算,但是高斯型函数与斯莱特型函数在r=0处的行为差异较大,直接使用高斯型函数构成基组会使得量子化学计算的精度下降。
3.压缩高斯型基组
压缩高斯基组是用压缩高斯型函数构成的量子化学基组。
为了弥补高斯型函数与r=0处行为的巨大差异,量子化学家使用多个高斯型函数进行线性组合,以组合获得的新函数作为基函数参与量子化学计算,这样获得的基组一方面可以较好地模拟原子轨道波函数的形态,另一方面可以利用高斯型函数在数学上的良好性质,简化计算。
压缩高斯型基组是目前应用最多的基组,根据研究体系的不同性质,量子化学家会选择不同形式的的压缩高斯型基组进行计算。
4.最小基组
最小基组又叫STO-3G基组,STO是斯莱特型原子轨道的缩写,3G表示每个斯莱特型原子轨道是由三个高斯型函数线性组合获得。
STO-3G基组是规模最小的压缩高斯型基组。
STO-3G基组用三个高斯型函数的线性组合来描述一个原子轨道,对原子轨道列出HF方程进行自洽场计算,以获得高斯型函数的指数和组合系数。
STO-3G基组规模小,计算精度相对差,但是计算量最小,适合较大分子体系的计算。
5.劈裂价键基组
根据量子化学理论,基组规模越大,量化计算的精度就越高,当基组规模趋于无限大时,量化计算的结果也就逼近真实值,为了提高量子化学计算精度,需要加大基组的规模,即增加基组中基函数的数量,增大基组规模的一个方法是劈裂原子轨道,即使用多于一个基函数来表示一个原子轨道。
劈裂价键基组就是应用上述方法构造的较大型基组,所谓劈裂价键就是将价层电子的原子轨道用两个或以上基函数来表示。
常见的劈裂价键基组有3-21G、4-21G、4-31G、6-31G、6-311G等,在这些表示中前一个数字用来表示构成内层电子原子轨道的高斯型函数数目,“-”以后的数字表示构成价层电子原子轨道的高斯型函数数目。
如6-31G所代表的基组,每个内层电子轨道是由6个高斯型函数线性组合而成,每个价层电子轨道则会被劈裂成两个基函数,分别由3个和1个高斯型函数线性组合而成。
劈裂价键基组能够比STO-3G基组更好地描述体系波函数,同时计算量也比最小基组有显著的上升需要根据研究的体系不同而选择相应的基组进行计算。
6.极化基组
劈裂价键基组对于电子云的变型等性质不能较好地描述,为了解决这一问题,方便强共轭体系的计算,量子化学家在劈裂价键基组的基础上引入新的函数,构成了极化基组。
所谓极化基组就是在劈裂价键基组的基础上添加更高能级原子轨道所对应的基函数,如在第一周期的氢原子上添加p轨道波函数,在第二周期的碳原子上添加d轨道波函数,在过渡金属原子上添加f轨道波函数等等。
这些新引入的基函数虽然经过计算没有电子分布,但是实际上会对内层电子构成影响,因而考虑了极化基函数的极化基组能够比劈裂价键基组更好地描述体系。
极化基组的表示方法基本沿用劈裂价键基组,所不同的是需要在劈裂价键基组符号的后面添加*号以示区别,如6-31G**就是在6-31G基组基础上扩大而形成的极化基组,两个*符号表示基组中不仅对重原子添加了极化基函数,而且对氢等轻原子也添加了极化基函数
7.弥散基组
弥散基组是对劈裂价键基组的另一种扩大。
在高斯函数中,变量α对函数形态有极大的作用,当α的取值很大时,函数图像会向原点附近聚集,而当α取值很小的时候,函数的图像会向着远离原点的方向弥散,这种α很小的高斯函数被称为弥散函数。
所谓弥散基组就是在劈裂价键基组的基础上添加了弥散函数的基组,这样的基组可以用于非键相互作用体系的计算。
8.高角动量基组
高角动量基组是对极化基组的进一步扩大,它在极化基组的基础上进一步添加高能级原子轨道所对应的基函数,这一基组通常用于在电子相关方法中描述电子间相互作用。
在最常用的量化高斯软件当中大多数方法都需要定义基组;
如果在计算执行路径部分没有定义基组,则使用STO-3G基组。
例外的情况是在一些方法中,基组作为该方法的积分部分定义;
这些方法在下面列出:
所有的半经验方法,含计算激发态的ZINDO。
所有的分子力学方法
混合的模型化学:
所有的Gn,CBS以及W1方法。
下面的基组存储在Gaussian03程序的内部(完整的说明见引用的参考文献),按照它们在Gaussian03中的关键字列出(有两个例外):
STO-3G[309,310]
3-21G[311-316]
6-21G[311,312]
4-31G[317-320]
6-31G[317-326]
6-31G+:
Gaussian03还包含GeorgePetersson等人的6-31G+和6-31G++基组,定义为完全基组方[88,327]的一部分。
它们可以用关键字6-31G(d'
)和6-31G(d'
p'
)访问,并可以添加单个或两个弥散函数;
还可以添加f-函数:
例如,6-31H(d'
f),等。
6-311G:
对第一行原子定义6-311G基组,对第二行原子定义MacLean-Chandler的
(12s,9p)—>
(621111,52111)基组[328,329](注意P,S,和Cl的基组是被 MacLean和Chandler称作“负离子”的基组;
它们被认为能比中性原子基组给出更好的结果),Ca和K是Blaudeau等人的基组[322],对第一行过渡元素定义Wachters-Hay[330,331]全电子基组,并使用Raghavachari和Trucks[332]的换算因子,对第三行其它元素使用McGrath,Curtiss等人的6-311G基组[324,333,334]。
注意在对第一行过渡元素使用Wachters-Hay基组时,Raghavachari和Trucks推荐使用换算因子并包含弥散函数;
包含弥散函数需要使用6-311+G的形式。
MC-311G等价于6-311G。
D95V:
Dunning/Huzinaga价电子双zeta基组[335]。
D95:
Dunning/Huzinaga全电子双zeta基组[335]。
SHC:
第一行原子用D95V,第二行原子用Goddard/SmedleyECP[335,336]。
也就是SEC。
CEP-4G:
Stephens/Basch/Krauss的ECP最小基组[337-339]。
CEP-31G:
Stephens/Basch/Krauss的分裂价电子ECP基组[337-339]。
CEP-121G:
Stephens/Basch/Krauss的三重分裂价电子ECP基组[337-339]。
注意超过第二行以后,只有一种CEP基组,所有三个关键字对这些原子是等价的。
LanL2MB:
对第一行原子是STO-3G[309,310],对Na-Bi是LosAlamosECP加上MBS[340-342]。
LanL2DZ:
对第一行原子是D95V[335],对Na-Bi是LosAlamosECP加上DZ[340-342]。
SDD:
对一直到Ar的原子是D95V[335],对周期表其它的原子使用Stuttgart/Dresden
ECP[343-367]。
在Gen的基组输入中,可以用SDD,SHF,SDF,MHF,MDF,MWB的形式指定这些基组/芯势。
注意,芯电子的数量必须按照本章后面的格式指定(例如,MDF28表示替代28个芯电子的MDF势)。
SDDAll:
对Z>
2的原子选用Stuttgart势。
cc-pVDZ,cc-pVTZ,cc-pVQZ,cc-pV5Z,cc-V6Z:
Dunning的相关一致基组[368-372](分别为双-zeta,三-zeta,四-zeta,五-zeta,和六-zeta)。
为了提高计算效率,这些基组删除了多余的函数并进行了旋转[373]。
这些基组在定义中包含了极化函数。
下面的表列出了对于各种原子,基组中包含的价极化函数:
原子cc-pVDZcc-pVTZcc-pVQZcc-pV5Zcc-pV6Z
H2s,1p3s,2p,1d4s,3p,2d,1f5s,4p,3d,2f,1g6s,5p,4d,3f,2g,1h
He2s,1p3s,2p,1d4s,3p,2d,1f5s,4p,3d,2f,1g不可用
B-Ne3s,2p,1d4s,3p,2d,1f5s,4p,3d,2f,1g6s,5p,4d,3f,2g,1h7s,6p,5d,4f,3g,2h,1i
Al-Ar4s,3p,1d5s,4p,2d,1f6s,5p,3d,2f,1g7s,6p,4d,3f,2g,1h不可用
Ga-Kr5s,4p,1d6s,5p,3d,1f不可用不可用不可用
这些基组可以通过给基组关键字添加AUG-前缀(而不是使用+和++符号——见下),用弥散函数增大基组。
但是,He,Mg,Li,Be,和Na元素在这些基组中没有定义弥散函数。
Ahlrichs等人的SV,SVP,TZV和TZVP基组[374,375]。
Truhlar等人的MIDI!
基组[376]。
使用这个基组需要MidiX关键字。
Epr-II和EPRIII:
Barone[377]的基组,对DFT方法(特别是B3LYP)的超精细耦合常数计算进行了优化。
EPR-II是一套极化函数的双-zeta基组,并对s-部分进行了增强:
对H是(6,1)/[4,1],对B到F是(10,5,1)/[6,2,1]。
EPR-III是包含弥散函数(双d-极化函数和一套f-极化函数)的三-zeta基组。
同样,s-部分也进行了改善,以更好地描述芯区域:
对H是(6,2)/[4,2],对B到F是(11,7,2,1)/[7,4,2,1]。
UGBS,UGBS1P,UGBS2P,和UGBS3P:
deCastro,Jorge等人[378-386]的普适Gaussian基组。
后三种关键字形式对标准UGBS基组中的每个函数添加1个,2个或者3个极化函数(也即,UGBS1P对每个s-函数添加一个p-函数,对每个p-函数添加一个d-函数,等等;
UGBS2P对每个s-函数添加一个p-和d-函数,对每个p-函数添加一个d-和f-函 数,UGBS3P对每个s-函数添加一个p-,d-和f-函数,等等)。
Martin和deOliveira的MTSmall,定义为W1方法的一部分(参见W1U关键字)[94]。
DGauss中使用的DGDZVP,DGDZVP2,和DGTZVP基组[387,388]。
优化的构型如何消去虚频
我用gaussian03对一个稳定态构型进行优化时,正常结束,但是查看频率时有三个虚频,-607.-94,-92,如何消去虚频?
谢谢.这个构型还有必要要吗?
csfn
2009-09-27
虚频不小呀,如果结构不太大,可以放弃这个,多建立几个结构测试,注意多重度的影响
要改这个结构,可以看虚频朝哪个方向振动,然后朝该方向调整结构
zeoliters
在高斯的输出文件里,列出了虚频的振动模式(以XYZ三个方向的坐标列出的),你可以把这些坐标同时乘以一个0~1之间的系数(建议使用0.5),再加到优化后的坐标里,再重新优化,一般都消掉。
例如:
1
A'
Frequencies--
-707.3181
Red.masses--
1.0573
Frcconsts
--
0.3117
IRInten
--
24.5978
AtomAN
X
Y
Z
6
0.05
0.02
0.00
2
0.03
0.01
3
-0.21
0.40
4
-0.61
-0.53
5
-0.19
-0.32
对应的虚频-707.3181的振动模式就是
X
-0.21
-0.61
-0.19
将这些坐标乘以一个0~1之间的系数,再加到Standardorientation的坐标里面去,再重新优化就可以了
Standardorientation:
---------------------------------------------------------------------
Center
Atomic
Coordinates(Angstroms)
Number
Type
Z
1
6
0
0.047915
0.719192
0.000000
2
-0.586151
3
-0.880556
-1.164400
4
-0.663020
1.532329
5
0.968595
-1.166173
[Lasteditedbyzeoliterson2009-9-27at12:
49]
只能一个一个消。
如果都想消去,就加那个绝对值最大的那个虚频对应的振动模式,一般后面的都可以自动消去。
自由基正离子的计算
我想算一下自由基正离子的轨道系数,比如苯酚自由基正离子的,我输入的结构就是苯酚结构吗?
用的方法是b3lyp/6-31g**电荷1多重度2,这样算对吗?
请高手多帮忙。
[Lasteditedbymingdongon2009-9-22at10:
24]
2009-09-22
能不能给出结构啊?
好久不用,竟然不知道苯酚自由基正离子是什么东西了,帮你问了几个人,大家都不确定,呵呵
liuyanzhu_79
2009-09-23
就是苯酚结构做质谱时它会失去一个电子带一个正电和一个自由基,我也不知道具体要输入什么结构
fegg7502
Originallypostedbyliuyanzhu_79at2009-9-2317:
10:
就是苯酚结构做质谱时它会失去一个电子带一个正电和一个自由基,我也不知道具体要输入什么结构
把你的结构图给出来
2009-09-24
就是[Ar-OH]+.
不好意思chemdraw画的结构式黏贴不过来
这样啊,你的输入应该是对的
2009-09-25
非常谢谢
自旋态4A2u代表什么意思?
是不是多重度为4?
A2u代表什么意思?
如果要用高斯计算这个自旋态,应该怎么设置?
jghe
2009-09-21
是4重度,A2u更对称性有关
高斯设自旋多重度为4就可以了
zdfbb
4是电子态多重度,A2u是电子态对称性
zhaolixia0814
哦,:
):
)
906980368
2009-10-04
2s+1=4,4表示四重态;
A2u表示电子态,其对称性性是A2u
psfan
2009-10-06
A2u对称性在几何空间上是如何伸展的呢?
怎么用gaussian来计算反应焓变
反应焓变的公式是
dH=dE298+d(PV)
dE298=dEe0+d(dEe)298+dEv0+d(dEv)298+dEr298+dEt298
这其中
1.dEe0是0k时产物与反应物的额能量差,在gaussian输出文件里是指哪一项的值?
2.d(dEe)298指0到298k电子能量变化是不是个固定值?
3.d(dEv)298指0到298k振动能量变化是不是个固定值?
4.dEr298旋转能在gaussian输出文件中是哪一项?
5.dEt298平动能在gaussian输出文件中是哪一项?
如果有其他方法,请指点。
[Lasteditedbybergkamp888818on2009-8-5at10:
53]
ZDBWHZ
2009-08-05
这是理论的原理,在Gaussian中是这样实现的。
先对反应中的每种物质进行结构优化,再算频率,输出文件中就有每种分子的焓值,再相减就得到焓变了。
bergkamp888818
2009-08-06
gaussian只算了内能,没有算焓值吧。
gongxd325
2009-08-07
结构优化之后作振动频率分析,得到SumofelectronicandthermalEnthalpies为焓
即可求焓变
chinanhui
帮楼主顶一下,我是不会
ADF计算过渡金属的过渡态结果怎么样?
相对于高斯,ADF寻找稳定点和过渡态都要快些和容易些.当然了,有老师告诉我ADF有一个问题就是容易收敛,就是说容易收敛到一个虚假的稳定结构上,跟高斯结果相比,高斯更可信.个人觉得虽然有这类问题,使用时注意就行了.因为ADF在做过渡金属体系方面的功能太强了,高斯对于大一些的体系都不可行.我尝试过一个双金属体系,在高斯中花了很大的时间才找到的稳定点,
- 配套讲稿:
如PPT文件的首页显示word图标,表示该PPT已包含配套word讲稿。双击word图标可打开word文档。
- 特殊限制:
部分文档作品中含有的国旗、国徽等图片,仅作为作品整体效果示例展示,禁止商用。设计者仅对作品中独创性部分享有著作权。
- 关 键 词:
- 收集 问题