Neutral Nickel Ethylene Oligo‐ and Polymerization Catalysts: Towards Computational Catalyst Prediction and Design

DFT calculations have been used to elucidate the chain termination mechanisms for neutral nickel ethylene oligo‐ and polymerization catalysts and to rationalize the kind of oligomers and polymers produced by each catalyst. The catalysts studied are the (κ²‐O,O)‐coordinated (1,1,1,5,5,5‐hexafluoro‐2,4‐acetylacetonato)nickel catalyst I , the (κ²‐P,O)‐coordinated SHOP‐type nickel catalyst II , the (κ²‐N,O)‐coordinated anilinotropone and salicylaldiminato nickel catalysts III and IV , respectively, and the (κ²‐P,N)‐coordinated phosphinosulfonamide nickel catalyst V . Numerous termination pathways involving β‐H elimination and β‐H transfer steps have been investigated, and the most probable routes identified. Despite the complexity and multitude of the possible termination pathways, the information most critical to chain termination is contained in only few transition states. In addition, by consideration of the propagation pathway, we have been able to estimate chain lengths and discriminate between oligo‐ and polymerization catalysts. In agreement with experiment, we found the Gibbs free energy difference between the overall barrier for the most facile propagation and termination pathways to be close to 0 kcal mol^?1 for the ethylene oligomerization catalysts I and V , whereas values of at least 7 kcal mol^?1 in favor of propagation were determined for the polymerization catalysts III and IV . Because of the shared intermediates between the termination and branching pathways, we have been able to identify the preferred cis/trans regiochemistry of β‐H elimination and show that a pronounced difference in σ donation of the two bridgehead atoms of the bidentate ligand can suppress hydride formation and thus branching. The degree of rationalization obtained here from a handful of key intermediates and transition states is promising for the use of computational methods in the screening and prediction of new catalysts of the title class.     

Computational Study of the Mechanism and Selectivity of Palladium‐Catalyzed Propargylic Substitution with Phosphorus Nucleophiles

The mechanism and sources of selectivity in the palladium‐catalyzed propargylic substitution reaction that involves phosphorus nucleophiles, and which yields predominantly allenylphosphonates and related compounds, have been studied computationally by means of density functional theory. Full free‐energy profiles are computed for both H‐phosphonate and H‐phosphonothioate substrates. The calculations show that the special behavior of H‐phosphonates among other heteroatom nucleophiles is indeed reflected in higher energy barriers for the attack on the central carbon atom of the allenyl/propargyl ligand relative to the ligand‐exchange pathway, which leads to the experimentally observed products. It is argued that, to explain the preference of allenyl‐ versus propargyl‐phosphonate/phosphonothioate formation in reactions that involve H‐phosphonates and H‐phosphonothioates, analysis of the complete free‐energy surfaces is necessary, because the product ratio is determined by different transition states in the respective branches of the catalytic cycle. In addition, these transition states change in going from a H‐phosphonate to a H‐phosphonothioate nucleophile.     

A density functional theory study on the electrocyclization of 1,2,4,6-heptatetraene analogues: converting a pericyclic to a pseudopericyclic reaction

A comprehensive B3LYP/6-31+G* study on the electrocyclization of 1,2,4,6-heptatetraene analogues was conducted. Starting from the cyclization of (2Z)-2,4,5-hexatrienal, a pericyclic disrotatory process favored by the assistance of a electron lone pair, we incorporated small modifications in its molecular structure to obtain a truly pseudopericyclic process. To this purpose electronegative atoms (fluorine and nitrogen) were added to give a more electrophilic character on the carbon atom which is attacked by the electron lone pair of the oxygen atom. The complete pathway for each reaction was determined, and changes in magnetic properties were monitored with a view to estimating the aromatization associated with each process. This information, together with the energetic and structural results, allowed us to classify the reactions as pseudopericyclic or pericyclic. Among all studied reactions only one was a truly pseudopericyclic process and another was a borderline case. The features of this unequivocally pseudopericyclic case were analyzed in depth.     

Modeling a halogen dance reaction mechanism: A density functional theory study

