Theoretical study on two types of weak interactions between methylenecyclopropane and XY (X,Y = H,F, Cl,and Br)

We present theoretical evidence that the two types of interactions exist in the complexes formed between methylenecyclopropane (MECP) and XY (X, Y = H, F, Cl, and Br). Two seats of XY interacted with MECP are located: (a) is via the pseudo‐π bonding electron pair associated with a C? C bond of the cyclopropane ring and (b) is via the typical‐π bonding of electron pair of the C?C bond of MECP. These two types of weak interactions are compared based on the calculated geometries, interaction energies, frequency changes, and topological properties of electron density. The integration of electron density over the interatomic surface is found to be a good measure for the strength of weak interaction. Furthermore, the total electron density and separated σ and π electron densities are also computed and discussed in this article. The separated electron density shows σ electron density determined the strength and π electron density influenced the direction of the hydrogen/halogen bond. © 2010 Wiley Periodicals, Inc. Int J Quantum Chem, 2011     

Activation of the aromatic system by the SO2CF3 group: Kinetics study and structure–reactivity relationships

The rate constants for the reaction of 2,6‐bis(trifluoromethanesulfonyl)‐4‐nitroanisole with some substituted anilines have been measured by spectrophotometric methods in methanol at various temperatures. The data are consistent with the S_NAr mechanism. The effect of substituents on the rate of reaction has been examined. Good linear relationships were obtained from the plots of log k₁ against Hammett σ_para constants values at all temperature with negative ρ values (?1.68 to ?1.11). Activation parameters Δ^≠H varied from 41.6 to 54.3 kJ mol^?1 and Δ^≠S from ?142.7 to ?114.6 J mol^?1 K^?1. The δΔ^≠H and δΔ^≠S reaction constants were determined from the dependence of Δ^≠H and Δ^≠S activation parameters on the σ substituent constants, by analogy with the Hammett equation. A plot of Δ^≠H versus Δ^≠S for the reaction gave good straight line with 177°C isokinetic temperature. © 2010 Wiley Periodicals, Inc. Int J Chem Kinet 42: 203–210, 2010     

Theoretical investigations on the weak nonbonded CS···CH2 interactions: Chalcogen‐bonded complexes with singlet carbene as an electron donor

In this article, we explored the noncovalent bonding interactions between O?C?S, S?C?S, F₂C?S, Cl₂C?S, and singlet carbene. Six chalcogen‐bonded complexes were obtained. It is found that all the vibrational frequencies of C?S bond presented a red shift character. Interaction energy, topology property of the electron density and its Laplacian, and the donor–acceptor interaction have been investigated. All these results show that there exists a weak nonbonded interaction between the chalcogen bond donor and CH₂. An energy decomposition analysis was performed to disclose that the electrostatic interaction is the main stabilized factor in these nonbonded complexes. © 2010 Wiley Periodicals, Inc. Int J Quantum Chem, 2011     

The effect of the functional,basis set,and solvent in the simulation of the geometry and spectroscopic properties of VIVO2+ complexes. chemical and biological applications

