Theoretical studies in palladium and platinum molecular chemistry

Palladium-catalyzed activation of HnA–AHn bonds (AHn = CH3, NH2, OH, F)

We have quantum chemically studied activation of HnA–AHn bonds (AHn = CH3, NH2, OH, F) by PdLn catalysts with Ln = no ligand, PH3, (PH3)2, using relativistic density functional theory at ZORA-BLYP/TZ2P. The activation energy associated with the oxidative addition step decreases from H3C–CH3 to H2N–NH2 to HO–OH to F–F, where the activation of the F–F bond is barrierless. Activation strain and Kohn–Sham molecular orbital analyses reveal that the enhanced reactivity along this series of substrates originates from a combination of (i) reduced activation strain due to a weaker HnA–AHn bond; (ii) decreased Pauli repulsion as a result of a difference in steric shielding of the HnA–AHn bond; and (iii) enhanced backbonding interaction between the occupied 4d atomic orbitals of the palladium catalyst and σ* acceptor orbital of the substrate.

C−X Bond Activation by Palladium: Steric Shielding versus Steric Attraction

The C−X bond activation (X = H, C) of a series of substituted C(n°)−H and C(n°)−C(m°) bonds with C(n°) and C(m°) = H₃C− (methyl, 0°), CH₃H₂C− (primary, 1°), (CH₃)₂HC− (secondary, 2°), (CH₃)₃C− (tertiary, 3°) by palladium were investigated using relativistic dispersion‐corrected density functional theory at ZORA‐BLYP‐D3(BJ)/TZ2P. The effect of the stepwise introduction of substituents was pinpointed at the C−X bond on the bond activation process. The C(n°)−X bonds become substantially weaker going from C(0°)−X, to C(1°)−X, to C(2°)−X, to C(3°)−X because of the increasing steric repulsion between the C(n°)‐ and X‐group. Interestingly, this often does not lead to a lower barrier for the C(n°)−X bond activation. The C−H activation barrier, for example, decreases from C(0°)−X, to C(1°)−X, to C(2°)−X and then increases again for the very crowded C(3°)−X bond. For the more congested C−C bond, in contrast, the activation barrier always increases as the degree of substitution is increased. Our activation strain and matching energy decomposition analyses reveal that these differences in C−H and C−C bond activation can be traced back to the opposing interplay between steric repulsion across the C−X bond versus that between the catalyst and substrate.

Predicting Pt-195 NMR Chemical Shift and 1J(195Pt-31P) Coupling Constant for Pt(0) Complexes Using the NMR-DKH Basis Sets

Pt(0) complexes have been widely used as catalysts for important reactions, such as the hydrosilylation of olefins. In this context, nuclear magnetic resonance (NMR) spectroscopy plays an important role in characterising of new structures and elucidating reaction mechanisms. In particular, the Pt-195 NMR is fundamental, as it is very sensitive to the ligand type and the oxidation state of the metal. In the present study, quantum mechanics computational schemes are proposed for the theoretical prediction of the Pt-195 NMR chemical shift and 1J(195Pt–31P) in Pt(0) complexes. The protocols were constructed using the B3LYP/LANL2DZ/def2-SVP/IEF-PCM(UFF) level for geometry optimization and the GIAO-PBE/NMR-DKH/IEF-PCM(UFF) level for NMR calculation. The NMR fundamental quantities were then scaled by empirical procedures using linear correlations. For a set of 30 Pt(0) complexes, the results showed a mean absolute deviation (MAD) and mean relative deviation (MRD) of only 107 ppm and 2.3%, respectively, for the Pt-195 NMR chemical shift. When the coupling constant is taken into account, the MAD and MRD for a set of 33 coupling constants in 26 Pt(0) complexes were of 127 Hz and 3.3%, respectively. In addition, the models were validated for a group of 17 Pt(0) complexes not included in the original group that had MAD/MRD of 92 ppm/1.7% for the Pt-195 NMR chemical shift and 146 Hz/3.6% for the 1J(195Pt–31P).

In silico rationalisation of selectivity and reactivity in Pd-catalysed C-H activation reactions

