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.     

CH2F自由基与HNCO反应机理的理论研究

The reaction mechanism of CH2F radical with HNCO was investigated by density functional theory (DFT)at the B3LYP/6-311++G(d,p) level. The geometries of the reactants, the intermediates, the transition states and the products were optimized. The transition states were verified through the vibration analysis.The relative energies were calculated at the QCISD(T)/6-311++G**//B3LYP/6-311++G(d,p) level. Seven feasible reaction pathways of the reaction were studied. The results indicate that the pathway (5) is the most favorable to occur, so it is the main pathway of the reaction.     

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     

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.     

The Dowd-Beckwith ring expansion: a theoretical study

UB3LYP/6-311++G(d,p) and ROMP2/6-311++G(d,p)//UB3LYP/6-311++G(d,p) calculations including the effect of benzene solvent through the PCM-UAHF method render a concerted mechanism without fragmentation as the most favourable one for the Dowd-Beckwith radical ring expansion of the bromomethyl adduct of methyl cyclopentanone-2-carboxylate to yield methyl cyclohexanone-3-carboxylate. The corresponding concerted TS is a bicyclic alcoxy radical.     

Sugars in the gas phase: The conformational properties of erythrose, threose, and erythrulose characterized by quantum chemistry calculations

The structural properties of the three open chain C₄H₈O₄ sugars, i.e. two aldoses (erythrose and threose) and one ketose (erythrulose), have been investigated by DFT and ab initio calculations to get accurate structures and relative energies. The structure of all the conformers predicted within 10 kJ/mol has been optimized at the B3LYP/6-311++G(d,p) level of the theory. Two types of intramolecular hydrogen bonds have been clearly identified. They are related to the hydroxyl and to the carbonyl oxygen atoms and are of weak and middle strength, respectively. The most stable structures have been optimized at the B3LYP/6-311++G(2df,p) and at the MP2/6-311++G(2df,p) levels of the theory in order to calculate accurate rotational parameters and dipole moment for their future detection in the microwave range in the gas phase. Their corresponding harmonic IR spectra have also been calculated and their fingerprint signature is discussed in the region of the OH stretching vibrations, of the torsion of the C–O bonds and of the deformation of the C–C skeleton.     

Inductive effects in radicals calculated from DFT energies; substituted bicyclo[2.2.2]octan-1-yloxy radicals

Energies of a series of 4-substituted 1-oxybicyclo[2.2.2]octan-1-yloxy radicals with 18 various substituents were calculated within the framework of the DFT theory at the levels UB3LYP/6-311+G(d,p)//UB3LYP/6-311+G(d,p) and UB3LYP/6-311++G(2df,p)//UB3LYP/6-311+G(d,p) and compared with similar series of the parent alcohols, their deprotonated and protonated forms calculated at the levels B3LYP/6-311+G(d,p)//B3LYP/6-311+G(d,p) and B3LYP/6-311++G(2df,p)//B3LYP/6-311+G(d,p). The two levels are of the same performance and both are sufficient for molecules of this type according to comparison with scarce experimental gas-phase acidities and basicities. The substituent effects were analyzed in terms of isodesmic equations. In addition to strong dependence on the substituent inductive effect, a slight dependence on the electronegativity of the first atom of the substituent was proven in certain cases. In all aspects, there is no qualitative difference between the effects on radicals and on similar closed shell species. Radicals behave as slightly electron deficient; the substituent effect is weaker than that on the ions but stronger than on neutral molecules.     

Nucleofugality of phenyl and methyl carbonates

