详细内容

Gaussian 讲解06-基组

大多数方法需要定义基组;如果在执行路径部分没有包含基组,则使用 STO-3G 基组。

例外的情况是在一些方法中,基组为该方法不可分割的一部分;这些方法在下面列出:

u 所有的半经验方法,含用于激发态的 ZINDO

u 所有的分子力学方法。

u 复合的模型化学:所有的 GnCBS 以及 W1 方法。 以下是 Gaussian 09 程序内部存储的基组(完整的说明见引用的参考文献),下面按照 它们对应的 Gaussian 09 关键词列出(有两个例外):

u STO-3G [Hehre69, Collins76]

u 3-21G [Binkley80a, Gordon82, Pietro82, Dobbs86, Dobbs87, Dobbs87a]

u 6-21G [Binkley80a, Gordon82]

u 4-31G [Ditchfield71, Hehre72, Hariharan74, Gordon80]

u 6-31G [Ditchfield71, Hehre72, Hariharan73, Hariharan74, Gordon80, Francl82, Binning90, Blaudeau97, Rassolov98, Rassolov01]

u 6-31GGaussian 09 还包含 George Petersson 及合作者的 6-31G 6-31G基组,定 义为完备基组方法[Petersson88, Petersson91]的一部分。它们可以用关键词 6-31G(d') 6-31G(d',p')调用,并可以添加单个或两个弥散函数;还可以添加 f 函数:例如,6-31G(d'f),等。

u 6-311G:对第一行原子定义 6-311G 基组,对第二行原子定义 MacLean-Chandler 的(12s,9p)à(621111,52111)基组[McLean80, Raghavachari80b] (注意 P,S,和 Cl 的基组是被 MacLean 和 Chandler 称作“负离子”的基组;它们被认为对中性分子能给出更好的结果),Ca 和 K 是 Blaudeau 等人的基组[Blaudeau97],对第一行过渡元素定义 Wachters-Hay [Wachters70, Hay77]全电子基组,并使用 Raghavachari 和 Trucks [Raghavachari89]的换算因子, 对第三行其它元素使用 McGrath,Curtiss 等人的 6-311G 基组[Binning90, McGrath91, Curtiss95]。注意在对第一行过渡元素使用 Wachters-Hay 基组时,Raghavachari 和 Trucks推荐使用换算因子并包含弥散函数;若包含弥散函数需使用 6-311+G 的形式。MC-311G 6-311G 的同义词。

u D95VDunning/Huzinaga 价电子双 zeta 基组[Dunning76]

u D95Dunning/Huzinaga 完全双 zeta 基组[Dunning76]

u SHC:第一行原子用 D95V,第二行原子用 Goddard/Smedley ECP [Dunning76, Rappe81]。也就是 SEC

u CEP-4GStephens/Basch/Krauss 的 ECP 最小基组[Stevens84, Stevens92, Cundari93]

u CEP-31GStephens/Basch/Krauss 的 ECP 分裂价基组[Stevens84, Stevens92, Cundari93]

u CEP-121GStephens/Basch/Krauss 的 ECP 三分裂基组[Stevens84, Stevens92, Cundari93]

注意超过第二行以后,只定义了一种 CEP 基组,所有三个关键词对这些原子都是等价的。

u LanL2MB:对第一行原子是 STO-3G [Hehre69, Collins76],对 Na-La,Hf-Bi 是 Los Alamos ECP 加上 MBS [Hay85, Wadt85, Hay85a]

u LanL2DZ:对第一行原子是 D95V [Dunning76],对 Na-La,Hf-Bi 是 Los Alamos ECP 加上 DZ [Hay85, Wadt85, Hay85a]

u SDDD95V 一直到 Ar 原子[Dunning76],对周期表的其它原子使用 Stuttgart/Dresden ECP [Fuentealba82, Szentpaly82, Fuentealba83, Stoll84, Fuentealba85, Wedig86, Dolg87, Igel-Mann88, Dolg89, Schwerdtfeger89, Dolg89a, Andrae90, Dolg91, Kaupp91, Kuechle91, Dolg92, Bergner93, Dolg93, Haeussermann93, Dolg93a, Kuechle94, Nicklass95, Leininger96, Cao01, Cao02]。在 Gen 的基组输入中,可以用 SDD, SHF, SDF, MHF, MDF, MWB 等种类指定这些基组/芯势。注意,种类之后必须指定芯电子数(例如,MDF28 表示替代 28 个芯电子的 MDF 势)。