A computational approach has been developed to automatically generate and analyse the structures of the intermediates of palladium-catalysed carbon-hydrogen (C-H) activation reactions as well as to predict the final products. Implemented as a high-performance computing cluster tool, it has been shown to correctly choose the mechanism and rationalise regioselectivity of chosen examples from open literature reports. The developed methodology is capable of predicting reactivity of various substrates by differentiation between two major mechanisms - proton abstraction and electrophilic aromatic substitution. An attempt has been made to predict new C-H activation reactions. This methodology can also be used for the automated reaction planning, as well as a starting point for microkinetic modelling.

Multi gas sensors using one nanomaterial, temperature gradient, and machine learning algorithms for discrimination of gases and their concentration

In this work, four identical micro sensors on the same chip with noble metal decorated tin oxide nanowires as gas sensing material were located at different distances from an integrated heater to work at different temperatures. Their responses are combined in highly informative 4D points that can qualitatively (gas recognition) and quantitatively (concentration estimate) discriminate all the tested gases. Two identical chips were fabricated with tin oxide (SnO₂) nanowires decorated with different metal nanoparticles: one decorated with Ag nanoparticles and one with Pt nanoparticles. Support Vector Machine was used as the "brain" of the sensing system. The results show that the systems using these multisensor chips were capable of achieving perfect classification (100%) and good estimation of the concentration of tested gases (errors in the range 8-28%). The Ag decorated sensors did not have a preferential gas, while Pt decorated sensors showed a lower error towards acetone, hydrogen and ammonia. Combination of the two sensor chips improved the overall estimation of gas concentrations, but the individual sensor chips were better for some specific target gases.

Hierarchy of Commonly Used DFT Methods for Predicting the Thermochemistry of Rh-Mediated Chemical Transformations

The accuracy and reliability of 17 commonly used density functionals in conjunction with Poisson-Boltzmann finite solvation model were gauged for predicting the free energy of Rh(I)- and Rh(III)-mediated chemical transformations such as ligand exchange, hydride elimination, dihydrogen elimination, chloride affinity, and silyl hydride bond activation reactions. In total, six Rh-mediated reactions were examined, and the computed density functional theory results were then subjected to comparison with the experimentally reported values. For reaction A, involving replacement of N₂ with η²-H₂ over Rh(I), MPWB1K-D3, B3PW91, B3LYP, and BHandHYLP emerged to be the best functionals of all the tested methods in terms of their deviations ≤2 kcal mol^-1 from experimental data. For reaction B, in which exchange of η²-C₂H₄ with N₂ over Rh(I) takes place, MPWB1K-D3 and M06-2X-D3 functionals performed the best, while as for reaction C (hydride elimination reaction in Rh(III) complex), it is PBE functional that showed impressive performance. Similarly, for reaction D (H₂ elimination reaction in Rh(III) complex), PBE0-D3 and PBE-D3 showed exceptional results compared to other functionals. For reaction E (H₂O/Cl^- exchange), the PBE0 again shows impressive performance as compared to other functionals. For reaction F (Si-H activation), M06-2X-D3, PBE0-D3, and MPWB1K-D3 functionals are undoubtedly the best functionals. Overall, PBE0-D3 and MPWB1K-D3 functionals were impressive in all cases with lowest mean unsigned errors (3.2 and 3.4 kcal mol^-1, respectively) with respect to experimental Gibbs free energies. Thus, these two functionals are recommended for studying Rh-mediated chemical transformations.

Ab initio calculations on the ground and excited electronic states of neutral and charged palladium monoxide, PdO0,+,−

Multi-reference configuration interaction and coupled cluster calculations were carried out for the ground and several low-lying excited electronic states for PdO, PdO⁺, and PdO⁻. The photoelectron spectrum peaks of PdO were assigned.

Platinum for hydrogen sensing: surface and grain boundary scattering antagonistic effects in Pt@Au core-shell nanoparticle assemblies prepared using a Langmuir-Blodgett method

