Ultra-high Seebeck coefficient of a thermal sensor through entropic optimisation of ligand length of Fe(ii) spin-crossover (SCO) materials

In this work, we present a spin-crossover (SCO) complex molecular formulation [Fe(Ln)2](BF4)2 in an electrochemical single couple solution. A Seebeck voltage arises when an electrochemical single couple solution is subjected to a temperature difference, resulting in a single couple reaction at either terminal of the electrochemical cell. The ultrahigh Seebeck coefficients were obtained due to a number of molecular optimisation strategies. The [Fe(L16)2](BF4)2 complex demonstrated a maximum Seebeck coefficient of 8.67 mV K−1, achieved through a six-pronged approach to maximise entropy during the transition from low spin (LS) to high spin (HS) through: (i) a change in spin state, (ii) a change in physical liquid crystalline state, (iii) the spin Seebeck effect, (iv) the kosmotropic and chaotropic effect, (v) the fastener effect and (vi) thermal heat absorbance. A reduction of the Seebeck coefficient to 1.68 mV K−1 during the HS–LS transition at higher temperatures is related to the single spin state transition entropy change. In summary, this paper presents a systematic study to identify the contributing factors in the production of a sensor with an ultrahigh Seebeck coefficient for energy harvesting through the optimisation of its molecular entropy elements.


