Tautomerism of protonated imidazoles: A perspective from ab initio valence bond theory

The tautomerism of protonated imidazoles concerns one aromatic (1H-Imi+) and one nonaromatic (4H-Imi+) tautomers. Both experiments and computations have shown that substituents to the imidazole ring can change the relative stability of tautomers. A detailed theoretical study of the inherent mechanism would benefit the rational design of experimental syntheses related to imidazoles. In this work, we used the block-localized wavefunction (BLW) method to explore the factors governing the tautomerism between 1H-Imi+ and 4H-Imi+. While π resonance always favors the aromatic tautomer 1H-Imi+, the aromaticity noticeably reduces with electron donating groups (EDGs) as substituents due to the increased π-π repulsion, leading to the stability of 4H-Imi+ over 1H-Imi+ with EDGs such as NH₂ and OH. For electron withdrawing groups (EWGs), the reduced π-π repulsion promotes the aromatic stability and favors 1H-Imi+. DFT computations were also performed to study the tautomerism mechanisms. Results show that tautomerism can hardly occur in gaseous phase, but in aqueous solution, water molecules can build hydrogen bonding network with 1H-Imi+ and facilitate the hydrogen transfers.     

XMVB: a program for ab initio nonorthogonal valence bond computations

An ab initio nonorthogonal valence bond program, called XMVB, is described in this article. The XMVB package uses Heitler-London-Slater-Pauling (HLSP) functions as state functions, and calculations can be performed with either all independent state functions for a molecule or preferably a few selected important state functions. Both our proposed paired-permanent-determinant approach and conventional Slater determinant expansion algorithm are implemented for the evaluation of the Hamiltonian and overlap matrix elements among VB functions. XMVB contains the capabilities of valence bond self-consistent field (VBSCF), breathing orbital valence bond (BOVB), and valence bond configuration interaction (VBCI) computations. The VB orbitals, used to construct VB functions, can be defined flexibly in the calculations depending on particular applications and focused problems, and they may be strictly localized, delocalized, or bonded-distorted (semidelocalized). The parallel version of XMVB based on MPI (Message Passing Interface) is also available.     

Three new phosphoric triamides with a [C(O)NH]P(O)[N(C)(C)]2 skeleton: a database analysis of C—N—C and P—N—C bond angles

In N,N,N′,N′‐tetraethyl‐N′′‐(4‐fluorobenzoyl)phosphoric triamide, C₁₅H₂₅FN₃O₂P, (I), and N‐(2,6‐difluorobenzoyl)‐N′,N′′‐bis(4‐methylpiperidin‐1‐yl)phosphoric triamide, C₁₉H₂₈F₂N₃O₂P, (II), the C—N—C angle at each tertiary N atom is significantly smaller than the two P—N—C angles. For the other new structure, N,N′‐dicyclohexyl‐N′′‐(2‐fluorobenzoyl)‐N,N′‐dimethylphosphoric triamide, C₂₁H₃₃FN₃O₂P, (III), one C—N—C angle [117.08 (12)°] has a greater value than the related P—N—C angle [115.59 (9)°] at the same N atom. Furthermore, for most of the analogous structures with a [C(=O)NH]P(=O)[N(C)(C)]₂ skeleton deposited in the Cambridge Structural Database [CSD; Allen (2002). Acta Cryst. B 58 , 380–388], the C—N—C angle is significantly smaller than the two P—N—C angles; exceptions were found for four structures with the N‐methylcyclohexylamide substituent, similar to (III), one structure with the seven‐membered cyclic amide azepan‐1‐yl substituent and one structure with an N‐methylbenzylamide substituent. The asymmetric units of (I), (II) and (III) contain one molecule, and in the crystal structures, adjacent molecules are linked via pairs of N—H...O=P hydrogen bonds to form dimers.     

Agostic Interactions in Early Transition-Metal Complexes: Roles of Hyperconjugation,Dispersion, and Steric Effect

The agostic interaction is a ubiquitous phenomenon in catalytic processes and transition-metal complexes, and hyperconjugation has been well recognized as its origin. Yet, recent studies showed that either short-range London dispersion or structural constraints could be the driving force, although proper evaluation of the role of hyperconjugation therein is needed. Herein, a simple variant of valence bond theory was employed to study a few exemplary Ti complexes with α- or β-agostic interactions and interpret the agostic effect in terms of the steric effect, hyperconjugation, and dispersion. For the complexes [MeTiCl₃(dmpe)] and [MeTiCl₃(dhpe)] with α-agostic interactions, hyperconjugation plays the dominant role with comparable magnitudes in both systems, but dispersion is solely responsible for the stronger agostic interaction in the former compared with the latter. For the complexes [EtTiCl₃(dmpe)] and [EtTiCl₃(dhpe)] with β-agostic interactions, however, hyperconjugation and dispersion play comparable roles, and the weaker steric repulsion leads to a stronger agostic effect in the former than in the latter. Thus, the present study clarifies the variable and sensitive roles of steric, hyperconjugative, and dispersion interactions in the agostic interaction.     