A series of X,Y-substituted benzhydryl phenyl carbonates 1 and X,Y-substituted benzhydryl methyl carbonates 2 were subjected to solvolysis in different methanol/water, ethanol/water, and acetone/water mixtures at 25 degrees C. The LFER equation, log k = sf(Ef + Nf), was used to derive the nucleofuge-specific parameters (Nf and sf) for phenyl carbonate (1LG) and methyl carbonate (2LG) leaving groups in a given solvent in SN1 type reaction. Kinetic measurements showed that phenyl carbonates solvolyze one order of magnitude faster than methyl carbonates. Optimized geometries of 1LG and 2LG at B3LYP/6-311G(d,p), B3LYP/6-311++G(d,p), and MP2(full)/6-311++G(d,p) levels revealed that negative charge delocalization in carbonate anions to all three oxygen atoms occurs due to negative hyperconjugation. Phenyl carbonate (1LG) is a better leaving group (Nf = -0.84 +/- 0.07 in 80% v/v aq EtOH) than methyl carbonate 2LG (Nf = -1.84 +/- 0.07 in 80% v/v aq EtOH) because of more pronounced negative hyperconjugation, which is characterized with a more elongated RO-C bond and more increased RO-C-CO angle in 1LG than in 2LG. Calculated affinities of benzhydryl cation toward methyl and phenyl carbonate anions (DeltaDeltaEaff = 11.7 kcal/mol at the B3LYP/6-311++G(d,p) level and DeltaDeltaEaff = 2.7 kcal/mol at the PCM-B3LYP/6-311++G(d,p) level in methanol, respectively) showed that 1LG is more stabilized than 2LG, which is in accordance with greater solvolytic reactivity of 1 than 2.     

Protonated nitro group: structure, energy and conjugation

Structure of protonated nitro compounds was investigated by calculations at the levels MP2(FC)/6-311++G(2d,2p)//MP2(FC)/6-311++G(2d,2p)(nitromethane and reference compounds) or B3LYP/6-311+G(d,p)//B3LYP/6-311+G(d,p)(nitrobenzene and its 18 meta- and para-substituted derivatives). The group NO2H+ reveals many similarities with the isoelectronic group CO2H as the preferred conformation, conformational equilibrium, and stabilization by interaction (resonance) within the group quantified by means of isodesmic reactions. However, there is a difference in the interaction with donor groups (for instance in 4-nitroaniline) that is much stronger with NO2H+ than with CO2H. This interaction may be called resonance and may be described by standard resonance formulas, but these formulas predict only partially the geometry and cannot explain the great interaction energy.     

Theoretical studies on the intramolecular hydrogen bond and tautomerism of 8-mercaptoquinoline in the gaseous phase and in solution using modern DFT methods

A theoretical quantum chemical study of the intramolecular hydrogen bonding interactions in 8-mercaptoquinoline has been carried out. Special attention has been paid to the rotation of S-H bond and intramolecular proton-transfer reactions. Therewith, the B3LYP/6-311++G(d,p), B3LYP/6-31+G(2d,2p), MPW1K/6-311++G(d,p), MPW1K/6-31+G(2d,2p), BH&HLYP/6-311++G(d,p), and G96LYP/6-311++G(d,p) methods have been used. By means of the Onsager and PCM reaction field methods, the effects of solvent on hydrogen-bond energies, conformational equilibria, rotational barriers, and tautomerism in aqueous solution have been studied. These simulations were done at the MPW1K/6-311++G(d,p) and B3LYP/6-311++G(d,p) levels. Natural-bond orbital analysis has been performed to study the intramolecular hydrogen bond (IHB) in the gaseous phase and in aqueous medium. The stability of forms under consideration in solution does not coincide with that in the gaseous phase, underlining a great importance of the electrostatic influence of solvent. Double-proton transfer in the prototropic tautomerization of 8-mercaptoquinoline, one water molecule complex in the gaseous phase and in solution, has been systematically studied. The double-proton transfer occurs concertedly and synchronously. The water-assisted tautomerization is kinetically less, but thermodynamically more favorable, compared to that of the single-proton transfer. As in the case with single-proton transfer, for water-assisted reaction, the tautomerization energies and barrier heights decrease with the increase in dielectric constant, which implies faster and more complete tautomerization of 8-mercaptoquinoline in a polar solvent.     

FT-IR and FT-Raman spectroscopic investigation, computed vibrational frequency analysis and IR intensity and Raman activity peak resemblance analysis on 2-nitroanisole using HF and DFT (B3LYP and B3PW91) calculations