Introduction
The thermoelectrochemical effect offers an effective way of converting heat directly into electrical potential, and thus can become a generator of electricity. 1 However, there is still a challenge in the thermoelectrochemical generator (TEG) system that should be overcome from the low current, as well as the electrical potential due to the material properties, which ultimately affect their performances. 1 Therefore, this work explores an alternative application of the thermoelectrochemical effect as thermal sensors by using spincrossover (SCO) complexes. The thermal sensor studied in this work has the same operating principle as the TEG, but the ultrahigh Seebeck (S e ) coefficient is achieved using the SCO complex by primarily exploiting its intrinsic phase and magnetic transitions, namely the liquid crystal (LC) effect and spin Seebeck effect (SSE). [2][3][4][5] The sensitivity of the SCO complex to a thermal stimulus has driven entropy changes, which subsequently drives a Seebeck voltage. The entropy change (DS) and the transport properties of the charge carrier affect the magnitude and sign of S e . The performance of the thermal sensor based on the Seebeck effect can be dened according to the gure of merit (ZT) below: Here, T is the absolute temperature, S e is the Seebeck coef-cient, s is the ionic conductivity and k is the thermal conductivity. 6 This work introduces three new novel spin-crossovers (SCO) based Fe(II) complexes which are [Fe(L 12 ) 2 ](BF 4 ) 2 (complex 1), [Fe(L 14 ) 2 ](BF 4 ) 2 (complex 2) and [Fe(L 16 ) 2 ](BF 4 ) 2 (complex 3) where L 12,14,16 were N 3 Schiff bases appended with linear C 12,14,16 carbon chains at the N atoms. These Fe(II) SCO complexes are characterized for their Seebeck voltage as a function of temperature in order to identify the effect of ligand length on the Seebeck voltage. Previously, Ibrahim et al. 7,8 have reported the use of similar Fe(II) metal centre with different structural complexes, and provided the structural correlations between molecular and ionic structures with spin state behaviour and on thermoelectric performance. The results exhibited the Seebeck coefficient with À0.51 mV K À1 for [Fe(L) 3 ](BF 4 ) 2 $2H 2 O (Fe 2+ -mono), À0.57 mV K À1 for [Fe 2 (CH 3 COO) 4 (L) 2 ]$2H 2 O (Fedinuc) and À0.52 mV K À1 for [Fe 2 (CH 3 COO) 4 (L) 2 . 7,8 SCO complexes have indicated great potential in thermal sensor applications, owing to their high S e values between À0.088 mV K À1 and 1.8 mV K À1 , reported by other researchers, listed in Table 1.
The switchable behaviour of the SCO material between highspin (HS) to low-spin (LS) or vice versa as shown in Fig. 1 distinguishes this material for its capacity to convert electrochemical energy to electrical potential, in response to temperature gradients. 12 Fe(II) complexes are the most extensively studied SCO materials, 13,14 as the change of spin state transition during electron transferring is capable of producing high entropy. Theoretically, Fe(II) complexes have good thermoelectric behaviour due to its large DS value. The DS spin for Fe(II) ¼ R [ln(2S + 1) HS À ln(2S + 1) LS ] is 13.4 J K À1 mol À1 , is higher than Fe(III) of 9.1 J K À1 mol À1 and Co(II) of 5.8 J K À1 mol À1 . 15 This is due to the electronic transition of Fe(II) complexes occurring between LS with no unpaired electrons (S ¼ 0), and HS with four unpaired electrons (S ¼ 2), as referred to the diagram in Fig. 1. 16,17 In terms of metal-ligand (M-L) bonding (between Fe and ligand), the LS state has a shorter bond which induces the N-donor on the ligand to easily donate electrons to the Fe centre. The positive charge at the Fe centre becomes unbalanced as the number of electrons are donated to the metal centre. The coulombic attraction between cation and anion becomes weak and easy for bonds to break, thus increasing entropy. 7 In contrast to the HS state which has a longer bond between Fe-L, the N-donor ligand has difficulty in contributing electrons to the Fe centre. The number of charges positive (at the metal centre) and negative (at N-donor ligand) cause the interaction between cation and anion to become stronger, thus requiring high energy to break the coulombic attraction, and hence lower the entropy. This is the reason why the LS state of Fe(II) complexes can produce higher entropies compared to the HS state, and in turn result in a higher Seebeck coefficient. Entropy is one of the important parameters in determining the Seebeck coefficients of thermal sensors. [18][19][20][21] Seebeck coefficients, S e can be described by the following equation: 22 In eqn (2), parameter n refers to the number of electrons involved in the redox reaction, F is the Faraday constant and DV is the potential generated by the temperature difference, DT. If the S e is in a positive value, the entropies of the cations are increased. 11 The advances in SCO material development can be seen through the combination of this material with its liquid crystal properties, otherwise known as metallomesogens (thermotropic liquid crystals). The liquid crystal phase transition results from the transition process between the crystal phase and liquid phase which can be observed through polarizing optical microscopy (POM). Metallomesogens complexes provide advantages for processing SCO materials in the form of thin lms, enhanced spin transition signals, switching and sensing in different temperature regimes and achievement of photoand thermochromism. 13,23 This complex structure change associated with the crystalline (Cr) 4 liquid crystal induces the spin state transition, thus contributing to an increase in entropy and Seebeck values, as reported by Gaspar et al. 24 The spin Seebeck effect (SSE) refers to the generation of spin voltage as a result of a temperature gradient in magnetic materials. SSE have been introduced as one of the novelties that can be induced for high entropy and Seebeck coefficient through charge transport properties in SCO materials. Megat Hasnan et al. has reported that SSE of polymer-based electrolyte in a thermoelectrochemical cell study has induced the spin S e and can act as a thermal booster to enhance ionic mobility and diffusivity. 5 Besides that, the spin state transition can cause two other effects such as the kosmotropic and chaotropic effect. The kosmotropic effect (inducing structure) is the ordered structure of complexes by stabilizing intramolecular interactions between Fe-ligand, whilst the chaotropic effect (inducing disorder) produces the opposite effect by disrupting and destabilizing the complex structure. 25 Hence, the kosmotropic and chaotropic effects involve the structure of complexes from an ordered structure to a disordered structure, brought to increase the entropy and consequently enhance the Seebeck effect. 26 The alkyl ligand lengths that are attached to the metal centre can provide variation in entropy and Seebeck measurements. A long alkyl ligand has a fastener effect where it forms a exible and so space eld around the central metal complex moiety through the thermal motions. 27 The exible alkyl chains cause the emergence of a strong cooperative effect called interchain interaction, formed between the adjacent chains of the surrounding ligands. This exibility can contribute to enhance the entropy and Seebeck values by facilitating the formation of spin state transition of the centre metal ions. 28,29 The thermal heat absorbance can contribute to maximise the entropy. This thermal heat can be measured using differential scanning calorimetry (DSC) where it records the phase change of SCO complexes during the heating and cooling process by observing the endothermic and exothermic peaks. The SCO complexes will absorb the heat to break the bond since ligands will begin to dissociate at certain temperatures, hence promoting the heat capacity and increasing the entropy. Throughout this work, we have focused on a six-pronged strategy which contributes to maximizing entropy during the transition from low spin (LS) to high spin (HS): (i) change in spin state, (ii) change in physical liquid crystalline state, (iii) spin Seebeck effect (iv) kosmotropic and chaotropic effect, (v) fastener effect, (vi) thermal heat absorbance and its impact on the Seebeck coefficient, which is especially apparent in the heating range between 298 K and 373 K.

Methodology
This section describes the synthesis of each complexes by the facile one pot reaction method, the preparation of Fe(II) complex electrolytes, and the characterization of each of the complexes on physical, electrochemical and thermal sensor properties.

Materials
Tetrabutylammonium tetrauoroborate (TBATFB) with a molecular weight of 327.29 g mol À1 , purchased from Merck, was used as a charge supplier in the Fe(II) complexes electrolyte. Dimethyl sulfoxide (DMSO), anhydrous 99.9% with a molecular weight of 78.13 g mol À1 , purchased from Sigma Aldrich, was used as a solvent to dissolve the solute. All the chemicals were of analytical grade and used without further purication.

Physical characterization of Fe(II) SCO complexes
The characterizations of Fe(II) SCO complexes with different ligand lengths are as listed below: (a) Fe(II) SCO complexes (i) Magnetic susceptibility was used to indicate the condition of the complex either in low spin or high spin at room temperature.
(ii) Mesomorphism was observed during the transition of liquid crystal properties of complexes by using polarizing optical microscopy.
(iii) A thermal study was used to evaluate the mass loss by thermogravimetric analysis.
(iv) Determination of enthalpy and predicted peak for endothermic and exothermic reactions were used to observe the spin state transition peak by differential scanning calorimetry.
(b) Fe(II) SCO electrolytes (v) Cyclic voltammetry was used to observe the single couple potential reaction and estimate the diffusion, mobility and electrochemical bandgap.
(vi) Electrochemical impedance spectroscopy was used to measure the ionic conductivity.
(vii) Thermal sensor properties were studied to measure the Seebeck value and power output.
2.4.1. Thermal studies. The thermogravimetric analysis (TGA) traces were recorded using a Perkin-Elmer Pyris Diamond  4 ] was used to calibrate the instrument.
2.4.3. Mesomorphism. Photomicrographs were viewed under Nikon-H600L Eclipse Microscopes, equipped with a Mettler Toledo FP90 central processor and a Linkam THMS 600 hot stage with a magnication of 50Â. The samples were kept in an oven at 333 K for two days prior to the analysis. This was to remove absorbed atmospheric moisture. The heating and cooling rates were 10 K min À1 and 3 K min À1 respectively in the temperature range between 298 K and 473 K. This is to observe the LC phase change.
2.4.4. Cyclic voltammetry. The cyclic voltammetry for Fe(II) SCO electrolytes under the temperature gradient were measured using Gamry Instruments. The measurement was performed between 2.0 V and À2.0 V with a 0.05 V s À1 scan rate. The current resulting from the oxidation and reduction was then measured. The oxidation peak (i p ) was calculated using the Randles-Sevcik's equation as follows: Here, n is the number of electrons exchanged per molecule, A is the electrode area, C is the compound concentration, v is the potential scan rate, F is the Faraday constant and R is the gas constant. The electrochemical bandgap, E e was then calculated using the relationship as shown below: In eqn (4), HOMO ¼ (onset oxidation peak voltage + 4.4) eV and LUMO ¼ (onset reduction peak voltage + 4.4) eV. The mobility (m) of the charge carrier was estimated using the Nernst-Einstein relationship.
The diffusion (D) value in eqn (5) was obtained from Randles-Sevcik's equation. Parameter k b refers to the Boltzmann's constant (1.38 Â 10 À23 m 2 kg s À2 K À1 ), where e is the electron charge (1.602 Â 10 À19 coulomb) and T is the absolute temperature.
2.4.5. Electrochemical impedance spectroscopy. In order to determine the impedance of the Fe(II) SCO electrolytes, the solution was placed into a Teon holder. The conductivity of the electrolyte was calculated using the equation below: Here, A is the electrode-electrolyte contact area, and d is the distance between two electrodes in the Teon holder. Parameter R b is the bulk resistance of the Fe(II) SCO electrolyte which was obtained from the intersection of the complex impedance plot with the real impedance axis.
2.4.6. Seebeck coefficient (S e ). The measurements were made in a non-isothermal cell containing approximately 1 ml of solution, as shown in Fig. 2. The cell consists of two platinum electrodes of which one side is a cold electrode while the other side is a hot electrode. The hot electrode was heated by a hotplate. The potential difference between the cold and hot electrodes was measured using a KEYSIGHT 34461A high impedance electrometer. Initially, both compartments were kept at room temperature, and the temperature of hot electrode was increased by 10 K. The potential difference between two electrodes was measured. The heating was continued until a temperature difference (DT) of 20 K was obtained. A graph of DT against the potential difference was plotted. From the plot, a gradient was obtained from the straight line which gives the S e values of the corresponding Fe(II) SCO electrolyte system.
2.4.7. Power output. The power output of the thermal sensor with Fe(II) SCO electrolytes at various temperature ranges were measured at a constant resistance of 1 kU. The voltage value was collected aer 30 minutes to allow the voltage and power output to stabilize. 30,31 The power output reached maximum, P max when the external resistance is equal to the overall internal resistance of the cell.
Fe(II) SCO compounds oen exhibit thermochromism (dark purple 4 white colour change) due to the rearrangement in the molecules of the complexes. The LS state absorbs visible light more strongly than the HS state and usually at distinctly different wavelengths, making UV-vis spectroscopy a useful tool for observing SCO. The HS state has more electrons in the antibonding e g orbitals compared to the LS state, thus has longer metal-ligand bond lengths. The electronic absorption spectral data in CHCl 3 tabulated in Table 3 shows d-d bands at approximately 700 nm assigned to the overlapping of 1 A 1g (F) / 1 T 1g and 1 A 1g / 1 T 2g electronic transitions respectively. Three other bands that were observed were metal-to-ligand charge transfer (MLCT) and intraligand charge transfers. It is interesting to note that the MLCT band occurred at low energies for all complexes.

