Screening substituent and backbone effects on the properties of bidentate P,P-donor ligands (LKB-PP(screen))

Bulky, electron-rich, renewable: analogues of Beller's phosphine for cross-couplings

In recent years, considerable progress has been made in the conversion of biomass into renewable chemicals, yet the range of value-added products that can be formed from biomass remains relatively small. Herein, we demonstrate that molecules available from biomass serve as viable starting materials for the synthesis of phosphine ligands, which can be used in homogeneous catalysis. Specifically, we prepared renewable analogues of Beller's ligand (di(1-adamantyl)-n-butylphosphine, cataCXium® A), which is widely used in homogeneous catalysis. Our new renewable phosphine ligands facilitate Pd-catalysed Suzuki-Miyaura, Stille, and Buchwald-Hartwig coupling reactions with high yields, and our catalytic results can be rationalized based on the stereoelectronic properties of the ligands. The new phosphine ligands generate catalytic systems that can be applied for the late-stage functionalization of commercial drugs.

Building a Toolbox for the Analysis and Prediction of Ligand and Catalyst Effects in Organometallic Catalysis

Computers have become closely involved with most aspects of modern life, and these developments are tracked in the chemical sciences. Recent years have seen the integration of computing across chemical research, made possible by investment in equipment, software development, improved networking between researchers, and rapid growth in the application of predictive approaches to chemistry, but also a change of attitude rooted in the successes of computational chemistry-it is now entirely possible to complete research projects where computation and synthesis are cooperative and integrated, and work in synergy to achieve better insights and improved results. It remains our ambition to put computational prediction before experiment, and we have been working toward developing the key ingredients and workflows to achieve this.The ability to precisely tune selectivity along with high catalyst activity make organometallic catalysts using transition metal (TM) centers ideal for high-value-added transformations, and this can make them appealing for industrial applications. However, mechanistic variations of TM-catalyzed reactions across the vast chemical space of different catalysts and substrates are not fully explored, and such an exploration is not feasible with current resources. This can lead to complete synthetic failures when new substrates are used, but more commonly we see outcomes that require further optimization, such as incomplete conversion, insufficient selectivity, or the appearance of unwanted side products. These processes consume time and resources, but the insights and data generated are usually not tied to a broader predictive workflow where experiments test hypotheses quantitatively, reducing their impact.These failures suggest at least a partial deviation of the reaction pathway from that hypothesized, hinting at quite complex mechanistic manifolds for organometallic catalysts that are affected by the combination of input variables. Mechanistic deviation is most likely when challenging multifunctional substrates are being used, and the quest for so-called privileged catalysts is quickly replaced by a need to screen catalyst libraries until a new "best" match between the catalyst and substrate can be identified and the reaction conditions can be optimized. As a community we remain confined to broad interpretations of the substrate scope of new catalysts and focus on small changes based on idealized catalytic cycles rather than working toward a "big data" view of organometallic homogeneous catalysis with routine use of predictive models and transparent data sharing.Databases of DFT-calculated steric and electronic descriptors can be built for such catalysts, and we summarize here how these can be used in the mapping, interpretation, and prediction of catalyst properties and reactivities. Our motivation is to make these databases useful as tools for synthetic chemists so that they challenge and validate quantitative computational approaches. In this Account, we demonstrate their application to different aspects of catalyst design and discovery and their integration with computational mechanistic studies and thus describe the progress of our journey toward truly predictive models in homogeneous organometallic catalysis.

Mapping the properties of bidentate ligands with calculated descriptors (LKB-bid)

Ligand space for bidentates has been mapped, computationally, varying donors, substituents and backbones, to give a new database, LKB-bid.

Computational Ligand Descriptors for Catalyst Design

Ligands, especially phosphines and carbenes, can play a key role in modifying and controlling homogeneous organometallic catalysts, and they often provide a convenient approach to fine-tuning the performance of known catalysts. The measurable outcomes of such catalyst modifications (yields, rates, selectivity) can be set into context by establishing their relationship to steric and electronic descriptors of ligand properties, and such models can guide the discovery, optimization, and design of catalysts. In this review we present a survey of calculated ligand descriptors, with a particular focus on homogeneous organometallic catalysis. A range of different approaches to calculating steric and electronic parameters are set out and compared, and we have collected descriptors for a range of representative ligand sets, including 30 monodentate phosphorus(III) donor ligands, 23 bidentate P,P-donor ligands, and 30 carbenes, with a view to providing a useful resource for analysis to practitioners. In addition, several case studies of applications of such descriptors, covering both maps and models, have been reviewed, illustrating how descriptor-led studies of catalysis can inform experiments and highlighting good practice for model comparison and evaluation.

