槲皮素及其异构体分子的结构性质研究

槲皮素是一类有着重要药理活性的黄酮类化合物。本文在B3LYP/6-31G**水平上对槲皮素及其14种异构体进行了几何优化,得到各种异构体的结构与性质。由电荷布居数分析出B环是此类多羟基化合物的活性部分,3位-OH的存在抑制了A、C环的活性。同时,表明此类化合物中必有分子内氢键形成。并给出了最稳构型、槲皮素和标准结构三类分子的振动光谱。探讨了它们的构效关系。     

Theoretical study on stability and properties of NC2O isomers

The structures, energetics, spectroscopies, and stabilities of the doublet NC(2)O radical are explored at density functional theory and ab initio levels. Nine minimum isomers are located connected by 22 interconversion transition states. At the CCSD(T)/6-311+G(2df)//QCISD/6-311G(d)+ZPVE level, the lowest-lying isomer is bent NCCO 1 (0.0 kcal/mol) with (2)A' state followed by bent isomer CNCO 2 (16.7). Two isomers (1 and 2) and another high-lying species CCNO 4 (99.4) with bent structure are considerably stabilized by a barrier of at least 20 kcal/mol. All of the three isomers should be experimentally or astrophysically observable. This result is consistent with their indication of neutralization-reionization mass spectrometry experiments. Also, the calculated spectroscopic properties and bond distances of known NCCO 1 are consistent with recent experimental observations and theoretical studies. The bonding natures of the isomers 1, 2, and 4 are analyzed. Their molecular properties including the heats of formation, adiabatic ionization potentials, and adiabatic electronic affinities are calculated at the higher levels G3//B3LYP, G3(MP2)//B3LYP, QCISD, and CCSD(T) (single-point). Possible formation strategies of the isomers 1, 2, and 4 in laboratory and space are also discussed in detail.     

TBHQ及其衍生物结构和热力学性质的理论研究

采用Materials Studio软件中的DMol3程序对特丁基对苯二酚(TBHQ)及衍生物二特丁基对苯二酚(DTBHQ)的几何结构和性质(振动频率、反应活性及热力学性质)进行了理论研究,得到了分子的稳定几何构型、各原子上的电荷分布、热力学性质、Fukui指数和前线分子轨道.计算结果表明:TBHQ与DTBHQ分子结构中的2个—OH是影响其抗氧化活性的主要部位和亲电反应发生的活性位点,整体表现为比较理想的稳定抗氧化性能.     

Theoretical study on structures and stability of Si2CP isomers

The structures, energetics, spectroscopies, and isomerization of various doublet Si2CP species are explored theoretically. In contrast to the previously studied SiC2N and SiC2P radicals that have linear SiCCN and SiCCP ground states, the title Si2CP radical has a four-membered-ring form cSiSiPC 1 (0.0 kcal/mol) with Si-C cross-bonding as the ground-state isomer at the CCSD(T)/6-311G(2df)//B3LYP/6-311G(d)+ZPVE level, similar to the Si2CN radical. The second low-lying isomer 2 at 11.6 kcal/mol has a SiCSiP four-membered ring with C-P cross-bonding, yet it is kinetically quite unstable toward conversion to 1 with a barrier of 3.5 kcal/mol. In addition, three cyclic species with divalent carbene character, i.e., cSiSiCP 7, 7' with C-P cross-bonding and cSiCSiP 8 with Si-Si cross-bonding, are found to possess considerable kinetic stability, although they are energetically high lying at 44.4, 46.5, and 41.4 kcal/mol, respectively. Moreover, a linear isomer SiCSiP 5 at 44.3 kcal/mol also has considerable kinetic stability and predominantly features the interesting cumulenic /Si=C=Si=P/* form with a slight contribution from the silicon-phosphorus triply bonded form /Si=C*-Si[triple bond]P/. The silicon-carbon triply bonded form *Si[triple bond]C-Si[triple bond]P/ has negligible contribution. All five isomers are expected to be observable in low-temperature environments. Their bonding nature and possible formation strategies are discussed. For relevant species, the QCISD/6-311G(d) and CCSD(T)/6-311+G(2df) (single-point) calculations are performed to provide more reliable results. The calculated results are compared to those of the analogous C3N, C3P, SiC2N, and Si2CN radicals with 17 valence electrons. Implications in interstellar space and P-doped SiC vaporization processes are also discussed.     

