Ligand Effect on Bonding in Gold(III) Carbonyl Complexes

Quantification of the Donor-Acceptor Character of Ligands by the Effective Fragment Orbitals

Metal-ligand interactions are at the heart of transition metal complexes. The Dewar-Chat-Duncanson model is often invoked, whereby the metal-ligand bonding is decomposed into the simultaneous ligand→metal electron donation and the metal→ligand back-donation. The separate quantification of both effects is not a trivial task, neither from experimental nor computational exercises. In this work we present the effective fragment orbitals (EFOs) and their occupations as a general procedure beyond the Kohn-Sham density functional theory (KS-DFT) framework for the identification and quantification of donor-acceptor interactions, using solely the wavefunction of the complex. Using a common Fe(II) octahedral complex framework, we quantify the σ-donor, π-donor, and π-acceptor character for a large and chemically diverse set of ligands, by introducing the respective descriptors σd, πd, and πa. We also explore the effect of the metal size, coordination number, and spin state on the donor/acceptor features. The spin-state is shown to be the most critical effect, inducing a systematic decrease of the sigma donation and π-backdonation going from low spin to high spin. Finally, we illustrate the ability of the EFOs to rationalize the Tolman electronic parameter and the trans influence in planar square complexes in terms of these new descriptors.

Infrared Photodissociation Spectroscopy of Mass-Selected Dinuclear Transition Metal Boride Carbonyl Cluster Cations

The transition-metal-boron bonding interactions and geometric structures of heterodinuclear transition metal carbonyl cluster cations BM(CO)n+ (M = Co, Ni, and Cu) are studied by a combination of the infrared photodissociation spectroscopy and density functional theory calculations at the B3LYP/def2-TZVP level. The BCu(CO)5+ and BCo(CO)6+ cations are characterized as an (CO)2B-M(CO)3/4+ structure involving an σ-type (OC)2B → M(CO)3,4+ dative bonding with end-on carbonyls, while for BNi(CO)5,6+ complexes with a bridged carbonyl, a 3c-2e bond involving the 5σ electrons of the bridged carbonyl and an electron-sharing bond between the B(CO)2 fragment and the Ni(CO)2,3+ subunits were revealed. Moreover, the fundamental driving force of the exclusive existence of a bridged carbonyl group in the boron-nickel complexes has been demonstrated to stem from the desire of the B and Ni centers for the favorable 8- and 18-electron structures.

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.

Bond-Dissociation Energies to Probe Pyridine Electronic Effects on Organogold(III) Complexes: From Methodological Developments to Application in π-Backdonation Investigation and Catalysis

In this work, we report on the synthesis of several organogold(III) complexes based on 4,4'-diterbutylbiphenyl (C^C) and 2,6-bis(4-terbutylphenyl)pyridine (C^N^C) ligands and bond with variously substituted pyridine ligands (pyrR). Altogether, 33 complexes have been prepared and studied with mass spectrometry using higher-energy collision dissociation (HCD) in an Orbitrap mass spectrometer. A complete methodology including the kinetic modeling of the dissociation process based on the Rice-Ramsperger-Kassel-Marcus (RRKM) statistical method is proposed to obtain critical energies E0 of the pyrR loss for all complexes. The capacity of these E0 values to describe the pyridine ligand effect is further explored, at the same time as more classical descriptors such as 1H pyridinic NMR shift variation upon coordination and Au-NpyrR bond length measured by X-ray diffraction. An extensive theoretical work, including density functional theory (DFT) and domain-based local pair natural orbital coupled-cluster theory (DLPNO-CCSD(T)) methods, is also carried out to provide bond-dissociation energies, which are compared to experimental results. Results show that dissociation energy outperforms other descriptors, in particular to describe ligand effects over a large electronic effect range as seen by confronting the results to the pyrR pKa values. Further insights into the Au-NpyrR bond are obtained through an energy decomposition analysis (EDA) study, which confirms the isolobal character of Au+ with H+. Finally, the correlation between the lability of the pyridine ligands toward the catalytic efficiency of the complexes could be demonstrated in an intramolecular hydroarylation reaction of alkyne. The results were rationalized considering both pre-catalyst activation and catalyst reactivity. This study establishes the possibility of correlating dissociation energy, which is a gas-phase descriptor, with condensed-phase parameters such as catalysis efficiency. It therefore holds great potential for inorganic and organometallic chemistry by opening a convenient and easy way to evaluate the electronic influence of a ligand toward a metallic center.

