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-GIAO calculations of 19F NMR chemical shifts for perfluoro compounds

The (19)F NMR shieldings for 53 kinds of perfluoro compounds were calculated by the B3LYP-GIAO method using the 6-31G(d), 6-31+G(d), 6-31G(d,p), 6-31++G(d,p), 6-311G(d,p), 6-311++G(d,p), 6-311G(2d,2p), 6-311++G(2d,2p), 6-311++G(2df,2p), 6-311++G(3d,2p), and 6-311++G(3df,2p) basis sets. The diffuse functions markedly reduce the difference between the calculated and experimental chemical shifts. The calculations using the 6-31++G(d,p) basis set give the chemical shifts within 10 ppm deviations from experimental values except for the fluorine nuclei attached to an oxygen atom, a four- and a six-coordinated sulfur atom, and FC(CF(3))(2) attached to a sulfur atom.     

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.     

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.     

Conformational inversion-topomerization processes of ethylcyclohexane and 1,2-dimethylcyclohexane: A computational investigation

To deeply reveal the impact of the substituents and their special orientations in ring on conformational behaviors for substituted cyclohexanes, a comprehensive study of ethylcyclohexane, cis-, and trans-1,2-dimethylcyclohexanes has been carried out. All conformational structures for them were captured by the accurate ab intio method, that is, B3LYP/6-311++G(d,p) method was used for geometry optimizations, and MP2/6-311++G(d,p), G4, and CCSD(T)/6-311++G(d,p) methods were applied for the high-level single point energy refinements. Based on CCSD(T)/6-311++G(d,p) quantum results, the conformational populations of minima for these three substituted cyclohexanes were calculated by Boltzmann distribution over 300-2500 K. Additionally, the conformational inversion-topomerization pathways for them were thoroughly investigated. The complete characterization involved in their potential energy surfaces are clearly presented by three or two-dimensional schemes.     

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.     

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.     

FT-IR, FT-Raman vibrational spectra and molecular structure investigation of 2-chloro-4-methylaniline: A combined experimental and theoretical study

In this work, the experimental and theoretical vibrational spectra of 2-chloro-4-methylaniline (2Cl4MA, C₇H₈NCl) were studied. FT-IR and FT-Raman spectra of 2Cl4MA in the liquid phase have been recorded in the region 4000–400 cm⁻¹ and 3500–50 cm⁻¹, respectively. The structural and spectroscopic data of the molecule in the ground state have been calculated by using Hartree-Fock (HF) and density functional method (B3LYP) with the 6-31G(d), 6-31G(d,p), 6-31+G(d,p), 6-31++G(d,p) and 6-311G(d), 6-311G(d,p), 6-311+G(d,p), 6-311++G(d,p) basis sets. The vibrational frequencies have been calculated and scaled values have been compared with experimental FT-IR and FT-Raman spectra. The observed and calculated frequencies are found to be in good agreement. The complete assignments were performed on the basis of the total energy distribution (TED) of the vibrational modes, calculated with scaled quantum mechanics (SQM) method. The DFT-B3LYP/6-311++G(d,p) calculations have been found more reliable than the ab initio HF/6-311++G(d,p) calculations for the vibrational study of 2Cl4MA. The optimized geometric parameters (bond lengths and bond angles) were compared with experimental values of aniline and p-methylaniline molecules.     

FT-IR, FT-Raman spectra and quantum chemical calculations of some chloro substituted phenoxy acetic acids

In this work, we will report a combined experimental and theoretical study on molecular and vibrational structure of o-chlorophenoxy acetic acid (OCPAA) and p-chlorophenoxy acetic acids (PCPAA). The FT-IR and Fourier transform-Raman spectra of both the compounds was recorded in the solid phase. The optimized geometry was calculated by HF and B3LYP methods with 6-311++G(d,p) basis set and harmonic vibrational frequencies, infrared intensities and Raman scattering activities were calculated by density functional B3LYP method with the 6-311++G(d,p) basis set. The scaled theoretical wavenumbers showed very good agreement with the experimental values. The thermodynamic functions of the title compounds were also performed at B3LYP/6-311++G(d,p) level of theory. A detailed interpretation of the infrared and Raman spectra of o-chloro and p-chlorophenoxy acetic acid is reported. The theoretical FT-IR spectrograms for the title molecules have been constructed.     

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 studies on dynamics and thermochemistry of the reactions CF(3)CHCl(2)+Cl-->CF(3)CCl(2)+HCl and CF(3)CHFCl+Cl-->CF(3)CFCl+HCl

