Quantum chemical study on excited states and electronic coupling matrix element in a catechol-bridge-dicyanoethylene system

The excited states of a donor-bridge-acceptor (DBA) model system have been investigated using time-dependent density-functional theory (TD-DFT) in vacuo and in solution. It is found that the MPW1PW91 functional always gives higher excitation energies than those with a B3LYP functional. Results from both TD-B3LYP and TD-MPW1PW91 are found consistent with the experimental observations. The two most intense absorptions of the DBA system, one resulting from the local excitation of catechol moiety and the other from that of dicyanoethylene, possess the pipi* transition feature. It seems that the solvent polarity does not remarkably influence the positions of absorption peaks. The spectroscopic properties of isolated donor, acceptor, and bridge and the donor-bridge compound have been investigated at the TD-B3LYP/6-31+G* and TD-MPW1PW91/6-31+G* levels. Results indicate that the donor and the acceptor are weakly coupled with the bridge. Therefore, it is more likely that the electron transfer takes place through a superexchange mechanism. In addition, we calculate the electronic coupling matrix elements according to the generalized Mulliken-Hush theory, and the detailed analyses also predict that the strong absorptions are due to the local excitation of the DBA system.     

A computational study of the ground and excited state structure and absorption spectra of free-base N-confused porphine and free-base N-confused tetraphenylporphyrin

Computational investigations into the ground and singlet excited-state structures and the experimental ground-state absorption spectra of N-confused tetraphenylporphyrin tautomers 1e and 1i and N-confused porphines (NCP) 2e and 2i have been performed. Structural data for the ground state, performed at the B3LYP/6-31G(d), B3LYP/6-31+G(d)//B3LYP/6-31G(d), and B3LYP/6-311+G(d)//B3LYP/6-31G(d) levels, are consistent with those performed at lower levels of theory. Calculations of the gas-phase, ground-state absorption spectrum are qualitatively consistent with condensed phase experiments for predicting the relative intensities of the Q(0,0) and Soret bands. Inclusion of implicit solvation in the calculations substantially improves the correlation of the energy of the Soret band with experiment for both tautomers (1e, 435 nm predicted, 442 nm observed in DMAc; 1i, 435 nm predicted, 437 nm observed in CH2Cl2). The x- and y-polarized Q-band transitions were qualitatively reproduced for 1e in both the gas phase and with solvation, although the low-energy absorption band in 1i was predicted at substantially higher energy (646 nm in the gas phase and 655 nm with solvation) than observed experimentally (724 nm in CH2Cl2). Franck-Condon state and equilibrated singlet excited-state geometries were calculated for unsubstituted NCP tautomers 2e and 2i at the TD-B3LYP/SVP and TD-B3LYP/TZVP//TD-B3LYP/SVP levels. Electronic difference density plots were calculated from these geometries, thereby indicating the change of electron density in the singlet excited states. Adiabatic S1 and S2 geometries of these compounds were also calculated at the TD-B3LYP/SVP level, and the results indicate that while 2i is a more stable ground-state molecule by approximately 7.0 kcal mol-1, the energy difference for the S1 excited states is only approximately 1.0 kcal mol-1 and is 6.1 kcal mol-1 for the S2 excited states.     

三种席夫碱-Ni(Ⅱ)配合物的电子结构和吸收光谱的理论计算

采用B3LYP/6-31+G(d)-LANL2DZ方法优化了三种Ni(Ⅱ)的席夫碱配合物基态的几何构型,并在相同水平下进行了频率分析以确认稳定点的性质;利用含时密度泛函理论和极化连续介质模型(PCM),按TDB3LYP/6-31+G(d)-LANL2DZ水平计算了目标配合物在N,N-二甲基甲酰胺溶剂中的电子结构和吸收光谱.计算结果表明,配体中间位甲氧基的存在使配合物A具有较大的HOMO-LUMO能级差;且三种Ni(Ⅱ)配合物的S0→S1态的跃迁能按照A→B→C的顺序依次降低.     

DFT study of conjugational electronic structures of aminoalkyl end-capped oligothiophenes up to octamers

The molecular geometries and electronic properties of a series of bis(aminoalkyl) end-capped oligothiophenes (BRnTs) were investigated by means of the density functional theory (DFT). The calculations were performed on dimers up to octamers in the neutral and ionic species using the B3LYP/6-31G(d,p) level of theory. The results obtained show that the conjugated systems in the p- and n-doped oligomers had more aromaticity, with expanded and planar chains. The calculated energy gap values between the frontier molecular orbitals for the end-capped oligomers were larger than those for the unsubstituted oligomers, in which with increase in the oligomer chain length, the conduction band gap decreased. The calculated first excitation energies of BRnTs at the TD-B3LYP/6-31G(d,p) level indicated that both doped oligomers (p- and n-type) had lower excitation energies than the neutral states, and that they displayed red shifts in their absorption spectra. Moreover, the results obtained for the natural bond orbital (NBO) analysis showed that closing the end-side oligothiophene chains with the aminoalkyl groups eased the hole or electron transfer, owning to better charge delocalization through the backbone structures of BRnTs.     