Fourier-transform Raman and infrared spectra of 2-nitroanisole are recorded (4000-100 cm(-1)) and interpreted by comparison with respective theoretical spectra calculated using HF and DFT method. The geometrical parameters with C(S) symmetry, harmonic vibrational frequencies, infrared and Raman scattering intensities are determined using HF/6-311++G (d, p), B3LYP/6-311+G (d, p), B3LYP/6-311++G (d, p) and B3PW91/6-311++G (d, p) level of theories. A detailed vibrational spectral analysis has been carried out and assignments of the observed fundamental bands have been proposed on the basis of peak positions and relative intensities. The results of the calculations have been used to simulate IR and Raman spectra for the molecule that showed good agreement with the observed spectra. The SQM method, which implies multiple scaling of the DFT force fields has been shown superior to the uniform scaling approach. The vibrational frequencies and the infrared intensities of the C-H modes involved in back-donation and conjugation are also investigated.     

On the flexibility of beta-peptides

The full conformational space was explored for an achiral and two chiral beta-peptide models: namely For-beta-Ala-NH2, For-beta-Abu-NH2, and For-beta-Aib-NH2. Stability and conformational properties of all three model systems were computed at different levels of theory: RHF/3-21G, B3LYP/6-311++G(d,p)//RHF/3-21G, B3LYP/6-311++G(d,p), MP2//B3LYP/6-311++G(d,p), CCSD//B3LYP/6-311++G(d,p), and CCSD(T)//B3LYP/6-311++G(d,p). In addition, ab initio E = E(phi, micro, psi) potential energy hypersurfaces of all three models were determined, and their topologies were analyzed to determine the inherent flexibility properties of these beta-peptide models. Fewer points were found and assigned than expected on the basis of Multidimensional Conformational Analysis (MDCA). Furthermore, it has been demonstrated, that the four-dimensional surface, E = E(phi, mu, psi), can be reduced into a three-dimensional one: E = E[phi, f(phi), psi]. This reduction of dimensionality of freedom of motion suggests that beta-peptides are less flexible than one would have thought. This agrees with experimental data published on the conformational properties of peptides composed of beta-amino acid residues.     

Reaction of ketenyl radical with acetylene: a promising route for cyclopropenyl radical

The reaction of the ketenyl radical (HCCO) with acetylene (C(2)H(2)) is very relevant to the oxygen-acetylene flames and fuel-rich combustion process for nitrogen-containing compounds. Unfortunately, except for several rate constant measurements, the mechanism is completely unknown for this reaction. In this paper, detailed theoretical investigations are performed for the HCCO + C(2)H(2) reaction at the G3B3 level using the B3LYP/6-31G(d), B3LYP/6-311++G(d,p), and QCISD/6-31G(d) geometries. The exclusive fragmentation channel is the formation of the cyclopropenyl radical (c-C(3)H(3)) and carbon monoxide (CO) via the chainlike OCCHCHCH and three-membered ring OC-cCHCHCH intermediates. Thus, the mass spectroscopic peak of C(3)H(3)(+) in a previous experiment can be explained. The calculated overall reaction barrier is 4.4, 4.4, and 5.3 kcal/mol at the G3B3//B3LYP/6-31G(d), G3B3//B3LYP/6-311++G(d,p), and G3B3//QCISD/6-31G(d) levels, respectively. The title reaction may provide an effective route for generating the long-sought cyclopropenyl radical in the laboratory, which has been the long-standing subject of numerous theoretical studies as the simplest cyclic conjugate radical, and its bulky derivatives were already known. Future experimental investigations for the HCCO + C(2)H(2) reaction are greatly desired to test the predicted fragmentation channel. The implication of the present study in combustion and interstellar processes is discussed.     

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.     

Theoretical and experimental study of the absolute configuration of helical structure of (2R,3S)-Rubiginone A2 analog

