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.     

Study of conformational and optical rotation for the alaninamide

Six stationary points of alaninamide have been located on the potential surface energy (PES) at the B3LYP/6‐311++G(2d,2p) level of theory both in the gas phase and in aqueous solution. In the aqueous solution, to take the water solvent effect into account, the polarizable continuum model (PCM) method has been used. Accurate geometric structures and their relative stabilities have been investigated. The results show that the intramolecular hydrogen bond plays a very important role in stabilizing the global minimum of the alaninamide. Moreover, the consistent result in relative energy using high‐level computations, including the MP2 and MP3 methods with the same basis set [6‐311++G(2d,2p)], indicates that the B3LYP/6‐311++G(d,p) level may be applied to the analogue system. More importantly, the optical rotation of the optimized conformers (both in the gas phase and in aqueous solution) of alaninamide have been calculated using the density functional theory (DFT) and Hartree–Fock (HF) method at various basis sets (6‐31+G*, 6‐311++G(d,p), 6‐311++G(2d,2p) and aug‐cc‐pvdz). The results show that the selection of the computation method and the basis set in calculation has great influence on the results of the optical rotations. The reliability of the HF method is less than that of DFT, and selecting the basis set of 6‐311++G(2d,2p) and aug‐cc‐pvDZ produces relative reliable results. Analysis of the computational results of the structure parameters and the optical rotations yields the conclusion that just the helixes in molecules caused the chiral molecules to be optical active. The Boltzmann equilibrium distributions for the six conformers (both in the gas phase and in the aqueous solution) are also carried out. © 2006 Wiley Periodicals, Inc. Int J Quantum Chem, 2007     

The Role of Intramolecular Hydrogen Bonds vs. Other Weak Interactions on the Conformation of Hyponitrous Acid and Its Mono- and Dithio-Derivatives

High level ab initio and density functional theory calculations have been carried out to investigate the relative stability of the different conformers of hyponitrous acid and its mono- and dithio-derivatives. Geometries and vibrational frequencies were obtained at the B3LYP/6-311+G(d,p) level and final energies through B3LYP/6-311++G(3df,2pd) single point calculations. The reliability of this theoretical scheme has been assessed by comparing these DFT results with those obtained at the G3 level of theory, for some suitable cases. The cis conformers of hyponitrous acid and its mono- and dithio-derivatives are systematically more stable than the trans ones because in the cis conformation a dative interaction between the nitrogen-lone pairs and the σ_NX^* antibonding orbital is significantly favored. Quite interestingly, in general, the conformers presenting an intramolecular hydrogen bond (IHB) are not the global minima of the corresponding potential energy surfaces and only for hyponitrous acid the conformer with a OH ⋅s O IHB is slightly more stable than the cis conformer without IHB. The low stability of the tautomers with IHB is closely related with another weak intramolecular interaction which involves the lone-pairs of the chalcogen atoms and the π_NN^* antibondig orbital, and which is significantly perturbed when the IHB is formed.     

Conformational and tautomeric preferences in 3‐aminoacrylaldehyde: A theoretical study

The conformational study of 3‐aminoacrylaldehyde were performed at various theoretical levels and the equilibrium conformations were determined. Furthermore, to have more reliable energies, the total energies of all forms recomputed at G2MP2 and CBS‐QB3. Theoretical calculations clearly show that the intramolecular hydrogen bond (IHB) is the origin of conformational preference and the resultant IHB order at HF/6‐311++G(d,p), MP2/6‐311++G(d,p), and CBS‐QB3 levels is different from the order which obtained from the B3LYP/6‐311++G(d,p), G2MP2, the geometrical parameters, AIM, and NBO analyses. Furthermore, our theoretical results reveal that the ketoamine (KA) tautomeric group is more stable than the enolimine (EI) and ketoimine (KI) ones. The IHB and tautomeric process could not rationalize the irregular stability of KA group with respect the others. But the population analyses of the possible conformations by NBO predict that the π‐electron delocalization, especially unusual π → π charge transfer, is the origin of tautomeric preference. © 2009 Wiley Periodicals, Inc. Int J Quantum Chem, 2010     