The geometry of 32 V^IVO²⁺ complexes with different donor set, electric charge, geometry, arrangement of the ligands with respect to the V?O bond and type of ligand was calculated by density functional theory methods. 32 V?O, 45 V? O, 16 V? OH, 40 V? N, 24 V? S, and 14 V? Cl bonds were examined. The performance of several functionals (B3LYP, B3P86, B3PW91, HCTH, TPSS, PBE0, and MPW1PW91), keeping constant the Pople triple‐zeta basis sets 6‐311g, was tested. The order of accuracy of the functional in the prediction of the bond distances, expressed in terms of mean of the deviation Δd (Δd = d^calcd ? d^exptl) and absolute deviation |Δd| (|Δd| = |d^calcd ? d^exptl|) from the experimental values and of the corresponding standard deviations (SD(Δd) and SD(|Δd|)), is: B3P86 ～ PBE0 ～ MPW1PW91 > B3PW91 ? TPSS > B3LYP ? HCTH. In the gas phase the prediction of V?O, V? O, V? N bond lengths is rather good, but that of V? OH, V? S and V? Cl distances is by far worse. An improvement in the optimization of V? S and V? Cl lengths is reached by adding polarization and diffuse functions on the sulfur and chlorine atoms. Finally, a general improvement in the prediction of all the calculated bond lengths and angles is obtained by simulating the structures in the solvent where they are isolated within the framework of the polarizable continuum model. The last choice allows also to improve the prediction of structural (the deviation of a penta‐coordinate geometry toward the trigonal bipyramid) and spectroscopic parameters (⁵¹V and ¹⁴N hyperfine coupling constants and ¹⁴N nuclear quadrupolar coupling constant). In most of the cases, the structures optimized in solution closely approach the experimental ones and this can be of great help in the simulations of naturally occurring vanadium compounds and metal site of V‐proteins, like amavadin and the reduced form of vanadium bromoperoxidase (VBrPO). © 2011 Wiley Periodicals, Inc. Int J Quantum Chem, 2011     

Mechanisms of norbornadiene dimerization to Binor‐S using cationic CoI,RhI, and IrI catalysts

We investigated the transition metal‐catalyzed reaction mechanisms of NBD dimerization to Binor‐S using cationic Co^I, Rh^I, and Ir^I catalysts, using mPW1PW91, mPW1K, and B3LYP density functional methods. Our results indicate that the monomeric metal center has the ability to bind with four double bonds of two NBD molecules with a syn spatial geometry to form a penta‐coordinated complex. We designed three possible pathways, but found two of them blocked. The favored pathway involves three steps from the reactant precursor to the product precursor: the first step is the formation of a single bond to connect two NBD units, the second is the alkene insertion leading to the formation of the three‐membered ring structure, and the final step is the formation of the final product precursor. Orbital analysis showed metal…C? C σ agostic interaction in the product precursor, which is in agreement with the previous experimental findings. In addition, we found that the solvent and counter‐ions had significant effects on the dimerization reactions. © 2010 Wiley Periodicals, Inc. J Comput Chem 2010     

Nature of the ring-closure process along the rearrangement of octa-1,3,5,7-tetraene to cycloocta-1,3,5-triene from the perspective of the electron localization function and catastrophe theory

We analyze the behavior of the energy profile of the ring‐closure process for the transformation of (3Z,5Z)‐octa‐1,3,5,7‐tetraene 5 to (1Z,3Z,5Z)‐cycloocta‐1,3,5‐triene 6 through a combination of electron localization function (ELF) and catastrophe theory (CT). From this analysis, concepts such as bond breaking/forming processes, formation/annihilation of lone pairs, and other electron pair rearrangements arise naturally through the reaction progress simply in terms of the different ways of pairing up the electrons. A relationship between the topology and the nature of the bond breaking/forming processes along this rearrangement is reported. The different domains of structural stability of the ELF occurring along the intrinsic reaction path have been identified. The reaction mechanism consists of six steps separated by fold and cusp catastrophes. The transition structure is observed in the third step, d(C1? C8) = 2.342 Å, where all bonds have topological signature of single bonds (C? C). The “new” C1? C8 single bond is not formed in transition state and respective catastrophe of the ELF field (cusp) is localized in the last step, d(C1? C8) ≈ 1.97 Å, where the two monosynaptic nonbonding basins V(C1) and V(C8) are joined into single disynaptic bonding basin V(C1,C8). The V(C1,C8) basin corresponds to classical picture of the C1? C8 bond in the Lewis formula. In cycloocta‐1,3,5‐triene 6 the single C1? C8 bond is characterized by relatively small basin population 1.72e, which is much smaller than other single bonds with 2.03 and 2.26e. © 2011 Wiley Periodicals, Inc. J Comput Chem, 2011     

Transition Metal Complexation of σ Bonds

