DFT Calculations and NBO Analysis of 2-Chloroethylethyldichlorosilane Unimolecular Elimination Kinetics in the Gas Phase

Ab initio molecular orbital calculations have been used to investigate the thermal decomposition kinetics of 2-chloroethylethyldichlorosilane at the B3LYP/6-311+G**,B3PW91/6-311+G**,and MPW1PW91/6-311+G** levels of theory.Among these methods,the results(activation parameters) obtained using the B3LYP/6-311+G** level are in good agreement with the available experimental data.The calculated data imply that in the unimolecular β-elimination reactions of the studied compound in the gas phase,the polarization of C(1)-Cl(3) and C(1)-H(4) bonds in the sense of C(1)δ+-Cl(3)δ-and C(1)δ+-H(4)δ-,respectively,is a determining factor in the gas phase elimination reactions 1,2 and 3.Analysis of bond order,natural bond orbital charges,bond indexes,synchro-nicity parameters,and IRC calculations suggest the elimination of 2-chloroethylethyldichlorosilane via reactions 1～3 can be described as concerted and slightly asynchronous.The transition state structures of these reactions are a four-membered cyclic structure.     

Theoretical Study and NBO Analysis of the Decomposition Kinetics of Oxetane, 2-Methyloxetane and 2,2-Dimethyloxetane

A theoretical study of the thermal decomposition kinetics of oxetane (1), 2-methyloxetane (2), and 2,2-dimethyloxetane (3) has been carried out at the B3LYP/6-311+G**, B3PW91/6-311+G**, and MPW1PW91/6-311+G** levels of theory. The MPW1PW91/6-311+G** method was found to give a reasonable good agreement with the experimental kinetics and thermodynamic parameters. The decomposition reaction of compounds 1～3 yields formaldehyde and the corresponding substituted olefin. Based on the optimized ground state geometries using MPW1PW91/6-311+G** method, the natural bond orbital (NBO) analysis of donor-acceptor (bond-antibond) interactions revealed that the stabilization energies associated with the electronic delocalization from σC3-C4 bonding to σ*O1-C2 antibonding orbitals decrease from compounds 1 to 3. The σC3-C4→σO1-C2 resonance energies for compounds 1～3 are 2.63, 2.59 and 2.45 kcal mol-1, respectively. Further, the results showed that the energy gaps between σC3-C4 bonding and σ*O1-C2 antibonding orbitals decrease from compounds 1 to 3. Also, the decomposition process in these compounds are controlled by σ→σ* resonance energies. Moreover, the obtained order of energy barriers could be explained by the number of electron-releasing methyl groups substituted to the Csp3 atom (which is attached to oxygen atom). NBO analysis shows that the occupancies of σCsp3-O bonds decrease for compounds 1～3 as 3<2<1, and those of σCsp3-O bonds increase in the opposite order (3 > 2 > 1). This fact illustrates a comparatively easier thermal decomposition of the sCsp3-O bond in compound 3 compared to compound 2, and in compound 2 compared to compound 1. NBO results indicate that these reactions are occurring through a concerted and asynchronous four-membered cyclic transition state type of mechanism.     

Density Functional Calculations of C–NO_2 Bond Dissociation Energies for Nitroalkanes Molecules

Bond dissociation energies for the removal of nitrogen dioxide group in some nit- roalkane energetic materials have been calculated by using the three hybrid density functional theory (DFT) methods B3LYP, B3PW91 and B3P86 with 6-31g** and 6-311g** basis sets. The computed BDEs have been compared with the available experimental results. It is found that the B3P86 method with 6-31g** and 6-311g** basis sets can obtain satisfactory bond dissociation energies (BDEs), which are in extraordinary agreement with the experimental data. Considering the smaller mean absolute deviation and maximum difference, the reliable B3P86/6-311g** method was recommended to compute the BDEs for the removal of nitrogen dioxide group in the nitroalkane energetic materials. Using the method, the BDEs of 8 other nitroalkane energetic materials have been calculated and the maximum difference from experimental value is 1.76 kcal·mol-1 (for the BDE of tC4H9–NO2), which further proves the reliability of B3P86/6-311g** method. In addition, it is noted that the BDEs of C–NO2 bond change slightly for main chain nitroalkane compounds with the maximum difference of only 3.43 kcal mol-1.     