Blue‐shifted and red‐shifted hydrogen bonds: Theoretical study of the CH3CHO· · ·HNO complexes

The blue‐shifted and red‐shifted H‐bonds have been studied in complexes CH₃CHO…HNO. At the MP2/6‐31G(d), MP2/6‐31+G(d,p) MP2/6‐311++G(d,p), B3LYP/6‐31G(d), B3LYP/6‐31+G(d,p) and B3LYP/6‐311++G(d,p) levels, the geometric structures and vibrational frequencies of complexes CH₃CHO…HNO are calculated by both standard and CP‐corrected methods, respectively. Complex A exhibits simultaneously red‐shifted C? H…O and blue‐shifted N? H…O H‐bonds. Complex B possesses simultaneously two blue‐shifted H‐bonds: C? H…O and N? H…O. From NBO analysis, it becomes evident that the red‐shifted C? H…O H‐bond can be explained on the basis of the two opposite effects: hyperconjugation and rehybridization. The blue‐shifted C? H…O H‐bond is a result of conjunct C? H bond strengthening effects of the hyperconjugation and the rehybridization due to existence of the significant electron density redistribution effect. For the blue‐shifted N? H…O H‐bonds, the hyperconjugation is inhibited due to existence of the electron density redistribution effect. The large blue shift of the N? H stretching frequency is observed because the rehybridization dominates the hyperconjugation. © 2006 Wiley Periodicals, Inc. Int J Quantum Chem, 2006     

DFT studies on the multi‐channel reaction of CH3S+NO2

The mechanisms for the reaction of CH₃S with NO₂ are investigated at the QCISD(T)/6‐311++G(d,p)//B3LYP/6‐311++G(d,p) on both single and triple potential energy surfaces (PESs). The geometries, vibrational frequencies, and zero‐point energy (ZPE) correction of all stationary points involved in the title reaction are calculated at the B3LYP/6‐311++G(d,p) level. More accurate energies are obtained at the QCISD(T)/6‐311++G(d,p). The results show that 5 intermediates and 14 transition states are found. The reaction is more predominant on the single PES, while it is negligible on the triple PES. Without any barrier height for the whole process, the main channel of the reaction is to form CH₃SONO and then dissociate to CH₃SO+NO. © 2006 Wiley Periodicals, Inc. Int J Quantum Chem, 2007     

Theoretical calculations on the mechanisms of the gas‐phase elimination kinetics of 2‐chloro‐1‐phenylethane, 3‐chloro‐1‐phenylpropane, 4‐chloro‐1‐phenylbutane, 5‐chloro‐1‐phenylpentane,and their corresponding chloroalkanes: The effect of the phenyl ring

The kinetics and mechanisms of the dehydrochlorination of 2‐chloro‐1‐ phenylethane, 3‐chloro‐1‐phenylpropane, 4‐chloro‐1‐phenylbutane, 5‐chloro‐1‐phenylpentane, and their corresponding chloroalkanes were examined by means of electronic structure calculation using density functional theory methods B3LYP/6–31G(d,p), B3LYP/6–31++G(d,p), MPW1PW91/6–31G(d,p), MPW1PW91/6–31++G(d,p), PBEPBE/6–31G(d,p), and PBEPBE/6–31++G(d,p). The potential energy surface was investigated for the minimum energy path. Calculated enthalpies and energies of activation are in good agreement with experimental values using the MPW1PW91 and B3LYP methods. The transition state of these reactions is a four‐centered cyclic structure. The reported experimental results proposing neighboring group participation by the phenyl group was not supported by theoretical calculations. The rate‐determining process in these reactions is the breaking of Cl? C bond. The reactions are described as concerted moderately polar and nonsynchronous. © 2011 Wiley Peiodicals, Inc. Int J Chem Kinet 43: 292–302, 2011     

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.     

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.     

Substituent effect on intramolecular hydrogen bonding in 2‐hydroxybenzaldehyde