The first σ complexes were found in the 1960s and 1970s, but they did not attract more than passing attention. Only now are we beginning to recognize their key role in the chemical reactions of σ bonds, and this has encouraged more detailed study. In contrast with the more familiar π-donor complexes such as M? (CH₂?CH₂) and complexes like M? NH₃, in which the one pair of electrons on the N atom is bound to the metal atom, in a σ complex an X? H group binds to the transition metal atom; the X? H σ-bonding electron pair acts as a 2e donor to give an (X-H)-M type complex. Dihydrogen complexes (X = H) are one important group of σ complexes. C-H-M complexes (X = R₃C) with an agostic C-H-M interaction have not only been found in the ground state but also implicated in the transition states of many important organometallic transformations such as Ziegler–Natta catalysis and sigma bond metathesis. The importance of X? H bond activation will encourage continued growth in this field.     

New insight into the electronic structure of iron(IV)‐oxo porphyrin compound I. A quantum chemical topological analysis

The electronic structure of iron‐oxo porphyrin π‐cation radical complex Por^·+Fe^IV?O (S? H) has been studied for doublet and quartet electronic states by means of two methods of the quantum chemical topology analysis: electron localization function (ELF) η(r) and electron density ρ(r). The formation of this complex leads to essential perturbation of the topological structure of the carbon–carbon bonds in porphyrin moiety. The double C?C bonds in the pyrrole anion subunits, represented by pair of bonding disynaptic basins V_i=1,2(C,C) in isolated porphyrin, are replaced by single attractor V(C,C)_i=1–20 after complexation with the Fe cation. The iron–nitrogen bonds are covalent dative bonds, N→Fe, described by the disynaptic bonding basins V(Fe,N)_i=1–4, where electron density is almost formed by the lone pairs of the N atoms. The nature of the iron–oxygen bond predicted by the ELF topological analysis, shows a main contribution of the electrostatic interaction, Fe^δ+···O^δ?, as long as no attractors between the C(Fe) and C(O) core basins were found, although there are common surfaces between the iron and oxygen basines and coupling between iron and oxygen lone pairs, that could be interpreted as a charge‐shift bond. The Fe? S bond, characterized by the disynaptic bonding basin V(Fe,S), is partially a dative bond with the lone pair donated from sulfur atom. The change of electronic state from the doublet (M = 2) to quartet (M = 4) leads to reorganization of spin polarization, which is observed only for the porphyrin skeleton (?0.43e to 0.50e) and S? H bond (?0.55e to 0.52e). © 2012 Wiley Periodicals, Inc.     

A density functional theory study on the reactions of Y atom and Y+ cation with carbonyl sulfide

The mechanism of the title reactions have been studied by using the DFT (B3LYP/ECP/6‐311+G*) level of theory. Both ground and excited state potential energy surfaces are discussed. It is found the reaction mechanism is insertion mechanism both along the C? S and C? O bond activation branches, but the C? S bond activation is much more favorable in energy than the C? O bond activation. The reaction of Y atom with SCO was shown to occur preferentially on the ground state (doublet) PES throughout the reaction process, and the experimentally observed species, have been explained according to the mechanism revealed in this work. Different from that of Y + SCO system, the reaction between Y⁺ cation and SCO involves potential energy curve‐crossing which dramatically affects reaction mechanism. Due to the intersystem crossing existing in the reaction process of Y⁺ with SCO, the intermediates SY⁺(η²CO) and OY⁺(η²CS) may not form. All our theoretical results not only support the existing conclusions inferred from early experiment, but also complement the pathway and mechanism for this reaction. © 2009 Wiley Periodicals, Inc. Int J Quantum Chem, 2010     

Carbon Monoxide Bonding With BeO and BeCO3: Surprisingly High CO Stretching Frequency of OCBeCO3