How to properly compute the resonance energy within the ab initio valence bond theory: a response to the ZHJVL paper

Valence bond (VB) theory describes a conjugated system by a set of electron-localized Lewis resonance structures. VB assumes that the magnitude of the intramolecular electron delocalization can be measured in terms of a resonance energy (RE), taken to be the energy difference between the real conjugated system (delocalized) and the corresponding most stable virtual resonance structure (localized). Proper RE estimates within VB theory require both delocalized and localized states to be defined at the same theoretical level, and the definition of the localized state to closely correspond to the intuitive picture of the corresponding VB structure. In contrast, the VB-delocal and VB-local computational approaches adopted by Zielinski, et al. [preceding paper in this issue] used definitions for either the delocalized or the localized states which, in our view, depart from the intuitive chemical picture. Consequently, their RE estimates are much lower than seemingly appropriate experimental evaluations with which they strongly disagree. Very large basis sets approaching completeness blur the boundaries among resonance structures and result in “basis set artifact” problems within any variant of VB theory. However, block-localized wavefunction (BLW) computations with mid-size basis sets not only exhibit insignificant variations with theoretical levels, but the resulting RE estimates also are justified by comparisons with those employing experimental data and MO computations. We stress that RE differs from the aromatic stabilization energy (ASE). The RE measures the total stabilization of an aromatic system, whereas ASE measures only the part of the RE that exceeds that of appropriate conjugated (but non-aromatic) reference molecules.     

Origin of Stability in Branched Alkanes

The potential origins of stability in branched alkanes are investigated, paying close attention to two recent hypotheses: geminal steric repulsion and protobranching. All alkane isomers through C₆H₁₄ along with heptane and octane were investigated at the MPW1B95/6‐311++G(d,p) level. Their geminal steric repulsion, total steric repulsion, and orbital interactions were evaluated by using natural bond orbital analysis. All measures of steric repulsion fail to explain the stability of branched alkanes. The extra stability of branched alkanes and protobranching, in general, is tied to stabilizing geminal σ→σ* delocalization, particularly of the type that involves adjacent C? C bonds and, thus, preferentially stabilizes branched alkanes. This picture is corroborated by valence bond calculations that attribute the effect to additional ionic structures (e.g., CH₃⁺ :CH₂ :CH₃^? and CH₃:^? CH₂: CH₃⁺ for propane) that are not possible without protobranching.     

重叠式乙烷构象是过渡态吗*

从化学和数学的角度指出乙烷的重叠式构象是过渡态而不是稳定构象,剖析分子旋转异构过程中的构象变化过渡态的本质。同时,采用Python语言编程实现了过渡态搜索算法并开展了乙烷旋转构象变化过渡态的自动搜索,成功得到旋转过渡态-重叠式乙烷结构,加深了对重叠式乙烷构象是乙烷旋转构象变化过程过渡态的理解和认识。     

Birch reduction of hexaphenyl- and pentaphenylbenzene and an X-ray crystallography and NMR spectroscopy study of cis- and epi-1,2,3,4,5,6-hexaphenylcyclohexane and of 2,3,5,6-tetraphenyl-1,1'-bicyclohexylidene: Cannizzaro's conundrum revisited

The Birch reduction of hexaphenylbenzene yields two isomers of 1,2,3,4,5,6-hexaphenylcyclohexane. The X-ray crystal structure of the all-cis isomer, 1, reveals that the severe steric crowding among the three axial phenyls is alleviated by a marked splaying out of those three aryl substituents relative to the positioning in a conventional chair structure. A second product, 2, was identified crystallographically and by NMR spectroscopy as the 1,3-diaxial-2,4,5,6-tetraequatorial (epi) isomer of hexaphenylcyclohexane, in which only five of the six additional hydrogen atoms are positioned on the same face of the C(6)Ph(6) precursor. A variable-temperature NMR study of the all-cis isomer 1 yielded a chair-to-chair inversion barrier of approximately 19 kcal mol(-1), which is somewhat higher than the previously reported values for all-cis-1,2,3,4,5,6-C(6)H(6)R(6) in which R=Me or CO(2)Me. The possible relevance to Cannizzaro's 1854 report of a product with the formula (C(7)H(6))(n) is discussed. By contrast, Birch reduction of pentaphenylbenzene led to the formation of 2,3,5,6-tetraphenyl-1,1'-bicyclohexylidene.     

Steric repulsions, rotation barriers, and stereoelectronic effects: a real space perspective