Magnetic susceptibilities
The magnetic data (c M T) obtained by the Gouy method at 298 K is used to deduce the spin state of Fe(II) atom either at LS or HS. In the case of Fe(II), the HS state (S ¼ 2) is paramagnetic, while the LS state (S ¼ 0) is diamagnetic. Thus, the Fe(II) SCO switches the paramagnetism to "on" and "off". The room temperature magnetic susceptibilities (c) for complex 1 to 3 were initially measured at 298 K and are tabulated in Table 4.
From Table 4, complexes 1 and 3 gave almost identical c M T values. Complex 1 gave a c M T value of 0.31 cm 3 K mol À1 made up of 10.3% of HS and 89.7% LS. Complex 3 has a c M T value of 0.27 cm 3 K mol À1 which consists of 9.0% HS and 91.0% LS. The c M T value for complex 2 obtained 0 cm 3 K mol À1 with 100% LS. It is reported that a full HS Fe(II) SCO complex (S ¼ 2) typically exhibits at c M T z 3.00 cm 3 K mol À1 and a LS Fe(II) SCO complex (S ¼ 0) is visible at c M T z 0 cm 3 K mol À1 . [29][30][31][32][33] Based on this study, it may be inferred that the prepared complexes are mostly LS Fe(II) at this temperature. Hence, this is in line with the data obtained from electronic absorption spectra. The SQUID magnetometry elucidated a low spin state for Fe(II) complexes at room temperature, revealing a diamagnetic electronic structure. 34