Hydrogen resistive sensors are fabricated through the synthesis of a series of Pt@Au core-shell nanoparticles showing various Pt shell thicknesses. Resulting colloids are assembled as hexagonal close-packed 2D monolayers of various dimension characteristics using a simple Langmuir-Blodgett method. The fabricated sensors show attractive hydrogen sensing performances with reversible responses in extended sensing ranges, a good specificity towards H₂, short response and recovery times… Sensing measurements and data analyses allow the demonstration of the associated sensing mechanisms. The dissociative chemisorption of H₂ and O₂ on the Pt surface through a Langmuir-Hinshelwood mechanism leads to the formation of chemisorbed hydrogen and hydroxyl groups. This surface nature change induces the modification of the scattering of the conduction electrons at both the grain surface and intercontacts, tuned by the extent of hydrogen and hydroxyl group coverages. In assemblies made of particles showing thin Pt shells, the predominance of the surface scattering described by the Fuchs-Sondheimer model accounts for the observed conductive responses as the number of involved grain boundaries is limited. In contrast, in assemblies made of particles with thick Pt shells, the scattering at the grain boundaries described by the Mayadas-Shatzkes model mostly contributes to the observed resistive responses. The sensor behavior is balanced by these two antagonistic effects.

Modeling σ-Bond Activations by Nickel(0) Beyond Common Approximations: How Accurately Can We Describe Closed-Shell Oxidative Addition Reactions Mediated by Low-Valent Late 3d Transition Metal?

Accurate modelings of reactions involving 3d transition metals (TMs) are very challenging to both ab initio and DFT approaches. To gain more knowledge in this field, we herein explored typical σ-bond activations of H-H, C-H, C-Cl, and C-C bonds promoted by nickel(0), a low-valent late 3d TM. For the key parameters of activation energy (ΔE^‡) and reaction energy (ΔE_R) for these reactions, various issues related to the computational accuracy were systematically investigated. From the scrutiny of convergence issue with one-electron basis set, augmented (A) basis functions are found to be important, and the CCSD(T)/CBS level with complete basis set (CBS) limit extrapolation based on augmented double-ζ and triple-ζ basis pair (ADZ and ATZ), which produces deviations below 1 kcal/mol from the reference, is recommended for larger systems. As an alternative, the explicitly correlated F12 method can accelerate the basis set convergence further, especially after its CBS extrapolations. Thus, the CCSD(T)-F12/CBS(ADZ-ATZ) level with computational cost comparable to the conventional CCSD(T)/CBS(ADZ-ATZ) level, is found to reach the accuracy of the conventional CCSD(T)/A5Z level, which produces deviations below 0.5 kcal/mol from the reference, and is also highly recommendable. Scalar relativistic effects and 3s3p core-valence correlation are non-negligible for achieving chemical accuracy of around 1 kcal/mol. From the scrutiny of convergence issue with the N-electron basis set, in comparison with the reference CCSDTQ result, CCSD(T) is found to be able to calculate ΔE^‡ quite accurately, which is not true for the ΔE_R calculations. Using highest-level CCSD(T) results of ΔE^‡ in this work as references, we tested 18 DFT methods and found that PBE0 and CAM-B3LYP are among the three best performing functionals, irrespective of DFT empirical dispersion correction. With empirical dispersion correction included, ωB97XD is also recommendable due to its improved performance.

Metal Ion Modeling Using Classical Mechanics

Metal ions play significant roles in numerous fields including chemistry, geochemistry, biochemistry, and materials science. With computational tools increasingly becoming important in chemical research, methods have emerged to effectively face the challenge of modeling metal ions in the gas, aqueous, and solid phases. Herein, we review both quantum and classical modeling strategies for metal ion-containing systems that have been developed over the past few decades. This Review focuses on classical metal ion modeling based on unpolarized models (including the nonbonded, bonded, cationic dummy atom, and combined models), polarizable models (e.g., the fluctuating charge, Drude oscillator, and the induced dipole models), the angular overlap model, and valence bond-based models. Quantum mechanical studies of metal ion-containing systems at the semiempirical, ab initio, and density functional levels of theory are reviewed as well with a particular focus on how these methods inform classical modeling efforts. Finally, conclusions and future prospects and directions are offered that will further enhance the classical modeling of metal ion-containing systems.

A Square-Planar Complex of Platinum(0)

