外电场对甲醛分子的结构与电子光谱的影响

采用密度泛函理论(DFT)在B3LYP/6-311++G(d,p)基组水平上,计算了不同外加电场(-8.22×10~9~8.22×10~9 V/m)下甲醛分子基态稳定构型、分子键长、电荷分布、能级分布、能隙、红外光谱、拉曼光谱和分子的总能量.在此基础上利用TDDFT/B3LYP/6-311++G(d,p)方法研究了甲醛分子由基态跃迁到前25个激发态的激发能E、谐振强度f、吸收波长λ受外电场的影响.结果表明:随着C=O连线方向外电场的增加,C=O键键长、氢原子电荷、偶极矩和能隙递增;C—H键键长、C,O原子电荷递减,总能量降低.振动频率与红外强度及拉曼强度由于不同振动有不同变化.甲醛分子UV-Vis光谱随外电场的增加,不同的吸收峰发生了不同程度的蓝移或者红移;外电场对甲醛分子的激发能、谐振强度和吸收波长的强度有一定影响,但随电场变化比较复杂.     

电子转移的外电场效应──苯分子与苯正离子自由基间电子转移的从头算研究

基于半经典电子转移理论,结合量子化学计算,在HF/DZP水平上,研究外电场作用下平行的苯分子-苯正离子自由基体系(C₆H₄)₂⁺的分子内电子转移问题.在给体和受体几何构型优化的基础上,用线性反应坐标确定电子转移过渡态,分别用两态变分方法和基于Koopmans定理的分子轨道跃迁能方法计算电子转移矩阵元V_AB,讨论了V_AB对给体和受体中心距d的指数衰减关系.取中心距为0.6nm,研究了外电场对反应热的影响,计算得到在不同外电场强度下分子内气相电子转移的速率常数k.     

外电场作用下甲基乙烯基硅酮分子结构和电子光谱

采用密度泛函B3P86方法在6-311++G(d, p)基组水平上优化得到了沿分子轴方向不同外电场(0-0.04 a.u.)作用下, 甲基乙烯基硅酮分子的基态电子状态、几何结构、电偶极矩和分子总能量. 在优化构型下利用杂化CIS-DFT方法(CIS-B3P86)研究了同样外电场条件下对甲基乙烯基硅酮的激发能和振子强度的影响. 计算结果表明, 分子几何构型与电场大小呈现强烈的依赖关系, 分子偶极矩μ随电场的增加先减小后急剧增大. 电场为零时, 分子总能量为-483.5532137 a.u., 随着电场增加, 能量升高, 在F=0.02 a.u.时达到最大值-483.5393952 a.u., 此后, 继续增大电场系统总能量则开始降低. 激发能随电场增加急剧减小, 表明在电场作用下, 分子易于激发和离解.     

聚丙烯酸在纳米TiO2表面吸附行为的研究

讨论了聚丙烯酸在纳米TiO₂水悬浮体系中的吸附行为.红外光谱分析和吸附实验结果表明,纳米TiO₂通过氢键吸附PAA.PAA吸附量随着浓度的升高而增大直至饱和吸附量,且分子量越大,饱和吸附量越大.pH值增大,则饱和吸附量减小.在相同条件下,表面吸附层的厚度随PAA分子量、浓度和pH值增大而增大.这是由PAA在颗粒表面构型的变化所致.吸附PAA后的纳米TiO₂的表面电荷密度和ζ电位发生变化,pHiep值向低值方向移动.表面吸附自由能的计算结果说明,PAA在纳米TiO₂表面的吸附是自发过程.     

线性响应和两种态特定方法模拟AF350电子跃迁的溶剂效应

在极化连续模型框架下比较了线性响应与两种不同态特定方法计算的溶液中Alexa Fluor 350(AF350)分子激发能和光谱移动值的差异. AF350的第一激发态S0→S1电子跃迁属于π→π^*跃迁, 主要对应于最高占据分子轨道(HOMO)到最低空轨道(LUMO)的跃迁. 该分子激发态偶极矩大于基态偶极矩, 激发态时溶质溶剂相互作用比基态时更强, 随着溶剂极性增大, 会发生光谱红移的现象. 与实验值相比, 线性响应和两种态特定方法均高估了激发能, 其中以IBSF(Improta-Barone-Scalmani-Frisch)方法得到的激发能最小, 矫正的基态反应场方法(cGSRF)得到的激发能最大. 对于光谱移动值, 3种方法与实验值相比都偏小, 线性响应方法(LR)计算出的误差最大, 而IBSF方法得到的结果与实验值最吻合, 是预测溶液中AF350分子激发能和光谱移动值最准确的方法. 对比了Marcus传统理论和基于约束平衡的非平衡溶剂化理论的结果, 发现后者得到的激发能和光谱移动值更接近于实验值.     