Widely used chemical concepts like Pauli repulsion or hyperconjugation, and their role in determining rotation barriers or stereoelectronic effects, are analyzed from the real space perspective of the interacting quantum atoms approach (IQA). IQA emerges from the quantum theory of atoms in molecules (QTAIM), but is free from the equilibrium geometry constraint of the former. A framework with both electronically unrelaxed and relaxed wavefunctions is presented that leads to an approximate correspondence between the IQA concepts and those used in the EDA (energy decomposition analysis) or NBO (natural bond orbital) procedures. We show that no net force acts upon the electrons in an electronically relaxed system, so that any reasonable definition of Pauli repulsion must involve unrelaxed state functions. Using antisymmetrized fragments clarifies that Pauli repulsions are energetically connected to the IQA deformation energies, leaving footprints in the finally relaxed states. Similarly, EDA or NBO hyperconjugative stabilizations are found to be naturally related to the IQA electron delocalization patterns. Applications to the rotation barrier of ethane and other simple systems are presented, and the very often forgotten role of electrostatic contributions in determining preferred conformations is highlighted.     

Electronic structure of NH+: An ab initio study

We report ab initio self‐consistent field MRSD‐CI electronic structure calculations of the NH⁺ cation. A basis set of DZ + POL quality augmented with Rydberg and bond functions was employed together with an extensive treatment of electron correlation. More than 50 electronic states of NH⁺ are reported, including doublets, quartets, and sextets. Leading configurations, vertical ionization energies of NH, vertical excitation energies of NH⁺, and potential energy curves are reported. Spectroscopic properties calculated for the known bound electronic states of NH⁺ are found in good agreement with experiment. © 2004 Wiley Periodicals, Inc. Int J Quantum Chem, 2005     

Restricted rotation in unbridged sandwich complexes: rotational behavior of closo-[Co(eta 5-NC4H4)(C2B9H11)] derivatives

Rotation about the centroid/metal/centroid axis in ferrocene is facile; the activation energy is 1-5 kcal mol(-1). The structurally similar sandwich complexes derived from closo-[3-Co(eta5-NC4H4)-1,2-C2B9H11] (1) have a different rotational habit. In 1, the cis rotamer in which the pyrrolyl nitrogen atom bisects the carboranyl cluster atoms is 3.5 kcal mol(-1) more stable in energy than the rotamer that is second lowest in energy. This cis rotamer is wide, spanning 216 degrees , and may be split into three rotamers of almost equal energy by substituting the N and the carboranyl carbon atoms adequately. To support this statement, closo-[3-Co(eta5-NC4H4)-1,2-(CH3)2-1,2-C2B9H9] (2), closo-[3-Co(eta5-NC4H4)-1,2-(mu-CH2)3-1,2-C2B9H9] 3, 2-->BF3, and 3-->BF3 have been prepared. Two rotamers are found at low temperature for 2-->BF(3) and 3-->BF3. Compounds 2, 3, and 1-->BF3 behave similarly to 1. Rotational energy barriers and the relative populations of the different energy states are calculated from 1H DNMR spectroscopy (DNMR, dynamic NMR). These results agree with those of semiempirical calculations. Without exception, the cis rotamer is energetically the more stable. The fixed conformation of 1 assists in elucidating the rotational preferences of the [3,3'-Co(1,2-C2B9H11)2]- ion in the absence of steric hindrance; the [3,3'-Co(1,2-C2B9H11)2]- ion is commonly accepted to present a cisoid orientation. Complex 1 is electronically similar to the [3,3'-Co(1,2-C2B9H11)2]- ion. Both have heteroatoms in the pi ligands, and they have the same electronegativity difference between the constituent atoms. This leads to a view of the [NC4H4]- as [7,8-C2B9H11]2- ion, with no steric implications. Therefore the [3,3'-Co(1,2-C2B9H11)2]- ion should be considered to have a cisoid structure, and the different rotamers observed to be the result of steric factors and of the interaction of the counterion with either B-H groups and/or ancillary ligands. The rotamer adopted is the one with the atoms holding the negative charges furthest apart.     

Comparative study of the electronic structure of pregnanolones by ab initio theory