Theoretical study on structures and stability of C4P isomers

The structures, energetics, spectroscopies, and stabilities of doublet C(4)P isomeric species are explored at the DFT/B3LYP, QCISD, and CCSD(T) (singlet-point) levels. A total of 12 minimum isomers and 27 interconversion transition states are located. At the CCSD(T)/6-311G(2df)//QCISD/6-311G(d)+ZPVE level, the lowest-lying isomer is a floppy CCCCP 1 (0.0 kcal/mol) mainly featuring a cumulenic structure |C=C=C=C=P*|, which differs much from the analogous C4N radical (|*C-C[triple bond]C-C[triple bond]N|). The quasi-linearity and the low bending mode of 1 are in contrast to the previous prediction. The second energetically followed isomer PC-cCCC 3 (14.9 kcal/mol) possesses a CCC ring-bonded to CP. The two low-lying isomers are separated by a high-energy ring-closure/open transition state (26.5 kcal/mol) and thus are very promising candidates for future laboratory and astrophysical detection. Furthermore, four high-energy isomers, that is, two bent isomers CCPCC 2 (68.4 kcal/mol) and CCPCC 2' (68.5 kcal/mol) and two cagelike species 10 (56.0 kcal/mol) and 11 (67.9 kcal/mol), are also stabilized by considerable barriers. The present work is the first detailed potential energy survey of CnP clusters and can provide useful information for the investigation of larger CnP radicals and for understanding the isomerism of P-doped C vaporization processes.     

三重态SiCP2异构体的结构和稳定性的理论研究

采用DFT,QCISD和CCSD(T)等理论计算方法对三重态SiCP2异构体的结构和稳定性进行了理论研究.在B3LYP/6-311G(d)水平下,共计算得到由17个过渡态相连接的15个异构体.在CCSD(T)/6-311 +G(2df)//QCISD/6-311G(d)水平下,考虑重点振动能相对能量最低的三元环状异构体P-cCSiP 8(0.0 kJ/mol)及四元环状结构的cPCSiP 4具有相当大的动力学稳定性,在一定的实验室和星际条件下可能被检测到.另外,对它们的成键性质也进行了分析.     

Theoretical study on structures and stability of HC2P isomers

 The structures and isomerization pathways of various HC₂P isomers in both singlet and triplet states are investigated at the B3LYP/6-311G(d,p), QCISD/6-311G(d,p) (for isomers only) and single-point CCSD(T)/6-311G(d,p)//B3LYP/6-311G(d,p) levels. At the CCSD(T)/6-311G(d,p)//B3LYP/6-311G(d,p) level, the lowest-lying isomer is a linear HCCP structure ³ 1 in the ³∑⁻ state. The second low-lying isomer has a CPC ring with exocyclic CH bonding ¹ 5 in a singlet state at 10.5 kcal/mol. The following third and fourth low-lying isomers are a singlet bent HCCP structure ¹ 1 at 20.9 kcal/mol and a bent singlet HPCC structure ¹ 3 at 35.8 kcal/mol, respectively. Investigation of the HC₂P potential-energy surface indicates that in addition to the experimentally known isomer ³ 1, the other isomers ¹ 1, ¹ 3 and ¹ 5 also have considerable kinetic stability and may thus be observable. However, the singlet and triplet bent isomers HCPC ¹ 2 and ³ 2 as well as the triplet bent isomer HPCC ³ 3 are not only high-lying but are also kinetically unstable, in sharp contrast to the situation of the analogous HCNC and HNCC species that are both kinetically stable and that have been observed experimentally. Furthermore, the reactivity of various HC₂P isomers towards oxygen atoms is briefly discussed. The results presented here may be useful for future identification of the completely unknown yet kinetically stable HC₂P isomers ¹ 1, ¹ 3 and ¹ 5 either in the laboratory or in interstellar space. Received: 5 November 2000 / Accepted: 25 November 2001 / Published online: 8 April 2002     

A quantum-chemical study of the thermodynamic properties of the CNOH<Subscript>3</Subscript> isomers

