Koosha Ehtiatiab,
Ilya Anufrievac,
Christian Friebed,
Ivan A. Volodinab,
Christian Stolzeab,
Simon Muenchab,
Grit Festagac,
Ivo Nischangabde,
Martin D. Hagerabcd and
Ulrich S. Schubert*abcd
aLaboratory of Organic and Macromolecular Chemistry (IOMC), Friedrich Schiller University Jena, Humboldtstrasse 10, Jena 07743, Germany. E-mail: ulrich.schubert@uni-jena.de
bCenter for Energy and Environmental Chemistry Jena (CEEC Jena), Friedrich Schiller University Jena, Philosophenweg 7a, Jena 07743, Germany
cJena Center for Soft Matter (JCSM), Friedrich Schiller University Jena, Philosophenweg 7, Jena 07743, Germany
dHelmholtz Institute for Polymers in Energy Applications Jena (HIPOLE Jena), Lessingstr. 12 – 14, 07743 Jena, Germany
eHelmholtz-Zentrum Berlin für Materialien und Energie GmbH (HZB), Hahn-Meitner-Platz 1, 14109 Berlin, Germany
First published on 18th October 2024
Application of redox-active polymers (RAPs) in redox flow batteries (RFBs) can potentially reduce the stack cost through substitution of costly ion-exchange membranes by cheap size-exclusion membranes. However, intermolecular interactions of polymer molecules, i.e., entanglements, particularly in concentrated solutions, result in relatively high electrolyte viscosities. Furthermore, the large size and limited mobility of polymers lead to slow diffusion and more sluggish heterogeneous electron transfer rates compared to quickly diffusing small molecules. Although a number of RAPs with varying electrolyte viscosities have been reported in the literature, the relation between the RAP structure and the hydrodynamic properties has not been thoroughly investigated. Herein, hyperbranched 2,2,6,6-tetramethylpiperidinyloxyl (TEMPO)-based polymers intended for application as low-viscosity catholytes for RFBs are presented and the influence of the structure and the molar mass distribution on the hydrodynamic properties is investigated. A new synthesis approach for TEMPO-based polymers is established based on step-growth polymerization of a TEMPO-containing monomer using an aza-Michael addition followed by a postpolymerization modification to improve solubility in aqueous solutions. The compact structure of hyperbranched polymers was demonstrated using size-exclusion chromatography (SEC) with viscometric detection and the optimum molar mass was found based on the results of viscometric and crossover investigations. The resulting RAP revealed a viscosity of around 21 mPas at a concentration corresponding to around 1 M TEMPO-containing units, according to the calculated mass of the repeating unit, showing potential for high capacity polymer-based catholytes for RFBs. Nevertheless, possible partial deactivation of TEMPO units lowered the active TEMPO concentration of the hyperbranched RAPs. A faster diffusion and higher charge transfer rate were observed for the hyperbranched polymer compared to the previously reported linear polymers. However, in RFB tests, a poor performance was observed, which is attributed to the side reactions of the oxidized TEMPO moieties. Finally, pathways for overcoming the main remaining challenges, i.e., high loss of material during dialysis as an indication of being prone to crossover, the partial deactivation of TEMPO moieties, and the subsequent side reactions under battery conditions, are suggested.
One common approach to reduce the stack cost is to increase power density of individual cells, through reduction of the internal resistance and/or increasing the cell voltage.4–6
Alternatively, substitution of expensive cell parts has been examined. With regard to this, applying redox-active polymers (RAPs) instead of small molecules has been demonstrated to enable the substitution of expensive ion-exchange membranes by relatively inexpensive size-exclusion membranes via introduction of the size-exclusion mechanism for crossover prevention.7–9 The size-exclusion mechanism can also be achieved through application of micellar structures, slurries, and microgels.10–13
On the downside, concentrated polymer solutions are highly viscous at high concentrations and the use of lower polymer concentrations limits the battery energy efficiency and the electrolyte capacity.7,8,14,15 The first reported polymer-based RFB comprised a poly(2,2,6,6-tetramethylpiperidinyloxyl) (TEMPO)-based catholyte in which a TEMPO-containing monomer was copolymerized with a cationic comonomer to improve solubility in aqueous solutions.8 The applicable capacity of this polymer solution was 10 A h L−1 with a viscosity of 17 mPas. Later, Hagemann et al. and Fu and Zhang et al. were able to reduce the electrolyte viscosity via copolymerization with a zwitterionic comonomer, and terpolymerization using anionic comonomers, respectively.7,14 However, the relation between the polymer structure and hydrodynamic properties remains under-investigated.
The viscosity of polymer solutions is highly dependent on the polymer concentration, chain length, and salt concentration (in case of polyelectrolytes).16–18 Scaling relations for an uncharged polymer describes that while the specific viscosity scales with polymer concentration with scaling power of 1.25 in the unentangled semi-dilute regime, the scaling becomes much stronger after entanglements are formed with scaling power of 3.75 in the entangled semi-dilute regime. As a consequence, the viscosity of highly concentrated polymer solutions is strongly related to their ability to form entanglements. It is known that hyperbranched polymers possess a more compact structure compared to linear polymers at the same molar mass.19–22 As a result of such compact structure, hyperbranched polymers are less prone to form entanglements and, thus, enabling a weaker dependency of viscosity on the polymer concentration than that of linear polymers.23–25
We herein report hyperbranched TEMPO-based polymers as potentially low-viscosity catholytes for RFBs aiming to improve the electrolyte capacity. We employ an oxygen-tolerant synthesis of hyperbranched TEMPO-based polymers and thoroughly investigate the hydrodynamic properties of the resulting hyperbranched polymers with varying molar masses and concentrations. Furthermore, the electrochemical properties of the hyperbranched polymers are investigated and their potential application in RFBs is assessed.
Time limitation was due to the sensitivity of the 100 to 500 Da dialysis tubes which were already very swollen after 6 h of dialysis against deionized water. In the case of HPT-I0 which was dialyzed for around 15 minutes longer, a damage occurred to the top of one of the two dialysis tubes causing uncertainties with the amount of obtained material. For this practical limitation, only the 1000 Da tube was chosen to dialyze the HPT-II1 sample. In this case after a round of dialysis against NaCl solution and deionized water, a trace amount of iodide was observed in the sample (detected by CV, for more details see the ESI S6†) and further dialysis against NaCl solution and deionized water was applied to completely remove the trace of iodide. The polymer products were finally obtained by freeze drying. The masses of the obtained polymer (data of Fig. 5a) were obtained as the ratio of (mass of obtained polymer)/(expected polymer mass with Cl−) where the expected polymer mass with Cl− is calculated based on the mass of the polymer used before dialysis (with I−) and the iodide content of each sample (measured by elemental analysis).
Nuclear magnetic resonance (NMR) was performed on a Bruker Avance I (300 MHz) spectrometer. Elemental analysis (halogen analysis) was measured with a TLalpha20 (Si Analytics GmbH, Germany) titrator.26 UV absorption of the polymer samples were obtained using a PerkinElmer Lambda 750 UV-Vis spectrometer (1 cm quartz cuvette). The radical content of the TEMPO-based polymers was measured via electron paramagnetic resonance spectroscopy (EPR). X-band EPR spectra were acquired on an EMXmicro CW-EPR spectrometer from Bruker, Germany. The SpinCountQ software module was used for the determination of the spin count. A known PTMA polymer (radical content of 80%, determined through redox titration) was used as a reference. The radical contents of the TEMPO-containing compounds were determined from the mean values derived from the EPR spectra of three samples per compound.
The molar mass values reported in Fig. 1 and 2 are based on light scattering using the dn/dc and the instrument factor. The apparent contraction factor values were obtained as the ratio of the intrinsic viscosity of the hyperbranched polymers relative to the intrinsic viscosity of linear TEMPO-based polymer (PTMA) at the same molar mass:
(1) |
Fig. 1 Cumulative molar mass distribution of the synthesized hyperbranched polymers based on SEC-MALLS. |
Polymer sample | AEP:TEMPO | Mwa (g mol−1) | [η]ga (cm3 g−1) | αKMHS a | k a | Ms,fb (g mol−1) | [η]c (cm3 g−1) | kH c | vd (cm3 g−1) |
---|---|---|---|---|---|---|---|---|---|
a Measured using SEC with viscometric detection.b Measured using AUC and based on the intrinsic viscosity values obtained from viscometry using capillary ball combination.c Obtained from viscometry using capillary ball combination.d Obtained from densimetry. | |||||||||
HPT-I0 | 0 | 3100 ± 700 | 3.2 ± 0.1 | 0.24 ± 0.05 | 0.55 ± 0.23 | 1100 | 4.8 | 1.7 | 0.74 |
HPT-I1 | 0.1:1 | 4100 ± 300 | 3.7 ± 0.1 | 0.19 ± 0.04 | 0.83 ± 0.24 | 1700 | 6.3 | 1.0 | 0.73 |
HPT-I2 | 0.2:1 | 7600 ± 900 | 3.9 ± 0.3 | 0.23 ± 0.05 | 0.57 ± 0.27 | 2900 | 6.1 | 1.5 | 0.72 |
HPT-I3 | 0.3:1 | 18300 ± 2400 | 4.4 ± 0.2 | 0.28 ± 0.04 | 0.32 ± 0.09 | 11300 | 10.2 | 0.9 | 0.73 |
A density meter DMA 4500 M (Anton Paar, Graz, Austria) was used for the density measurements of the polymer solutions at 20 °C. The partial specific volume of the polymers was determined according to the procedure described elsewhere.28
CV was conducted by scanning from 0.2 V up to 1.1 V versus the reference electrode for six consecutive cycles with scanning rates from 5 mV s−1 up to 1000 mV s−1. Voltammetry experiments were conducted with a random order of scan rates to avoid systematic errors of electrode surface changes on the scan rate dependency and the solution was stirred shortly between the tests to ensure the same boundary conditions for every test. Determination of the peak-to-peak potential and the peak current was based on the first cycle. To determine the diffusion coefficient, the range of scan rates was selected where the peak-to-peak potential was independent of the scan rate. The diffusion coefficient for a one-electron reaction was obtained based on the eqn (2):30
(2) |
CA was conducted by applying a step potential of 1.0 V for 60 s. Before application of this step potential the electrode was held at the potential of 0.2 V for 60 s. Another period of 60 s at 0.2 V was applied after the step potential to ensure all the oxidized molecules were reduced again so the bulk concentration remains constant. The test was conducted three times, and the solution was shortly stirred in between. The diffusion coefficient was calculated based on the Cottrell equation for a semi-infinite diffusion:30
(3) |
(4) |
(5) |
The kinetic currents were then plotted as a function of overpotential (here defined as the difference between the potential and the formal potential) to determine the heterogenous charge transfer constant (k0) and the transfer coefficient (α) as described for a one-step one-electron reaction:
(6) |
As a model system, cyclohexylamine was first employed instead of 4-amino TEMPO, due to its availability and facility of NMR analysis compared to 4-amino TEMPO, to study the polymerization kinetics and molar mass variation over time. Methanol was used as a protic solvent to increase the rate of aza-Michael addition. Since acrylates are more reactive than acrylamides in aza-Michael addition, first a triacrylate (trimethylolpropane triacrylate) was used as B3.
Around 50% conversion of NH groups was obtained quickly but the reaction slowly progressed further reaching a conversion of around 90% after around 50 h corresponding to a theoretical degree of polymerization of only five (ESI S1†). This is in line with the higher reactivity of the primary amines compared to formed secondary amines in an aza-Michael addition. Longer polymer chains were aimed by having polymerization times of around four days which resulted in a conversion of 97% corresponding to a degree of polymerization around 16. However, transesterification of ester groups and subsequent shortening of polymer chains were observed as a drawback of ester groups in the backbone (see ESI S2†).
As an alternative, a triacrylamide (1,3,5-triacryloylhexahydro-1,3,5-triazine) was employed, which is expected to be more stable in polymerization reaction and less prone to hydrolysis under conditions of the final application compared to acrylates. In case of triacrylamide the reaction mixture was slightly heated to 40 °C to improve the solubility of the reactants and to increase the reaction rate of the formed secondary amines. In such conditions, around 90% conversion of the NH groups was already obtained after around 1 h of reaction (ESI S3†). Following the kinetics of the reaction, the initial polymers were obtained by variation of the polymerization time between around 4 to around 8 h (ESI S4†). At certain times, polymerization was terminated by addition of an excess amount of 1-methylpiperazine to consume all the remaining double bonds through aza-Michael addition (Scheme 1). With increasing polymerization time, the molar mass distribution became broader with an increasing fraction of long chains (Mw > 40000 g mol−1). This behavior is expected for an A2 + B3 chemistry because longer chains, possessing higher number of functionalities, have higher probability for further growth compared to shorter chains leading to broadening of the molar mass distribution as the polymerization progresses. Since very long chains can potentially increase the viscosity of the polymer solutions which is undesirable, a short polymerization time of 4 h was chosen for the TEMPO-based samples in this study.
Polymerization of 4-amino TEMPO with 1,3,5-triacryloylhexahydro-1,3,5-triazine was performed under similar conditions. However, in this case molar masses of the obtained polymers were lower with narrower distributions compared to the cyclohexylamine polymerization (ESI S4†). Further attempts to increase molar mass by varying the stoichiometry, temperature, and time did not lead to an increase indicating that the molar mass was limited by the monomer functionality. To this end, various amounts of a trifunctional amine (AEP) were added to increase the overall amine functionality and, thus, to increase the molar mass. Fig. 1 and Table 1 show the properties of the polymers obtained by this strategy. The incorporation of both monomers and the chemical structure was supported by 1H NMR spectroscopy and UV-Vis absorption spectroscopy (ESI S5†).
After the polymerization, the polymers were further modified through methylation using iodomethane, as shown in Scheme 1 followed by ion exchange to substitute I− by Cl−. Methylation and subsequent ion exchange were performed to improve the solubility of the redox-active polymers in aqueous solutions and to avoid any interference with the redox activity of the TEMPO moiety as tertiary amines and I− ions can be oxidized within the potential window of TEMPO oxidation (see ESI S6† for more details). The methylation reaction was verified qualitatively by comparing the 1H NMR spectra of the polymers before and after amine quaternization and a high conversion (>90%) was quantitatively obtained using elemental analysis and acid titration (see ESI S7† for more details). Finally, the exchange of counterions from I− to Cl− was confirmed by cyclic voltammetry.
Using SEC coupled with viscometric detection enabled estimation of the intrinsic viscosity, [η], of each molar mass elution slice, Mi, of the disperse polymer population. The relationship of [η] and Mi can be described by an exponential Kuhn–Mark–Houwink–Sakurada (KMHS) scaling relationship, . A double logarithmic plot of [η] against Mi from SEC with viscometric detection enables the estimation of the exponent αKMHS. Generally, for linear polymer systems it assumes a value of 0.5 or larger. For ideal hyperbranched systems [η] is independent of the molar mass.38 Practically, for branched or hyperbranched systems, αKMHS assumes values below 0.5.39,40
The acquired viscosity data of the hyperbranched polymers appear noisy, as observed in Fig. 2. This can be attributed to the compact nature of the hyperbranched polymers, low molar mass, relatively broad distribution, and overall low [η]-values. From the molar mass distribution of polymer samples, the middle 80% fraction of each sample is selected for investigating the KMHS relation. The resulting values of the exponent for the newly synthesized polymers determined over a limited molar mass range of around one order of magnitude are in the range of 0.19 to 0.28, indicating more compact polymer architectures compared to the linear polymer. Furthermore, the apparent contraction factor is obtained from the ratio of the intrinsic viscosity of the hyperbranched to that of the linear polymer at the same molar masses (calculated based on the obtained KMHS relation for linear PTMA), as shown in Fig. 2b.19,40 The contraction factor reflects the compactness of the hyperbranched polymers compared to linear chains with the same chemical structures.19,41 Here, due to synthesis and solubility limitations, a linear polymer with the same chemical structure was not synthesized and linear PTMA, containing TEMPO moieties but possessing a different backbone than that of the hyperbranched polymers, was used to compare the intrinsic viscosities. Hence, we have used the term “apparent contraction factor” to distinguish the ratio of intrinsic viscosities in our case from the contraction factor, which should ideally be obtained from polymers with the same chemical structure. For a wide range of molar masses, from around 1000 to 30000 g mol−1, the apparent contraction factor of all the polymer samples is between 0.8 to 0.3 and decreases with increasing molar mass, indicating the compactness of the hyperbranched polymers as compared to PTMA with similar molar mass. The much more compact structure of the hyperbranched polymers manifests itself at larger molar masses; at small molar masses the [η] values tend to converge such that there is no particular benefit in the low molar mass region in terms of intrinsic viscosity.
While SEC with viscometric detection could verify the more compact structure of polymers, although over a small and limited molar mass range, there are technical issues associated with the reported molar mass distribution. First, the intensity of the light scattering signal is low even at a rather high concentration of polymer samples used here (10 mg mL−1) which is attributed to the low molar mass of the polymer samples and the low value of refractive index increment. Additionally, the functional groups from the polymer backbone can interact with the column material interfering with the size-exclusion mechanism for separation of various polymer fractions. Such effect has previously been observed for branched polymers with functional groups with larger molecules being more likely to interact with the column.42 Here we have added 1% v/v of DEEDA to the eluent to screen such interactions. While addition of the amine enables proper elution of the sample, the polymer–column interactions cannot be completely excluded. We utilized another hydrodynamic technique known as analytical ultracentrifugation (AUC) to independently determine the molar masses, Ms,f. We conducted sedimentation velocity experiments and sedimentation-diffusion analysis (Table 1 and S3†). More details can be found in the ESI S8.† The resulting values of the molar masses, Ms,f, of the individual polymer populations are smaller than those from SEC-MALLS, an aspect that aligns with observations made with nonrelated polymer populations and when comparing SEC-MALLS and AUC.43 The trend of molar masses based on AUC follows the amount of AEP added in the polymerization. It is also confirmed that by addition of AEP, the polymer samples become more disperse. While the majority of chains are still small, larger chains are formed by addition of AEP.
Fig. 3 Dynamic solution viscosity of the synthesized hyperbranched polymers as a function of polymer concentration in ethanol measured with a capillary ball combination. |
To investigate the influence of the molar mass on the intrinsic viscosity and the entanglement formation, reduced viscosity () and specific viscosity (ηsp = ηr − 1), where ηr is the ratio of the solution viscosity to the solvent viscosity known as relative viscosity and cp is the polymer concentration, are plotted in Fig. 4. The Huggins extrapolations were performed to determine the values of intrinsic viscosity, [η], in a range of relative viscosities of ηr = 1.2–2.5, which corresponds to the dilute solution regime (Fig. 4a). The Huggins parameters for HPT-I0 to HPT-I3 are in the range of 0.9 to 1.7, which corresponds to previously observed values for hyperbranched polymer systems.38
Overall, the intrinsic viscosity values appear to be higher compared to those obtained from SEC with viscometric detection. This difference is attributed to the way those values are obtained. The global intrinsic viscosity obtained from SEC with viscometric detection is an average value of the intrinsic viscosity of polymer fractions of various molar masses, while viscometry is conducted on the entire sample solution in which the larger macromolecules have a larger influence on the viscosity than the smaller ones.
Comparing the highest and the lowest molar mass samples, it is observed that while the molar mass has increased with a factor of 6 or 10 (derived from SEC results and AUC results, respectively), the intrinsic viscosity has increased with a factor of around 2. This moderate dependency of the intrinsic viscosity on the molar mass is expected for hyperbranched polymer solutions as was discussed in the previous section with respect to the KMHS analysis. By plotting the specific viscosity as a function of polymer concentration in a double logarithmic scale, as shown in Fig. 4b, information about entanglement concentration of the polymer solutions is obtained. The dashed lines show a slope of 1.25 corresponding to the description of scaling relations for unentangled semi-dilute solutions. A positive deviation from this slope can be considered as an indication of the entanglement formation as the scaling power is higher in the entangled semi-dilute solution. While such deviation has already started from the concentration of around 11 g dL−1 for HPT-I3, the specific viscosities of HPT-I1 and HPT-I2 start to deviate at the concentration of around 15 g dL−1 and HPT-I0 shows such deviation starting from 20 g dL−1. Alternatively, the degree of dilution of a polymer solution is the product c[η], also known as the Debye parameter.44 This parameter describes the volume of the solution that is occupied by the polymer. Values of cp[η] < 1 indicate that the polymers occupy a smaller volume than the volume of the solution. A value of cp[η] of 1 is the polymer concentration at which the polymer occupies the entire solution volume. Apparently, this will result in similar concentration values as obtained from the double-logarithmic plots shown in Fig. 4b, not requiring knowledge of [η]. Values of cp[η] = 1 are found for polymer HPT-I3 being cp = 9 g dL−1, for HPT-I2 and HPT-I1 cp = 16 g dL−1, and for HPT-I0 cp = 21 g dL−1. Those values appear to correlate with the ones above and derived from the double logarithmic plots in Fig. 4b.
Overall, it can be concluded that increasing the molar mass only moderately influences the intrinsic viscosity of hyperbranched polymers due to their compact structure. The viscosity of highly concentrated solutions is, however, still significantly affected by the molar mass. This can be due to a not ideally hyperbranched molecular architecture and the higher chance of entanglement formation for larger macromolecules resulting in a significant influence of the molar mass on the overall viscosity of entangled semi-dilute solutions.
As was discussed in the synthesis section, in the last step an ion-exchange process was performed on the polymer samples using dialysis in a concentrated NaCl solution. After these processes, the samples corresponding to various molar masses are named as HPT-IxCl with x = 0 to 3. The size exclusion of the polymer samples was evaluated by measuring the amount of obtained material after dialysis. Fig. 5a shows the results of the obtained materials for HPT-I0Cl to HPT-I3Cl using size-exclusion membranes with MWCOs of lower values (100 to 500 Da) and higher values (1000 Da).
First, in all the cases, a significant amount of material (at least 22%) is lost after two days of dialysis. Considering the molar mass distribution from SEC and AUC (Fig. 1 and S13†), this fact can be understood based on the presence of a large fraction of low molar mass polymers (or oligomers) in all the samples. Results obtained from AUC reveal that even HPT-I3 with the highest average molar mass, has a significant amount of material with small sedimentation coefficients (Fig. S13†). The smaller species can pass the dialysis membrane, rationalizing the high amount of material loss in all cases. Using the membrane with the lower MWCO is shown to significantly improve the amount of obtained material (by 22% on average), indicating the presence of a significant fraction of polymers with sizes in the range of the pore size of the membrane with a MWCO of 1000 Da. Regarding the effect of molar mass, by comparing HPT-I0 with HPT-I1, a significant increase in the amount of obtained polymer is observed with introduction of the trifunctional amine (29% and 56% for HPT-I0 compared to 51% and 78% for HPT-I1 using membranes with MWCO of 1000 Da and 100 to 500 Da, respectively). However, with further addition of trifunctional amine in the synthesis (going from HPT-I1 to HPT-I3) the amount of obtained material shows insignificant changes.
Besides, dynamic viscosities of polymer solutions obtained after ion exchange and dialysis using the membrane with MWCO of 1000 Da was measured in 1.5 M NaCl solutions to evaluate final applicability of the hyperbranched polymers in the redox flow battery electrolyte, as shown in Fig. 5b.
HPT-I2Cl and HPT-I3Cl were highly soluble in 1.5 M NaCl aqueous solutions enabling viscosity measurements in the entire range of concentrations. For HPT-I1Cl, however, small particles were found floating in the solution which were not dissolved with further dilution. While these particles did not impede the movement of the ball in the larger capillary of the viscometer, the data of which are shown in Fig. 5b, they interfered with the moving ball in case of the smaller capillary. The latter has a smaller gap between ball and capillary walls. In case of HPT-I0Cl (which was exceptionally obtained from dialysis with the membrane with MWCO of 100 to 500 Da) similar particles were observed with a higher content and, thus, HPT-I0Cl was considered insoluble and the corresponding viscosity data are not presented.
With increasing molar mass from HPT-I1Cl to HPT-I3Cl, the increase of viscosity with concentration becomes steeper, similar to the previously observed trend of viscosity for polymers in the neutral state. Overall, it can be concluded that using a trifunctional amine for chain extension has improved the solubility and suppressed crossover by comparing HPT-I1Cl and HPT-I0Cl. However, further increase in molar mass, from HPT-I1Cl to HPT-I3Cl, does not significantly impact crossover while causing a significant increase in solution viscosity. In fact, additional trifunctional amine molecules have a higher chance to react with large molecules possessing relatively high number of reactive groups. As a result of such reactivity, larger molecules grow even faster compared to the smaller molecules by addition of trifunctional amine causing a broader molar mass distribution with increased trifunctional amine content. Such behavior is observed in the molar mass distributions obtained by both SEC and sedimentation coefficient distribution from AUC measurements. As a result of such broad distributions, a small fraction of large macromolecules is formed which dominate the viscosity behavior through introduction of macromolecular entanglements. On the other hand, the major fraction of small molecules is not significantly affected by additional amounts of trifunctional amine. This fraction can transport through the membrane causing a high material loss via crossover.
Based on these results, HPT-I1Cl is chosen as the optimum point of viscosity and crossover behavior. A new sample (HPT-II1) was synthesized with a similar amount of trifunctional amine but terminated with 4-amino TEMPO instead of 1-methylpiperazine which resulted in an increased TEMPO content in each polymer chain compared to HPT-I1 (ESI S9†). For HPT-II1Cl, the ion exchange and dialysis were continued until no trace of iodide was observed in the sample (based on cyclic voltammetry). With extension of the dialysis time even more material was lost, which resulted in only about 10% of the obtained material. HPT-II1Cl revealed excellent solubility in 1.5 M aqueous NaCl. For comparison, examples of linear TEMPO-based polymers from literature are also plotted (cyan squares correspond to datapoints of a linear zwitterionic copolymer reported by Hagemann et al. and cyan triangles correspond to those of a linear anionic terpolymer reported by Fu and Zhang et al.).7,14 These data are calculated based on the concentrations provided in the respective publication (in A h L−1), the mass of the repeating unit for the provided structure, and the reported degree of oxidation. These two examples are selected among other polymer samples in the literature for multiple reasons. First, to our best knowledge, they represent the lowest reported viscosity values for TEMPO-based polymers which are soluble in aqueous solutions. Second, the selected supporting electrolyte solution (1.5 M NaCl) and the utilized method of viscosity measurement are similar to our case, enabling quantitative comparison. However, from the viewpoint of electrostatic interactions the selected examples are different. While HPT-II1Cl is a cationic polymer, the copolymer reported by Hagemann et al. is zwitterionic and the terpolymer reported by Fu and Zhang et al. is anionic. In the dilute region (1.2 ≤ ηr ≤ 2.5), the viscosity behavior of all new systems is comparable to the linear polymers reported earlier. With the increase of the concentration, differences in the viscosity become more prominent.
The terpolymer shows the lowest viscosities among the polymers revealing a viscosity of 8.1 mPas at a concentration of around 40 g dL−1. Next, the hyperbranched polymer HPT-II1Cl reveals a superior viscosity behavior compared to the zwitterionic copolymer. While the zwitterionic copolymer shows a viscosity of 20 mPas at a mass concentration of 29 g dL−1, HPT-II1Cl shows a lower viscosity (14 mPas) at a higher concentration of 34 g dL−1 and a slightly higher viscosity of 21 mPas at a mass concentration of 38 g dL−1. Furthermore, the molar mass of the repeating unit (one TEMPO monomer + one comonomer) in case of the linear zwitterionic copolymer is around 520 g mol−1 while in the case of HPT-II1Cl the calculated mass containing 1 mole of TEMPO is around 385 g mol−1 (for details of the calculations see ESI S10†) which is expected to further improve the applicable TEMPO concentration and corresponding electrolyte capacity compared to the previously reported linear polymer.
However, the mass of the polymer samples equivalent to 1 mole TEMPO was also measured experimentally using EPR and compared with the calculated values, as shown in Fig. 6. In all cases the mass obtained from the EPR was higher than that of the calculated one based on the polymer repeating units by a factor of 1.6 to 2.4. One source of error can arise from the calculation of the incorporated AEP. In the calculations it is assumed that each AEP molecule reacts with three acrylamide groups but if the reactivity of the formed secondary amine is too low to consume the third acrylamide, the remaining double bond will be consumed in the termination step with 1-methylpiperazine or 4-amino TEMPO. In that case, the overall mass of polymer for each TEMPO unit would increase, in case of HPT-I1Cl to HPT-I3Cl, compared to the currently calculated value. However, we argue that this cannot be the only source of this difference because even in the sample without AEP, i.e., HPT-I0Cl, where the repeating unit is exactly known, the experimental mass equivalent to one active TEMPO is higher than the mass of the repeating unit by a factor of 1.6. This can be explained by deactivation of the TEMPO moieties through a side reaction. To further investigate the cause of deactivation, the monomer (4-amino TEMPO) and the polymer HPT-II1 (before methylation) were also evaluated using EPR. The mass of 224 ± 32 g mol−1 corresponding to a ratio of EPR-based mass to the calculated mass of 1.31 was obtained for the monomer. The deviation from 1 can be explained by a combination of error margin in the EPR measurement and the presence of impurities in the monomer. For HPT-II1 before methylation a mass of 385 ± 19 g mol−1 corresponding to a ratio of 1.18 was obtained. This indicates that the major deactivation had not occurred during the polymerization as the TEMPO moieties are active after polymerization. One could, therefore, hypothesize that side reactions during the methylation reaction could potentially deactivate the TEMPO moieties. To this end, for future works, one direction of improvements can be to optimize the methylation reaction.
Fig. 6 Calculated (blue) and EPR-based (red) mass of polymer samples equivalent to 1 mole active TEMPO. |
Based on the experimentally obtained equivalent masses, the viscosity of polymer solutions is also plotted as a function of active TEMPO molarity in Fig. 7a. In this case a superior behavior, i.e., significantly higher accessible concentration range within the relevant viscosity window, of HPT-II1Cl is observed in comparison to HPT-I1Cl, HPT-I2Cl, and HPT-I3Cl. The highest concentration of HPT-II1Cl used herein is 0.69 M corresponding to a theoretical capacity of 18.5 A h L−1 which is superior compared to previously reported capacities of TEMPO-based polymer catholytes (ESI S16†). However, at a lower concentration, the viscosity value of HPT-II1Cl appears between those of the zwitterionic copolymer and the anionic terpolymer, i.e., at a concentration of around 0.56 M active TEMPO (15 A h L−1), viscosities of 20 mPas, 14 mPas, and 8 mPas were observed for the zwitterionic copolymer, HPT-II1Cl, and the anionic terpolymer, respectively. Overall, it can be concluded that partial deactivation of TEMPO moieties limited the viscosity advantage of HPT-II1Cl by decreasing the concentration of active TEMPOs.
Fig. 7 (a) Dynamic solution viscosity as a function of molarity of active TEMPO. Molarity of active TEMPO is based on the experimentally determined masses from Fig. 6. (b) Dynamic solution viscosity as a function of calculated molarity of TEMPO. Calculated molarity assumes that all the incorporated TEMPO units are active. For hyperbranched polymer, it corresponds to the calculated value of molar mass equivalent to 1 TEMPO in Fig. 6 and in case of the zwitterionic copolymer and the anionic terpolymer such value corresponds to an oxidation degree of 100%. |
Fig. 7b displays an alternate view of such viscosity plots, in which it is assumed that all TEMPO moieties would be active. For HPT-II1Cl such case corresponds to the calculated molar masses equivalent to 1 mole TEMPO as in Fig. 6 (blue bar), while in cases of the zwitterionic copolymer and the anionic terpolymer it corresponds to an oxidation degree of 100%. By this view, we can reduce the discussion on pure effects of the viscosity behavior. In that case, HPT-II1Cl reveals similar properties to that of the anionic terpolymer, and a significantly improved behavior compared to the zwitterionic copolymer. At a calculated molarity of TEMPO of 0.57 mol L−1 the solution viscosity reached 20 mPas for the linear zwitterionic copolymer, while a solution viscosity of 21 mPas is obtained at a TEMPO molarity of 0.98 mol L−1 for the hyperbranched polymer system. In other words, if both polymer systems would contain 100% activity of TEMPO, the branched system allows for approximately 70% more active TEMPO at a dynamic viscosity value of around 20 mPas (Fig. 7b).
Fig. 8 displays the behavior of HPT-II1Cl in cyclic voltammetry at a scan rate of 100 mV s−1 for six consecutive cycles. Here, to prevent any interference of the remaining tertiary amine groups in the electrochemical activity, we have conducted cyclic voltammetry at pH = 2.0 to ensure that any remaining tertiary amine groups of the polymer are protonated and electrochemically inactive. To prevent the influence of migration on the electrochemical behavior of HPT-II1Cl, an electrolyte concentration of 1.5 M NaCl was used similar to the targeted solution for viscosity and battery experiments. Furthermore, a standard potential of E0 = 0.72 V vs. Ag/AgCl was obtained, which is within the range of other TEMPO derivatives (examples are provided in the ESI S14†). The redox potential of a redox-active moiety can be shifted by the effect of electron-donating or electron-withdrawing groups. In most previously reported TEMPO polymers for redox flow batteries, a methacrylate-based monomer was used for the polymerization in which the TEMPO moiety is attached to the backbone via an ester group.7,8,14 Two of them reported an E0 value of 0.7 vs. Ag/AgCl, which is slightly lower compared to our case.7,8 This is attributed to the difference in electron-withdrawing effect of the group on the 4-position of the TEMPO ring (quaternary ammonium in HPT-II1Cl compared to the ester group in the methacrylate-based polymers). N,N,N-2,2,6,6-Heptamethylpiperidinyloxy-4-ammonium chloride (TMATEMPO) small molecule with a quaternary ammonium reveals a higher redox potential of E0 = 0.79 V vs. Ag/AgCl.45 A peak potential separation (ΔE) of 72 mV and an anodic to cathodic peak current ratio of represents an electrochemically quasi-reversible behavior at this scan rate.
Moreover, the properties of HPT-II1Cl are evaluated as a function of cycling scan rate from 5 mV s−1 to 1000 mV s−1, as shown in Fig. 8b. ΔE is plotted as a function of scan rates in Fig. 8c. Here, ΔE fluctuates around a value of 80 mV at scan rates up to 100 mV s−1. However, with further increase ΔE increases reaching up to around 200 mV at a scan rate of 750 mV s−1. In Fig. 8d, the peak oxidation current is plotted as a function of square root of the scan rate. If electrochemical reversibility is assumed, eqn (2) can be applied to this relation from which a diffusion coefficient of D = 1.6 × 10−6 cm2 s−1 is estimated. This value is three times lower than the diffusion coefficient reported for the TMATEMPO small molecule (4.8 × 10−6 cm2 s−1), but higher than the previously reported values for other TEMPO-based polymers (see the ESI S15† for comparison).
This faster diffusion arises from the same origin as the discussed viscosity advantage, namely the more compact structure of hyperbranched polymers and the lower molar mass compared to the previously reported linear polymer. Herein, it should be noted that the diffusion coefficient obtained based on the electrochemical methods are derived based on the theories which are mainly developed for molecules with one redox-active moiety. In the case of polymers, each molecule has multiple redox-active moieties and, thus, when the molecule is transported through diffusion to the electrode surface, there are multiple redox-active moieties that can react at the electrode surface. Therefore, the diffusion coefficient obtained here is not necessarily the diffusion coefficient of polymer molecules but the result of polymer chain diffusion and multiple reactions at the surface. In that sense the reported value is an apparent diffusion coefficient reflecting the transport of each redox-active moiety.
However, considering the higher ΔE of the polymer (around 80 mV) compared to a reversible process (around 59 mV), the behavior of HPT-II1Cl may deviate from eqn (2). Thus, the diffusion coefficient obtained herein is later compared to those values based on other electrochemical experiments (CA and RDE), which do not rely on electrochemical reversibility. Fig. 9 displays the result of potential step chronoamperometry for HPT-II1Cl at the similar concentration of 5.7 mM in 1.5 M NaCl solution at pH = 2.0, where different colors show the results of the repetitions of the same experiment. Application of a step potential of 1.0 V starts the oxidation of the TEMPO moieties in the polymers at the electrode surface, which is described by the Cottrell equation (eqn (3)). The measured current at the electrode surface decreases over time due to an increasing thickness of the diffusion layer. In case of semi-infinite diffusion, the current scales inversely with the square root of time. Under most experimental conditions, natural convection limits the thickness of the diffusion layer and, thus, a deviation from the semi-infinite diffusion is observed for long step durations. Anyhow, as observed in Fig. 9b, a nearly linear behavior is observed in the time range from 1 to 60 seconds. Therefore, a linear regression is built for this range and from its slope, the diffusion coefficient is estimated as D = 1.2 × 10−6 cm2 s−1, which is close to the value obtained from cyclic voltammetry.
Comparing diffusion coefficients obtained from CV and RDE, it is worth mentioning that application of eqn (2) to CV data requires electrochemical reversibility in the investigated scan rates, while with RDE such conditions are not required. Furthermore, in our experiments, a blank subtraction was performed for RDE data evaluation, while for CV, the current of the blank solution was not subtracted, which might slightly influence the peak currents due to the presence of capacitive currents.
Furthermore, the inverse current is plotted as a function of inverse square root of rotation rates, which shows a linear relationship, following the Koutecký–Levich equation (eqn (5)), as shown in Fig. 10c. From the intercept of this linear relation, the mass-transfer-independent kinetic currents (Ik) can be obtained at various potentials with respect to standard potential. Finally, from the logarithmic plot of Ik versus the potential difference from the standard potential, as shown in Fig. 10d, the electrochemical reaction rate constant (k0 = 1.9 × 10−3 cm s−1) and the transfer coefficient (α = 0.71) are obtained. This value of the electrochemical rate constant is lower than that of the small molecule TMATEMPO (k0 = 4.2 × 10−3 cm s−1) but higher than those of the previously reported TEMPO polymers (see ESI S15†) including the zwitterionic copolymer (9.7 × 10−4 cm s−1) and the anionic terpolymer (7.8 × 10−4 cm s−1).7,14,45 The fact that polymers reveal a lower heterogenous charge transfer rate than that of small molecules has been attributed to their limited Brownian motion.46 In addition, it was reported that long and flexible alkyl chains could improve the charge transfer rate. Based on that, the higher charge transfer rate of HPT-II1Cl compared to methacrylate-based linear polymers can similarly be attributed to the higher mobility of the TEMPO moieties in HPT-II1Cl, which are attached to the backbone via three sp3 atoms before the carbonyl group, compared to methacrylate-based polymers where TEMPO is directly attached to the ester group (1 sp3 atom between the TEMPO moiety and the carbonyl group).
Overall, using cyclic voltammetry, chronoamperometry, and rotating disc electrode voltammetry, a quasi-reversible oxidation/reduction for the TEMPO-containing hyperbranched polymer with a higher diffusion coefficient and a faster heterogenous charge transfer rate compared to previously reported TEMPO-based linear polymers is demonstrated.
The first scenario is the overreduction of the TEMPO moieties through side reactions in the methylation reaction. As was discussed based on EPR results, the TEMPO moieties were partially deactivated, and the deactivation could be due to overreduction through reaction with iodide in the methylation reaction. To mitigate this issue, a new batch of polymer, similar to HPT-II1, was synthesized. This time chloromethane was utilized for the methylation reaction instead of iodomethane. The resulting methylated polymer was characterized using EPR from which a mass equivalent to 1 mole TEMPO of 557 ± 33 g mol−1 was obtained, which is slightly lower compared to HPT-II1Cl. A battery test under similar conditions was performed utilizing this polymer (results are presented in the ESI, S13†). However, even this polymer revealed a similar behavior with an overcharging in the first cycle followed by a low discharge capacity and a fast capacity decay like the HPT-II1Cl material. Therefore, one can conclude that possible side reactions of iodide during the methylation reaction were not the main cause of the overoxidation of the TEMPO moieties under battery conditions.
The second possible cause of the overcharging could be oxidation of the remaining tertiary amines in the polymer chain after methylation. Although the methylation reaction was quite efficient (as described in ESI S7† for HPT-I0 to HPT-I3), a fraction of tertiary amines is expected to remain. Those can undergo irreversible oxidation. As discussed in the ESI S6,† the deprotonated tertiary amines of the polymer start to get oxidized at potentials higher than 0.8 V vs. Ag/AgCl, while the protonation at low pH can prevent such side reactions. Although the battery tests were started at pH = 2.0 to prevent such oxidation, the pH value of the system could increase during cycling if hydrogen evolution reaction had consumed the hydronium ions on the anolyte side. Sampling of the catholyte after cycling revealed the pH value still being 2.0 and, thus, the possibility of deprotonation and oxidation of the tertiary amines can be ruled out.
The third possibility is the interaction of the polymer backbone with the oxidized TEMPO. Previous reports of small molecules and polymers containing TEMPO and amide moieties have shown reasonable stability and, thus, the amide groups of the polymer backbone are not expected to interfere with the redox activity of the TEMPO moieties.48,49 Therefore, it is expected that other chemical interactions with, for example, protonated amine groups, triazine rings, etc. are responsible for such TEMPO degradation. Understanding the underlying mechanism requires further investigations.
Another possible parasitic reaction can be the disproportionation of TEMPO moieties in acidic medium.50 However, comparing the voltammogram of the TEMPO polymer kept at a pH = 2.0 without cycling to a similar solution after cycling shows that cycling was the main cause of the TEMPO degradation but not keeping the TEMPO polymer at pH = 2.0 (ESI S13†).
Overall, the battery performance of the hyperbranched TEMPO polymer indicates the presence of a side reaction which degrades the TEMPO moieties during cycling. Here we analyzed three possible scenarios, from which we excluded two: (i) overreduction of the TEMPO moieties by iodide during methylation and (ii) interference of the oxidation of the remaining tertiary amine groups. However, to understand the mechanism of such degradation, designing symmetrical battery experiments with further chemical characterizations are suggested which is out of the scope of this work.
Additionally, potential applicability of the polymers with size-exclusion membranes was evaluated by measuring the amount of material crossover in ex situ dialysis experiments using membranes with various MWCOs. Overall, a high amount of material loss was observed (between around 20 to 70% w/w depending on the sample molar mass and the membrane MWCO) after only two days of dialysis followed by an even higher loss when dialysis was continued longer. Therefore, material crossover remained as one of the main challenges of the hyperbranched polymers, which is mainly attributed to the broad distribution of the molar mass in the synthesized polymers. Interestingly, a significant decrease in crossover was obtained through introduction of the trifunctional amine in polymerizations. Based on the results of crossover and viscosity, an optimum amount of added trifunctional amine was found and a hyperbranched polymer with the selected amount of trifunctional amine and additionally incorporated TEMPO moieties was synthesized.
Although the selected polymer possessed an overall high concentration of TEMPO with relatively moderate viscosities (around 1 M TEMPO, according to the calculated mass of the repeating unit, at a viscosity of 21 mPas), possible partial deactivation of TEMPO moieties decreased the concentration of active TEMPO moieties leading to a viscosity behavior in the range of previously reported linear polymers.
The electrochemical properties of the selected polymer were investigated via cyclic voltammetry, chronoamperometry, and rotating disc electrode voltammetry. An electrochemically quasi-reversible redox activity with faster diffusion and higher heterogenous charge transfer rate compared to previously reported linear polymers was observed. In the galvanostatic battery cycling test, however, the resulting polymer revealed a poor performance with an irreversible charging in the first cycle, followed by a fast-decaying capacity and low coulombic efficiencies, which are mainly attributed to the side reactions of the oxidized TEMPO moieties.
To sum up, a new synthesis and thorough hydrodynamic and electrochemical characterization of hyperbranched TEMPO-based polymers were presented, which enabled decreased viscosity and faster diffusion and charge transfer rate compared to linear polymers. Nevertheless, the main challenges of being prone to crossover, partial deactivation of the TEMPO moieties, and subsequent side reactions represent remaining challenges with regard to their application in redox flow batteries.
In future studies, these challenges can be addressed through various strategies. To mitigate the crossover problem, one idea can be to choose multifunctional amines (e.g., diethylenetriamine) instead of trifunctional amine to further increase the molar mass at the same molar ratio. Alternatively, multifunctional alkyl halides (e.g., 1,3-dichloropropane) can be used to further increase molar mass during the methylation reaction by linking the polymer molecules together on their amine groups. Deactivation of TEMPO molecules should also be further investigated. Alternative methylation procedures using, e.g., dimethyl carbonate can be considered to further optimize the methylation reaction.
Footnote |
† Electronic supplementary information (ESI) available. See DOI: https://doi.org/10.1039/d4ra03925d |
This journal is © The Royal Society of Chemistry 2024 |