Thermal properties
The decomposition temperatures and thermal properties of Fe(II) SCO complexes with different ligand length [Fe(L 12 ) 2 ](BF 4 ) 2 (1) and [Fe(L 16 ) 2 ](BF 4 ) 2 $H 2 O (3) were studied by thermogravimetric analysis (TGA) and differential scanning calorimetry (DSC) respectively. Table 5 tabulated the major mass losses and strongest peaks obtained for both samples. From the TGA, the major mass losses for each complex agreed with the values calculated from the loss of two BF 3 and two ligands. 2 Besides that, complex 1 suffered initial mass losses of 0.8% (temperature ranged from 357 K to 513 K) due to the evaporation of lattice H 2 O. Thus, the thermal stability of complex 1 (T dec ¼ 515 K) is higher than complex 3 (T dec ¼ 463 K) due to a higher HS of complex 1. This produced a strong coulombic attraction in the ionic complex, and as a result, longer metal-to-ligand bonds. 7 Based on Table 5, it can be seen that the measured mass loss (percentage) is in good agreement with the theoretical results (percentage) obtained from their structural formula, as it indicates the purity of the complexes and supports the proposed structural formula.
From the DSC, the strongest peak for all complexes in the temperature range of 315 K to 350 K was assigned to Cr-to-Col transitions. It is in agreement with the polarizing optical microscopy (POM). During the heating process, the endothermic peak corresponds to the rst order phase transition that accompanies the spin-crossover. The enthalpy of transition from LS to HS state for complex 3 (DH ¼ +79.2 kJ mol À1 ) is higher than the enthalpy during the endothermic peak of complex 1. This can be attributed to the ligand length effect which was aided to promote the reaction enthalpy formed for this complex. 35 Therefore, the reaction entropy produced for  a l is wavelength and 3 max is molar extinction coefficient. complex 3 (0.25 kJ mol À1 K À1 ) is higher compared to complex 1 (0.19 kJ mol À1 K À1 ). During the cooling process, complex 3 exhibited the highest exothermic peak of DH ¼ À32.2 kJ mol À1 and showed a phase transition from isotropic to solid. Complex 3 was expected to produce a better Seebeck performance as the thermal sensor material of this complex can generate a higher entropy reaction than complex 1. As a result, a spin-crossover behaviour occurred in the former complex. Moreover, complex 3 can undergo a reversible spin-state transition during the heating and cooling process. This is of great importance its application as the complex can absorb more heat easily due to the convection process that takes place in the thermal sensor cell on the cold electrode. This is to promote high reaction entropy which contributes to the Seebeck effect. It is noted that there were no corresponding peaks on cooling for complex 1, indicating the slow rate of phase formation.

