Interactions of Molecules with cis and trans Double Bonds: A Theoretical Study of cis‐ and trans‐2‐Butene

Noncovalent interactions of cis‐ and trans‐2‐butene, as the smallest model systems of molecules with cis and trans double bonds, were studied to find potential differences in interactions of these molecules. The study was performed using quantum chemical methods including very accurate CCSD(T)/CBS method. We studied parallel and displaced parallel interactions in 2‐butene dimers, in butane dimers, and between 2‐butene and saturated butane. The results show the trend that interactions of 2‐butene with butane are the strongest, followed by interactions in butane dimers, whereas the interaction in 2‐butene dimers are the weakest. The strongest calculated interaction energy is between trans‐2‐butene and butane, with a CCSD(T)/CBS energy of ?2.80 kcal mol^?1. Interactions in cis‐2‐butene dimers are stronger than interactions in trans‐2‐butene dimers. Interestingly, some of the interactions involving 2‐butene are as strong as interactions in a benzene dimer. These insights into interactions of cis‐ and trans‐2‐butene can improve understanding of the properties and processes that involve molecules with cis and trans double bonds, such as fatty acids and polymers.     

Mechanistic Insights into Palladium‐Catalyzed Silylation of Aryl Iodides with Hydrosilanes through a DFT Study

The catalytic cycles of palladium‐catalyzed silylation of aryl iodides, which are initiated by oxidative addition of hydrosilane or aryl iodide through three different mechanisms characterized by intermediates R₃Si−Pd^II−H (Cycle A), Ar−Pd^II−I (Cycle B), and Pd^IV (Cycle C), have been explored in detail by hybrid DFT. Calculations suggest that the chemical selectivity and reactivity of the reaction depend on the ligation state of the catalyst and specific reaction conditions, including feeding order of substrates and the presence of base. For less bulky biligated catalyst, Cycle C is energetically favored over Cycle A, through which the silylation process is slightly favored over the reduction process. Interestingly, for bulky monoligated catalyst, Cycle B is energetically more favored over generally accepted Cycle A, in which the silylation channel is slightly disfavored in comparison to that of the reduction channel. Moreover, the inclusion of base in this channel allows the silylated product become dominant. These findings offer a good explanation for the complex experimental observations. Designing a reaction process that allows the oxidative addition of palladium(0) complex to aryl iodide to occur prior to that with hydrosilane is thus suggested to improve the reactivity and chemoselectivity for the silylated product by encouraging the catalytic cycle to proceed through Cycles B (monoligated Pd⁰ catalyst) or C (biligated Pd⁰ catalyst), instead of Cycle A.     

DFT Mechanistic Study on Diene Metathesis Catalyzed by Ru‐Based Grubbs–Hoveyda‐Type Carbenes: The Key Role of π‐Electron Density Delocalization in the Hoveyda Ligand

The catalytic activity and catalyst recovery of two heterogenized ruthenium‐based precatalysts ( H and NO₂(4) ) in diene ring‐closing metathesis have been studied by means of density functional calculations at the B3LYP level of theory. For comparison and rationalization of the key factors that lead to higher activities and higher catalyst recoveries, four other Grubbs–Hoveyda complexes have also been investigated. The full catalytic cycle (catalyst formation, propagation, and precatalyst regeneration) has been considered. DFT calculations suggest that either for the homogeneous and heterogenized systems the activity of the catalysts mainly depends on the ability of the precursor to generate the propagating carbene. This ability does not correlate with the traditionally identified key factor, the Ru???O interaction strength. In contrast, precatalysts with lower alkoxy‐dissociation energy barriers and lower stabilities compared with the propagating carbene also present larger C1? C2 bond length (i.e., lower π character of the C? C bond that exists between the metal–carbene (Ru?C) and the phenyl ring of the Hoveyda ligand). Catalyst recovery, regardless of whether a release–return mechanism occurs or not, is also mainly determined by the π delocalization. Therefore, future Grubbs–Hoveyda‐type catalyst development should be based on fine‐tuning the π‐electron density of the phenyl moiety, with the subsequent effect on the metalloaromaticity of the ruthenafurane ring, rather than considering the modification of the Ru???O interaction.     

