Redox Switches for Single‐Molecule Magnet Activity: An Ab Initio Insight

A dinuclear Co^II complex ( 1 ) featuring unprecedented anodic and cathodic switches for single‐molecule magnet (SMM) activity has been recently investigated (J. Am. Chem. Soc. 2013 , 135, 14670). The presence of sandwiched radicals in different oxidation states of this compound mediates magnetic coupling between the high‐spin (S=3/2) cobalt ions, which gives rise to SMM activity in both the oxidized ([ 1 (OEt₂)]⁺) and reduced ([ 1 ]^?) states. This feature represents the first example of a SMM exhibiting fully reversible, dual ON/OFF switchability. Here we apply ab initio and broken‐symmetry DFT calculations to elucidate the mechanisms responsible for magnetic properties and magnetization blocking in these compounds. It is found that due to the strong delocalization of the magnetic molecular orbital, there is a strong antiferromagnetic interaction between the radical and cobalt ions. The lack of high axiality of the cobalt centres explains why these compounds possess slow relaxation of magnetization only in an applied dc magnetic field.     

Ab Initio Molecular Dynamics with Explicit Solvent Reveals a Two‐Step Pathway in the Frustrated Lewis Pair Reaction

The role solvent plays in reactions involving frustrated Lewis pairs (FLPs)—for example, the stoichiometric mixture of a bulky Lewis acid and a bulky Lewis base—still remains largely unexplored at the molecular level. For a reaction of the phosphorus/boron FLP and dissolved CO₂ gas, first principles (Born–Oppenheimer) molecular dynamics with explicit solvent reveals a hitherto unknown two‐step reaction pathway—one that complements the concerted (one‐step) mechanism known from the minimum‐energy‐path calculations. The rationalization of the discovered reaction pathway—that is, the stepwise formation of P?C and O?B bonds—is that the environment (typical organic solvents) stabilizes an intermediate which results from nucleophilic attack of the phosphorus Lewis base on CO₂. This finding is significant because presently the concerted reaction‐path paradigm predominates in the rationalization of FLP reactivity. Herein we point out how to attain experimental proof of our results.     

Modulation of Stacking Interactions by Transition‐Metal Coordination: Ab Initio Benchmark Studies

A series of ab initio calculations are used to determine the C? H???π and π???π‐stacking interactions of aromatic rings coordinated to transition‐metal centres. Two model complexes have been employed, namely, ferrocene and chromium benzene tricarbonyl. Benchmark data obtained from extrapolation of MP2 energies to the basis set limit, coupled with CCSD(T) correction, indicate that coordinated aromatic rings are slightly weaker hydrogen‐bond acceptors but are significantly stronger hydrogen‐bond donors than uncomplexed rings. It is found that π???π stacking to a second benzene is stronger than in the free benzene dimer, especially in the chromium case. This is assigned, by use of energy partitioning in the local correlation method, to dispersion interactions between metal d and benzene π orbitals. The benchmark data is also used to test the performance of more efficient theoretical methods, indicating that spin‐component scaling of MP2 energies performs well in all cases, whereas various density functionals describe some complexes well, but others with errors of more than 1 kcal mol^?1.     

Synthesis and Reactivity of Six‐Membered Oxa‐Nickelacycles: A Ring‐Opening Reaction of Cyclopropyl Ketones

