二元碳氮化合物CN_8的理论研究

采用B3LYP,B3PW91和MP2 3种理论方法优化得到了13种稳定的CN8异构体.对比这些异构体的总能发现分别具有六元杂环结构和平面链状结构的2种异构体是最稳定的结构.进一步计算这2种异构体的生成能和解离掉一个N2的过渡态,发现它们都具有与实验合成的CN12相近的热力学稳定性,并且具有六元杂环结构的异构体展现出更好的动力学稳定性.这些研究结果表明具有六元杂环结构的CN8异构体是可能在实验上合成出来的.     

铝原子团簇Al5、Al5-和Al5+稳定结构的密度泛函理论研究

采用密度泛函理论的四种方法:杂化密度泛函B3LYP与B3PW91、Perdew-Wang91交换与相关泛函WP91PW91、局域自旋密度近似SVWN,研究了A15、Al5-和Al5+团簇的多种可能结构,找到了它们稳定的结构与自旋态,与已有的理论结果作了比较,并计算了Al5-的绝热与垂直电子离解能、Al5的绝热与垂直电离势,同有关的实验数据比较,符合较好.同时对四种密度泛函方法的计算结果作了一些比较与讨论.     

C_(24)和B_(12)N_(12)团簇的结构与稳定性

本文采用量子化学密度泛函理论的B3LYP/6-31G*方法,对C24和B12N12团簇的12种异构体进行了优化,并对它们的几何构型、振动频率、核独立化学位移(NICS)和结合能进行了理论探讨, 比较了C24和B12N12团簇结构的稳定性。研究表明:C24团簇的最稳定几何构型为类石墨结构d,B12N12团簇的最稳定结构为4/6笼状结构g。C24异构体的稳定性大小顺序为d > b > f > c > a > e。B12N12团簇异构体稳定性大小顺序为a > f> c> d > e >b。     

膦配体对金团簇[Au@Au8(PR3)8]3+(R=Me,OMe,H,F,Cl,CN)稳定化作用的理论研究

采用ab initio HF, MP2方法和密度泛函理论方法, 对具有D2h和D4d构型的膦配体稳定的过渡金属团簇[Au@Au8(PR3)8]3+(R=Me, OMe, H, F, Cl, CN)进行了几何结构、 电子结构及团簇稳定性等方面的研究. 计算表明, 与D2h构型相比, D4d构型更稳定, 两者能量相差约5~10 kJ/mol. SVWN局域泛函能够对团簇的几何结构给予较准确的描述, MP2方法对团簇的结构参数有所低估, 而离域和杂化泛函则过高地估计了团簇的结构参数. 电子结构分析表明, 中心Au原子与外围的Au原子之间通过 d 电子的成键作用构成团簇内核[Au@Au8]3+, [Au@Au8]3+与PR3配体则通过"σ给予/π反馈"模式成键. PR3配体与[Au@Au8]3+的结合能够加强内核-外围Au原子间的成键作用, 缩小外围Au原子在成键上的差异, 增大前线轨道能级间隙, 从而提高团簇的稳定性. PR3配体中R基团供、 吸电子能力的变化对[Au@Au8(PR3)8]3+结构影响较小, 但对[Au@Au8]3+-PR3结合能影响较大. 能量分析显示, 不同PR3配体与[Au@Au8]3+之间具有相近的轨道作用能, 与R基团供、 吸电子能力相关的非轨道作用能成为影响两者连接牢固程度的决定因素.     

Au12M (M=Na,Mg, Al,Si, P,S, Cl)团簇的结构、稳定性和电子性质