The Pt⁰ complex [Pt(PPh₃ )(Eind₂ -BPEP)] with a pyridine-based PNP-pincer-type phosphaalkene ligand (Eind₂ -BPEP) has a highly planar geometry around Pt with ∑(Pt)=358.6°. This coordination geometry is very uncommon for formal d¹⁰ complexes, and the Pd and Ni homologues with the same ligands adopt distorted tetrahedral geometries. DFT calculations reveal that both the Pt and Pd complexes are M⁰ species with nearly ten valence electrons on the metals whereas their atomic orbital occupancies are evidently different from one another. The Pt complex has a higher occupancy of the atomic 6s orbital because of strong s-d hybridization due to relativistic effects, thereby adopting a highly planar geometry reflecting the shape and orientation of the partially unoccupied dx2-y2 orbital.

Shielding the chemical reactivity using graphene layers for controlling the surface properties of carbon materials

Graphene can efficiently shield chemical interactions and gradually decrease the binding to reactive defect areas. In the present study, we have used the observed graphene shielding effect to control the reactivity patterns on the carbon surface. The experimental findings show that a surface coating with a tiny carbon layer of 1-2 nm thickness is sufficient to shield the defect-mediated reactivity and create a surface with uniform binding ability. The shielding effect was directly observed using a combination of microscopy techniques and evaluated with computational modeling. The theoretical calculations indicate that a few graphene layers can drastically reduce the binding energy of the metal centers to the surface defects by 40-50 kcal mol(-1). The construction of large carbon areas with controlled surface reactivity is extremely difficult, which is a key limitation in many practical applications. Indeed, the developed approach provides a flexible and simple tool to change the reactivity patterns on large surface areas within a few minutes.

Allylic amination reactivity of Ni, Pd, and Pt heterobimetallic and monometallic complexes

Transition metal heterobimetallic complexes with dative metal–metal interactions have the potential for novel fast reactivity.

Reversible [4 + 2] cycloaddition reaction of 1,3,2,5-diazadiborinine with ethylene

Under ambient conditions, a [4 + 2] cycloaddition reaction of 1,3,2,5-diazadiborinine 1 with ethylene afforded a bicyclo[2.2.2] derivative 2, which was structurally characterized.

Under ambient conditions, a [4 + 2] cycloaddition reaction of 1,3,2,5-diazadiborinine 1 with ethylene afforded a bicyclo[2.2.2] derivative 2, which was structurally characterized. The cyclization process was found to be reversible, and thus retro-[4 + 2] cycloaddition reproduced 1 quantitatively, concomitant with the release of ethylene. Compound 1 reacted regio-selectively and stereo-selectively with styrene derivatives and norbornene, respectively, and these processes were found to be reversible too. Computational studies determined the reaction pathways which were consistent with the regio-selectivity observed in the reaction of styrene, and the reaction was suggested to be essentially concerted but highly asynchronous.

Highly selective and sensitive fluorescence chemosensor for the detection of palladium species based on Tsuji-Trost reaction

A new chemosensor 7-nitro-2,1,3-benzoxadiazole-4-allyl-N-(thiophen-2-ylmethyl)carbamate (NBDTC) was synthesized and utilized for palladium detection based on the Tsuji-Trost reaction. NBDTC displayed specific and ratiometric fluorescent responses toward palladium species. The chemosensor showed more than 50-fold enhancement in fluorescence intensity with the presence of PEG400 and palladium because NBDTC can be transformed to NBDT under palladium-catalyzing Tsuji-Trost reaction. NBDTC displayed high selectivity and sensitivity for palladium species with the detection limit of 1.13×10(-9) M.

Activation of Strong Boron-Fluorine and Silicon-Fluorine σ-Bonds: Theoretical Understanding and Prediction

The oxidative addition of BF3 to a platinum(0) bis(phosphine) complex [Pt(PMe3)2] (1) was investigated by density functional calculations. Both the cis and trans pathways for the oxidative addition of BF3 to 1 are endergonic (ΔG°=26.8 and 35.7 kcal mol(-1), respectively) and require large Gibbs activation energies (ΔG°(≠)=56.3 and 38.9 kcal mol(-1), respectively). A second borane plays crucial roles in accelerating the activation; the trans oxidative addition of BF3 to 1 in the presence of a second BF3 molecule occurs with ΔG°(≠) and ΔG° values of 10.1 and -4.7 kcal mol(-1), respectively. ΔG°(≠) becomes very small and ΔG° becomes negative. A charge transfer (CT), F→BF3, occurs from the dissociating fluoride to the second non-coordinated BF3. This CT interaction stabilizes both the transition state and the product. The B-F σ-bond cleavage of BF2Ar(F) (Ar(F)=3,5-bis(trifluoromethyl)phenyl) and the B-Cl σ-bond cleavage of BCl3 by 1 are accelerated by the participation of the second borane. The calculations predict that trans oxidative addition of SiF4 to 1 easily occurs in the presence of a second SiF4 molecule via the formation of a hypervalent Si species.