The complexes OCBeCO₃ and COBeCO₃ have been isolated in a low‐temperature neon matrix. The more stable isomer OCBeCO₃ has a very high C? O stretching mode of 2263 cm^?1, which is blue‐shifted by 122 cm^?1 with respect to free CO and 79 cm^?1 higher than in OCBeO. Bonding analysis of the complexes shows that OCBeO has a stronger OC? BeY bond than OCBeCO₃ because it encounters stronger π backdonation. The isomers COBeCO₃ and COBeO exhibit red‐shifted C? O stretching modes with respect to free CO. The inverse change of C? O stretching frequency in OC? BeY and CO? BeY is explained with the reversed polarization of the σ and π bonds in CO.     

Theoretical study of methoxy group influence in the gas‐phase elimination kinetics of methoxyalkyl chlorides

The unimolecular gas‐phase elimination kinetics of 2‐methoxy‐1‐chloroethane, 3‐methoxy‐1‐chloropropane, and 4‐methoxyl‐1‐chloroburane has been studied by using density functional theory (DFT) methods to propose the most reasonable mechanisms of decomposition of the aforementioned compounds. Calculation results of 2‐methoxy‐1‐chloroethane and 3‐methoxy‐1‐chloropropane suggest dehydrochlorination through a concerted nonsynchronous four‐centered cyclic transition state (TS) to give the corresponding olefin. In the case of 4‐methoxyl‐1‐chloroburane, in addition to the 1,2‐elimination mechanism, the anchimeric assistance by the methoxy group, through a polar five‐centered cyclic TS, provides additional pathways to give 4‐methoxy‐butene, tetrahydrofuran and chloromethane. The bond polarization of the C? Cl, in the direction of C^δ+ ··· Cl^δ?, is the limiting step of these elimination reactions. The significant increase in rate together with the formation of a cyclic product tetrahydrofuran in the gas‐phase elimination of 4‐methoxyl‐1‐chloroburane is attributed to neighboring group participation of the oxygen of the methoxy group in the TS. The theoretical calculations show a good agreement with the reported experimental results. © 2011 Wiley Periodicals, Inc. Int J Quantum Chem, 2011     

Oxidation of cellulose fibers mediated by nonpersistent nitroxyl radicals

Three different nonpersistent radicals bearing >NO^? moiety have been used to oxidize the viscose fibers at room temperature and alkaline pH. The generation of the free radical species was achieved in situ, from their corresponding ? OH precursors: 1‐hydroxybenzotriazole, violuric acid, and N‐hydroxy‐3,4,5,6‐tetraphenylphthalimide. Three different routes were used: (i) in the presence of metallic cocatalyst (lead tetraacetate), (ii) under metal‐free conditions (anthraquinone as organic cocatalyst), and (iii) a homolytically scission of ? OH bond through a 365‐nm UV irradiation. The oxidized fibers were subjected to attenuated total reflection FTIR characterization, potentiometric titration, wide angle X‐ray, energy dispersive X‐ray analyses, microscopic investigations, and solid‐state ¹³C‐NMR. The patterns of CP/MAS ¹³C‐NMR spectra revealed that the oxidation occurred at the C6 primary hydroxyl groups of cellulose. Notably, as a result of the introduction of hydrophilic carboxylate groups, the water retention values of the oxidized fibers increased by 35% as compared to the original nonoxidized sample. © 2010 Wiley Periodicals, Inc. J Polym Sci Part A: Polym Chem, 2010     

NMR and X‐ray structural studies on 3‐benzyl‐8‐bromoadenine

8‐Bromoadenine was benzylated in the presence of base to give a mixture of two regioisomers. One was easily recognized as 9‐benzyl‐8‐bromoadenine, but the other structure could not be determined with absolute certainty by NMR. Therefore, X‐ray crystallography was used to prove that the benzyl group was attached to N‐3. Furthermore, it is shown that the 3‐benzyl adenine derivative exists as the amine tautomer both in the crystalline state as well as in solution (DMSO‐d₆), with restricted rotation around the N⁶? C6 bond. J. Heterocyclic Chem., (2011).     

