Inequivalence of substitution pairs in hydroxynaphthaldehyde: A theoretical measurement by intramolecular hydrogen bond strength,aromaticity, and excited‐state intramolecular proton transfer reaction

The inequivalence of substitution pair positions of naphthalene ring has been investigated by a theoretical measurement of hydrogen bond strength, aromaticity, and excited state intramolecular proton transfer (ESIPT) reaction as the tools in three substituted naphthalene compounds viz 1‐hydroxy‐2‐naphthaldehyde (HN12), 2‐hydroxy‐1‐naphthaldehyde (HN21), and 2‐hydroxy‐3‐naphthaldehyde (HN23). The difference in intramolecular hydrogen bond (IMHB) strength clearly reflects the inequivalence of substitution pairs where the calculated IMHB strength is found to be greater for HN12 and HN21 than HN23. The H‐bonding interactions have been explored by calculation of electron density ρ(r) and Laplacian ?²ρ(r) at the bond critical point using atoms in molecule method and by calculation of interaction between σ* of OH with lone pair of carbonyl oxygen atom using NBO analysis. The ground and excited state potential energy surfaces (PESs) for the proton transfer reaction at HF (6‐31G**) and DFT (B3LYP/6‐31G**) levels are similar for HN12, HN21 and different for HN23. The computed aromaticity of the two rings of naphthalene moiety at B3LYP/6‐31G** method also predicts similarity between HN12 and HN21, but different for HN23. © 2010 Wiley Periodicals, Inc. J Comput Chem, 2010     

Quantitative Correlation of Chromatographic Retention and Acute Toxicity for Alkyl（1- phenylsulfonyl） Cycloalkane Carboxylates and Their Structural Parameters by DFT

1 INTRODUCTION Quantitative structure-activity relationship (QSAR)is one of the necessary methods that could be emp- loyed to evaluate the hazards of organic chemicals. QSAR equation could be applied to predict the biological activity of unknown compounds, espe- cially for initial screening and evaluation of toxic compounds[1]. Moreover, the quantitative relation- ship between molecular structure and chromatogra- phic retention (capacity factor lgKW) could also bedeveloped to explain …     

偏最小二乘和量子化学参数预测酯类物质的闪点

基于化合物定量结构与性质的关系,对其中88个酯类化合物进行了B3LYP/6-31G*水平上的结构优化,并在优化结构基础上进行了量子化学结构参数的计算。应用偏最小二乘(PLS)方法对酯类化合物闪点温度这个理化性质与量子化学参数进行了关联。结果表明,分子热力学自由能(ΔG)结合偶极距(μ)可以很好地表达酯类化合物闪点温度与量子化学参数间的定量关系。所建立的QSPR模型具有较强的稳健性和预测能力。     

羰基硼化合物(BCO)n(n==1～12)的理论研究

运用HF/3-21G方法和密度泛函理论(DFT)的B3LYP/6-31G*方法, 对羰基硼化合物(BCO)n(n=1~12)的各种可能结构进行了优化, 对在B3LYP/6-31G*水平上得到的几何构型、电子态、结合能、振动频率、核独立化学位移(NICS)、能量二次差分和热力学性质进行了理论研究, 得到了(BCO)n(n=1~12)结构的稳定性信息. 十二种基态结构都是端配位(μ1-CO), (BCO)n(n=1~3, 5, 6) 的基态是线型或平面结构, (BCO)n(n=4, 7~12)的基态是笼状结构; B—C平均键能呈现奇偶交替现象, 偶数的结构比奇数稳定; 能量二次差分得到同样的结论；羰基的振动频率与实验值非常吻合; 热力学性质的研究对实验具有重要的指导意义.     

DFT Study of substituent effects of 2‐substituted alkyl ethyl methylcarbonates in homogeneous,unimolecular gas phase elimination kinetics