The B3LYP/6-311++G(3df,3pd) DFT method was used to determine the geometry and vibrational frequencies for the formamide, formaldoxime, nitrosomethane, oxaziridine, and formimide. The potential energy functions of the hindered internal rotations were calculated. The conformers of formaldoxime and formimide were determined. For all the molecules, including the conformers, the thermodynamic characteristics Δ_f H° (298 K), S°(298 K), C _p(298 K), and Δ_f G° (298 K) were determined. The temperature dependences C _p(T) within 298–1500 K were represented in the form of cubic polynomials.     

Theoretical and experimental NMR study of protopine hydrochloride isomers

The 1H and 13C NMR chemical shifts of cis- and trans-protopinium salts were measured and calculated. The calculations of the chemical shifts consisted of conformational analysis, geometry optimization (RHF/6-31G** method) and shielding constants calculations (B3LYP/6-31G** method). Based on the results of the quantum chemical calculations, two sets of experimental chemical shifts were assigned to the particular isomers. According to the experimental results, the trans-isomer is more stable and its population is approximately 68%.     

多溴代芴热力学性质的理论预测研究

为了预测环境污染物多溴代芴的热力学性质,计算了60种多溴代芴的分子连接性指数、电性距离矢量指数、电性拓扑状态指数以及分子形状指数.采用多元线性回归的方法建立了多溴代芴化合物的标准生成焓(Δ_fH~0)、标准熵(S~0)、标准生成吉布斯自由能(Δ_fG~0)及等容热容等热力学性质的QSPR相关性模型,方程的相关系数值均在0.99以上.与相关文献的相关系数R和标准误差S进行比较,该模型具有良好的预测能力和稳定性.     

Theoretical study on the structures, isomerization and stability of SiC4 isomers

The structures, energetics, dipole moments, vibrational spectra, rotational constants, and isomerization of singlet SiC₄ isomers were explored using ab initio methods. Five types of isomers, a total of 11 minima, connected by 11 interconversion transition states, were located on the potential energy surface at the MP2/6-311G(d, p) level. More accurate energies were obtained at the G3(MP2) level. With the highest isomerization barrier, a C2v tetra-angular cone possesses the largest kinetic stability. The lowest-lying structure, linear SiCCCC is also highly kinetically stabilized. Besides, D2d bicyclic c-Si(CC)₂, C2v five-membered ring c-SiCCCC, another C2v tetra-angular cone isomer and C3v trigonal bipyramid isomer are also considered to be kinetically stable, because their isomerization barriers are all over 10 kcal/mol. Other isomers cannot be kinetically stabilized with considerably low isomerization barriers. Investigation on the vibrational spectra, dipole moments, and rotational constants for SiC₄ isomers are valuable for their detections in the interstellar space and laboratory.     

Theoretical study of the structures and properties of cyclic nitramines: Tetranitrotetraazadecalin (TNAD) and its isomers

Density Functional Theory (DFT) was employed to study the geometries, electronic structures, infrared vibrational spectra, and thermodynamic properties of seven isomeric cyclic nitramines of C₆H₁₀N₈O₈ (i.e., TNAD and its six isomers) at the B3LYP/6‐31G* level of theory. The experimental results available for TNAD were used to determine the reliability of the DFT method for generating structural and IR spectroscopic values for these molecular systems. The relative stabilities of the conformers were evaluated from the energy differences of the structures. Detonation properties of various conformers were evaluated using the Kamlet‐Jacobs equations, and it was found that all the calculated results are comparable to the available experimental data. In addition, the calculated results demonstrate that all title compounds can be used as excellent propellant ingredients. © 2005 Wiley Periodicals, Inc. Int J Quantum Chem, 2005     

Theoretical study in gas phase of linear and nonlinear optical properties of the ortho-, meta- and para-nitrophenol isomers