采用基于密度泛函理论的第一性原理方法系统地研究了Au12M(M=Na,Mg,Al,Si,P,S,Cl)团簇的结构、稳定性和电子性质.对团簇的平均结合能、镶嵌能、垂直离化势、最高占据分子轨道(HOMO)和最低未占据分子轨道(LUMO)的能级差、电荷布居分析、自然键轨道(NBO)进行了计算和讨论.对于Au12M(M=Na,Mg,Al)团簇,它们形成了内含M原子的最稳定的笼状结构.然而对于Au12M(M=Si,P,S,Cl)团簇,它们却形成了以M元素为顶点的稳定锥形结构.在这些团簇中发现Au12S团簇相对是最稳定的,这是由于Au12S团簇形成了稳定的满壳层的电子结构.自然电荷布居分析表明:对于所有的Au12M(M=Na,Mg,Al,Si,P,S,Cl)团簇电荷总是从Au原子转向M原子.自然键轨道和HOMO分析表明Au12M团簇中发生了Au原子的s-d轨道和M原子的p轨道间的杂化现象.     

B24N24团簇的结构与稳定性

采用密度泛函理论,在B3LYP/6-31G^*水平下,对B₂₄N₂₄笼状团簇的12种异构体进行了优化,并对它们的几何构型、化学键性质、振动光谱和稳定性进行了探讨.研究表明:具有S₈对称的含有2个八元环、8个四元环和16个六元环的结构h是B₂₄N₂₄笼状团簇最稳定的异构体,只存在B-N键,而无N-N和B-B键.含有五元环结构的稳定性最低.B-B和N-N键对的数目越多,结构的稳定性越低.12种异构体的稳定性顺序为h＞a＞b＞I＞g＞l＞c＞k＞j＞d＞e＞f.     

(Mg3N2)n(n=1～4)团簇结构与性质的密度泛函研究

用密度泛函理论(DFT)的杂化密度泛函B3LYP方法在6-31G*基组水平上对(Mg3N2)n(n=1～4)团簇各种可能的构型进行几何结构优化,预测了各团簇的最稳定结构.并对最稳定结构的振动特性、成键特性、电荷特性和稳定性等进行了理论分析.结果表明:(Mg3N2)n=1～4团簇易形成笼状结构,其最稳定构型中N原子配位数以3、4较多见;团簇主要由Mg-N键组成,Mg-N键长为0.194～0.218nm,Mg-Mg 键长为0.262～0.298 nm;N原子的平均自然电荷为-2.06 e,Mg原子的平均自然电荷为 1.37 e;(Mg3N2)2团簇有相对较高的动力学稳定性.     

Al掺杂Sn12-团簇结构及电子性质的密度泛函理论研究

采用密度泛函理论(DFT)的B3PW91方法,在混合基组水平上对Al掺杂Sn12-团簇几何结构和电子结构进行了计算分析.结果表明,Al内掺杂Sn12-团簇能量更低更稳定,但LU-MO-HOMO能隙较小.外掺杂多面体簇中,电荷从Al原子移向Sn12-笼,趋向形成[Al+Sn122-]结构;内掺杂多面体簇中,电荷从Sn12-笼移向Al原子,趋向形成[Al-@Sn12]结构.     

金基二元合金团簇Au_(12)M(M=Cu,Pt,Ni)催化水煤气变换反应的理论研究

利用密度泛函理论(DFT)对Au12M(M=Cu,Pt,Ni)3种合金团簇的结构稳定性、热力学稳定性和反应活性进行研究,并对金基二元合金团簇催化水煤气变换反应(WGSR)的反应机理进行探讨.研究发现,Au12Ni合金团簇的稳定性及电子活性最优.考察了WGSR在金基二元合金团簇上的氧化还原机理和羧基机理,表明Au12Cu合金团簇上WGSR按照氧化还原机理A进行,水解离后产生的OH*会继续解离为O*和H*(*代表吸附物质);Au12Pt及Au12Ni合金团簇上按照氧化还原机理B进行,2个OH*发生歧化反应.比较3种团簇上的最佳反应路径发现,Au12Cu团簇对WGSR表现出较好的催化活性.     

(BN)_n团簇的结构和稳定性