Theoretical studies on the gas phase elimination of 2‐substituted alkyl ethyl methylcarbonates were performed at the B3LYP/6‐31G* and B3LYP/6‐31+G** level of theory. The results of these calculations provide additional evidence that the mechanism of carbonates with a C_β H bond proceeds through a concerted nonsynchronous six‐membered cyclic transition state to produce methylcarbonic acid and the corresponding olefin. The unstable intermediate, methylcarbonic acid, rapidly decomposes through a four‐membered cyclic transition state to methanol and carbon dioxide. The correlation of the logarithm of theoretical rate coefficients against Hancock's steric parameters E gave an approximate straight line (δ = 0.30, r = 0.996 at 400°C). An additional fact is that when experimental log k_re.l is plotted against the theoretical log k_re.l. for 2‐alkyl ethyl methylcarbonates an approximate straight line (r = 0.997 at 400°C) is obtained. © 2006 Wiley Periodicals, Inc. Int J Chem Kinet 38: 184–193, 2006     

Theoretical Study of Heterocyclic Nitro Compounds for Use as Novel Explosives and Propellants

The heats of formation (HOFs) of heterocyclic nitro compounds were obtained by using a density functional theory B3LYP method with 6‐31G* and 6‐311+G** basis sets. The isodesmic reactions designed for the evaluation of HOFs keep most of the basic ring structures of the title compounds and thus ensure the credibility of the results. The values of HOFs are 567.90, 874.29 and 975.83 kJ/mol at the B3LYP/6‐31G* level for hexanitrohexazaadamantane ( A ), nonanitrononaza‐tetracyclo[7.3.1.1^3,7.1^5,11] pentadecane ( B ) and tetranitrotetrazacubane ( C ) respectively. The predicted detonation velocities of the title compounds are larger than, and detonation pressures are much larger than that of the widely used 1,3,5,7‐tetranitro‐1,3,5,7‐tetraazacyclooctane (HMX). The dissociation energy for the weakest C‐N bonds in the cage skeleton of the title compounds are 137‐144 kJ/mol at the B3LYP/6‐31G* level.     

Quantitative Correlation of the Acute Toxicity of Phenylthio-carboxylates with Their Structural and Thermodynamic Parameters by DFT Calculation

1 INTRODUCTION As industrial raw materials and pesticide, thio-ether, sulfoxide and sulfone have been used moreand more extensively. On the other hand, they havealso attracted chemists’ much attention due to theirpotential environmental hazards. Quantitative struc-ture-activity relationship (QSAR) is an effective me-thod evaluating the hazards of organic chemicals[1].QSAR equation could be employed to forecast thebiological activity of unknown compounds, which issignificant for initia…     

Theoretical Study on Gas—phase Pyrolytic Reactions of N—Ethyl,Isopropyl and N—t—Butyl Substituted 2—Aminopyrazine

Density functional theory (DFT) and ab initio methods were used to study gas‐phase pyrolytic reaction mechanisms of iV‐ethyl, N‐isopropyl and N‐t‐butyl substituted 2‐aminopyrazine at B3LYP/6–31G* and MP2/6–31G*, respectively. Single‐point energies of all optimized molecular geometries were calculated at B3LYP/6–311 + G(2d,p) level. Results show that the pyrolytic reactions were carried out through a unimolecular first‐order mechanism which were caused by the migration of atom H(17) via a six‐member ring transition state. The activation energies which were verified by vibrational analysis and correlated with zero‐point energies along the reaction channel at B3LYP/6–311 + G(2d,p) level were 252.02 kJ. mo^?1 (N‐ethyl substituted), 235.92 kJ‐mol^?1 (N‐t‐isopropyl substituted) and 234.27 kJ‐mol^?1 (N‐t‐butyl substituted), respectively. The results were in good agreement with available experimental data.     

Benchmarking of DFT functionals for the kinetics and mechanisms of atmospheric addition reactions of OH radicals with phenyl and substituted phenyl‐based organic pollutants