Vibrational Properties of CO Adsorbed on Au Single Atom Catalysts on TiO2(101), ZrO2(101), CeO2(111), and LaFeO3(001) Surfaces: A DFT Study

The nature and local environment of Au single atoms supported and stabilized on four different oxides is studied by means of DFT + U calculations using CO as probe molecule and its stretching frequency, ωe, as a fingerprint of the site where the Au atom is bound. Four oxides are considered, anatase TiO2, tetragonal ZrO2, cubic CeO2, and a perovskite LaFeO3. In this latter case a recently reported experimental study has detected a stretching mode for CO adsorbed on Au1/LaFeO3 of 2215 cm−1, with a large blue shift, ∆ω(CO) = 72 cm−1 with respect to free CO. In order to identify the Au adsorption site that can give rise to this large blue-shift we have considered five cases: (a) Au replacing a lattice cation, (Au)subM; (b) Au replacing a lattice O anion, (Au)subO; (c) Au adsorbed on the surface, (Au)ads; (d) Au bound to an extra O atom on the surface, (AuO)ads, or (e) Au bound to two extra O atoms on the surface, (AuO2)ads. It turns out that the correct reproduction of ∆ω for CO adsorbed on positively charged gold, Auδ+, is challenging for DFT. Therefore, we have performed a comparative study of Auδ+-CO molecular compounds for which ωe(CO) is known experimentally using various kinds of DFT functionals and accurate CCSD and CCSD(T) quantum chemistry methods. Also based on this comparison we propose a tentative assignment for the observed frequency of CO adsorbed on Au1/LaFeO3 single atom catalyst.

Graphic Abstract

Experimental and theoretical studies on gold(III) carbonyl complexes: reductive C,H- and C,C bond formation