用HF方法、密度泛函理论的B3LYP以及微扰理论的MP2方法 ,在 6 3 1G(d)基组水平上 ,对 (BN) n(n =1～ 16)团簇的各种可能结构进行了优化 .讨论了环状与笼状稳定团簇的几何构型、自然键轨道 (NBO)、振动频率、结合能、核独立化学位移 (NICS)和能量二次差分 ,得到了 (BN) n(n =1～ 16)团簇结构的稳定性信息 .比较了HF ,B3LYP以及MP2三种理论方法对(BN) n 团簇的适应性所表现出的差异 .     

Density functional theory and the correlation consistent basis sets: the tight d effect on HSO and HOS

The HSO and HOS isomers have been revisited using the DFT functionals, B3LYP, B3PW91, and PBE, in combination with tight d-augmented correlation consistent basis sets, cc-pV(x+d)Z and aug-cc-pV(x+d)Z for second-row atoms. Structures, vibrationally averaged structures, relative energies, harmonic and anharmonic frequencies, enthalpies of formation of HSO and HOS, and the barrier for the HSO/HOS isomerization have been determined. These results were compared with results from previous DFT and ab initio studies in which the standard correlation consistent basis sets were used. The relative energies of the two isomers converge more rapidly and smoothly with respect to increasing basis set size for the tight d-augmented sets than for the standard basis sets. Our best calculations, B3PW91/aug-cc-pV(5+d)Z, for the relative energy of the isomers are in excellent agreement with previous CCSD(T) results given by Wilson and Dunning.     

Ab initio calculation of bowl, cage, and ring isomers of C20 and C20-

High-level ab initio calculations have been carried out to reexamine relative stability of bowl, cage, and ring isomers of C(20) and C(20)(-). The total electronic energies of the three isomers show different energy orderings, strongly depending on the hybrid functionals selected. It is found that among three popular hybrid density-functional (DF) methods B3LYP, B3PW91, PBE1PBE, and a new hybrid-meta-DF method TPSSKCIS, only the PBE1PBE method (with cc-pVTZ basis set) gives qualitatively correct energy ordering as that predicted from ab initio CCSD(T)/cc-pVDZ [CCSD(T)-coupled-cluster method including singles, doubles, and noniterative perturbative triples; cc-pVDZ-correlation consistent polarized valence double zeta] as well as from MP4(SDQ)/cc-pVTZ [MP4-fourth-order Moller-Plesset; cc-pVTZ-correlation consistent polarized valence triple zeta] calculations. Both CCSD(T) and MP4 calculations indicate that the bowl is most likely the global minimum of neutral C(20) isomers, followed by the fullerene cage and ring. For the anionic counterparts, the PBE1PBE calculation also agrees with MP4/cc-pVTZ calculation, both predicting that the bowl is still the lowest-energy structure of C(20)(-) at T=0 K, followed by the ring and the cage. In contrast, both B3LYP/cc-pVTZ and B3PW91/cc-pVTZ calculations predict that the ring is the lowest-energy structure of C(20)(-). Apparently, this good reliability in predicting the energy ordering renders the hybrid PBE method a leading choice for predicting relative stability among large-sized carbon clusters and other carbon nanostructures (e.g., finite-size carbon nanotubes, nano-onions, or nanohorns). The relative stabilities derived from total energy with Gibbs free-energy corrections demonstrate a changing ordering in which ring becomes more favorable for both C(20) and C(20)(-) at high temperatures. Finally, photoelectron spectra (PES) for the anionic C(20)(-) isomers have been computed. With binding energies up to 7 eV, the simulated PES show ample spectral features to distinguish the three competitive C(20)(-) isomers.     

<Emphasis Type="Italic">Ab initio</Emphasis> quantum-chemical calculations of the energies and structures of 1,2-acetylenedithiol isomers

