乙硼烷离子和自由基结构的量子化学研究

用B₃LYP/6-311G(d,p)密度泛函方法对B₂H₅⁺阳离子和B₂H₅^·自由基的几何异构体的空间构型进行了优化,并在此基础上用QCISD(T)/6-311++G(3df,2p)偶合簇法进行了单点能计算和零点能校正.结果表明,B₂H₅⁺单态有2种稳定的几何构型(D_3h,C₁),其中C₁构型是新发现的.B₂H₅⁺三重态阳离子除已知C_s构型外,又发现两种稳定构型(C₁).对于B₂H₅^·自由基体系,共有4种异构体(包括两种新发现的构型C_s),其中,具有单桥结构的C_2v最稳定.用二级多体微扰理论和密度泛函方法对前人所认为稳定的B₂H₅⁺单态的C_2v构型进行了全优化,结果发现该构型始终具有一个虚频,不是稳定构型.对B₂H₅^-阴离子体系的单态和三重态进行的全优化,理论上得出单态时具有C_2v和C_s两种稳定构型,而三重态只有C_2v一种稳定构型.     

(NO)2+复合物阳离子的结构和性质研究

用B3LYP/6-311++G**方法对NO二聚体阳离子(NO)2+进行了研究.几何优化并结合振动分析表明,该复合物存在5种可能的稳定构型.能量最低的是N-N相连的反式异构体,具有C2h对称性.分析了各稳定构型的相对稳定性及成键特征.建立了态态相关并给出异构化过渡态,分析了各构型之间转化的途径.     

原子簇Ni2Fe2P结构稳定性的密度泛函理论研究

利用DFT(密度泛函理论)方法在B3LYP/Lan12dz水平下对原子簇模型Ni₂Fe₂P的二十余种可能构型分别在二、四重态下进行优化计算,分析比较了优化结果的能量、结合能、吉布斯自由能变化和过渡态.结果表明:原子簇Ni₂Fe₂P十种异构体没有虚频,有可能稳定存在于非晶态合金中;其中以具有二重态的构型1的能量最低,结合能、吉布斯自由能变化及过渡态能垒最大,最为稳定;四重态中异构体1',2',3'和4'共存的可能性比较大.     

反式和顺式HOOOH的电子光谱的理论研究

用密度泛函方法(DFT)和全活化空间自洽场方法(CASSCF)以及耦合簇理论(CCSD)优化了反式和顺式HOOOH的平衡几何构型, 用DFT计算了HOOOH顺反异构化反应的势能曲线和谐振动频率. 用含时密度泛函理论(TD-DFT)和二阶全活化空间微扰理论(CASPT2)计算了反式和顺式HOOOH垂直激发能. 计算结果表明: (1)反式异构体比顺式异构体稳定; (2)两种稳定构型的异构化反应有两种路径; (3)对于垂直跃迁能最低的单态和叁态, 反式的垂直跃迁能比顺式的低; (4)在单激发态中, CASPT2方法预测的顺式HOOOH寿命最长的激发态为21A′′, 其跃迁能是167.43 nm, 寿命为 1.44×10−5 s; 反式HOOOH寿命最长的激发态为21A, 其跃迁能是165.52 nm, 寿命为 2.07×10−5 s.     

SiCS分子结构及其稳定性的理论研究

采用DFT,QCISD及CCSD(T)方法对单重态、三重态SiCS的分子体系势能面进行理论计算,在QCISD/6-311G(d)水平上得到由3个过渡态连接的5个稳定构型.经动力学及热力学分析均是稳定的三重态线型分子SiCS(³1)、单重态线型分子SiCS(¹1)以及单重态的环状分子cSiCS(¹2).     

苯在Si(111)-7×7表面化学吸附的理论研究

采用两种大小不同的原子簇模型Si_(30)H_(28)和Si_(13)H_(16)，分别用两层 ONIOM方法（对较大原子簇）和普通量子化学方法（对较小原子族）考察了苯分子 在Si(111)-7 * 7表面的化学吸附。对三种可能的吸附物种分别用DFT或HF方法进行 了计算。通过大小原子簇吸附物种的吸附能以及几何构型优化参数的比较发现，对 于稳定的吸附物种，较小的原子簇基本上可以代替较大的原子簇进行计算，而对于 不太稳定的吸附物种，就不得不考虑周边原子的影响。计算结果表明苯在Si(111)- 7 * 7表面的主要吸附种是双σ成键的1，4加成产物，不稳定的单吸附物种可能是 1，4加成物种的前驱态。     

大环炔基噻吩衍生物的结构和UV-Vis光谱的理论研究