Since the discovery of the halogen dance (HD) reaction more than 60 years ago, numerous insights into the mechanism have been unveiled. To date however, the reaction has not been investigated from a theoretical perspective. Density functional theory (DFT) was used to model the potential energy surface linking the starting reagents to the lithiated products for each step in the mechanism using a thiophene substrate. It was found that the lithium‐halogen exchange mechanism is critical to understand the HD mechanism in detail and yielded the knowledge that S_N2 transition states (TS) are favored over the four‐center type for the lithium‐bromine exchange steps. The overall driving force for the HD is thermodynamics, while the kinetic factors tightly control the reaction path through temperature. The S_N2 lithium‐bromide TS are barrierless, except the second, which is the limiting step. Finally, the model for the HD is discovered to be a pseudo‐clock type, due to a highly favorable bromide catalysis step and the reformation of 2‐bromothiophene. © 2016 Wiley Periodicals, Inc.     

Photophysical Properties of Benzoylgermane and para‐Substituted Derivatives: Substituent Effects on Electronic Transitions

In the present study, a selection of basic substitution patterns on benzoyl(trimethyl)germane was investigated using time‐dependent density‐functional theory (TDDFT) to explore the influence on the stability and on the relative order of the lowest excited electronic states. The theoretical results are in agreement with absorption and fluorescence measurements. We show that electron‐withdrawing groups decrease the energetic level of the lowest singlet and triplet state relative to the electron‐pushing systems resulting in red‐shifted radiative transitions (fluorescence). In the first triplet state electron‐withdrawing groups lead to an increased dissociation barrier and a close approach with the singlet ground state before the transition state in the triplet state is reached, favoring radiationless ground‐state recovery. The results are also in good agreement with empirical concepts of organic chemistry, therefore providing simple rules for synthetic strategies towards tuning the excited‐state properties of benzoylgermanes.     

A quantum chemical study of racemization pathways in substituted chrysene derivatives

The potential-energy surfaces of 5,11-disubstituted 6,12-dimethoxychrysene and chrysene-6,12-dione derivatives were investigated by means of density functional calculations. We report relative energies of all conformers and an identification of the racemisation pathways of the chiral equilibrium structures. By analysis of homodesmotic reactions we were able to obtain an estimate for the strain energy of the substituted compounds. This strain energy can be used as a means of measuring the steric effects exerted by the substituents.     

Hydrolysis of tetrafluorosilane under neutral conditions: A quantum chemical study

A density functional quantum chemical study of the first step of SiF₄ hydrolysis under neutral conditions in the presence of one or two H₂O molecules was carried out. The reaction is endothermic and can follow three different pathways that involve the formation of penta-coordinated (pathway A) and hexacoordinated (pathways B and C) intermediates and transition states as the key steps. Pathway B is the most energetically favorable. All three pathways of the hydrolysis reaction lead to a product with formal retention of the configuration of substituents at the silicon atom. Published in Russian in Izvestiya Akademii Nauk. Seriya Khimicheskaya, No. 5, pp. 749–753, May, 2006.     

Criteria for the elucidation of the pseudopericyclic character of the cyclization of (Z)-1,2,4,6-heptatetraene and its heterosubstituted analogues: magnetic properties and natural bond orbital analysis

The electrocyclization of heterosubstituted derivatives of (Z)-1,2,4,6-heptatetraene, (2Z)-2,4,5-hexatrien-1-imine and (2Z)-2,4,5-hexatrienal exhibit some features which suggest a pseudopericyclic mechanism. In order to examine this, a comprehensive study including the determination of magnetic properties to estimate aromaticity and an NBO analysis throughout the reaction path was conducted. The cyclization of 5oxo-2,4-pentadienal, a process of unequivocal pseudopericyclic nature, was studied for comparison. The results suggest that, although the lone electron pair on the heteroatom in the heptatetraene derivatives seemingly plays a crucial role in the reaction mechanism, it does not suffice to deprive the reaction from the essential features of a pericyclic disrotatory electrocyclization.     

Enantioselective Synthesis of Piperidines through the Formation of Chiral Mixed Phosphoric Acid Acetals: Experimental and Theoretical Studies