Cyclopropanecarboxaldehyde ( 1 a ), cyclopropyl methyl ketone ( 1 b ), and cyclopropyl phenyl ketone ( 1 c ) were reacted with [Ni(cod)₂] (cod=1,5‐cyclooctadiene) and PBu₃ at 100 °C to give η²‐enonenickel complexes ( 2 a – c ). In the presence of PCy₃ (Cy=cyclohexyl), 1 a and 1 b reacted with [Ni(cod)₂] to give the corresponding μ‐η²:η¹‐enonenickel complexes ( 3 a , 3 b ). However, the reaction of 1 c under the same reaction conditions gave a mixture of 3 c and cyclopentane derivatives ( 4 c , 4 c′ ), that is, a [3+2] cycloaddition product of 1 c with (E)‐1‐phenylbut‐2‐en‐1‐one, an isomer of 1 c . In the presence of a catalytic amount of [Ni(cod)₂] and PCy₃, [3+2] homo‐cycloaddition proceeded to give a mixture of 4 c (76 %) and 4 c′ (17 %). At room temperature, a possible intermediate, 6 c , was observed and isolated by reprecipitation at ?20 °C. In the presence of 1,3‐bis(2,6‐diisopropylphenyl)imidazol‐2‐ylidene (IPr), both 1 a and 1 c rapidly underwent oxidative addition to nickel(0) to give the corresponding six‐membered oxa‐nickelacycles ( 6 ai , 6 ci ). On the other hand, 1 b reacted with nickel(0) to give the corresponding μ‐η²:η¹‐enonenickel complex ( 3 bi ). The molecular structures of 6 ai and 6 ci were confirmed by X‐ray crystallography. The molecular structure of 6 ai shows a dimeric η¹‐nickelenolate structure. However, the molecular structure of 6 ci shows a monomeric η¹‐nickelenolate structure, and the nickel(II) 14‐electron center is regarded as having “an unusual T‐shaped planar” coordination geometry. The insertion of enones into monomeric η¹‐nickelenolate complexes 6 c and 6 ci occurred at room temperature to generate η³‐oxa‐allylnickel complexes ( 8 , 9 ), whereas insertion into dimeric η¹‐nickelenolate complex 6 ai did not take place. The diastereoselectivity of the insertion of an enone into 6 c having PCy₃ as a ligand differs from that into 6 ci having IPr as a ligand. In addition, the stereochemistry of η³‐oxa‐allylnickel complexes having IPr as a ligand is retained during reductive elimination to yield the corresponding [3+2] cycloaddition product, which is consistent with the diastereoselectivity observed in Ni⁰/IPr‐catalyzed [3+2] cycloaddition reactions of cyclopropyl ketones with enones. In contrast, reductive elimination from the η³‐oxa‐allylnickel having PCy₃ as a ligand proceeds with inversion of stereochemistry. This is probably due to rapid isomerization between syn and anti isomers prior to reductive elimination.     

An Improved Model for Predicting Proton NMR Shielding by Alkenes Based on Ab Initio GIAO Calculations

In a strong magnetic field, nuclei located over a carbon-carbon double bond experience NMR shielding effects that are the net result of the magnetic anisotropy of the nearby double bond and various other intramolecular shielding effects. We have used GIAO, a subroutine in Gaussian 4, to calculate isotropic shielding values and to predict the proton NMR shielding increment for a simple model system: methane held in various orientations and positions over ethene. The average proton NMR shielding increments of several orientations of methane have been plotted versus the Cartesian coordinates of the methane protons relative to the center of ethene. A single empirical equation for predicting the NMR shielding experienced by protons over a carbon-carbon double bond has been developed from these data. The predictive capability of this equation has been validated by comparing the shielding increments for several alkenes calculated using our equation to the experimentally observed shielding increments. This equation predicts the NMR shielding effects more accurately than a previous model that was based on only one orientation of methane over ethene. Deshielding is predicted by this equation for protons over the center and within about 3 Å of a carbon-carbon double bond. This result is in contrast to predictions made by the long-held shielding cone model based on the McConnell equation found in nearly every textbook on NMR, but is consistent with experimental observations.     

Magnetic Relaxation in Single‐Electron Single‐Ion Cerium(III) Magnets: Insights from Ab Initio Calculations