The linear (α), and nonlinear (β, γ) optical NLO properties of ortho-, meta- and para-nitrophenol (ONP, MNP and PNP) isomers have been calculated in gas phase by using ab initio (HF, MP2 and MP4) and density functional theory (DFT) (B3LYP, CAM-B3LYP) methods, with the 6-31+G(d,p) and 6-311+G(3d,3p) standard and the Sadlej specialized basis sets. These properties were evaluated both at static and at dynamic regime within the finite field FF numerical techniques and the time-dependent-Hartree–Fock approach at 1,910 nm, respectively. Additional calculations were performed for the β static hyperpolarizability of these isomers in presence of p-dioxane solvent with the Onsager Model and the SCRF-PCM approach, using B3LYP/6-31+G(d,p) and MP2/6-31+G(d,p) levels of theory. Additionally, CCSD/6-31+G(d,p) calculations were performed for the α, β and γ properties of PNP isomer. The B3LYP and MP2 α _ave results of the nitrophenol isomers are comparable to the experimental α _ave reports; while the tendency for the β _v calculated values (β _v PNP > β _v MNP > β _v ONP), that can be explained in terms of the O _x atomic charge of the –NO₂ group, does not follow exactly the experimental ones. The B3LYP γ _ave results are in correspondence to the experimental measurements, the correlation of which is r ² = 0.99. The use of FF methodology in conjunction with the B3LYP and MP2 methods and the 6-31+G(d,p) basis set show to be appropriate approaches to predict qualitative optical properties of Push–Pull like organic molecules, provided are considered the solvent effects or frequency dependence. However, to have a clear picture of the NLO properties of an isolated molecule, higher order correlation effects combined with specialized basis sets, frequency and solvent effects should be employed. We have demonstrated that MP4/Sadlej level of theory is able to reproduce NLO properties that can be considered equivalent to those from more sophisticated approaches, such as CCSD together with extended basis sets.     

Theoretical study on the structures and aromaticity of N‐confused porphyrazine isomers

The energy and nucleus‐independent chemical shift (NICS) of 95 isomers of N‐confused porphyrazine (NCPz), including normal porphyrazine (N₀CPz), N‐confused porphyrazine (N₁CPz), doubly N‐confused porphyrazine (N₂CPz), triply N‐confused porphyrazine (N₃CPz), and fully N‐confused porphyrazine (N₄CPz), have been calculated by the density‐functional theory (DFT) method. The stability of NCPz decreased by increasing the number of confused pyrrole rings and the macrocycle tend to be destabilized stepwise by approximately +21 kcal/mol. The relative energies of the most stable isomers in confusion level are 0 kcal/mol (N₀CPz‐1), +23.481 (N₁CPz‐5), +41.849 (N₂CPz‐a4), +61.738 (N₃CPz‐b3), and +84.596 (N₄CPz‐b13), respectively. The most stable isomers of N₂CPz, N₃CPz, and N₄CPz are not necessarily the most aromatic but rather nonaromatic, especially in the case of N₃CPz and N₄CPz. On the other hand, the magnitude of the aromaticity estimated by NICS for these isomers does not differ largely. The NICS values of the most aromatic isomers are ?15.5411 (N₀CPz‐1), ?14.0458 (N₁CPz‐2), ?12.8171 (N₂CPz‐d1), ?11.5961 (N₃CPz‐b6), and ?12.8012 ppm (N₄CPz‐a6), respectively. © 2009 Wiley Periodicals, Inc. Int J Quantum Chem, 2010     

Theoretical study on the structures, isomerization and stability of SiC3H isomers

The calculations of the geometry optimizations, energies, dipole moments, vibrational spectra, rotational constants, and isomerization of doublet SiC₃H species were performed using density functional theory and ab initio methods. Four types of isomers, a total of 18 minima, connected by 16 interconversion transition states, were located on the potential energy surface (PES) at the B3LYP/6-311G (d, p) level. More accurate energies were obtained at the CCSD(T)/6-311G(2df, 2p), and G3(MP2) levels. With the highest isomerization barrier, the lowest lying structure, linear A1 possesses the largest kinetic stability. Besides, the isomerization barriers of A2, A4, C2, F1, F4 and F5 are over 10 kcal/mol, and these isomers are also considered to be higher kinetically stable. Other isomers cannot be kinetically stabilized with considerably low isomerization barriers. Investigation on the bonding properties and the computations of vibrational spectra, dipole moments, and rotational constants for SiC₃H isomers are helpful for understanding their structures and also valuable for their detections in the interstellar space and laboratory.     