用DFT 和TDDFT方法对大环炔基噻吩衍生物的结构和UV-Vis光谱进行了理论研究. 对分子不同的对称性结构(C₁, C₅, C_5v)进行了优化, 得到了稳定的几何构型.以优化构型为基础计算了分子的UV-Vis光谱, 结果表明, C₅,C_5v对称性下丁基取代的分子构型(C[3T_DA]₅-Bu)都是较稳定的; 当分子构型具有C₅对称性时, 得到的光谱数据与实验值符合的较好. 对于大环噻吩C[3T_DA]₅衍生物, 性质相同取代基的体积大小及分子对称性都将影响结构的稳定性.     

C40H2各种可能异构体稳定性的理论研究

采用AM1和PM3两种半经验方法,对D_5d对称性的C₄₀及C₄₀H₂所有可能异构体的几何构型进行了非限制对称性全优化,得到51种稳定异构体,在此基础上研究了氢的加成反应规律及本体C₄₀和最稳定及最不稳定C₄₀H₂异构体的红外光谱,讨论了影响C₄₀(D_5d)氢加成异构体稳定性及加成位置选择性的三种主要因素:(1)C₄₀本体几何结构;(2)共轭效应;(3)电荷分布影响.     

L-缬氨酸旋光异构的两种光反应可能途径

利用密度泛函理论(DFT)和从头算分子轨道理论对L-缬氨酸的旋光异构光反应机理进行了研究. 分别用B3LYP和MP2方法在6-311++G(d, p)基组级别上全优化得到了S₀和T₁态反应路径上的反应物、产物、中间体以及过渡态结构的几何构型, 给出了反应能垒, 利用含时密度泛函理论(TD-DFT)中的B3LYP/6-311++G(d, p)方法优化得到了S₁态反应路径上的平衡态几何构型. 通过分析反应途径上各个驻点的几何构型特征, 确定了L-缬氨酸在激发态可能通过手性碳上的氢原子以羰基氧或氨基氮为中转媒介发生质子迁移来完成旋光异构反应. 进一步用自洽反应场理论中的极化连续模型(PCM)方法探讨了溶剂化效应对旋光异构反应机理的影响.     

C60亮氨酸衍生物的合成及其理论研究

用1,3-偶极环加成方法合成了含吡咯环C₆₀衍生物C₆₆NH₁₃,并以FTIR、UV-Vis、¹HNMR和LD-TOFMS进行表征.用AM1方法对两种可能的加成产物-[6,6]和[6,5]异构体进行几何构型优化.结果表明,[6,6]异构体更稳定.以优化构型为基础,用INDO/CI方法计算两种加成产物的UV谱,结果表明,[6,6]异构体的特征吸收与实验值一致.本文对这两种异构体的电子跃迁进行理论指认,并分析了光谱红移的原因.     

Excitation of Electronic States of N<Subscript>2</Subscript> in Radio-Frequency Electric Field by Electron Impact

Excitation of electronic states of the N₂ molecule by electron impact is recognized as an essential process in nitrogen plasmas that strongly impacts their chemical reactivity and other properties. Many surface and coating technologies are based on radio-frequency plasma discharges in nitrogen. In this paper the electron impact excitation rate coefficients for singlet and triplet electronic states of the N₂ molecule have been calculated in non-equilibrium conditions in the presence of a radio-frequency electric field. A Monte Carlo simulation has been performed in order to determine non-equilibrium electron energy distribution functions within one period of the electric field. By using these distribution functions, the excitation rate coefficients have been obtained in the frequency range from 13.56 up to 500 MHz, at reduced electric field values from 200 to 700 Td.     

Relative stability of the 3A2, 1A2, and 1A1 states of phenylnitrene: A difference-dedicated configuration interaction calculation

The relative stability of the ³A₂, ¹A₂, and ¹A₁ states of phenylnitrene is evaluated by means of ab initio calculations followed by difference-dedicated configuration interaction (DDCI). This approach is based on effective Hamiltonian theory at a low order of perturbation to select rationally the determinants which contribute to the energy difference. The CI space built on this criterion is then treated variationally. The method allows a considerable reduction of the CI space compared with a complete CAS*SDCI calculation (where CAS stands for complete active space). Depending on the concerned energy difference, different model spaces may be chosen, as illustrated in the ³A₂ → ¹A₂ and the ³A₂ → ¹A₁ transitions in phenylnitrene. Since the CI space may reach considerable dimensions, a direct CI algorithm for selected CI spaces, the SCIEL algorithm, has been used to perform the calculations. The results are in excellent agreement with previous calculations and with available experimental data. © 1996 by John Wiley & Sons, Inc.     

MRD CI study of the photodissociation of HO2 into OH(X2II) + O(3P, 1D)

