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.     

Cyclic vinyl p-tolyl sulfilimines as chiral dienophiles: theoretical study on the stereoselectivity, Lewis acid catalysis, and solvent effects in their Diels-Alder reactions

The theoretical study reported in the present work deals with chiral cyclic vinyl sulfilimines and their reactivity as dienophiles in [4 + 2] cycloaddition reactions, using B3LYP/6-31G(d)//AM1 and B3LYP/6-31G(d)//B3LYP/6-31G(d) model chemistries. Consideration of Lewis acid catalysis, illustrated by BF(3), decreases the activation energies of the cycloaddition process while the charge transfer from the diene to the sulfilimine is augmented. The [4 + 2] cycloaddition reactions of sulfilimines with both furan and cyclopentadiene occur in the gas phase with endo stereoselectivity, which is more pronounced with the latter diene. Endo-exo energy differences in the gas phase with the B3LYP/6-31+G(d)//B3LYP/6-31+G(d), B3LYP/6-31G(d)//B3LYP/6-31G(d), and B3LYP/6-31G(d)//AM1 model chemistries are almost the same. Solvent effects are responsible for the inversion of the stereoselectivity in the reactions of sulfilimines with furan because of the great difference in the dipole moments in endo and exo approaches.     

Structural elucidation of nitro-substituted five-membered aromatic heterocycles utilizing GIAO DFT calculations

The GIAO (Gauge Including Atomic Orbitals) DFT (Density Functional Theory) method is applied at the B3LYP/6-31+G(d,p)//B3LYP/6-31+G(d), B3LYP/6-311++G(d,p)//B3LYP/6-31+G(d), B3LYP/6-311+G (2d,p)//B3LYP/6-31+G(d) and B3LYP/6-311++G(d,p)//B3LYP/6-311++G(d,p) levels of theory for the calculation of proton and carbon chemicals shifts and coupling constants for 25 nitro-substituted five-membered heterocycles. Difference (1D NOE) spectra in combination with long-range gHMBC experiments were used as tools for the structural elucidation of nitro-substituted five-membered heterocycles. The assigned NMR data (chemical shifts and coupling constants) for all compounds were found to be in good agreement with theoretical calculations using the GIAO DFT method. The magnitudes of one-bond (1JCH) and long-range (nJCH, n>1) coupling constants were utilized for unambiguous differentiation between regioisomers of nitro-substituted five-membered heterocycles.     

Quantum chemical calculations based on ONIOM and the DFT methods in the inclusion complex: doxycycline/2-O-Me-β-cyclodextrin

A computational study of inclusion complexes of 2-methyl-βCD with Doxycycline tautomeric (enol and keto form) has been performed with several combinations of ONIOM hybrid calculations. The reliability of the ONIOM2 calculations at the integrated level, ONIOM2 (M05-2X/6-31G(d): M05-2X/3-21G*), ONIOM2 (M05-2X/6-31G(d):HF/3-21G*), ONIOM2 (B3LYP/6-31G(d):HF/3-21G*), ONIOM2 (B3LYP/6-31G(d):B3LYP/3-21G*) and ONIOM2 (B3PW91/6-31G(d):B3PW91/3-21G*) was examined. Their complexation, binding, deformation and stabilization energies, and geometrical data were compared with those of the target geometry structure optimized at the M05-2X/6-31G(d) level of theory. Mixed combinations ONIOM2 (M05-2X 6-31G(d):HF 3-21G*) and ONIOM2 (B3LYP 6-31G(d):HF 3-21G*) reproduces nearly the target geometry structure and provides realistic energetic results at a relatively low computational cost.     

Comparison of Some Computational Methods for Geometry Optimisation of Phosphorus Acid Derivatives