Theoretical studies on the structures, thermodynamic properties, detonation properties, and pyrolysis mechanisms of spiro nitramines

Density function theory (DFT) has been employed to study the geometric and electronic structures of a series of spiro nitramines at the B3LYP/6-31G level. The calculated results agree reasonably with available experimental data. Thermodynamic properties derived from the infrared spectra on the basis of statistical thermodynamic principles are linearly correlated with the number of nitramine groups as well as the temperature. Detonation performances were evaluated by the Kamlet-Jacobs equations based on the calculated densities and heats of formation. It is found that some compounds with the predicted densities of ca. 1.9 g/cm3, detonation velocities over 9 km/s, and detonation pressures of about 39 GPa (some even over 40 GPa) may be novel potential candidates of high energy density materials (HEDMs). Thermal stability and the pyrolysis mechanism of the title compounds were investigated by calculating the bond dissociation energies (BDE) at the B3LYP/6-31G level and the activation energies (E(a)) with the selected PM3 semiempirical molecular orbital (MO) based on the unrestricted Hartree-Fock model. The relationships between BDE, E(a), and the electronic structures of the spiro nitramines were discussed in detail. Thermal stabilities and decomposition mechanisms of the title compounds derived from the B3LYP/6-31G BDE and the UHF-PM3 E(a) are basically consistent. Considering the thermal stability, TNSHe (tetranitrotetraazaspirohexane), TNSH (tetranitrotetraazaspiroheptane), and TNSO (tetranitrotetraazaspirooctane) are recommended as the preferred candidates of HEDMs. These results may provide basic information for the molecular design of HEDMs.     

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.     

Theoretical investigation on the structural and thermodynamic properties of FeSe at high pressure and high temperature

A theoretical investigation on structural and thermodynamic properties of 11-type iron-based superconductor FeSe at high pressure and high temperature was performed by employing the first-principles method based on the density functional theory. Some structural parameters of FeSe in both tetragonal and hexagonal phases are reported. According to the fourth-order Birch-Murnaghan equation of states, the transition pressure P(t) of FeSe from the PbO-type phase to the NiAs-type phase was determined. The calculated results are found to be in good agreement with the available experimental data. Based on the quasi-harmonic Debye model, the pressure and temperature dependence of the thermodynamic properties for hexagonal phase FeSe were investigated. Our theoretical calculations suggest that the pressure and temperature have significant effects on the heat capacity, vibrational internal energy, vibrational entropy, vibrational Helmholtz free energy, thermal expansion coefficient and Debye temperature. Even though few theoretical reports on the structural properties of FeSe are found in the current literature, to our knowledge, this is a novel theoretical investigation on the structural and thermodynamic properties of FeSe at high temperature. We hope that the theoretical results reported here can give more insight into the structural and thermodynamic properties of other iron-based superconductors at high temperature.     

Theoretical studies on the vibrational spectra, thermodynamic properties, detonation properties, and pyrolysis mechanisms for polynitroadamantanes

To look for high energy density materials (HEDM), the relationships between the structures and the performances of polynitroadamantanes (PNAs) were studied. The assigned infrared spectra of PNAs obtained at the density functional theory (DFT) B3LYP/6-31G level were used to compute the thermodynamic properties on the basis of the principle of statistical thermodynamics. The thermodynamic properties are linearly related with the number of nitro groups as well as with the temperatures. Detonation properties of PNAs were evaluated by using the Kamlet-Jacobs equation based on the calculated densities and heats of formation for titled compounds, and it is found that only when the number of nitro groups of PNA is equal to or more than eight can it be possible for PNAs to be used as HEDMs. The relative stabilities of PNAs were studied by the pyrolysis mechanism using the UHF-PM3 method. The homolysis of the C-NO2 bond is predicted to be the initial step of thermal decomposition. The activation energies (Ea) for the homolysis decrease with the number of nitro groups being increased on the whole. The stability order of dinitroadamantane isomers derived from the interactions among nitro groups is consistent with what is determined by Ea. The relations between the Ea's and the electronic structure parameters were discussed. In combination with the stability, PNA (1,2,3,4,5,6,7,8,9,10-) is recommended as the target of HEDM with insensitivity.