取代基对corrole二阶非线性光学性质的影响

在6-31G(d,p)水平上用B3LYP方法对三苯基corrole(TPC)以及氟苯基corrole系列分子进行几何构型优化.用ZINDO/SOS方法计算了这些分子的一阶超极化率.计算结果表明:TPC分子由于其结构非中心对称使其具有大小为2.46×10-50 C3·m3·J-2的一阶超极化率.引入氟苯基后,氟取代苯基corrole系列分子的一阶超极化率均有显著提高,其中5位及10,15位氟苯基取代的corrole分子具有较大的一阶超极化率(3.71×10-49 C3·m3·J-2).氟苯基corrole系列分子一阶超极化率提高主要起源于分子的二能级分量及三能级分量显著增大,且二能级分量与三能级分量大小相当、符号相同.计算表明,二能级分量的增大主要来源于第三激发态(B1)与基态之间的偶极矩的变化.跃迁通道分析表明三能级分量的增大主要通过改变与结构敏感的激发态跃迁通道,即由TPC分子B带与Q带的耦合变成了氟苯基corrole分子B带内2个不同激发态的耦合.     

戊烯自由基阳离子的密度泛函理论研究

使用密度泛函理论B3LYP方法和6-31G(d,p)、6-31＋G(d,p)、6-311G(d,p)及6-311＋G(d,p)基组，分别对2-C5H10＋和1-C5H10＋的各种构象进行了几何构型优化，并用B3LYP/6-311G(d,p)进行了频率分析计算.计算预言1-C5H10＋具有非平面构型，与以往报导的从头算计算结论相反.在两个自由基阳离子的各种构象的B3LYP几何构型上，进行了B3LYP和UMP2(full)方法的超精细偶合常数计算，得到了比以往更好的结果.     

利用密度泛函理论研究α-联噻吩体系H(C_4H_2S)_nH的结构和电子光谱

应用密度泛函理论,在B3LYP/6-31G(d)和CAM-B3LYP/6-31G(d)水平上优化了α-联噻吩体系H(C4H2S)n H(n=2¹3)的基态几何构型;与此同时,利用TD-B3LYP方法计算了H(C4H2S)n H的吸收光谱,得到了其垂直激发能和体系大小n的解析表达式;并采用TD-CAM-B3LYP方法研究了其发射光谱.研究表明,H(C4H2S)n H的基态结构呈现三种构型:螺旋上升型、环型和稍有弯曲的带状结构;其中前两种构型是顺式结构,最后一种构型是反式结构;环型结构具有C2对称性,其他两种结构具有C1对称性.此外,其吸收光谱和发射光谱计算值与实验值吻合.     

双环金属铱(Ⅲ)异腈配合物的二阶非线性光学性质

采用密度泛函理论(DFT)方法对双环金属Ir(III)异腈配合物的非线性光学(NLO)性质进行计算研究。用B3PW91(UB3PW91)(金属原子采用LANL2DZ基组,非金属原子采用6-31G*基组)方法对配合物进行几何结构优化。在优化构型基础上,采用B3PW91(UB3PW91)和B3LYP(UB3LYP)方法计算了配合物的第一超极化率(βtot),并用CAM-B3LYP(UCAM-B3LYP)(金属原子采用LANL2DZ基组,非金属原子采用6-31G**基组)方法模拟配合物的吸收光谱。结果表明,主配体的取代基(R1)和副配体的取代基(R2)对第一超极化率值贡献不大。配合物发生氧化还原反应,电荷转移方式增多,电荷转移程度增大,使βtot值显著增加,其中1a+([(C∧N)2Ir(CNR)2]+(R=CH3))发生氧化反应和还原反应的βtot值分别增大了75倍和144倍。因此,这类双环金属铱(III)异腈配合物的氧化还原反应可以有效地调节其二阶NLO性质。     

端基取代的长链硅烷二阶超极化率的量子化学研究