Pregnanedione (5β‐pregnane, 3,20‐dione), pregnanolone (3β‐hydroxy‐5β‐pregnan‐20‐one), and epipregnanolone (3α‐hydroxy‐5β‐pregnan‐20‐one) result from the 5β‐reduction of progesterone [4‐pregnene, 3‐20‐dione (P)]. These P metabolites induce anesthesia and smooth muscle relaxation (nongenomic actions). In the present study, geometries and electronic structure of these steroids were assessed by ab initio calculations using the 6‐31G* basis set. Consequently, bond distances, valence angles, and dihedral angles were measured. In addition total energy, frontier orbitals, i.e., highest occupied molecular orbital (HOMO), lowest unoccupied molecular orbital (LUMO), dipole moment, and electrostatic potentials were calculated. Total energy was higher for P, followed by pregnanedione. Pregnanolones, the hydroxylated progestins, showed the lower energies. Concerning frontier orbitals, P showed the highest HOMO energy and the lowest LUMO energy. Pregnanedione showed lower HOMO and LUMO energy values than pregnanolone and epipregnanolone. P showed both HOMO and LUMO located at the A ring, including the π bond at C4, C5, and the carbonyl at C3. The HOMO in pregnanedione was included mostly in the A ring and the C3 carbonyl group, while the LUMO was shared by the carbonyl groups at C3 and C20. The frontier orbitals of pregnanolone and epipregnanolone were quite similar. The HOMO in both steroids included the B, C, and D rings and the carbonyl at C20. The LUMO was also similar in both pregnanolones including mostly the carbonyl at C20. The dipole moment was shorter for P and pregnanedione and directed toward the acetyl side chain at C17. Pregnanolone and epipregnanolone showed the dipole moment vector larger and directed toward the A ring. The electrostatic potentials were related mostly with the lone pairs of electrons from the oxygens. By the total energy and frontier orbitals energies of the hormones studied, it is concluded that the metabolism of progesterone toward its 5β‐reduced metabolites might be rationalized from the theoretical chemistry point of view. Besides, the importance of the A/B ring cis configuration, dipole moment, and electrostatic potential are highlighted as possible improving elements of molecular interactions to explain the nongenomic biological action of 5β‐reduced progestins. ©1999 John Wiley & Sons, Inc. Int J Quant Chem 71: 433–440, 1999     

The control of the nitrogen inversion in alkyl-substituted diaziridines

For the first time the nitrogen inversion barriers in 3,3-unsubstituted trans-diaziridines, such as 1,2-di-tert-butyldiaziridine (1) and 1,2-di-n-butyldiaziridine (2) were determined. Enantioselective stopped-flow multidimensional gas chromatography was used to investigate the enantiomerization barrier of 1 between 126.2 and 171.0 degrees C (DeltaG ++ gas (150.7 degrees C) = 135.8+/-0.2 kJ mol(-1), DeltaH++ gas = 116.1+/-2.5 kJ mol(-1), DeltaS ++ gas == -46+/-2 J K(-1) mol(-1)). The separation of the enantiomers has been achieved in presence of the chiral stationary phase (CSP) Chirasil-beta-Dex with a high separation factor (alpha = 1.44 at 80 degrees C). In a complementary approach, the enantiomerization barriers of 1,2-di-tert-butyldiaziridine (1), 1,2-di-n-butyldiaziridine (2), 1-n-butyl-3,3-dimethyldiaziridine (3), and 1,2,3,3-tetramethyldiaziridine (4) were determined for comparison by enantioselective dynamic chromatography (DGC) and computer simulation of the dynamic elution profiles. The enantiomerization barrier of 2 was shown to be the highest among the nonsterically hindered diaziridines studied so far, whereas 1 exhibited the highest value found for strained nitrogen-containing rings, that is, aziridines, diaziridines and oxaziridines.     

Structure and conformation of 1-Bromopropane: A gas-phase electron-diffraction investigation using microwave-spectroscopy data and results from ab initio calculations as constraints

1-Bromopropane has been studied by gas-phase electron diffraction (ED) at 24C. Earlier published values for rotational constants from microwave spectroscopy (MW), together with results from ab initio molecular-orbital calculations, have been included in the ED analysis. Two conformers with C-C-C-Br torsion angles of 180 (anti) or 66.0(17) (gauche) have been observed. The results obtained for the bond distances (r _g) and valence angles () from this combined ED/MW analysis, with the ab initio results used as constraints are r (C-H)=1.114(9) å,r(C₁-C₂)=1.521(5) å,r(C₂-C₃)=1.535(5) å,r (C₁-Br)=1.962(6) å, <(C-C-C)anti,=110.0(11), <(C-C-C)gauche=113.3(11), < (C-C-Br)anti=111.1(6), < (C-C-Br)gauche=112.1(6), <C₂-C₁-H=112.1 (ab initio value), <C₂-C₃-H=111.4 (ab initio value), <H-C₂-H=107.0 (ab initio value). Error limits are given as 2, where (standard deviation) includes estimates of uncertainties in voltage/height measurements and correlation in the experimental data. The observed amountof gauche conformer was 64(14)%. Using the entropy difference between conformers obtained in the ab initio calculations, this composition corresponds to an energy difference of E=E _anti —E _gauche=0.03(36) kcal/mol. The results are compared with those earlier obtained for other 1-halopropanes.