Detailed ab initio calculations were performed on two structurally different cerium(III) single‐molecule magnets (SMMs) to probe the origin of magnetic anisotropy and to understand the mechanism of magnetic relaxations. The complexes [Ce^III{Zn^II(L)}₂(MeOH)]BPh₄ ( 1 ) and [Li(dme)₃][Ce^III(cot′′)₂] ( 1 ; L=N,N,O,O‐tetradentate Schiff base ligand; 2 ; DME=dimethoxyethane, COT′′=1,4‐bis(trimethylsilyl)cyclooctatetraenyldianion), which are reported to be zero‐field and field‐induced SMMs with effective barrier heights of 21.2 and 30 K respectively, were chosen as examples. CASSCF+RASSI/SINGLE_ANISO calculations unequivocally suggest that m_J|±5/2〉 and |±1/2〉 are the ground states for complexes 1 and 2 , respectively. The origin of these differences is rooted back to the nature of the ligand field and the symmetry around the cerium(III) ions. Ab initio magnetisation blockade barriers constructed for complexes 1 and 2 expose a contrasting energy‐level pattern with significant quantum tunnelling of magnetisation between the ground state Kramers doublet in complex 2 . Calculations performed on several model complexes stress the need for a suitable ligand environment and high symmetry around the cerium(III) ions to obtain a large effective barrier.     

Molecular Weight Development in Free‐Radical Polymerization with Polyfunctional Chain‐Transfer Agents, 1. Equal Reactivity Model

The analytic expression for the weight‐average molecular weight development in free‐radical polymerization that involves a polyfunctional chain‐transfer agent is proposed. Free‐radical polymerization is kinetically controlled; therefore, the probability of chain connection with a polyfunctional chain‐transfer agent as well as the primary chain‐length distribution changes during the course of polymerization. We consider the primary chains formed at different times as different types of chains, and the heterochain branching model is used to obtain the weight‐average chain length at a given conversion level in a matrix formula, described as P_w = W { D _w + ( I + T ) SP ( I – TSP )^–1 D_f }. Because the primary chains are formed consecutively, the number of chain types N is extrapolated to infinity, but such extrapolation can be conducted with the calculated values for only three different N values. The criterion for the onset of gelation is simply described as a point at which the largest eigenvalue of the product of matrixes, TSP reaches unity, i. e., det  ( I – TSP ) = 0. The present model can readily be extended for the star‐shaped polyfunctional initiators, and the relationships between the model parameters and kinetic rate expression for such reaction systems are also shown.     

Solvent‐Dependent Structure of the I3− Ion Derived from Photoelectron Spectroscopy and Ab Initio Molecular Dynamics Simulations

Ab initio molecular dynamics (MD) simulations of the solvation of LiI₃ in four different solvents (water, methanol, ethanol, and acetonitrile) are employed to investigate the molecular and electronic structure of the I₃^? ion in relation to X‐ray photoelectron spectroscopy (XPS). Simulations show that hydrogen‐bond rearrangement in the solvation shell is coupled to intramolecular bond‐length asymmetry in the I₃^? ion. By a combination of charge analysis and I 4 d core‐level XPS measurements, the mechanism of the solvent‐induced distortions has been studied, and it has been concluded that charge localization mediates intermolecular interactions and intramolecular distortion. The approach involving a synergistic combination of theory and experiment probes the solvent‐dependent structure of the I₃^? ion, and the geometric structure has been correlated with the electronic structure.     

Slide Fastener Reduction of Graphene‐Oxide Edges by Calcium: Insight from Ab Initio Molecular Dynamics

The reduction of graphene oxide can be used as a simple way to produce graphene on a large scale. However, the numerous edges produced by the oxidation of graphite seriously degrade the quality of the graphene and its carrier transport property. In this work, the reduction of oxygen‐passivated graphene edges and the subsequent linking of separated graphene sheets by calcium are investigated by using first‐principles calculations. The calculations show that calcium can effectively remove the oxygen groups from two adjacent edges. The joining point of the edges serves as the starting point of the reduction and facilitates the reaction. Once the oxygen groups are removed, the crack is sutured. If the joining point is lacking, it becomes difficult to zip the separated fragments. A general electron‐reduction model and a random atom‐reduction model are suggested for these two situations. The present study sheds light on the reduction of graphene‐oxide edges by using reactive metals to give large‐sized graphene through a simple chemical reaction.     

Correlation of Quantitative Structure and Inhibition Phytotoxicity for Aromatic Compounds Using Ab Initio Method