The effect of some substituents on intramolecular hydrogen bonding of 5‐X‐2‐hydroxybenzaldehyde (5‐X‐2‐HBA) has been studied by B3LYP and MP2 methods using 6‐311++G** and AUG‐cc‐PVTZ basis sets. The relationship between hydrogen bond energy E_HB and electron donation (or withdrawal) of substituents has been investigated. An approximately good linear relationship has been detected between Hammett coefficients and hydrogen bond formation energy (R² = 0.98). Herein, population analysis has been performed by atoms in molecules (AIM) and natural bond orbital (NBO) methods. The results of AIM and NBO analyses are in a good agreement with calculated energy values. Furthermore, correlation between ring aromaticity and hydrogen bonding has been investigated by nucleus‐independent chemical shift (NICS) at GIAO/B3LYP/6‐311++G** level of theory. © 2009 Wiley Periodicals, Inc. Int J Quantum Chem, 2009     

复合物HNO···H2O2内N-H···O蓝移氢键的理论研究

A theoretical study on the blue-shifted H-bond N-H…O and red-shifted H-bond O-H…O in the complexHNO…H_2O_2 was conducted by employment of both standard and counterpoise-corrected methods to calculate thegeometric structures and vibrational frequencies at the MP2/6-31G(d),MP2/6-31 G(d,p),MP2/6-311 q G(d,p),B3LYP/6-31G(d),B3LYP/6-31 G(d,p) and B3LYP/6-311 G(d,p) levels.In the H-bond N-H…O,the calcu-lated blue shift of N-H stretching frequency is in the vicinity of 120 cm~(-1) and this is indeed the largest theoreticalestimate of a blue shift in the X-H…Y H-bond ever reported in the literature.From the natural bond orbital analy-sis,the red-shifted H-bond O-H…O can be explained on the basis of the dominant role of the hyperconjugation.For the blue-shifted H-bond N-H…O,the hyperconjugation was inhibited due to the existence of significant elec-tron density redistribution effect,and the large blue shift of the N-H stretching frequency was prominently due tothe rehybridization of sp~n N-H hybrid orbital.     

Optical rotation and two chiral C atoms of aliphatic cyclic dipeptides

Accurate geometries structures and total energies have been determined for the conformers of cyclo(L-Pro-Gly), cyclo(L-Ala-L-Ala), and cyclo(L-Pro-Ala) in the gaseous phase, using HF and B3LYP correlation methods at 6−31++G(d), 6−311++G(d, p), 6−311++G(2d, 2p) and aug-cc-pvdz basis sets. High level computations MP2 with 6−311++G(2d, 2p) basis set indicate that the relative stabilities of the available conformers can be determined correctly at the B3LYP/6−311++g(2d, 2p) level of theory. We have also described the implementation of DFT and HF theory for calculations of the optical rotation at 589.3 nm. In L-Ala-L-Ala, and L-Pro-Ala molecules, they have two chiral C (C*), so we discuss the different effect of two chiral C to optical activity of cydo(L-Pro-Gly), cyclo(L-Ala-L-Ala), and cyclo(L-Pro-Ala).     

Intramolecular hydrogen bond in 3‐imino‐propenylamine isomers: AIM and NBO studies

The molecular structure and intramolecular hydrogen bond energy of 18 conformers of 3‐imino‐propenyl‐amine were investigated at MP2 and B3LYP levels of theory using the standard 6‐311++G** basis set. The atom in molecules or AIM theory of Bader, which is based on the topological properties of the electron density (ρ), was used additionally and the natural bond orbital (NBO) analysis was also carried out. Furthermore calculations for all possible conformations of 3‐imino‐propenyl‐amin in water solution were also carried out at B3LYP/6‐311++G** and MP2/6‐311++G** levels of theory. The calculated geometrical parameters and conformational analyses in gas phase and water solution show that the imine–amine conformers of this compound are more stable than the other conformers. B3LYP method predicts the IMA‐1 as global minimum. This stability is mainly due to the formation of a strong N? H···N intramolecular hydrogen bond, which is assisted by π‐electrons resonance, and this π‐electrons are established by NH2 functional group. Hydrogen bond energies for all conformers of 3‐imino‐propenyl‐amine were obtained from the related rotamers methods. © 2009 Wiley Periodicals, Inc. Int J Quantum Chem, 2010     

Factors influencing the C?ON bond strength of the alkoxyamines in the styrene–acrylonitrile–TEMPO copolymerization system