Ab initio quantum-chemical calculations of 1,2-acetylenedithiol isomers were carried out. The MP2(full), DFT(B3PW91, MPW1PW91), G3, G3B3, and CBS-Q methods were used. According to the calculations, the most stable isomers were 1,2-dithiet, thiiranethione, and trans-1,2-dithioglyoxal. The necessity of including basis set functions with a large angular momentum in calculations was confirmed. The relatively high stability of 1,2-dithiet was attributed to the aromaticity of its four-membered ring. It was noted that the carbon-carbon bond in the three-membered rings of the cis- and trans-isomers of thiirenethiols was unusually short.     

Quantum chemical studies on structure, conformation and isomerization of nitroso, nitro substituted benzene and 1,3-cyclopentadiene

The molecular structure, conformational stability and isomerization of nitroso, nitro substituted benzene and 1,3-cyclopentadiene in gas phase have been investigated using ab initio and density functional theory methods. The molecular geometries and energetics of possible conformers were obtained by employing MP2, B3LYP and B3PW91 levels of theory implementing 6-31G* basis set. The relative stabilities of the conformations were evaluated from the energy differences of the structure. Chemical hardness (η) and chemical potential (μ) were calculated at HF/6-31G* level of theory for all the positional and geometrical isomers to study the maximum hardness principle. Each optimized structure has been tested against the imaginary frequencies at MP2/6-31G* level of theory in order to be sure they are located at energy minimum.     

Structure and stability of B8 clusters

We investigated various isomers of B₈ clusters with ab initio (MP2) and density function theory (DFT) methods (B3LYP and B3PW91). Nineteen B₈ isomers were determined to be local minima on their potential energy hypersurfaces by the B3LYP, B3PW91, and MP2 methods. Fifteen of these structures are first reported. The most stable neutral B₈ cluster is the regular heptagon, with another boron atom at the center (D_7h, triplet), in agreement with results reported previously. The natural bond orbital (NBO) analysis and nucleus‐independent chemical shifts (NICS) further reveal that the most stable species have delocalized π bond and multicentered σ bonds and therefore exhibit multiple‐fold aromaticity. © 2004 Wiley Periodicals, Inc. Int J Quantum Chem, 2005     

Structure and energetics of small gold nanoclusters and their positive ions

We performed density functional theory calculations at the PW91PW91/LANL2DZ, PW91PW91/Stuttgart 1997, PW91PW91/CRENBL, B3LYP/LANL2DZ, and SVWN5/LANL2DZ levels of theory to attain the minimum-energy structures of neutral and cationic gold clusters of up to nine gold atoms. We locate the 2D-to-3D (two-dimensional to three-dimensional) transition in cationic clusters as occurring between Au(8) (+) and Au(9) (+). We also demonstrate that we can correlate the 2D-to-3D transition in cationic clusters with a linear extrapolation of the energy differences of the lowest-lying 2D and 3D structures of cluster sizes below the transition. We then use the same approach to predict that the 2D-to-3D transition occurs in neutral clusters at Au(11); this is confirmed by locating 3D Au(11) structures that are lower in energy than the best 2D structures reported to date. We examine the effects of choice of basis sets and exchange-correlation functionals on the relative stabilities and other properties of the calculated structures. Finally we find that there is good agreement between calculated and experimental data for clusters with up to six constituent atoms. For clusters with more than seven atoms, there are significant differences observed between the calculated and experimental properties using SVWN5/LANL2DZ, but there is still good agreement for the other levels of theory used.     

Fourier transforms infrared spectra and structure of triformylmethane. A density functional theoretical study

Molecular structure and vibrational frequencies of triformylmethane have been investigated by means of density functional theory (DFT) calculations. The geometrical parameters and vibrational frequencies obtained in the B3LYP, B3PW91, BLYP, BPW91, G96LYP and G96PW91 levels of DFT and compared with the corresponding parameters of malonaldehyde (MA). Fourier transform infrared spectra of triformylmethane and its deuterated analogue were clearly assigned. Theoretical calculations show that the hydrogen bond strength of triformylmethane is stronger than that of MA, which is in agreement with spectroscopic results.     