Computational assessment on the Tolman cone angles for P-ligands

P-ligand cone angles have been properly recomputed in different transition-metal coordination linear, tetrahedral and octahedral environments by combining molecular mechanics and DFT methodologies.

Towards operando computational modeling in heterogeneous catalysis

An increased synergy between experimental and theoretical investigations in heterogeneous catalysis has become apparent during the last decade. Experimental work has extended from ultra-high vacuum and low temperature towards operando conditions. These developments have motivated the computational community to move from standard descriptive computational models, based on inspection of the potential energy surface at 0 K and low reactant concentrations (0 K/UHV model), to more realistic conditions. The transition from 0 K/UHV to operando models has been backed by significant developments in computer hardware and software over the past few decades. New methodological developments, designed to overcome part of the gap between 0 K/UHV and operando conditions, include (i) global optimization techniques, (ii) ab initio constrained thermodynamics, (iii) biased molecular dynamics, (iv) microkinetic models of reaction networks and (v) machine learning approaches. The importance of the transition is highlighted by discussing how the molecular level picture of catalytic sites and the associated reaction mechanisms changes when the chemical environment, pressure and temperature effects are correctly accounted for in molecular simulations. It is the purpose of this review to discuss each method on an equal footing, and to draw connections between methods, particularly where they may be applied in combination.

Exploiting coordination geometry to selectively predict the σ-donor and π-acceptor abilities of ligands: a back-and-forth journey between electronic properties and spectroscopy

Coordination geometry switches the carbonyl stretching frequency into a selective probe of the σ-donor and π-acceptor abilities of ligands.

On the relation between carbonyl stretching frequencies and the donor power of chelating diphosphines in nickel dicarbonyl complexes

The relation between spectroscopic observables and the detailed metal-ligand bonding features in chelation complexes is addressed using both experimental and state-of-the-art theoretical and computational methods. We synthesized and characterized a set of six nickel dicarbonyl complexes of general formula [Ni(CO)2(PP)], where PP is an atropoisomeric chelating diphosphine ligand. The analysis of the obtained experimental data and the basicity and oxidative potentials of the free ligands suggests a close relation between the donor ability of the chelating ligand and the carbonyl stretching frequencies observed in the complexes. We then use theory to unravel the detailed mechanisms of chelation-bond formation in terms of partial charge flows between the molecular orbitals of the fragments. By extending the promising, recently published natural orbitals for chemical valence/charge displacement (NOCV/CD) analysis scheme we provide a thorough, quantitative description of the several charge fluxes following the metal-ligand bond formation and demonstrate that the carbonyl stretching frequencies in the considered complexes selectively respond to the σ-donation charge flow from the phosphorus lone pairs of the ligands, with the frequency shift being in quantitative correlation with the extent of the ligand-to-metal charge transfer.

Palladium-catalyzed carbonylative transformation of aryl chlorides and aryl tosylates

The developments in the carbonylative transformations of aryl chlorides and aryl tosylates have been collected and discussed.

Lost in chemical space? Maps to support organometallic catalysis

Descriptors calculated from molecular structures have been used to map different areas of chemical space. A number of applications for such maps can be identified, ranging from the fine-tuning and optimisation of catalytic activity and compound properties to virtual screening of novel compounds, as well as the exhaustive exploration of large areas of chemical space by automated combinatorial building and evaluation. This review focuses on organometallic catalysis, but also touches on other areas where similar approaches have been used, with a view to assessing the extent to which chemical space has been explored.

Graphical abstract

Impact of Coordination Geometry, Bite Angle, and Trans Influence on Metal-Ligand Covalency in Phenyl-Substituted Phosphine Complexes of Ni and Pd

