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.     

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.     

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.     

Modern valence-bond description of chemical reaction mechanisms: the 1,3-dipolar addition of diazomethane to ethene

The electronic mechanism for the gas-phase concerted 1,3-dipolar cycloaddition of diazomethane (CH2N2) to ethene (C2H4) is described through spin-coupled (SC) calculations at a sequence of geometries along the intrinsic reaction coordinate obtained at the MP2/6-31G(d) level of theory. It is shown that the bonding rearrangements occurring during the course of this reaction follow a heterolytic pattern, characterized by the movement of three well-identifiable orbital pairs, which are initially responsible for the pi bond in ethene and the C-N pi bond and one of the N-N pi bonds in diazomethane and are retained throughout the entire reaction path from reactants to product. Taken together with our previous SC study of the electronic mechanism of the 1,3-dipolar cycloaddition of fulminic acid (HCNO) to ethyne (C2H2) (Theor. Chim. Acc. 1998, 100, 222), the results of the present work suggest strongly that most gas-phase concerted 1,3-dipolar cycloaddition reactions can be expected to follow a heterolytic mechanism of this type, which does not involve an aromatic transition state. The more conventional aspects of the gas-phase concerted 1,3-dipolar cycloaddition of diazomethane to ethene, including optimized transition structure geometry, electronic activation energy, activation barrier corrected for zero-point energies, standard enthalpy, entropy and Gibbs free energy of activation, have been calculated at the HF/6-31G(d), B3LYP/6-31G(d), MP2/6-31G(d), MP2/6-31G(d,p), QCISD/6-31G(d) and CCD/6-31G(d) levels of theory. We also report the CCD/6-311++G(2d, 2p)//CCD/6-31G(d), MP4(SDTQ)/6-311++G(2d,2p)//CCD/6-31G(d) and CCSD(T)/6-311++G(2d, 2p)//CCD/6-31G(d) electronic activation energies.     

Structure and conformational stability of CH2CHCH2X (X=F, Cl and Br) molecules—post Hartree–Fock and density functional theory methods

The molecular structure and conformational stability of CH₂CHCH₂X (X=F, Cl and Br) molecules were studied using ab initio and density functional theory (DFT) methods. The molecular geometries of 3-fluoropropene were optimized employing BLYP and B3LYP levels of theory of DFT method implementing 6-311+G(d,p) basis set. The MP2/6-31G^*, BLYP and B3LYP levels of theory of ab initio and DFT methods were used to optimize the 3-chloropropene and 3-bromopropene molecules. The structural and physical parameters of the molecules are discussed with the available experimental values. The rotational potential energy surface of the above molecules were obtained at MP2/6-31G^* and B3LYP/6-311+G(d,p) levels of theory. The Fourier decomposition of the rotational potentials were analyzed. The HF/6-31G^* and MP2/6-31G^* levels of theory have predicted the cis conformer as the minimum energy structure for 3-fluoropropene, which is in agreement with the experimental values, whereas the BLYP/6-311+G(d,p) and B3LYP/6-311+G(d,p) levels of theory reverses the order of conformation. The ΔE values calculated for 3-chloropropene at MP2/6-31G^*, BLYP/6-311+G(d,p) and B3LYP/6-311+G(d,p) levels of theory show that the gauche form is more stable than the cis form, which is in agreement with the experimental value. The same levels of theory have also predicted that the gauche form is stable than cis for 3-bromopropene molecule. The maximum hardness principle has been able to predict the stable conformer of 3-fluoropropene at HF/6-31G^* level of theory, but the same level of theory reverses the conformational stability of 3-chloropropene and 3-bromopropene molecules and MP2/6-31G^* level of theory predicted the stable conformer correctly.     

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.     

Addition reaction of adamantylideneadamantane with Br2 and 2Br2: a computational study

