Dioxygen insertion into the gold(i)-hydride bond: spin orbit coupling effects in the spotlight for oxidative addition

Dioxygen Activation by Gold(I)-Distorted Porphyrin Dinuclear Complexes

Interactions between gold-based materials and dioxygen (O2) have motivated researchers to understand reaction mechanisms for O2 activation by homo- and heterogeneous gold catalysts. In this work, gold(I) porphyrin dinuclear complexes were synthesized with a saddle-distorted porphyrin ligand. The gold(I) porphyrin complexes showed unprecedented O2 activation in the presence of protic solvents to form gold(III) tetradentate porphyrin complexes. Mechanistic insights into the O2 activation by the gold(I) center were elucidated by spectroscopic measurements and theoretical calculations, revealing that dissociation of halides on the gold(I) center by alcohol solvents and hydrogen bonding of an N-H proton in the distorted porphyrin with dioxygen played important roles in establishing the unique reactivities of gold(I) complexes.

Diffusion Nuclear Magnetic Resonance Measurements on Cationic Gold (I) Complexes in Catalytic Conditions: Counterion and Solvent Effects

The amount of free ions, ion pairs, and higher aggregate of the possible species present in a solution during the gold(I)-catalyzed alkoxylation of unsaturated hydrocarbon, i.e., ISIP (inner sphere ion pair) [(NHC)AuX] and OSIP (outer sphere ion pairs) [(NHC)Au(TME)X] [NHC 1,3-bis(2,6-di-isopropylphenyl)-imidazol-2-ylidene; TME = tetramethylethylene (2,3-bis methyl-butene); X- = Cl-, BF4-, OTf-; and OTs- BArF4- (ArF = 3,5-(CF3)2C6H3)], has been determined. The 1H and 19F DOSY NMR measurements conducted in catalytic conditions indicate that the dissociation degree (α) of the equilibrium ion pair/free ions {[(NHC)Au(TME)X] [(NHC)Au(TME)]+ + X-} depends on the nature of the counterion (X-) when chloroform is the catalytic solvent: while the compounds containing OTs- and OTf- as the counterion gave a low α (which means a high number of ion pairs) of 0.13 and 0.24, respectively, the compounds containing BF4- and BArF4- showed higher α values of 0.36 and 0.32, respectively. These results experimentally confirm previous deductions based on catalytic and theoretical data: the lower the α value, the greater the catalytic activity because the anion that can activate methanol during a nucleophilic attack, although the lower propensity to activate methanol of BF4- and BArF4-, as suggested by the DFT calculations, cannot be completely overlooked. As for the effect of the solvent, α increases as the dielectric constant increases, as expected, and in particular, green solvents with high dielectric constants show a very high α (0.90, 0.84, 0.80, and 0.70 for propylene carbonate, γ-valerolactone, acetone, and methanol, respectively), thus confirming that the moderately high activity of NHC-Au-OTf in these solvents is due to the specific effect of polar functionalities (O-H, C=O, O-R) in activating methanol. Finally, the DOSY measurements conducted in p-Cymene show the formation of quadrupole species: under these conditions, the anion can better exercise its 'template' and 'activating' roles, giving the highest TOF.

Monomeric gold hydrides for carbon dioxide reduction: ligand effect on the reactivity

We analyzed the ligand electronic effect in the reaction between a [LAu(I)H]0/- hydride species and CO2, leading to a coordinated formate [LAu(HCOO)]0/-. We explored 20 different ligands, such as carbenes, phosphines and others, carefully selected to cover a wide range of electron-donor and -acceptor properties. We included in the study the only ligand, an NHC-coordinated diphosphene, that, thus far, experimentally demonstrated facile and reversible reaction between the monomeric gold(I) hydride and carbon dioxide. We elucidated the previously unknown reaction mechanism, which resulted to be concerted and common to all the ligands: the gold-hydrogen bond attacks the carbon atom of CO2 with one oxygen atom coordinating to the gold center. A correlation between the ligand σ donor ability, which affects the electron density at the reactive site, and the kinetic activation barriers of the reaction has been found. This systematic study offers useful guidelines for the rational design of new ligands for this reaction, while suggesting a few promising and experimentally accessible potential candidates for the stoichiometric or catalytic CO2 activation.

Metal-Free Activation of Molecular Oxygen by Quaternary Ammonium-Based Ionic Liquid: A Detail Mechanistic Study