Several economical methods for geometry optimization, that should be applicable to larger molecules, have been evaluated for 19 phosphorus acid derivatives. MP2/cc-pVDZ geometry optimizations are used as reference points and the geometries obtained from the other methods are evaluated with respect to deviations in bond lengths and angles, from the reference geometries. The geometry optimization methods are also compared to the much used B3LYP/6-31G(d) method. Single point energies obtained by subsequent EDF1/6-31+G(d) or B3LYP/6-31+G(d,p) calculations on the respective equilibrium geometries are also reported relative to the energies obtained from the reference geometries. The geometries from HF/MIDI! optimizations were closer to those of the references than the geometries of the HF/3-21G(d), HF/6-31G(d), and B3LYP/MIDI! optimizations. The EDF1/6-31+G(d) or B3LYP/6-31+G(d,p) single point energies obtained from the HF/3-21G(d), HF/6-31G(d), and B3LYP/MIDI! geometries gave a mean absolute deviation (MAD) from that of the reference geometries of 1.4-3.9 kcal mol m 1 . The HF/MIDI! geometries, however, gave EDF1/6-31+G(d) and B3LYP/6-31+G(d,p) energies with a MAD of only about 0.5 and 0.55 kcal mol m 1 respectively from the energies obtained with the reference geometries. Thus, use of HF/MIDI! for geometry optimization of phosphorus acids is a method that gives geometries of near-MP2 quality, resulting in a fair accuracy of energies in subsequent single point calculations, at a much lower computational cost other methods that give similar accuracies.     

Neutral homoaromaticity in some heterocyclic systems

[reaction: see text] Neutral homoaromaticity has been evaluated in heterocyclic systems related to the bicyclo[3.2.1]octane skeleton with replacement of CH(2) at C-2 in bicyclo[3.2.1]octa-3,6-diene with X = BH, AlH, Be, Mg, O, S, PH, NH (12); replacement of CH at C-3 in bicyclo[3.2.1]octa-3,6-dien-2-yl anion with PH, S, NH, O (13); and replacement at C-2 and C-3 with N and O (14). Stabilization energies (SE) are evaluated using density functional theory and homodesmotic equations at the B3LYP/6-311+G(3df,2p)//B3LYP/6-31G(d) level for series 12-14. Stabilization energies are compared with diamagnetic susceptibility exaltations, Lambda, CSGT-B3LYP/6-31G(d)//B3LYP/6-31G(d), and nucleus-independent chemical shifts (NICS), GIAO-B3LYP/6-311+G(2d,p)//B3LYP/6-31G(d). Analysis of frontier orbitals and geometries, B3LYP/6-31G(d)//B3LYP/6-31G(d), and proton affinities of 2-azabicyclo[3.2.1]octa-3,6-diene, pyrrole, and related model systems, B3LYP/6-311+G(2d,2p)//B3LYP/6-31G(d), provide complementary evidence supporting the division of the substrates evaluated into antihomoaromatic (12, X = BH, AlH, and Be), nonhomoaromatic (12, X = O, S, NH, PH), and homoaromatic (13, X = S, PH, NH, O and 14 X = ON), with 12 (X = Mg) appearing as transitional between anti- and nonhomoaromatic.     

Theoretical study of the conformational energy hypersurface of cyclotrisarcosyl

The multidimensional Conformational Potential Energy Hypersurface (PEHS) of cyclotrisarcosyl was comprehensively investigated at the DFT (B3LYP/6-31G(d), B3LYP/6-31G(d,p) and B3LYP/6-311++G(d,p)), levels of theory. The equilibrium structures, their relative stability, and the Transition State (TS) structures involved in the conformational interconversion pathways were analyzed. Aug-cc-pVTZ//B3LYP/6-311++G(d,p) and MP2/6-31G(d)//B3LYP/6-311++G(d,p) single point calculations predict a symmetric cis-cis-cis crown conformation as the energetically preferred form for this compound, which is in agreement with the experimental data. The conformational interconversion between the global minimum and the twist form requires 20.88 kcal mol-1 at the MP2/6-31G(d)//B3LYP/6-311++G(d,p) level of theory. Our results allow us to form a concise idea about the internal intricacies of the PEHSs of this cyclic tripeptide, describing the conformations as well as the conformational interconversion processes in this hypersurface. In addition, a comparative analysis between the conformational behaviors of cyclotrisarcosyl with that previously reported for cyclotriglycine was carried out     