Electronic states and nature of the chemical bond in the molecule CrC by all-electron ab initio calculations

All-electron ab initio Hartree–Fock (HF ), valence configuration interaction (CI ), and multiconfiguration self-consistent-field (CASSCF ) calculations have been applied to investigate the electronic states of the CrC molecule. The molecule is predicted as having four low-lying electronic states, ³∑^?, ⁵∑^?, ⁷∑^?, and ⁹∑^?, separated by an energy gap of 0.55 eV from the next higher-lying state, ¹∑^?, which is followed by the states ⁵Π and ⁷Π. The four lowest-lying electronic states are due to the coupling of the angular momenta of the ⁶S_g Cr⁺ ion with those of the ⁴S_u C^? anion. The chemical bond in the ³∑^? ground state can be viewed as a quadruple bond composed of two σ and two π bonds. One σ bond is due to the formation of a molecular orbital that is doubly occupied. The remaining bonds, i.e., one σ and two π bonds, arise from valence-bond couplings. The π bonds originate from the valence-bond couplings of the electrons in the C 2pπ orbitals with those in the Cr 3dπ orbitals. The σ bond originates from the valence-bond coupling of the C 2pσ electron with an electron in the Cr 4s, 4p hybrid that is polarized away from the C atom.     

Probing the weakest bond and the cleavage of p‐chlorobenzaldehyde diperoxide energetic molecule via quantum chemical calculations and theoretical charge density analysis

A high‐level ab initio Hartree‐Fock/Møller‐Plesset 2 and density functional theory quantum chemical calculations were performed on p‐chlorobenzaldehyde diperoxide energetic molecule to understand its bond topological, electrostatic, and energetic properties. The optimized molecular geometry for the basis set 6‐311G** exhibit chair diperoxide ring and planar aromatic side rings. Although the diperoxide ring bear same type of side rings, surprisingly, both the rings are almost perpendicular to each other, and the dihedral angle is 96.1°. The MP2 method predicts the O? O bond distance as ～1.466 Å. The charge density calculation reveals that the C? C bonds of chlorobenzaldehyde ring have rich electron density and the value is ～2.14 e Å^?3. The maximum electron density of the O? O bonds does not lie along the internuclear axes; in view of this, a feeble density is noticed in the ring plane. The high negative values of laplacian of C? C bonds (approximately ?22.4 e Å^?5) indicate the solidarity of these bonds, whereas it is found too small (approximately ?1.8 e Å^?5 for MP2 calculation) in O? O bonds that shows the existence of high degree of bond charge depletion. The energy density in all the C? C bonds are found to be uniform. A high electronegative potential region is found at the diperoxide ring which is expected to be a nucleophilic attack area. Among the bonds, the O? O bond charge is highly depleted and it also has high bond kinetic energy density; in consequence of this, the molecular cleavage is expected to happen across these bonds when the material expose to any external stimuli such as heat or pressure treatment. © 2010 Wiley Periodicals, Inc. Int J Quantum Chem, 2011     

CH⋅⋅⋅N Hydrogen‐Bonding Interaction in 7‐Azaindole:CHX3 (X=F,Cl) Complexes