Most oxidation processes in common organic synthesis and chemical biology require transition metal catalysts or metalloenzymes. Herein, we report a detailed mechanistic study of a metal-free oxygen (O2) activation protocol on benzylamine/alcohols using simple quaternary alkylammonium-based ionic liquids to produce products such as amide, aldehyde, imine, and in some cases, even aromatized products. NMR and various control experiments established the product formation and reaction mechanism, which involved the conversion of molecular oxygen into a hydroperoxyl radical via a proton-coupled electron transfer process. Detection of hydrogen peroxide in the reaction medium using colorimetric analysis supported the proposed mechanism of oxygen activation. Furthermore, first-principles calculations using density functional theory (DFT) revealed that reaction coordinates and transition state spin densities have a unique spin conversion of triplet oxygen leading to formation of singlet products via a minimum energy crossing point. In addition to DFT, domain-based local pair natural orbital coupled cluster, (DLPNO-CCSD(T)), and complete active space self-consistent field, CASSCF(20,14) methods complemented the above findings. Partial density of states analysis showed stabilization of π* orbital of oxygen in the presence of ionic liquid, making it susceptible to hydrogen abstraction in a mild, metal-free condition. Inductively coupled plasma atomic emission spectroscopic (ICP-AES) analysis of reactant and ionic liquids clearly showed the absence of any significant transition metal contamination. The current results described the origin of O2 activation within the context of molecular orbital (MO) theory and opened up a new avenue for the use of ionic liquids as inexpensive, multifunctional and high-performance alternative to metal-based catalysts for O2 activation.

Insertion of Molecular Oxygen into a Gold(III)-Hydride Bond

The use of molecular oxygen as an oxidant in chemical synthesis has significant environmental and economic benefits, and it is widely used as such in large-scale industrial processes. However, its adoption in highly selective homogeneous catalytic transformations, particularly to produce oxygenated organics, has been hindered by our limited understanding of the mechanisms by which O2 reacts with transition metals. Of particular relevance are the mechanisms of the reactions of oxygen with late transition metal hydrides as these metal centers are better poised to release oxygenated products. Homogeneous catalysis with gold complexes has markedly increased, and herein we report the synthesis and full characterization of a rare AuIII-H, supported by a diphosphine pincer ligand (tBuPCP = 2,6-bis(di-tert-butylphosphinomethyl)benzene). [(tBuPCP)AuIII-H]+ was found to cleanly react with molecular oxygen to yield a stable AuIII-OOH complex that was also fully characterized. Extensive kinetic studies on the reaction via variable temperature NMR spectroscopy have been completed, and the results are consistent with an autoaccelerating radical chain mechanism. The observed kinetic behavior exhibits similarities to that of previously reported PdII-H and PtIV-H reactions with O2 but is not fully consistent with any known O2 insertion mechanism. As such, this study contributes to the nascent fundamental understanding of the mechanisms of aerobic oxidation of late metal hydrides.

Experimental and Theoretical Investigation of Ion Pairing in Gold(III) Catalysts

The ion pairing structure of the possible species present in solution during the gold(III)-catalyzed hydration of alkynes: [(ppy)Au(NHC)Y]X2 and [(ppy)Au(NHC)X]X [ppy = 2-phenylpyridine, NHC = NHCiPr = 1,3-bis(2,6-di-isopropylphenyl)-imidazol-2-ylidene; NHC = NHCmes = 1,3-bis(2,4,6-trimethylphenyl)-imidazol-2-ylidene X = Cl-, BF4-, OTf-; Y = H2O and 3-hexyne] are determined. The nuclear overhauser effect nuclear magnetic resonance (NMR) experimental measurements integrated with a theoretical description of the system (full optimization of different ion pairs and calculation of the Coulomb potential surface) indicate that the preferential position of the counterion is tunable through the choice of the ancillary ligands (NHCiPr, NHCmes, ppy, and Y) in [(ppy)Au(NHC)(3-hexyne)]X2 activated complexes that undergo nucleophilic attack. The counterion can approach near NHC, pyridine ring of ppy, and gold atom. From these positions, the anion can act as a template, holding water in the right position for the outer-sphere attack, as observed in gold(I) catalysts.

A general method for locating stationary points on the mixed-spin surface of spin-forbidden reaction with multiple spin states