A MRD CI procedure has been used to calculate several electronic states of the hydroperoxyl radical. The basis set is of double-zeta plus polarization quality augmented with s- and p-type bond and Rydberg functions. The vertical excitation energies of the lowest eight doublet and six quartet states are reported. Oscillator strengths for transitions form the ground to upper doublet states were calculated. A cut of the potential energy surfaces along the OOH fragmentation pathway is used to discuss the mechanisms of HO₂ photodissociation below 6.4 eV. Arguments are presented which indicate O(¹D) rather than O(³P) is the primary dissociation product, and so support the experimental findings rather than theory in the conflict raised earlier on this matter. Ostensibly the dissociation proceeds diabatically on the surface of the initially populated ²A″(1a″ → 2a″) state yielding OH(X²II) + O(¹D).     

Jahn-Teller effect in tetrahedral symmetry: large-amplitude tunneling motion and rovibronic structure of CH4+ and CD4+

The energy level structures of the ground vibronic states of 12CH4+, 13CH4+, and 12CD4+ have been measured by pulsed-field-ionization zero-kinetic-energy photoelectron spectroscopy. The nuclear spin symmetries of the tunneling-rotational levels have been determined in double-resonance experiments via selected rotational levels of the v3=1 and v3=2 vibrational levels of the X 1A1 ground state of CH4. The energy level structures of 12CH4+, 13CH4+, and 12CD4+ have been analyzed with an effective tunneling-rotational Hamiltonian. The analysis together with a group theoretical treatment of the Tx(e+t2) Jahn-Teller effect in the Td(M) group prove that the equilibrium geometry of 12CH4+, 13CH4+, and 12CD4+ has C2v symmetry and characterize the pseudorotational dynamics in these fluxional cations. The tunneling behavior is discussed in terms of the relevant properties of the potential energy surface, some of which have been recalculated at the CCSD(T)/cc-pVTZ level of ab initio theory.     

Theoretical study of the ArH+ photodissociation

The multireference Spin-Orbit (SO) Configuration Interaction (CI) method in its Lambda-S Contracted SO-CI (LSC-SO-CI) version is employed to calculate potential energy curves for the ground and low-lying excited states of the ArH(+) cation. For the first time, electric dipole moments are also computed in the approach, including SO coupling for transitions to the states responsible for the first absorption continuum (A-band) of ArH(+). On this basis, the partial and total absorption spectra in this energy range are obtained. It is shown that absorption in the A-band is dominated by the parallel A(1)Sigma(+)<--X(1)Sigma(+) transition. In the low-energy part of the band (<95 x 10(3) cm(-1)) the absorption is caused by the perpendicular B(1)Pi<--X(1)Sigma(+) excitation, but transitions to the b(3)Pi(0(+),1) states are also not negligible. The branching ratio Gamma for the final photodissociation products is calculated and it is shown to increase smoothly from 0 in the red tail of the band to 1 at E>or= 10(5) cm(-1). The latter value corresponds to the exclusive formation of the spin-excited Ar(+)((2)P(1/2)) ions, and thus leads to the inverse population of the Ar(+)((2)P(1/2)-(2)P(3/2)) ion states.     

The configuration interaction method,and the triplet-singlet splitting in CH2

A configuration-interaction (CI) method in which the interaction matrix is never constructed has been investigated, following the original suggestion of Roos. Two methods have been used (1) for singlet states, which can be represented by a one determinant configuration of doubly occupied orbitals, CI with all singly and doubly excited configurations, (2) for states for which the restricted self-consistent field approximation is a single determinant, CI with all singly and doubly excited determinants. In case (2), the wavefunction may not be exactly an eigenfunction of S ². The methods were investigated using a double-zeta plus polarisation basis for CH₂. Both methods must give the same result for the lowest singlet ground state. Keeping the bond length fixed at 2.10 and 2.04 bohr respectively the bond angle for the singlet and triplet were found to be 100.8 ° and 132.0 °, with energies ?39.0312 a.u. and ?39.0563 a.u. respectively. These are the lowest variational energies obtained for these systems; the singlet-triplet splitting is thus predicted to be 15.4 kcal/mol.     

Electron energy distribution functions and rate and transport coefficients of H2/H/CH4 reactive plasmas for diamond film deposition

Electron energy distribution functions (eedf) and rate and transport coefficients for H₂/H/CH₄ mixtures have been calculated by solving a stationary Boltzmann equation as a function of reduced electric field E/N, of molar fraction, and of different concentrations of electronically excited states. Superelastic electronic collisions superimpose structures to eedf especially for E/N values < 40 Td.     

Optimal Geometrical Configuration of Cobalt Cations in Spinel Oxides to Promote Oxygen Evolution Reaction