In-depth quantum chemical investigation of electro-optical and charge-transport properties of trans-3-(3,4-dimethoxyphenyl)-2-(4-nitrophenyl)prop-2-enenitrile

The structural, electro-optical and charge-transport properties of compound trans-3-(3,4-dimethoxyphenyl)-2-(4-nitrophenyl)prop-2-enenitrile (DMNPN) were studied using quantum chemical methods. The neutral, cation and anion molecular geometries were optimized in the ground state using density functional theory (DFT) at the restricted and unrestricted B3LYP/6-31G** level of theory. The excited state geometries were optimized by applying time-dependent DFT at the TD-B3LYP/6-31G** level of theory. The absorption and fluorescence wavelengths were calculated at the TD-CAM-B3LYP/6-31G** and TD-LC-BLYP/6-31G** levels of theory. The distribution pattern of the charge densities on the highest occupied molecular orbitals (HOMOs) and lowest unoccupied molecular orbitals (LUMOs) are discussed. Intramolecular charge transfer was observed from the dimethoxyphenyl to (nitrophenyl)prop-2-enenitrile moieties. The detailed charge-transport behavior of the DMNPN molecule is investigated based on its ionization potential, electron affinity, hole and electron reorganization energies, hole and electron-transfer integrals, and hole and electron intrinsic mobilities. The total/partial densities of states and structure–property relationship are discussed in detail. The higher computed hole intrinsic mobility than electron intrinsic mobility reveals that DMNPN is an efficient hole-transport material.     

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.     

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的吡啶环上; 吸收光谱和发射光谱的计算值与实验值基本符合. 该类配合物都是优良的电子传输材料, 改变金属离子和取代基均可以调控发光材料的光谱波段.     

Electronic, optical, and charge transfer properties of donor–bridge–acceptor hydrazone sensitizers

The ground state geometries have been computed by using density functional theory (DFT) at B3LYP/6-31G*, B3LYP/6-31G**, and PCM-B3LYP/6-31G* level of theories. The highest occupied molecular orbitals (HOMOs) are delocalized on whole of the molecule and the lowest unoccupied molecular orbitals (LUMOs) are localized on the tricarbonitrile. The lowest HOMO and LUMO energies have been observed for Dye1 while highest for Dye4. The LUMO energies of Dye1–Dye4 are above the conduction band of TiO₂ and HOMOs are below the redox couple. The absorption spectra have been computed in solvent (methanol) and without solvent by using time-dependant DFT at TD-B3LYP/6-31G*, TD-B3LYP/6-31G**, and PCM-TD-B3LYP/6-31G* level of theories. The calculated maximum absorption wavelengths of the spectra in methanol are in good agreement with experimental evidences. The maximum absorption wavelengths of new designed sensitizers are red shifted compared to parent molecule. The electronic coupling constant and electron injection have been computed by first principle investigations. The improved electronic coupling constant and electron injection revealed that new modeled systems would be efficient sensitizers.     

Formation of sugar radicals in RNA model systems and oligomers via excitation of guanine cation radical

In previous work, we have shown that photoexcitation of guanine cation radical (G*+) in frozen aqueous solutions of DNA and its model compounds at 143 K results in the formation of neutral sugar radicals with substantial yield. In this report, we present electron spin resonance (ESR) and theoretical (DFT) evidence regarding the formation of sugar radicals after photoexcitation of guanine cation radical (G*+) in frozen aqueous solutions of one-electron-oxidized RNA model compounds (nucleosides, nucleotides and oligomers) at 143 K. Specific sugar radicals C5'*, C3'* and C1'* were identified employing derivatives of Guo deuterated at specific sites in the sugar moiety, namely, C1'-, C2'-, C3'- and C5'-. These results suggest C2'* is not formed upon photoexcitation of G*+ in one-electron-oxidized Guo and deuterated Guo derivatives. Phosphate substitution at C5'- (i.e., in 5-GMP) hinders formation of C5'* via photoexcitation at 143 K but not at 77 K. For the RNA-oligomers studied, we observe on photoexcitation of oligomer-G*+ the formation of mainly C1'* and an unidentified radical with a ca. 28 G doublet. The hyperfine coupling constants of each of the possible sugar radicals were calculated employing the DFT B3LYP/6-31G* approach for comparison to experiment. This work shows that formation of specific neutral sugar radicals occurs via photoexcitation of guanine cation radical (G*+) in RNA systems but not by photoexcitation of its N1 deprotonated species (G(-H)*). Thus, our mechanism regarding neutral sugar formation via photoexcitation of base cation radicals in DNA appears to be valid for RNA systems as well.     