Theoretical investigation on the GaCl3-catalyzed ring-closing metathesis reaction of N-2,3-butadienyl-2-propynyl-1-amine: three-membered ring versus four-membered ring mechanism

The gallium chloride (GaCl(3))-catalyzed ring-closing metathesis reaction mechanism of N-2,3-butadienyl-2-propynyl-1-amine has been studied at the Becke three-parameter hybrid functional combined with Lee-Yang-Parr correlation functional (B3LYP)/6-31G(d), B3LYP/6-31+G(d,p), B3LYP/6-311++G(d,p)//B3LYP/ 6-31G(d) and the second-order M?ller-Plesset perturbation (MP2)/6-311++G(d,p)//B3LYP/6-31+G(d,p) levels. It was found that the final metathesis product can be yielded via a three-membered or four-membered ring mechanism. The three-membered ring pathway is favorable due to its low energy barrier at the rate determining step. The whole reaction is stepwise and strongly exothermic.     

Conformational space of 4,4′-methoxypropylstilbene molecule

Quantum-chemical simulation of the ground state [the density function B3LYP/6-31G, B3LYP/6-31G(d), and B3LYP/6-31+G(d,p) and the perturbation theory MP2/6-31G(d) methods] and the transition states [the B3LYP/6-31G(d) method] of 4,4′-methoxypropylstilbene molecule has been performed. Using the Ellinger MM2 force field method, the potentials of internal rotation have been obtained for each rotational degree of freedom of the molecule. The B3LYP simulation has revealed the planarity of the conjugated system and the orthogonal position of the alkyl substituent, whereas the benzene rings have deviated by about 20° with respect to the double bond plane according to the MP2 data. Three transition states of the molecule corresponding to the saddle points of the first and the second orders have been revealed. The stationary points have been identified by means of vibrational analysis.     

DFT study on mechanism of the classical Biginelli reaction

The condensation of benzaldehyde, urea, and ethyl acetoacetate according to the procedure described by Biginelli was investigated at the B3LYP/6-31G（d）, B3LYP/6-31 ＋G（d,p）, and B3LYP/6-311＋G（3df,2p）//B3LYP/6-31＋G（d,p） levels to explore the reaction mechanism. According to the mechanism proposed by Kappe, structures of five intermediates were optimized and four transition states were found. The calculation results proved that the mechanism proposed by Kappe is right.     

A density functional theory study of the structure and vibrational spectra of beta-carotene, capsanthin, and capsorubin

A theoretical study of the structure and the vibrational spectra of the beta-carotene molecule and its derivatives capsanthin and capsorubin is carried out. We first investigate systematically the theoretical method which provides the best results for beta-carotene by performing ab initio calculations at the HF/6-31G(d), SVWN/6-31G(d), PBE0/6-31G(d), BLYP/6-31G(d), B3LYP/6-31G(d), B3LYP/6-31G(d,p), B3LYP/6-311G(d), and B3LYP/6-311G(d,p) levels and by using previous theoretical results available in the literature obtained at the AM1 and BPW91/6-31G(d) levels. The influence of both the level of calculation and the size of the basis set used in the geometry optimization and in the determination of the IR and Raman spectra of this molecule is thus analyzed. It is confirmed that the hybrid functional B3LYP with the basis 6-31G(d) is the method that gives the best results as a whole. By use of this level of calculation, we next optimize the molecular geometries of related molecules of capsanthin and capsorubin, which to the best of our knowledge have only been studied at the semiempirical AM1 level. In addition we calculate the IR and Raman spectra of these molecules at the B3LYP/6-31G(d) level of theory. The results obtained for capsanthin show on the one hand that the double bond of the beta-ionone ring is outside the polyene chain plane, due to the repulsion between the hydrogen atoms of the ring methyl groups and the hydrogen atoms of the polyene chain, and on the other hand that the carbonyl double bond in the other headgroup is very close to planarity with the polyene chain, since in this case such a repulsion does not exist. For the molecule of capsorubin the two carbonyl groups also take the same coplanar orientation relative to the polyene chain. The IR and Raman spectra theoretically computed for these two molecules are finally compared with their experimental spectra and the vibrational normal modes of the main signals are interpreted.     