Theoretical Studies on Intermolecular Hydrogen-bond Interactions between Hexamethylenetetramine and Nitric Acid

The structures of the complexes generated by hexamethylenetetramine and nitric acid have been fully optimized by B3LYP method at the 6-311++G** and aug-cc-pVTZ levels. The intermolecular hydrogen-bonding interactions have been calculated by the B3LYP/6-311++G**, B3LYP/aug-cc-pVTZ, MP2(full)/6-311++G** and CCSD(T)/6-311++G** methods, respectively. The NBO (nature bond orbital), AIM (atom in molecule), temperature effect and solvation effect have been analyzed to reveal the origin of the interactions. The results indicate that the stable hydrogen-bonded complexes could be generated by hexamethylenetetramine and nitric acid. The interactions follow the order of (a)>>(e)>(b)>(c)>(d)>(f)>(g). The C-N bonds which are adjacent to the methylene involving the hydrogen bonds tend to break in the chemical reaction. Due to the exothermic process, low temperature is conducive to the formation of the composition, which tallies with the experimental result.     

Structures and Aromaticity of Three-membered BeXP （X = C,Si, Ge） and CYP （Y = O,S, Se） Rings

Three-membered BeXP （X = C, Si, Ge） and CYP （Y = O, S, Se） rings are theoretically investigated using density functional theory （DFT） methods at the B3LYP/6-311＋G^＊ and B3PW91/6-311＋G^＊ levels of theory. The research results show that the size of atoms has a great influence on the structural stability of these species. The wiberg bond indexs （WBIs） suggest the existence of delocalization in these structures. Negative nucleus-independent chemical shift （NICS） values for these species indicate that a strong ring current exists in these three-membered structures （Cs symmetry）. A detailed molecular orbital （MO） analysis further reveals that a delocalized π MO strengthens the structural stability and makes these species show strong aromatic characters.     

Ab Initio and RRKM Study of the Unimolecular Decomposition of Five-membered Ring Heterocyclic （X：O,S,N-H）

Quantum mechanical and Rice-Ramsperger-Kassel-Marcus (RRKM) calculations are carried out to study the thermal unimolecular decomposition of 2,5-dihydrofuran (1),2,5dihydrothiophene (2),and 3-pyrroline (3) at the MPW1PW91/6-31++G level of theory,and the results are in good agreement with the experimental values.The predicted high pressure limit rate constants (k(T)) in various states of activation energy and pre-exponential (S1:(A(calc.),E a(calc.)),S2:(A (calc.),E a(exp.)),and S3:(A (exp.),E a(calc.))) for the thermal decomposition processes 1-3 were evaluated.Also,the fall-off pressures (P1/2) for compounds 1-3 in the states 1-3 are found to be (1.24×10-2,1.09×10-3,and 4.19×10-2mmHg),(1.24×10-2,1.63×10-3,and 2.79×10-2mmHg),and (1.24×10-2,1.63×10-3,and 4.19×10-2 mmHg),respectively.As the fall-off pressure of thermal decomposition process of compounds 1-3 is in the following order:P1/2(3)> P1/2(1)> P1/2 (2),the decomposition rates are as below:rate(3) 相似文献   

Density Functional Calculations of C-NO2 Bond Dissociation Energies for Nitroalkanes Molecules

