构型与构象对苯双自由基体系基态自旋多重度及其稳定性的影响

利用AM1-CI方法计算了构象对邻、间、对二取代苯双自由基体系基态自旋多重度及其稳定性的影响. 结果表明单 - 三重态能量差(△ES－T)和部分占据分子轨道的能量劈裂(△EPOMO)随自由基与苯环间的二面角而变化. 当二面角接近90°时, 分子具有平面或近平面构象时强的铁磁或反铁磁耦合单元, 由于具有近简并的高自旋和低自旋基态, 而变成弱的反铁磁或铁磁耦合单元. 由此提出为获得具有稳定高自旋基态的高自旋分子, 实验上应尽量避免选用强烈扭曲的邻、对苯分子构象.     

以喹啉为耦合单元双自由基体系的理论研究

设计用4种自由基自旋中心连接在耦合单元喹啉的不同位置上的双自由基体系,用AM1-CI方法计算的结果表明:双自由基连接的位置不同对体系耦合作用的影响符合双自由基之间磁性耦合的拓扑规则,即共轭体系中,两个自由基之间以偶数个C(或N)原子耦合,体系具有低自旋基态,表现反铁磁耦合;两个自由基之间以奇数个C(或N)原子耦合,体系具有高自旋基态,表现铁磁耦合.当双自由基连接在喹啉的相邻奇数个C或N原子位置时,体系具有高自旋基态,表现铁磁耦合.     

对嘧啶的邻、间、对位双自由基体系的自旋耦合规律的研究

采用量子化学abinitio法和密度泛函方法对不同取代位置的嘧啶自旋耦合规律进行研究 .两种方法比较 ,用UHF方法计算导致自旋污染严重 ,而用UB3LYP方法计算 ,自旋污染则减少了许多 .计算结果得到了双自由基之间磁性耦合的拓扑规则 :共轭体系中 ,两个自由基之间以偶数个碳 (或氮 )原子耦合 ,则有效交换积分Jij<0 ,体系具有低自旋基态 ,表现为反铁磁耦合 ;两个自由基之间以奇数个碳 (或氮 )原子耦合 ,则有效交换积分Jij>0 ,体系具有高自旋基态 ,表现为铁磁耦合 .自由基性质和铁磁耦合单元的不同位置对自旋耦合的影响较大 ,这些结论为有机磁性材料的分子设计与实验合成提供了理论依据 .     

以-•N-N-作自旋中心杂环作端基的高自旋有机分子的理论研究

本文设计了三种新颖的具有不同排列方式的高自旋有机分子，它们是-•N-N-作自旋中心、苯作铁磁耦合单元、杂环（如苯、吡啶、哒嗪、嘧啶、吡嗪、三嗪）作端基得到的。通过AM1-CI方法计算可知不同的端基对体系的高自旋基态的稳定性影响不同，同时发现由苯(FC)、杂环(EG)和-•N-N- (SC)组成的双自由基体系由于自旋密度在双中心的部分离域，导致-•N-N-自由基的特殊稳定性.而且由于SC三种不同的排列方式，导致其三重态的稳定性随主要SC(-N-)原子间距离的增大而降低。从三种体系的三重态稳定性顺序可知以杂环作端基的三重态的稳定性高于以苯作端基的三重态的稳定性，而且在杂环中，以三嗪作端基的三重态的稳定性最高，嘧啶和吡嗪次之，吡啶和哒嗪最低。     

以杂环为铁磁耦合单元的双中心双自由基高自旋有机分子的理论设计

利用－·N－S－为自旋中心（SC），苯为端基（EG），苯、吡啶、哒嗪、嘧啶 、吡嗪、三嗪为耦合单元（FC），设计三种不同排列方式的新型稳定高自旋分子。 由于自旋密度在杂环（FC）和－·N－S－（SC）组成的体系中自由基双中心的部分 离域，导致－·N-S－自由基的特殊稳定性。三种不同的排列方式中，其三重态的 稳定性随主要SC（－N－）原子间距离的增大而降低。从三个系统八个体系三重态 的稳定性来看，FC上的杂原子位于取代基的间位能提高体系的铁磁耦合作用，而位 于邻位和对位则不利于铁磁耦合作用。     

邻、间、对-二甲氧基亚甲基苯及衍生物双自由基体系的自旋耦 合规律

