Francesco Pedrolia,
Alessio Marranib,
Minh-Quyen Lea,
Olivier Sanseauc,
Pierre-Jean Cottineta and
Jean-Fabien Capsal*a
aUniv Lyon, INSA-Lyon, LGEF, EA682, F-69621, Villeurbanne, France. E-mail: jean-fabien.capsal@insa-lyon.fr
bSolvay Specialty Polymers, viale Lombardia 20, 20021 Bollate, Italy
cP2D, CNRS/Rhodia-Solvay, UMR 5268, 85 avenue des Frères Perret, F-69192 Saint Fons, France
First published on 26th April 2019
Electroactive polymers (EAPs) such as P(VDF-TrFE-CTFE) are very promising in the field of flexible sensors and actuators. Their advantages in smart electrical devices are due to their low cost, elastic properties, low density, and ability to be manufactured into various shapes and thicknesses. In earlier years, terpolymer P(VDF-TrFE-CTFE) attracted a lot of research due to its relaxor-ferroelectric property that exhibits high electrostriction phenomena. While widely used in flexible actuation, this class of material is still limited by the high electric fields required (≥30 V μm−1) to achieve sufficient strain levels (>2%). This inevitably leads to high levels of leakage current and thus a short lifetime. This paper proposes a new approach based on electro-annealing thermal treatment for a pure terpolymer P(VDF-TrFE-CTFE) matrix in order to limit the conduction mechanisms. This in turn reduces the dielectric losses at a high level of electric fields. The experimental results demonstrate that a huge decrease in leakage current of 80% is achieved for a wide range of electric fields (i.e. up to 90 V μm−1) with a 4-fold extension in time-to-breakdown at high voltage excitations of 40 V μm−1.
EAPs can be classified into two main categories, i.e. ionic and dielectric. For both electro-active polymer types, actuation is driven by an applied electric field, but the material deformation is steered by very different physical mechanisms. The working principle of the ionic EAPs is based on ionic exchange between an electrolyte and a polymer matrix upon the appliance of the electric field; in the case of dielectrics, polymer deformation is driven by the electrostatic force between the two electrodes generated by the external electric field.5 Here we focus only on this second category of polymers because they have wider versatility and a faster response to electrical stimuli—this makes them more promising in the development of active actuators. Unlike ionic polymers, dielectrics do not need an environment rich in ionic species to operate. Thus, they can be implemented in several kinds of atmospheres.6
Of the existing dielectric EAPs (silicones, acrylates, polyurethanes and PVDF-based polymers), the most performing one in terms of electromechanical conversion is fluorinated electrostrictive terpolymer P(VDF-TrFE-CTFE).7–10 The small value of Young's modulus of silicones and acrylates such as polydimethylsiloxane (PDMS) and polyurethane (PU) allows these polymers to reach very large strains.11 At the same time, it severely limits the electro-to-mechanical energy conversion.12,13 The semicrystalline terpolymer P(VDF-TrFE-CTFE) exhibits extremely large elastic energy density.14,15 Its peculiar semicrystalline morphology leads to strongly enhanced polarization levels thanks to the presence of cooperative nanopolar regions in crystalline domains.16,17
Despite the augmented electrostrictive performances of fluorinated dielectric polymers, high electrical fields are still required to attain satisfactory strain levels.1 The driving voltage is proportional to dielectric thickness, and electroactive polymer films remain in the range of tens of micrometers. They did not go beyond those levels of driving voltage that can represent a strict limitation to wide industrial and commercial applications.3,18 Thin polymer films with high quality and homogeneity are now achievable and scalable (even at very low thicknesses on the order of micrometers).19–23 However, the use of high electric fields with low-glass-transition polymers inevitably results in a high level of leakage current.24
The control of dielectric losses and their limits in terms of an exponential increase at operating voltage is fundamental to the development of high performing and long-life polymeric actuators.25 High levels of leakage current also increase energy dissipation26 and underlies self-heating.27,28 This promotes material degradation processes29–31 and leads to a dramatic drop in material dielectric strength or electrical breakdown. This severely reduced electrical breakdown does not allow EAPs to fully assess their potential for target applications such as energy storage or actuation.
The main objective of this work involves developing long-life polymeric actuators by perfectly controlling the dielectric mechanisms governing ionic conduction in pure P(VDF-TrFE-CTFE) terpolymers. To enhance its intrinsic electrical properties, we propose here a new electro-thermal annealing process for semicrystalline terpolymer. A 2-fold reduced ionic conductivity at low voltage input has been recorded at 0.1 Hz along with a 5-fold reduced leakage current at high voltage input (up to 90 V μm−1). Structural analysis, dielectric spectroscopy, and electrical characterization can identify the physical phenomena behind this decrease in conductivity. In addition, aging tests and electromechanical measurements can better assess the material performance—especially in terms of smart actuator design.
Films were prepared by casting a dissolution of 25 %wt. terpolymer powder dissolved in 2-butanone (also known as methyl ethyl ketone, or MEK). Before casting, the dissolution was filtered via a pressure column system equipped with polytetrafluoroethylene (PTFE) membrane filters with sub-micron pore dimensions. The polymeric dissolution was then cast on tempered glass plates and dried at room temperature. Prior to thermal treatment, the film was peeled away from the glass plate and laid down on PTFE foils to not induce any internal residual stress due to the different thermal expansion coefficients of terpolymer and tempered glass. To ensure complete solvent evaporation even in the case of thick film deposition (i.e. 90 μm), the prepared film underwent a first low-temperature thermal treatment 60 °C in a convection oven. We note that such a temperature did not lead to modifications in crystal morphology nor did it promote crystal growth. Finally, circular gold electrodes 20 mm in diameter and 17 nm-thick were coated on both side of the terpolymer film to obtaining a metal/polymer/metal capacitor-like architecture. Deposition used gold sputtering (Cressington 208 HR).
Fig. 1 Temperature and voltage profile for the standard thermal and the novel electro-thermal annealing treatments. |
Two samples were prepared based on the two different thermal processes: samples prepared with the standard annealing (“TH”) and samples prepared with novel electro-thermal annealing (“E-TH”).
DSC used a SETARAM DSC131 EVO calorimeter. The thermal characterization determined the potential effects on crystal morphology induced by the two different annealing treatments mentioned previously. The first series of post-annealed samples was heated from room temperature to 140 °C at 10 °C min−1. Subsequently, a second series was analyzed with a faster heating rate of 20 °C min−1 to ensure no modifications to the crystal morphology during the heating ramp. The thermograph for the second melting of the terpolymer sample showed that the optimal annealing temperature was identified as the temperature onset at melting peak.39 Finally, the crystallinity degree (or χc) of the post-annealed samples was calculated from integration of the melting peak divided by the enthalpy fusion of a hypothetic 100% crystal of terpolymer (i.e. 42 J g−1).40
The XRD was performed via an X'Pert Pro MPD Panalytical diffractometer using Cu-Kα radiation (λ = 1.5406 Å) of 45 kV and an electrical input of 40 mA in tandem with an incident-beam monochromator (Inc. beam Johansson 1xGe111 Cu/Co) and a X'Celerator detector. The diffraction patterns were recorded over an angular range of 10–30° (2θ) where the characteristic peak relative to P(VDF-TrFE-CTFE) crystalline phase was localized.40 A step length of the angular (2θ) equaled 0.017° with a counting time of 120 s per step. The extraction of the peak positions for indexing was performed via the X'Pert High Score.
The IR spectroscopy used “attenuated total reflectance” (ATR) because the spectra recorded via transmission of IR rays across the sample thickness were saturated. Measurement were achieved via a Spectro IR Alpha analyzer (Bruker) equipped with a ATR Diamant tip from 4000–400 cm−1 with a resolution of 4 cm−1 and 32 scans.
We emphasize that this setup can monitor the resulting current and the actuation displacement at the same time. The current amplification used a Stanford SR-570 amplifier connected to the moving electrode; the displacement was achieved via a high-precision capacitive sensor (FOGALE MC 940). The measured displacement referred to the thickness variation ΔdE driven by a input voltage excitation (longitudinal displacement) because it parallels the electric field direction (i.e. the 33-direction).42 Accordingly, the longitudinal strain S33 was given by:
(1) |
Finally, all signals including displacement, voltage, and current were simultaneously recorded in real-time using DEWE software (Sirius 8XSGT). Post-data treatment was performed thanks to Origin software.
To better understand the mechanical property of the terpolymer P(VDF-TrFE-CTFE), it is necessary to empirically determine its Young's modulus. The measurement was performed by recording the longitudinal tensile force under a given uniaxial displacement at 100 mHz. The 50 × 10 mm rectangular specimens with 90 μm-thick terpolymer films were used and a full description of the developed test bench was detailed in our previous works.43 Fig. 3 shows that determining the slope of the stress-versus-strain curve, an estimation of the Young's modulus values can be obtained equal to 148 MPa for both TH sample and E-TH sample. Consequently, the electro-annealing method does not modify the mechanical property of polymer. This is consistent with values previously reported in the literature.44,45 A mechanical model was implemented by fitting it to the experimental curves as shown in Fig. 3 demonstrating good coherence.
Fig. 3 Stress versus strain: experimental curves for P(VDF-TrFE-CTFE) films prepared via both standard annealing (black) and electro-annealing (red), and fitting curves (blue). |
The DSC results (Fig. 4) confirmed no modifications in crystalline morphology from the two different annealing methods. The thermographs had no variation in melting peak position, width, or appearance of bimodal shapes. This result demonstrated that the electro-annealing did not introduce any change in crystal formation to the crystal size or crystal size distribution. Fig. 4a shows the thermograph for the two differently annealed samples with an heating rate of 10 °C min−1. Integration of melting peaks results in values of melting enthalpy of 12.2 J g−1 corresponding to 29% degree of crystallinity.40 To verify that the heating rate does not distort the measurements—especially in case of fast crystallization dynamic of material—thermographs with faster heating ramp (20 °C min−1) were taken as illustrated in Fig. 4b. There were no differences between samples. The results of thermal analysis are summarized in Table 1.
Fig. 4 DSC thermographs for the two samples with two different heating rates: (a) 10 °C min−1 and (b) 20 °C min−1. |
TH sample | E-TH sample | |||
---|---|---|---|---|
10 °C min−1 | 20 °C min−1 | 10 °C min−1 | 20 °C min−1 | |
Tmelting (°C) | 114.6 | 114.3 | 114.3 | 114.0 |
ΔHmelting (J g−1) | 12.2 | 13.4 | 12.1 | 14.1 |
FWHM (°C) | 12.7 | 13.0 | 11.5 | 14.5 |
χc% | 29.3 | 31.9 | 28.8 | 33.6 |
A deeper investigation of crystal structure was carried out based the XRD analysis. The diffraction spectra (Fig. 5) shows diffraction peaks perfectly superimposable for the two samples because there was no variation in peak position (2θ = 18.4°) or peak width (FWHM = 0.83°). Therefore, no modification in crystal lattice space was found confirming that the electro-annealing treatment did not alter the crystalline phase conformation. Similar to the DSC thermographs, the XRD diffraction peaks do not exhibit any variation in crystal size distribution for both samples.
Fig. 6a depicts the ATR-IR spectra of the TH and E-TH samples based on two different annealing procedures. Both spectra were very similar affirming that the novel electro-annealing technique did not induce any structural modification to the P(VDF-TrFE-CTFE) films. These results are of particular interest and suggest that the improvements in electrical behavior of EAPs in terms of leakage current reduction and dielectric losses are mainly due to the interfacial polarization of ionic impurities. Fig. 6b shows IR analysis of the E-TH samples before and after 22 h of aging. Again, no structural modifications were observed.
Fig. 6 (a) IR spectra for TH and E-TH samples. (b) IR spectra for E-TH sample before and after 22 h of aging. |
Fig. 7 Current versus electric field for unipolar sinusoidal voltages. Black and red lines refer to experimental measurement of the TH-sample and E-TH sample, respectively. Blue curves represent the conduction current component modelled by eqn (3). |
The blue slope of Fig. 7 described the conduction (or leakage) current that can be estimated based on the following Hopping model:46,47
(2) |
For the sake of simplicity, the Hopping model can be rewritten as follows:
(3) |
(4) |
The evaluation of A0 and A1 parameters of each sample has been performed by fitting the curve of experimental current versus electric field. Table 2 summarizes the modeling results. As expected, only A0 parameters varied whereas A1 was constant for the both TH and E-TH terpolymers. This is consistent to previous works.33,48–51 Interestingly, A0 showed a 5-fold reduction for the E-TH sample with respect to the TH one, and this parameter is inversely proportional to the activation energy of conduction mechanism.
A0 (A) | A1 (A m2 V−1) | a (nm) | ATH0/AE-TH0 | |
---|---|---|---|---|
TH sample | 7.4 × 10−3 | 2.89 × 107 | 0.89 | |
E-TH sample | 1.49 × 10−3 | 2.89 × 107 | 0.89 | 5 |
Accordingly, the fitting results revealed that the proposed electro-thermal annealing did not cause any morphological modification of polymer amorphous phase;52 the remarkable 80% decrease in leakage current was inferred to a 5-fold increased energy barrier of the electronic conduction mechanisms.
The spectrum in Fig. 8a shows that under low frequency of 100 mHz, the dielectric losses tan(δ) of all E-TH samples dropped dramatically achieving a value of 0.078 instead of 0.15 as in the case of the conventional TH sample. This property leads to substantially decreased dielectric losses of almost 50%. This was caused by limited ionic conductivity of the E-TH samples. Similar results under very low frequencies were obtained for the E-TH discharged sample. Nonetheless, at higher frequency ranges (i.e., around 106 Hz where the relaxation peak occurs relating to the polymer chain α-relaxations (so-called dielectric glass transition)) there were no difference in dielectric losses behavior of the all samples. No structural modifications were observed regardless of which thermal treatment was selected: standard or electro-annealing process. These characteristics lead to an unchanged dielectric constant in the terpolymer films across the entire frequency range as demonstrated in Fig. 8b.
To conclude, high-voltage electronic conduction can be controlled via a simple polarization of ionic species presented in the polymer matrix. The polymerization agents used in the polymer synthesis of P(VDF-TrFE-CTFE) terpolymer contain ionic impurities that, although only on the order of ppm, contribute to interfacial polarization.24,53–55 These ionic impurities represent heterocharges that can be driven toward the oppositely charged electrode and accumulate under a constant DC voltage excitation.56 A schematic representation of this ionic polarization process is shown in Fig. 9.
In our case, during the electro-thermal annealing, the polarization process of ionic species was further promoted by the high temperatures at which the ion mobility in the polymer matrix is enhanced24 especially at longer times. Subsequently, the ionic charges turned out to be “locked” at the polymer/electrode interface during the cool-down step to room temperature when the molecular mobility was strongly reduced. This formed an electric double layer that had a dual effect on dielectric loss. First, their contribution to ionic conductivity was strongly reduced because the ionic species were constrained at the interface.53 This was observable in dielectric losses spectrum at low-frequency (Fig. 8a). Secondly, the accumulation of hetero-charges build up a local electric field (renamed Eions) opposed to the applied electric field that lead to a smaller internal electric fields and thus reduced driving force for electronic conduction57,58 as measured at high-voltage. This was confirmed by the modeled increased activation energy for electronic conduction in Subsection 4.2.
(5) |
Fig. 10 Current vs. AC electric field for (a) TH sample and (b) E-TH sample. Blue lines represent the conduction of Ohm's law, and slope values are dedicated to the polymer volume resistance. |
The ρ of the TH samples was estimated to be 6.3 × 109 ohm m. The E-TH sample (Fig. 10b) exhibited an asymmetric current that depends on the sign of the applied electric field. Based on the fitting results of eqn (5), two different values of the polymer volume resistivity were obtained, i.e. ρ = 6.5 × 109 ohm m for E < 0, and ρ = 12.0 × 109 ohm m for E > 0. Interestingly, the resistivity value of the E-TH sample was similar to that of the TH sample under a negative input voltage—this was not the case for the positive electric field.
According to the observation drawn in the previous Subsection 4.3, such an asymmetrically resistive behavior to electronic conduction of the E-TH sample is due to the electric double layer formed at the polymer/metal interfaces by ions polarization achieved during the electro-thermal annealing. Under appliance of positive electric field, as represented in Fig. 9, the ionic charges accumulated at the polymer/electrode interfaces during the electro-thermal annealing are oppositely charged with respect to the adjacent electrode. The so-built electric double layers can represent a limitation for injection of homocharges56 (i.e. electrons and holes) resulting in the reduced leakage current through the sample,60,61 which was not the case for the standard annealed TH sample. Moreover, overall internal electric field turns out to be reduced by the interfacial polarization of the ionic hetero-charges. Anions and cations separation builds up a local electric field that, as depicted in Fig. 9, opposes to the positive applied electric field leading to lower driving forces for electron motion,57,58 reducing the electronic conduction through the sample.
Consequently, using positive voltage excitation nearly doubled the resistivity of the E-TH sample resulting in considerably improved electrical properties like reduced leakage current as well as dielectric losses. In case of negative voltage appliance, the electrode polarity was reversed; thus, the electric double layer no longer holds resulting in unchanged resistivity values that is similar to the one of TH sample.
Fig. 10 highlights that the experimental current corresponding to a total intensity (itot) consists of a capacitive current (icap) and a conduction current (ileakage); the latter is given by eqn (5). Subsequently, the capacitive current can be easily estimated (eqn (6)) by subtracting the leakage current from the total measured current.
icap = itot − ileakage | (6) |
Eventually, the electric displacement D(E) is determined by time integration of the capacitive current density that is written by:
(7) |
Fig. 11 displays the electric displacement versus the applied electric field for the TH and E-TH samples. As expected, the polarization loop D(E) of the E-TH sample was asymmetric—this is the opposite of the case of the TH sample. This result is consistent with the previous analysis relating to the leakage current behavior.33
Fig. 11 Experimental electric displacement vs. electric field curves for (a) TH sample and (b) E-TH sample. Blue lines represent the polarization calculated from the Debye/Langevin model. |
The electric displacement loop of the two samples were then fitted by the theoretical Debye/Langevin model. Further details of this approach have been fully described in our previous works.33,62,63 The two-component Debye/Langevin formalism considers the capacitive current icap to constitute the sum of contributions generated by the dielectric response of two inherent morphological phases (i.e. crystal and amorphous) present in the semi-crystalline P(VDF-TrFE-CTFE):
icap = icrystalcap + iamorphcap | (8) |
(9) |
The relative dielectric permittivity εrθ(E) can be estimated based on the following formula:
(10) |
The experimental polarization curves as well as the corresponding fitting models are depicted in Fig. 11. The fitting results are in Table 3:
TH sample | E-TH sample | |||
---|---|---|---|---|
εE=0r | Esat (V μm−1) | εE=0r | Esat (V μm−1) | |
Crystal | 59 | 13 | 59 | 20 |
Amorphous | 10 | 400 | 10 | 400 |
The Debye/Langevin fitting result of the TH and E-TH samples shows that their intrinsic dielectric permittivity is identical in both the crystal and amorphous phases in contrast to the saturation electric field that is higher for the E-TH sample. For a better comparison, experimental loops and modeling curves for the two samples are superimposed in Fig. 12; characteristic polarization parameters are listed in Table 4.
Fig. 12 (a) Superposition of experimental electric displacement curves and (b) modeled polarization for the two samples. Direction of applied electric field is indicated by blue arrows. |
TH sample | E-TH sample | |||
---|---|---|---|---|
E < 0 | E > 0 | E < 0 | E > 0 | |
Pmax (C m−2) | 1.40 × 10−2 | −1.40 × 10−2 | 1.65 × 10−2 | −1.64 × 10−2 |
Premnant (C m−2) | 3.02 × 10−3 | −3.03 × 10−3 | 2.10 × 10−3 | −1.00 × 10−3 |
The relevant asymmetry in the electric displacement loop of the E-TH sample is clearly observed; the absolute value of its remnant polarization under a positive electric field is two-fold lower than under a negative input excitation. This effect is largely caused by the non-relaxed ionic polarization (Eions in Fig. 9) that counterbalances the intrinsic terpolymer remnant polarization induced by ferroelectric domains.
Fig. 13 Current density vs. aging time for one pair of specimens including the TH sample and the E-TH sample. Dashed blue lines represent the level of leakage current density. |
Fig. 13 shows that the current versus time measurement. This measurement was characterized by an initial transitory regime that represents the capacitive and ionic drift current components due to dipole polarization and ionic polarization, respectively. Next, a steady regime occurred that was merely related to the DC conduction effect. At longer testing times, the current intensity again became unstable and began to increase. This final regime related to a run-away process due to the self-heating of polymer created by electric conduction in a thermally-activated mechanism.29 This mechanism corresponds to the material degradation phase.
It has been observed that the E-TH film leads to a drastically decreased DC leakage current: the value dropped to 4.5 mA m−2 as opposed to 6.5 mA m−2 for the standard material. This result demonstrates a significant improvement in resistive properties of around 30% for the novel developed polymer. Moreover, the last regime of the aging experiment showed that the E-TH sample had current-induced self-heating effects after 10 hours of electrical solicitation. This effect occurred after only 2 h of testing for the standard TH sample. Similar behavior was observed where the E-TH polymer led to 4-fold longer time-to-breakdown as opposed to the standard sample. In conclusion, the E-TH polymer is clearly much more resistant than the conventional one—particularly under a high voltage excitation.
Fig. 14 Strain S33 vs. time under 20V μm−1 sinusoidal AC electrical input for (a) TH sample and (b) E-TH sample. Blue lines correspond to the applied electric field. |
To better assess the actuation performance of the EAPs, we provide here mathematical expressions of the electrostrictive coefficient (M33) and the electromechanical coupling constant (Q33).
The Q33 is defined by the relationship between the longitudinal strain S33 and the sample macroscopic polarization D(E) based on the following formula:42,62,64
S33 = Q33D(E)2 | (11) |
Hence, Q33 can be expressed as a function of the Young's modulus Y and the dielectric permittivity ε:
(12) |
(13) |
(14) |
Regarding the electrostrictive coefficient (M33), the quadratic relationship between the longitudinal strain S33 and the applied electric field E allows to identify its value:
(15) |
Importantly, the Young's modulus used in eqn (12)–(15) formally refers to the polymer's mechanical behavior under compression strains. In this study, it can be assumed that the Young's modulus is unchanged in both compression and tensile configurations considering the small deformation in the linear elasticity area and the low crystallinity of terpolymer.
Fig. 15a depicts the experimental electromechanical characterization as well as the corresponding theoretical model of the conventional TH materials. Longitudinal strain S33 versus applied electric field curve shown in Fig. 15a represents the typical electro-mechanical behavior of a ferro-relaxor material65,66 where the residual or remnant polarization, due to the not completely relaxed ferroelectricity, induces hysteresis in the longitudinal strain S33 versus applied electric field curves.67,68 Both theoretical models (13) and (14) perfectly fit the experimental curves. Eqn (12) with a Young's modulus of 148 MPa (Subsection 3.4) yielded an estimate of 13.2–14.5 m4 C−2 for the electromechanical coupling constant Q33 under an applied electric field varying from 0 to 20 V μm−1.
The electrostrictive coefficient M33, as defined by eqn (15), turns out to be not constant and dependent on the applied electric field. More precisely, the value of permittivity tends to somewhat decrease with the input electric field due to dipole saturation. However, this variation in permittivity is relatively low and it can be negligible under electric field not excess 20 V μm−1. Indeed, as observed in Fig. 12b, the dielectric response of both TH and E-TH samples is quasi-linear at such an input voltage range. This leads to unchanged electrostrictive coefficient M33 for an electric field range ≤ 20V μm−1.
Accordingly, from eqn (15), the estimate values of the quasi-constant electrostrictive coefficient M33 were 3.6 × 10−18 m2 V−2 and 4.4 × 10−18 m2 V−2 for the TH and E-TH polymers, respectively. The estimates of both Q33 and M33 were consistent with literature reports for polar electrostrictive polymers.13,17,18,62
As demonstrated in Subsection 4.4, the Debye/Langevin model was not valid for the E-TH sample because of its asymmetrical electrical polarization that was inverse to the standard TH material. Indeed, such asymmetrical dipolar properties cannot be correlated by a unique model parameter that describes the strain response under both negative and positive voltages. Thereby, the experimental strain of the E-TH polymer were fitted by eqn (11) with different parameters depending on the sign of the input excitation. Fig. 15b shows the empirical curves of the E-TH sample under negative electric fields. These curves were well-fitted based on the same values of M33 and Q33 used for the TH sample. Nonetheless, under positive electrical inputs, the E-TH polymer leads to increased electrostrictive coefficient M33 of 22% and enhanced electromechanical coupling constant Q33 of 14% with respect to the former TH terpolymer. For a better comparison, Table 5 summarizes the values of both coefficients characterizing the electromechanical performances of the E-TH and TH materials.
TH sample | E-TH sample | |
---|---|---|
M33 (m2 V−2) | 3.6 × 10−18 | 4.4 × 10−18 |
Q33 (m4 C−2) | 13.2–14.5 | 16.5 |
This journal is © The Royal Society of Chemistry 2019 |