Some chemical reactions proceed on multiple potential energy surfaces and are often accompanied by a change in spin multiplicity, being called spin-forbidden reactions, where the spin-orbit coupling (SOC) effects play a crucial role. In order to efficiently investigate spin-forbidden reactions with two spin states, Yang et al. [Phys. Chem. Chem. Phys. 20, 4129-4136 (2018)] proposed a two-state spin-mixing (TSSM) model, where the SOC effects between the two spin states are simulated by a geometry-independent constant. Inspired by the TSSM model, we suggest a multiple-state spin-mixing (MSSM) model in this paper for the general case with any number of spin states, and its analytic first and second derivatives have been developed for locating stationary points on the mixed-spin potential energy surface and estimating thermochemical energies. To demonstrate the performance of the MSSM model, some spin-forbidden reactions involving 5d transition elements are calculated using the density functional theory (DFT), and the results are compared with the two-component relativistic ones. It is found that MSSM DFT and two-component DFT calculations may provide very similar stationary-point information on the lowest mixed-spin/spinor energy surface, including structures, vibrational frequencies, and zero-point energies. For the reactions containing saturated 5d elements, the reaction energies by MSSM DFT and two-component DFT agree very well within 3 kcal/mol. As for the two reactions OsO+ + CH4 → OOs(CH2)+ + H2 and W + CH4 → WCH2 + H2 involving unsaturated 5d elements, MSSM DFT may also yield good reaction energies of similar accuracy but with some counterexamples. Nevertheless, the energies may be remarkably improved by a posteriori single point energy calculations using two-component DFT at the MSSM DFT optimized geometries, and the maximum error of about 1 kcal/mol is almost independent of the SOC constant used. The MSSM method as well as the developed computer program provides an effective utility for studying spin-forbidden reactions.

Reactivity of NHC/diphosphene-coordinated Au(I)-hydride

We report the reactivity of isolable Au(i)-hydride stabilized by an NHC-coordinated diphosphene towards substrates containing C-C and N-N multiple bonds (NHC = N-heterocyclcic carbene). Reactions with dimethyl acetylenedicarboxylate and azobenzene lead to a trans-addition of the Au(i)-H across the C-C triple bond and the N-N double bond, respectively. In contrast, the reaction with ethyl diazoacetate affords a gold(i)-hydrazonide as the 1,1-addition product to the terminal nitrogen atom. With phenyl acetylene, the corresponding Au(i)-alkynyl complex is obtained under the elimination of dihydrogen. Strikingly, diphosphene-containing Au(i)-hydride is more reactive - affording different products in some cases - than a related NHC-stabilized Au(i)-hydride without the mediating diphosphene moiety.

Constructing Spin-Adiabatic States for the Modeling of Spin-Crossing Reactions. I. A Shared-Orbital Implementation

In the modeling of spin-crossing reactions, it has become popular to directly explore the spin-adiabatic surfaces. Specifically, through constructing spin-adiabatic states from a two-state Hamiltonian (with spin-orbit coupling matrix elements) at each geometry, one can readily employ advanced geometry optimization algorithms to acquire a "transition state" structure, where the spin crossing occurs. In this work, we report the implementation of a fully-variational spin-adiabatic approach based on Kohn-Sham density functional theory spin states (sharing the same set of molecular orbitals) and the Breit-Pauli one-electron spin-orbit operator. For three model spin-crossing reactions [predissociation of N₂O, singlet-triplet conversion in CH₂, and CO addition to Fe(CO)₄], the spin-crossing points were obtained. Our results also indicated the Breit-Pauli one-electron spin-orbit coupling can vary significantly along the reaction pathway on the spin-adiabatic energy surface. On the other hand, due to the restriction that low-spin and high-spin states share the same set of molecular orbitals, the acquired spin-adiabatic energy surface shows a cusp (i.e. a first-order discontinuity) at the crossing point, which prevents the use of standard geometry optimization algorithms to pinpoint the crossing point. An extension with this restriction removed is being developed to achieve the smoothness of spin-adiabatic surfaces.

Spin-Forbidden Reactivity of Transition Metal Oxo Species: Exploring the Potential Energy Surfaces