采用量子化学abinitio法对具有甲氧基的碳、氧双自由基邻、间、对二甲氧基亚甲基苯及衍生物体系基态自旋耦合规律进行研究,得到非平面共轭体系中自由基之间磁性耦合的拓朴规则:共轭体系中,两个自由基之间以偶数个碳原子耦合,则有效交换积分J~i~j<0,体系具有低自旋基态;两个自由基之间以奇数个碳原子耦合,则J~i~j>0,体系具有高自旋基态。自由基性质对自旋耦合的影响较大,正离子自由基间磁性耦合能力较强,这些结论为有机磁性材料的分子设计与实验合成提供了理论依据。     

TTF+1-π共轭桥-6-氧四联氮阳离子双自由基NLO性质的理论研究

采用量子化学UPBE1PBE结合有限场(FF)方法,对系列TTF+?-π共轭桥-6-氧四联氮阳离子双自由基体系的稳定性,极化率αs和第一超极化率βtot进行研究.结果表明,TTF+?-6-氧四联氮阳离子双自由基引入共轭桥后,随体系共轭性增强,αs和βtot值均增大(体系2S的βtot值除外).自旋多重度和构象对双自由基体系的极化率和第一超极化率都有影响,双自由基体系由单重态转变为三重态时,极化率减小,而第一超极化率明显增大.以体系2为例,在单重态时αs和βtot值随构象变化较小,而三重态时αs和βtot值随二面角θ1和θ2的增加而减小.     

异自旋中心三自由基分子的磁性设计

采用量子化学DFT UB3LYP方法,在6-31G(d)基组水平上讨论了以1,3,5-苯为 铁磁耦合单元的异自旋中心三自由基体系。结果表明：对于异自旋三自由基体系, 形成异自旋后对称性降低使部分占据的近简并轨道能级劈裂值增大,反铁磁耦合作 用普遍增强。由两个同自旋和一个异自旋中心构成的三自由基若保持较好的对称性 (C_(2v))可望表现铁磁耦合。在通常情况下这类分子服从已有规律：ΔE_(POMO)小 ,ΔE_(D-Q)大,表现为稳定的四重态基态,且其稳定性与苯环2,4,6碳位置及自 由基原子上的自旋密度相关。     

双卡宾、双氮宾体系基态自旋的从头算研究

采用量子化学从头算UHF方法, 对平面型双卡宾及双氮宾体系的基态自旋情况进行研究。结合前面的分析结果, 进一步探讨了多自由基体系基态自旋的耦合规律, 为有机磁性体的分子设计提供了可靠的理论依据。     

2,2'-联吡啶同自旋双自由基的量子化学研究

采用量子化学abinitio方法,讨论了2,2'-联吡啶同自旋双自由基体系构象变化对铁磁耦合的影响。结果表明,在各种构象下,体系的磁性耦合符合自旋极化规则;对于·CH2,·MH2^+两种自由基磁性耦合性质是机同的,只影响到体系磁性耦合的强度,这一结论为有机磁性材料的分子设计提供了有益的信息。     

Studies on Effect of Conformation on Spin Multiplicity of m－Phenylene Type of Biradicals by Different Approaches

The effect of the conformation on the spin multiplicity of the ground state and the stability of the ground state were investigated for m-phenylene type of biradicals by means of comparative study with DFT, CASSCF and AM1-CI approaches. It was found that AM1-CI approach is reliable in dealing with the stability of the high-spin ground state with the change of conformation; DFT method can give the reasonable results of the spin density of the high-spin state. Furthermore, when one or two radical centers are twisted sufficiently out of the coniugation with the benzene ring, m-phenylene turns into weak ferromagnetic and weak antiferro-magnetic coupling units, respectively.     

Effect of configuration and conformation on the spin multiplicity in xylylene type biradicals