9,10-二氰基蒽和杜烯间光诱导电子转移的电荷分离态和电子耦合

采用从头算方法,讨论了9,10-二氰基蒽(DCA)和杜烯(DUR)间光诱导电子转移反应的态-态跃迁.考虑基组重叠误差(BSSE)对相互作用能的校正,用MP₂方法优化得到重叠式[DCA…DUR]配合物的稳定构型.用单激发组态相互作用(CIS)方法讨论了[DCA…DUR]配合物的光诱导电荷分离和电荷复合过程.根据广义Mulliken-Hush(GMH)模型,计算了电荷复合过程的电子耦合矩阵元.结果表明,[DCA…DUR]配合物的S₀→S₁和S₀→S₂跃迁产生了两个强的局域激发态,S₀→S₃跃迁直接导致电荷分离态,小的振子强度预测该电荷转移(CT)跃迁是一弱跃迁,电荷分离态S₃衰变到低局域激发态或基态的电荷复合是可能的.     

4-氨基-1,8-萘二酰亚胺的结构及光谱性质的理论研究

采用HF和密度泛函理论中的B3LYP和PBE0方法,在7个不同的基组下优化得到了4-氨基-1,8-萘二酰亚胺(ANI)的基态几何构型,用CIS/6-31+G(d)方法得到第一激发态几何构型,频率分析无虚频.在此基础上运用HF-CIS,TD-B3LYP和TD-PBE0方法研究了在气相及DMSO,DMF,MeCN,THF,CHCl3和EtOH溶剂中ANI的前线轨道及电子光谱.结果表明,HOMO→LUMO的跃迁是π→π*跃迁.随溶剂极性的增加,其最大吸收和荧光波长红移.用TD-B3LYP/6-31+G(d)方法得到的溶剂中ANI的吸收光谱计算值与实验值吻合性较好,但荧光光谱计算值与实验值有较大差异.进一步经线性拟合校正,ANI在非质子溶剂中的计算值与实验值能较好地吻合.计算显示激发态ANI具有较大的偶极矩,与解释相关荧光分子探针的光诱导电子转移方向选择性现象的光生电场理论一致.     

BPh-2(mqp)的电子结构和光谱性质的含时密度泛函理论研究

采用abinitioHF和DFTB3LYP方法,对配合物BPh₂(mqp)基态结构进行优化,分析了前线分子轨道特征和能级分布.用abinitioCIS方法优化体系激发态结构.用含时密度泛函理论(TD-DFT)对BPh₂(mqp)的电子光谱进行了研究.结果发现,该物质是配体发光配合物,其发光源于mqp配体内π^*→π的电子跃迁.这表明在mqp配体上进行修饰,可有效地影响配合物前线分子轨道分布,达到调整发光波段的目的.     

2-乙酰基-1-甲基吡咯激发态结构动力学及溶剂效应的共振拉曼光谱和密度泛函理论计算

获取了覆盖紫外光谱中A带和B带吸收的共7个不同激发波长的共振拉曼光谱, 并结合密度泛函理论方法研究了2-乙酰基-1-甲基吡咯(2-Ac-NMP)的A带和B带电子激发和Franck-Condon区域结构动力学. 在TD-B3LYP/6-311++G(d,p)计算水平上, A带和B带吸收的跃迁主体为π→π^* . A带和B带共振拉曼光谱分别指认为13个振动模式和8个振动模式的基频、泛频和组合频, 其中C＝O伸缩振动(ν₈)、C3-C4-C5不对称伸缩振动+C2-C6伸缩振动(ν₁₄)及环上CH面内摇摆(ν₁₈)对拉曼光谱强度贡献最大, 表明2-Ac-NMP的S_π激发态结构动力学主要沿反应坐标展开. 考察了溶剂对共振拉曼光谱强度模式的影响, 结果表明, 在同一溶剂中, 随激发波长由长变短, C＝O伸缩振动模(ν₈)的强度呈现出由强变弱再变强的现象. 这种变化规律与Franck-Condon区域S_n/S_π态混合或势能面交叉相关, 并受溶剂的有效调控.     