Spin-forbidden reactions are frequently encountered when transition metal oxo species are involved, particularly in oxygen transfer reactivity. The computational study of such reactions is challenging, because reactants and products are located on different spin potential energy surfaces (PESs). One possible approach to describe these reactions is the so-called minimum energy crossing point (MECP) between the diabatic reactants and products PESs. Alternatively, inclusion of spin-orbit coupling (SOC) effects allows to locate a saddle point on a single adiabatic PES (TS SOC). The TS SOC approach is rarely applied because of its high computational cost. Recently evidence for a TS SOC impact on significantly lowering the activation barrier in dioxygen addition to a carbene-gold(I)-hydride complex reaction (Chem. Sci. 2016, 7, 7034-7039) or even on predicting a qualitatively different reaction mechanism in mercury methylation by cobalt corrinoid (Angew. Chem. Int. Ed. 2016, 55, 11503-11506) has been put forward. Using MECP and TS SOC approaches a systematic analysis is provided here of three prototypical transition metal oxo spin-forbidden processes to investigate their implications on reactivity. Cycloaddition of ethylene to chromyl chloride (CrO₂ Cl₂ +C₂ H₄ ), iron oxide cation insertion into the hydrogen molecule (FeO⁺ +H₂ ) and H-abstraction from toluene by a Mn^V -oxo-porphyrin cation (MnOP(H₂ O)⁺ +C₆ H₅ CH₃ ) are case studies. For all these processes the MECP and TS SOC results are compared, which show that the spin-forbidden reactivity of transition metal oxo species can be safely described by a MECP approach, at least for the first-row transition metals investigated here, where the spin-orbit coupling is relatively weak. However, for the Mn-oxo reactivity, the MECP and TS SOC have been found to be crucial for a correct description of the reaction mechanism. In particular, the TS SOC approach allows to straightforwardly explore detailed features of the adiabatic potential energy surface which in principle could affect the overall reaction rate in cases where the involved diabatic PESs are tricky.

Experimental and Computational Investigation of the Aerobic Oxidation of a Late Transition Metal-Hydride

The rational development of homogeneous catalytic systems for selective aerobic oxidations of organics has been hampered by the limited available knowledge of how oxygen reacts with important organometallic intermediates. Recently, several mechanisms for oxygen insertion into late transition metal-hydride bonds have been described. Contributing to this nascent understanding of how oxygen reacts with metal-hydrides, a detailed mechanistic study of the reaction of oxygen with the Ir^III hydride complex (^dmPhebox)Ir(OAc)(H) (1) in the presence of acetic acid, which proceeds to form the Ir^III complex (^dmPhebox)Ir(OAc)₂(OH₂) (2), is described. The evidence supports a multifaceted mechanism wherein a small amount of an initially formed metal hydroperoxide proceeds to generate a metal-oxyl species that then initiates a radical chain reaction to rapidly convert the remaining Ir^III-H. Insight into the initiation step was gained through kinetic and mechanistic studies of the radical chain inhibition by BHT (butylated hydroxytoluene). Computational studies were employed to contribute to a further understanding of initiation and propagation in this system.

Alkyne Activation with Gold(III) Complexes: A Quantitative Assessment of the Ligand Effect by Charge-Displacement Analysis

A quantitative assessment of the Dewar-Chatt-Duncanson components of the Au(III)-alkyne bond in a series of cationic and dicationic bis- and monocyclometalated gold(III) complexes with 2-butyne via charge-displacement (CD) analysis is reported. Bonding between Au(III) and 2-butyne invariably shows a dominant σ donation component, a smaller, but significant, π back-donation, and a remarkable polarization of the alkyne CC triple bond toward the metal fragment. A very large net electron charge transfer from CC triple bond to the metal fragment results, which turns out to be unexpectedly insensitive to the charge of the complex and more strictly related to the nature of the ancillary ligand. The combination of σ donation, π back-donation, and polarization effects is in fact modulated by the different ligand frameworks, with ligands bearing atoms different from carbon in trans position with respect to the alkyne emerging as especially interesting for both imparting Au(III)-alkyne bond stability and inducing a more effective alkyne activation. A first attempt to figure out a rationale on the bonding/reactivity relationship for Au(III)-alkyne is made by performing a comparative study in a model nucleophilic attack of water to the alkyne triple bond. Smaller π back-donation facilitates alkyne slippage in the transition states, which is energetically less demanding for Au(III) than for Au(I), and suggests a greater propensity of Au(III) to facilitate the nucleophilic attack.

Hydride-induced ligand dynamic and structural transformation of gold nanoclusters during a catalytic reaction