Redox potentials of dopamine and its supramolecular complex with aspartic acid

Dopamine (DA) can be oxidized to dopamine quinone (DAquinone) through a one-step, two-electron redox reaction. The electron transfer property of DA and its supramolecular complex with aspartic acid (Asp) has been investigated by the theoretical calculations. We calculated the standard redox potentials (E ^o) of DA/DAquinone at the MP2/6-31G(d,p)//B3LYP/6-31G(d,p), MP2/6-31+G(d,p)//B3LYP/6-31+G(d,p), MP2/6-31G(d,p)//B3LYP/6-311G(d,p), and MP2/6-311+G(d,p)//B3LYP/6-311+G(d,p) levels. Comparing the experimental value, the redox potentials of DA/DAquinone obtained at MP2//B3LYP/6-311G(d,p) and MP2//B3LYP/6-311+G(d,p) levels can be considered as the upper and lower estimates. DA can form supramolecular complex (DA-Asp) with Asp through hydrogen bond (H-bond). Therefore, the values of 0.631 and 0.628 V obtained at MP2//B3LYP/6-311G(d,p) and MP2//B3LYP/6-311+G(d,p) levels for DA-Asp/DAquinone-Asp can be proposed as the upper and lower estimates of a probable (about 0.630 V) value of the corresponding redox potential. The calculated E ^o values of DA-Asp/DAquinone-Asp at the four theoretical levels are upper than those of DA/DAquinone, which indicates that the formation of H-bonds weaken the electron-donating ability of DA.     

苯并氧化呋咱稳定性和异构化的DFT和ab initio研究

运用B3LYP/6-31G(d)密度泛函理论（DFT）方法对苯并氧化呋咱、邻二亚硝基苯及其间的异构化反应进行了计算研究。结果表明,苯并氧化呋咱的分子总能量比邻二亚硝基苯的低;由苯并氧化呋咱异构为邻二亚硝基苯的正向反应活化能（Ea+=51.0kJ/mol）,与文献实测值（58.6kJ/mol）较接近,而其逆向反应活化能（Ea-=4.6kJ/mol）很小,从而揭示了苯并氧化呋咱比邻二亚硝基苯更稳定·此外,进行了HF/3-21G、HF/6-31G(d)和MP2/6-31G(d)//6-31G(d)水平下相应的计算,发现B3LYP-DFT的结果较abinitio为优。谐振动频率的B3LYP/6-31G(d)计算还支持了邻二亚硝基苯为苯并氧化呋咱“自-自”互变重排反应的中间体。     

An experimental and theoretical study of the basicity of tetra-tert-butyltetrahedrane