Mesomorphic properties
The mesomorphic properties of complex 1 to 3 were studied by polarizing optical microscopy (POM). Under POM, complex 1 was observed to melt at about 343 K, in agreement with the DSC data, and cleared to an isotropic liquid (I iso ) at 402 K during the heating process. Upon cooling from the isotropic liquid phase, the dendritic-like optical textures [36][37][38] was developed at 398 K, which corresponds to a columnar mesophase (Col) as shown in Fig. 4(a). This texture gradually transformed to the mosaic character of a crystalline phase when further cooled to room temperature. Similarly, complex 2 melted at about 348 K and cleared to an isotropic liquid at 437 K. Aer further cooling, an optical texture like complex 1 (Col) was developed at 404 K. This texture then transformed to a broken-fan texture at 390 K. Subsequently, complex 3 melted at approximately 362 K, and cleared to an isotropic liquid at 472 K. Upon cooling, a focal conic fan-like (Col) optical texture [39][40][41] was formed at 397 K ( Fig. 4(b)). The results showed that the melting point and the phase transition from clear crystalline to isotropic liquid point became higher with an increase in the ligand length of the complexes. These behaviours were ascribed to the fastening effect of the long alkyl chains. The weaker van der Waals interactions of the long alkyloxy chains of ligand bonded to the LS Fe(II) is due to the less exible geometry at this site. It possesses a stronger Fe-L bond compared to those bonded to the HS Fe(II). As a result, L at the HS site melted at a lower temperature, which enabled the Fe-L bonds to become shorter (stronger), leading to a change in its electronic conguration to LS. 3 Hence, complex 3 revealed a high Seebeck generation owing to the fastening effect of the long alkyl chain that affects the melting point and phase transition point.

