Tyler M.
Porter
,
Andrew L.
Ostericher
and
Clifford P.
Kubiak
*
Department of Chemistry and Biochemistry, University of California San Diego, 9500 Gilman Drive, La Jolla, California, USA. E-mail: ckubiak@ucsd.edu
First published on 11th July 2019
Synthetic control of the influence of steric and electronic factors on the ultrafast (picosecond) isomerization of penta-coordinate ruthenium dithietene complexes (Ru((CF3)2C2S2)(CO)(L)2, where L = a monodentate phosphine ligand) is reported. Seven new ruthenium dithietene complexes were prepared and characterized by single crystal X-ray diffraction. The complexes are all square pyramidal and differ only in the axial vs. equatorial coordination of the carbonyl ligand. Fourier Transform Infrared (FTIR) spectroscopy was used to study the ν(CO) bandshapes of the complexes in solution, and these reveal rapid exchange between two or three isomers of each complex. Isomerization is proposed to follow a Berry psuedorotation-like mechanism where a metastable, trigonal bipyramidal (TBP) intermediate is observed spectroscopically. Electronic tuning of the phosphine ligands L = PPh3, P((p-Me)Ph)3, ((p-Cl)Ph)3, at constant cone angle is found to have little effect on the kinetics or thermodynamic stabilities of the axial, equatorial and TBP isomers of the differently substituted complexes. Steric tuning of the phosphine ligands over a range of phosphine cone angles (135 < θ < 165°) has a profound impact on the isomerization process, and in the limit of greatest steric bulk, the axial isomer is not observable. Temperature dependence of the FTIR spectra was used to obtain the relative thermodynamic stabilities of the different isomers of each of the seven ruthenium dithietene complexes. This study details how ligand steric effects can be used to direct the solution state dynamics on the picosecond time scale of discrete isomers energetically separated by <2.2 kcal mol−1. This work provides the most detailed description to date of ultrafast isomerization in the ground states of transition metal complexes.
Fig. 1 Isomerization process for 1 (top).1 The eight structurally characterized ruthenium complexes and their respective phosphine cone angles (bottom). |
Through the use of 2D IR spectroscopy, the equatorial (1a) and axial (1c) isomers of complex 1 were observed to exchange in picoseconds (kex ≈ 10−12) through an observable, trigonal bipyramidal (TBP) intermediate (1b).1 Additional spectroscopic and DFT analyses suggested that the mechanism of isomerization followed a Berry pseudorotation-like (BPR) pathway (1a ↔ 1b ↔ 1c). To further probe the steric and electronic effects, that govern this low barrier isomerization, we investigated seven new ruthenium dithietene complexes and report herein on their synthesis, characterization, and exchange dynamics as a function of the steric and electronic character of the phosphine ligands.
Here, the steric and electronic effects of the monodentate phosphine ligands are readily accounted for by consideration of the phosphine cone angles and their relative σ-donor strengths.6 Complexes 1–3 address the influence of phosphine electronic effects, as para-substitution on the phenyl rings perturbs the σ-donor strength while the cone angle remains fixed. In complexes 4–8 however, the σ-donor strength of the phosphine ligand remains relatively constant, while the cone angle can be varied across a 30° span (135° < θ < 165°).
Single crystal X-ray diffraction studies revealed that all complexes were square-pyramidal and differed only in the orientation of the carbonyl ligand about the metal center. The orange crystals (2–6) were found to have the CO ligand in the apical site while the violet crystals (7 and 8) were found to have the CO ligand in the equatorial plane (Fig. 2). While each complex is predicted to have additional stable isomers,1–5 their crystallization is remarkably condition dependent and any additional isomers have yet to be crystallized. It is rather surprising however, that 7 and 8 crystallize as the COequatorial isomer as previous studies indicated it is the thermodynamically less favorable isomer.1,2 Nevertheless, this does reveal initial effects of phosphine substitution, and suggests that increased ligand sterics favor the sterically less encumbered COequatorial isomer.3 Minor variations in the bonding parameters are observed (Tables 1 and S9†), suggesting a similar electronic structure for all complexes. Comparison of the C–C and C–S bond distances of the 1,2-dithioketone ligand allow determination of the dithioketone (and metal) oxidation state.7–9 For complexes 1–8, average C–C and C–S bond distances of 1.354 Å and 1.733 Å, respectively, are in excellent agreement with a doubly reduced dithietene ligand (C–C: 1.337 Å and C–S: 1.761 Å),7,9 suggesting a Ru(II), 16 e− complex. This assignment is additionally supported through preparation of the hexa-coordinate, trisphosphine complex, Ru((CF3)2C2S2)(CO)(P(Ph)2(Me))3 (9, Fig. 3).
Fig. 2 ORTEP structures of complexes 2–8. Thermal ellipsoids are shown at 50% probability. Hydrogen atoms and co-crystallized solvent molecules have been omitted for clarity and a single position is shown where rotational disorder of the CF3 groups exists. Notes on refinement and crystallographic data is included in the ESI.† |
Bond (Å) | 1 | 2 | 3 | 4 | 5 | 6 | 7 | 8 | 9 | Average | |
---|---|---|---|---|---|---|---|---|---|---|---|
COap | COeq | COap | COap | COap | COap | COap | COeq | COeq | (P(Ph)2Me)3 | ||
a Values were adapted from previously published structures.3–5 b Crystals were ground from acetonitrile solutions and the structure was found to contain a coordinating acetonitrile molecule. | |||||||||||
C(2)–C(3) | 1.350 | 1.359 | 1.357 | 1.359 | 1.354 | 1.356 | 1.343 | 1.347 | 1.350 | 1.350 | 1.353 |
C(2)–S(1) | 1.735 | 1.718 | 1.728 | 1.727 | 1.735 | 1.733 | 1.762 | 1.732 | 1.741 | 1.765 | 1.735 |
C(3)–S(2) | 1.726 | 1.717 | 1.741 | 1.737 | 1.740 | 1.737 | 1.756 | 1.737 | 1.727 | 1.752 | 1.735 |
Based upon simple principles of electron counting, complex 9 is expected to have a doubly reduced dithiolate ligand, giving the 18 e−, Ru(II) complex. This is in excellent agreement with observed C–C and C–S bond distances of 1.350 Å and 1.765/1.752 Å respectively, supporting a Ru(II) assignment for 9, and 1–8 as well. It is important to note that the tris-phosphine adduct has only been isolated for complex 4, where the relatively small cone angle (θ = 135°) is believed to allow the coordination of an additional phosphine ligand. While complex 9 does not appear to undergo dynamic exchange on the picosecond timescale (as judged by the narrow ν(CO) band shape, Fig. S1†), significantly broadened, partially coalesced NMR signals do suggest dynamics on the NMR timescale (kex < 104 s−1, Fig. S29–S30†). Clearly, additional studies are required, however this remains outside of the scope of the current investigation and here 9 aids only in providing a structural benchmark for Ru(II)/dithiolate(2-) oxidation states.
Upon dissolution of 2–8 into dichloromethane (DCM) a unique, significantly broadened ν(CO) band is observed for each complex, behavior indicative of exchange dynamics occurring on the vibrational timescale (Fig. 4, S2 and S3).1,10–19 Owing to the similar line shape observed for complexes 1–3 (Fig. S2†) phosphine electronic effects appear to have little impact on the isomerization reaction. However, variation of the phosphine cone angles is observed to have substantial impact (Fig. 4). At small ligand cone angles (θ < 145°, 4 and 5) two significantly broadened, partially coalesced ν(CO) bands are observed (Fig. 4) in the vicinity of 1936 and 1970 cm−1. In comparison to complex 1, the observed ν(CO) bands are in excellent agreement with the COequatorial (1a: ν(CO) ≈ 1960 cm−1) and COaxial (1c: ν(CO) ≈ 1940 cm−1) isomers, while the partially coalesced spectra suggests relatively slow exchange dynamics (in comparison to the highly coalesced ν(CO) bands for 1–3 and 6–8). Interestingly, two-component spectral deconvolutions are less than sufficient to reproduce the observed FTIR line shape (Fig. S4†) and instead, at-least three Voigt functions are required to fit the ν(CO) bands in the vicinity of 1975 cm−1, 1965 cm−1, and 1932 cm−1. These findings are again in excellent agreement with the reported frequencies for complex 1 (ν(CO) ≈ 1a: 1950 cm−1, 1b: 1980 cm−1, 1c: 1960 cm−1),1 and additionally suggests the presence of the COTBP intermediate.
Fig. 4 Variable temperature Fourier transform infrared (VT-FTIR) spectra of complexes 2 and 4–8 in DCM across a temperature range of 20 °C to −70 °C. |
As the ligand cone angle is increased (θ ≈ 145°, 2, 3, and 6), the ν(CO) band in solution is observed to coalesce into a broad absorbance (FWHM ≈ 50 cm−1) centered near 1958 cm−1, where the nearly identical ν(CO) line shapes (in comparison to 1) suggest very similar solution state dynamics and energetics. Thus, similarly sized phosphine ligands display similar isomerization dynamics, even when their electronic differences are significant. These observations are confirmed through spectral deconvolution, where again, three Voigt functions in the vicinity of 1978 cm−1, 1960 cm−1, and 1940 cm−1 are required to fit the experimental line shape (Fig. S5†). These findings in addition to the nearly identical variable temperature response (Fig. S2†) support very similar solution state dynamics for complexes 1–3 and 6.
At large cone angles (θ > 145°, 7 and 8), a near complete coalescence of the ν(CO) band is observed, and unlike complexes 1–6, only two-component deconvolutions were required (Fig. S6†). While this behavior would initially indicate faster exchange rates in comparison to complexes 1–6, a more likely explanation is inhibition of one of the isomerization pathways. Here, the deconvoluted center line frequencies (ν(CO) ≈ 1960 cm−1 and 1975 cm−1) are in excellent agreement with observation of the COequatorial and COTBP isomers, with a notable lack of intensity in the vicinity of 1940 cm−1. This suggests an absence of the COaxial isomer, which is well supported by crystallographic and density functional theory (DFT) studies (vide infra), where the COaxial isomer is believed to be too sterically encumbered with high cone angle phosphine ligands.
The thermodynamics of the isomerization reaction were investigated through the use of variable temperature Fourier transform infrared (VT-FTIR) spectroscopy. The VT-FTIR spectra of complexes 2–8 (Fig. 4, S2 and S3†) were recorded using a SPECAC flow-through optical cryostat (Model: GS21525-C) and a Bruker Equinox 55 FTIR spectrometer. The samples were enclosed in a CaF2 crystal windowed FTIR flow cell, contained in a vacuum jacketed housing, that was cooled using a liquid nitrogen/ethanol slurry to within 1 °C of the desired temperature using a built-in temperature controller. FTIR spectra for each complex were collected across a 90 °C span from 20 °C to −70 °C. After solvent subtraction, the VT-FTIR line shapes for each complex were fit using constrained Voigt functionals in accord with the aforementioned spectral deconvolutions (Fig. S10–S13†). The Voigt functionals were constrained by ±5 cm−1 about the isomers center line frequencies, with FWHM constraints of 1–20 cm−1, and peak heights of 0.001–0.05 a.u. Using the determined spectral areas the population ratios of the isomers for complexes 4–8 were then determined (Tables S11–S17†), allowing application of a Van 't Hoff analysis (Fig. S14† and Table 2).
Complex | Pathway | (kcal mol−1) | ΔHi (kcal mol−1) | ΔSi (e.u.) |
---|---|---|---|---|
4 | COeq to COTBP | 0.9 (0.2) | 1.4 (0.3) | 1.9 (0.6) |
COTBP to COax | −0.8 (0.2) | −2.2 (0.2) | 4.7 (0.9) | |
COeq to COax | 0.094 (0.007) | −0.76 (0.08) | −2.9 (0.4) | |
5 | COeq to COTBP | 1.6 (0.2) | 1.96 (0.09) | 1.1 (0.4) |
COTBP to COax | −1.4 (0.2) | −2.1 (0.1) | −2.2 (0.5) | |
COeq to COax | 0.26 (0.01) | −0.08 (0.04) | −1.2 (0.2) | |
7 | COeq to COTBP | 1.29 (0.04) | 0.95 (0.03) | −1.2 (0.1) |
8 | COeq to COTBP | −0.40 (0.01) | 0.62 (0.03) | 3.4 (0.1) |
In accord with previous studies,1 the trigonal bipyramidal isomer was found to be the highest in energy (ΔHi), lying on average 1.2 kcal mol−1 and 2.1 kcal mol−1 higher in energy than the COequatorial and COaxial isomers, respectively. Overall, isomerization from the COequatorial to COTBP was found to be an endergonic process while isomerization from COTBP to COaxial was exergonic (Table 2). While isomerization between the COequatorial and COTBP isomers was observed for complexes 7 and 8; interestingly, the isomerization from COequatorial to COTBP in 8 was found to be exergonic, suggesting an entropically driven isomerization.
While ligand sterics likely play a major role in the thermal isomerization process, the changes in entropy across the entire series would suggest that solvation dynamics may also contribute to the overall reaction. While solvent studies are currently underway, this concept is not surprising as the reordering of solvent about each isomer would be expected to play a key role. This is especially true for an isomerization reaction occurring on the same dynamic timescale as solvent reorientation.20–24 It is important to note however, that the thermochemical values are highly dependent on the peak locations of the spectral deconvolutions, and therefore should be only taken as estimates. Higher ordered 2D IR spectroscopy will be required to better resolve the role of steric effects in setting the relative stabilities of the different isomers, and the dynamics of exchange between them.
Despite the limitations of linear spectroscopy, additional support for these observations is obtained through the use of density functional theory (DFT). Due to computational time constraints however, DFT calculations were only performed on complexes 5 and 7 using the ORCA computational software package (version 3.0.3).25–29 At the BP86 level of theory, all hydrogen, carbon, and fluorine atoms were treated with Ahlrichs' def2-SVP basis set while all sulfur, phosphorus, oxygen, and ruthenium atoms were treated with the def2-TZVP basis set.30–38 Relativistic effects were applied through the zero-order regular approximation (ZORA), dispersion corrections were applied with a Becke–Johnson damping scheme (D3BJ), and solvation was accounted for using the conductor-like screening model (COSMO) in DCM.39–41 Initial geometries were adapted from the reported crystal structures and geometry optimizations were confirmed to be minima by frequency calculations performed at the same level of theory.
In agreement with previous studies, all three isomers for 5 were found to be minima while only the COequatorial and COTBP isomers were favorable for 7 (Fig. S7†). The COaxial isomer for complex 7 consistently converged to a structure matching the COequatorial isomer (Fig. S8†). These findings support that the COaxial isomer is too sterically encumbered at large cone angles and are in agreement with experimental observations. Upon comparison of the bonding parameters from crystallography to DFT, only minor variations are observed, suggesting DFT well represents the electronic structure complexes 5 and 7 (Tables S9 and S10†). This is further confirmed by comparison of the DFT predicted FTIR with experiment (Fig. S9†). Here the DFT calculated ν(CO) bands are in excellent agreement with experiment suggesting the DFT calculated structures match those observed experimentally. This is well supported by previous studies,1 and offers further support for the presence of the described isomers in these systems.
While previous studies have proposed the isomerization shown in Fig. 1, where isomerization is mediated by the trigonal bipyramidal intermediate, the direct isomerization between COaxial and COequatorial could not be ruled out.1 Owing to the scope of the current investigation, this question is partly answered through superposition of the DFT structures. When the Ru, CO, and P atoms of each isomer are superimposed, one isomerization pathway is revealed. Starting from the COaxial isomer shown in orange, structural rearrangement between the COaxial (orange) and COequatorial (violet) occurs by a slight, 45° twist of the dithietene ligand.
From there, isomerization between the COequatorial (violet) and COTBP (green) isomers occurs by a second dithietene twist (<35°), which opens a coordination site on the ruthenium atom that was formerly occupied by the dithietene sulfur atom. This newly opened coordination site allows the phosphine ligand to “slide” into place forming the COTBP isomer (green). This suggests that the COTBP isomer cannot be accessed directly from the COaxial isomer without a complete dithietene twist. While the thermochemical data initially suggests that the COTBP is acting as an intermediate, this is only one slice along a multi-dimensional reaction pathway, and all three isomers may in fact be interconverting directly. The magnitude required for the second dithietene twist (COequatorial to COTBP) decreases as phosphine cone angle is increased. This is most clearly seen in complex 7 (Fig. 5, right), where very little movement of the dithietene ligand is needed. Consistent with previous findings,1 several low-frequency normal modes (ν < 1000 cm−1) have been identified for complexes 5 and 7 with nuclear displacements in alignment with these proposed rearrangement pathways.
Footnote |
† Electronic supplementary information (ESI) available: A collection of additional NMR, FTIR, DFT, and crystallographic data. See DOI: 10.1039/c9sc02359c |
This journal is © The Royal Society of Chemistry 2019 |