Bond dissociation energies for the removal of nitrogen dioxide group in some nitroalkane energetic materials have been calculated by using the three hybrid density functional theory （DFT） methods B3LYP, B3PW91 and B3P86 with 6-31g^＊＊ and 6-311g^＊＊ basis sets. The computed BDEs have been compared with the available experimental results. It is found that the B3P86 method with 6-31g^＊＊ and 6-311g^＊＊ basis sets can obtain satisfactory bond dissociation energies （BDEs）, which are in extraordinary agreement with the experimental data. Considering the smaller mean absolute deviation and maximum difference, the reliable B3P86/6-311g^＊＊ method was recommended to compute the BDEs for the removal of nitrogen dioxide group in the nitroalkane energetic materials. Using the method, the BDEs of 8 other nitroalkane energetic materials have been calculated and the maximum difference from experimental value is 1.76 kcal.mo1^-1 （for the BDE of tC4Hg-NOz）, which further proves the reliability of B3P86/6-311g^＊＊ method. In addition, it is noted that the BDEs of C-NO2 bond change slightly for main chain nitroalkane compounds with the maximum difference of only 3.43 kcal mo1^-1.     

Theoretical Study of the N-NO2 Bond Dissociation Energies for Energetic Materials with Density Functional Theory

The N-NO2 bond dissociation energies （BDEs） for 7 energetic materials were computed by means of accurate density functional theory （B3LYP, B3PW91 and B3P86） with 6-31G＊＊ and 6-311G＊＊ basis sets. By comparing the computed energies and experimental results, we find that the B3P86/6-311G＊＊ method can give good results of BDE, which has the mean absolute deviation of 1.30kcal/mol. In addition, substituent effects were also taken into account. It is noted that the Hammett constants of substituent groups are related to the BDEs of the N-NO2 bond and the bond dissociation energies of the energetic materials studied decrease when increasing the number of NO2 group.     

Theoretical Studies on the Intermolecular Interactions of Aza-calix[2]arene[2]-triazines with RDX

Six fully optimized structures of the aza-calix[2]arene[2]-triazines/RDX supramo-lecular complexes have been obtained at the DFT-B3LYP/6-311++G** level,and the corresponding intermolecular interactions have been investigated using the B3LYP,mPWPW91 and MP2 methods at the 6-311++G** level,respectively.The natural bond orbital(NBO) and atoms in molecules(AIM) analyses have been performed to reveal the origin of interactions.To our interest,the result indicates that the strongest interaction is up to -22.34 kJ/mol after basis set superposition error(BSSE) and zero point energy(ZPE) correction at the MP2/6-311++G** level.Furthermore,the intermolecular interactions between aza-calix[2]arene[2]-triazines with the substituted amidos and RDX are stronger than those of other complexes.Thus,the complexes with amidos can be used as the candidates to increase the stability of explosive and eliminate the explosive wastewater.     

Toxicity （-lgEC50） Measurement of the Fluorobenzene Derivants against Vibrio Qinghaiensis （Q67） and Their 2D,3D-QSAR Study

Toxicities (-lgEC50) of 16 fluorobenzene derivants against vibrio qinghaiensis (Q67) were measured systematically,and their quantum chemistry parameters were calculated at the B3LYP/6-311G** level. Based on the experimental toxicity data and quantum chemistry parameters,2D-QSAR model was proposed,which was validated by variance inflation factors (VIF),t-value and cross-validation method. At the mean time,comparative molecular force field (CoMFA) based on molecular simulation was used to investigate the toxicity of fluorobenzene derivants. Furthermore,the intoxicating mechanism of fluorobenzene derivants was discussed. To our interest,2D-QSAR and CoMFA models exhibit good prediction ability,with which the toxicity of similar compounds can be predicted. Finally,toxicities (-lgEC50) of 12 fluorobenzene derivants against vibrio qinghaiensis (Q67) were predicted with these models.     

NMR Parameters and Vibrational Investigation of 2-Dicyanovinyl-5-（4-ethoxyphenyl） thiophene by ab initio HF and DFT Calculations

Optimized geometry and harmonic vibrational frequency of 2-dicyanovinyl-5-(4- ethoxyphenyl)thiophene (C16H12N2OS) are calculated at the HF/6-31++G(d,p) and B3LYP/6- 311++G(d,p) levels. Mulliken charges in the ground state are also calculated. The research shows the presence of intermolecular interaction in the title compound. The scaled harmonic vibrational frequencies have been compared with experimental FT-IR spectra. A detailed interpretation of the infrared spectra of the title compound is reported. The theoretical spectrograms for IR spectra of the title compound have been constructed. The isotropic chemical shift computed by 13C and 1H NMR analyses also shows good agreement with the experimental observations.     

