取代基对3(5)-(9-蒽基)吡唑光学性质的影响

使用密度泛函理论(DFT)B3LYP/6-31G(d)方法优化得到了3(5)-(9-蒽基)吡唑及其衍生物的基态(S₀)分子结构, 使用单激发组态相互作用(CIS)/6-31G(d)方法优化得到这些分子的第一单重激发态(S₁)的几何结构, 并使用含时密度泛函理论(TD-DFT)B3LYP/6-311++G(d,p)方法计算了它们的吸收和发射光谱. 计算结果表明, 与3(5)-(9-蒽基)吡唑相比, 无论取代基是吸电子基团还是供电子基团, 衍生物的吸收和发射峰均发生红移, 并且当取代基―R=―BH₂, ―CCl₃, ―CHO, ―NH₂时衍生物有较长的吸收波长和发射波长.     

8-巯基喹啉阴离子的锌配合物及其衍生物的电子光谱性质的含时密度泛函理论研究

采用从头算(ab initio)和密度泛函理论(DFT B3LYP)方法, 对配合物8-巯基喹啉锌Zn(tq)₂及其5种衍生物基态结构进行优化, 用含时密度泛函理论(TD-DFT/B3LYP)及6-31＋G(d)基组计算吸收光谱; 同时用ab initio HF 单激发组态相互作用(CIS)法在6-31G(d)基组上优化其最低激发单重态几何结构, 用含时密度泛函理论计算发射光谱. 结果表明, 电子在基态与激发态间的跃迁, 主要是在配体8-巯基喹啉(tq)环内的电荷转移, 电子从含S的苯硫酚环转移至含N的吡啶环上; 吸收光谱和发射光谱的计算值与实验值基本符合. 该类配合物都是优良的电子传输材料, 改变金属离子和取代基均可以调控发光材料的光谱波段.     

CnBδ(δ=0, ±1; n =1～6)团簇的结构、稳定性和光谱

采用密度泛函理论(DFT)的B3LYP方法,在6-31G*和6-311＋G(3df)水平上对C_nB(n=1～6)团簇及其阴离子和阳离子的几何构型和电子结构进行了优化和振动频率计算.得到了C_nB(n=1～6)团簇的电离能,绝热电子亲合势以及C_nB^δ(δ=0,±1)团簇的能隙.结果表明C_nB^－(n=1～6)团簇的基态构型均为线形,这与等电子的C_n簇合物的结构是一致的; C_nB(n=1～6)团簇的基态构型中,除C₂B为不对称的三角形,C₆B为具有C_2v对称性的环状结构外,其余均为线形结构.阳离子团簇中n=2、3、6的基态结构具有C_2v对称性外,其它几个均为线形结构.从几何参数和振动频率上发现,采用密度泛函B3LYP方法在6-311＋G(3df)和6-31G*两种基组上计算得到的键长参数和振动频率非常接近,说明B3LYP方法在计算C_nB簇合物结构参数上对于基组的选择是不太敏感的.通过对C_nB(n=1～6)的光电子能谱性质的研究发现,C₄B容易获得一个电子形成阴离子团簇,但失去一个电子是很困难的,这与实验上观测到的结果非常吻合.     

2，7′-(乙烯基)-二-8-羟基喹啉及其金属有机配合物的密度泛函和自然键轨道理论计算

采用密度泛函理论(DFT),在B3LYP/6-31G水平上对2,7′-(乙烯基)-二-8-羟基喹啉(2,7′-Ethq₂)及其3种金属M(M=Zn,Mg,Be)有机配合物M(2,7′-Ethq₂)₂的结构进行了全优化,并用ZINDO和TDDFT方法计算了它们的吸收光谱。同时,利用自然键轨道理论(NBO)对分子内氢键进行了分析。结果表明,光谱计算值与实验值基本符合,该类化合物均具有较大的电子亲和能,改变中心金属原子对配合物吸收光谱性质影响不大。和2,7′-Ethq₂相比,M(2,7′-Ethq₂)₂的吸收光谱产生明显红移。2,7′-Ethq₂及其M(2,7′-Ethq₂)₂分子内存在较强的氢键,氢键与环上的碳原子形成五元环,分子内氢键的存在使分子的稳定性增加。     

芴-硫芴共聚物的电子结构和光谱性质的理论研究