Addition reaction of adamantylideneadamantane with Br2 and 2Br2: a computational study

Ab initio calculations were carried out for the reaction of adamantylideneadamantane (Ad=Ad) with Br2 and 2Br2. Geometries of the reactants, transition states, intermediates, and products were optimized at HF and B3LYP levels of theory using the 6-31G(d) basis set. Energies were also obtained using single point calculations at the MP2/6-31G(d)//HF/6-31G(d), MP2/6-31G(d)//B3LYP/6-31G(d), and B3LYP/6-31+G(d)//B3LYP/6-31G(d) levels of theory. Intrinsic reaction coordinate (IRC) calculations were performed to characterize the transition states on the potential energy surface. Only one pathway was found for the reaction of Ad=Ad with one Br2 producing a bromonium/bromide ion pair. Three mechanisms for the reaction of Ad=Ad with 2Br2 were found, leading to three different structural forms of the bromonium/Br3- ion pair. Activation energies, free energies, and enthalpies of activation along with the relative stability of products for each reaction pathway were calculated. The reaction of Ad=Ad with 2Br2 was strongly favored over the reaction with only one Br2. According to B3LYP/6-31G(d) and single point calculations at MP2, the most stable bromonium/Br3- ion pair would form spontaneously. The most stable of the three bromonium/Br3- ion pairs has a structure very similar to the observed X-ray structure. Free energies of activation and relative stabilities of reactants and products in CCl4 and CH2ClCH2Cl were also calculated with PCM using the united atom (UA0) cavity model and, in general, results similar to the gas phase were obtained. An optimized structure for the trans-1,2-dibromo product was also found at all levels of theory both in gas phase and in solution, but no transition state leading to the trans-1,2-dibromo product was obtained.     

AN EPR STUDY OF π-CATION RADICALS IN DINUCLEOSIDE PHOSPHATES AND DNA PRODUCED BY PHOTOIONIZATION

In this work we have produced the π-cation radicals of a number of nucleotides, dinucleoside phosphates, and DNA in aqueous glasses (8M NaCIO₄) by photoionization and investigated these species by EPR spectroscopy. Results found for nucleotides and dinucleoside phosphates containing one type of DNA base, e.g. TpT, GpG, or dApdA, were used in the analysis of spectra found for mixed dinucleoside phosphates, e.g. TpdG. For TpdG and TpdA in neutral glasses photoionization takes place from the purine base and no transfer of charge to the pyrimidine base is found. In basic conditions both the adenine and thymine π-cations are observed in TpdA. In both neutral and basic conditions the results found for mixed dinucleoside phosphates containing guanine show that the guanine cation is formed preferentially by photolysis. This result was found to extend to DNA. Photolysis of DNA in 8M NaC1O₄ produced principally the guanine cation. Computer simulations using parameters determined by other workers from a study of γ-irradiated oriented DNA closely match the spectrum found in this work attributed to the guanine cation in dinucleoside phosphates and DNA. This work thus confirms the presence of the guanine cation in γ-irradiated DNA.     

线型碳链PC2nP的结构和电子光谱的密度泛函理论研究

应用密度泛函理论,在B3LYP/6-31G**和B3LYP/6-311G**水平上优化得到了线型簇合物PC2nP(n=1-10)的基态平衡几何构型,计算了它们的谐振动频率.在基态平衡构型下,利用含时密度泛函理论,计算得到了簇合物PC2nP(n=1-10)的垂直激发能和相应的振子强度,导出了激发能与体系大小n的解析关系式.     

On the origin of cis selectivity in the cyclization of N-protected 2-substituted 3-aza-5-hexenyl radicals: a density functional study

Cyclization of the N-dimethylphosphinoyl-2-methyl-3-aza-5-hexenyl radical has been studied at the UB3LYP/6-31+G(d)//UB3LYP/6-31G(d) hybrid density functional level. The corresponding radical precursor has been synthesized and found to give cis/trans ratios of up to 10/1 in reductive radical cyclizations. The relative energies of reactant and transition state conformers were determined. In discord with the Beckwith-Houk model, it has been found that chair-axial transition states, which lead to cis products, are lowest in energy, rationalizing the observed experimental diastereoselectivity.     

尿酸质子转移反应的理论研究I——分子内质子转移研究