Electrochemical reversibility of the Fe(II) SCO complexes
The ligands that surround the Fe 2+ ion will create a ligand eld (spin state). This causes a redistribution of the electrons over the orbitals. This condition results in different spin states depending on the strength of the ligand eld, thus inuencing the single couple potential. In order to investigate the electrochemical properties of complex 1 to 3, a cyclic voltammetry (CV) was carried out. The measurements were conducted at various temperature ranges and repeated by three cycles for each temperature. The complexes showed good electrochemical reversibility. The charge supplier (TBATFB) was added into the solution as it easier to undergo oxidation and reduction processes due to the lower single couple potential. As electron transfer occurs at the electrodes, the charge supplier will migrate to balance the charge and complete the electrical circuit. 42 The CV scan for the complexes shows the anodic and cathodic peaks, as indicated respectively in ESI 1. † All complexes showed the anodic and cathodic potential peaks in the range outside À1 V and +1 V. These peaks can be assigned to the oxidation of L to L + and the reduction of L + to L. 43 The Fe ion peak is absent during the oxidation and reduction for all complexes. This is attributed to the spin state transition either in LS or HS, where the HS state have longer bonds between Fe-L. Hence, the Fe metal centre was shielded by the long alkyl ligands which acted as an insulating layer, hindering from further participation of the Fe in the single couple reaction at the electrode. 11,44 In contrast, complexes in LS have shorter bonds between Fe-L, allowing donation of electron to take part easily. Many negative charges surrounded the Fe metal centre and made it more difficult for the Fe ion to oxidize due to repulsion towards the BF À anion. The aromatic structure of complexes have shielded the solvent to oxidize at the metal centre, thus limiting the solvent reorganisation of the entropy. 45 Referring to Table 6 and ESI 1, † the increase of temperature gave the anodic peak shi towards positive voltage and the cathodic peak shi towards negative voltage. This behaviour is characteristic of a diffusion-controlled reaction.
Subsequently, ESI 2 † showed the comparison of the CV graph changes between three complexes at temperatures of 303 K, 333 K and 363 K. At temperatures of 303 K and 333 K, complex 1 and complex 3 indicated ligand peaks for oxidation and reduction, but no peak was observed for complex 2. At 363 K, all complexes are observed to have cleared oxidation and reduction peaks. This indicates that the complexes are easier to undergo oxidation and reduction at high temperatures. 46 The diffusion, electrochemical band gap and mobility are calculated by using eqn (3)-(5), as summarised in Table 7.
Referring to Table 7, the values of diffusion, mobility and electrochemical bandgap for all complexes increased with temperature. The diffusion coefficient is related to the movement of ions in the electrolyte caused by an external forced eld such as the temperature. This occurred when the ion is removed from the electrolyte by the discharge at an electrode. Besides that, the diffusion coefficient is related to the ionic mobility where the higher diffusion affects the ionic mobility of the charge carrier. This is based on the Einstein relation in eqn (5) that showed the ionic mobility is proportional to diffusion. The diffusion coefficient for complex 3 is higher than complex 1 and complex 2 at 333 K due to weaker Fe-L bonds of the HS. A weaker coulombic interaction enabled the ion to move easily, this resulted in high diffusion as well as high mobility. 27,47 Moreover, complex 1-3 experienced a spin Seebeck effect due to the increasing value of diffusion and ion mobility of the charge carrier that have arisen aer the spin state transition. In Table  7, it can be seen that the electrochemical bandgap for complex 3