Homolytic bond dissociation energy (BDE) of the (C? ON) bond for several N‐alkoxyamines derived from 2,2,6,6‐tetramethylpiperidine‐1‐oxyl (TEMPO) and the corresponding (C? H) bonds were determined from quantum‐mechanical calculations including the B3‐LYP6‐31G(d), B3‐LYP/6‐311++G(2df,p), UB3‐LYP/6‐311+G(3df,2p), and integrated IMOMO (G3:ROMP2/6‐31G(d)) method. The investigated N‐alkoxyamines were considered as models for dormant forms of propagating chains in the radical copolymerization process of styrene with acrylonitrile in the presence of TEMPO according to the terminal and penultimate model. The substituent effect on BDE was investigated. Radical stabilization energies (RSE) for radicals created from homolysis of the investigated N‐alkoxyamines were calculated according to Rüchardt's method. Polar, steric, and stabilization effects on C? ON alkoxyamine bond homolysis were studied. A dramatically weakened C? ON bond in the alkoxyamine‐containing two consecutive styrene units in the propagating chain was ascribed to geometric parameters characterizing energetically unfavorable conformation of the substituents. These phenomena can be regarded as the penultimate effect in the radical living/controlled copolymerization system of styrene with acrylonitrile. © 2008 Wiley Periodicals, Inc. J Polym Sci Part A: Polym Chem 46: 1165–1177, 2008     

Calculation of dissociation constants and chemical hardness of some biologically important molecules: A theoretical study

The gas‐phase geometries of neutral, protonated, and deprotonated forms of some biologically important molecules, alanine (Ala), glycine (Gly), phenylalanine (Phe), and tyrosine (Tyr), were optimized using density functional theory at B3LYP/6‐311++G(d) and the ab initio HF/6‐311++G(d) level of theories. The neutral and different stable ionic states of Ala, Gly, Phe, and Tyr have also been solvated in aqueous medium using polarizable continuum model for the determination of solvation free energies in the aqueous solution. The gas‐phase acidity constants of above four molecules have been also calculated at both levels of theories and found that the values calculated at HF/6‐311++G(d) method are in good agreement with experimental results. A thermodynamic cycle was used to determine the solvation free energies for the proton dissociation process in aqueous solution and the corresponding pK_a values of these molecules. The pK_a values calculated at B3LYP/6‐311++G(d) method are well supported by the experimental data with a mean absolute deviation 0.12 pK_a units. Additionally, the chemical hardness and the ionization potential (IP) for these molecules have been also explored at both the level of theories. The Tyr has less value of chemical hardness and IP at both levels of theories compared with other three molecules, Ala, Gly, and Phe. The calculated values of chemical hardness and IP are decreasing gradually with the substitution of the various functional groups in the side chain of the amino acids. © 2010 Wiley Periodicals, Inc. Int J Quantum Chem, 2011     

Contemplation on Some Prismatic Polynitrogen Structures – A DFT Treatment

Polynitrogen compounds, are rare molecules having only nitrogen atoms. In recent years, they have been considered as promising candidates of clean (green) high energy density materials. They possess high energy content and their sole decomposition product is N₂. Presently some prismatic polynitrogen structures (N₆ – N₁₄) are considered within the limitations of density functional theory at the levels of B3LYP/6‐311++G(d,p) and B3LYP/cc‐PVTZ. The calculations reveal that they are all highly endothermic but stable. Certain quantum chemical properties, IR and UV/Vis spectra are reported. Homolytic bond cleavage of top rings are considered. Then, the transition state and activation energies, and also homolytic bond dissociation energies for the top rings have been calculated at the level of UB3LYP/6‐311++G(d,p). Also NICS(0) values have been calculated at the level of B3LYP/6‐311++G(d,p). The results indicate that N₆ and N₁₀ have aromatic and nonaromatic top (and also base) rings, respectively. All the rest of the structures have antiaromatic rings (all the structures have antiaromatic N₄ rings).     

Isomerizations of CH3SO2: Quantum Chemistry and Topological Study

