DOI:
10.1039/D0SC04967K
(Edge Article)
Chem. Sci., 2020,
11, 12436-12445
Structure and dynamics of catalytically competent but labile paramagnetic metal-hydrides: the Ti(III)-H in homogeneous olefin polymerization†
Received
8th September 2020
, Accepted 24th September 2020
First published on 24th September 2020
Abstract
Metal hydride complexes find widespread application in catalysis and their properties are often understood on the basis of the available crystal structures. However, some catalytically relevant metal hydrides are only spontaneously formed in situ, cannot be isolated in large quantities or crystallised and their structure is therefore ill defined. One such example is the paramagnetic Ti(III)-hydride involved in homogeneous Ziegler–Natta catalysis, formed upon activation of CpTi(IV)Cl3 with modified methylalumoxane (MMAO). In this contribution, through a combined use of electron paramagnetic resonance (EPR), electron-nuclear double resonance (ENDOR) and hyperfine sublevel correlation (HYSCORE) spectroscopies we identify the nature of the ligands, their bonding interaction and the extent of the spin distribution. From the data, an atomistic and electronic model is proposed, which supports the presence of a self-assembled ion pair between a cationic terminal Ti-hydride and an aluminate anion, with a hydrodynamic radius of ca. 16 Å.
Introduction
Transition metal hydrides play crucial roles in many catalytic processes of importance in biology, materials science and industrial chemistry. Indeed, M–H bonds are involved in the biosynthesis of dihydrogen by hydrogenase,1 in nitrogen fixation by MoFe cofactor2 as well as in different catalytic reactions of industrial relevance such as hydrogenation of multiple carbon–carbon and carbon–heteroatom bonds,3 hydrocyanation,4 hydroformylation,5 hydrosilylation6 and in olefin isomerisation and polymerisation.7,8 The structure and chemistry of metal-hydride species involved in homogeneous catalysis has been comprehensively reviewed by Butler and Crimmin.9 However, only recently thorough characterisation of paramagnetic systems has gained momentum.10 Focusing on titanium, Ti-hydride complexes usually belong to the class of heterobimetallic bridged hydrides with general formula Ti-H-M (M = Mg, Zn, Al, Ga), these are active in a number of relevant synthetic reactions11 and are stable enough to be crystallized. Ti/Al heterobimetallic hydrides are usually obtained by salt-metathesis from Ti(IV) precursors whereby the formation of the Ti-H-Al motif often leads to the Ti(III) paramagnetic state. Such open-shell hydrides are notoriously challenging to characterise in particular when formed in situ owing to their short-lived and labile nature, which hampers standard structural characterisation. One such example is the Ti(III)-hydride complex formed during the activation reaction of titanium half-metallocenes with methylalumoxane (MAO) which is the subject of this work. In the early 1980s, the development of single site homogenous olefin polymerization catalysts based on group 4 metallocenes or half-metallocenes activated with MAO pivoted the mechanistic understanding of various features of olefins polymerization. To date, this reaction is one of the most investigated and well-understood processes where most mechanistic details concerning the catalytic cycle are well-defined and modelled.12–14 MAO is an alkyl aluminum compound of general formula [–Al(CH3)O–]n coordinating a variable amount of AlMe3 and featuring three-dimensional cage-like cluster structures with tetrahedral aluminum centers bridged by three-coordinate oxygen atoms.15–17 In modified MAO (MMAO), bulkier isobutyl or isooctyl groups partially replace the methyl groups in the MAO structure.18,19 During the in situ activation of group 4 metallocenes, MAO/MMAO acts as an alkylating agent leading to a metal–alkyl bond, via ligand exchange reactions. Moreover, in the case of titanocenes and half-titanocenes, the Ti(IV) pre-catalysts can be readily reduced to lower oxidation states, namely Ti(III) and Ti(II), due to the lower redox potential of this metal as compared to Zr and Hf. This reduction process is thought to be at the origin of the extremely high polymerization activity of the CpTiX3 catalysts (Cp = μ5-C5H5; X= halide; alkyl) in syndiospecific styrene polymerization,20–23 chemoselective polymerization of conjugated dienes24–26 and the synthesis of stereoregular styrene-butadiene27,28 copolymers where the role of the Ti(III) species as active catalysts has been firmly assessed.29 The formation of a number of open-shell Ti(III) species upon MAO/MMAO activation was demonstrated through conventional continuous wave (CW) electron paramagnetic resonance (EPR) spectroscopy in solution.30–32 Furthermore, quantitative 1H NMR analysis of the syndiotactic polystyrene produced under quasi living polymerization conditions (polydispersity index (PDI) = 1.3–1.6) by these catalysts demonstrated that Ti(III)-H species are involved in the initiation reaction and that a stationary concentration (∼20% of the paramagnetic fraction) is produced under polymerization conditions through a reversible insertion–elimination of the growing polymer chain.33
Despite their relevance, the chemical structure of such open-shell Ti(III) hydride complexes and the origin of these species are still an open question of great interest in MAO activated half-titanocene catalysis. NMR spectroscopy has been instrumental in elucidating the activation and the structure of Ti(IV) cationic complexes and their stabilization through formation of ion pairs formed upon MAO/MMAO activators.34–37 Similar detailed spectroscopic studies of the corresponding Ti(III) open-shell complexes are completely lacking. EPR is especially suited for this task as the magnetic hyperfine couplings between the electron spin and the surrounding magnetically active nuclei encode geometrical information such as the distance between the magnetic nuclei and the electron spin-carrying centre and their relative orientation. So far only conventional X-band CW EPR techniques have been used to characterize the paramagnetic products generated upon activation of the precursors, providing key insight into the reaction mechanism.30–32 However, pulsed EPR offers a whole arsenal of specific experiments (ENDOR, ESEEM, HYSCORE) capable of investigating the local coordination environment of paramagnetic species,38–41 affording resolution comparable to that of advanced NMR techniques.
In this contribution we employ X- and Q-band EPR spectroscopy to elucidate the electronic and geometrical structure of Ti(III) species obtained by reaction of CpTiCl3 with MMAO. We provide compelling evidence for a terminal Ti-H bond and for the presence of a 27Al in the second coordination sphere of the Ti(III) cationic complex. This is the first detailed structural description of the metal hydride formed in situ upon reductive activation of a titanium half-metallocene precursor, which is stabilized by the formation of an ion pair with the alumoxane moiety.
Results and discussion
Solution X-band CW EPR
The toluene solution of the CpTiCl3 precursor is EPR silent, as expected for a diamagnetic Ti(IV) system. When such solution was treated with MMAO (see also Fig. S1 in the ESI†) under argon atmosphere, the resulting reaction mixture ([Ti] = 6.0 mM; Al/Ti molar ratio = 250; 25 °C) readily turned from a red to a brown colour. The corresponding room temperature X-band EPR CW spectra recorded at different reaction times after preparation and their relative simulations are given in Fig. 1. The 5 minute spectrum (Fig. 1a) consists of a doublet centred at giso = 1.9917 ± 0.0002 – characteristic of a Ti(III) ion – and separated by 0.65 mT. Minor satellite peaks are visible at either side of the spectrum and account for 47,49Ti (47Ti: 7.44%, I = 5/2; 49Ti: 5.41%, I = 7/2) hyperfine couplings, if a |aisoTi| = 17.5 ± 0.2 MHz is imposed in the simulation. This |aisoTi| value is consistent with literature reports on similar systems42 and corresponds to a spin density at Ti of ∼95%. Given that 1H and 35,37Cl are the only magnetically active nuclei present in significant natural abundances and based on 1H ENDOR experiments (vide infra), the doubled splitting is assigned to the hyperfine coupling of the unpaired electron with a 1H nucleus with |aisoH| = 18.16 ± 0.20 MHz, corresponding to a spin density at the proton of 1.3% (see ESI†). These values are in agreement with previous studies on related systems.32–34 When the sample is left to evolve at room temperature, the intensity of the doublet spectrum smoothly decreases over time as a second species gradually appears, while the total paramagnetic fraction stays constant. Besides the already discussed doublet, the 95 minute spectrum (Fig. 1b and S2 in the ESI†) shows this second species characterized by a giso = 1.9821 ± 0.0002 and line shape with a characteristic inflection point in the middle that can be simulated by imposing a 27Al hyperfine coupling of |aisoAl| = 6.2 ± 0.3 MHz (27Al: 100%, I = 5/2). Both these signals have been previously attributed32,43 to a Ti-hydride adduct (Ti-H hereafter) and a Ti-alkyl adduct (Ti-R) stabilized by MAO.
|
| Fig. 1 Experimental (black) and simulated (red) X-band EPR spectra of CpTiCl3 in presence of MMAO ([Ti] = 6.0 mM; Al/Ti molar ratio = 250) recorded (a) after 5 minutes and (b) after 95 minutes from the mixing of the reagents. The simulated spectrum in b is the sum of the Ti-H and Ti-R components, with relative weight 30:70. On the right hand side the X-band CW-EPR spectra recorded during a polymerization reaction performed in situ in the EPR cell are shown on the g-scale: (c) activated catalyst (CpTiCl3 in presence of MMAO, [Ti] = 6.0 mM; Al/Ti molar ratio = 250); (d) after 10 minutes of reaction with styrene; (e) after 45 minutes of reaction and heating at 50 °C. The Ti-H signal is highlighted in red while asterisks indicate extra species formed during the polymerisation reaction. Experimental conditions: temperature = 298 K, mwfrequency = 9.8 GHz, mwpower = 10 mW. Simulation parameters are summarized in Table 1. | |
One crucial point at this stage is the demonstration of the catalytic relevance of the Ti-H species. To this end, X-band CW-EPR spectra were recorded at different stages of the polymerization reaction with styrene, which was chosen as a benchmark substrate. The EPR spectrum of the activated sample is shown in Fig. 1c, it displays the characteristic doublet of the Ti-H and a small fraction of the Ti-R species. Polymerization of styrene was initiated by adding 0.3 mL of anhydrous styrene to the catalyst. In the EPR spectrum recorded after 10 minutes from the mixing (Fig. 1d) a drastic decrease of the Ti-H signal was observed as the polymer started forming in the EPR tube. After 45 min of reaction followed by heating at about 50 °C, the Ti-H signal grows again (Fig. 1e) indicating that the Ti-H species is restored under these conditions. The decrease of the Ti-H signal during the initial stages of the polymerization points to a direct involvement in the initiation reaction, whereby a styrene monomer is inserted in the Ti-H bond producing a Ti-R species. On the other hand, restoring of the Ti-H signal intensity at high polymerization temperature can be explained considering the generation of Ti-H species in the course of the polymerization reaction via β-hydride transfer from the last monomer unit of the growing polymer chain (see Scheme 2). The oscillating nature of the Ti-H signal thus firmly demonstrates that this species is actively involved in the catalytic cycle, participating in both initiation and termination reactions.33 In the following, we use pulsed EPR techniques to determine the electronic and geometrical structure of the Ti-H species.
Solid state Q-band CW EPR
An understanding of the electronic ground state and the coordination geometry of the Ti(III) ion can be obtained if the full anisotropy of the g tensor is resolved. To this end, Q-band CW spectroscopy at cryogenic temperature is advantageous as the suppression of molecular reorientation with respect to the applied magnetic field allows the resolution of the full g tensor, and the higher applied magnetic field permits greater accuracy in the determination of its principal values. Fig. 2 shows a representative experimental CW Q-band spectrum measured at 50 K and its simulation (see also Fig. S2 in the ESI†). The spectrum shows two partially overlapping signals, one with a rhombic g tensor [gx, gy, gz] = [1.9812, 1.9916, 2.0023] and the other with an almost axial g tensor [gx, gy, gz] = [1.9737, 1.9770, 1.9956]. On the basis of the giso measured at room temperature at X-band and the relationship giso = 1/3(gx + gy + gz), it can be concluded that the rhombic tensor belongs to the Ti-H signal, whereas the Ti-R species accounts for the nearly axial tensor. The largest principal g value (gz) of Ti-H corresponds, within the experimental error, to the free electron value (ge = 2.0023) unambiguously indicating a singly occupied molecular orbital (SOMO) comprised predominantly of a Ti dz2 orbital. In this picture the directions of the principal g value gz (2.0023) can be set along the cylindrical axis of the dz2 orbital, whereas the directions of gx and gy cannot be unequivocally determined at this stage.
|
| Fig. 2 Experimental (black) and simulated (red) Q-band EPR spectra of CpTiCl3 after reaction with MMAO. (a) CW spectrum recorded with experimental parameters: temperature = 50 K, mwfrequency = 33.7 GHz, mwpower = 1 mW. The simulated spectrum is the sum of the Ti-H and Ti-R components, with relative weight 35:65. (b) Two-pulse ESE-EPR spectrum recorded with experimental parameters: temperature = 50 K, mwfrequency = 33.7 GHz, interpulse delay = 200 ns. The simulated spectrum only accounts for the Ti-H component. Simulation parameters are summarized in Table 1. Arrows indicate the magnetic field settings at which Davies ENDOR and HYSCORE experiments were performed. | |
|
| Scheme 1 (a) Summary of all constraints and structural parameters derived from the analysis of the EPR experiments (for clarity, only one out of the five Cp protons is shown), dashed lines represent trough–space interaction, whereas continuous lines represent direct bonds; (b) chemical structure derived from the EPR data, dashed lines represent the ion pair interaction between the [CpTi(III)H]+ cation and the MMAO anion. The hydrodynamic diameter (dh) of the ion pair is also reported. | |
Electron Spin Echo (ESE)-detected EPR
Further information on the electronic and geometrical structure of the paramagnetic Ti(III) complexes can be obtained through pulse EPR measurements.
The Q-band ESE-detected EPR spectrum recorded at 50 K is shown in Fig. 2b. The echo detected EPR signal corresponds to the absorption of the CW EPR spectrum, however, inspection of the signal in Fig. 2b reveals that the axial component resonating at 1220 mT is missing in the ESE spectrum. This clearly indicates that such species are characterized by fast relaxation times and are filtered out in the ESE spectrum, which can be well simulated considering only the parameters relative to the Ti-H species (Fig. 2b and Table 1). The simulation accounts for 90% of the observed signal and displays a poorer agreement with the experimental spectrum only in the high field region where the major contribution of the Ti-R species would be expected. Attempts to recover the Ti-R component by lowering the temperature at 20 K were not successful. The very different relaxation dynamics of Ti-H and Ti-R allow the selective study of the former via hyperfine spectroscopy in order to define the coordination environment of the Ti(III) centre.
Table 1 Summary of the EPR parameters for Ti-H and Ti-R derived in this work, errors were estimated from spectral simulations. All hyperfine and nuclear quadrupole interactions are given in units of MHz, while angles are in degrees
|
Ti-H |
Ti-R |
g
x
, gy, gz |
[1.9812, 1.9916, 2.0023] ± 0.0002 |
[1.9737, 1.9770, 1.9956] ± 0.0002 |
g
iso
|
1.9917 ± 0.0002 |
1.9821 ± 0.0002 |
|aisoTi| |
17.5 ± 0.2 |
|
a
iso
H
|
−18.16 ± 0.10 |
|
T
H
|
[−10.16, −0.67, +10.82] ± 0.20 |
|
[α, β, γ]H |
[40, 70, 90]° ± 10° |
|
a
iso
Al
|
3.6 ± 1.0 |
6.2 ± 0.3 |
T
Al
|
[−1.7, +4.3, −2.7] ± 1.0 |
|
[α, β, γ]Al |
[0, 0, 10]° ± 5° |
|
e
2
qQ/hAl |
18 ± 2 |
|
H Q-band ENDOR
The 1H hyperfine interaction is no longer resolved in the solid state Q-band CW- and ESE-EPR spectra due to inhomogeneous line broadening. In order to recover such vital information for the structural characterisation of Ti-H, field-dependent Q-band Davies ENDOR spectra were collected at the three canonical positions of the EPR absorption envelope at 20 K, these are indicated by arrows in Fig. 2b. ENDOR spectroscopy provides the NMR spectrum of a paramagnetic centre. Pairs of transitions are expected at νENDOR = |νI ± 1/2a|, where νI is the Larmor frequency of the magnetic isotope under investigation and a is the (orientation-dependent) sought hyperfine coupling. At 33.7 GHz, with an applied field of ∼1200 mT ENDOR lines belonging to protons with a < 2νH, are expected to be centred at ∼51 MHz.
The experimental spectra and their relative simulations are collected in Fig. 3. Two sets of lines are readily discriminated: the outer features centred at about ±9 MHz are attributed to the 1H giving rise to the doublet of lines in solution; whereas inner features corresponding to interacting protons with a maximum coupling of the order of 8 MHz (±4 MHz) can be assigned to the protons of the cyclopentadienyl ring and alkyl groups in MMAO.45 The relatively low signal-to-noise ratio of the ENDOR experiments stems from the elusive nature of the Ti-H complex which spontaneously forms in solution at a low steady-state concentration.
|
| Fig. 3
1H orientation selective Davies ENDOR spectra of Ti-H recorded at Q-band and at the three canonical positions of the EPR spectrum (black traces) and their simulations (red traces) considering only the hydride ligand with A1H = [−8, −17.5, −29] and Euler angles = [40, 70, 90]°. Experimental parameters: temperature = 20 K, mwfrequency = 33.7 GHz, π pulse = 60 ns, π/2 pulse = 30 ns, τ = 400 ns, RF pulse 14 μs. In the x-axis the frequency scale gives the deviation of the resonance lines from the Larmor frequency, νH, at each field position. | |
Strongly coupled 1H.
The outer region of the ENDOR spectra, corresponding to large hyperfine couplings, can be satisfactorily simulated by one 1H hyperfine tensor with components A1H = [−8, −17.5, −29] MHz, where a negative sign is assumed based on a direct Ti-H bond model, as a result of the through-bond spin polarization mechanism, in analogy with the pπ C–H fragment.44 The best fit was obtained with a relative orientation of the A tensor with respect to the g tensor expressed by the Euler angles [α, β, γ] = [40, 70, 90]° ± 10°. A slightly less convincing fitting is achieved for the gx orientation at which a residual contribution from Ti-R may still be present. To further validate the determination of the principal components of A1H, the simulation of the ENDOR spectrum obtained as the sum of the spectra reported in Fig. 3 was also evaluated (Fig. S3 in ESI†). By removing the orientation selectivity, this procedure simplifies the ENDOR spectrum and allows to assess the overall lineshape of the ENDOR envelope.46
The 1H hyperfine tensor can be decomposed into an isotropic component aiso = −18.16 MHz, fully consistent with the value determined from CW EPR, and an anisotropic – almost fully rhombic – component T1H = [−10.16, −0.67, 10.82] MHz.
In general, since the anisotropic component of the hyperfine tensor depends on the through space electron dipole–nuclear dipole interaction it can be used to derive bond length and relative orientation between the spin bearing centre (in the present case Ti(III)) and 1H. For short distances (<2.5 Å), as in the case of a nucleus directly bound to the paramagnetic metal ion, the commonly used point-dipole approximation fails and is of no meaningful use. Moreover, the point-dipole approximation would predict an axial anisotropic tensor [T, T, −2T], which is in open contradiction with the experimentally derived (fully) rhombic tensor. A three-point dipole model has been put forward to account for the origin of 1H hyperfine anisotropy for a hydride bound to a metal ion with ndz2 ground state. The model, derived as an extension of the case of pπ C–H fragment,44 has been validated by successfully reproducing the anisotropies of the Ni–C state of the [NiFe] regulatory hydrogenase of Ralstonia eutropha1 and to the Ti(III) complex Cp2TiH3AlCTMS3 (TMS = Si(CH3)3; Cp = μ5-C5H5).42 By applying the same model to Ti-H (Scheme 1) and by setting the angle between gz and the Ti⋯H vector equal to the Euler angle β, considering a spin density on the Ti 3dz2 orbital of 0.95, the displacement of the 3dz2 lobes equal to 1.40 Å and a Ti⋯H bond length of r = 1.8 ± 0.1 Å, the experimental rhombicity derived by ENDOR spectroscopy can be reasonably reproduced Tcalc = [−10.4, −5.7, 16.1] MHz, providing an independent validation of the hyperfine components obtained by fitting.
Weakly coupled 1H.
The central region of the ENDOR spectra shows a number of unresolved proton hyperfine couplings spanning the ∼0 to ±4 MHz. The sharp central peak at νI is the matrix signal due to distant solvent protons. The two wings spanning the regions ∼±2–4 MHz are consistent with the hyperfine couplings of the Cp ring protons at ∼3 Å from Ti, as reported in previous ENDOR studies on the Ti(III) pseudo sandwich compound CpTicot.45 Moreover 1H hyperfine couplings of the order of 4.2 MHz were observed in the case of CpTiMe2 and attributed to directly ligated Ti(III)–CH3 moieties41 (Fig. S3 in the ESI†).
To summarise, since ENDOR peaks depend on the nuclear Larmor frequency, the present experiments unequivocally prove that the doublet signal originates from a 1H hyperfine coupling. Furthermore, analysis of the isotropic and anisotropic components and comparison with literature data confirm that the Ti(III) ion is coordinated by a single hydride ligand placed at ∼1.8 Å and at an angle of ∼70° with respect to the cylindrical axis of the dz2 orbital and that the hydride ligand bears 1.3% of the total spin density. Additional signals are observed and are consistent with the Cp protons at about 3 Å from the Ti centre and further alkyl groups in close proximity (Scheme 1 and Fig. S3 in the ESI†).
27Al Q-band HYSCORE
The data collected through CW EPR at X- and Q-bands and Davies ENDOR experiments support that two of the four coordinating sites surrounding the Ti(III) ion are occupied by a hydride and the Cp ligand, respectively. In order to elucidate the nature of the remaining sites, whether the chlorine ligands of the CpTiCl3 precursor are still bound to the central metal and to understand the fate of the MMAO fragment, Q-band HYSCORE spectra were recorded at the three canonical positions of the EPR envelope (indicated by arrows in Fig. 2b). HYSCORE spectra present correlation peaks only in the (+,+) quadrant centred at approximately the 27Al nuclear Larmor frequency, indicating the presence of nearby Al nuclei (Fig. 4 and full spectra in Fig. S4 in the ESI†). The maximum coupling of about 8 MHz is observed at a magnetic field position corresponding to the gy tensor component. Computer simulation (red trace in Fig. 4) indicates that the experimental pattern can be satisfactorily reproduced with a slightly rhombic hyperfine tensor A27Al = [2.0, 8.0, 1.0] ± 1.0 MHz with the Euler angles [α, β, γ] = [0, 0, 10]° ± 5° in the g tensor reference frame and a relatively large nuclear quadrupole coupling (e2qQ/h = 18 ± 2 MHz).
|
| Fig. 4
27Al orientation selective Q-band HYSCORE spectra (black traces) of Ti-H at the three canonical positions of the EPR spectrum (as indicated by arrows in Fig. 2). The simulations (red traces) correspond to a |A|27Al = [2.0, 8.0, 1.0] MHz with Euler angles = [0, 0, 10]° and nuclear quadrupole parameters e2Qq/h = 18 MHz, η = 0 with Euler angles = [0, 0, 0]°. Experimental conditions: temperature = 50 K, mwfrequency = 33.7 GHz, π/2 pulse = 16 ns, π pulse = 32 ns. Both experimental and simulated spectra are the sum of two spectra recorded with two different τ values, which were summed together after Fourier transform. The following interpulse delays were used: τ = 100 and 150 ns (gz), τ = 100 and 148 ns (gy) and τ = 110 and 148 ns (gx). | |
The latter value is in line with values reported for Ti(III)–Al heterobimetallic complexes.42 From the isotropic part of the hyperfine coupling aisoAl = 1/3(Ax + Ay + Az) = 3.6 ± 1.0 MHz and considering the value of a0 = 3367.76 MHz for unit spin density in the 27Al 3s orbital, one can calculate the electron spin population in the Al 3s orbital to be 0.1%. From the anisotropic hyperfine tensor, a lower limit for the Ti–Al distance of about 2 Å can be estimated (Scheme 1a). Both aiso and the dipolar matrix correlate with values observed in the case of VO2+–O–Al linkages in V doped aluminumphosphate47 and Ti3+–O–Al linkages in Ti doped aluminosilicate materials.48 Somewhat larger values have been reported for Ti3+ ions coordinated to an anionic aluminum hydride (alanate) ligand42 and for Ti3+ ions on the surface of chlorinated aluminum oxide.49 Finally, we note that no evidence of Cl interactions could be detected. Chlorine has two magnetically active isotopes 35Cl and 37Cl both with nuclear spin I = 3/2 and relative natural abundance of 75.77% and 24.24%, respectively. However, the HYSCORE spectra presented herein (see also Fig. S4 in the ESI†) lack any sign of magnetic coupling between the electron spin and the pristine chloride ligands, suggesting that they have been removed from the first coordination sphere of Ti. This is in stark contrast with heterogeneous Ziegler–Natta catalysts50 and Ti(III) molecular complexes,51 where intense and well defined Cl cross peaks demonstrated the direct Ti–Cl bonding interactions.
Dynamic properties of Ti-H
The A1H tensor derived from Q-band ENDOR measurements (Fig. 3), along with the g tensor obtained from the simulations of Q-band CW of the frozen solution spectra (Fig. 2a), was used to simulate the motionally averaged (fluid solution) CW spectrum at X-band, assuming an isotropic tumbling.
Considering an intrinsic linewidth 0.3 mT (≈8.4 MHz), which is compatible with the hyperfine values for weakly coupled protons (see Fig. S3 in ESI†) and aAliso derived in this work, a good fit is achieved with a rotational diffusion coefficient Dr = 4 × 108 s−1 (see Fig. S5 in ESI†). We note that the successful simulation of the motionally averaged spectrum permits to validate at once the g-components and the hyperfine parameters derived in the previous sections. Moreover, through the Stoke–Einstein–Debye equation Dr can be translated into a hydrodynamic radius (rh) of ≈16 Å (for an estimation on the uncertainty on τcorr and rh see Fig. S5 in ESI†). Given that the calculated radius is orders of magnitude larger than what it would be expected for a small metal complex, this analysis further supports the idea that the Ti-H fragment is ion paired to a negatively charged MMAO anion. The estimated radius is consistent with NMR derived values for diamagnetic Ti(IV) complexes forming ion pairs with modified MAO.52 From a methodological perspective, this approach provides thus a viable way of estimating hydrodynamic parameters for self-assembled catalytically-relevant paramagnetic complexes, which are not accessible through the more common NMR techniques (e.g. PGSE).
Electronic and geometrical structure of Ti-H
The EPR data presented show thus that activation of CpTiCl3 with MMAO in toluene leads to the formation of a structurally and electronically well-defined Ti-H complex, featuring a Ti centred SOMO with dominant dz2 character. 1H Davies ENDOR experiments allow the obtainment of compelling evidence for a 1H hyperfine coupling characterized by aiso = −18.16 MHz and a dipolar tensor T = [−10.16, −0.67, 10.82] MHz. The experimental 1H dipolar tensor can be satisfactorily reproduced adopting a three-point dipole model42,44 with a Ti-H distance of the order of r = 1.8 ± 0.1 Å and a tilting angle of ∼70° with respect to the cylindrical axis of the dz2 orbital. M–H distances in early transition metals of the first row range typically in the interval 1.55–1.77 Å,59 while longer distances are usually observed for bridging hydride species. A Ti-H distance of 1.84 Å has been reported for a terminal Ti(IV) hydride.59b
The aiso value of open-shell metal hydrides critically depends on both M–H distance and orientation, reflecting the degree of overlap between the metal and the hydride orbitals and it is found to span a wide interval ranging from 3.5 up to 132 MHz, as reported in Table 2 for selected examples. Although no simple correlation can be established, aiso values of the order of 18.5 MHz, similar to those measured here correlate with M–H distances of the order of 1.8 Å.
Table 2 EPR and structural parameters of selected open-shell metal hydride species
Compound name |
Metal |
S
|
M–H distance (Å) |
a
iso (MHz) |
T (MHz) |
Reference |
[NiFe] hydrogenase R. eutropha Ni–C state |
Ni(III)–Fe(II)ls |
1/2 |
|
−3.5 |
[21.9, −7.3, −14.5] |
1
|
α-70Ile MoFe (2H) |
|
1/2 |
|
24.3 |
[−13.3, 0.7, 12.7] |
2
|
22.3 |
[10.7, −12.3, 1.7] |
Cp2TiH3AlCTMS3 (2 equivalent H) |
Ti(III) |
1/2 |
1.892 |
−9.7 |
[16.5, −5.5, −11.0] |
42
|
[LtBu,iPr2FeH]− |
Fe(I) |
3/2 |
1.76(8) (X-ray) |
−18.5 |
[28.9, −10.7, −16.1] |
53
|
1.85 (EPR) |
{[PhBP3]Fe}2(μ-H)(μ-NH) |
Mixed valence Fe(II)–Fe(III) |
1/2 |
1.68(3) (Fe1, X-ray) |
−38.9 |
[+19.4, −17.4, −2] |
54
|
1.64(3) (Fe2, X-ray) |
FeIII(H)(N2)(thiolate) |
Fe(III) |
1/2 |
1.54(4) |
±43 |
±[−28, 13, 15] |
55
|
{N2–Fe(μ-H)2Fe–N2}+ |
Mixed valence Fe(II)–Fe(I) |
1/2 |
1.73 (average distance) |
−73 |
[11, −7, −4] |
56
|
Cp2Ti(μ-H)MgR (2 equivalent H) |
Ti(III) |
1/2 |
|
−19.6 |
— |
57
|
Co(CO)4H− |
Co(0) |
1/2 |
|
−132 |
[2.2, −1.1, −1.1] |
58
|
Ti-H |
Ti(III) |
1/2 |
1.8 (EPR) |
−18.16 |
[−10.16, −0.67, 10.82] |
This work |
In the case of Ti(III)–Al linked by two μ-H bridging hydrides with a Ti(III)-H distance of 1.892 Å, 1H couplings with aiso = 9.7 MHz were observed accompanied by a 27Al isotropic hyperfine coupling aiso = 9.4 MHz.42 Comparison to our values shows that our system is characterized by a larger 1H and a smaller 27Al hyperfine coupling, implying a reduced spin density on the Al and a larger spin density at the proton. This fact suggests a stronger Ti-H interaction with respect to bridging H.42 Even though we cannot definitely exclude an H–Al interaction, this would be rather loose, and all evidence point toward a terminal character of the Ti-H hydride. On the other hand, the aiso and T values extracted from the analysis of the 27Al HYSCORE experiments imply a Ti-L-Al linkage and an intimate interaction between the alumoxane and the transition metal.
It is generally accepted that the MAO activation of group 4 metallocenes proceeds through alkylation steps, which occurs via halide/alkyl exchange equilibria possibly up to complete substitution.34 In this work no evidence of Cl ions coordination to titanium is obtained pointing to methyl or alkyl bridges, even though oxygen coordination cannot be definitely excluded.14,18 Recent experimental and theoretical studies17,60,61 on the activation of transition metal catalysts with MAO proved that the strong Lewis acid [AlMe2]+ cation, resulting from the heterolytic dissociation of AlMe3 coordinated to the surface of the MAO cage, is responsible for the in situ generation of the alkylmetallocenium cation from the neutral dialkylmetallocene. The resulting MAO anions can in turn coordinate the metallocenium cation via methyl bridged structures.16 Our EPR results are consistent with this picture pointing to the formation of an ion pair between a MMAO anion linked via methyl bridges to a highly unsaturated [CpTi(III)H]+ cation where the H− ligand occupies a terminal position (Scheme 1b). This is at odds with a number of heterobimetallic hydride complexes used in C–H activation reactions, where a bridged TM-H-M motif (TM = transition metal, M = Al, Ga, Zn) has been established through crystallographic data.9,62 Finally, from the analysis of the motionally averaged EPR spectrum and the derived rotational diffusion coefficient, an estimate can be made of the hydrodynamic radius of the paramagnetic complex, rh ≈ 16 Å, which directly points to a MMAO supported species (ESI Fig. S5†).
Our results provide thus the first electronic and structural determination of the labile Ti-H species. They demonstrate that Ti-H is active in syndiospecific styrene polymerization and give further confirmation that the oxidation state of the titanium catalyst is +3, a paramagnetic state. A catalytic cycle involving the Ti(III)-H species can be proposed and is reported in Scheme 2. The active Ti-H species spontaneously forms upon MMAO activation (Scheme 2 top) and exhibits the characteristic two-line EPR spectrum. In presence of styrene, Ti-H undergoes monomer insertion and chain growth, where both stages are associated with the decrease of the two-line EPR signal (Fig. 1d). β-hydrogen elimination then restores the Ti-H species (Fig. 1e). This latter process may represent either a resting state or an intermediate in equilibrium with the chain growth step.
|
| Scheme 2 Proposed Ti-H catalytic cycle in the syndiospecific polymerization of styrene. | |
Conclusions
In this work, we report a thorough EPR characterization of the organometallic complex obtained by reaction of CpTiCl3 with MMAO. A combination of CW and pulse EPR/ENDOR experiments provides compelling evidence for a structurally well-defined open-shell Ti(III) with 3dz2 ground state featuring a characteristic Ti-H bond and an intimate interaction between the alumoxane and the transition metal. ENDOR and HYSCORE spectra offer unprecedented detail on the bonding interaction of the metal centre with the hydride ligand in the first coordination sphere and an Al ion on the second coordination sphere. This data can be combined to yield a structural model as the one reported in Scheme 1b pointing towards an ion pair between an unsaturated [CpTi(III)H]+ cation and a MMAO anion, with no participation of chloride ligands. It is worth noting that catalytically competent metal-hydrides are usually reactive and difficult to isolate, in this respect EPR spectroscopy and the collection of related techniques provide powerful means to define structural models. To the best of our knowledge, this is one of the few reports of advanced EPR spectroscopy applied to a labile terminal Ti-hydride complex with a well-defined catalytic function. The data presented also establish a wealth of constraints that shall guide future calculations on the reaction pathways and intermediates of the homogenous Ziegler–Natta reaction catalysed by half-titanocenes. The approach described here, with mention to the determination of dynamic parameters, may be of general applicability for structural determination of paramagnetic ion pairs.
Experimental
Materials
All manipulations of air- and/or moisture-sensitive compounds were performed under inert atmosphere (dinitrogen or argon) using standard Schlenk techniques or a MBraun drybox. MMAO-3A (7% w/w heptane solution; approximate molecular formula [(CH3)0.7(i-C4H9)0.3AlO]n and CpTiCl3 were purchased from Nouryon and Merck, respectively and used without further purification. Commercial grade toluene (Merck) was dried over anhydrous CaCl2, refluxed 48 h over Na, and distilled before use. Styrene (Merck) was stirred overnight over CaH2, distilled in vacuum and stored in the refrigerator prior to use.
EPR samples preparation
CpTiCl3 (65 mg; 300 μmole) was dissolved in toluene (10 mL); 0.10 mL of this stock solution were transferred into a glass vial and premixed with 0.40 mL of the commercial solution of MMAO to produce a toluene solution of [Ti] = 6.0 mM (Al/Ti molar ratio = 250) which was finally introduced into the EPR tube and rapidly analysed (after 5–10 min at room temperature). In the polymerization test, 0.30 mL of styrene were added in a vial containing said solution and rapidly transferred into the EPR tube and analysed.
EPR spectroscopy
All measurements were collected on a Bruker ElexSys E500 spectrometer operating either a X- or Q-band (∼9.5 and ∼34 GHz respectively). Further details for each set of experiments are given below.
CW-EPR spectroscopy
Room temperature X-band CW EPR spectra were with a super-high Q resonator (ER4122 SHQE, operating at ∼9.8 GHz). Low temperature Q-band CW EPR spectra were collected using the EN 5107D2 Bruker resonator and an Oxford CF935 cryostat.
Q-band pulsed spectroscopy
All low temperature Q-band pulsed EPR experiments were collected using the Bruker EN 5107D2 resonator housed in an Oxford CF935 cryostat and a 10 W solid-state microwave amplifier.
Electron-spin-echo (ESE) detected EPR experiments were carried out in the range 20–50 K with the pulse sequence: π/2–τ–π–τ-echo, with microwave pulse lengths π/2 = 16 ns and π = 32 ns and a τ value of 200 ns.
ENDOR measurements were performed at 20 K by employing the Davies pulse sequence (π-RF–π/2–τ–π–τ-echo), with π/2 = 30 ns, π = 60 ns and an RF pulse πRF = 14 μs. ENDOR spectra were collected in stochastic mode by randomly hopping the RF excitation frequency over a 36 MHz sweep with a resolution of 0.07 MHz (512 points).
HYSCORE spectra were recorded at 50 K using the pulse sequence π/2–τ–π/2–t1–πinv–t2–π/2–τ-echo. The pulse length for inversion pulse (πinv) and the π/2 and π pulses were 28 ns 16 and 32 ns, respectively. The first interpulse delay, τ, was adjusted at each field position either to maximise or minimise the 27Al Larmor contribution. Eight-step phase cycling was used to eliminate unwanted echoes. The experimental dataset comprised a 300 by 300 points. Time-domain spectra were baseline-corrected (third-order polynomial), apodized with a hamming window function, zero-filled to 2048 points, and fast Fourier-transformed to yield the frequency-domain spectra.
Both ENDOR and HYSCORE spectra were recorded at all three canonical positions to fully characterise the angular dependence of the coupling tensors.
Further spectrometer and experimental settings are given in the corresponding figure captions.
EPR simulations
Simulations of both CW and PULSED EPR spectra were performed using EasySpin 5.2.28 toolbox63 within the Matlab 2019b software suite (The Mathworks Inc., Natick, MA).
Conversion between magnetic field and frequency units
Much of the existing literature on Ti(III) complexes involved in Ziegler–Natta catalysis reports the hyperfine couplings in unit of Gauss (G). This unit of measurement is not a SI unit and its use is discouraged in favour of Tesla (T) or, more commonly milli-Tesla (mT). However, given that hyperfine couplings (a) express an interaction energy, Mega-Hertz (MHz) should be used. Useful conversion factors are: 1 mT = 10 G (i.e. 1 T = 10000 G) and . It is worth noting that the conversion between G and MHz is inversely proportional to the g-factor of the species considered.
Conflicts of interest
Authors declare no conflicts of interest.
Notes and references
- M. Brecht, M. van Gastel, T. Buhrke, B. Friedrich and W. Lubitz, J. Am. Chem. Soc., 2003, 125, 13075–13083 CrossRef CAS.
- R. Y. Igarashi, M. Laryukhin, P. C. Dos Santos, H. Lee, D. R. Dean, L. C. Seefeldt and B. M. Hoffman, J. Am. Chem. Soc., 2005, 127, 6231–6241 CrossRef CAS.
- H. U. Blaser, C. Malan, B. Pugin, F. Spindler, H. Steiner and M. Studer, Adv. Synth. Catal., 2003, 345, 103–151 CrossRef CAS.
- C. A. Tolman, W. C. Seidel, J. D. Druliner and P. J. Domaille, Organometallics, 1984, 3, 33–38 CrossRef CAS.
- F. Agbossou, J.-F. Carpentier and A. Mortreux, Chem. Rev., 1995, 55, 2485–2506 CrossRef.
- X. Du and Z. Huang, ACS Catal., 2017, 7, 1227–1243 CrossRef CAS.
- B. J. Burger, M. E. Thompson, W. D. Cotter and J. E. Bercaw, J. Am. Chem. Soc., 1990, 112, 1566–1577 CrossRef CAS.
- L. H. Shultz and M. Brookhart, Organometallics, 2001, 20, 3975–3982 CrossRef CAS.
- M. J. Butler and M. R. Crimmin, Chem. Commun., 2017, 53, 1348 RSC.
-
(a) H. Yang, J. Rittle, A. R. Marts, J. C. Peters and B. M. Hoffman, Inorg. Chem., 2018, 57(19), 12323–12330 CrossRef CAS;
(b) M. W. Drover, D. J. Schild, P. H. Oyala and J. C. Peters, Angew. Chem., Int. Ed., 2019, 58, 15504 CrossRef CAS;
(c) N. X. Gu, P. H. Oyala and J. C. Peters, J. Am. Chem. Soc., 2020, 142(17), 7827–7835 CrossRef CAS.
-
(a) R. Sun, J. Liu, S. Yang, M. Chen, N. Sun, H. Chen, X. Xie, X. You, S. Li and Y. Liu, Chem. Commun., 2015, 51, 6426–6429 RSC;
(b) J. Thomas, M. Klahn, A. Spannenberg and T. Beweries, Dalton Trans., 2013, 42, 14668–14672 RSC.
- H. H. Brintzinger, D. Fischer, R. Mulhaupt, B. Rieger and R. M. Waymouth, Angew. Chem., Int. Ed. Engl., 1995, 34, 1143–1170 CrossRef CAS.
- G. W. Coates, Chem. Rev., 2000, 100, 1223–1252 CrossRef CAS.
- W. Kaminsky, A. Funck and H. Hahnsen, Dalton Trans., 2009, 8803–8810 RSC.
- L. Negureanu, R. l. W. Hall, L. G. Butler and L. A. Simeral, J. Am. Chem. Soc., 2006, 128, 16816–16826 CrossRef CAS.
-
(a) H. S. Zijlstra and S. Harder, Eur. J. Inorg. Chem., 2015, 1, 19–43 CrossRef; and references therein
(b) G. Theurkauff, A. Bondon, V. Dorcet, J. F. Carpentier and E. Kirillov, Angew. Chem., Int. Ed., 2015, 54, 6343–6346 CrossRef CAS.
- H. S. Zijlstra, M. Linnolhati, S. Collins and J. S. McIndoe, Organometallics, 2017, 36, 1803–1809 CrossRef CAS.
- C. J. Harlan, S. G. Bott and A. R. Barron, J. Am. Chem. Soc., 1995, 117, 6465–6474 CrossRef CAS.
- M. Galimberti, M. Destro, O. Fusco, F. Piemontesi and I. Camurati, Macromolecules, 1999, 32, 258–263 CrossRef CAS.
- N. Ishihara, T. Seimiya, M. Kuramoto and M. Uoi, Macromolecules, 1986, 19, 2464–2465 CrossRef CAS.
- N. Hishihara, M. Kuramoto and M. Uoi, Macromolecules, 1988, 21, 3356–3360 CrossRef.
- A. Zambelli, C. Pellecchia, L. Oliva, P. Longo and A. Grassi, Makromol. Chem., 1991, 192, 223–231 CrossRef CAS.
- C. Pellecchia and A. Grassi, Top. Catal., 1999, 7, 125–132 CrossRef CAS.
- L. Oliva, P. Longo, A. Grassi, P. Ammendola and C. Pellecchia, Makromol. Chem., Rapid Commun., 1990, 11, 519–524 CrossRef CAS.
- G. Ricci, S. Italia, A. Giarrusso and L. Porri, J. Organomet. Chem., 1993, 451, 67–72 CrossRef CAS.
- A. Miyazawa, T. Kase and K. Soga, Macromolecules, 2000, 33, 2796–2800 CrossRef CAS.
- A. Zambelli, M. Caprio, A. Grassi and D. E. Bowen, Macromol. Chem. Phys., 2000, 201, 393–400 CrossRef CAS.
- M. Caprio, M. C. Serra, D. E. Bowen and A. Grassi, Macromolecules, 2002, 35, 9315–9322 CrossRef CAS.
- J. C. W. Chien, Z. Salajka and S. Dong, Macromolecules, 1992, 25, 3199–3203 CrossRef CAS.
- A. Grassi, A. Zambelli and F. Laschi, Organometallics, 1996, 15, 480–482 CrossRef CAS.
- A. Grassi, S. Saccheo, A. Zambelli and F. Laschi, Macromolecules, 1998, 31, 5588–5591 CrossRef CAS.
- L. Bonoldi, L. Abis, L. Fiocca, R. Fusco, L. Longo, F. Simone and S. Spera, J. Mol. Catal. A: Chem., 2004, 219, 47–56 CrossRef CAS.
- R. Pastorino, C. Capacchione, R. Ferro, S. Milione and A. Grassi, Macromolecules, 2009, 42, 2480–2487 CrossRef CAS.
- F. Zaccaria, L. Siana, C. Zuccaccia and A. Macchioni, Adv. Organomet. Chem., 2020, 73, 1–78 CrossRef.
- D. E. Babushkin, N. V. Semikolenova, V. A. Zakharov and E. P. Talsi, Macromol. Chem. Phys., 2000, 201, 558 CrossRef CAS.
- K. P. Bryliakov, E. P. Talsi, A. Z. Voskoboynikoc, S. J. Lancaster and M. Bochmann, Organometallics, 2008, 27, 6333–6342 CrossRef CAS.
- I. Tritto, R. Donetti, M. C. Sacchi, P. Locatelli and G. Zannoni, Macromolecules, 1997, 30, 1247–1252 CrossRef CAS.
-
S. Van Doorslaer and D. M. Murphy, EPR Spectroscopy in Catalysis, in EPR Spectroscopy. Topics in Current Chemistry 2011, ed. M. Drescher, G. Jeschke, Springer, Berlin, Heidelberg, vol. 321 Search PubMed.
- M. M. Roessler and E. Salvadori, Chem. Soc. Rev., 2018, 47, 2534–2553 RSC.
- E. Morra, E. Giamello and M. Chiesa, J. Magn. Reson., 2017, 280, 89–102 CrossRef CAS.
- S. Van Doorslaer, J. J. Shane, S. Stoll, A. Schweiger, M. Kranenburg and R. J. Meier, J. Organomet. Chem., 2001, 634, 185–192 CrossRef CAS.
- G. Rao, A. B. Altman, A. C. Brown, L. Tao, T. A. Stich, J. Arnold and R. D. Britt, Inorg. Chem., 2019, 58, 7978–7988 CrossRef.
- K. Bryliakov, N. Semikolenova, V. A. Zhakarov and E. Talsi, J. Organomet. Chem., 2003, 683, 23 CrossRef CAS.
-
W. Gordy, Theory and applications of Electron Spin Resonance, John Wiley & Sons Inc, 1980, ch 6, p. 207 Search PubMed.
- D. Gourier and E. Samuel, J. Am. Chem. Soc., 1987, 109, 4571–4578 CrossRef CAS.
- A. V. Astashkin, E. L. Klein and J. H. Enemark, J. Inorg. Biochem., 2007, 101, 1623–1629 CrossRef CAS (and references therein).
- S. Maurelli, G. Berlier, M. Chiesa, F. Musso and F. Corà, J. Phys. Chem. C, 2014, 118, 19879–19888 CrossRef CAS.
- E. Morra, S. Maurelli, M. Chiesa and E. Giamello, Top. Catal., 2015, 58, 783–795 CrossRef CAS.
- A. Piovano, K. S. Thushara, E. Morra, M. Chiesa and E. Groppo, Angew. Chem., Int. Ed., 2016, 55, 11203–11206 CrossRef CAS.
- E. Morra, E. Giamello, S. Van Doorslaer, G. Antinucci, M. D'Amore, V. Busico and M. Chiesa, Angew. Chem., Int. Ed., 2015, 54, 4857–4860 CrossRef CAS.
- S. Maurelli, E. Morra, S. Van Doorslaer, V. Busico and M. Chiesa, Phys. Chem. Chem. Phys., 2014, 16, 19625–19633 RSC.
- F. Zaccaria, C. Zuccaccia, R. Cipullo, P. H. M. Budzelaar, A. Macchioni, V. Busico and C. Ehm, ACS Catal., 2019, 9(4), 2996–3010 CrossRef CAS.
- K. P. Chiang, C. C. Scarborough, M. Horitani, N. S. Lees, K. Ding, T. R. Dugan, W. W. Brennessel, E. Bill, B. M. Hoffman and P. L. Holland, Angew. Chem., Int. Ed., 2012, 51, 3658–3662 CrossRef CAS.
- R. A. Kinney, C. T. Saouma, J. C. Peters and B. M. Hoffman, J. Am. Chem. Soc., 2012, 134, 12637–12647 CrossRef CAS.
- N. X. Gu, P. H. Oyala and J. C. Peters, J. Am. Chem. Soc., 2018, 140, 6374–6382 CrossRef CAS.
- J. Rittle, C. C. L. McCrory and J. C. Peters, J. Am. Chem. Soc., 2014, 136, 13853–13862 CrossRef CAS.
- S. I. Troyanov, V. Varga and K. Mach, J. Organomet. Chem., 1993, 461, 85–90 CrossRef CAS.
- S. A. Fairhurst, J. R. Morton and F. Preston, J. Magn. Reson., 1983, 55, 453–459 CAS.
-
(a)
R. H. Crabtree, “Hydride Complexes of the Transition Metals” in Encyclopedia of Inorganic and Bioinorganic Chemistry, John Wiley & Sons, Ltd, 2011 Search PubMed;
(b) K. Ma, W. E. Piers, Y. Gao and M. Parvez, J. Am. Chem. Soc., 2004, 126, 5668–5669 CrossRef CAS.
- T. K. Trefz, M. A. Henderson, M. Y. Wang, S. Collins and J. S. McIndoe, Organometallics, 2013, 32, 3149–3152 CrossRef CAS.
- F. Ghiotto, C. Pateraki, J. Tanskanen, J. R. Sever, N. Luehmann, A. Kusmin, J. Stellbrink, M. Linnolahti and M. Bochmann, Organometallics, 2013, 32, 3354–3362 CrossRef CAS.
- M. J. Butler, A. J. P. White and M. R. Crimmin, Organometallics, 2018, 37, 949–956 CrossRef CAS.
- S. Stoll and A. Schweiger, J. Magn. Reson., 2006, 178, 42–55 CrossRef CAS.
Footnote |
† Electronic supplementary information (ESI) available. See DOI: 10.1039/d0sc04967k |
|
This journal is © The Royal Society of Chemistry 2020 |
Click here to see how this site uses Cookies. View our privacy policy here.