8-羟基喹啉铍及其衍生物电子光谱性质的含时密度泛函理论研究

运用密度泛函理论(DFT)B3LYP方法和abinitioHF单激发组态相互作用(CIS)法分别优化了有机金属配合物8-羟基喹啉铍(BeQ₂)及其3种衍生物分子的基态及最低激发单重态几何结构.系统分析了分子结构、前线分子轨道特征和能级分布规律以探索电子跃迁机理.应用含时密度泛函理论(TD-DFT)计算分子的电子光谱,揭示了BeQ₂及其衍生物的发光源于配体中π→π^*电子跃迁,指出通过配体修饰可以有效地影响配合物前线分子轨道分布,调整发光波段,并有效提高电荷转移量.     

CH2＝CHCl与O(3P)反应的理论研究

用量子化学密度泛函理论和QCISD (Quadratic configuration interaction calculation)方法, 对O(³P)与CH₂CHCl的反应进行了理论研究. 在UB3LYP/6-311＋＋G(d,p), UB3LYP/6-31＋＋G(3df, 3pd)计算水平上, 优化了反应物、产物、中间体和过渡态的几何构型, 并在UQCISD(T)/6-311＋＋G(2df,2pd)水平上计算了单点能量. 为了确证过渡态的真实性, 在UB3LYP/6-311＋＋G(3df,3pd)水平上进行了内禀坐标(IRC)计算和频率分析, 并确定了反应机理. 研究结果表明, 反应主要产物为CH₂CHO和Cl.     

Theoretical study on structures and stability of SiCP2 isomers

The structures, energetics, spectral parameters and stability of the singlet SiCP₂ isomers are explored at the density functional theory and ab initio levels. Eight isomers connected by ten interconversion transition states are located at the CCSD(T)/6-311G(2d)//B3LYP/6-311G(d)level. The kinetically stable isomers and their relevant interconversion transition states are further refined at CCSD(T)/6-311+G(2df)//QCISD/6-311G(d) level. At QCISD/6-311G(d) level, one four-membered ring isomer cSiPCP and two linear structures PSiCP, SiCPP possess considerable kinetic stability (more than 15 kcal/mol). The valence bond structures of three kinetically stable SiCP₂ isomers are analyzed. The similarities and discrepancies in structure, energy and stability between SiCP₂ and its analogous C₂P₂, Si₂P₂, SiCN₂ and CSiNP molecules are also discussed. The predicted structures and spectroscopic properties are expected to be informative for the identification of the SiCP₂ in the laboratory and space.     

Excited state proton transfer in guanine in the gas phase and in water solution: a theoretical study

Theoretical investigations were performed to study the phenomena of ground and electronic excited state proton transfer in the isolated and monohydrated forms of guanine. Ground and transition state geometries were optimized at both the B3LYP/6-311++G(d,p) and HF/6-311G(d,p) levels. The geometries of tautomers including those of transition states corresponding to the proton transfer from the keto to the enol form of guanine were also optimized in the lowest singlet pipi* excited state using the configuration interaction singles (CIS) method and the 6-311G(d,p) basis set. The time-dependent density function theory method augmented with the B3LYP functional (TD-B3LYP) and the 6-311++G(d,p) basis set was used to compute vertical transition energies using the B3LYP/6-311++G(d,p) geometries. The TD-B3LYP/6-311++G(d,p) calculations were also performed using the CIS/6-311G(d,p) geometries to predict the adiabatic transition energies of different tautomers and the excited state proton transfer barrier heights of guanine tautomerization. The effect of the bulk aqueous environment was considered using the polarizable continuum model (PCM). The harmonic vibrational frequency calculations were performed to ascertain the nature of potential energy surfaces. The excited state geometries including that of transition states were found to be largely nonplanar. The nonplanar fragment was mostly localized in the six-membered ring. Geometries of the hydrated transition states in the ground and lowest singlet pipi* excited states were found to be zwitterionic in which the water molecule is in the form of hydronium cation (H3O(+)) and guanine is in the anionic form, except for the N9H form in the excited state where water molecule is in the hydroxyl anionic form (OH(-)) and the guanine is in the cationic form. It was found that proton transfer is characterized by a high barrier height both in the gas phase and in the bulk water solution. The explicit inclusion of a water molecule in the proton transfer reaction path reduces the barrier height drastically. The excited state barrier height was generally found to be increased as compared to that in the ground state. On the basis of the current theoretical calculation it appears that the singlet electronic excitation of guanine may not facilitate the excited state proton transfer corresponding to the tautomerization of the keto to the enol form.     