Ab initio calculations were carried out for the reaction of adamantylideneadamantane (Ad=Ad) with Br2 and 2Br2. Geometries of the reactants, transition states, intermediates, and products were optimized at HF and B3LYP levels of theory using the 6-31G(d) basis set. Energies were also obtained using single point calculations at the MP2/6-31G(d)//HF/6-31G(d), MP2/6-31G(d)//B3LYP/6-31G(d), and B3LYP/6-31+G(d)//B3LYP/6-31G(d) levels of theory. Intrinsic reaction coordinate (IRC) calculations were performed to characterize the transition states on the potential energy surface. Only one pathway was found for the reaction of Ad=Ad with one Br2 producing a bromonium/bromide ion pair. Three mechanisms for the reaction of Ad=Ad with 2Br2 were found, leading to three different structural forms of the bromonium/Br3- ion pair. Activation energies, free energies, and enthalpies of activation along with the relative stability of products for each reaction pathway were calculated. The reaction of Ad=Ad with 2Br2 was strongly favored over the reaction with only one Br2. According to B3LYP/6-31G(d) and single point calculations at MP2, the most stable bromonium/Br3- ion pair would form spontaneously. The most stable of the three bromonium/Br3- ion pairs has a structure very similar to the observed X-ray structure. Free energies of activation and relative stabilities of reactants and products in CCl4 and CH2ClCH2Cl were also calculated with PCM using the united atom (UA0) cavity model and, in general, results similar to the gas phase were obtained. An optimized structure for the trans-1,2-dibromo product was also found at all levels of theory both in gas phase and in solution, but no transition state leading to the trans-1,2-dibromo product was obtained.     

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.     

Molecular Dipole Moment Computed with Ab Initio MKS Charges

1 INTRODUCTION Dipole moment (DM) can be measured by a variety of experimental methods[1] or computed with atomic charge distribution directly derived from molecular orbital calculations[2~10]. The average errors in DM computed with AM1, PM3 and MNDO methods for 125 organic compounds are 0.35 D, 0.45 D and 0.45 D, respectively[2]. Feller et al. reported that the average errors in the computed DM at the HF level are 0.65 D, 0.49 D and 0.30 D with the basis set STO-3G, 3-21G and…     

Computational study of the deamination reaction of cytosine with H2O and OH-

The mechanism for the deamination reaction of cytosine with H(2)O and OH(-) to produce uracil was investigated using ab initio calculations. Optimized geometries of reactants, transition states, intermediates, and products were determined at RHF/6-31G(d), MP2/6-31G(d), and B3LYP/6-31G(d) levels and for anions at the B3LYP/6-31+G(d) level. Single-point energies were also determined at B3LYP/6-31+G(d), MP2/GTMP2Large, and G3MP2 levels of theory. Thermodynamic properties (DeltaE, DeltaH, and DeltaG), activation energies, enthalpies, and free energies of activation were calculated for each reaction pathway that was investigated. Intrinsic reaction coordinate analysis was performed to characterize the transition states on the potential energy surface. Two pathways for deamination with H(2)O were found, a five-step mechanism (pathway A) and a two-step mechanism (pathway B). The activation energy for the rate-determining steps, the formation of the tetrahedral intermediate for pathway A and the formation of the uracil tautomer for pathway B, are 221.3 and 260.3 kJ/mol, respectively, at the G3MP2 level of theory. The deamination reaction by either pathway is therefore unlikely because of the high barriers that are involved. Two pathways for deamination with OH(-) were also found, and both of them are five-step mechanisms. Pathways C and D produce an initial tetrahedral intermediate by adding H(2)O to deprotonated cytosine which then undergoes three conformational changes. The final intermediate dissociates to product via a 1-3 proton shift. Deamination with OH(-), through pathway C, resulted in the lowest activation energy, 148.0 kJ/mol, at the G3MP2 level of theory.     

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.     

Investigation of structural and biological properties of N-heterocyclic carbene silver(I) and palladium(II) complexes

abstract

Computational investigations were done on bis(1-allyl-3-benzyl-2,3-dihydro-1H-benzo[d]imidazol-2-yl)silver(I), bis(1-benzyl-3-butyl-2,3-dihydro-1H-benzo[d]imidazol-2-yl)silver(I), bis(1-allyl-3-benzyl-2,3-dihydro-1H-benzo[d]imidazol-2-yl)dibromidepalladium(II), and bis(1-benzyl-3-butyl-2,3-dihydro-1H-benzo[d]imidazol-2-yl)dibromidepalladium(II) complexes. Related complexes were optimized at different six calculation levels which are HF/6-31G(LANL2DZ), HF/6-31G(d,p)(LANL2DZ), B3LYP/6-31G(LANL2DZ), B3LYP/6-31G(d,p)(LANL2DZ), M062X/6-31G(LANL2DZ) and M062X/6-31G(d,p)(LANL2DZ) levels in vacuo. IR and NMR spectrum are calculated and examined in detail. Energy diagram of molecular orbitals, contour diagram of frontier molecular orbitals, molecular electrostatic potential maps and the harmonic surface of related molecules are examined in detail. Finally, interactions between mentioned complexes and related proteins (1BNA, 1JNX, and 2ING) are investigated in detail. As a result, it is found that biological and anti-cancer properties of silver N-heterocyclic carbene complexes are higher than those of palladium complexes.     