Selective C-H and C-C Bond Activation: Electronic Regimes as a Tool for Designing d(10) MLn Catalysts

We wish to understand how a transition-metal catalyst can be rationally designed so as to selectively activate one particular bond in a substrate, herein, C-H and C-C bonds in ethane. To this end, we quantum chemically analyzed the activity and selectivity of a large series of model catalysts towards ethane and, for comparison, methane, by using the activation strain model and quantitative molecular orbital theory. The model catalysts comprise d(10) MLn complexes with coordination numbers n=0, 1, and 2; metal centers M=Co(-), Rh(-), Ir(-), Ni, Pd, Pt, Cu(+), Ag(+), and Au(+); and ligands L=NH3, PH3, and CO. Our analyses reveal that rather subtle electronic differences between bonds can be exploited to induce a lower barrier for activating one or the other, depending, among other factors, on the catalysts electronic regime (i.e., s-regime versus d-regime catalysts). Interestingly, the concepts and design principles emerging from this work can also be applied to the more challenging problem of differentiating between activation of the C-H bonds in ethane versus those in methane.

Computational aspects of hydroformylation

This review is to focus on computational studies on hydroformylation and theoretical coordination chemistry results related to hydroformylation.

Pd-catalyzed aerobic oxidative coupling of arenes: evidence for transmetalation between two Pd(II)-aryl intermediates

Pd-catalyzed aerobic oxidative coupling of arenes provides efficient access to biaryl compounds. The biaryl product forms via C-H activation of two arenes to afford a Pd(II)ArAr' intermediate, which then undergoes C-C reductive elimination. The key Pd(II)ArAr' intermediate could form via a "monometallic" pathway involving sequential C-H activation at a single Pd(II) center, or via a "bimetallic" pathway involving parallel C-H activation at separate Pd(II) centers, followed by a transmetalation step between two Pd(II)-aryl intermediates. Here, we investigate the oxidative coupling of o-xylene catalyzed by a PdX2/2-fluoropyridine catalyst (X = trifluoroacetate, acetate). Kinetic studies, H/D exchange experiments, and kinetic isotope effects provide clear support for a bimetallic/transmetalation mechanism.

Mechanism, reactivity, and selectivity in palladium-catalyzed redox-relay Heck arylations of alkenyl alcohols

The enantioselective Pd-catalyzed redox-relay Heck arylation of acyclic alkenyl alcohols allows access to various useful chiral building blocks from simple olefinic substrates. Mechanistically, after the initial migratory insertion, a succession of β-hydride elimination and migratory insertion steps yields a saturated carbonyl product instead of the more general Heck product, an unsaturated alcohol. Here, we investigate the reaction mechanism, including the relay function, yielding the final carbonyl group transformation. M06 calculations predict a ΔΔG^⧧ of 1 kcal/mol for the site selectivity and 2.5 kcal/mol for the enantioselectivity, in quantitative agreement with experimental results. The site selectivity is controlled by a remote electronic effect, where the developing polarization of the alkene in the migratory insertion transition state is stabilized by the C–O dipole of the alcohol moiety. The enantioselectivity is controlled by steric repulsion between the oxazoline substituent and the alcohol-bearing alkene substituent. The relay efficiency is due to an unusually smooth potential energy surface without high barriers, where the hydroxyalkyl-palladium species acts as a thermodynamic sink, driving the reaction toward the carbonyl product. Computational predictions of the relative reactivity and selectivity of the double bond isomers are validated experimentally.

Structural, spectroscopic and theoretical studies of a vapochromic platinum(ii) terpyridyl complex

A combination of spectroscopic methods, crystallographic analyses, and theoretical studies provides a rationale to understand the nature of the structural and electronic response of the simplest platinum(ii) terpyridyl complex to acetonitrile.