Despite the long-standing use of phosphine and diphosphine ligands in coordination chemistry and catalysis, questions remain as to their effects on metal-ligand bonding in transition metal complexes. Here we report ligand K-edge XAS, DFT, and TDDFT studies aimed at quantifying the impact of coordination geometry, diphosphine bite angle, and phosphine trans influence on covalency in M-P and M-Cl bonds. A series of four-coordinate NiCl2 and PdCl2 complexes containing PPh3 or Ph2P(CH2)nPPh2, where n = 1 (dppm), 2 (dppe), 3 (dppp), and 4 (dppb), was analyzed. The XAS data revealed that changing the coordination geometry from tetrahedral in Ni(PPh3)2Cl2 (1) to square planar in Ni(dppe)Cl2 (2) more than doubles the intensity of pre-edge features assigned to Ni-P and Ni-Cl 1s → σ* transitions. By way of comparison, varying the diphosphine in Pd(dppm)Cl2 (4), Pd(dppp)Cl2 (6), and Pd(dppb)Cl2 (7) yielded Pd-P 1s → σ* transitions with identical intensities, but a 10% increase was observed in the P K-edge XAS spectrum of Pd(dppe)Cl2 (5). A similar observation was made when comparing Ni(dppe)Cl2 (2) to Ni(dppp)Cl2 (3), and DFT and TDDFT calculations corroborated XAS results obtained for both series. Comparison of the spectroscopic and theoretical results to the diphosphine structures revealed that changes in M-P covalency were not correlated to changes in bite angles or coordination geometry. As a final measure, P and Cl K-edge XAS data were collected on trans-Pd(PPh3)2Cl2 (8) for comparison to the cis diphosphine complex Pd(dppe)Cl2 (5). Consistent with phosphine's stronger trans influence compared to chloride, a 35% decrease in the intensity of the Pd-P 1s → σ* pre-edge feature and a complementary 34% increase in Pd-Cl 1s → σ* feature was observed for 8 (trans) compared to 5 (cis). Overall, the results reveal how coordination geometry, ligand arrangement, and diphosphine structure affect covalent metal-phosphorus and metal-chloride bonding in these late transition metal complexes.

Tuneable reactivity with PPh3 and SnX2 of four- and five-coordinate Pd(II) and Pt(II) complexes containing polyphosphines

The reactivity of the unusual d(8) trigonal-bipyramidal systems [MX(PP3)]X (X = Cl: M = Pd(1a), Pt(2a); X = Br: M = Pd(3a), Pt(4a); X = I: M = Pd(5a), Pt(6a); PP3 = tris[2-(diphenylphosphino)ethyl]phosphine) in CHCl3-CH3OH, the square-pyramidal compounds [MCl(NP3)]Cl (M = Pd(7a); Pt(8a); NP3 = tris[2-(diphenylphosphino)ethyl]amine) in CD3OD-DMF and the distorted square-planar mononuclear [MX(PNP)]X (M = Pd: X = Cl(10a); M = Pt: X = I(10b); PNP = bis[2-(diphenylphosphino)ethyl]amine) and the heteronuclear [PdAu2X4(PP3)] [X = I(9a), Cl(14a), Br(15a)] and [MAuX2(PP3)]X [M = Pd: X = Cl(16a); M = Pt: X = Cl(17a), Br(18a)] species in CDCl3 with PPh3 + SnX2 has been explored to establish the factors that influence the nature of the products. With the mononuclear precursors the course of the reaction is strongly dependent on the tripodal or linear arrangement of the polydentate ligand and in the former case on the halogen. Thus, while for chlorides (1a-2a, 7a-8a) and bromides (3a-4a) the reaction led to the trigonal-bipyramidal compounds [M(SnCl3)(AP3)][SnCl3] [A = P: M = Pd(1), Pt(2); A = N: M = Pd(7), Pt(8)], [MBr(PP3)][SnBr3] [M = Pd(4), Pt(6)] containing M-Sn and M-Br bonds, respectively, for iodides (5a-6a) resulted in the unknown neutral square-planar compounds [MI2(PP(PO)2)(SnI2)2] [M = Pd(9) and Pt(10)] bearing two dangling P=O-SnI2 units and P2MI2 environments. However, complexes of the type [PtCl(PP2PO)X]X' [X = SnCl2, X' = [SnCl3](-)(11)] and [M(PP(PO)2)2X4]X'2 [X = SnCl2, X' = [SnCl3](-): M = Pd(12), Pt(13)] showing P=O-SnCl2 arms were obtained by direct reaction of [PtCl(PP2PO)]Cl (11a) and [M(PP(PO)2)2]Cl2 [M = Pd(12a), Pt(13a)] with SnCl2 in CH3OH. Although complex 9 was also prepared by interaction of the heteronuclear iodide 9a with PPh3 + SnI2 in CDCl3, the use of the neutral and ionic heteronuclear chlorides and bromides (14a-18a) as starting materials afforded the distorted square-planar ionic systems [MAuX'(PP3)(PPh3)][SnX3]2 [M = Pd: X = Cl, X' = SnCl3(-)(14); X = Br, X' = SnBr3(-)(15); M = Pt: X = Cl, X' = SnCl3(-)(17); X = Br, X' = SnBr3(-)(18)] containing M-SnX3 and P-Au-PPh3 functionalities. It was found that these reactions where the heteronuclear species are the precursors proceed via the trigonal-bipyramidal halides not only with X = Cl and Br(1a-4a) but also I(5a). When the precursors were 10a and 10b the reaction occurred with formation of [Pd(PNP)(PPh3)][SnCl3]2 (23) and [Pt(PNP)(PPh3)][SnCl2I]2 (24) showing M-PPh3 units and trihalostannato counter anions.