Javad
Jeddi
a,
Jukka
Niskanen
b,
Benoît H.
Lessard
*cd and
Joshua
Sangoro
*a
aDepartment of Chemical and Biomolecular Engineering Ohio State University, Columbus, Ohio 43210, USA. E-mail: sangoro.1@osu.edu
bDepartment of Chemical and Metallurgical Engineering, Aalto University, Kemistintie 1, 02150 Espoo, Finland
cDepartment of Chemical and Biological Engineering, University of Ottawa, 161 Louis Pasteur, Ottawa, Ontario K1N 6N5, Canada. E-mail: benoit.lessard@uottawa.ca
dSchool of Electrical Engineering and Computer Science, University of Ottawa, 800 King Edward Ave, Ottawa, ON K1N 6N5, Canada
First published on 17th May 2024
The correlation among chemical structure, mesoscale structure, and ion transport in 1,2,3-triazole-based polymerized ionic liquids (polyILs) featuring comparable polycation and polyanion backbones is investigated by wide-angle X-ray scattering (WAXS), differential scanning calorimetry, and broadband dielectric spectroscopy (BDS). Above the glass transition temperature, Tg, higher ionic conductivity is observed in polycation polyILs compared to their polyanion counterparts, and ion conduction is enhanced by increasing the counterion volume in both polycation or polyanion polyILs. Below Tg, polyanions show lower activation energy associated with ion conduction. However, the validity of the Barton–Nakajima–Namikawa relation indicates that hopping conduction is the dominant charge transport mechanism in all the polyILs studied. While a significant transition from a Vogel–Fulcher–Tammann to Arrhenius type of thermal activation is observed below Tg, the decoupling index, often used to quantify the extent to which segmental dynamics and ion conduction are correlated, remains unaltered for the polyILs studied, suggesting that this index may not be a general parameter to characterize charge transport in polymerized ionic liquids. Furthermore, detailed analyses of the WAXS results indicate that both the mobile ion type and the structure of the pendant groups control mesoscale organization. These findings are discussed within the framework of recent models, which account for the subtle interplay between electrostatic and elastic forces in determining ion transport in polyILs. The findings demonstrate the intricate balance between the chemical structure and interactions in polyILs that determine ion conduction in this class of polymer electrolytes.
One of the peculiar characteristics of organic ionic salts such as ionic liquids and polyILs is the tendency to form distinct structural features due to the mesoscopic organization of ionic, polar, and non-polar moieties.27–30 Wide-angle X-ray scattering is instrumental in examining the structure of amorphous ionic liquids and polyILs.30,31 In the polyILs based on alkyl-imidazolium cations, different morphologies of polar and nonpolar domains can form, and the backbone-to-backbone correlation peak is observed to shift to lower q values with increasing alkyl chain length. These structural features significantly influence ion transport in polycation polyILs.22,30,32 However, it is unclear if such structural signatures exist and influence charge transport in polyILs with anions tethered on the polymer backbones.
In this study, four model 1,2,3-triazole-based polyILs with either polycation or polyanion backbones comprising comparable chemical structures, are investigated to understand the impact of mobile ion types and chemical structure of the backbone on the correlation between mesoscopic organization and ionic conductivity. A combination of broadband dielectric spectroscopy (BDS) and X-ray scattering techniques are employed to elucidate the effects of the chemistry of anionic and cationic backbone structures and ion dynamics. It is found that a subtle interplay between the chemical structure and interactions in polyILs determines the ion conduction in both classes of polyILs.
The model 1,2,3-triazole-based polymerized ionic liquids (polyILs) studied in this work include poly(1-(4-vinyl benzyl)-1H-1,2,3-triazole-4-yl)methanesulfonate) with 1-butyl-3-methyl-imidazolium (C4MIm+) and Na+ as the cations (PVB-SO3C4MIm and PVB-SO3Na), poly(4-butyl-3-methyl-1-(4-vinyl benzyl)-1H-1,2,3-triazol-3-ium chloride) (PVB-MeCl), and poly(4-butyl-3-methyl-1-(4-vinylbenzyl)-1H-1,2,3-triazol-3-ium TFSI) (PVB-MeTFSI). These systems were synthesized using a copper-catalyzed azidealkyne cycloaddition (CuAAC) ‘click’ reaction as reported previously.33 The thermal properties were characterized using differential scanning calorimetry (DSC). DSC measurements of polyIL samples were carried out using a TA Instruments DSC2500 unit calibrated with indium (temperature and enthalpy) and sapphire (heat capacity). The samples were dried in a vacuum oven at 350 K for 24 hours before sealing in the Tzero aluminum hermetic pans. Standard heat–cool–heat cycles at 10 K min−1 were applied to determine the glass transition temperatures, Tgs. The Tg was considered as the midpoint of the step in the second heating cycle.
A Novocontrol High-Resolution Alpha Dielectric Analyzer equipped with a Quatro Cryosystem for temperature control was employed for the BDS measurements in the frequency range 0.1–107 Hz. The samples were prepared by thermomechanically pressing between two Teflon sheets to obtain a film of about 10 mm in diameter. The films were sandwiched between two polished brass electrodes in a parallel plate configuration. The sample thickness was kept constant at 100 μm using silica spacers, then annealed at 350 K for 24 h in a vacuum oven to remove moisture before starting the BDS measurements. In the dielectric analyzer, the samples were further annealed at 375 K under dry nitrogen gas while the dielectric response was continuously monitored until a plateau in the real part of the complex dielectric function versus time at various frequency points, was reached. The Quatro Cryosystem was employed to control the temperature of the samples within ±0.1 K accuracy of the set point.
The X-ray scattering experiments were conducted in the transmission mode using a Xenocs Xeuss 3.0 X-ray scattering instrument equipped with a Cu Kα 50 kV Xenocs Genix ULD SL X-ray source (λ = 1.5418 Å) at room temperature. Pressed sample films (100 μm) were dried in a vacuum oven at 350 K for 24 h before starting the measurement to minimize the moisture effect on the results. The two-dimensional scattering patterns were registered in a PILATUS3 R 300K detector from Dectris, Ltd. The scattered intensity was obtained as a function of the scattering vector q = (4πsinθ)/λ, where 2θ is the scattering angle through azimuthal averages of the 2D data in the given angular range.
The chemical structures and molecular properties of the polyILs investigated in this work are shown in Fig. 1 and Table 1, respectively. Fig. 2 shows the differential scanning calorimetry (DSC) thermograph of the polyILs studied in the current work. PVB-MeTFSI shows the lowest Tg among all systems investigated and this result can be understood at a qualitative level. The glass transition temperature, Tg, of polyILs is thought to be controlled by cooperative dynamics determined by intermolecular interactions, which can be influenced by free volume or packing effects, coulombic interactions between anions and cations, and molecular weight.30 For the TFSI− anion, the charge delocalization between S–N–S bonds along with steric hindrances from sulfonyl and trifluoromethane groups, results in the reduction of the strength of coulombic interactions with the triazolium cation, leading to a lower Tg.22,34,35 However, the high charge density of the Cl− anion results in stronger ionic interactions with the triazolium cation, leading to the higher Tg of PVB-MeCl. Based on similar arguments, the higher Tg values of the polyanions compared to their polycation counterparts can be rationalized. In both polyanion and polycation systems, larger counterions result in lower Tg. This observation is consistent with previous studies that explored the counterion charge density in determining Tg of the polyILs.4 However, rationalization of the higher Tg in the polyanions requires a more in-depth study of the structure, as discussed later based on the wide-angle X-ray scattering (WAXS) results.
Fig. 2 Differential scanning calorimetry thermograph of the polymerized ionic liquids investigated in the current work. |
Fig. 3 shows the wide-angle X-ray scattering (WAXS) spectra presented in terms of the intensity versus scattering vector (q). The characteristic X-ray scattering peaks for the cationic polyILs indicate the rich mesoscale organization associated with polar, non-polar, and ionic moieties.30,36–39 Through a combination of molecular dynamics simulations and X-ray scattering experiments to analyze the mesoscopic structure in polyILs, three general structural features corresponding to the backbone-to-backbone (low q), ion–ion (intermediate q), and pendant group (high q) correlation distances have been identified.19,22,30,32,36,40,41
To quantitatively determine the peak position, asymmetric Lorentzian functions, are used to fit the X-ray scattering data.42,43 This equation represents an asymmetric distribution with the center at q0 = 2π/d, a full-width at half maximum (FWHM) of Δ, and an asymmetric parameter b. The fit parameters yield characteristic correlation lengths, presented in Table 2.
PVB-MeCl | PVB-MeTFSI | PVB-SO3C4MIm | PVB-SO3Na | |
---|---|---|---|---|
d backbone (Å) | 14.85 ± 0.1 | 12.69 ± 0.1 | — | 13.68 ± 0.1 |
d ion (Å) | 7.22 ± 0.2 | 8.16 ± 0.2 | 8.24 ± 0.2 | — |
d pendant (Å) | 4.65 ± 0.1 | 4.67 ± 0.1 | 4.51 ± 0.1 | 4.65 ± 0.1 |
Despite the similarity of WAXS spectra of polyanion and polycation samples, there are some unique differences. The low-q and intermediate-q peaks of PVB-SO3Na appear at lower q in comparison with polycations. The sharp low-q peak of PVB-SO3Na corresponds to the center-to-center spacing between ion aggregates (daggregate), which is associated with the characteristic length scale of 21.7 Å.44–47 The intermediate-q peak of correlation distance 13.7 Å may be related to the backbone-to-backbone distance, which is comparable to the polycation counterparts. On the other hand, the broad low-q peak of PVB-SO3C4MIm that appears in the q = 0.31 Å−1, daggregate = 20.3 Å, may be due to the merging of the backbone-to-backbone and ion aggregate peaks. The second intermediate-q peak of PVB-SO3C4MIm shows a correlation length scale of 8.24 Å, corresponding to the ion–ion correlation.
The length of the fully stretched pendant groups of polyanion and polycation polyILs were estimated to be 12.25 and 12.5 Å using the open-source Avogadro software. Thus, idealized, fully stretched, non-interdigitated systems should yield nearly double backbone-to-backbone correlation distances. Although this is a crude estimate, some insights may be gleaned from these values. The comparable pendant group sizes of the polycation and polyanion polyILs agree with dbackbone values from WAXS, implying interdigitation of the pendant groups. However, the association of SO3− and the free cations at the end of the pendant group, likely results in large aggregates while the butyl pendant group of the polycations presumably disrupts the ion aggregate formation.48–50 These ion aggregates can act as physical cross-linking points resulting in the segmental dynamics restriction and higher Tg in polyanions.
The dielectric spectra of PVB-MeTFSI at different temperatures are shown in Fig. 4. A combination of the empirical Havriliak–Negami (HN) and the random barrier model (RBM) functions were used to describe the dielectric spectra:51 where Δε is the dielectric relaxation strength, τHN is the relaxation time, α and γ are symmetric and asymmetric stretching parameters of the HN function in the first term. The second term represents the RBM part where σ0 is the dc ionic conductivity, τe is the characteristic charge transport time, ε0 is the vacuum permittivity, and ε∞ is the limiting high-frequency permittivity. The Aωn term is used to account for electrode polarization. The empirical HN model describes the orientational polarization effects associated with ion conduction.18,19 On the other hand, the RBM model accounts for charge transport. Within the RBM framework, ion hopping occurs by overcoming the spatially randomly distributed energy barriers.52 Above Tg, the counterion diffusion in the polyILs is coupled to some degree to the segmental dynamics or α-relaxation. However, the dielectric relaxations in this spectral range are not detectable due to the dominant contributions of the ionic conductivity to ε′′. To suppress the dc conductivity contributions, the derivative representation of the imaginary part of the complex dielectric function, defined as based on the Kramers–Kronig relation, is used. Fig. 5a shows variation in the dc ionic conductivity with the inverse temperature. The high dc conductivity of PVB-MeTFSI and PVBSO3C4MIm compared to their counterparts correlates with their lower Tg. The bulky structure of the counterions likely lowers the charge density and coulombic interaction strengths, resulting in faster segmental dynamics of the polymer backbone. To analyze the results at comparable timescales, the dc conductivity is plotted versus Tg-normalized temperature. Above Tg (Tg/T < 1), the conductivity of PVB-MeTFSI is slightly higher while PVBSO3C4MIm has the lowest conductivity.
This suggests that above Tg, ion transport is more effective in the PVB-MeTFSI at the same value of Tg/T. However, we could not study the σdc conductivity of PVB-SO3Na above Tg due to the possibility of sample degradation at higher temperatures. The temperature dependence of σ0 above Tg can be represented by the Vogel–Fulcher–Tammann (VFT) equation , where B, σ∞, and T0 are fitting parameters (Table 3), below Tg, the Arrhenius equation is employed. Fig. 5b shows the variation of the dc conductivity versus ωc(=1/τe) which follows the empirical Barton, Nakajima, and Namikawa (BNN) relation.51 The BNN relation correlates the dc conductivity to the dielectric relaxation strength through σ0 = pε0Δεωc, where p is a dimensionless empirical number, σ0 is the apparent dc conductivity, Δε is the dielectric relaxation strength, and ωc is from the RBM. Fig. 5b suggests that the underlying mechanisms of ion conduction and dielectric relaxation in polyanion and polycation polyILs are similar.21
Sample | σ ∞ | B | T 0 | E σ (kJ mol−1) | Decoupling index (m) |
---|---|---|---|---|---|
PVB-MeCl | 64 | 2927 | 207 | 165 ± 1 | 6.7 |
PVB-MeTFSI | 0.24 | 1013 | 227 | 123 ± 1 | 6.7 |
PVB-SO3-C4MIm | 18 | 2702 | 217 | 43 ± 1 | 5.9 |
PVB-SO3Na | — | — | — | 95 ± 1 | 6.1 |
The degree of decoupling, m, determined at Tg using m = 15 + log10(σ0), was proposed to quantify the coupling between segmental and ion conduction.53 The extent of decoupling of the ion conduction from segmental dynamics is dictated by the backbone and pendant group stiffness.54 The close decoupling indices suggest that the chain stiffness of polycations and polyanions are comparable. In the current study, the effect of the counterions of either the cation or anion on the degree of decoupling is negligible (see Table 3), indicating that m is likely controlled by the backbone structure. However, a close examination of the remarkable differences below Tg implies that m may not be an appropriate parameter for characterizing ion transport in polyILs.
The Eσ values can be calculated from the temperature dependence of ionic conductivity below the Tg using the Arrhenius equation , where σ∞ is a pre-exponential factor and R is the universal gas constant. More accurate estimates are obtained by the fitting data in a broad range of temperatures by combining the VFT and Arrhenius equations. Despite the comparable decoupling indices, polyanion polyILs show lower ion conduction activation energy, Eσ, than their polycation counterparts. Moreover, smaller counterions result in higher Eσ in polycation and polyanion polyILs. Deeper insights into the ion conduction properties of the polyILs can be gleaned by attempting to find correlations between their ion transport properties below Tg and their chemical structure.
In the polyanion systems, the strong coulombic interactions within the ion aggregates result in slower segmental dynamics, as reflected by the Tg values. Below Tg, polyanion systems show lower ion activation energy than their polycation counterparts. This result is presumably associated with the percolation of the ion aggregates. We conjecture that the larger characteristic lengths associated with the ion aggregates of the polyanions lead to frustrations in the packing of the polymeric chains and form percolated ion paths resulting in facilitated ion mobility below Tg. On the other hand, the interdigitation of butyl pendant groups leads to better-packed structures and steric hindrances that limit the available sites for ion hopping as reflected by the higher Eσ.
The counterion size is another factor that controls the Eσ. The overall effect of the counterion size and charge density on the Eσ can be interpreted using the Anderson–Stuart Model. Within the framework of this model, the activation energy is given by:55,56 where the first term describes the coulombic interactions between mobile ions with charge, q1 and radius R1, and counterions with charge q2 and radius R2. ε0 and ε are vacuum permittivity and medium static dielectric permittivity, and β is the Madelung constant. The second term defines the elastic contribution of the ion jump process. The G term is the shear modulus, RD is the radius of free volume for ion hopping, and λ is the ion jump length. In the polycation systems, it is observed that the Eσ of PVB-MeCl is higher than PVB-MeTFSI. The small size of Cl− results in stronger coulombic interactions with the triazolium ring. So, according to this model, the coulombic term dominates the Eσ. It was previously shown that the elastic modulus of polyILs increases with the ionic volume of the counteranion.57 However, the lower Eσ of PVB-MeTFSI correspond to the weaker coulombic interactions of the large TFSI− and backbone despite higher shear modulus.
The higher Eσ of PVB-SO3Na than PVB-SO3C4MIm can also be understood using the Anderson–Stuart Model. The strong coulombic interaction of the Na+ and polyanion backbone enhances Eσ. On the other hand, C4MIm+ can increase the shear modulus and decrease the coulombic interaction. The low Eσ indicates the effect of the weaker coulombic interaction dominates the higher shear modulus role in controlling ion conduction activation energy. This microscopic understanding of ion dynamics in polymerized ionic liquids will be valuable in designing solid polymer electrolytes. These results suggest that a subtle counterbalance between electrostatic and elastic forces determines the overall ionic conduction. A rational design of polyILs for practical applications should involve carefully selecting the target ionic groups that optimize the total contribution of these forces in a wide range of temperatures.
In summary, differential scanning calorimetry (DSC), broadband dielectric spectroscopy (BDS), and wide-angle X-ray scattering (WAXS) were employed to understand the correlation between chemical structure, nanoscale organization, and ion transport properties of 1,2,3-triazole based polymerized ionic liquids. Above Tg, polycation polyILs showed higher ion conduction compared to anionic backbone counterparts due to the coupling of the ion conduction and segmental dynamics. Below Tg, polyanions show lower Eσ of ion conduction than their polycation counterpart. The strong coulombic interactions between the backbone and the counterions result in higher Eσ while increasing the volume of the counterion leads to higher shear modulus and reduction of the Eσ, qualitatively consistent with the Anderson–Stuart model. The findings of this work provide insights into the rational design of solid electrolytes.
This journal is © The Royal Society of Chemistry 2024 |