对端基取代的一维无限长反位硅烷H₂N-(SiH₂-SiH₂)_n-NO₂的二阶超极化率进行了系统的量子化学研究. 通过仔细检验和选择外场强度, 采用9个外场强度(0.0000, ±0.0008, ±0.0012, ±0.0016, ±0.0020 a.u.)计算的体系能量来确定4阶场强展开式中的5个系数, 从而得到可靠的二阶超极化率. 建议数据拟合时用二阶超极化率单元值的平均值形式γ(n)/n作为拟合对象, 同时用1/n的2阶多项式作为拟合函数, 以得到无限长链的二阶超极化率极限值. 拟合数据范围的选择应该使该数据范围得到的极限值与其临近数据范围得到的极限值的均方偏差最小. 分子构型的优化使计算的二阶超极化率增加大约20%, 在基组中增加极化函数使二阶超极化率在无限长链时的极限值减少大约15%. 相关效应的影响最大, MP2的结果比RHF的结果增加近一倍. 根据本文最高水平MP2/6-31G(d)//RHF/6-31G的计算, 端基取代的一维无限长反位硅烷H₂N-(SiH₂-SiH₂)_n-NO₂的二阶超极化率的每单元极限值为0.8364×10⁶ a.u.     

供体-锯齿型碳纳米管-受体的非线性光学响应对管径的依赖性

用密度泛函 B3LYP 方法得到一系列不同管径的氨基-硝基双取代纳米管NH₂-(n,0)CNT-NO₂ (n＝5～10)的几何结构. 用CAM-B3LYP 方法计算了体系的第一超极化率(β₀). 研究表明, β₀ 对管径大小有着极强的依赖性. 特别是对偶数纳米管(n＝6, 8 和10), 其β₀ (5.5×10³～9.8×10⁴ au)是奇数纳米管β₀ (1.0×10³～2.8×10³ au)的5～35 倍.     

间位取代卟啉及其锌卟啉的电子吸收光谱

取代的卟啉类衍生物在气敏传感器方面具有广泛的应用前景.本文采用了密度泛函理论(DFT)和含时密度泛函理论(TD-DFT)研究了四种不同取代基的卟啉衍生物(meso位四硝基苯基/四氨基苯基卟啉(NO2PP,NH2PP)及其相应的锌金属卟啉衍生物(NO2ZnPP,NH2ZnPP))的紫外和近紫外光谱特征.利用两种不同的交换相关泛函(广义梯度近似泛函(PBE)和杂化密度泛函(B3LYP))优化了上述四种物质的结构,并应用TD-DFT计算了相应的电子激发能量和振动强度.结果表明,取代卟啉的吸收光谱与大量的电子跃迁有关;与B3LYP泛函预测的光谱相比,PBE泛函所得B带以及Q带的波长位置与实验值更为接近.另外,计算所得硝基取代基卟啉的B带相对于氨基取代基卟啉的B带发生了红移,这与实验现象也保持一致.由于卟啉衍生物的三重激发态在电子转移中有很重要的应用,因此在PBE/6-31G(d)水平上计算了四种物质的最低三重激发态能量,分别为1.426、1.469、1.608和1.581eV.     

1,3,5-三硝基苯在高岭石表面吸附的理论研究

构造了高岭石硅氧层和铝氧层的团簇模型(分别为Si₁₃O₃₇H₂₂和Al₆O₂₄H₃₀), 并分别在B3LYP/6-31G(d), MP2/6-31G(d)//B3LYP/6-31G(d)和B3LYP/6-311++G(d,p)//B3LYP/6-31G(d)理论水平上对1,3,5-三硝基苯(TNB)在高岭石表面的吸附性质(如优化的几何构型、 结构参数、 吸附能、 振动频率、 静电势和分子轨道等)进行了研究. 结果表明, TNB和硅氧层表面间的相互作用以静电和范德华相互作用为主; TNB与铝氧层间的相互作用以氢键为主, 且TNB和铝氧层间相互作用的能量更低, 结构更稳定.     

C_3O_2分子结构和光谱的密度泛函理论研究

使用密度泛函理论,在B3LYP/6-31G(d)和B3LYP/6-311G(2d)水平上,研究了C_3O_2分子的可能几何构型,并在6-31G(d)水平上计算了其中2种总能量最小的构型的振动频率,同时与实验观察值进行了比较, 计算结果当C_3O_2分子具有C2v对称性的W型弯曲结构（键角C－C－C和C－C－O分别为162.3°和178.8°）时,振动频率的计算值和实验观察值非常吻合。     

双噻唑苯二聚体自由基分子的极化率与二阶超极化率的理论研究