^•OH addition reactions play a pivotal role in the atmospheric transformation of a number of phenyl and substituted phenyl‐based persistent and toxic organic pollutants. Here, we screened appropriate DFT functionals to predict reaction mechanisms and rate constants (k_OH) of the ^•OH additions by taking benzene and substituted benzenes (C₆H₅F, C₆H₅Cl, C₆H₅Br, C₆H₅CH₃, C₆H₅OH) as model compounds. By comparing the k_OH values calculated with DFT methods to experimental values, we found that the ωB97 functional is the best among the 18 functionals considered (using the basis sets 6‐31 + G(d,p) for optimizations and 6‐311++G(3df,2pd) for single point energy calculations) in the temperature range of 230‐330 K. In addition, we found that some other functionals performed well in specific conditions, e.g., BMKD3 is good for benzene, halogenated benzenes and C₆H₅CH₃, and CAM‐B3LYP is good for the reaction of C₆H₅OH at room temperature. Based on the diversity of the electronic structures of the selected model compounds and the frequent occurrence of certain substituents ( CH₃,  OH,  F,  Cl, and  Br) in the target compounds, the functionals recommended here can be used for future study of the reaction mechanisms and k_OH values for ^•OH addition to phenyl and substituted phenyl‐based persistent and toxic organic pollutants.     

Structure,conformation and NMR studies on 1,2-dioxane and halogen substituted 1,2-dioxane molecules

Molecular structure and conformational stability of chair and twist conformers of 1,2-dioxane and halogen substituted compounds of the 1,2-dioxane have been studied using ab initio and density functional theory (DFT) methods. The molecular geometries of 1,2-dioxane, 3,6-difluoro, 3,6-dichloro, 3,3,6,6-tetrafluoro and 3,3,6,6-tetrachloro 1,2-dioxane compounds were optimized at HF, MP2, B3LYP and B3PW91 levels of theory by implementing 6-31G* basis set. To study the effect of polar medium, self-consistent reaction field theory is used to optimize the conformers at B3LYP/6-31G* level of theory. The geometrical parameters of chair and twist conformers have been discussed in the light of interaction between lone pair electrons present in the oxygen and substituted halogen atoms. The relative stability of the conformers have been studied using relative energy, maximum hardness principle and thermodynamical quantities. The 13C-NMR chemical shift study for carbon atoms in the title compounds are calculated and the results have been discussed.     

Theoretical Estimate of Hydride Affinities of Aromatic Carbonyl Compounds

Aromatic carbonyl compounds are one type of the most important organic compounds, and the reductions ofthem by hydride agents such as LiAlH4 or NaBH4 are widely used in organic synthesis. The reactivity of carbonyl compounds generally increases in the following order: ketone ＜ aldehyde, and amide ＜ acid ＜ ester ＜ acid halide, which could be related to their hydride affinities (HA). In the previous paper, Robert[1] calculated the absolute HAof a series of small non-aromatic carbonyl compounds. In this paper, we use DFT method at B3LYP/6-311 + + G (2d, 2p)∥B3LYP/6-31 + G* level to estimate hydride affinities of five groups of aromatic carbonyl compounds. The detailed results are listed in Table 1.     

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     

Controlling the Conformational Changes in Donor–Acceptor [4]‐Dendralenes through Intramolecular Charge‐Transfer Processes

The synthesis of two [4]‐dendralene compounds incorporating thiophene‐(p‐nitrophenyl) donor–acceptor units is presented. The dendralenes adopt two different conformers in solution and solid state and the transformation between the structures can be controlled by light and heat. The electron‐donating components of the dendralenes are represented by bromothienyl (in 13 ) and ethylenedioxythiophene(EDOT)‐thienyl (in 15 ) end‐groups. The most facile transformation involves the isomerisation of donor–acceptor conjugated systems ( a conformers) into structures in which only the thiophenes are conjugated ( b conformers), and this process is driven by ambient light. The structures of the two conformers of compound 13 are confirmed by single‐crystal X‐ray diffraction studies and the structural changes in both compounds have been monitored by ¹H NMR spectroscopy and absorption studies. The transformations were found to be first‐order processes with rate constants of k=0.0027 s^?1 and k=0.00022 s^?1 for 13 and 15 , respectively. Density functional theory calculations at the B3LYP/6‐31G* level give credence to the proposed mechanism for the a → b conversion, which involves photoinduced intramolecular charge transfer (ICT) as the key step. The EDOT derivative ( 15 ) can be polymerised by electrochemical oxidation and a combination of cyclic voltammetry and UV/Vis spectroelectrochemical experiments indicate that the a conformer can be trapped and stabilised in the solid state.