使用量子化学中的Hartree-Fock方法和密度泛函理论中的B3LYP方法，分别在3－21G^*和6－31G（d)水平上，计算了尿酸分子从三羰基异构体向三羟基异构体的转化。结果表明，转化过程经历了单羟基和双羟基异构体2种中间物和3种过渡态时的分子内质子转移（IPT），转移中的H原邻近的N，O和C原子形成了具有四元环结构的过渡态。随着IPT的进行，N－H键逐渐被削弱和断裂，O－H键则逐渐生成。3个反应的活化能分别为190.3kJ/mol,181.4kJ/mol和249.9kJ/mol(B3LYP/6-31G(d))。较高的活化能表明在室温下，无催化剂的IPT难以进行。     

Computational study of the deamination reaction of cytosine with H2O and OH-

The mechanism for the deamination reaction of cytosine with H(2)O and OH(-) to produce uracil was investigated using ab initio calculations. Optimized geometries of reactants, transition states, intermediates, and products were determined at RHF/6-31G(d), MP2/6-31G(d), and B3LYP/6-31G(d) levels and for anions at the B3LYP/6-31+G(d) level. Single-point energies were also determined at B3LYP/6-31+G(d), MP2/GTMP2Large, and G3MP2 levels of theory. Thermodynamic properties (DeltaE, DeltaH, and DeltaG), activation energies, enthalpies, and free energies of activation were calculated for each reaction pathway that was investigated. Intrinsic reaction coordinate analysis was performed to characterize the transition states on the potential energy surface. Two pathways for deamination with H(2)O were found, a five-step mechanism (pathway A) and a two-step mechanism (pathway B). The activation energy for the rate-determining steps, the formation of the tetrahedral intermediate for pathway A and the formation of the uracil tautomer for pathway B, are 221.3 and 260.3 kJ/mol, respectively, at the G3MP2 level of theory. The deamination reaction by either pathway is therefore unlikely because of the high barriers that are involved. Two pathways for deamination with OH(-) were also found, and both of them are five-step mechanisms. Pathways C and D produce an initial tetrahedral intermediate by adding H(2)O to deprotonated cytosine which then undergoes three conformational changes. The final intermediate dissociates to product via a 1-3 proton shift. Deamination with OH(-), through pathway C, resulted in the lowest activation energy, 148.0 kJ/mol, at the G3MP2 level of theory.     

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-hydroxyquinoline 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.     

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.     

单线态卡宾衍生物与臭氧反应机理的密度泛函数理论研究

用密度泛函理论中B3LYP方法,在6-31G(d)水平上优化了单线态卡宾衍生物(CX~2,X=F,Cl,Br)与臭氧反应过程中驻点的几何构型,通过振动分析对反应过渡态和中间体构型进行确认,对单点进行了CCSD(T)/6-31G(d)计算,并进行了零点能(ZPE)校正,同时对反应绝对速率常数进行了理论计算,研究结果表明：相对CF~2和CCl~2而言CBr~2对臭氧的损耗较大,卡宾衍生物与臭氧反应过程均为强放热反应。     

TD-CI simulation of the electronic optical response of molecules in intense fields II: comparison of DFT functionals and EOM-CCSD

Time-dependent configuration interaction (TD-CI) simulations can be used to simulate molecules in intense laser fields. TD-CI calculations use the excitation energies and transition dipoles calculated in the absence of a field. The EOM-CCSD method provides a good estimate of the field-free excited states but is rather expensive. Linear-response time-dependent density functional theory (TD-DFT) is an inexpensive alternative for computing the field-free excitation energies and transition dipoles needed for TD-CI simulations. Linear-response TD-DFT calculations were carried out with standard functionals (B3LYP, BH&HLYP, HSE2PBE (HSE03), BLYP, PBE, PW91, and TPSS) and long-range corrected functionals (LC-ωPBE, ωB97XD, CAM-B3LYP, LC-BLYP, LC-PBE, LC-PW91, and LC-TPSS). These calculations used the 6-31G(d,p) basis set augmented with three sets of diffuse sp functions on each heavy atom. Butadiene was employed as a test case, and 500 excited states were calculated with each functional. Standard functionals yield average excitation energies that are significantly lower than the EOM-CC, while long-range corrected functionals tend to produce average excitation energies slightly higher. Long-range corrected functionals also yield transition dipoles that are somewhat larger than EOM-CC on average. The TD-CI simulations were carried out with a three-cycle Gaussian pulse (ω = 0.06 au, 760 nm) with intensities up to 1.26 × 10(14) W cm(-2) directed along the vector connecting the end carbons. The nonlinear response as indicated by the residual populations of the excited states after the pulse is far too large with standard functionals, primarily because the excitation energies are too low. The LC-ωPBE, LC-PBE, LC-PW91, and LC-TPSS long-range corrected functionals produce responses comparable to EOM-CC.