Quasi-homogeneous ligand-protected gold nanoclusters (Au NCs) with atomic precision and well-defined structure offer great opportunity for exploring the catalytic nature of nanogold catalysts at a molecular level. Herein, using real-time electrospray ionization mass spectrometry (ESI-MS), we have successfully identified the desorption and re-adsorption of p-mercaptobenzoic acid (p-MBA) ligands from Au25(p-MBA)18 NC catalysts during the hydrogenation of 4-nitrophenol in solution. This ligand dynamic (desorption and re-adsorption) would initiate structural transformation of Au25(p-MBA)18 NC catalysts during the reaction, forming a mixture of smaller Au NCs (Au23(p-MBA)16 as the major species) at the beginning of catalytic reaction, which could further be transformed into larger Au NCs (Au26(p-MBA)19 as the major species). The adsorption of hydrides (from NaBH4) is identified as the determining factor that could induce the ligand dynamic and structural transformation of NC catalysts. This study provides fundamental insights into the catalytic nature of Au NCs, including catalytic mechanism, active species and stability of Au NC catalysts during a catalytic reaction.

Ligand Effect on Bonding in Gold(III) Carbonyl Complexes

We quantitatively assess the Dewar-Chatt-Duncanson (DCD) components of the Au(III)-CO bond and the charge density polarization at the CO, in a series of neutral, cationic, and dicationic bis- and monocyclometalated gold(III) complexes via charge-displacement (CD) analysis. A striking feature concerns the very small net electron charge flux from CO to the metal fragment which is unexpectedly stable toward both the charge of the complex and the oxidation state of gold (I, III). All systems exhibit a similar trend for the σ charge rearrangement in the region of the carbonyl bond, where, by contrast, the π back-donation trend variation is large, which is strictly correlated to the change in CO bond distance and the shift in CO stretching frequencies, in close analogy with the gold(I) carbonyl complexes. In the whole series of gold(III) compounds, a large Au(III) ← CO σ donation is measured (from 0.19 to 0.31 electrons), as well as a significant Au(III) → CO π back-donation (from -0.09 to -0.22 electrons), which however is not generally able to completely balance the polarization of the CO π electrons in the direction from oxygen to carbon (C ← O) induced by the presence of the metal fragment [LAu(III)]^0/+1/+2. Surprisingly, all the gold(III) complexes in the series are characterized by a very small anisotropy in the Au(III) → CO in-plane and out-of-plane π back-donation components, in sharp contrast with the marked anisotropy found before for the experimentally characterized [(C^N^C)Au(III)CO]⁺ complex. A first attempt to figure out a rationale on the bonding/reactivity relationship for Au(III)-CO is made by performing a comparative study with an isostructural [(N^N^C)Pt(II)CO]⁺ complex in a model water-gas shift (WGS) reaction.

Automated reaction path searches for spin-forbidden reactions

Many catalytic and biomolecular reactions containing transition metals involve changes in the electronic spin state. These processes are referred to as "spin-forbidden" reactions within nonrelativistic quantum mechanics framework. To understand detailed reaction mechanisms of spin-forbidden reactions, one must characterize reaction pathways on potential energy surfaces with different spin states and then identify crossing points. Here we propose a practical computational scheme, where only the lowest mixed-spin eigenstate obtained from the diagonalization of the spin-coupled Hamiltonian matrix is used in reaction path search calculations. We applied this method to the ^6,4 FeO⁺  + H₂ → ^6,4 Fe⁺  + H₂ O, ^6,4 FeO⁺  + CH₄ → ^6,4 Fe⁺  + CH₃ OH, and ⁷ Mn⁺  + OCS → ⁵ MnS⁺  + CO reactions, for which crossings between the different spin states are known to play essential roles in the overall reaction kinetics. © 2018 Wiley Periodicals, Inc.

Transition states of spin-forbidden reactions

New approximation method for locating stationary points on lowest spin-coupled potential energy surface (PES) using density functional calculations.

Mechanistic insights into the light-driven hydrogen evolution reaction from formic acid mediated by an iridium photocatalyst

Electronic structure calculations shed important mechanistic light on light-driven hydrogen evolution from formic acid mediated by an iridium photocatalyst.

The gold(iii)–CO bond: a missing piece in the gold carbonyl complex landscape

Surprisingly, charge-displacement (CD) analysis reveals a significant π back-donation in the gold(iii)–CO bond.

The ligand effect on the oxidative addition of dioxygen to gold(i)–hydride complexes

The activation energy barriers of the O₂ to [LAuH] oxidative addition, calculated by including spin–orbit coupling (SOC) effects, quantitatively correlate with the σ donation component of the L–AuH bond.