An Ab Initio Study of the Complexes Formed by the Positive Acetylene Ion with the Hydrogen and Nitrogen Molecules and the Argon Atom

The complexes formed by the positive acetylene ion with the hydrogen molecule, the nitrogen molecule, and the argon atom are investigated with ab initio calculations using the 6-311G** and the 6-31+G(2df,2pd) basis sets. MP2/6-311G** energies and optimum geometries are obtained, as well as single-point MP3, MP4, and QCISD(T) energies with the MP2/6-311G** optimized geometries. Single-point calculations are performed with the 6-31+G(2df,2pd) basis set at MP2/6-311G** optimized geometries.     

NH0-1+2-3离解能等的高级ab initio计算与评价

A recent experimental determination[1] of the dissociation energies (D0) for H2N-H, H2N+-H and H2N-H+, the ionization energies for NH3 and NH2 resulted in large deviations when compared with those of the earlier values and the QCISD(T)/6-311+G(3df,2p) ab initio calculations. We have performed some higher level ab initio calculations on these data by utilizing the Gaussian 92/DFT and Gaussian 94 pakages of programs and have assessed the available experimental values. Our calculations were carried out at the QCISD (TQ)/aug-cc-pVDZ, G2(QCI), QCISD(T)/6-311 ++G(3df,3pd) and QCISD(T)/aug-cc-pVTZ levels of theory. Geometries were optimized at both of the MP2(full)/6-31G(d) and the MP2(full)/6-31(d,p) levels, and were compared with those of the experiments if available. The MP2(full)/6-31G(d,p) tight-optimized geometries for the neutrals are closer to those of the experiments than those of the MP2 (full)/6-31G(d), and are in excellent agreement with the experimental results as shown in Table 1. In this case, we assumed that the optimized geometries for the cations would be better if p polarization functions are added to the hydrogen atoms. We firstly noted that the symmetry of the NH3+ cation was D3h, other than Cs. as reported in ref.[1]. All of the zero-point energies and the final geometries are calculated at the MP2(full)/6-31G(d,p) level of theory. We have also repeated the QCISD(T )/6-311 + G(3df,2p) calculations of ref. [1], because we could not identify their level of goemetry optimization. It is found that the total energy, -55.244 19 Hartrees, for NH2+(1A1 ) in ref.[1] might be in error. Our result is -55.336 29 Hartrees at the same level of theory. At our highest level [QCISD(T)/aug-cc-pVTZ] of calculations as shown in Table 3, the D0 (temperature at zero Kelvin) values of H2N-H, H2N+-H(3B1for NH2+ ) and H2N- H+ are 4.51, 5.49 and 8.00 eV, respectively. These data reported in re f.[1] were 4.97, 5.59 and 8.41 eV, respectively. Our result on D0(H2N-H) supports the work of ref.[2,3,5,6]. The ionization energies (IE) for NH3 and NH2 (3B1 for NH2+) at our highest level are 10.11 and 11.09 eV while in ref.[1] were 10.16 and 10.78 eV, respectively. For the latter, our result supports the experiment of ref.[3]. Our predicted D0 for HN2+-H and IE for NH2 (1A1 for each NH2+) are 6.80 and 12.39 eV, respectively. These values differ greatly from the predicted values (9.29 and 14.88 eV) of ref.[1] where the total energy of NH2+(1A1) might be in error. The D0 value for HN-H has not been found in ref.[1]. Our result supports the work of ref.[3]. We have also derived all of these values at the temperature of 298K and under the pressure of 101kPa at several levels of thoery as shown in Table 3. On examining the experiment of ref.[1] in detail, it is easy to find that all of the larger deviations might be from a too high value of the appearance potential of proton AP(H+). Indeed, ref.[1] has mentioned that the determintion of AP(H+), due to kinetic shift, would lead to a hihger value for the dissociation energy as has been pointed out by Berkowitz and Ruscic. In this work, we concluded that, besides some mistakes in the theoretical calculations of ref.[1], the dissociation energies for H2N-H and H2N-H+,the IE for NH2 (3B1 for NH2+) might also be unreliable and need to be re-examined.
     

四氢二嵌萘衍生物在不同环境中的电子结构及光谱特征的理论研究

