(CH2)2N和(CH3)2NH+的密度泛函理论计算

用密度泛函理论UB3LYP／6-31G（d,p）方法研究了二甲胺自由基（CH3）2N及其质子化离子（CH3）2NH 的构象和超精细结构．在由构象研究得到的两种自由基的最稳定结构上,用密度泛函的UB3LYP和UB3PW91方法及从头计算的UHF、UMP2（FULL）方法计算了α-质子、β-质子和N核上的超精细偶合常数A（Hα）、A（Hβ）和A（N）结果表明：两种自由基中甲基内旋转的位垒均很低,分别为0．46kJ·mol-1（（CH3）2NH ）和2．6kJ·mol-1（（CH3）2N）．UB3LYP／6-31G（d,p）和UB3PW91／6-31G（d,p）计算的A（Hα）、A（Hβ）和A（N）与ESR实验结果符合得很好,UMP2／6-31G（d,p）方法的计算值与实验值符合得也较好．     

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.     

Quantum mechanical investigation of the inner-sphere reorganization energy of cyclooctatetraene/cyclooctatetraene radical anion. Part I

The inner-sphere reorganization energy of the electron self-exchange of the couple cyclooctatetraene/cyclooctatetraene radical anion has been investigated by quantum mechanical calculations. The more stable Jahn Teller distorted B2g conformation of the radical anion has been used in this study. Two different theories have been applied in this first part. The harmonic approximation in the classical Marcus scheme has been modified by using projected force constants, which are obtained from the complete force constant matrix and the geometry changes of the molecule during the ET (introduced by Mikkelsen). A different approach (introduced by Nelsen) combines the different energies of the neutral and radical anion with and without relaxation corresponding to the vertical ionization potential and the vertical electron affinity. The electronic energies of the neutral molecule and the radical anion differ dramatically applying three different levels of quantum mechanical calculations (UAM1, UB3LYP, PMP2 with three different basis sets with and without diffuse functions). Nevertheless the Nelsen method gives almost consistent results for the inner-sphere reorganization energies: 120.1 kJ/mol for semiempirical UAM1 method, 159.3 kJ/mol, 156.4 kJ/mol and 158.3 kJ/mol for density functional UB3LYP/6-31G*, UB3LYP/6-31++G* and UB3LYP/AUG-cc-pVDZ calculations and 192.5 kJ/mol for ab-initio PMP2/6-31G* investigations, respectively. These values are in agreement with earlier experimental work supposing the total reorganization energy to be larger than 38 kcal/mol assuming an electron self-exchange rate of 10(4) M(-1) s(-1). The simple harmonic approximation of Marcus relation has not yet been applied for a molecule like cyclooctatetraene with large torsional geometry changes. Using the projected force constants after scaling, considerably different results for the inner-sphere reorganization energy have been calculated: 738.1 kJ/mol for the UB3LYP/6-31G*, 743.3 kJ/mol for UB3LYP/6-31++G* and 759.1 kJ/mol for UB3LYP/AUG-cc-pVDZ level of theory. Comparison with our concentration dependent EPR experiments are controversial to the earlier experimental results, but the latter supports the assumption that the electron self-exchange occurs in a time scale so that the molecules cannot complete their vibrational motions. Therefore the projected Marcus relation is not valid for cyclooctatetraene/cyclooctatetraene radical anion including a large torsional change during the electron transfer.     

Inductive effects in radicals calculated from DFT energies; substituted bicyclo[2.2.2]octan-1-yloxy radicals

Energies of a series of 4-substituted 1-oxybicyclo[2.2.2]octan-1-yloxy radicals with 18 various substituents were calculated within the framework of the DFT theory at the levels UB3LYP/6-311+G(d,p)//UB3LYP/6-311+G(d,p) and UB3LYP/6-311++G(2df,p)//UB3LYP/6-311+G(d,p) and compared with similar series of the parent alcohols, their deprotonated and protonated forms calculated at the levels B3LYP/6-311+G(d,p)//B3LYP/6-311+G(d,p) and B3LYP/6-311++G(2df,p)//B3LYP/6-311+G(d,p). The two levels are of the same performance and both are sufficient for molecules of this type according to comparison with scarce experimental gas-phase acidities and basicities. The substituent effects were analyzed in terms of isodesmic equations. In addition to strong dependence on the substituent inductive effect, a slight dependence on the electronegativity of the first atom of the substituent was proven in certain cases. In all aspects, there is no qualitative difference between the effects on radicals and on similar closed shell species. Radicals behave as slightly electron deficient; the substituent effect is weaker than that on the ions but stronger than on neutral molecules.     