Vibrational spectrum and assignments of 2-(4-methoxyphenyl)-1H-benzo[d]imidazole by ab initio Hartree-Fock and density functional methods

The room temperature attenuated total reflection Fourier transform infrared spectrum of the 2-(4-methoxyphenyl)-1H-benzo[d]imidazole has been recorded with diamond/ZnSe prism. The conformational behaviour, structural stability of optimized geometry, frequency and intensity of the vibrational bands of the title compound were investigated by utilizing ab initio calculations with 6-311G** basis set at HF, B3LYP, BLYP, B3PW91 and mPW1PW91 levels. The harmonic vibrational frequencies were calculated and scaled values have been compared with experimental IR spectrum. The observed and the calculated frequencies are found to be in good agreement. The theoretical vibrational spectra of the title compound were interpreted by means of potential energy distributions using VEDA 4 program. Furthermore, the optimal uniform scaling factors calculated for the title compound are 0.9120, 0.9596, 0.9660, 0.9699, and 0.9993 for HF, mPW1PW91, B3PW91, B3LYP and BLYP methods, respectively.     

Pd8簇催化氢化乙炔的反应机理

采用密度泛函理论(DFT)中的B3PW91方法研究了过渡金属Pd₈簇催化氢化乙炔的反应机理. 研究表明: H₂进入Pd₈簇后解离成H原子, 并且只有当H原子吸附在Pd₈簇表面上时, 催化氢化反应才能发生. 过渡金属Pd₈簇催化氢化乙炔的反应机理的研究证实该催化反应从两种反应物出发经过两条不同的反应途径完成催化氢化反应, 两种反应物分别为吸附在Pd₈簇上的乙炔(Pd₈(2H)－CH＝CH)和其同分异构体亚乙烯基吸附物(Pd₈(2H)－C＝CH₂). 两条途径均为多步连续的加氢反应, 不同之处在于从Pd₈(2H)－CH＝CH出发的为单一路径, 解离后的H原子分步依次加成到吸附在Pd₈簇上的乙炔中的C原子上, 直到反应完成生成乙烷. 而从Pd₈(2H)－C＝CH₂出发的路径较为复杂, 分别经过两个不同的过渡态和中间体生成次乙基中间体, 该过程相对应的反应位垒相差约12.552 kJ·mol^-1, 说明这两个过渡态同时存在, 无先后次序. 然后继续加成H原子直到生成乙烷完成反应. 同时, 两条路径分别形成一系列具有应用价值的C₂有机化合物中间体, 其中一些中间体通过分子内质子转移相互转化, 使得原本独立的两条反应路径联系在一起, 成为网状路径.     

Isomers of NCO2: IR-absorption spectra of ONCO in solid Ne

Irradiation of a Ne matrix sample containing NO and CO near 4 K with an ArF excimer laser at 193 nm yielded new lines at 2045.1 and 968.0 cm(-1) that were depleted upon secondary photolysis at 308 nm. These lines are assigned to C=O stretching and mixed stretching modes of ONCO, based on results of 15N-, 13C-, and 18O-isotopic experiments and quantum-chemical calculations. These calculations using density-functional theory (B3LYP and PW91PW91/aug-cc-pVTZ) predict five stable isomers of NCO2: ONCO, NCOO, N-cyc-CO2, CNOO, and cyc-CNOO, listed in order of increasing energy. According to B3LYP calculations, ONCO has a trans configuration, with bond angles of angleONC approximately 136.3 degrees and angleNOC approximately 160.7 degrees. Calculated vibrational wave numbers, IR intensities, 15N-, 13C-, and 18O-isotopic shifts for ONCO agree satisfactorily with experimental results. ONCO was formed from reaction of CO with NO in its excited state.