采用密度泛函理论的B3LYP泛函和单激发组态相关方法(CIS)在6-31+G*基组水平下分别优化了四氢二嵌萘衍生物(THP)基态和第一单重激发态的结构,并在此基础上利用含时密度泛函理论计算了THP在不同环境中的电子结构和光谱性质.溶剂效应采用连续极化介质模型(PCM)计算.结构分析发现在酸性环境中,基态时THP二甲氨基上的N原子容易质子化.计算结果表明,尽管THP的第一单重激发态是一个电荷转移态,但是它本身并不发射双荧光.在酸性溶剂中,THP的双峰发射来源于两种物质:跃迁能为2.71eV的发射峰由激发的THP发射,而3.69eV的发射峰则来源于质子化的THP.根据计算结果,我们建议了THP在酸性溶剂中的激发与驰豫途径.理论预测的吸收和发射光谱与实验结果一致.     

DFT法研究3-羟基丙烯醛的双键旋转异构反应机理

利用密度泛函理论(DFT)分别在B3LYP/6-31G**和B3LYP/6-311＋＋G**的计算水平上优化了基态3-羟基丙烯醛分子在双键旋转异构反应过程中的平衡态以及过渡态的几何构型，分析了反应过程中键参数的变化，计算了该反应的内禀反应坐标(IRC)，发现在重排反应途径上存在一个四元环骨架的中间体.通过振动分析对平衡态和过渡态进行了确认，并得到了零点能.计算结果表明，基态3-羟基丙烯醛分子的双键旋转异构反应经过两步完成，第一步反应位垒稍高，第二步反应位垒较低，存在着发生重排反应的可能性.     

Structure and Stability of Isomers of the Promising Interstellar Molecule PC3O

DFT/B3LYP/6-311G(d) and CCSD(T)/6-311G(2d) single-point calculations are carried out for exploring the doublet potential energy surface (PES) of PC₃O, a molecule of potential interest in interstellar chemistry. A total of 29 minima connected by 65 interconversion transition states are located. The structures of the most relevant isomers and transition states are further optimized at the QCISD level followed by CCSD(T) single-point energy calculations. At the CCSD(T)/6-311G(2df)//QCISD/6-311G(d)+ZPVE level, the global minimum is the quasi-linear structure PCCCO 1 (0.0 kcal/mol) with a great kinetic stability of 47.9 kcal/mol, and the cumulenic form features largely in its resonance structures. Moreover, the chainlike isomer OPCCC 3 (64.5) and five-membered-ring species cPCCCO 19 (77.8) possess considerable kinetic stability of about 18.0 kcal/mol. All these three isomers are very promising candidates for future experimental and astrophysical detection. Additionally, a three-membered-ring isomer CC-cCOP 10 (69.6) has slightly lower kinetic stability of around 15 kcal/mol and may also be experimentally observable. Possible formation mechanisms of the four stable isomers in interstellar space are discussed. The present research is the first attempt to study the isomerization and dissociation mechanisms of PC _n O series. The predicted spectroscopic properties, including harmonic vibrational frequencies, dipole moments and rotational constants for the relevant isomers, are expected to be informative for the identification of PC₃O in laboratory and interstellar medium.     

O（~1D）与CF_2HCl反应的理论研究

用量子化学密度泛函理论（DFT）和G3B3方法，对为（~1D）与CF_2HCl的反应 进行了研究，在B3LYP/6-311+G(d)，B3LYP/6-311+G(2df,2pd)和G3B3计算水平上， 优化了反应热能面上各驻点的几何结构，通过内禀反应坐标（IRC）计算和振动分 析，对反应过渡态进行了确认，并确定了反应机理。     

Theoretical Studies on Structures and Stabilities of C4H2+ Isomers

The structures, energies, stabilities and spectroscopies of doublet C₄H₂⁺ cations were explored at the DFT/B3LYP/6-311G(d,p), CCSD(T)/6-311+G(2df,2pd)(single-point), and G3B3 levels. Ten minimum isomers including the chainlike, three-member-ring, and four-member-ring structures are interconverted by means of 15 interconversion transition states. The potential energy surface was investigated. At the CCSD(T)/6-311+G(2df,2pd) and G3B3 levels, the global minimum isomer was found to be a linear HCCCCH. The structures of the stable isomer and its relevant transition state are further optimized at the QCISD/6-311G(d,p) level. The bonding nature and structure of isomer HCCCCH were analyzed.