The reactivity of cationic (C^C)gold(iii) carbonyl complexes was investigated. While the in situ-formed IPrAu(bph)CO+ complex (bph = biphenyl-2,2'-diyl) does not undergo a migratory insertion of CO into the neighboring gold-carbon bond, nucleophiles can attack the coordinated CO moiety intermolecularly. Water as a nucleophile initiates a CO2 extrusion combined with a reductive C,H bond formation. The rapid formation of a gold(i) species from an intermediary gold(iii) carbonyl has not been observed before and shows a significant difference in reactivity between (C^C) and (C^N^C)gold(iii) carbonyls. The latter have been reported to form stable gold(iii) hydrides via the WGS reaction. In the case of methanol acting as a nucleophile attacking the gold(iii) carbonyl, no extrusion of CO2 is observed. Instead an intermediary gold(iii) carboxyl complex forms an aryl carboxylate via reductive C-C bond elimination. Experimental and theoretical studies on the mechanism explain the observed selectivities and give new insights into the reactivity of elusive gold(iii) carbonyls.

Monitoring of the Pre-Equilibrium Step in the Alkyne Hydration Reaction Catalyzed by Au(III) Complexes: A Computational Study Based on Experimental Evidences

The coordination ability of the [(ppy)Au(IPr)]²⁺ fragment [ppy = 2-phenylpyridine, IPr = 1,3-bis(2,6-di-isopropylphenyl)-imidazol-2-ylidene] towards different anionic and neutral X ligands (X = Cl⁻, BF₄⁻, OTf⁻, H₂O, 2-butyne, 3-hexyne) commonly involved in the crucial pre-equilibrium step of the alkyne hydration reaction is computationally investigated to shed light on unexpected experimental observations on its catalytic activity. Experiment reveals that BF₄⁻ and OTf⁻ have very similar coordination ability towards [(ppy)Au(IPr)]²⁺ and slightly less than water, whereas the alkyne complex could not be observed in solution at least at the NMR sensitivity. Due to the steric hindrance/dispersion interaction balance between X and IPr, the [(ppy)Au(IPr)]²⁺ fragment is computationally found to be much less selective than a model [(ppy)Au(NHC)]²⁺ (NHC = 1,3-dimethylimidazol-2-ylidene) fragment towards the different ligands, in particular OTf⁻ and BF₄⁻, in agreement with experiment. Effect of the ancillary ligand substitution demonstrates that the coordination ability of Au(III) is quantitatively strongly affected by the nature of the ligands (even more than the net charge of the complex) and that all the investigated gold fragments coordinate to alkynes more strongly than H₂O. Remarkably, a stabilization of the water-coordinating species with respect to the alkyne-coordinating one can only be achieved within a microsolvation model, which reconciles theory with experiment. All the results reported here suggest that both the Au(III) fragment coordination ability and its proper computational modelling in the experimental conditions are fundamental issues for the design of efficient catalysts.

Tuning the hyperconjugative aromaticity in Au(III)-substituted indoliums

As a fundamental concept in chemistry, aromaticity has been extended from traditional organics to organometallics. Similarly, hyperconjugative aromaticity (HCA) has also been developed from main group to transition metal systems through the hyperconjugation of the substituents. However, it remains unclear that how the oxidation state of transition metal in the substituents affects the HCA. Herein, we demonstrate via density functional theory calculations that HCA could disappear in indoliums when the Au(i) substituents are changed to the Au(iii) ones. By tuning the ligand or cis-trans isomerization, HCA could be regained or enhanced in indoliums containing Au(iii) substitutents.

Spectroscopic/Bond Property Relationship in Group 11 Dihydrides via Relativistic Four-Component Methods

Group 11 dihydrides MH2- (M = Cu, Ag, Au, Rg) have been much less studied than the corresponding MH compounds, despite having potentially several interesting applications in chemical research. In this work, their main spectroscopic constants (bond lengths, dissociation energies, and force constants) have been evaluated by means of highly accurate relativistic four-component coupled cluster (4c-CCSD(T)) calculations in combination with large basis sets. Periodic trends have been quantitatively explained by the charge-displacement/natural orbitals for chemical valence (CD-NOCV) analysis based on the four-component relativistic Dirac-Kohn-Sham method, which allows a consistent picture of the nature of the M-H bond to be obtained on going down the periodic table in terms of Dewar-Chatt-Duncanson bonding components. A strong ligand-to-metal donation drives the M-H bond and it is responsible for the heterolytic (HM···H-) dissociation energies to increase monotonically from Cu to Rg, with RgH2- showing the strongest and most covalent M-H bond. The "V"-shaped trend observed for the bond lengths, dissociation energies, and stretching frequencies can be explained in terms of relativistic effects and, in particular, of the relativistically enhanced sd hybridization occurring at the metal, which affects the metal-ligand distances in heavy transition-metal complexes. The sd hybridization is very small for Cu and Ag, whereas it becomes increasingly important for Au and Rg, being responsible for the increasing covalent character of the bond, the sizable contraction of the Au-H and Rg-H bonds, and the observed trend. This work rationalizes the spectroscopic/bond property relationship in group 11 dihydrides within highly accurate relativistic quantum chemistry methods, paving the way for their applications in chemical bond investigations involving heavy and superheavy elements.

Spin-orbit coupling is the key to unraveling intriguing features of the halogen bond involving astatine

The effect of spin-orbit coupling (SOC) on the halogen bond involving astatine has been investigated using state-of-the-art two- and four-component relativistic calculations. Adducts between Cl-X (X = Cl, Br, I and At) and ammonia have been selected to establish a trend on going down the periodic table. The SOC influence has been explored not only on the geometric and energetic features that can be used to characterize the halogen bond strength but also on the three main contributions to it that are the charge transfer, the "σ-hole" (i.e. the localized region with a net positive electrostatic potential at the halogen site) and the "polar flattening" (which is related to the effective shape of the halogen site). A surprisingly large increase of the Cl-At dipole moment, due to the inclusion of SOC, has been worked out using four-component CCSD(T) reference calculations, indicating that this bond is significantly more ionic than one may predict. Due to the SOC effect, which induces a peculiar charge accumulation on the At side in the Cl-At dimer, a weakening of the astatine-mediated halogen bond occurs arising from the (i) reduced amount of charge transfer, (ii) decrease of the polar flattening and (iii) lowering of the short-range Coulomb potential. The analysis of the electronic structure of the Cl-At moiety allows for a rationalization of the SOC effects on all the considered features of the halogen bond, including an unprecedented unsymmetrical charge back-donation from Cl-At to ammonia.

Charge Displacement Analysis-A Tool to Theoretically Characterize the Charge Transfer Contribution of Halogen Bonds

Theoretical bonding analysis is of prime importance for the deep understanding of the various chemical interactions, covalent or not. Among the various methods that have been developed in the last decades, the analysis of the Charge Displacement function (CD) demonstrated to be useful to reveal the charge transfer effects in many contexts, from weak hydrogen bonds, to the characterization of σ hole interactions, as halogen, chalcogen and pnictogen bonding or even in the decomposition of the metal-ligand bond. Quite often, the CD analysis has also been coupled with experimental techniques, in order to give a complete description of the system under study. In this review, we focus on the use of CD analysis on halogen bonded systems, describing the most relevant literature examples about gas phase and condensed phase systems. Chemical insights will be drawn about the nature of halogen bond, its cooperativity and its influence on metal-ligand bond components.

Catalytic Au(i)/Au(iii) arylation with the hemilabile MeDalphos ligand: unusual selectivity for electron-rich iodoarenes and efficient application to indoles

The ability of the hemilabile (P,N) MeDalphos ligand to trigger oxidative addition of iodoarenes to gold has been thoroughly studied. Competition experiments and Hammett correlations substantiate a clear preference of gold for electron-enriched substrates both in stoichiometric oxidative addition reactions and in catalytic C-C cross-coupling with 1,3,5-trimethoxybenzene. This feature markedly contrasts with the higher reactivity of electron-deprived substrates typically encountered with palladium. Based on DFT calculations and detailed analysis of the key transition states (using NBO, CDA and ETS-NOCV methods in particular), the different behavior of the two metals is proposed to result from inverse electron flow between the substrate and metal. Indeed, oxidative addition of iodobenzene is associated with a charge transfer from the substrate to the metal at the transition state for gold, but opposite for palladium. The higher electrophilicity of the gold center favors electron-rich substrates while important back-donation from palladium favors electron-poor substrates. Facile oxidative addition of iodoarenes combined with the propensity of gold(iii) complexes to readily react with electron-rich (hetero)arenes prompted us to apply the (MeDalphos)AuCl complex in the catalytic arylation of indoles, a challenging but very important transformation. The gold complex proved to be very efficient, general and robust. It displays complete regioselectivity for C3 arylation, it tolerates a variety of functional groups at both the iodoarene and indole partners (NO2, CO2Me, Br, OTf, Bpin, OMe…) and it proceeds under mild conditions (75 °C, 2 h).

Cationic Gold(I) Diarylallenylidene Complexes: Bonding Features and Ligand Effects

Using computational approaches, we qualitatively and quantitatively assess the bonding components of a series of experimentally characterized Au(I) diarylallenylidene complexes (N.Kim, R.A.Widenhoefer, Angew. Chem. Int. Ed. 2018, 57, 4722-4726). Our results clearly demonstrate that Au(I) engages only weakly in π-backbonding, which is, however, a tunable bonding component. Computationally identified trends in bonding are clearly correlated with the substitution patterns of the aryl substituents in the Au(I) diarylallenylidene complexes and good agreement is found with the previously reported experimental data, such as IR spectra, 13 C NMR chemical shifts and rates of decomposition together with their corresponding barrier heights, further substantiating the computational findings. The description of the bonding patterns in these complexes allow predictions of their spectroscopic features, their reactivity and stability.

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.

Ubiquity of cis-Halide → Isocyanide Direct Interligand Interaction in Organometallic Complexes

We recently reported a density functional theory (DFT) analysis of the Nb(V)-C bond in various NbCl₅(L) complexes, discovering that the carbon ligand L receives electronic density from the metal (classical back-donation) and from the chlorides in the cis position (direct interligand interaction). Here we report the synthesis and the structural characterization of two new coordination compounds of niobium pentahalides, i.e., NbX₅(CNXyl) (X = Cl, Br; Xyl = 2,6-C₆H₃Me₂), and the corresponding DFT analyses of the Nb(V)-C bond using the Natural Orbitals for Chemical Valence-Charge Displacement (NOCV-CD) approach, confirming the presence of a cis-halide → isocyanide direct interligand interaction. To verify whether the latter is limited to Nb complexes or not, we performed a NOCV-CD analysis on a series of several organometallic complexes based on Ti(IV), Nb(V), Ta(V), Rh(III), Pd(II), and Au(III), all of which bear one halide ligand and m-xylyl-isocyanide in a mutual cis position, revealing that the cis-halide → isocyanide interaction is always present.

Lewis Base Activation of Lewis Acid: A Detailed Bond Analysis

The effect of a Lewis base (LB) on the nucleophilic attack on chalcogeniranium (chalcogen = sulfur and selenium) cations, the so-called LB activation of a Lewis acid, has been studied coupling natural orbital for chemical valence decomposition of the orbital interaction energy with charge displacement analysis. This methodology provides a detailed and accurate description of all the interactions (LB···chalcogen, chalcogen···olefin and olefin···ammonia) present in the system and leads to a deeper understanding of how they influence each other at all stages of the reaction: reactant complex, transition state, and product complex. In particular, the bond between the chalcogen and the olefin has been decomposed in terms of σ donation/π back-donation and the bond components quantified. This allowed determination of a linear relationship between the activation barrier of the nucleophilic attack and the net amount of charge donated by the olefin to the chalcogen.