用DFT-B3LYP方法对低聚物(PF30T)_n [n(芴)∶n(硫芴)＝2∶1, 物质的量之比, n＝1～4], (PF50T)_n [n(芴)∶n(硫芴)＝1∶1, 物质的量之比, n＝1～4]体系全优化, 得到两系列低聚物的电离能(IP_(a,v))、电子亲和势(EA_(a,v))、空穴抽取能(HEP)、电子抽取能(EEP). 在此基础上用ZINDO和TD-DFT方法计算吸收光谱, 分析了两个系列的HOMO-LUMO能隙随着n递增的变化趋势及硫芴含量对低聚物电子结构和光谱性质的影响, 推断了高聚物的电子和光谱性质. 用ab initio CIS方法优化了低聚物的S₁激发态结构并分析了其与发射光谱的关系. 研究显示: 2,8位引入的硫芴基团, 破坏了链的共轭, 而且随着硫芴含量的增加, HOMO-LUMO能隙变大, 光谱蓝移; 激发态结构趋于共面化.     

苯乙烯基吡啶类化合物几何结构与光谱的理论研究

采用量子化学密度泛函理论(DFT)方法分别在B3LYP/6-31G*, 6-31G**, 6-31＋G*水平上对苯乙烯基吡啶类化合物进行计算研究. 通过在相同水平下的振动频率分析发现苯乙烯基吡啶类化合物具有C₁对称性, 酯基的碳氧原子与苯环形成不同的离域大π键, 空间位阻和共轭效应使得两苯环处于两个不同平面, 二面角在60°与62°之间. 使用含时密度泛函理论(TD-DFT)方法计算第一激发态的电子垂直跃迁能, 得到最大吸收波长λ_max. 计算结果表明末端烷基链的长度对该类化合物的几何结构与振动光谱、电子光谱无影响.     

藏药五脉绿绒蒿碱结构和性质的理论研究

五脉绿绒蒿碱是一种从藏药五脉绿绒蒿中提取并已确认结构的新的生物碱. 采用密度泛函理论(DFT)和从头算(ab initio)方法, 在HF/6-31G*和B3LYP/6-31G*水平下全优化计算了该化合物的分子几何构型和电子结构; 依据Onsager自恰反应场(SCRF)模型考察了五脉绿绒蒿碱在氯仿、丙酮、二甲亚砜及水等溶剂中的溶剂化作用; 基于气相优化结构进行了B3LYP/6-31G*振动分析与红外光谱计算, 进一步按照统计力学原理求得了298～1500 K温度范围内该化合物的热力学性质. 此外, 还讨论了五脉绿绒蒿碱的分子结构与药效的关系.     

6-羟基-5,12-萘并萘醌及其CH3, C6H5取代衍生物的光致变色理论研究

用密度泛函理论的B3LYP/6-31G方法和从头算的CIS/6-31G方法分别研究了6-羟基-5,12-萘并萘醌及其CH₃, C₆H₅取代衍生物基态和激发态的异构化反应，.对反应势能面的研究发现, 在光异构化反应中化合物M21和M21、M31的基态和激发态虽然都可以构成四能级反应过程, 但由于M21异构化过程的活化能较高, 使其所构成的四能级反应难以进行, 这就从理论上解释了迁移基团为甲基的M21变色性能低于迁移基团为苯基的M31的实验结果. 此外用TD/B3LYP方法在溶剂存在下计算了上述化合物的紫外吸收光谱和荧光发射光谱, 计算所得到的光谱数据与实验值基本一致, 与光异构化反应的光激发条件相符合.     

6-羟基-5,12-萘并萘醌及其CH3, C6H5取代衍生物的光致变色理论研究

用密度泛函理论的B3LYP/6-31G方法和从头算的CIS/6-31G方法分别研究了6-羟基-5,12-萘并萘醌及其CH₃, C₆H₅取代衍生物基态和激发态的异构化反应，.对反应势能面的研究发现, 在光异构化反应中化合物M21和M21、M31的基态和激发态虽然都可以构成四能级反应过程, 但由于M21异构化过程的活化能较高, 使其所构成的四能级反应难以进行, 这就从理论上解释了迁移基团为甲基的M21变色性能低于迁移基团为苯基的M31的实验结果. 此外用TD/B3LYP方法在溶剂存在下计算了上述化合物的紫外吸收光谱和荧光发射光谱, 计算所得到的光谱数据与实验值基本一致, 与光异构化反应的光激发条件相符合.     

锂在共轭双键高分子中的电化学嵌入反应III. 锂嵌入聚噻吩的量子化学计算†

采用Gaussian软件和HF方法, 通过从头计算(ab initio)法选取4-31G基组计算锂离子嵌入聚噻吩过程中结构与结合能的变化关系. 发现噻吩聚合时主要生成三或四聚合物. 聚合物在Li原子(或Li^＋离子)嵌入后, 聚噻吩间距离明显变小, 同时发生电荷转移, 形成稳定嵌合物; 并使噻吩环的C-α—C-β键级变小. 同时, 研究了锂离子(或原子)嵌入后体系的HOMO, LUMO能级. 聚噻吩在嵌入锂离子时LUMO轨道能级变为负值, 成为电池反应得电子的正极. 而金属Li₂ 释放Li^＋后的Li^－的HOMO能级为＋0.7427 eV, 则成为给电子的负极. 由此, 可以完成由锂/聚噻吩在高氯酸锂电解质中组成的放电过程, 并提出嵌合键级概念用来表征锂在聚噻吩间的结合程度.     