N2O在焦炭表面异相生成和分解机理的密度泛函理论研究

选用合理简化的焦炭模型,对煤焦燃烧过程中N₂O的异相生成和分解机理进行了分子水平上的研究。采用UB3LYP/6-31G(d)密度泛函理论方法优化得到了反应路径上反应物、产物、中间体和过渡态的几何构型和各中间反应的活化能和反应焓变。NO与其预先吸附在焦炭表面解离生成的表面氮组分反应生成N₂O的路径有两个,需要克服的势垒分别为69.3kJ/mol和200.0kJ/mol;NO亦可直接与焦炭中的吡啶氮结合释放出N₂O,该反应路径所需克服的最大势垒为418.0kJ/mol。N₂O可在焦炭表面分解释放出N₂,异相分解反应为一步反应,计算所得活化能为100.8kJ/mol。N₂O的异相生成和异相分解反应均为放热反应。采用经典过渡态理论计算得到了各路径中速率控制步骤的反应速率常数。低温条件下,N₂O的异相分解反应速率略低于其异相生成速率,随着温度的升高,两者逐渐接近,说明高温条件有利于N₂O的异相分解。     

Theoretical calculations on the cycloreversion of oxetane radical cations

The molecular mechanism for the cycloreversion of oxetane radical cations has been studied at the UB3LYP/6-31G* level. Calculations support that the cycloreversion takes place via a concerted but asynchronous process, where C-C bond breaking at the transition state is more advanced than O-C breaking. This allows a favorable rearrangement of the spin electron density from the oxetane radical cation (with the spin density located mainly on the oxygen atom) to the alkene radical cation which is one of the final products. Inclusion of solvent effects does not modify the gas-phase results.     

12]Annulene radical anions revisited: evaluation of structure assignments based on computed energetic and electron spin resonance data

We report density functional and coupled cluster calculations on numerous monocyclic and bicyclic (CH)12(*-) isomers. At the RCCSD(T)/cc-pVDZ//UB3LYP/6-31+G* level, a nearly planar, bond-equalized radical anion of 1,7-di-trans-[12]annulene (4a(*-)) is lowest in energy; several other isomers and conformations lie within 3 kcal/mol of 4a(*-). RCCSD(T)/AUG-cc-pVDZ//UB3LYP/6-31+G* results place the all-cis isomer 3(*-) slightly below 4a(*-) in energy. Validation studies on the heptalene radical anion, [16]annulene radical anion, and tri-trans-[12]annulene radical anion indicate that electron spin resonance (ESR) hyperfine coupling constants (aH values) computed at the BLYP/EPR-III level on DFT geometries give much better agreement with experimental values than those computed using B3LYP/6-31G*. We were unable to locate any C12H12(*-) isomer that could account for the ESR spectrum previously attributed to a highly twisted structure for the 1,7-di-trans-[12]annulene radical anion. Our computed energetic and ESR data for [12]annulene radical anions and their valence isomers suggest that 4a(*-) may have been made, yet its ESR spectrum was incorrectly assigned to the bicyclic isomer 6b(*-). Finally, the computed (1)H NMR shift values of the dianion of 4 reveal a distinct diatropic ring current that should aid in its characterization.     

Theoretical elucidation of kinetic and thermodynamic control of radical addition regioselectivity