Mechanistic Insights into the Rhenium‐Catalyzed Alcohol‐To‐Olefin Dehydration Reaction

Rhenium‐based complexes are powerful catalysts for the dehydration of various alcohols to the corresponding olefins. Here, we report on both experimental and theoretical (DFT) studies into the mechanism of the rhenium‐catalyzed dehydration of alcohols to olefins in general, and the methyltrioxorhenium‐catalyzed dehydration of 1‐phenylethanol to styrene in particular. The experimental and theoretical studies are in good agreement, both showing the involvement of several proton transfers, and of a carbenium ion intermediate in the catalytic cycle.     

Mechanistic insights into iridium-catalyzed asymmetric hydrogenation of dienes

Hydrogenation of 2,3-diphenylbutadiene (1) with the chiral carbene-oxazoline-iridium complex C has been studied by means of a combined experimental and computational approach. A detailed kinetic profile of the reaction was obtained with respect to consumption of the substrate and formation of the intermediate half-reduction products, 2,3-diphenylbut-1-ene (2) and the final product, 2,3-diphenylbutane (3). The data generated from these analyses, and from NMR experiments, revealed several facets of the reaction. After a brief induction period (presumably involving reduction of the cyclooctadiene ligand on C), the diene concentration declines in a zero-order process primarily to give monoene intermediates. When all the diene is consumed, the reaction accelerates and compound 3 begins to accumulate. Interestingly, the prevalent enantiomer of the monoene intermediate 2 is converted mostly to meso-3 so the enantioselectivity of the reaction appears to reverse. The reaction seems to be first-order with respect to the catalyst when the catalyst concentration is less than 0.0075 M; diffusion of hydrogen across the gas-liquid interface complicates the analysis at higher catalyst concentrations. Similarly, these diffusion effects complicated measurements of reaction rate versus applied pressure of dihydrogen; other factors like stir speed and flask geometry come into play under some, but not all, the conditions examined. Density functional theory (DFT) calculations, using the PBE method, were used to probe the reaction. These studies indicate a transoid-eta(4)-diene-dihydride complex forms in the first stages of the catalytic cycle. Further reaction requires dissociation of one alkene ligand to give a eta(2)-diene-dihydride-dihydrogen intermediate. A catalytic cycle that features Ir(3+)/Ir(5+) seems to be involved thereafter.     

Kinetic Stabilization and Reactivity of π Single‐Bonded Species: Effect of the Alkoxy Group on the Lifetime of Singlet 2,2‐Dialkoxy‐1,3‐diphenyloctahydropentalene‐1,3‐diyls

Kinetic stabilization and reactivity of π single‐bonded species have been investigated in detail by generating a series of singlet 2,2‐dialkoxy‐1,3‐diphenyloctahydropentalene‐1,3‐diyls ( DR s). The lifetime at 293 K in benzene was found to increase when the carbon chain length of the alkoxy groups was increased; 292 ns ( DRb ; OR=OR′=OCH₃) <880 ns ( DRc ; OR=OR′=OC₂H₅) <1899 ns ( DRd ; OR=OR′=OC₃H₇) ≈2292 ns ( DRe ; OR=OR′=OC₆H₁₃) ≈2146 ns ( DRf ; OR=OR′=OC₁₀H₂₁). DRh (OR=OC₃H₇, OR′=OCH₃; 935 ns) with the mixed‐acetal moiety is a longer‐lived species than another diastereomer DRg (OR=OCH₃, OR′=OC₃H₇; 516 ns). Activation parameters determined for the first‐order decay process reveal that the enthalpy factor plays a crucial role in determining the energy barrier of the ring‐closing reaction, that is, from the π‐bonding to the σ‐bonding compounds. Computational studies using density functional theory provided more insight into the structures of the singlet species with π single‐bonded character and the transition states for the ring‐closing reaction, thereby clarifying the role of the alkoxy group on the lifetime and the stereoselectivity of the ring‐closing reaction.     