A novel analogue of (2R,3S)-Rubiginone A₂ was synthesized as a chiral helical model compound via an eight-step procedure (2.7% overall yield). Quantum methods, such as density functional theory (DFT) at different basis sets of 6-311+(d), 6-311++G(2d,p), were used to compute its optical rotation and electronic circular dichroism at the B3LYP/6-311++G(2d,p) level in the gas phase and in solution using PCM model, respectively. UV corrections were performed in electronic circular dichroism (ECD) simulations to match the experimental ECD well. The suitable computational methods, e.g., B3LYP/6-311++G(2d,p)//B3LYP/6-311++G(2d,p) in the gas phase using zero-point energy in Boltzmann statistics, were found and suggested for optical rotation and circular dichroism computations that can be used for absolute configuration determination of chiral helical compounds.     

Molecular structure and vibrational spectra of 3-chloro-4-fluoro benzonitrile by ab initio HF and density functional method

In this work, the experimental and theoretical spectra of 3-chloro-4-fluoro benzonitrile (3C4FBN) were studied. The Fourier transform infrared and Fourier transform Raman spectra of 3C4FBN were recorded in the solid phase. The optimized geometry was calculated by HF and B3LYP methods with 6-311++G(d,p) basis set. The harmonic-vibrational frequencies, infrared intensities and Raman scattering activities of the title compound were performed at and HF/B3LYP/6-311++G(d,p) level of theories. The scaled theoretical wave number showed very good agreement with the experimental values. The thermodynamic functions of the title compound was also performed at HF/6-31G(d,p) and B3LYP/6-311++G(d,p) level of theories. A detailed interpretation of the infrared and Raman spectra of 3C4FBN was reported. The theoretical spectrograms for FT-IR and FT-Raman spectra of the title molecule have been constructed.     

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.     

OH-与CH2ClF反应的阴离子产物通道

理论研究了羟基负离子(OH^-)与氟氯代甲烷(CH₂ClF)反应的阴离子产物通道. 分别在B3LYP/6-31+G(d,p)和B3LYP/6-311++G(2d,p)水平上得到反应势能面上各关键物种的优化构型, 进而计算得到谐振频率和零点能. 基于CCSD(T)/6-311+G(3df,3dp)水平的相对能量, 描述了由质子转移和双分子亲核取代(S_N2)过程生成各阴离子产物的途径. 各阴离子产物途径势垒的计算结果表明质子转移过程是实验中的主要产物通道, 与以往实验测量的结论相符. 此外, 计算还显示双分子亲核取代过程得到了非典型的阴离子产物, 其中动力学效应可能会导致F^-的生成.     

Topological study of intramolecular hydrogen bonding in beta-hydroxyethylperoxy radical and beta-hydroxyethoxy radical along its dissociation pathway

The presence of intramolecular hydrogen bond (IHB) in beta-hydroxyethylperoxy and beta-hydroxyethoxy radicals was investigated using the QTAIM topological study on UB3LYP/6-311++G(2d,2p) charge densities. Only one of the two conformers of beta-hydroxyethylperoxy radical which were previously considered to present IHB displays a bond critical point (BCP) associated to an IHB. Furthermore, the atomic energies and electron populations indicate no evidence of IHB in the second conformer. Nevertheless, very small differences in molecular energies were obtained using several one-step and multi-step methods (G3, G3B3) between both conformers. No BCP is found between the hydroxyl hydrogen and the oxygen in the most stable conformer or in the transition state for the dissociation path of beta-hydroxyethoxy radical. However, a BCP is formed in the last steps of this path, thereby yielding H-bonded products.     

戊烯自由基阳离子的密度泛函理论研究

使用密度泛函理论B3LYP方法和6-31G(d,p)、6-31＋G(d,p)、6-311G(d,p)及6-311＋G(d,p)基组，分别对2-C5H10＋和1-C5H10＋的各种构象进行了几何构型优化，并用B3LYP/6-311G(d,p)进行了频率分析计算.计算预言1-C5H10＋具有非平面构型，与以往报导的从头算计算结论相反.在两个自由基阳离子的各种构象的B3LYP几何构型上，进行了B3LYP和UMP2(full)方法的超精细偶合常数计算，得到了比以往更好的结果.