u SDDAll:对 Z>2 的原子选用 Stuttgart 势。

u cc-pVDZ, cc-pVTZ, cc-pVQZ, cc-pV5Z, cc-pV6ZDunning 的关联一致基组[Dunning89, Kendall92, Woon93, Peterson94, Wilson96](分别为双-zeta,三-zeta,四-zeta,五-zeta,和 六-zeta)。为了提高计算效率,这些基组删除了冗余函数并进行了旋转[Davidson96]

这些基组在定义中已包含了极化函数。下面的表列出了各种原子的基组中包含的价极化

函数:

image.png

这些基组可以通过给基组关键词添加 AUG-前缀(而不是使用+++符号——见下),用 弥散函数增大基组。

u Ahlrichs 等人的 SV, SVPTZVTZVP [Schaefer92, Schaefer94],和 QZVP [Weigend05]基组。

u Truhlar 等人的 MIDI!基组[Easton96]。使用这个基组需要 MidiX 关键词。

u Epr-II EPR-IIIBarone [Barone96a]的基组,对 DFT 方法(特别是 B3LYP)的超精细耦 合常数计算进行了优化。EPR-II 是具有一套极化函数的双-zeta 基组,并对 s-部分进 行了加强:对 H 是(6,1)/[4,1],对 B 到 F 是(10,5,1)/[6,2,1]。EPR-III 是三-zeta 基组,包括弥散函数,双 d-极化和一套 f-极化函数。同样,s-部分也进行了改善,以更好地描述芯区域:对 H 是(6,2)/[4,2],对 B 到 F 是(11,7,2,1)/[7,4,2,1]。

u UGBSde Castro,Jorge 等人[Silver78, Silver78a, Mohallem86, Mohallem87, daCosta87, daSilva89,

Jorge97, Jorge97a, deCastro98]的通用 Gaussian 基组。给关键词加上后缀,可以添加附加的

极化函数:

UGBSnP|V|O

其中的 n 是整数,表示是否对常规的 UGBS 基组中的每个函数加上 1 个,2 个或者 3 个 极化函数。第二项是一个字母代码,表示给哪个函数加上极化函数:P 表示给所有函数 加上极化函数,V 给所有价函数加上极化函数,O 使用 Gaussian 03 所用的方案(见下)。

例如,UGBS1P 关键词请求对该基组的所有轨道添加一个附加的极化函数,而 UGBS2V 所有价轨道添加两个附加的极化函数。

Gaussian 03 的 UGBSnP 关键词一样,O 后缀添加同样的函数。UGBS1O 对每个 s-函数添加一个 p-函数,对每个 p-函数添加一个 d-函数,等;UGBS2O 对每个 s-函数添加一 个 p-和 d-函数,对每个 p-函数添加一个 d-和 f-函数,等;UGBS3O 对每个 s-函数添加 一个 p-,d-和 f-函数,等。可以照常用+++添加弥散函数;第一个可以指定为 2+,表示对重原子加上两个弥散函 数。

u Martin 和 de Oliveira 的 MTSmall,定义为 W1 方法的一部分(参见 W1U 关键词)[Martin99]

u DGauss 中使用的 DGDZVPDGDZVP2,和 DGTZVP 基组[Godbout92, Sosa92]

u CBSB7:选择 6-311G(2d,d,p)基组,用于 CBS-QB3 高精度能量方法[Montgomery99]。这一名称对第二行原子指定两个附加的 d 极化函数,对第一行原子指定一个 d 函数,对氢指定一个 p 函数(注意,Gaussian 09 不支持这种三段的极化函数语法)。 


最新评论
请先登录才能进行回复登录
技术支持: CLOUD | 管理登录
seo seo