The structures, stabilities and the isomerization reactions of CH₃SO₂ isomers in a doublet electronic state have been studied at B3LYP/6‐311+ +G (d,p), MP2/6‐311++G (d,p) and CCSD(T)/6‐311++G (d,p) levels. The three different levels of calculation give the similar results: thirteen minimum isomers were located and they were connected by eleven transition states. Among the thirteen isomers, cis‐CH₃OSO, trans‐CH₃OSO and CH₃SO₂ are the most stable species, and they should be detected easily in experiment. This is well consistent with the experimental result. These isomers could isomerize to each other by chemical bond vibration, chemical bond rotation and atom migration. The non‐planar ring structure transition state (STS), which was found in this paper, extended the concept of ring STS to the non‐planar systems.     

Rate constants calculation of hydrogen abstraction reactions CH3CHBr + HBr and CH3CBr2 + HBr

Dual‐level direct dynamics method is used to study the kinetic properties of the hydrogen abstraction reactions of CH₃CHBr + HBr → CH₃CH₂Br + Br (R1) and CH₃CBr₂ + HBr → CH₃CHBr₂ + Br (R2). Optimized geometries and frequencies of all the stationary points and extra points along the minimum‐energy path are obtained at the MPW1K/6‐311+G(d,p), MPW1K/ma‐TZVP, and BMK/6‐311+G(d,p) levels. Two complexes with energies less than that of the reactants are located in the entrance of each reaction at the MPW1K/6‐311+G(d,p) and MPW1K/ma‐TZVP levels, respectively. The energy profiles are further refined with the interpolated single‐point energies method at the G2M(RCC5)//MPW1K/6‐311+G(d,p) level of theory. By the improved canonical variational transition‐state theory with the small‐curvature tunneling correction (SCT), the rate constants are evaluated over a wide temperature range of 200–2000 K. Our calculations have shown that the radical reactivity decreases from CH₃CHBr to CH₃CBr₂. Finally, the total rate constants are fitted by two modified Arrhenius expression. © 2012 Wiley Periodicals, Inc.     

Adsorption of benzyl alcohol and 1-phenylethanol on graphitized thermal carbon black

The adsorption of benzyl alcohol and 1-phenylethanol on graphitized thermal carbon black at 300 K is investigated. Thermodynamic characteristics of the adsorption of conformers stable in the gas phase and in the adsorbed state are determined by a molecular statistics method. The geometry of the molecules is optimized using the MP2 and B3LYP methods with a 6-311++G(d,p) basis set. Using the PBE0/6-31G(d) method, it is established that, in the adsorption of benzyl alcohol, two types of interactions, CH…O and OH…π_GTCB, occur. It is shown that, in 1-phenylethanol, intramolecular interactions are maintained.     

Intramolecular hydrogen bond in the hydroxycyclohexadienyl peroxy radicals

The hydroxycyclohexadienyl peroxy radicals (HO? C₆H₆? O₂) produced from the reaction of OH‐benzene adduct with O₂ were studied with density functional theory (DFT) calculations to determine their characteristics. The optimized geometries, vibrational frequencies, and total energies of 2‐hydroxycyclohexadienyl peroxy radical II_s and 4‐hydroxycyclohexadienyl peroxy radical III_s were calculated at the following theoretical levels, B3LYP/6‐31G(d), B3LYP/6‐311G(d,p), and B3LYP/6‐311+G(d,p). Both were shown to contain a red‐shifted intramolecular hydrogen bond (O? H … O? H bond). According to atoms‐in‐molecules (AIM) analysis, the intramolecular hydrogen bond in the 2‐hydroxycyclohexadienyl peroxy radical II_s is stronger than that one in 4‐hydroxycyclohexadienyl peroxy radical III_s, and the former is the most stable conformation among its isomers. Generally speaking, hydrogen bonding in these radicals plays an important role to make them more stable. Based on natural bond orbital (NBO) analysis, the stabilization energy between orbitals is the main factor to produce red‐shifted intramolecular hydrogen bond within these peroxy radicals. The hyperconjugative interactions can promote the transfer of some electron density to the O? H antibonding orbital, while the increased electron density in the O? H antibonding orbital leads to the elongation of the O? H bond and the red shift of the O? H stretching frequency. © 2006 Wiley Periodicals, Inc. Int J Quantum Chem, 2007