Mismatch base pairing of the mutagen 8‐oxoguanine and its derivatives with adenine: A theoretical search for possible antimutagenic agents

Molecular geometries of 8‐oxoguanine (8OG), those of its substituted derivatives with the substitutions CH₂, CF₂, CO, CNH, O, and S in place of the N7H7 group, adenine (A), and the base pairs of 8OG and its substituted derivatives with adenine were optimized using the RHF/6‐31+G* and B3LYP/6‐31+G* methods in gas phase. All the molecules and their hydrogen‐bonded complexes were solvated in aqueous media employing the polarized continuum model (PCM) of the self‐consistent reaction field (SCRF) theory using the RHF/6‐31+G* and B3LYP/6‐31+G* methods. The optimized geometrical parameters of the 8OG‐A base pair at the RHF/6‐31+G* and B3LYP/6‐31+G* levels of theory agree satisfactorily with those of an oligonucleotide containing the base pair found from X‐ray crystallography. The pattern of hydrogen bonding in the CF₂‐ and O‐substituted 8OG‐A base pair is of Watson–Crick type and that in the unsubstituted and CH₂‐, CNH‐, and S‐substituted base pairs is of Hoogsteen type. In the CO‐substituted base pair, the hydrogen bonding pattern is of neither Watson–Crick nor Hoogsteen type. The CF₂‐substitution appears to introduce steric hindrance for stacking of DNA bases. On the basis of these results, it appears that among all the substituted 8OG molecules considered here, the O‐substituted derivative may be useful as an antimutagenic drug. It is, however, subject to experimental verification. © 2004 Wiley Periodicals, Inc. Int J Quantum Chem, 2005     

A computational mechanistic study of the deamination reaction of melamine

A detailed computational study of the deamination reaction of melamine by OH^–, n H₂O/OH^–, n H₂O (where n = 1, 2, 3), and protonated melamine with H₂O, has been carried out using density functional theory and ab initio calculations. All structures were optimized at M06/6‐31G(d) level of theory, as well as with the B3LYP functional with each of the basis sets: 6‐31G(d), 6‐31 + G(d), 6‐31G(2df,p), and 6‐311++G(3df,3pd). B3LYP, M06, and ω B97XD calculations with 6‐31 + G(d,p) have also been performed. All structures were optimized at B3LYP/6‐31 + G(d,p) level of theory for deamination simulations in an aqueous medium, using both the polarizable continuum solvation model and the solvation model based on solute electron density. Composite method calculations have been conducted at G4MP2 and CBS‐QB3. Fifteen different mechanistic pathways were explored. Most pathways consisted of two key steps: formation of a tetrahedral intermediate and in the final step, an intermediate that dissociates to products via a 1,3‐proton shift. The lowest overall activation energy, 111 kJ mol^?1 at G4MP2, was obtained for the deamination of melamine with 3H₂O/OH^?.     

Effects of Epoxidation and Nitration on Ballistic Properties of FOX‐7 – A DFT Study

1, 1‐Diamino‐2, 2‐dinitroethylene (FOX‐7) has received increasing attention since it was industrialized in the late 1990s. It has lower sensitivity and comparable performance to RDX. This paper presents ballistic properties of FOX‐7, its mono and dinitro derivatives and their epoxide derivatives computationally. The structures were optimized at the B3LYP/6‐31G(d, p) level and the bond lengths were calculated. The calculated data for FOX‐7 are compatible with the literature one. We have investigated the bond dissociation energies of the molecules. Mulliken electro negativities (χ_M) and chemical hardness (η) were reviewed using Frontier Molecular Orbitals at HF/6‐31G(d, p)//B3LYP/6‐31G(d, p) theoretical level. The detonation performance analyses were done using empirical Kamlet‐Jacobs equations. Additionally, power index values were calculated. All the compounds considered in the present article are powerful candidates for high energy materials.     

Computational studies on tetrazole derivatives as potential high energy materials