The gas-phase basicity (GB) of tetra-tert-butyltetrahedrane (tBu4THD) was determined by FT-ICR mass spectrometry and comparison with reference compounds of known basicity. Its GB, 1035+/-10 kJ x mol(-1), makes tetra-tert-butyltetrahedrane one of the strongest bases reported so far. Ab initio calculations [B3LYP/6-31G(d) and B3LYP/6-311 + G(d,p)//6-31G(d)] have been carried out in order to compare the high experimental basicity of tBu4THD with that estimated theoretically. Both B3LYP/6-31G(d) and QCISD(T) calculations were used to determine the reaction path which connects the initial tetrahedrane-ammonium complex with the final products, protonated cyclobutadiene (CBDH+) and ammonia.     

1,3,5-三硝基苯在高岭石表面吸附的理论研究

构造了高岭石硅氧层和铝氧层的团簇模型(分别为Si₁₃O₃₇H₂₂和Al₆O₂₄H₃₀), 并分别在B3LYP/6-31G(d), MP2/6-31G(d)//B3LYP/6-31G(d)和B3LYP/6-311++G(d,p)//B3LYP/6-31G(d)理论水平上对1,3,5-三硝基苯(TNB)在高岭石表面的吸附性质(如优化的几何构型、 结构参数、 吸附能、 振动频率、 静电势和分子轨道等)进行了研究. 结果表明, TNB和硅氧层表面间的相互作用以静电和范德华相互作用为主; TNB与铝氧层间的相互作用以氢键为主, 且TNB和铝氧层间相互作用的能量更低, 结构更稳定.     

DFT study of the electronic and charge transfer properties of perfluoroarene–thiophene oligomers

The ground state structures of 5,5″-diperfluorophenyl-2,2′:5′,2″:5″,2‴-quaterthiophene (1), 5,5′-bis{1-[4-(thien-2-yl)perfluorophenyl]}-2,2′-dithiophene (2), 4,4′-bis[5-(2,2′-dithiophenyl)]-perfluorobiphenyl (3), 5,5″-diperfluorophenyl-2,2′:5′,2″-tertthiophene (4), 5,5′-diperfluorophenyl-2,2′-dihiophene (5), and 5,5-diperfluorophenylthiophene (6) have been optimized at the B3LYP/6-31G(d), B3LYP/6-31G(d,p), PBE0/6-31G(d), and PBE0/6-31G(d,p) level of theories. The B3LYP/6-31+G(d) and PBE0/6-31+G(d) level of theories have been applied to investigate the absorption spectra. The PBE0 functional is good to predict the C–S bond lengths while the C–F bond lengths are good envisaged with B3LYP functional. The increment of thiophene rings between two perfluoroarene rings leads to red shift in absorption spectra. The electron affinities are energetically destabilized while energetic stabilization of the radical-cation increases by decreasing the thiophene rings from four to one. The perfluoroarene rings leads to enhance the electron injection.     

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°）时,振动频率的计算值和实验观察值非常吻合。     

A new scheme for determining the intramolecular seven-membered ring N-H...O=C hydrogen-bonding energies of glycine and alanine peptides

In this paper a new scheme was proposed to calculate the intramolecular hydrogen-bonding energies in peptides and was applied to calculate the intramolecular seven-membered ring N-H...O=C hydrogen-bonding energies of the glycine and alanine peptides. The density-functional theory B3LYP6-31G(d) and B3LYP6-311G(d,p) methods and the second-order Moller-Plesset perturbation theory MP26-31G(d) method were used to calculate the optimal geometries and frequencies of glycine and alanine peptides and related structures. MP26-311++G(d,p), MP26-311++G(3df,2p), and MP2/aug-cc-pVTZ methods were then used to evaluate the single-point energies. It was found that the B3LYP6-31G(d), MP26-31G(d), and B3LYP6-311G(d,p) methods yield almost similar structural parameters for the conformers of the glycine and alanine dipeptides. MP2/aug-cc-pVTZ predicts that the intramolecular seven-membered ring N-H...O=C hydrogen-bonding strength has a value of 5.54 kcal/mol in glycine dipeptide and 5.73 and 5.19 kcal/mol in alanine dipeptides, while the steric repulsive interactions of the seven-membered ring conformers are 4.13 kcal/mol in glycine dipeptide and 6.62 and 3.71 kcal/mol in alanine dipeptides. It was also found that MP26-311++G(3df,2p) gives as accurate intramolecular N-H...O=C hydrogen-bonding energies and steric repulsive interactions as the much more costly MP2/aug-cc-pVTZ does.