The C?H???Y (Y=hydrogen‐bond acceptor) interactions are somewhat unconventional in the context of hydrogen‐bonding interactions. Typical C?H stretching frequency shifts in the hydrogen‐bond donor C?H group are not only small, that is, of the order of a few tens of cm^?1, but also bidirectional, that is, they can be red or blue shifted depending on the hydrogen‐bond acceptor. In this work we examine the C?H???N interaction in complexes of 7‐azaindole with CHCl₃ and CHF₃ that are prepared in the gas phase through supersonic jet expansion using the fluorescence depletion by infra‐red (FDIR) method. Although the hydrogen‐bond acceptor, 7‐azaindole, has multiple sites of interaction, it is found that the C?H???N hydrogen‐bonding interaction prevails over the others. The electronic excitation spectra suggest that both complexes are more stabilized in the S₁ state than in the S₀ state. The C?H stretching frequency is found to be red shifted by 82 cm^?1 in the CHCl₃ complex, which is the largest redshift reported so far in gas‐phase investigations of 1:1 haloform complexes with various substrates. In the CHF₃ complex the observed C?H frequency is blue shifted by 4 cm^?1. This is at variance with the frequency shifts that are predicted using several computational methods; these predict at best a redshift of 8.5 cm^?1. This discrepancy is analogous to that reported for the pyridine‐CHF₃ complex [W. A. Herrebout, S. M. Melikova, S. N. Delanoye, K. S. Rutkowski, D. N. Shchepkin, B. J. van der Veken, J. Phys. Chem. A­ 2005 , 109, 3038], in which the blueshift is termed a pseudo blueshift and is shown to be due to the shifting of levels caused by Fermi resonance between the overtones of the C?H bending and stretching modes. The dissociation energies, (D₀), of the CHCl₃ and CHF₃ complexes are computed (MP2/aug‐cc‐pVDZ level) as 6.46 and 5.06 kcal mol^?1, respectively.     

The effect of substitution on the intramolecular hydrogen bonding in 3‐hydroxy‐propenethial

The intramolecular hydrogen bond strength of 3‐hydroxy‐propenethial (HPT) as well as the fluoro, chloro, bromo, and methyl derivatives were investigated at the B3LYP/6‐311++G^** level of theory. Solvent‐based calculations (in water) for HPT and derivatives were also carried out. The nature of the intramolecular hydrogen bond existing within the molecular under investigation has been studied by means of the Bader theory of atoms in molecules (AIM) that is based upon the use topological properties in terms of the electron density. © 2008 Wiley Periodicals, Inc. Int J Quantum Chem, 2008     

Carbon?hydrogen versus carbon?heteroatom activation by a high‐valent zirconium‐imido complex

A density functional theory (DFT) study of carbon? hydrogen versus carbon? heteroatom bond activation is presented. Heteroatom groups (X) investigated include X = F, Cl, OH, SH, NH₂, PH₂. The activating model complex is a prototypical d⁰ zirconium‐imide. While C? X activation has a thermodynamic advantage over C? H activation, the former has been found to have a kinetic advantage. Implications for catalytic hydrocarbon functionalization and phosphine–ligand degradation are discussed. The present results for a high‐valent metal complex are compared/contrasted with low‐valent bond activating complexes. © 2006 Wiley Periodicals, Inc. Int J Quantum Chem, 2006     

Density functional theory study of the oxidation of methanol to formaldehyde on a hydrated vanadia cluster

Density functional theory was used to study the mechanism for the oxidation of methanol to formaldehyde. A vanadium oxide cluster O?V(OH)₃ has been utilized to represent the catalytic system under hydrated conditions, i.e., in the presence of V? OH hydroxyl groups. Two types of methoxy‐intermediates have been considered: a penta‐coordinate methoxy‐intermediate (OH)₄V(OCH₃) and a tetrahedral methoxy‐intermediate (OH)₂VO(OCH₃)(H₂O). The most plausible reaction pathway corresponds to the process involving first the formation of the tetrahedral methoxide, and a subsequent rate‐limiting step where hydrogen is transferred from the methoxy groups toward the oxygen atom of the vanadyl V?O site. The reaction mechanism is a typical two‐state reactivity process due to a change of the multiplicity (reactive singlet → product triplet) along the reaction coordinate accompanied by a reduction of the vanadium center from V^V (d⁰) to V^III (d²). Minimum energy crossing points were localized and possible spin inversion processes are discussed by means of the intrinsic reaction coordinate approach to find the most favorable reaction pathways. The hydration effect is found to be mainly the destabilization of the methoxy intermediates. An alternative reaction pathway with a lower apparent barrier is presented. © 2010 Wiley Periodicals, Inc. J Comput Chem, 2010