An enantioselective intramolecular chiral phosphoric acid‐catalyzed cyclization of unsaturated acetals has been utilized for the synthesis of functionalized chiral piperidines. The chiral enol ether products of these cyclizations undergo subsequent in situ enantioenrichment through acetalization of the minor enantiomer. A new computational reaction exploration method was utilized to elucidate the mechanism and stereoselectivity of this transformation. Rather than confirming the originally postulated cyclization proceeding directly through a vinyl oxocarbenium ion, simulations identified an alternative two‐step mechanism involving the formation of a mixed chiral phosphate acetal, which undergoes a concerted, asynchronous S_N2′‐like displacement to yield the product with stereoselectivity in agreement with experimental observations.     

Carbon Dioxide Capture and Release by Anions with Solvent‐Dependent Behaviour: A Theoretical Study

Recently, Clyburne and co‐workers [Science, 2014 , 344, 75–78] reported the novel synthesis of the elusive cyanoformate anion, NCCO₂^?. The stability of this anion is dependent on the dielectric constant of the local environment (polarity‐switchable solvent): it is stable in low‐polarity media and unstable in high‐polarity solvents; hence, capturing and then releasing CO₂. The possibility of extending such behaviour to other anions is explored herein. Specifically, the CO₂ capture process is studied for 26 anions in the gas phase and 3 distinct solvents (water, tetrahydrofuran, and toluene) by using the polarisable continuum model. Calculations are performed with the M06‐2X and B3LYP‐D3 density functionals and the aug‐cc‐pVTZ basis set. The design of new CO₂ complexes with the anion, which can be formed or destroyed on demand by changing the solvent, is possible; the results for the alkoxylate and thiolate anions are especially promising. The nature of the substituents connected to the atom that bonds to CO₂ in the anion is crucial in modulating the relative stability of the products—a key point for reversibility in the CO₂ capture process. A moderate interaction for the anion–CO₂ adduct—about 10 kcal mol^?1 relative free energy with respect to the isolated reactants in the gas phase—and a relevant effect in the dielectric constant of the local environment are also key ingredients to achieve solvent dependency.     

The Catalytic Mechanism of Fluoroacetate Dehalogenase: A Computational Exploration of Biological Dehalogenation

The biological dehalogenation of fluoroacetate carried out by fluoroacetate dehalogenase is discussed by using quantum mechanical/molecular mechanical (QM/MM) calculations for a whole‐enzyme model of 10 800 atoms. Substrate fluoroacetate is anchored by a hydrogen‐bonding network with water molecules and the surrounding amino acid residues of Arg105, Arg108, His149, Trp150, and Tyr212 in the active site in a similar way to haloalkane dehalogenase. Asp104 is likely to act as a nucleophile to attack the α‐carbon of fluoroacetate, resulting in the formation of an ester intermediate, which is subsequently hydrolyzed by the nucleophilic attack of a water molecule to the carbonyl carbon atom. The cleavage of the strong C? F bond is greatly facilitated by the hydrogen‐bonding interactions between the leaving fluorine atom and the three amino acid residues of His149, Trp150, and Tyr212. The hydrolysis of the ester intermediate is initiated by a proton transfer from the water molecule to His271 and by the simultaneous nucleophilic attack of the water molecule. The transition state and produced tetrahedral intermediate are stabilized by Asp128 and the oxyanion hole composed of Phe34 and Arg105.     

Theoretical Investigation of the OH.‐Initiated Oxidation of Benzaldehyde in the Troposphere

A mechanistic and kinetic study of the OH^.‐initiated oxidation of benzaldehyde is carried out using quantum chemical methods and classical transition state theory. We calculate the rate constant for this reaction within the temperature range of 200–350 K at atmospheric pressure. All possible hydrogen abstraction and OH^. addition channels are considered and branching ratios are obtained. Tunneling corrections are taken into account for abstraction channels, assuming unsymmetrical Eckart barriers. The aldehydic abstraction is by far the most important reaction channel within the entire range of temperatures studied, especially at room temperature and lower—the temperatures relevant to atmospheric chemistry. The relative importance of all the other possible channels increases slightly with temperature. Branching ratios show that addition at the ring and abstraction of an ortho hydrogen contribute about 1 % each at about 300 K, while the branching ratio for the main reaction decreases from 99 % at 200 K to 93 % at 350 K. The results are compared with available experimental measurements.