Theoretical Investigation on the Geometries and Properties of Guanine-BX3 （X = F, Cl） Complex

Geometries and binding energies were predicted at the B3LYP/6-311+G* level for the guanine-BX3 (X = F, Cl) systems and four isomers with no imaginary frequencies have been obtained for both guanine-BF3 and guanine-BCl3, respectively. Single energy calculations using much larger basis sets (6-311+G(2df,p) and aug-cc-pVDZ were carried out as well. It was found that the most stable isomer of guanine-BF3 is BF3 connected to N3 of guanine with the stabilization energy of –19.93 kcal/mol (BSSE corrected), while that of guanine-BCl3 is BCl3 connected to O10 of guanine having stabilization energy of –15.02 kcal/mol at the same level. The analyses for the combining interaction between BX3 and guanine with the atom-in-molecules theory (AIM) and natural bond orbital (NBO) methods have been performed. The results indicated that all the isomers are formed with σ-p type interactions between guanine and BX3, in which pyridine-type nitrogen or carbonyl oxygen or nitrogen atom of amino group offers its lone pair electrons to the empty p orbital of boron atom and the concomitance of charge transfer from guanine to BX3 has occurred. Still, one or two hydrogen bonds exist in some isomers of guanine-BX3 system and contribute to the stability of complex systems. Frequency analysis suggested that the stretching vibration of BX3 undergoes a red shift in complexes. Guanine-BF3 complex is more stable than guanine-BCl3 although the B–Y (Y=N, O) bond distance in the latter is shorter.     

Crystal Structure and DFT Studies of a Triazole Derivative： 4- （ 2-Hydrobenzylideneamino ） -3- （ 1,2,4-triazol-4-ylmethyl） 1H-1, 2,4-triazole-5 （4H） -thione

A novel triazole derivative 4-(2-hydrobenzylideneamino)-3-(1,2,4-triazol-4-ylmethyl)-1H-1,2,4-triazole-5(4H)-thione(1) was synthesized and characterized using elemental analysis, FTIR, and 1H NMR, and its crystal structure was determined via X-ray single crystal diffraction analysis. Crystal data: monoclinic, P2(1)/c, a=0.83335(9) nm, b=1.49777(16) nm, c=1.14724(12) nm, β =107.990(2)°, D=1.470 Mg/m3, and Z=4. The geometries and the vibrational frequencies were determined using the density functional theory(DFT) method at the B3LYP/ 6-31G level. To demonstrate the accuracy of the reaction route of compound 1, one of the important intermediates was also tested using the same method. The structural parameters of the two compounds calculated using the DFT study are close to those of the crystals, and the harmonic vibrations of the two compounds computed via the DFT method are in good agreement with those in the observed IR spectral data. The thermodynamic properties of the title compound were calculated, and the compound shows a good structural stability at normal temperature. The test results of biological activities show that it has a certain bactericidal ability.     

Theoretical Investigations into the Intermolecular Hydrogen-bonding Interactions between Azacyclopentane-2-one and N-Methylol Ethanone

The structures of the complexes formed between N-methylol ethanone(model molecule of ceramide) and azacyclopentane-2-one(the model molecule of azone) have been fully optimized at the B3LYP/6-311++G** level.The intermolecular hydrogen bonding interaction energies have been calculated by using the B3LYP/6-311++G**,B3LYP/6-311++G(2df,2p),MP2(full)/6-311++G** and MP2(full)/6-311++G(2df,2p) methods,respectively.The results show that strong O–H···O=C,N–H···O=C and C–H···O=C hydrogen bonds could exist between azacyclopentane-2-one and N-methylol ethanone.The formation of the complexes might change the conformation of ceramide molecule and thus cause better percutaneous permeation for the drugs.This is perhaps the origin of the permeation enhances the activity of azone for medicament,as is in accordance with the experimental results.The hydrogen-bonding interactions follow the order of(a) (c) (b) (d) (g) ≈(e) ≈(i) (h) (f).The analyses of frequency,NBO,AIM and electron density shift are used to further reveal the nature of the complex formation.In the range of 263.0~328.0 K,the complex is formed via an exothermic reaction,and the solvent with lower temperature and dielectric constant is favorable to this process.     

3-单酰胺喹唑啉酮衍生物的表征及密度泛函理论研究

Several 3-aminoquinazolin-4-(3H)-one derivatives were synthesized and characterized.Using proton nuclear magnetic resonance (NMR) spectra,we have investigated the barriers to rotation around the N ― N bond as a function of temperature.Changes in the NMR spectra at high temperatures are explained in terms of hindered rotations of the N ― N bond.Free energies of activation for the rate determining stereochemical process were calculated to be as high as 67-75 kJ · mol-1.Ground state molecular geometries and vibrational frequencies were calculated using the HF/6-31G and B3LYP/6-31G level of theories.The optimized bond lengths and bond angles are in good agreement with experimental values at both theoretical levels.     

Theoretical Study on Hydrogen Interaction between 5-Fluorouracil and Glycine

The geometries,electronic structure,IR spectrum and other properties of hydrogen interaction between 5-fluorouracil and glycine were studied at the B3LYP/6-31+G* level.Single point energy calculations were executed at the B3LYP/6-311++G** and B3LYP/aug-cc-pvdz levels,and natural bond orbital (NBO) analysis was carried out at the B3LYP/6-31+G* level.Finally,the hydrogen bonds were discussed via AIM electronic density topology analysis.     

FT-IR, Raman and NMR Spectra, Molecular Geometry, Vibrational Assignments, ab initio and Density Functional Theory Calculations for Diethyl Phthalate

Phthalate acid esters (PAEs) possess endocrine disruptive effects and can produce reproductive and developmental toxicities. In this paper, both experimental and theoretical studies on FT-IR, Raman and 1H NMR spectra of diethyl phthalate (DEP) have been carried out. The geometrical structure of DEP was optimized at the HF/6-31G*, HF/6-311G**, B3LYP/6-31G*, and B3LYP/6-311G** levels, respectively. The harmonic vibrational frequencies, IR intensity, Raman activity and 1H NMR chemical shifts have been computed at the B3LYP/6-31G* and B3LYP/6-311G** levels. Anharmonic corrections to frequencies were obtained by means of second-order perturbation theory (PT2) at the B3LYP/6-31G* level. Based on potential energy distribution (PED), the vibrational assignments have also been performed. The theoretical calculation values were compared with the experimental observations and the results indicate they are in excellent agreement.     

（E）-1-[4-（2-hydroxy-5-methoxybenzylideneamino）-phenyl]ethanone and （E）-1-[4-（2-hydroxy-4-methoxybenzylideneamino） phenyl]ethanone：X-ray and DFT-calculated structures

The isomeric structures of(E)-1-[4-(2-hydroxy-5-methoxybenzylideneamino)-phenyl] ethanone(I) and(E)-1-[4-(2-hydroxy-4-methoxybenzylideneamino) phenyl]ethanone(Ⅱ) ,both C16H15NO3,have been determined using X-ray diffraction techniques and characterized by IR,and their molecular structures have also been optimized at the B3LYP/6-31G(d,p) level using density functional theory(DFT) . The energetic behaviors of the title compounds in solvent media have been examined using B3LYP method with the 6-31G(d,p) basis set by applying the polarizable continuum model(PCM) . The total energies of the title compounds decrease with the increasing polarity of the solvent. In addition,DFT calculations of the title compounds' molecular electrostatic potentials(MEP) were performed at the B3LYP/6-31G(d,p) level of theory. X-ray study shows that the title compounds both have strong intramolecular O-H…N hydrogen bonds. The molecules of Ⅰ are linked into a one-dimensional framework structure by C-H…π interactions,while in Ⅱ,intermolecular π···π interactions result in the formation of infinite chains running along the [010].