Cooperative Catalysis of Metal and O?H⋅⋅⋅O/sp3‐C?H⋅⋅⋅O Two‐Point Hydrogen Bonds in Alcoholic Solvents: Cu‐Catalyzed Enantioselective Direct Alkynylation of Aldehydes with Terminal Alkynes

Catalyst–substrate hydrogen bonds in artificial catalysts usually occur in aprotic solvents, but not in protic solvents, in contrast to enzymatic catalysis. We report a case in which ligand–substrate hydrogen‐bonding interactions cooperate with a transition‐metal center in alcoholic solvents for enantioselective catalysis. Copper(I) complexes with prolinol‐based hydroxy amino phosphane chiral ligands catalytically promoted the direct alkynylation of aldehydes with terminal alkynes in alcoholic solvents to afford nonracemic secondary propargylic alcohols with high enantioselectivities. Quantum‐mechanical calculations of enantiodiscriminating transition states show the occurrence of a nonclassical sp³‐C? H???O hydrogen bond as a secondary interaction between the ligand and substrate, which results in highly directional catalyst–substrate two‐point hydrogen bonding.     

Intermolecular Interactions in Weak Donor–Acceptor Complexes from Symmetry‐Adapted Perturbation and Coupled‐Cluster Theory: Tetracyanoethylene–Benzene and Tetracyanoethylene–p‐Xylene

The interactions in the complexes of tetracyanothylene (TCNE) with benzene and p‐xylene, often classified as weak electron donor–acceptor (EDA) complexes, are investigated by a range of quantum chemical methods including intermolecular perturbation theory at the DFT‐SAPT (symmetry‐adapted perturbation theory combined with density functional theory) level and explicitly correlated coupled‐cluster theory at the CCSD(T)‐F12 level. The DFT‐SAPT interaction energies for TCNE–benzene and TCNE–p‐xylene are estimated to be ?35.7 and ?44.9 kJ mol^?1, respectively, at the complete basis set limit. The best estimates for the CCSD(T) interaction energy are ?37.5 and ?46.0 kJ mol^?1, respectively. It is shown that the second‐order dispersion term provides the most important attractive contribution to the interaction energy, followed by the first‐order electrostatic term. The sum of second‐ and higher‐order induction and exchange–induction energies is found to provide nearly 40 % of the total interaction energy. After addition of vibrational, rigid‐rotor, and translational contributions, the computed internal energy changes on complex formation approach results from gas‐phase spectrophotometry at elevated temperatures within experimental uncertainties, while the corresponding entropy changes differ substantially.     

Intramolecular Hydrogen Bonds in Low‐Molecular‐Weight Polyethylene Glycol

We used static DFT calculations to analyze, in detail, the intramolecular hydrogen bonds formed in low‐molecular‐weight polyethylene glycol (PEG) with two to five repeat subunits. Both red‐shifted O?H???O and blue‐shifting C?H???O hydrogen bonds, which control the structural flexibility of PEG, were detected. To estimate the strength of these hydrogen bonds, the quantum theory of atoms in molecules was used. Car–Parrinello molecular dynamics simulations were used to mimic the structural rearrangements and hydrogen‐bond breaking/formation in the PEG molecule at 300 K. The time evolution of the H???O bond length and valence angles of the formed hydrogen bonds were fully analyzed. The characteristic hydrogen‐bonding patterns of low‐molecular‐weight PEG were described with an estimation of their lifetime. The theoretical results obtained, in particular the presence of weak C?H???O hydrogen bonds, could serve as an explanation of the PEG structural stability in the experimental investigation.     

Synthesis,Electronic Properties,and Reactivity of Phospholes and 1,1′‐Biphospholes Bearing 2‐ or 3‐Thienyl C‐Substituents