采用量子化学UMP2/6-31G(d,p)方法优化双噻唑苯二聚体自由基分子的几何结构,以0.05nm为单位步长拉长与缩短2分子片之间的距离,选取5个点,采用DFTUB3LYP/6-31G(d,p)方法,对双噻唑苯二聚体自由基分子的极化率和二阶超极化率进行理论计算.结果表明,自由基体系的单重态为相对稳定状态.在完全重叠的体系中,在单、三重态时极化率都随着2分子片间距离的增大而增加;三重态时二阶超极化率的绝对值随着2分子片间距离的增大而增大.部分重叠的体系,单重态时极化率随2分子片距离的增大而减小;三重态时,二阶超极化率的绝对值随着2分子片间距离的增大而增大.     

OH+C2H2与OH+C2H4加成反应的热力学和动力学

OH基与乙炔(乙烯)的反应是控制大气中OH基浓度的重要化学反应。对反应OH+C_2H_2,1975年Davis等用FP-RF技术测定了反应的速率常数。1977年,Perry又用同样方法对该反应作出了研究;结果表明,此反应的速率常数强烈地依赖于压力,与Davis等人的实验结果不符。对反应OH+C_2H_4,Atkinson等人的实验研究表明,压力在30.0kPa以下,反应的速率常数随压力而改变;而在30.0—88.4 kPa之间,与压力无关.这与前人的结果不同。两个反应的产物也因温度不同而异.此外,上述反应一般是在近101kPa下发生的,而实验结果大多是在低压下获得的。且目前尚无直接的理论计算结果。为此,我们在从头算水平上用统计热力学方法,对反应     

Theoretical studies on vibrational spectra of mixed cyanide-halide complexes of platinum(IV) and palladium(IV)

The vibrational spectra of mixed cyanide-halide complexes, M(CN)4X 2 2- and M(CN)5X2- (M=Pt and Pd; X=F, Cl, Br and I), have been systematically investigated by ab initio RHF, B3LYP and MP2 methods with LanL2DZ and SDD basis sets. The calculated vibrational frequencies of platinum complexes are evaluated via comparison with the experimental values. In the infrared frequency region, the C--N stretching vibrational frequencies calculated at B3LYP level with two basis sets are in good agreement with the observed values with deviations, -16-4 cm(-1) for Pt(CN)4X 2 2- and -18 to -2 cm(-1) for Pt(CN)5X2-. However, in far-infrared region, the results obtained at RHF level are better than those calculated at B3LYP and MP2 levels. For RHF/SDD method, the deviations for Ptz.sbnd;X and Ptz.sbnd;C stretching vibrational frequencies are -14-1 and -12 to -2 cm(-1) in the complex Pt(CN)4X2 2-, -19 to -11 and -15-14 cm(-1) in the Pt(CN)5X2- complex, respectively. The vibrational frequencies of palladium(IV) and some platinum(IV) complexes that have not been experimentally reported are predicted.     

CH3SH…HOO开壳型氢键复合物的结构及其电子密度拓扑性质

在DFT-B3LYP/6-311++G**水平下求得CH3SH…HOO复合物势能面上的稳定构型. 计算结果表明, 在HOO以其O8—H7作为质子供体与CH3SH分子中的S5原子为质子受体形成的氢键复合物1和2中, O8—H7明显被“拉长”, 且其伸缩振动频率发生显著的红移, 红移值分别为330.1和320.4 cm-1; 在CH3SH分子以其S5—H6作为质子供体与HOO的端基O9原子为质子受体形成的氢键复合物3和4中, 也存在类似的情况, 但S5—H6伸缩振动频率红移不大. 经MP2/6-311++G**水平计算的4种复合物含BSSE校正的相互作用能分别为-20.81, -20.10, -4.46和-4.52 kJ/mol. 自然键轨道理论(NBO)分析表明, 在CH3SH…HOO复合物1和2中, 引起H7—O8键长增加的因素包括两种电荷转移, 即孤对电子n1(S5)→σ*(H7—O8)和孤对电子n2(S5)→σ*(H7—O8), 其中后者为主要作用. 在复合物3和4中也有相似的电荷转移情况, 但轨道间的相互作用要弱一些. AIM理论分析结果表明, 4个复合物中的S5…H7间和O9…H6间都存在键鞍点, 且其Laplacian量▽2ρ(r)都是很小的正值, 说明这种相互作用介于共价键和离子键之间, 偏静电作用为主.     

A Raman microprobe study of melt inclusions in kimberlites from Siberia, Canada, SW Greenland and South Africa

