Quantum chemical studies on structure, conformation and isomerization of nitroso, nitro substituted benzene and 1,3-cyclopentadiene

The molecular structure, conformational stability and isomerization of nitroso, nitro substituted benzene and 1,3-cyclopentadiene in gas phase have been investigated using ab initio and density functional theory methods. The molecular geometries and energetics of possible conformers were obtained by employing MP2, B3LYP and B3PW91 levels of theory implementing 6-31G* basis set. The relative stabilities of the conformations were evaluated from the energy differences of the structure. Chemical hardness (η) and chemical potential (μ) were calculated at HF/6-31G* level of theory for all the positional and geometrical isomers to study the maximum hardness principle. Each optimized structure has been tested against the imaginary frequencies at MP2/6-31G* level of theory in order to be sure they are located at energy minimum.     

异丙醇的O-H伸缩振动泛频光谱和分子构像

使用高灵敏的光腔衰荡光谱(CavityRingDownSpectroscopy)技术测出了异丙醇的O-H伸缩v=4、5振动泛频光谱,每个振动能级都有三个吸收峰,被归属为分子构像的O-H伸缩泛频吸收.给出了光腔衰荡光谱的振动泛频吸收的谱带强度公式,并求得分子不同构像在不同振动能级的O-H伸缩泛频吸收的谱带强度;同时利用局域模理论,求得分子各O-H伸缩局域模振子的机械频率(X1)、非谐性(X2)以及解离能(D).用密度泛函(DFT)B3LYP/6-31+G*理论方法优化了分子的各种可能构像,验证了分子存在反式(trans)和偏转(gauche)两种稳定构像,计算的分子的O-H伸缩频率及构像稳定性同实验结果是一致的.     

Theoretical study of the structure and vibrational spectra of the Ci(2) conformation of 18-crown-6

Optimized geometries of 18-crown-6 (18ce6) were calculated at the HF/6-31G* and B3LYP/6-31G* levels of theory for the D3d, Ci(1) and Ci(2) conformations. At the B3LYP level, the Ci(2) optimized geometry was higher in energy by 23.3 and 18.8 kcal mol(-1) than the Ci(1) and D3d optimized geometries, respectively. Harmonic force field, vibrational frequencies and IR absorption intensities were calculated at the corresponding optimized geometry at the B3LYP level for the Ci(2) conformation. Scaled Ci(2) frequencies were compared with the experimental frequencies of free 18ce6, Ci(1) conformation, and 18ce6-urea complex, Ci(2) conformation. This comparison showed possible misassignments in the fundamental vibrational frequencies of 18ce6.     

Proton affinity of proline and modified prolines using the kinetic method: role of the conformation investigated by ab initio calculations