The cyclizations of two structurally similar 2-oxo-5-hexenyl-type radicals have been investigated by ab initio and density functional (UB3LYP/6-31+G**//UHF/6-31G* and UB3LYP/6-31G*//UB3LYP/6-31G*) calculations. The origin of apparently contradictory reports of 6-endo and 5-exo cyclizations is determined. Kinetic control favors 6-endo cyclization, while thermodynamic control gives 5-exo cyclization, and the observation of different products from different research groups arises from the difference in experimental conditions used by the two groups. The outcome of a new cyclization reaction was predicted by using these theoretical techniques. Kinetic control is predicted to yield exclusively the products of 6-endo cyclization, while thermodynamic control would lead to an approximately equal mixture of one 6-endo and one 5-exo cyclized product. Experimental studies revealed that the reaction yields only the products of 6-endo cyclization through kinetic control.     

CH2FCF3与O(1D)反应机理的理论研究

用量子化学密度泛函理论(DFT)和G3(MP2)B3方法,对O(1↑D)与CH2FCF3的反应进行了研究．在UB3LYP/6-31G(d)计算水平上,优化了反应势能面上各驻点的几何结构,在G3(MP2)B3水平上进行了单点计算,并利用UB3LYP/6-311 G(3df,3pd)计算的波函数进行了电荷密度分析．通过内禀坐标(IRC)计算和振动分析,对反应过渡态进行了确认,并确定了反应机理．     

C_3O_2分子结构和光谱的密度泛函理论研究

使用密度泛函理论,在B3LYP/6-31G(d)和B3LYP/6-311G(2d)水平上,研究了C_3O_2分子的可能几何构型,并在6-31G(d)水平上计算了其中2种总能量最小的构型的振动频率,同时与实验观察值进行了比较, 计算结果当C_3O_2分子具有C2v对称性的W型弯曲结构（键角C－C－C和C－C－O分别为162.3°和178.8°）时,振动频率的计算值和实验观察值非常吻合。     

戊烯自由基阳离子的密度泛函理论研究

使用密度泛函理论B3LYP方法和6-31G(d,p)、6-31＋G(d,p)、6-311G(d,p)及6-311＋G(d,p)基组，分别对2-C5H10＋和1-C5H10＋的各种构象进行了几何构型优化，并用B3LYP/6-311G(d,p)进行了频率分析计算.计算预言1-C5H10＋具有非平面构型，与以往报导的从头算计算结论相反.在两个自由基阳离子的各种构象的B3LYP几何构型上，进行了B3LYP和UMP2(full)方法的超精细偶合常数计算，得到了比以往更好的结果.     

咪唑氧自由基分子NLO性质的密度泛函理论研究

采用密度泛函理论(DFT) UB3LYP方法, 在6-31g(d)水平上对2,2’-(1,2-乙炔基-4,1-亚苯基)双[4,4’,5,5’-四氢]咪唑氧自由基分子及其异构体的自旋耦合性质进行分析, 并结合有限场(FF)方法计算它们的非线性光学(NLO)系数, 以探讨咪唑氧环在共轭链不同位置时体系的自旋耦合规律和NLO系数. 结果表明, 所有体系基态自旋符合自旋极化规则, 它们的极化率随自旋多重度的增加而减小; 一阶超极化率因受分子对称性影响, 对称性不同其一阶超极化率的变化也不同; 二阶超极化率呈现随着自旋多重度的增加而增加的趋势. 从理论上探讨这些自由基分子自旋耦合规律与NLO活性的关系, 为有机自由基NLO材料的分子设计与实验研究提供一定的理论依据.     

NMR spectra, GIAO and charge density calculations of five-membered aromatic heterocycles

The B3LYP/6-31+G(d) molecular geometry optimized structures of 17 five-membered heterocycles were employed together with the gauge including atomic orbitals (GIAO) density functional theory (DFT) method at the B3LYP/6-31+G(d,p), B3LYP/6-311++G(d,p) and B3LYP/6-311+G(2d,p) levels of theory for the calculation of proton and carbon chemicals shifts and coupling constants. The method of geometry optimization for pyrrole (1), N-methylpyrrole (2) and thiophene (7) using the larger 6-311++G(d,p) basis sets at the B3LYP/6-31+G(d,p), B3LYP/6-311++G(d,p), B3LYP/6-31+G(2d,p) and B3LYP/cc-pVTZ levels of theory gave little difference between calculated and experimental values of coupling constants. In general, the (1)H and 13C chemical shifts for all compounds are in good agreement with theoretical calculations using the smaller 6-31 basis set. The values of nJHH(n=3, 4, 5) and rmnJ(CH)(n=1, 2, 3, 4) were predicted well using the larger 6-31+G(d,p) and 6-311++G(d,p) basis sets and at the B3LYP/6-31+G(d,p), B3LYP/6-311++G(d,p), B3LYP/6-31+G(2d,2p) levels of theory. The computed atomic charges [Mülliken; Natural Bond Orbital Analysis (NBO); Merz-Kollman (MK); CHELP and CHELPG] for the B3LYP/6-311++G(d,p) geometry optimized structures of 1-17 were used to explore correlations with the experimental proton and carbon chemical shifts.     

Quantum-Chemical Modeling of B<Subscript>32</Subscript> Complexes with Nitrogen: Endo or Exo?

Computational modeling of borospherene B₃₂ complexes with nitrogen has been performed by the density functional theory method (DFT UB3LYP*/6-311+G(d)). The geometry of isomers of these complexes with a different starting orientation of N has been found, their relative energies have been calculated, and the spin density distribution has been studied. Computation results predict that the structures with a nitrogen atom or ion built into the boron cage are more favorable than the endohedral isomers.     

Quantum chemistry study on the mechanism of the reaction between ozone and 2,3,7,8‐TCDD

Despite of its fundamental importance, the mechanism of the reaction between ozone and dioxins are still lack detailed investigation so far. It is well-known that quantum chemical calculation is a well-established method for investigating chemical reactions. In this article, quantum chemical calculation was employed to investigate the mechanism of the reaction between ozone and dioxins, as exemplified by 2,3,7,8-TCDD. The theoretical study showed that, 2,3,7,8-TCDD was gradually destructed by ozone via six cleavages of the CC bonds. All the six cleavages of the CC bonds were calculated and discussed in detail based on the theoretical calculations by the UB3LYP/6-31G(d) method. At the same time, the energies of stationary points along the reaction process were calculated by the UMP2/6-311g(d,p)//UB3LYP/6-31G(d) method and the activation energy was obtained. The obtained activation energy was 12.25 kcal/mol, which was lower than that of the reaction between benzene and O₃(16.64 kcal/mol). This indicated that, by comparison with benzene, 2,3,7,8-TCDD could be more efficiently destructed by O₃. The reason for this result was also discussed. © 2010 Wiley Periodicals, Inc. Int J Quantum Chem, 2011     

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.     

Si~2Cl~6分子的内旋转能垒和振动基频的理论研究

用从头算方法HF/6-31G^*^*和密度函方法B3LYP/6-31G^*^*,对Si~2Cl~6分子的平衡几何构型进行优化,优化的结果与实验结果吻合得较好.并用上述两种不同的方法计算Si~2Cl~6分子的内旋转能垒,结果分别为8.786和6.694kJ/mol,其中DFT方法的计算结果与实验结果4.18kJ/mol吻合得较好.对Si~2Cl~6分子的振动基频进行计算.用HF/6-31G^*^*SQM力场所计算的频率理论值与实验值的平均误差为7.3cm^-^1,用B3LYP/6-31G^*^*未标度的力场所计算的频率理论值与实验值的平均误差为6.0cm^-^1.该密度泛函方法(B3LYP/~6-31G^*^*)的理论计算值比用HF/6-31G^*^*标度后的SQM力场计算的频率与实验值(除Si--Si键扭转振动基频之外的11条振动基频)吻合得更好.并给出了Si--Si键扭转振动基频的预测值。     

HO3、HO+3、HO-3光谱和结构的MP2和DFT理论研究

用B3LYP和MP2量子化学理论,详细地计算了HO3、HO+3、HO-3可能的平衡构型、能量和光谱参数,分析了它们可能的解离通道和稳定性.还分析了HO3的电离势(IE)及电子亲和势(Ea),讨论了HO3、HO3-自由基的光电子能谱的特征,及用光电子能谱去探测HO3的可能性.     

Vibrational spectroscopy investigation using ab initio and density functional theory analysis on the structure of 3-aminobenzotrifluoride

In this work, we will report a combined experimental and theoretical study on molecular and vibrational structure of 3-aminobenzotrifluoride. The FT-Raman and Fourier transform infrared spectra of 3-aminobenzotrifluoride (3ABTF) were recorded in the liquid phase. The equilibrium geometry, harmonic vibrational frequencies, infrared intensities and Raman scattering activities, depolarization ratios, reduced masses were calculated by HF and density functional B3LYP method with the 6-31G(d,p) and 6-311G(d,p) basis sets. The scaled theoretical wavenumbers showed very good agreement with the experimental values. The thermodynamic functions of the title compound were also performed at HF/6-31G(d,p)/6-311G(d,p) and B3LYP/6-31G(d,p)/6-311G(d,p) levels of theory. A detailed interpretations of the infrared and Raman spectra of 3ABTF is reported. The theoretical spectrograms for FT-IR spectra of the title molecule have been constructed.