A dual-level direct dynamics study has been carried out for the two hydrogen abstraction reactions CF(3)CHCl(2)+Cl and CF(3)CHFCl+Cl. The geometries and frequencies of the stationary points are optimized at the BHLYP/6-311G(d,p), B3LYP/6-311G(d,p), and MP2/6-31G(d) levels, respectively, with single-point calculations for energy at the BHLYP/6-311++G(3df,2p), G3(MP2), and QCISD(T)/6-311G(d,p) levels. The enthalpies of formation for the species CF(3)CHCl(2), CF(3)CHFCl, CF(3)CCl(2), and CF(3)CFCl are evaluated at higher levels. With the information of the potential energy surface at BHLYP/6-311++G(3df,2p)//6-311G(d,p) level, we employ canonical variational transition-state theory with small-curvature tunneling correction to calculate the rate constants. The agreement between theoretical and experimental rate constants is good in the measured temperature range 276-382 K. The effect of fluorine substitution on reactivity of the C-H bond is discussed.     

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.     

Comparison of some representative density functional theory and wave function theory methods for the studies of amino acids

Energies of different conformers of 22 amino acid molecules and their protonated and deprotonated species were calculated by some density functional theory (DFT; SVWN, B3LYP, B3PW91, MPWB1K, BHandHLYP) and wave function theory (WFT; HF, MP2) methods with the 6-311++G(d,p) basis set to obtain the relative conformer energies, vertical electron detachment energies, deprotonation energies, and proton affinities. Taking the CCSD/6-311++G(d,p) results as the references, the performances of the tested DFT and WFT methods for amino acids with various intramolecular hydrogen bonds were determined. The BHandHLYP method was the best overall performer among the tested DFT methods, and its accuracy was even better than that of the more expensive MP2 method. The computational dependencies of the five DFT methods and the HF and MP2 methods on the basis sets were further examined with the 6-31G(d,p), 6-311++G(d,p), aug-cc-pVDZ, 6-311++G(2df,p), and aug-cc-pVTZ basis sets. The differences between the small and large basis set results have decreased quickly for the hybrid generalized gradient approximation (GGA) methods. The basis set convergence of the MP2 results has been, however, very slow. Considering both the cost and the accuracy, the BHandHLYP functional with the 6-311++G(d,p) basis set is the best choice for the amino acid systems that are rich in hydrogen bonds.     

Thermochemical properties of molecules and cations of MgNnHm (n=1,2 and m=1–6): enthalpies of formation, proton affinities and ionization energies

Ab initio molecular orbital calculations are reported for small neutral molecules and cations containing magnesium, nitrogen and hydrogen. Structures have been optimized using gradient techniques at B3LYP/6-31+G(d) and at MP2(full)/6-311++G(d,p). Single-point calculations are reported at QCISD(T)(full)/6-311++G(2df,p) and at CCSD(T)(full)/6-311++G(2df,p) levels using geometries optimized at MP2(full)/6-311++G(d,p). Standard enthalpies of formation at 298 K have been calculated at these two higher levels of theory. Other thermochemical properties calculated include ionization energies and proton affinities. The binding enthalpies of ammonia to Mg⁺, MgNH₂⁺ and MgNH₃⁺ are also reported.     

Theoretical Studies on the Aminolysis of Ester: Effects of the Catalysis and Substituent to the Reaction of Methyl Indole-3-acetate with Ammonia

A comprehensive exploration of the aminolysis mechanism for methyl indole-3-acetate with ammonia is carried out by employing the B3 LYP/6-311++G(d,p), M06-2 X/6-311++G(d,p) and MP2/6-311++G(d,p)//M06-2 X/6-311++G(d,p) levels. Two alterative reaction channels of the concerted and addition/elimination stepwise processes including the uncatalyzed, base-catalyzed reactions are taken into consideration. Subsequently, the substituent effects and solvent effects in methanol are also evaluated at the M06-2 X/6-311++G(d,p) level. The calculated results indicate that the calculated values of M06-2 X level are quite close to those of MP2, the stepwise pathway has more advantages to the concerted one for all of the reaction processes and the catalyst facilitates the proton migration and decreases the energy barriers as well. It is shown that the most preferred mechanism is the based-catalyzed stepwise process, the substituent of NH2 group slightly accelerates all the aminolysis reaction processes, and the solvent effect does not remarkably change the mechanism of the reaction.     