Effect of ligand length Fe(II) SCO complexes on the ionic conductivity
Results in Fig. 5 and Table 8 presented the ionic conductivity at various temperatures for each Fe(II) SCO electrolyte. The ionic conductivity of complex electrolytes were calculated using the eqn (6). From Fig. 5, it is observed that the ionic conductivity increased linearly with temperature for complex 1, complex 2 and complex 3. The values of conductivity are tabulated in Table  8. At room temperature (298 K), complex 2 gave the highest conductivity value compared to complex 1. As discussed before, complex 2 is 100% LS. Complexes in the LS state are more conductive compared to the HS state. This is due to shorter bonds between Fe-L which has compressed the structure, hence easing the donation of the electrons. At temperatures of 313 K and above, it can be seen that complex 1 revealed the highest conductivity value compared to other complexes. This observation can be attributed to the shorter 12-carbon linear alkyl chains that assisted in improving the ionic conductivity, as compared to complex 2 and complex 3 that have relatively lower ionic conductivity due to long alkyl chains. The graph presented a linear ionic conductivity trend for all complexes due to the spin Seebeck effect that has developed the charge transport properties aer the spin state transition. The ionic conductivity is highly sensitive to the changes induced by the spin-crossover transition. This is through changes in ionic conductivity via chemical compression, driven by the structural changes. In particular, it is the phase change of the liquid crystal order in solution that was thought to have somewhat retained and contributed to one of the causes for the increase in entropy. 49 Any perturbations will result in higher entropy change. This event will then be elucidated from the view of ligand length of the SCO complex. Fig. 6(a)-(c) depicted the results of potential difference as a function of the temperature of hot electrode (T hot ). The result indicates two sign gradients for the Seebeck coefficient of each complex, as tabulated in Table 9. A higher boiling point of DMSO (462 K) solvent allows the complex  and cobalt-based redox couples [Co 2+/3+ (2,2 0 -bipyridyl) 3 (NTf 2 ) 2/3 ] (2.19 mV K À1 ). 9 This is due to the dependency of the Seebeck coefficient on entropy, and in such SCO systems, the signicant enhancement in entropy of SCO complexes comes from a sixpronged approach, as shown in Fig. 6(c). Complex 3 shows the spin state transition from LS-HS occurred at 343 K which caused the structural changes and the enhancement of the entropy. During the LS state, complex 3 has short metal-ligand bonds that easily allowed the N-donor to donate electrons to the metal centre. The increase of charge density at the metal centre causes the metal-ligand bonds to break down due to a weak coulombic attraction. The kosmotropic effect occurred during the spin state transition of LS-HS. This effect increased the order structure of complex 3 by providing the pathway for electron transfer. 51 The fastener effect is another factor which contributed to the enhancement of entropy through the interchain interaction formed by long alkyl chains of complex 3. The exibility of long alkyl chains of complex 3 formed a strong cooperative effect between the adjacent chains of surrounding ligands that enhanced the entropy and Seebeck value. The spin state transition of centre metal ions was affected by the exibility of long alkyl chains of complex 3. 28,29 The endothermic peak at 317 K with 79.2 kJ mol À1 belongs to the rst order phase transition that promoted the enthalpy reaction, thus signicantly increasing its entropy. All four factors contribute to maximising the entropy activity, resulting in the highest Seebeck value of 8.67 mV K À1 for complex 3.
The second gradient showed the spin state transition which occurred from HS-LS at 343 K to 373 K. The metal-ligand bond during the HS state is longer and is difficult for the N-donor to donate the electron to the metal centre. The charge density was maintained at the metal centre and anion ligand due to the stronger coulombic attraction, making the breaking of the bond more difficult. Between the temperature ranges, the phase transition from Crystalline to Columnar also occurred at 362 K, as observed by the POM. Moreover, the enhancement of the chaotropic effect contributed to the entropy due to the disordered structure of complex 3. 11 The spin Seebeck effect has developed aer the spin transition and contributing to give the Seebeck coefficient of 1.68 mV K À1 which is lower than a second gradient of complex 2.
Referring to Fig. 6(a), the Seebeck coefficient of complex 1 was the lowest (0.51 mV K À1 ) compared to other complexes as the spin state transition from the HS-LS and the kosmotropic effect at 293 K to 333 K contributed to the entropy. Referring to the second gradient, complex 1 produced the highest Seebeck coefficient of 2.10 mV K À1 . This can be attributed to the spin Seebeck effect of complex 1 from the LS-HS, in which the    that the Seebeck coefficient of complex 2 is 0.88 mV K À1 which is higher than complex 1. This event is due to the spin state transition from LS-HS, as well as the kosmotropic effects which occurred from 298 K to 333 K. The second gradient for complex 2 produced the Seebeck coefficient of 1.93 mV K À1 which is higher than the second gradient of complex 3. This circumstance is related to the spin Seebeck effect from HS-LS, the chaotropic effect and the phase transition (Crystalline / Columnar) at 348 K. All six effects mentioned above have contributed to maximising the entropy, hence increasing the Seebeck value. were tested in a thermal sensor that were connected with two platinum electrodes and a resistor in parallel to the circuit. The hot sided electrode was heated from 298 K to 403 K. The highest power output was produced when the applied external resistance is equal to the internal resistance of the cell. 52 Referring to Fig. 7, it can be seen that complex 2 revealed the highest power density with the value (P max ¼ 0.166 mW m À2 ) at the temperature of 373 K compared to other SCO electrolytes. This is attributed to a combination of a higher Seebeck coefficient at the second gradient and spin state transition from HS-LS which occurred at 373 K. Complex 2 in a LS state is higher in charge density, indicating a higher capability to transport charge to the respective electrodes and accumulate on the electrode surface. However, higher ion concentrations have the effect of reducing the ion diffusivity (3.48 Â 10 À7 cm 2 s À1 ) and mobility (11.10 Â 10 À6 cm 2 V s À1 ). In addition, additional factors also come from the changes in physical liquid crystalline (Crystalline / Columnar state) and chaotropic effect. The spin Seebeck effect that has occurred aer the spin state transition contributed to the high charge transport properties. 5 During these occurrences, the potential in the thermal sensor increased dramatically. As a result, high cell power density is obtained. Although the highest power density obtained in this present study is 0.166 mW m À2 for thermal sensor with complex 2, this value is still considered low if compared to the reported value by Abraham et al. 9 of 499 mW m À2 with a cell consisting of 0.1 M Co II/III (bpy) 3 (NTf 2 ) 2/3 in MPN electrolyte. This can be due to the fact that power density is directly proportional to concentration. The low concentration of complex 2 (0.002 M) has restricted the performance of the thermal sensor device due to the diffusion limitations. 53 Another important factor that limits the device performance is through the observation of the relationship between the current and potential of electrochemical device.
In Fig. 8(a) and (b), the current and power density obtained are low for all complexes systems which led to the limitation of diffusion, thus revealing a smaller power output. This clearly showed that the power output and current-voltage relation have inuenced the overall performance of the device. The following are strategies to improve the performance of thermal sensor device; (i) optimization of the electrolyte concentration must be conducted, 9,53 (ii) the use of the redox couple instead of the single couple can increase the redox reaction, 11 (iii) the electrode surface area should be increased by replacing them with carbon nanotubes which increases the number of reaction sites, hence producing a higher power output. Baughman et al. 54 proved that high surface area carbon nanotube electrodes perform better than platinum in aqueous thermoelectrochemical systems. This strategy has been proved by Abraham et al. 9 where the thermocell power density increased from 499 to 522 mW m À2 for the MPN-based device and from 183 to 240 mW m À2 for [C 2 mim][B(CN) 4 ]. Furthermore, another strategy is to focus on the device level by introducing a thermal separator between hot and cold electrodes can increase the thermal gradient within the electrolyte. The optimization of the thermal separator thickness can contribute to maximizing the performance of the device. 55  with the longest carbon chain length of ligand yielded an ultrahigh Seebeck value of 8.67 mV K À1 due to a six-pronged strategy that increased the entropy of the molecular SCO system. This contributed to changes during spin state, physical phase (liquid crystal mesophase), spin Seebeck effect, kosmotropic and chaotropic effect, fastener effect and thermal heat absorbance. Additionally, in a real thermal sensor (where one terminal is hot whilst the other is cold), an enhancement of the thermal property of the SCO solution due to a strong endothermic peak (as illustrated by the DSC analysis) further enhances the entropy of the reaction. Ultimately, all these properties point towards a structured strategy for the optimization of stable SCO based electrolytes for ultrahigh Seebeck thermal sensor. These materials are electrochemically reversible and hence, demonstrate a long lifetime use of these solutions in thermal sensor.

Conflicts of interest
There are no conicts to declare.