1 INTRODUCTION Quantitative structure-activity relationship (QSAR) is one of the necessary methods to evaluate the ha- zards of organic chemicals. QSAR equation could be employed to forecast the biological activity of un- known compounds, which is significant for initial screening and evaluation of toxic compounds[1]. Aro- matic compounds are toxic organic compounds with relatively low water solubility, and their structure- activity relationship has been investigated with AM1 method[2]…     

Direct Ab Initio Dynamics Study on the Reaction of CH3CHF2 (HFC‐152a) with the Cl Atom

A direct ab initio dynamics method is used to investigate the hydrogen‐abstraction reaction CH₃CHF₂+Cl. One transition state is located for α‐H abstraction, and two are identified for β‐H abstraction. The potential‐energy surface (PES) is obtained at the G3(MP2)//MP2/6‐311G(d, p) level. Furthermore, the rate constants of the three channels are evaluated by using canonical variational transition‐state theory (CVT) with small‐curvature tunneling (SCT) contributions over a wide temperature range of 200–2500 K. The dynamic calculations show that the reaction proceeds mainly by α‐H abstraction over the whole temperature range. The calculated rate constants and branching ratios are both in good agreement with the available experimental values.     

A Reinvestigation of the Structure and Torsional Potential of N2O5 by Gas‐Phase Electron Diffraction Augmented by Ab Initio Theoretical Calculations

Gaseous N₂O₅ consists of two NO₂ groups bonded to a bridging O‐atom to form a nonlinear N−O−N moiety. The NO₂ groups undergo slightly hindered internal rotation around the bonds to the bridge so that instantaneous composition of the gaseous system is characterized by molecules with all combinations of torsion angles. In an earlier investigation, an attempt was made to determine the coefficients for an empirical form of the double‐rotor torsional potential, and the bond lengths and bond angles measured subject to assumptions that the structure of the O−NO₂ groups was invariant to torsion angle and that these groups had C_2v symmetry. The system has now been reinvestigated in terms of a more realistic model in which this symmetry restriction was relaxed, account was taken of structural changes in the NO₂ groups with torsion angle as predicted by ab initio theory at the B3LYP/6‐311+G* level, and a more convenient form of the torsional potential was assumed. The most stable conformation has C₂ symmetry with torsion angles τ₁ (defined as ∢(N−O−N=O₄)) equal to τ₂ (defined as ∢(N−O−N=O₆)) equal to 33.7°; because of the broad potential minimum in this region, the uncertainty in these angles is difficult to estimate, but is probably 3 – 4°. The results for the bond lengths and bond angles for the most stable conformation are r_g(N−O)=1.505(4) Å, r_g(N=O)=1.188(2) Å, ∢_α(N−O−N)=112.3(17)°, ∢_α(O=N=O)=134.2(4)°, 〈∢_α(O−N=O)〉=112.8(2)°. The difference between the symmetry‐nonequivalent O−N=O angles is estimated to be ca. 6.7° with the larger angle positioning the two N=O bonds on different NO₂ groups nearest each other. These average values are similar to those obtained in the original study. The main difference is found in the shape of the torsional potential, which at τ₁/τ₂=0/0 has a saddle point in the present work and a substantial peak in the earlier. The implication of the torsion‐angle findings for electron‐diffraction investigations of this type is discussed.     

Ab Initio MP2 and Density Functional Theory Computational Study of AcAlaNH2 Peptide Hydration: A Bottom‐Up Approach

AcAlaNH₂?n H₂O (n=1–13) complexes have been proposed as models to account for water solvent effects on the molecular properties of N‐acetyl‐L ‐alanine amide. Ab initio computations are planned to evaluate peptide–water interactions and to provide a means for approximating relative effects of the short‐range many‐body interactions arising in real solution without introducing any external parameters intended to quantify empirical or semiempirical potential‐energy functions. The present bottom‐up approach reveals the formation of compact ring clusters of water molecules strongly bonded to peptidic polar groups throughout hydrogen bonds. The explicit coordination of water molecules around the peptide renders the fully extended (FE) and polyproline II (PPII) conformers more stable with respect to the 3₁₀ helix or γ turn. The alternance of donor and acceptor groups on both sides of the FE and PPII conformers allows for synergy and extensive H‐bonding.