PS, I love you! Novel mixed phosphole/thiophene π‐conjugated systems were synthesized and their electronic properties have been studied both experimentally by UV/Vis spectroscopy and electrochemistry and by theoretical calculations. Exploiting the chemistry of both P‐ and S‐heteroles allows the generation of a diverse range of novel ring‐fused benzophosphole–thiophene derivatives.

     

Sulfur Ylide Promoted Synthesis of N‐Protected Aziridines: A Combined Experimental and Computational Approach

A range of N‐protected aziridines [N‐Tosyl (N‐Ts), N‐2‐trimethylsilylethanesulfonamide (N‐SES), N‐tert‐butoxycarbonylamido (N‐Boc), and N‐o‐nitrobenzenesulfonamide (oNs)] were prepared in moderate to good yield and with high enantiomeric excess of both isomers starting from N‐protected imines, using a sulfonium salt derived from Eliel’s oxathiane. The diastereoselectivities of the reactions are influenced by the imine N‐protecting group, the imine substituent, and the sulfide structure. An unusual cis selectivity was observed in the formation of N‐tosyl‐2‐phenyl‐3‐tert‐butylaziridine and N‐o‐trimethylsilylethanesulfonamide‐2‐phenyl‐3‐tert‐butylaziridine, which was explained by using computational models. The analysis suggests that betaine formation in the case of N‐tosyl‐tert‐butylaldimine aziridination using oxathiane benzyl sulfonium ylide 1 ′ is reversible and that the selectivity is determined at the rotation step, which is unusual for semistabilized ylide aziridination. We have shown herein that the steric bulk of an imine substituent, in combination with a sterically demanding sulfonium ylide, can also affect the reversibility of the reaction. This is the first example of this sort involving aziridinations using semistabilized ylides.     

Does Silica Surface Catalyse Peptide Bond Formation? New Insights from First‐Principles Calculations

The role that silica surface could have played in prebiotic chemistry as a catalyst for peptide bond formation has been addressed at the B3LYP/6-31+G(d,p) level for a model reaction involving glycine and ammonia on a silica cluster mimicking an isolated terminal silanol group present at the silica surface. Hydrogen-bond complexation between glycine and the silanol is followed by the formation of the mixed surface anhydride Si(surf)-O-C(=O)-R, which has been suggested in the literature to activate the C=O bond towards nucleophilic attack by a second glycine molecule, here simulated by the simpler NH3 molecule. However, B3LYP/6-31+G(d,p) calculations show that formation of the surface mixed anhydride Si(surf)-O-C(=O)-R is disfavoured (delta(r)G298 approximately 6 kcal mol(-1)), and that the surface bond only moderately lowers the free-energy barrier of the nucleophilic attack responsible for peptide bond formation (deltaG298(double dagger) approximately 48 kcal mol(-1)) in comparison with the uncatalysed reaction (deltaG298(double dagger) approximately 52 kcal mol(-1)). A further decrease of the free-energy barrier of peptide bond formation (deltaG298(double dagger) approximately 41 kcal mol(-1)) is achieved by a single water molecule close to the reaction centre acting as a proton-transfer helper in the activated complex. A possible role of strained silica surface defects on the formation of the surface mixed anhydride Si(surf)-O-C(=O)-R has also been addressed.     

Base‐Free Non‐Noble‐Metal‐Catalyzed Hydrogen Generation from Formic Acid: Scope and Mechanistic Insights