The tetrazole is an important functionality of the most of energetic materials due to 80% nitrogen content, stability, and high enthalpy of formation. The present structure–property relationship study focuses on the optimized geometries of tetrazole derivatives obtained from density functional theory (DFT) calculations at B3LYP/6-31G* levels. The heat of formation (HOF) of tetrazole derivatives have been calculated by designing the appropriate isodesmic reactions. The increase in nitro groups on azole rings shows the remarkable increase in HOF. Density has been predicted by using CVFF force field. Increase in the nitro group increases the density. Detonation properties of the designed compounds were evaluated by using the Kamlet–Jacobs equation based on predicted densities and HOFs. Designed tetrazole derivatives show detonation velocity (D) over 8 km/s and detonation pressure (P) of about 32 GPa. Thermal stability was evaluated via bond dissociation energies (BDE) of the weakest C–NO₂ bond at B3LYP/6-31G* level. Charge on the nitro group has been used to assess the sensitivity correlation. Overall, the study implies that designed compounds of this series are found to be stable and expected to be the novel candidates of high energy materials (HEMs).     

Structure and electronic properties of tris(4-hydroxy-1,5-naphthyridinato) aluminum (AlND3) and its methyl derivatives: a theoretical study

Molecular structures of the ground state (S₀) of tris(4-hydroxy-1,5-naphthyridinato) aluminum (AlND3) and its methyl derivatives have been optimized using B3LYP/6-31G(D) and the first singlet excited state (S₁) using CIS/6-31G(D) method, respectively. The frontier molecular orbitals characteristics of these Al-complexes have been analyzed systematically in order to understand the electronic transitions. It is seen from these results that in all these complexes, like in earlier reported mer-Alq3, the highest occupied molecular orbital (HOMO) is localized on the pyridine-4-ol ring of A-ligand while lowest unoccupied molecular orbital (LUMO) is on the pyridyl ring of B-ligand for S₀ states irrespective of the methyl substitution present on the ligands. The absorption and emission wavelengths have been evaluated at the TD-PBE0/6-31G(D) level and found to be comparable with the experiment. The charge transfer integrals have been calculated for AlND3, and results reveals that electron transport is larger than hole transport. The reorganization energies have been calculated at B3LYP/6-31G(D) level for these complexes, and the results show that the charge mobilities of these complexes are comparable with mer-Alq3.     

TD-DFT Study on the Electronic Spectrum Properties of 2,7＇-（Ethylene）- bis-8-hydroxyquinoline and Its Derivatives

The structures of 2,7＇-（ethylene）-bis-8-hydroxyquinoline and its derivatives were optimized at the ground states using ab initio HF and B3LYP methods. At the same time, the molecular structures of the first singlet excited state for 2,7＇-（ethylene）-bis-8-hydroxyqulnoline and its derivatives were optimized by CIS/6-31G（d）. The absorption and emission spectra based on the above structures were obtained by the time-dependent density functional theory （TD-DFT） by the B3LYP method with the 6-31G（d） basis set. The calculated results of luminescence originate from the electronic transition from the hydroxphenol ring of 8-hydroxyquinoline A to the pyridine ring of 8-hydroxyquinoline B. Their luminescence wave bands can be tuned by different substituents on the ligand of 8-hydroxyquinoline.     

Structures and stability of N9, N9 − and N9 + clusters

 Ab initio molecular orbital calculations for N₉, N⁻ ₉ and N⁺ ₉ isomers were carried out at the HF/ 6-31G*, B3PW91/6-31G*, B3LYP/6-31G* and MP2/ 6-31G* levels of theory. Stable equilibrium geometric structures were determined by harmonic vibrational frequency analyses at the HF/6-31G*, B3PW91/6-31G* and B3LYP/6-31G* levels of theory. The most stable free-radical N₉ cluster is structure 1 with C _{2 v} symmetry and that of anion N⁻ ₉ is structure 3 with C _s symmetry. Only one stable structure of the N⁺ ₉ cation with C _{2 v} symmetry was predicted. Their potential application as high-energy-density materials has been examined. Received: 15 June 1999 / Accepted: 11 October 1999 / Published online: 14 March 2000     

Computational DFT studies on a series of toluene derivatives as potential high energy density compounds

Based on the full optimized molecular geometric structures at B3LYP/6-31G**, B3LYP/6-31+G**, B3P86/6-31G**, and B3P86/6-31+G** levels, the densities (ρ), detonation velocities (D), and pressures (P) for a series of toluene derivatives, as well as their thermal stabilities, were investigated to look for high energy density compounds (HEDCs). The heats of formation (HOFs) are also calculated via designed isodesmic reactions. The calculations on the bond dissociation energies (BDEs) indicate that the BDEs of the initial scission step are between 48 and 59 kcal/mol, and pentanitrotoluene is the most reactive compound, while 2,4,6-trinitrotoluene is the least reactive compound for toluene derivatives studied. A good linear relationship between BDE/E and impact sensitivity is also obtained. The condensed phase HOFs are calculated for the title compounds. These results would provide basic information for the further studies of HEDCs. The detonation data of pentanitrotoluene show that it meets the requirement for HEDCs.