The singlet-triplet splitting energy gap ΔEs-T= Es - ET is calculated for the ortho-, meta-, and para-xylylenes and their heteroatomic analogous by means of AM1-CI approach. It is shownthat when the radical centers R(R=H2C- ,H2N - or HN'- ) are twisted sufficiently far out ofconjugation with the benzene ring, ΔEs-T tends to zero or is negative, i.e. ortho-, meta-, and para-phenylenes turn into weak ferromagnetic or antiferromagnetic coupling unit, while they are strong ferromagnetic (meta-isomers) or antiferromagnetic (ortho-, para-isomers) coupling units under planar conformation. It is suggested that serious twisted conformation is not recommended candidate for the design of novel high-spin molecules with stable high-spin ground states by ortho-or para-phenylene coupling unit.     

卤代氰基卡宾自由基XCCN(X=F,Cl,Br)的理论研究

在C_s对称性和ANO-S基组下, 使用全活化空间自洽场方法(CASSCF), 研究了卤代氰基卡宾自由基及其阴离子的低能电子激发态性质. 为了进一步考虑电子的动态相关效应,采用多组态二级微扰理论(CASPT2)获得更加精确的能量值. 计算结果表明, XCCN的基态是三重态. 单重态和三重态的能隙差ΔE_S-T(kJ/mol): 7.4(FCCN)<13.4(ClCCN)<16.6(BrCCN). 计算得到, XCCN(X=F, Cl, Br)最低垂直激发能分别为408.3, 385.4和 345.2 kJ/mol, 这归因于π(a′) →n_xy 的电子跃迁; XCCN的电子亲和势分别为235.7, 233.0和 237.2 kJ/mol, 与HCCN相比, 其电子亲和势变大.     

Ab initio diradical/zwitterionic polarizabilities and hyperpolarizabilities in twisted double bonds

We investigated the nature of the ground state and static response properties (mu, alpha, and beta) for a promising class of twisted pi-electron system nonlinear optical chromophores at the HF, B3LYP, MP2, and CASSCF levels. We report results for a substituted twisted ethylene and a larger tictoid analogue. Previous work has reported only a zwitterionic character for such tictoid species, however, (14,13) CASSCF calculations predict a ground-state diradical. At the HF, B3LYP, MP2, and (14,13) CASSCF levels, the values of beta are orders of magnitude smaller than those predicted by semiempirical methods.     

Hybrid density functional theory studies on the magnetic interactions and the weak covalent bonding for the phenalenyl radical dimeric pair

The phenalenyl radical (1) is a prototype of the hydrocarbon radical. Recently, the single crystal of 2,5,8-tri-tert-butylphenalenyl (2) was isolated and showed that the two phenalenyl radicals form a staggered dimeric pair, giving rise to strong antiferromagnetic interactions. The origin of the antiferromagnetic interactions and the nature of the chemical bond for the dimeric pair are challenging issues for chemists. First, spin-polarized hybrid DFT (Becke's half and half LYP (UB2LYP)) and CASSCF calculations were performed for 2 and its simplified model, the staggered-stacking phenalenyl radical dimeric pair (3a), to elucidate the origin of the strong antiferromagnetic coupling and the characteristics of the chemical bond. The calculated results showed that a SOMO-SOMO overlap effect was responsible for the strong antiferromagnetic interactions and weak or intermediate covalent bonding between phenalenyl radicals. The tert-butyl groups introduced at three beta-positions hardly affected the magnetic coupling, mainly causing steric hindrances in the crystalline state. Next, to obtain insight into ferromagnetic stacking, we investigated the stacking effect of staggered (3a)- and eclipsed (3b)-stacking phenalenyl radical dimeric pairs with a change of the SOMO-SOMO overlap on the basis of the extended McConnell model. We found that the stacking mode of the dimeric pair with both a small SOMO-SOMO overlap and a ferromagnetic spin polarization effect provided a ferromagnetic coupling.     

para‐Phenylene‐Bridged Spirobi(triarylamine) Dimer with Four Perpendicularly Linked Redox‐Active π Systems

para‐Phenylene‐bridged spirobi(triarylamine) dimer 2 , in which π conjugation through four redox‐active triarylamine subunits is partially segregated by the unique perpendicular conformation, was prepared and characterized by structural, electrochemical, and spectroscopic methods. Quantum chemical calculations (DFT and CASSCF) predicted that the frontier molecular orbitals of 2 are virtually fourfold degenerate, so that the oxidized states of 2 can give intriguing electronic and magnetic properties. In fact, the continuous‐wave ESR spectroscopy of radical cation 2 ^.+ showed that the unpaired electron was trapped in the inner two redox‐active dianisylamine subunits, and moreover was fully delocalized over them. Magnetic susceptibility measurements and pulsed ESR spectroscopy of the isolated salts of 2 , which can be prepared by treatment with SbCl₅, revealed that the generated tetracation 2 ⁴⁺ decomposed mainly into a mixture of 1) a decomposed tetra(radical cation) consisting of a tri(radical cation) moiety and a trianisylamine radical cation moiety (≈75 %) and 2) a diamagnetic quinoid dication in a tetraanisyl‐p‐phenylendiamine moiety and two trianisylamine radical cation moieties (≈25 %). Furthermore, the spin‐quartet state of the tri(radical cation) moiety in the decomposed tetra(radical cation) was found to be in the ground state lying 30 cal mol^?1 below the competing spin‐doublet state.     