The iron‐catalyzed dehydrogenation of formic acid has been studied both experimentally and mechanistically. The most active catalysts were generated in situ from cationic Fe^II/Fe^III precursors and tris[2‐(diphenylphosphino)ethyl]phosphine ( 1 , PP₃). In contrast to most known noble‐metal catalysts used for this transformation, no additional base was necessary. The activity of the iron catalyst depended highly on the solvent used, the presence of halide ions, the water content, and the ligand‐to‐metal ratio. The optimal catalytic performance was achieved by using [FeH(PP₃)]BF₄/PP₃ in propylene carbonate in the presence of traces of water. With the exception of fluoride, the presence of halide ions in solution inhibited the catalytic activity. IR, Raman, UV/Vis, and EXAFS/XANES analyses gave detailed insights into the mechanism of hydrogen generation from formic acid at low temperature, supported by DFT calculations. In situ transmission FTIR measurements revealed the formation of an active iron formate species by the band observed at 1543 cm^?1, which could be correlated with the evolution of gas. This active species was deactivated in the presence of chloride ions due to the formation of a chloro species (UV/Vis, Raman, IR, and XAS). In addition, XAS measurements demonstrated the importance of the solvent for the coordination of the PP₃ ligand.     

Mechanistic Investigations into the Enantioselective Conia‐Ene Reaction Catalyzed by Cinchona‐Derived Amino Urea Pre‐Catalysts and CuI

The enantioselective Conia‐ene cyclization of alkyne‐tethered β‐ketoesters is efficiently catalyzed by the combination of cinchona‐derived amino‐urea pre‐catalysts and copper(I) salts. The reaction scope is broad and a series of substrates can be efficiently cyclized with high yields and enantioselectivities. Herein, we present a detailed mechanistic study based on experimental considerations and quantum mechanical calculations. Several variables, such as the nature of the organic pre‐catalyst and the metal‐ion source, have been thoroughly investigated. Kinetic studies, as well as kinetic isotope effects and deuterium labeling experiments have been used to gain further insights into the mechanism and prove the cooperative nature of the catalytic system. Our studies suggest that the rate‐limiting step for the reaction involves the β‐ketoester deprotonation and that the active species responsible for the enantiodeterming step is monomeric in amino‐urea pre‐catalyst. Computational studies provide a quantitative understanding of the observed stereoinduction and identify hydrogen bonding from the urea group as a crucial factor in determining the observed enantioselectivity.     

Synthesis and Characterization of 6,13‐Diamino‐Substituted Pentacenes

A series of 6,13‐diamino‐substituted pentacenes 1 a – d has been prepared and characterized as a new class of pentacene derivatives with strong donor ability and enhanced solubility in common organic solvents. The spectroelectrochemical and DFT studies revealed that the two‐electron oxidation process was accompanied by the substantial structural change into a butterfly‐like conformation of the pentacene moiety. More importantly, the extent of deformation from the planar pentacene moiety in the dications of 6,13‐diaminopentacene is tunable by varying the N‐substituents.     

σ‐Hole Bond Versus Hydrogen Bond: From Tetravalent to Pentavalent N,P, and As Atoms

Ab initio calculations were performed on complexes of ZH₄⁺ (Z=N, P, As) and their fluoro derivatives, ZFH₃⁺ and ZF₄⁺, with a HCN (or LiCN) molecule acting as the Lewis base through the nitrogen electronegative center. It was found that the complexes are linked by the Z? H???N hydrogen bond or another type of noncovalent interaction in which the tetravalent heavy atom of the cation acts as the Lewis acid center, that is, when the Z???N link exists, which may be classified as the σ‐hole bond. The formation of the latter interaction is usually preferable to the formation of the corresponding hydrogen bond. The Z???N interaction may be also considered as the preliminary stage of the S_N2 reaction. This is supported by the observation that for a short Z???N contact, the corresponding complex geometry coincides with the trigonal‐bipyramidal geometry typical for the transition state of the S_N2 reaction. The Z???N interaction for some of complexes analyzed here possesses characteristics typical for covalent bonds. Numerous interrelations between geometrical, topological and energetic parameters are discussed. The natural bond orbital method as well as the Quantum Theory of “Atoms in Molecules” is applied to characterize interactions in the analyzed complexes. The experimental evidences of the existence of these interactions, based on the Cambridge Structure Database search, are also presented. In addition, it is justified that mechanisms of the formation of the Z???N interactions are similar to the processes occurring for the other noncovalent links. The formation of Z???N interaction as well as of other interactions may be explained with the use of the σ‐hole concept.