MgCo₂O₄, CoCr₂O₄, and Co₂TiO₄ were selected, where only Co³⁺ in the center of octahedron (Oh), Co²⁺ in the center of tetrahedron (Td), and Co²⁺ in the center of Oh, can be active sites for the oxygen evolution reaction (OER). Co³⁺(Oh) sites are the best geometrical configuration for OER. Co²⁺(Oh) sites exhibit better activity than Co²⁺(Td). Calculations demonstrate the conversion of O* into OOH* is the rate‐determining step for Co³⁺(Oh) and Co²⁺(Td). For Co²⁺(Oh), it is thermodynamically favorable for the formation of OOH* but difficult for the desorption of O₂. Co³⁺(Oh) needs to increase the lowest Gibbs free energy over Co²⁺(Oh) and Co²⁺(Td), which contributes to the best activity. The coexistence of Co³⁺(Oh) and Co²⁺(Td) in Co₃O₄ can promote the formation of OOH* and decrease the free‐energy barrier. This work screens out the optimal geometrical configuration of cobalt cations for OER and gives a valuable principle to design efficient electrocatalysts.     

Toward the observation of quartet states of the ozone radical cation: insights from coupled cluster theory

Since the discovery of ozone depletion, the doublet electronic states of the ozone radical cation have received much attention in experimental and theoretical investigations, while the low-lying quartet states have not. In the present research, viable pathways to the quartet states from the lowest three triplet states of ozone, (3)A(2), (3)B(2), and (3)B(1), and excitations from the (2)A(1) and (2)B(2) states of the ozone radical cation have been studied in detail. The potential energy surfaces, structural optimizations, and vibrational frequencies for several states of ozone and its radical cation have been thoroughly investigated using the complete active space self-consistent field, unrestricted coupled cluster theory from a restricted open-shell Hartree-Fock reference including all single and double excitations (UCCSD), UCCSD method with the effects of connected triple excitations included perturbatively, and unrestricted coupled cluster including all single, double, and triple excitations with the effects of connected quadruple excitations included perturbatively. These methods used Dunning's correlation-consistent polarized core-valence basis sets, cc-pCVXZ (X = D, T, Q, and 5). The most feasible pathways (symmetry and spin allowed transitions) to the quartet states are (4)A(1)<--(3)A(2), (4)A(2)<--(3)A(2), (4)A(1)<--(3)B(2), (4)A(2)<--(3)B(1), (4)B(2)<--(3)B(1), (4)A(2)<--(1)A(1), (4)B(2)<--(1)A(1), and (4)A(1)<--(1)A(1) with vertical ionization potentials of 12.46, 12.85, 12.82, 12.46, 12.65, 13.43, 13.93, and 14.90 eV, respectively.     

Ab initio potential energy surfaces, total absorption cross sections, and product quantum state distributions for the low-lying electronic states of N(2)O

Adiabatic potential energy surfaces for the six lowest singlet electronic states of N(2)O (X (1)A('), 2 (1)A('), 3 (1)A('), 1 (1)A("), 2 (1)A(") and 3 (1)A(")) have been computed using an ab initio multireference configuration interaction (MRCI) method and a large orbital basis set (aug-cc-pVQZ). The potential energy surfaces display several symmetry related and some nonsymmetry related conical intersections. Total photodissociation cross sections and product rotational state distributions have been calculated for the first ultraviolet absorption band of the system using the adiabatic ab initio potential energy and transition dipole moment surfaces corresponding to the lowest three excited electronic states. In the Franck-Condon region the potential energy curves corresponding to these three states lie very close in energy and they all contribute to the absorption cross section in the first ultraviolet band. The total angular momentum is treated correctly in both the initial and final states. The total photodissociation spectra and product rotational distributions are determined for N(2)O initially in its ground vibrational state (0,0,0) and in the vibrationally excited (0,1,0) (bending) state. The resulting total absorption spectra are in good quantitative agreement with the experimental results over the region of the first ultraviolet absorption band, from 150 to 220 nm. All of the lowest three electronically excited states [(1)Sigma(-)(1 (1)A(")), (1)Delta(2 (1)A(')), and (1)Delta(2 (1)A("))] have zero transition dipole moments from the ground state [(1)Sigma(+)(1 (1)A('))] in its equilibrium linear configuration. The absorption becomes possible only through the bending motion of the molecule. The (1)Delta(2 (1)A('))<--X (1)Sigma(+)((1)A(')) absorption dominates the absorption cross section with absorption to the other two electronic states contributing to the shape and diffuse structure of the band. It is suggested that absorption to the bound (1)Delta(2 (1)A(")) state makes an important contribution to the experimentally observed diffuse structure in the first ultraviolet absorption band. The predicted product rotational quantum state distribution at 203 nm agrees well with experimental observations.