N-Ethynylation of Anilides Decreases the Double-Bond Character of Amide Bond while Retaining trans-Conformation and Planarity

Activated amide bonds have been attracting intense attention; however, most of the studied moieties have twisted amide character. To add a new strategy to activate amide bonds while maintaining its planarity, we envisioned the introduction of an alkynyl group on the amide nitrogen to disrupt amide resonance by nN→C_sp conjugation. In this context, the conformations and properties of N-ethynyl-substituted aromatic amides were investigated by DFT calculations, crystallography, and NMR spectroscopic analysis. In contrast to the cis conformational preference of N-ethyl- and vinyl-substituted acetanilides, N-ethynyl-substituted acetanilide favors the trans conformation in the crystal and in solution. It also has a decreased double bond character of the C(O)−N bond, without twisting of the amide. N-Ethynyl-substituted acetanilides undergo selective C(O)−N bond or N−C(sp) bond cleavage reactions and have potential applications as activated amides for coupling reactions or easily cleavable tethers.     

Ab initio and density function theory computational studies of the CH4 + H → CH3 + H2 reaction

The activation barrier for the CH₄ + H → CH₃ + H₂ reaction was evaluated with traditional ab initio and Density Functional Theory (DFT) methods. None of the applied ab initio and DFT methods was able to reproduce the experimental activation barrier of 11.0-12.0 kcal/mol. All ab initio methods (HF, MP2, MP3, MP4, QCISD, QCISD(T), G1, G2, and G2MP2) overestimated the activation energy. The best results were obtained with the G2 and G2MP2 ab initio computational approaches. The zero-point corrected energy was 14.4 kcal mol⁻¹. Some of the exchange DFT methods (HFB) computed energies which were similar to the highly accurate ab initio methods, while the B3LYP hybrid DFT methods underestimated the activation barrier by 3 kcal mol⁻¹. Gradient-corrected DFT methods underestimated the barrier even more. The gradient-corrected DFT method that incorporated the PW91 correlational functional even generated a negative reaction barrier. The suitability of some computational methods for accurately predicting the potential energy surface for this hydrogen radical abstraction reaction was discussed.     

Effect of histidine on the sensitized generation of singlet oxygen in the complexes with chlorophyll

The energies of the intermolecular interactions of an O₂ molecule in the ground and excited states with the electron-excited and non-excited model complexes of chlorophyll were calculated using the DFT, CASSCF, SA-CASSCF, MCQDPT2, and XMCQDPT2 methods. The activation energies of formation and dissociation of the oxygen complexes were estimated. The radiative electric dipole moments of (0 → 0) spin-allowed S → S, T → T, and spin-forbidden S → T transitions were calculated taking into account the spin-orbit coupling, and rate constants of nonradiative transitions that determine the generation and deactivation of the O₂ molecule (¹Δ_g) were evaluated. The effect of histidine on the probability of singlet oxygen generation sensitized by the model chlorophyll complex was considered in detail.     

Exploring the potential energy surface of retinal, a comparison of the performance of different methods

The ground state structure of retinal has been investigated. We found that DFT and CASSCF produce different results for the bond length alternation in a model system of retinal. Quantum mechanics/molecular mechanics calculations including the closest surrounding amino acids have been performed, using DFT and CASSCF to calculate the structure of retinal in the protein cavity. The planarity of the retinal molecule is affected by the surrounding protein. DFT and CASSCF produce different twist angles. The difference between CASSCF and DFT appears to be related to the positively charged nitrogen of the Schiff base, which leads to different pi-bond orders produced by the two methods.