Open Access Article
Phakpoom Angpanitcharoen,
Jessica V. Lamb,
Jean-Charles Buffet
,
Zoë R. Turner
and
Dermot O'Hare
*
Chemistry Research Laboratory, Department of Chemistry, University of Oxford, OX1 3TA, UK. E-mail: Dermot.ohare@chem.ox.ac.uk
First published on 19th March 2021
A new family of zirconocene complexes of the type (3-RInd#)2ZrX2 (where Ind# = C6Me5H and R = Me, Et and Ph) have been synthesised and fully characterised. Six new crystal structures have been reported (meso-(3-EtInd#)2ZrBr2, rac-(3-EtInd#)2ZrCl2, rac-(3-EtInd#)2Zr(CH2Ph)2, meso-(3-EtInd#)2Zr(CH2Ph)2, meso-(3-MeInd#)2ZrBr2 and meso-(3-MeInd#)2Zr(CH2Ph)2). The complexes were studied for slurry-phase ethylene polymerisation when immobilised on solid polymethylaluminoxane (sMAO). Variation in the initiation group was found to have greater influence over polymerisation activity for meso-catalysts than rac-catalysts, with meso-alkyl catalysts showing higher polymerisation activities than meso-halide. Below 70 °C, polymerisation activity follows the order sMAO-meso-(3-EtInd#)2Zr(CH2Ph)2, sMAO-meso-(3-EtInd#)2ZrCl2 and sMAO-meso-(3-EtInd#)2ZrBr2 (activities of 657, 561, and 452 kgPE molM−1 h−1 bar−1, respectively). sMAO-meso-(3-EtInd#)2ZrBr2 produces HDPE with the highest molecular weight, followed by sMAO-meso-(3-EtInd#)2ZrCl2 and sMAO-meso-(3-EtInd#)2Zr(CH2Ph)2 (Mw of 503, 406, and 345 kg mol−1, respectively, at 50 °C). sMAO-meso-(3-MeInd#)2ZrBr2 produced HDPE with almost identical molecular weights to sMAO-meso-(3-EtInd#)2ZrCl2 (395 kg mol−1 at 50 °C).
The discovery of ‘single-site’ metallocene-based catalysts has revived α-olefin polymerisation chemistry as these catalyst systems enable the production of (co)polyolefins with tuneable molecular structures, stereochemistries and molecular weight distributions (MWD, Mw/Mn).7–14 Self-immobilisation of catalyst and co-catalysts for olefin polymerisation has shown that the high excess of MAO can reduced by more than 90%.15 Yet these catalysts are rarely used for industrial processes, with a few exceptions including Dowlex (Dow) and Sclairtech (Nova Chem. Corp.); due to the incompatibility with the existing gas or slurry phase olefin polymerisation processes. By immobilising metallocene catalysts onto support materials, the processability problem is solved and the advantages of metallocene catalysts, including narrow MWD, high activity, and precise control over polymer microstructure, are preserved to a great extent. Supporting metallocene complexes also allows for the reduction of the ratio of MAO to metallocene,16 and have an dramatic effect on the tacticity when polymerising α-olefins.8 Many materials were tested as olefin polymerisation supports including inorganic solids such as silica,17 silicate clays,18 MgCl2,19 and polymeric aluminoxanes.20 Various functionalised polymers have also been investigated as supports. Silica is the most commonly used support for heterogeneous metallocene catalysts. Precontacting the silica surface with MAO was reported to increase catalyst loadings.21–24 The use of layered double hydroxides (LDHs), a class of anionic clays consisting of positively charged Brucite-like layers with weakly bound anions intercalated between them, as catalyst supports for olefin polymerisation has been recently reported by the O'Hare group.25–27 Heat treatment of MAO for a prolonged period was reported to increase polymerisation activity of up to 25%. This is believed to be caused by the aggregation of MAO to form larger particles which are insoluble in hydrocarbon solvents. This insoluble ‘solid MAO’ can act simultaneously as both co-catalyst and support in slurry phase ethylene polymerisation.28,29
:
50 mixture of rac- and meso-(3-PhInd#)2ZrCl2 (1) as an orange yellow solid in 45% yield (Scheme 1). The two sets of resonances corresponding to each isomer are indistinguishable in the 1H NMR spectrum; resonances between 6.90 and 7.90 ppm correspond to the phenyl protons, while two singlets at 5.86 and 6.52 ppm correspond to the cyclopentadienyl protons and singlets between 1.60 and 2.80 ppm correspond to the methyl groups on the indenyl rings (Fig. S11†).
3-EtInd#Li was synthesised according to a previously reported procedure.30 Two equivalents of 3-EtInd#Li were reacted with one equivalent ZrBr4 to afford an orange solid comprising of a 40
:
60 mixture of rac- and meso-(3-EtInd#)2ZrBr2 after work-up. Recrystallisation of the isomeric mixture yielded meso-(3-EtInd#)2ZrBr2 (meso-2) as orange crystals in 6% yield. Similar to 1, meso-2 shows diagnostic resonances corresponding to the cyclopentadienyl proton and methyl groups of the indenyl rings. meso-2 also shows two doublet of quartets at 2.69 and 3.22 ppm and a triplet at 1.05 ppm corresponding to the diastereotopic methylene and methyl protons of the ethyl groups respectively (Fig. S13†).
rac-(3-EtInd#)2ZrCl2 (rac-3) and meso-(3-EtInd#)2ZrCl2 (meso-3) were prepared according to a literature procedure.30 They were reacted with two equivalents KCH2Ph to afford rac-(3-EtInd#)2Zr(CH2Ph)2 (rac-4) and meso-(3-EtInd#)2Zr(CH2Ph)2 (meso-4) as a yellow solids, both in 93% yield after work-up. The C2-symmetry of rac-4 results in magnetically equivalent indenyl rings and benzyl groups, as evidenced in the 1H NMR spectrum; the phenyl protons are seen as a doublet, doublet of doublets, and triplet at 6.62, 7.08, and 6.80 ppm, respectively, while two doublets at 0.67 and 0.87 ppm correspond to the diastereotopic benzylic protons (Fig. S15†). Similar to meso-Me2SB(3-EtI*)Zr(CH2Ph)2,31 the symmetry of meso-4 results in magnetically inequivalent benzyl groups in the 1H NMR spectrum; the ortho, meta, para, and benzylic protons of one benzyl group are observed as a doublet at 6.46 ppm, a doublet of doublets at 7.04 ppm, a triplet at 6.75 ppm, and a singlet at 0.24 ppm, while for the other benzyl group they appear at 6.58, 7.06, 6.79, and 0.79 ppm (Fig. S17†). 3-MeInd#Li was synthesised according to a previously reported procedure.32 Two equivalents 3-MeInd#Li were reacted with one equivalent ZrBr4 to afford an orange solid comprising of a 20
:
80 mixture of rac- and meso-(3-MeInd#)2ZrBr2, and impurities. Recrystallisation of the isomeric mixture in toluene at −30 °C yielded meso-(3-MeInd#)2ZrBr2 (meso-5) as orange crystals in 6% yield. The 1H NMR spectrum of meso-5 shows that the two indenyl rings are magnetically equivalent, with diagnostic resonances corresponding to the cyclopentadienyl protons and indenyl methyl groups (Fig. S19†). One equivalent meso-5 was further reacted with two equivalents of KCH2Ph to afford meso-(3-MeInd#)2Zr(CH2Ph)2 (meso-6) as a yellow solid in 88% yield. Similar to meso-4, the Cs-symmetry of meso-6 results in magnetically inequivalent benzyl groups in the 1H NMR spectrum (Fig. S21†).
Orange crystals of meso-2, rac-4, and meso-4 were grown from pentane at room temperature, while crystals of rac-3, meso-5 and meso-6 were grown from a toluene solution at −30 °C and room temperature. The solid-state molecular structures are depicted in Fig. 1, with selected bond lengths and angles presented in Table 1.
| Complex | Zr-X | Zr-Cpcent | α | δ | TA | |
|---|---|---|---|---|---|---|
| a X = Br, Cl or CH2Ph; α, δ and TA: crystallographic parameters are defining in Fig. S23. | ||||||
| meso-2 | meso-(3-EtInd#)2ZrBr2 | 2.6041(5) | 2.2425(11) | 52.44(14) | 131.93(4) | 81.49 |
| 2.5927(5) | 2.2430(11) | |||||
| rac-3 | rac-(3-EtInd#)2ZrCl2 | 2.4183(6) | 2.2453(8) | 51.72(10) | 132.44(3) | 173.14 |
| 2.4514(6) | 2.2374(8) | |||||
| meso-4 | meso-(3-EtInd#)2Zr(CH2Ph)2 | 2.307(2) | 2.3045(8) | 51.80(11) | 132.99(3) | 102.95 |
| 2.3284(19) | 2.2666(9) | |||||
| rac-4 | rac-(3-EtInd#)2Zr(CH2Ph)2 | 2.307(7) | 2.322(3) | 51.9(4) | 133.63(11) | 157.30 |
| 2.299(7) | 2.281(3) | |||||
| meso-5 | meso-(3-MeInd#)2ZrBr2 | 2.6133(4) | 2.2500(10) | 54.35(13) | 133.62(4) | 35.46 |
| 2.5745(4) | 2.2525(10) | |||||
| meso-6 | meso-(3-MeInd#)2Zr(CH2Ph)2 | 2.3302(18) | 2.2763(9) | 53.66(10) | 133.92(3) | 112.56 |
| 2.3132(18) | 2.3419(10) | |||||
For meso-2, the average Zr–Br bond length (2.598 Å) and Zr–Cpcent distance (2.243 Å) are within the range of those reported for other unbridged dibromobisindenyl zirconocene complexes.16–18 The ring tilt angle (α) (52.44°) and angle about the metal centre (δ) (131.93) are similar to meso-(3-MeInd#)2ZrCl2 (53.28° and 132.02°, respectively).32 The torsion angle (TA) measures the twist of the two indenyl rings relative to one another (Fig. S23†). The introduction of ethyl groups causes increased steric strain, resulting in meso-2 adopting ‘gauche’ conformation, which is reflected in the TA angle of 81.49° compared to 37.16° for meso-(3-MeInd#)2ZrCl2.32
The average Zr–Cl bond length and Zr–Cpcent distance of rac-3 (2.435 and 2.241 Å) are in good agreement with the 3-methyl analogue rac-(3-MeInd#)2ZrCl2 (2.421 and 2.228 Å).32 The α and δ values (Fig. S23†) of rac-3 (51.72° and 132.44°) are comparable to rac-(3-MeInd#)2ZrCl2 (51.58° and 133.59°) and rac-(3-EtInd#)2HfCl2 (51.80° and 132.56°), as are the torsion angles (173.14, 171.31, and 173.01°, respectively), implying that the introduction of 3-ethyl groups has little-to-no impact on the molecular frameworks. Structure comparisons of rac-3 with the unmethylated analog, rac-(1-EtInd)2ZrCl2, show similar average Zr–Cl bond lengths (2.435 and 2.442 Å) and Zr–Cpcent distances (2.241 and 2.233 Å).33 The α and δ values for rac-3 (51.72° and 132.44°) are slightly larger than for rac-(1-EtInd)2ZrCl2 (49.22° and 131.98°). However, the TA value for rac-3 (173.14°) is significantly larger than for rac-(1-EtInd)2ZrCl2 (120.02°) due to the increased sterics of methylation. For meso-4, the average Zr–CH2Ph bond length (2.318 Å) is similar to rac-4 (2.303 Å) and the hexamethylated analogue rac-(3-MeInd#)2Zr(CH2Ph)2 (2.320 Å),32 and is slightly longer than (4,7-FInd)2Zr(CH2Ph)2 (2.294 Å).34 Similar to meso-2, meso-4 has a TA value of 102.95°, adopting ‘gauche’ conformation to minimise the steric clash between the ethyl and benzyl groups. The average Zr–Br bond length and Zr–Cpcent distances of meso-5 (2.594 and 2.251 Å) are similar to meso-2 (2.598 and 2.243 Å). The α value for meso-5 (54.35°) is slightly larger than meso-(3-MeInd#)2ZrCl2 (53.28°) due to the larger Br ligands.32 The TA value of 35.46° for meso-5 is comparable to meso-(3-MeInd#)2ZrCl2 (37.16°),32 however, is significantly less than meso-2 (81.49°) due to the decrease in steric strain on exchanging ethyl for methyl. For meso-6, the average Zr–CH2Ph bond length (2.322 Å) is comparable to rac-4 (2.302 Å) and rac-(3-MeInd#)2Zr(CH2Ph) (2.320 Å).32 However, comparisons between the α values of meso-6 and rac-4 are challenging because the complexes adopt different conformations, reflected in the TA values of 112.56° and 157.30°, in order to minimise the steric interactions between the indenyl and benzyl ligands.
:
50 mixture of rac- and meso-(3-PhInd#)ZrCl2 (1), meso-(3-EtInd#)2ZrBr2 (meso-2), rac-(3-EtInd#)2ZrCl2 (rac-3), meso-(3-EtInd#)2ZrCl2 (meso-3), rac-(3-EtInd#)2Zr(CH2Ph)2 (rac-4), meso-(3-EtInd#)2Zr(CH2Ph)2 (meso-4) meso-(3-MeInd#)2ZrBr2 (meso-5sMAO) were immobilised on solid polymethylaluminoxane (sMAO)28,29,35 according to a literature procedure with an initial aluminium to zirconium ([AlsMAO]0/[Zr]0) loading of 200.30 Slurry-phase ethylene polymerisations were conducted in 150 mL ampoules with 50 mL hexane, 10 mg pre-catalyst and 2 bar ethylene at 50–90 °C for 30 minutes with triisobutylaluminium (TiBA, Al(CH2CH(CH3)2)3) scavenger with [Alscav]0/[M]0 of 1000 (Fig. 2).
1sMAO shows significantly lower polymerisation activities than rac-3sMAO and meso-3sMAO, likely due to the large steric bulk of the phenyl groups. rac-3sMAO displays higher polymerisation activities than meso-3sMAO below 70 °C, and almost identical polymerisation activities above 70 °C, while rac-4sMAO and meso-4sMAO display similar polymerisation activities across the entire temperature range. rac-3 shows significantly higher polymerisation activities when immobilised on sMAO compared to MAO modified silica (activities of 561 and approximately 75 kgPE molM−1 h−1 bar−1 respectively at 70 °C).30 rac-3sMAO and rac-4sMAO show lower polymerisation activities than the 3-methyl analogues (under similar polymerisation conditions with [AlsMAO]0/[Zr]0 = 300); activities of 858 and 870 kgPE molM−1 h−1 bar−1 at 60 °C for rac-3sMAO and sMAO-rac-(3-MeInd#)2ZrCl2 respectively, and activities of 864 and 1063 kgPE molM−1 h−1 bar−1 at 60 °C for rac-4sMAO and sMAO-rac-(3-MeInd#)2Zr(CH2Ph)2 respectively.35 However, meso-3sMAO shows polymerisation activities approximately double sMAO-meso-(3-MeInd#)2ZrCl2 (651 and 343 kgPE molM−1 h−1 bar−1 respectively at 60 °C).35
Looking at the meso-catalysts, below 70 °C polymerisation activity follows the order meso-4sMAO, meso-3sMAO, and meso-2sMAO (activities of 657, 561, and 452 kgPE molM−1 h−1 bar−1, respectively). This is likely due to faster formation of the active species, as catalysts containing halide initiation groups first need to be alkylated.36,37 meso-4sMAO shows less resilience to higher polymerisation temperatures and, as a result, the polymerisation activities converge to approximately 250 kgPE molM−1 h−1 bar−1 at 90 °C.
Reducing the steric bulk in the indenyl ligand leads to increases in polymerisation activities; the 3-methyl catalyst meso-5sMAO shows higher polymerisation activities than the corresponding 3-ethyl catalyst meso-2sMAO (maximum activities of 927 kgPE molM−1 h−1 bar−1 at 70 °C and 503 kgPE molM−1 h−1 bar−1 at 50 °C, respectively). Unlike the other catalysts, which show continual decreases in activity with increasing temperature, meso-5sMAO shows a peak in activity at 70 °C (927 kgPE molM−1 h−1 bar−1) which could be related to the fact that meso-5sMAO possesses the largest gap aperture (α = 54.35(13)°).38,39 meso-5sMAO displays lower polymerisation activity than the ethylene bridged analogue sMAO-meso-(EBI*)ZrCl2; activities of 874 and 1331 kgPE molM−1 h−1 bar−1 respectively at 60 °C (under similar polymerisation conditions with [Al]0/[Zr]0 = 300),35 and lower than solution phase of well-known metallocene complexes.40
The molecular weights (Mw) and molecular weights distribution (Mw/Mn) of the polyethylenes produced were analysed by gel permeation chromatography (GPC) (Fig. 3). rac-3sMAO and rac-4sMAO produce polymers with similar molecular weights; 391 and 381 kg mol−1 at 50 °C, which is expected due to their similar polymerisation activities. rac-3sMAO produces polymers with lower molecular weights than the dimethyl silyl bridged analog, Mw of 447 mol−1 at 50 °C for sMAO-rac-Me2SB(3-EtI*)ZrCl2. rac-3-HfsMAO produced polymers with the lowest molecular weights (297 kg mol−1 at 50 °C). For the meso-3-ethyl catalysts, meso-2sMAO produced polymers with the highest molecular weights, followed by meso-3sMAO and meso-4sMAO (Mw of 503, 406, and 345 kg mol−1, respectively, at 50 °C). meso-5sMAO produced polymers with almost identical molecular weights to meso-3sMAO (395 kg mol−1 at 50 °C). All polymers show the expected decreases in polymer molecular weights with increasing polymerisation temperature, attributed to stronger chain transfer reactions at elevated temperatures.41 No long chain branching was observed by GPC or NMR spectroscopy,42,43 but very large molecular weights distribution were observed (Mw/Mn around 5–8, Tables S4–S9†) which will point towards several catalytic species on the surface.44
When immobilised on solid polymethylaluminoxane (sMAO) these complexes produce HDPE in slurry-phase ethylene polymerisations. The supported meso-alkyl catalysts showing higher polymerisation activities than meso-halide catalysts; polymerisation activity follows the order meso-4sMAO, meso-3sMAO, and meso-2sMAO.
These metallocene complexes and their inorganic based supported catalysts show promising results for use as slurry-phase polymerisation catalysts. Further testing with hydrogen, co-monomers and α-olefins will be carried out.
:
50 mixture of rac-(3-PhInd#)2ZrCl2 (rac-1) and meso-(3-PhInd#)2ZrCl2 (meso-1) as an orange-yellow solid in 45% yield (1.15 g, 1.68 mmol). 1H NMR (chloroform-d1, 400 MHz, 298 K) δ (ppm): 7.75 (Ph-H, 2H, d, 3JHH = 7.5 Hz), 7.56 (Ph-H, 2H, d, 3JHH = 7.5 Hz), 7.28 (Ph-H, 4H, m), 7.18 (Ph-H, 6H, m), 7.14 (Ph-H, 2H, m), 6.91 (Ph-H, 2H, m), 6.88 (Ph-H, 2H, m), 6.42 (Ind-H, 2H, s), 5.76 (Ind-H, 2H, s), 2.63 (Ind-CH3, 6H, s), 2.19 (Ind-CH3, 6H, s), 2.17 (Ind-CH3, 12H, s), 2.11 (Ind-CH3, 12H, s), 2.08 (Ind-CH3, 6H, s), 1.92 (Ind-CH3, 6H, s), 1.88 (Ind-CH3, 6H, s), 1.54 (Ind-CH3, 6H, s). 13C{1H} NMR (chloroform-d1, 101 MHz, 298 K) δ (ppm): 136.1 (Ind-Ph), 136.0 (Ind-Ph), 134.9 (Ind), 134.8 (Ind), 134.0 (Ph-H), 133.7 (Ind), 133.6 (Ind), 132.5 (Ind), 132.4 (Ph-H), 132.2 (Ph-H), 132.1 (Ph-H), 131.8 (Ind), 130.7 (Ind), 130.5 (Ind), 130.5 (Ind), 127.9 (Ind), 127.7 (Ph-H), 127.6 (Ind), 127.6 (Ind), 127.3 (Ph-H), 127.3 (Ph-H), 127.1 (Ph-H), 127.0 (Ph-H), 127.0 (Ph-H), 125.8 (Ind), 125.4 (Ind), 125.1 (Ind), 124.4 (Ind), 98.6 (Ind-H), 94.4 (Ind-H), 18.8 (Ind-CH3), 18.5 (Ind-CH3), 18.3 (Ind-CH3), 17.1 (Ind-CH3), 16.9 (Ind-CH3), 16.8 (Ind-CH3), 16.7 (Ind-CH3), 16.7 (Ind-CH3), 16.6 (Ind-CH3), 15.1 (Ind-CH3). CHN analysis (%): calculated C 70.15, H 6.18; observed C 70.26, H 6.38. HRMS (ESI): expected m/z 682.1705; observed 682.1736 [M]+.
:
60 mixture of rac-(3-EtInd#)2ZrBr2 and meso-(3-EtInd#)2ZrBr2 as an orange solid. Recrystallisation of the isomeric mixture in DCM at −30 °C yielded meso-(3-EtInd#)2ZrBr2 (meso-2) as orange crystals in 6% yield (0.178 g, 0.263 mmol). Orange crystals of meso-(3-EtInd#)2ZrBr2 suitable for a single crystal X-ray diffraction study were grown in pentane at room temperature. meso-(3-EtInd#)2ZrBr2: 1H NMR (chloroform-d1, 400 MHz, 298 K) δ (ppm): 6.25 (Ind-H, 2H, s), 3.22 (Ind-CH2–CH3, 2H, dq, 2JHH = 15.0 Hz, 3JHH = 7.5 Hz), 2.69 (Ind-CH2–CH3, 2H, dq, 2JHH = 15.0 Hz, 3JHH = 7.5 Hz), 2.54 (Ind-CH3, 12H, s), 2.27 (Ind-CH3, 6H, s), 2.23 (Ind-CH3, 6H, s), 1.70 (Ind-CH3, 6H, s), 1.05 (Ind-CH2–CH3, 6H, t, 3JHH = 7.5 Hz). 13C{1H} NMR (chloroform-d1, 101 MHz, 298 K) δ (ppm): 134.2 (Ind), 133.4 (Ind), 130.9 (Ind), 130.8 (Ind), 127.8 (Ind), 126.0 (Ind), 125.6 (Ind), 125.2 (Ind), 94.6 (Ind-H), 23.0 (Ind-CH2–CH3), 17.6 (Ind-CH3), 17.1 (Ind-CH3), 16.8 (Ind-CH3), 16.6 (Ind-CH3), 15.6 (Ind-CH2–CH3), 14.6 (Ind-CH3). CHN analysis (%): calculated C 56.71, H 6.25; observed C 56.61, H 6.33. HRMS (ESI): expected m/z 674.0695; observed 674.0724 [M]+.
:
80 mixture of rac-(3-MeInd#)2ZrBr2 and meso-(3-MeInd#)2ZrBr2 and impurities. Recrystallisation of the isomeric mixture in toluene at −30 °C yielded meso-(3-MeInd#)2ZrBr2 (meso-5) as orange crystals, suitable for a single crystal X-ray diffraction study, in 6% yield (0.180 g, 0.277 mmol). 1H NMR (chloroform-d1, 400 MHz, 298 K) δ (ppm): 5.90 (Ind-H, 2H, s), 2.64 (Ind-CH3, 6H, s), 2.55 (Ind-CH3, 6H, s), 2.22 (Ind-CH3, 6H, s), 2.20 (Ind-CH3, 6H, s), 2.17 (Ind-CH3, 6H, s), 2.12 (Ind-CH3, 6H, s). 13C{1H} NMR (chloroform-d1, 101 MHz, 298 K) δ (ppm): 134.0 (Ind), 133.4 (Ind), 131.0 (Ind), 129.6 (Ind), 127.4 (Ind), 127.4 (Ind), 126.8 (Ind), 121.2 (Ind), 96.7 (Ind-H), 17.5 (Ind-CH3), 16.7 (Ind-CH3), 16.6 (Ind-CH3), 16.6 (Ind-CH3), 16.6 (Ind-CH3), 16.2 (Ind-CH3). CHN analysis (%): calculated C 55.46, H 5.90; observed C 55.50, H 6.02. HRMS (EI): expected m/z 646.0382; observed 646.0366 [M]+.
:
complex molar ratio of 200
:
1 were used unless stated otherwise. The reaction mixture was heated to 60 °C for 2 hours with occasional swirling. The coloured solid was allowed to settle from the clear, colourless solution which was decanted, and the solid was dried in vacuo to afford the product as follows:
Footnote |
| † Electronic supplementary information (ESI) available: NMR spectroscopy, X-ray crystallography, polymerisation data. CCDC 2058030–2058035, 2058038 and 2058039. For ESI and crystallographic data in CIF or other electronic format see DOI: 10.1039/d1ra01912k |
| This journal is © The Royal Society of Chemistry 2021 |