The proton affinities of proline, cis-3-methylproline and cis-3-ethylproline have been measured by the kinetic method using an ion trap instrument; the values obtained are 936, 940.5, and 943 kJ mol(-1), respectively. The experimental values are consistent with those obtained by high-level ab initio calculations (B3LYP/6-31+G*//B3LYP/6-31G* and B3P86/6-31+G*//B3LYP/6-31G*). Several conformations of neutral and protonated proline were considered, in particular the endo and exo ring structure and the position of the carboxyl group. These results show the importance of the position of the hydrogen atom of the carboxyl group in determining the most stable protonated proline structure.     

FC(O)NCS 分子振动光谱的理论研究

采用DFT(B3LYP)方法,以6-3G*为基组对FC(O)NCS的顺式和反式两种构型的几何结构,振动谐性力场和红外光谱进行了研究。B3LYP/6-31G*计算水平和大多数有机分子的一套固定标度因子进行标度。根据标度后的理论力场进行简正坐标分析得到的势能分布(PED)和红外光谱强度值对FC(O)NCS分子的顺式和反式两种构型的振动基频进行了理论归属。     

Non-aqueous solvation of n-octanol and ethanol: spectroscopic and computational studies

Raman spectroscopy was used to examine the interactions of the free O-H bonds in n-octanol and ethanol with the organic solvents carbon tetrachloride (CCl(4)), cyclohexane, and benzene. These spectra reveal that the solvents CCl(4) and cyclohexane have a small effect on the free O-H peak of alcohols, whereas benzene as a solvent significantly red-shifts the free O-H band. Calculated spectra were generated via MP2/6-31G* calculations and the B3LYP/6-31+G**//MP2/6-31G*-derived Boltzmann populations of each ethanol complex and are consistent with the experimental results. Additional spectra were calculated using Boltzmann populations derived from single-point energies at the polarizable continuum model (PCM) level with the B3LYP/6-31+G** level of theory to take overall solvent effects into account, and these simulated spectra are also largely consistent with the experimental results. Analysis of the computational results reveals a lengthening of the O-H bond from the O-H interaction with the delocalized electronic structure of benzene as well as a bimodal distribution of the free O-H peak of the alcohol/benzene mixtures due to two distinctly different types of alcohol/benzene complexes.     

Conformational study of the structure of 12-crown-4-alkali metal cation complexes

A conformational search was performed for the 12-crown-4 (12c4)-alkali metal cation complexes using two different methods, one of them is the CONFLEX method, whereby eight conformations were predicted. Computations were performed for the eight predicted conformations at the HF/6-31+G*, MP2/6-31+G*//HF/6-31+G*, B3LYP/6-31+G*, MP2/6-31+G*//B3LYP/6-31+G*, and MP2/6-31+G* levels. The calculated energies predict a C4 conformation for the 12c4-Na+, -K+, -Rb+, and -Cs+ complexes and a C(s) conformation for the 12c4-Li+ complex to be the lowest energy conformations. For most of the conformations considered, the relative energies, with respect to the C4 conformation, at the MP2/6-31+G*//B3LYP/6-31+G* are overestimated, compared to those at the MP2/6-31+G* level, the highest level of theory considerd in this report, by 0.2 kcal/mol. Larger relative energy differences are attributed to larger differences between the B3LYP and MP2 optimized geomtries. Binding enthalpies (BEs) were calculated at the above-mentioned levels for the eight conformations. The agreement between the calculated and experimental BEs is discussed.     

Spectroscopic analysis of 2-(2,3-dihydro-1,5-dimethyl-3-oxo-2-phenyl-1H-pyrazol-4-ylimino)-2-(4-nitro-phenyl) acetonitrile

The optimized geometry and vibrational frequencies of 2-(2,3-dihydro-1,5-dimethyl-3-oxo-2-phenyl-1H-pyrazol-4-ylimino)-2-(4-nitro-phenyl) acetonitrile (DOPNA) were obtained by ab initio DFT/B3LYP level with complete relaxation in the potential energy surface using 6-31G and 6-311G basis sets. The Fourier-transform infrared (FT-IR) spectrum of DOPNA has been recorded in the region 4000-400 cm(-1). The harmonic vibrational frequencies were calculated and the scaled values have been compared with experimental FT-IR spectrum. The calculated frequencies are in comparable agreement with the experimental frequencies. The calculated energy span between the HOMO and the LUMO of DOPNA is 2.94 and 2.87eV by B3LYP/6-31G and B3LYP/6-311G, respectively.     

2,3,7,8-四氯苯并二英分子结构与热力学性质的理论研究

用密度泛函理论(DFT)和从头算(ab initio)方法,在B3LYP/6-31G、 B3LYP/6-31G*、 B3LYP/6-311G*和MP2/6-31G*水平上全优化计算了2,3,7,8-四氯苯并二英(2,3,7,8-TCDD)的几何构型、电子结构和振动频率,并用校正后的频率计算了298～1500 K的标准热力学函数,同时用半经验的PM3 SCF-MO进行了同样的计算,计算结果与实验值及文献值较好地吻合.     

Molecular structure and vibrational spectra of three substituted 4-thioflavones by density functional theory and ab initio Hartree-Fock calculations

The vibrational frequencies of three substituted 4-thioflavones in the ground state have been calculated using the Hartree-Fock and density functional method (B3LYP) with 6-31G* and 6-31+G** basis sets. The structural analysis shows that there exists H-bonding in the selected compounds and the hydrogen bond lengths increase with the augment of the conjugate parameters of the substituent group on the benzene ring. A complete vibrational assignment aided by the theoretical harmonic wavenumber analysis was proposed. The theoretical spectrograms for FT-IR spectra of the title compounds have been constructed. In addition, it is noted that the selected compounds show significant activity against Shigella flexniri. Several electronic properties and thermodynamic parameters were also calculated.     

Experimental and theoretical study of the vibrational spectra of 12-crown-4-alkali metal cation complexes

The vibrational, Raman, and IR, spectra of the five 12-crown-4 (12c4) complexes with Li+, Na+, K+, Rb+, and Cs+ alkali metal cations were measured. Except for a small shift of the position of some bands in the vibrational spectra of the Li+ complex, the vibrational spectra of the five complexes are so similar that it is concluded that the five complexes exist in the same conformation. B3LYP/6-31+G* force fields were calculated for six of the eight predicted conformations in a previous report (J. Phys. Chem. A 2005, 109, 8041) of the 12c4-Li+, Na+, and K+ complexes that are of symmetries higher than the C1 symmetry. These six conformations, in energy order, are of C4, Cs, Cs, C(2v), C(2v), and Cs symmetries. Comparison between the experimental and calculated vibrational frequencies assuming any of the above-mentioned six conformations shows that the five complexes exist in the C4 conformation. This agrees with the fact that the five alkali metal cations are larger than the 12c4 ring cavity. The B3LYP/6-31+G* force fields of the C4 conformation of the Li+, Na+ and K+ complexes were scaled using a set of eight scale factors and the scale factors were varied so as to minimize the difference between the calculated and experimental vibrational frequencies. The root-mean-square (rms) deviations of the calculated frequencies from the experimental frequencies were 7.7, 5.6, and 5.1 cm(-1) for the Li+, Na+, and K+ complexes, respectively. To account for the earlier results of the Li+ complex that the Cs conformation is more stable than the C4 conformation by 0.16 kcal/mol at the MP2/6-31+G* level, optimized geometries of the complex were calculated for the C4 and Cs conformations at the MP2/6-311++G** level. The C4 conformation was calculated to be more stable than the Cs conformation by 0.13 kcal/mol.     

硝酸甲与硝酸乙酯的结构、振动频率和热力学性质的密度函数 理论(DFT)研究

用密度函数理论(DFT)的BLYP和B3LYP方法,取6-31G,6-31G^*,6-31G^*^*,6-311G,6-311G^*和6-311G^*^*六种基组,对硝酸甲酯和硝酸乙酯的几何构型和红外振动频率进行了计算研究.结果表明,B3LYP方法在采用极化基组(6-31G^*,6-31G^*^*,6-311G^*和6-311G^*^*)时计算得到的结果均较好,适用于硝酸酯类化合物的研究.而BLYP方法无论采用何种基组均不适用;运用校正后的B3LYP/6-31G^*频率(校正因子0.975)计算得到的热力学性质(C^o~p,H^o和S^o)与实验结果较吻合。     

Computational studies on carbohydrates: I. Density functional ab initio geometry optimization on maltose conformations

Ab initio geometry optimization was carried out on 10 selected conformations of maltose and two 2‐methoxytetrahydropyran conformations using the density functional denoted B3LYP combined with two basis sets. The 6‐31G* and 6‐311++G** basis sets make up the B3LYP/6‐31G* and B3LYP/6‐311++G** procedures. Internal coordinates were fully relaxed, and structures were gradient optimized at both levels of theory. Ten conformations were studied at the B3LYP/6‐31G* level, and five of these were continued with full gradient optimization at the B3LYP/6‐311++G** level of theory. The details of the ab initio optimized geometries are presented here, with particular attention given to the positions of the atoms around the anomeric center and the effect of the particular anomer and hydrogen bonding pattern on the maltose ring structures and relative conformational energies. The size and complexity of the hydrogen‐bonding network prevented a rigorous search of conformational space by ab initio calculations. However, using empirical force fields, low‐energy conformers of maltose were found that were subsequently gradient optimized at the two ab initio levels of theory. Three classes of conformations were studied, as defined by the clockwise or counterclockwise direction of the hydroxyl groups, or a flipped conformer in which the ψ‐dihedral is rotated by ∼180°. Different combinations of ω side‐chain rotations gave energy differences of more than 6 kcal/mol above the lowest energy structure found. The lowest energy structures bear remarkably close resemblance to the neutron and X‐ray diffraction crystal structures. © 2000 John Wiley & Sons, Inc. * J Comput Chem 21: 1204–1219, 2000     

1. 5-环辛二烯-3, 7-二炔(C~8H~4)结构和光谱性质的理论研究

采用量子化学从头算方法研究了1,5-环辛二烯-3,7-二炔(C~8H~4)的结构和光谱性质,根据等键反应分析和自然键轨道方法研究了它的稳定性、成键情况和共轭性。结果表明1,5-环辛二烯-3,7-二炔(C~8H~4)分子为平面刚性结构,可能稳定存在。分子中C≡C键与C=C键存在一定程度的共轭,可能具有芳香性。     

四唑互变异构反应的密度泛函理论(DFT)研究

运用11种密度泛函理论方法对四唑互变异构反应进行了计算研究。结果表明,B3LYP-DFT法与从头算的优化几何和能量最为吻合;在6-31^*基组下B3LYP计算的IR频率与MP2/6-311G^*^*计算结果相差很小;用未经校正的B3LYP计算频率求得的产物(2H-四唑)的热力学性质与实测结果也完全一致;由此推荐B3LYP-DFT法适合于对四唑化合物作系统研究。     

2-丁基-四氢噻吩亚砜13C-NMR的理论研究

在RHF/6-31G和B3LYP/6-31G水平上对顺式(Cis-)与反式(Trans-)2-丁基-四氢噻吩亚砜(BTHTO)进行几何优化,应用规范不变原子轨道法(GIAO)在6-31G、6-31+G、6-31++G和6-31+G(2d,p)水平上计算了Cis-和Trans-BTHTO的13C-NMR,对13C-NMR谱进行了归属。结果表明,BTHTO噻吩五元环的稳定构象呈半椅式,Cis-和Trans-BTHTO中与硫原子直接碳原子13C-NMR的显著差异主要是由于空间构型不同引起分子的静电势场对相应碳原子的屏蔽作用不同所致。     

二氯亚甲基锗烯与乙烯环加成反应机理的理论研究

The mechanism of a cycloaddition reaction between singlet dichloromethylene germylene and ethylene has been investigated with B3LYP/6-31G＊ method, including geometry optimization and vibrational analysis for the involved stationary points on the potential energy surface. Energies for the involved conformations were calculated by CCSD（T）//B3LYP/6-31G＊ method. On the basis of the surface energy profile obtained with CCSD（T）// B3LYP/6-31G＊ method for the cycloaddition reaction between singlet dichloromethylene germylene and ethylene, it can be predicted that the dominant reaction pathway is that an intermediate INT1 is firstly formed between the two reactants through a barrier-free exothermic reaction of 61.7 kJ/mol, and the intermediate INT1 then isomerizes to an active four-membered ring product P2.1 via a transition state TS2, an intermediate INT2 and a transition state TS2.1, in which energy barriers are 57.7 and 42.2 kJ/mol, respectively.     

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键扭转振动基频的预测值。     

Semiempirical PM5 molecular orbital study on chlorophylls and bacteriochlorophylls: comparison of semiempirical, ab initio, and density functional results

The semiempirical PM5 method has been used to calculate fully optimized structures of magnesium-bacteriochlorin, magnesium-chlorin, magnesium-porphin, mesochlorophyll a, chlorophylls a, b, c(1), c(2), c(3), and d, and bacteriochlorophylls a, b, c, d, e, f, g, and h with all homologous structures. Hartree-Fock/6-31G* ab initio and density functional B3LYP/6-31G* methods were used to optimize structures of methyl chlorophyllide a, chlorophyll c(1), and methyl bacteriochlorophyllides a and c for comparison. Spectroscopic transition energies of the chromophores and their 1:1 or 1:2 solvent complexes were calculated with the Zindo/S CIS method. The self-consistent reaction field model was used to estimate solvent shifts. The PM5 calculations predict planar structure of the porphyrin ring and central position of the four coordinated magnesium atoms in all pigments studied, in accord with the experimental, ab initio, and density functional results, a significant improvement as compared to the older semiempirical PM3 approach. Only small differences in PM5 and B3LYP/6-31G* or Hartree-Fock/6-31G* minimum energy geometries of the reference molecules were observed. Calculations show that in 1:1 solvent complexes, where the magnesium atom is five coordinated, the magnesium atom is shifted out of the plane of the porphyrin ring towards the solvent molecule, while the hexa coordinated 1:2 complexes are again planar. The PM5 method gives atomic charges that are comparable with those obtained from the Hartree-Fock/6-31G* and B3LYP/6-31G* calculations. The single point ZINDO/S CIS calculations with PM5 minimum energy structure gave excellent correlations between calculated and experimental transition energies of the chlorophylls and bacteriochlorophylls studied. Such correlations may be used for prediction of transition energies of the chromophores in protein binding sites. Calculations also predict existence of dark electronic states below the main Soret absorption band in all chromophores studied. The results suggest that the semiempirical PM5 method is a fairly reliable and computationally efficient method in predicting molecular parameters of porphyrin-like molecules.     

线性簇合物SC2nS2－(n =1～12)电子吸收光谱

应用密度泛函理论,在B3LYP/6-31G*水平上优化了线性簇合物SC2nS2－(n =1～12)的基态平衡几何结构,并计算了它们的谐振动频率.在基态平衡构型下,通过TD-B3LYP/cc-pvTZ和TD-B3LYP/cc-pvDZ计算,确定了簇合物SC2nS2－(n =1～10) 电子跃迁的垂直激发能和对应的振子强度.基于计算结果,导出了电子跃迁吸收波长与体系大小n的解析关系式,以及SC2nS2－体系第一电离能与体系大小n的解析表达式,并讨论了不同端位原子对碳链体系激发态性质的影响.