The geometries, energies and vibrational frequencies of various polyborates in both gaseous and aqueous phase were calculated at the B3LYP/aug-cc-pVDZ level. The calculated total symmetrical stretching Raman shifts of B(OH)(3), B(OH)(4)(-), B(2)O(OH)(4), B(2)O(OH)(5)(-), B(2)O(OH)(6)(2-), B(3)O(3)(OH)(3), B(3)O(3)(OH)(4)(-), B(3)O(3)(OH)(5)(2-), B(3)O(3)(OH)(6)(3-), B(4)O(5)(OH)(4)(2-) and B(5)O(6)(OH)(4)(-) were assigned to 877.40, 735.33, 785.22, 792.90, 696.79, 587.72, 599.06, 740.16, 705.01, 551.67 and 521.04cm(-1), respectively. The results can be used as the characteristic frequency for polyborates in aqueous phase at room temperature. At least six types of polyborates B(OH)(3), B(OH)(4)(-), B(3)O(3)(OH)(4)(-), B(3)O(3)(OH)(5)(2-), B(4)O(5)(OH)(4)(2-) and B(5)O(6)(OH)(4)(-), occur in aqueous solutions at ambient temperature. The chemical species distribution and the relevant interaction mechanisms among polyborates in the solutions were also suggested.     

Thermodynamic properties of molecular borane amines and the [BH4-][NH4+] salt for chemical hydrogen storage systems from ab initio electronic structure theory

The heats of formation for the borane amines BH3NH3, BH2NH2, and HBNH, tetrahedral BH4-, and the BN molecule have been calculated by using ab initio molecular orbital theory. Coupled cluster calculations with single and double excitations and perturbative triples (CCSD(T)) were employed for the total valence electronic energies. Correlation consistent basis sets were used, up through the augmented quadruple-zeta, to extrapolate to the complete basis set limit. Core/valence, scalar relativistic, and spin-orbit corrections were included in an additive fashion to predict the atomization energies. Geometries were calculated at the CCSD(T) level up through at least aug-cc-pVTZ and frequencies were calculated at the CCSD(T)/aug-cc-pVDZ level. The heats of formation (in kcal/mol) at 0 K in the gas phase are Delta Hf(BH3NH3) = -9.1, Delta Hf(BH2NH2) = -15.9, Delta Hf(BHNH) = 13.6, Delta Hf(BN) = 146.4, and Delta Hf(BH4-) = -11.6. The reported experimental value for Delta Hf(BN) is clearly in error. The heat of formation of the salt [BH4-][NH4+](s) has been estimated by using an empirical expression for the lattice energy and the calculated heats of formation of the two component ions. The calculations show that both BH3NH3(g) and [BH4-][NH4+](s) can serve as good hydrogen storage systems which release H2 in a slightly exothermic process. The hydride affinity of BH3 is calculated to be 72.2 kcal/mol, in excellent agreement with the experimental value at 298 K of 74.2 +/- 2.8 kcal/mol.     

Ab initio/GIAO-CCSD(T) study of structures, energies, and 13C NMR chemical shifts of C4H7(+) and C5H9(+) ions: relative stability and dynamic aspects of the cyclopropylcarbinyl vs bicyclobutonium ions

The structures and energies of the carbocations C 4H 7 (+) and C 5H 9 (+) were calculated using the ab initio method. The (13)C NMR chemical shifts of the carbocations were calculated using the GIAO-CCSD(T) method. The pisigma-delocalized bisected cyclopropylcarbinyl cation, 1 and nonclassical bicyclobutonium ion, 2 were found to be the minima for C 4H 7 (+) at the MP2/cc-pVTZ level. At the MP4(SDTQ)/cc-pVTZ//MP2/cc-pVTZ + ZPE level the structure 2 is 0.4 kcal/mol more stable than the structure 1. The (13)C NMR chemical shifts of 1 and 2 were calculated by the GIAO-CCSD(T) method. Based on relative energies and (13)C NMR chemical shift calculations, an equilibrium involving the 1 and 2 in superacid solutions is most likely responsible for the experimentally observed (13)C NMR chemical shifts, with the latter as the predominant equilibrating species. The alpha-methylcyclopropylcarbinyl cation, 4, and nonclassical bicyclobutonium ion, 5, were found to be the minima for C 5H 9 (+) at the MP2/cc-pVTZ level. At the MP4(SDTQ)/cc-pVTZ//MP2/cc-pVTZ + ZPE level ion 5 is 5.9 kcal/mol more stable than the structure 4. The calculated (13)C NMR chemical shifts of 5 agree rather well with the experimental values of C 5H 9 (+).     