FT-Raman and FT-IR spectra, ab initio and density functional studies of 2-amino-4,5-difluorobenzoic acid

The Fourier transform Raman and Fourier transform infrared spectra of 2-amino-4,5-difluorobenzoic acid (2A45DFBA) were recorded in the solid phase. The equilibrium geometry, harmonic vibrational frequencies, infrared intensities and Raman scattering activities, depolarization ratios were calculated by HF and density functional B3LYP method with the 6-31+G(d, p) and 6-311+G(d, p) basis sets. The scaled theoretical wavenumbers showed very good agreement with the experimental values. The thermodynamic functions of the title compound were also performed at HF/6-31+G(d, p)/6-311+G(d, p) and B3LYP/6-31+G(d, p)/6-311+G(d, p) levels of theory. A detailed interpretations of the infrared and Raman spectra of 2-amino-4,5-difluorobenzoic acid is reported. The theoretical spectrograms for FT-IR spectra of the title molecule have been constructed.     

Gas electron diffraction study of the geometric structure of triallylborane molecule

The main structural parameters of the triallylborane molecule having the C ₃ symmetry were determined by gas electron diffraction and quantum-chemical calculations at the MP2/6-31G(d,p) and B3LYP/6-31G(d,p) levels. The parameters calculated by the MP2/6-31G(d,p) method are in better agreement with the experimental data than those calculated by the B3LYP/6-31G(d,p) method.__________Published in Russian in Izvestiya Akademii Nauk. Seriya Khimicheskaya, No. 1, pp. 98–101, January, 2005.     

Vibrational spectroscopy investigation using ab initio and density functional theory analysis on the structure of 3-aminobenzotrifluoride

In this work, we will report a combined experimental and theoretical study on molecular and vibrational structure of 3-aminobenzotrifluoride. The FT-Raman and Fourier transform infrared spectra of 3-aminobenzotrifluoride (3ABTF) were recorded in the liquid phase. The equilibrium geometry, harmonic vibrational frequencies, infrared intensities and Raman scattering activities, depolarization ratios, reduced masses were calculated by HF and density functional B3LYP method with the 6-31G(d,p) and 6-311G(d,p) basis sets. The scaled theoretical wavenumbers showed very good agreement with the experimental values. The thermodynamic functions of the title compound were also performed at HF/6-31G(d,p)/6-311G(d,p) and B3LYP/6-31G(d,p)/6-311G(d,p) levels of theory. A detailed interpretations of the infrared and Raman spectra of 3ABTF is reported. The theoretical spectrograms for FT-IR spectra of the title molecule have been constructed.     

低硝酸盐离子N2O22-和亚硝酸HONO的异构体及其异构化反应的量化研究

采用Gaussian-98程序进行，在HF／6—31G(d)，B3LYP／6—31G(d)和MP2／6—31G(d)水平下优化分子结构并寻找过渡态，对于MP2／6—31G(d)结果在QCISD(T，E4T)，MP4／6—311 G(d，p)，MP4／6—311 G(2df，p)水平下重新计算能级．并用内禀反应坐标(IRC)法研究了N2O2^2-和亚硝酸HONO的异构化反应机理。     

The photoelectron spectrum and conformation of phenylphosphine and phenylarsine

He(I) and He(II) photoelectron spectra of phenylphosphine and phenylarsine have been investigated and assigned. The rotational barrier of the phosphino group has been investigated at the MP2/6-31G(d,p)//MP2/6-31G(d,p) and HF/6-31G(d,p)//HF/6-31G(d,p) levels of theory, and that of the arsino group at the HF/6-31G(d,p)//6-31G(d,p) levels of theory. The rotational barrier of the two molecules is nearly the same. The energy difference between the two possible conformers of the molecules is low (1.5 kJ/mol at the MP2/6-31G(d,p) level of theory), allowing nearly free rotation about the P-C bond. The photoelectron spectrum cannot be interpreted by considering the most stable rotamer, but all possible conformers should be taken into account. The present interpretation is consistent with the smalln _p - interaction concluded from other investigations. The rotational barrier ofo-phosphinophenol is significantly larger than for phenylphosphine, and the photoelectron spectrum of this compound can be interpreted by considering a single conformer, and no appreciable interaction between the -system of the ring and the phosphorus lone pair.     

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和铝氧层间相互作用的能量更低, 结构更稳定.