Evaluation of the enthalpy of formation, proton affinity, and gas-phase basicity of gamma-butyrolactone and 2-pyrrolidinone by isodesmic reactions

The knowledge of thermochemical parameters such as the enthalpy of formation, gas-phase basicity, and proton affinity may be the key to understanding molecular reactivity. The obtention of these thermochemical parameters by theoretical chemical models may be advantageous when experimental measurements are difficult to accomplish. The development of ab initio composite models represents a major advance in the obtention of these thermochemical parameters, but these methods do not always lead to accurate values. Aiming at achieving a comparison between the ab initio models and the hybrid models based on the density functional theory (DFT), we have studied gamma-butyrolactone and 2-pyrrolidinone with a goal of obtaining high-quality thermochemical parameters using the composite chemical models G2, G2MP2, MP2, G3, CBS-Q, CBS-4, and CBS-QB3; the DFT methods B3LYP, B3P86, PW91PW91, mPW1PW, and B98; and the basis sets 6-31G(d), 6-31+G(d), 6-31G(d,p), 6-31+G(d,p), 6-31++G(d,p), 6-311G(d), 6-311+G(d), 6-311G(d,p), 6-311+G(d,p), 6-311++G(d,p), aug-cc-pVDZ, and aug-cc-pVTZ. Values obtained for the enthalpies of formation, proton affinity, and gas-phase basicity of the two target molecules were compared to the experimental data reported in the literature. The best results were achieved with the use of DFT models, and the B3LYP method led to the most accurate data.     

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.     

Peptide models. XXXIII. Extrapolation of low-level Hartree-Fock data of peptide conformation to large basis set SCF,MP2, DFT,and CCSD(T) results. The Ramachandran surface of alanine dipeptide computed at various levels of theory

At the dawn of the new millenium, new concepts are required for a more profound understanding of protein structures. Together with NMR and X-ray-based 3D-structure determinations in silico methods are now widely accepted. Homology-based modeling studies, molecular dynamics methods, and quantum mechanical approaches are more commonly used. Despite the steady and exponential increase in computational power, high level ab initio methods will not be in common use for studying the structure and dynamics of large peptides and proteins in the near future. We are presenting here a novel approach, in which low- and medium-level ab initio energy results are scaled, thus extrapolating to a higher level of information. This scaling is of special significance, because we observed previously on molecular properties such as energy, chemical shielding data, etc., determined at a higher theoretical level, do correlate better with experimental data, than those originating from lower theoretical treatments. The Ramachandran surface of an alanine dipeptide now determined at six different levels of theory [RHF and B3LYP 3-21G, 6-31+G(d) and 6-311++G(d,p)] serves as a suitable test. Minima, first-order critical points and partially optimized structures, determined at different levels of theory (SCF, DFT), were completed with high level energy calculations such as MP2, MP4D, and CCSD(T). For the first time three different CCSD(T) sets of energies were determined for all stable B3LYP/6-311++G(d,p) minima of an alanine dipeptide. From the simplest ab initio data (e.g., RHF/3-21G) to more complex results [CCSD(T)/6-311+G(d,p)//B3LYP/6-311++G(d,p)] all data sets were compared, analyzed in a comprehensive manner, and evaluated by means of statistics.     

The azido-tetrazole and diazo-1,2,3-triazole tautomerism in six-membered heteroaromatic rings and their relationships with aromaticity: Azines and perimidine

The properties of 28 molecules together with 12 transition states belonging to the series of azido-azines, tetrazolo-azines, diazo-azines, and 1,2,3-triazolo-azines have been studied at the B3LYP/6-31G(d), B3LYP/6-311++G(d,p) and, for 26 cases, at the G3B3 level. Energies, NICS and bond critical points were used to discuss the ring-chain tautomerism of these compounds in relation with the aromaticity of the azines (pyridine, pyrimidine, quinazoline, 1,3-diazapyrene, and perimidine) and the azoles (tetrazoles and 1,2,3-triazoles).