InnNa和InnNa+(n=2-8)的团簇结构和电子性质

在密度泛函理论B3LYP水平上, 对InnNa和InnNa+(n=2-8)团簇进行了结构优化和振动频率计算. 计算结果表明, InnNa(n=2、3、4、6)最稳定结构中的对称性分别为C2v、C3v、C4v和C2v, 而InnNa(n=5、7、8)的最稳定结构的对称性为C1点群. 从InnNa(n=4-8)的最稳定结构可以看出, Na原子均位于四个In原子形成的四边形面上. 对于InnNa+(n=2-8), 除了In2Na+、In4Na+和In7Na+, 其它结构都与其中性结构相似. 进一步计算InnNa(n=2-8)团簇的平均结合能、能量的二阶差分以及绝热电离能表明, InnNa(n=2-8)团簇能量的二阶差分呈现奇偶交替特征, In4Na和In6Na较其它团簇更为稳定, 而且理论计算得到的绝热电离能和实验结果吻合得很好.     

Molecular structure of the azidoalane [Me2N(CH2)3]AltBu(N3). X-ray structural determination and ab initio calculations

A single crystal of the azidoalane [Me₂N(CH₂)₃]AltBu(N₃) (1a), grown in a capillary using a miniature zone melting procedure, was investigated by X-ray analysis. Compound 1a (C₉H₂₁AlN₄) is a monomeric species in the solid state, which crystallizes in the monoclinic space group P2₁ with a=6.8560(14) Å, b=12.251(3) Å, c=7.786(2) Å, β=108.51(3)° and Z=2. The results of the X-ray structural determination are compared with the calculated structure of 1a (HF/6-31G(d) and B3LYP/6-31G(d) level of theory). Whereas the overall agreement between the measured and calculated structure is good, the Al–N donor-bond length differs by 11 and 12 pm at the HF and B3LYP level, respectively. To evaluate the effects of a polar environment on the molecular structure of 1a self-consistent reaction field (SCRF) calculations at the HF and B3LYP level with the 6-31G(d) basis set were performed.     

SiH2自由基与HNCO反应机理的理论研究

采用密度泛函理论B3LYP方法研究了SiH₂自由基与HNCO的反应机理, 并在B3LYP/6-311＋＋G**水平上对反应物、中间体、过渡态进行了全几何参数优化, 通过频率分析和内禀反应坐标(IRC)确定了中间体和过渡态. 为了得到更精确的能量值, 又用QCISD(T)/6-311＋＋G**方法计算了在B3LYP/6-311＋＋G**水平优化后的各个驻点的相对能量. 计算结果表明SiH₂自由基与HNCO的反应有五条反应通道, 其中顺式反应通道SiH₂＋HNCO→IM3→ TS4→IM5→TS5→IM6→SiH₂NH＋CO反应能垒最低, 为主反应通道.     

Oxo-hydroxy tautomerism of 5-fluorouracil: water-assisted proton transfer

Post-Hartree-Fock ab initio quantum chemical calculations were performed for 5-fluorouracil in the gas phase and in a three-water cluster. Full geometry optimizations of the 5-fluorouracil-water complexes were carried out at the MP2/6-31+G(d,p) level of theory. MP4/6-31+G(d,p)//MP2/6-31+G(d,p) and MP4/6-31++G(d,p)//MP2/6-31+G(d,p) single-point calculations were performed to obtain more accurate energies. In water solution, 5-fluorouracil exists mainly in the 2,4-dioxo form (A). We propose that the populations of the 2-hydroxy-4-oxo (B) and 4-hydroxy-2-oxo (D) tautomers are 1 x 10(-4)% and 3.9 x 10(-8)%, respectively, on the basis of the relative stabilities of the tautomers calculated at the MP4/6-31++G(d,p)//MP2/6-31+G(d,p) level of theory. A profound difference between isolated and hydrated 5-fluorouracil is noted for the height of the tautomerization barrier. In the absence of water, the process of proton transfer is very slow. The addition of water molecules decreases the barrier by 2.3 times, making the process much faster. The minimum energy path (MP2/6-31+G(d,p)) for water-assisted proton transfer in trihydrated 5-fluorouracil was followed. CNDO/S-CI calculations predict singlet pi-pi(*) electron transitions at 312 nm for B and at 318 nm for D. The fluorescence spectrum of 5-fluorouracil in water confirms the presence of the hydroxy tautomer.