# Investigation of the Effect of Spin Crossover on the Static and Dynamic Properties of MEMS Microcantilevers Coated with Nanocomposite Films of [Fe(Htrz)2(trz)](BF4)@P(VDF-TrFE)

José Elías Angulo-Cervera, Mario Piedrahita-Bello, Fabrice Mathieu, Thierry Leichle, Liviu Nicu, Lionel Salmon, Gábor Molnár, Azzedine Bousseksou

► To cite this version:

José Elías Angulo-Cervera, Mario Piedrahita-Bello, Fabrice Mathieu, Thierry Leichle, Liviu Nicu, et al. Investigation of the Effect of Spin Crossover on the Static and Dynamic Properties of MEMS Microcantilevers Coated with Nanocomposite Films of [Fe(Htrz)2(trz)](BF4)@P(VDF-TrFE). Magnetochemistry, 2021, 7 (8), pp.114. 10.3390/magnetochemistry7080114. hal-03329761

HAL Id: hal-03329761

https://hal.science/hal-03329761

Submitted on 31 Aug 2021

HAL is a multi-disciplinary open access
archive for the deposit and dissemination of sci-
entific research documents, whether they are pub-
lished or not. The documents may come from
teaching and research institutions in France or
abroad, or from public or private research centers.

L'archive ouverte pluridisciplinaire HAL, est
destinée au dépôt et à la diffusion de documents
scientifiques de niveau recherche, publiés ou non,
émanant des établissements d'enseignement et de
recherche français ou étrangers, des laboratoires
publics ou privés.

<figure><img src="image_1.png" /></figure>

Distributed under a Creative Commons Attribution 4.0 International License
Article

# Investigation of the Effect of Spin Crossover on the Static and Dynamic Properties of MEMS Microcantilevers Coated with Nanocomposite Films of [Fe(Htrz)₂(trz)](BF₄)@P(VDF-TrFE)

José Elías Angulo-Cervera <fnref n="1" />,<fnref n="2" /> ID, Mario Piedrahita-Bello <fnref n="1" />,<fnref n="2" />, Fabrice Mathieu <fnref n="2" />, Thierry Leichle <fnref n="2" />,<fnref n="3" />, Liviu Nicu <fnref n="2" />, Lionel Salmon <fnref n="1" /> ID, Gábor Molnár <fnref n="1" />,* and Azzedine Bousseksou <fnref n="1" />,*

<figure><img src="image_2.png" /></figure>

check for
updates

**Citation:** Angulo-Cervera, J.E.; Piedrahita-Bello, M.; Mathieu, F.; Leichle, T.; Nicu, L.; Salmon, L.; Molnár, G.; Bousseksou, A. Investigation of the Effect of Spin Crossover on the Static and Dynamic Properties of MEMS Microcantilevers Coated with Nanocomposite Films of [Fe(Htrz)₂(trz)](BF₄)@P(VDF-TrFE). *Magnetochemistry* **2021**, *7*, 114.
https://doi.org/10.3390/magnetochemistry7080114

Academic Editor: Zoi G. Lada

Received: 28 June 2021

Accepted: 26 July 2021

Published: 9 August 2021

**Publisher's Note:** MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

<figure><img src="image_3.png" /></figure>

**Copyright:** © 2021 by the authors. Licensee MDPI, Basel, Switzerland. This article is an open access article distributed under the terms and conditions of the Creative Commons Attribution (CC BY) license (https://creativecommons.org/licenses/by/4.0/).

**Abstract:** We used a spray-coating process to cover silicon microcantilevers with ca. 33 wt% [Fe(Htrz)₂(trz)](BF₄)@P(VDF₇₀-TrFE₃₀) nanocomposite thin films of 1500 nm thickness. The bilayer cantilevers were then used to investigate the thermomechanical properties of the composites through a combined static and dynamic flexural analysis. The out-of-plane flexural resonance frequencies were used to assess the Young's modulus of the spray-coated films (3.2 GPa). Then, the quasi-static flexural bending data allowed us to extract the actuation strain (1.3%) and an actuation stress (7.7 MPa) associated with the spin transition in the composite.

**Keywords:** spin crossover; MEMS; mechanical properties

## 1. Introduction

The molecular spin crossover (SCO) phenomenon displayed by some transition metal complexes refers to the reversible switching between the low spin (LS) and high spin (HS) electronic configurations of the central metal ion—under the influence of an external stimulus (e.g., temperature, pressure, light irradiation, ... ) [1–4]. Due to the strong electron-lattice coupling, SCO materials are known to display a considerable volume change across the spin transition [5–7], which can be exploited in mechanical actuator devices [8–17]. In this context, an interesting asset of molecular SCO materials is that they can be operated at reduced size scales [18]. To explore the actuating properties of SCO complexes at the micro-nanometer scales, thin SCO films have been recently deposited onto moving parts of microelectromechanical systems (MEMS). The term of MEMS refers to a miniaturized device, which consists of both mechanical and electronic components with feature sizes at the microscale. In a broad sense, MEMS devices comprise microfabricated sensors, actuators and electronics that function as integrated systems. MEMS devices have significant impact in modern technology, for example in telecommunication and information technologies. Importantly for us, this class of sensitive micromechanical devices can provide quantitative information on the mechanical and actuating properties of SCO films.

Manrique et al. reported on the first SCO-MEMS [11], which consisted of a 200 nm film of the SCO complex [Fe(H₂B(pz)₂)₂(phen)] (pz = pyrazol-1-yl, phen = 1,10-phenanthroline) deposited by thermal evaporation on a silicon microcantilever. The latter was associated with a highly sensitive integrated piezoresistive detection system. The spin crossover was triggered by light irradiation (LIESST effect) at 40 K, which has led to a sizeable and

<footnote n="1">¹ Laboratoire de Chimie de Coordination, CNRS & University of Toulouse, 205 Route de Narbonne, F-31077 Toulouse, France; janguloc@laas.fr (J.E.A.-C.); mario.piedrahita@lcc-toulouse.fr (M.P.-B.); lionel.salmon@lcc-toulouse.fr (L.S.)</footnote>

<footnote n="2">² Laboratoire d'Analyse et d'Architecture des Systèmes, CNRS & University of Toulouse, 7 Avenue du Colonel Roche, F-31400 Toulouse, France; fmathieu@laas.fr (F.M.); leichle@laas.fr (T.L.); nicu@laas.fr (L.N.)</footnote>

<footnote n="3">³ Georgia Tech-CNRS International Research Laboratory, School of Electrical and Computer Engineering, Atlanta, GA 30332, USA</footnote>

<footnote>* Correspondence: gabor.molnar@lcc-toulouse.fr (G.M.); azzedine.bousseksou@lcc-toulouse.fr (A.B.)</footnote>
reversible shift in the resonance frequency of the cantilever. This resonance shift was quantitatively linked to the changes in thickness, density and Young's modulus of the film upon the SCO. In the continuity of these efforts, the microcantilevers were coated with a 140 nm-thin film of the [Fe(HB(tz)₃)₂] complex (tz = triazol-1-yl) [12]. Upon the SCO, changes in the resonance frequency and damping (dynamic response) as well as of the tip position (static response) of the cantilever could be clearly observed. In a follow-up work, these data were extracted and combined with mechanical measurements on bulk [Fe(HB(tz)₃)₂] samples, leading to the assessment of the full set of elastic constants of the compound [14]. In a similar approach, Urdampilleta and co-workers [15] drop-casted different SCO complexes over a polymer-based, soft MEMS cantilever, and witnessed an upshift of the resonance frequency when switching from the low to the high spin state. They attributed this behavior to surface stress effect. Finally, [Fe(Htrz)₂(trz)](BF₄) (Htrz = 1,2,4-4H-triazole, trz = 1,2,4-triazolato) coated silicon cantilevers were investigated by Dugay et al. [16], who reported on a thermal hysteresis of the cantilever resonance frequency, which closely matched the hysteresis loop associated with the SCO.

In addition to pure SCO materials, MEMS were also coated by SCO@polymer composite films. A dispersion of [Fe(Htrz)₂(trz)](BF₄) nanoparticles in the epoxy polymer SU8 was spray-coated over a silicon microcantilever [13]. Actuation induced by the thermal spin transition has led to a work density of 140 mJ/cm³ (see also Table 1), representing considerable performance metrics compared with other bending actuator technologies. The key interests of this SCO-composite spray-coating approach are: (1) the seamless integration of high topography micro/nanostructures, as well as (2) its versatility, enabling for deposition of various combinations of SCO particles and polymers. In the present work, we took advantage of the versatility of this method to spray-coat [Fe(Htrz)₂(trz)](BF₄) nanoparticles embedded in a poly(vinylidene fluoride—trifluoro-ethylene) P(VDF-TrFE) matrix over MEMS cantilevers. SCO@P(VDF-TrFE) composites have been recently investigated by some of us for the interesting electromechanical couplings between the components [19,20]. As a first step towards the exploration of these couplings at the micro- and nanometric scales, we present here an investigation of the thermomechanical properties of [Fe(Htrz)₂(trz)](BF₄)@P(VDF-TrFE) nanocomposite films, deposited on a MEMS cantilever.

<table><caption>**Table 1.** Mechanical and actuation properties of [Fe(Htrz)₂(trz)](BF₄)@polymer composite films.</caption><thead><tr><th></th><th>[Fe(Htrz)<sub>2</sub>(trz)](BF<sub>4</sub>)@P(VDF-TrFE)<br>[This Work]</th><th>[Fe(Htrz)<sub>2</sub>(trz)](BF<sub>4</sub>)@SU8 [13]</th></tr></thead><tbody><tr><td>Length, l (μm)</td><td>150</td><td>840</td></tr><tr><td>Si cantilever thickness, t<sub>s</sub> (μm)</td><td>2</td><td>20</td></tr><tr><td>Composite film thickness, t<sub>f</sub> (μm)</td><td>1.5</td><td>3.5</td></tr><tr><td>Width, w (μm)</td><td>50</td><td>100</td></tr><tr><td>Particle weight fraction, f (wt%)</td><td>33</td><td>30</td></tr><tr><td>Density, ρ (kg/m<sup>3</sup>)</td><td>1800</td><td>1400</td></tr><tr><td>Tip deflection at SCO δ (μm)</td><td>0.88</td><td>3.65</td></tr><tr><td>Cantilever curvature due to SCO k (m<sup>-1</sup>)</td><td>80</td><td>10</td></tr><tr><td>Actuation strain due to SCO, ε<sub>SCO</sub> (%)</td><td>1.3</td><td>1.0</td></tr><tr><td>Actuation stress due to SCO, σ<sub>SCO</sub> (MPa)</td><td>7.7</td><td>43 *</td></tr><tr><td>Young's modulus, E (GPa) at 40 °C</td><td>3.2</td><td>-</td></tr><tr><td>Young's modulus, E (GPa) at 95 °C</td><td>0.5</td><td>3.2</td></tr><tr><td>Work density due to SCO, w/v (mJ/cm<sup>3</sup>)</td><td>40</td><td>140</td></tr></tbody></table>

## 2. Results

In order to reduce particle aggregation effects during the spray-coating process, we opted for a reverse micelle-based particle synthesis method, which has been already described in [21]. The surfactant used in the synthesis enrobes the [Fe(Htrz)₂(trz)](BF₄)

<footnote>* Reevaluated (vs. [13]) using Equation (4).</footnote>
particles and reduces their aggregation [22]. TEM analysis of the particles revealed ca. 1000 nm long and 200 nm large rod-shaped particles. Films of the pure P(VDF-TrFE) copolymer (with 70%/30% monomer ratio) and nanocomposite films of 33 wt% [Fe(Htrz)₂(trz)](BF₄)@P(VDF-TrFE) with 1500 nm thickness were spray-coated on silicon microcantilevers, which were analyzed for their static and dynamical behavior both before and after the coating process.

Figure 1 shows a photograph and scanning electron microscopy (SEM) images of the SCO-MEMS device, revealing a smooth, continuous and defect free coating of the cantilevers. It is interesting to notice that the SEM images depict a homogeneous deposit even over the micro-wires that bond the MEMS system to the PCB (printed circuit board) and the PCB itself, proof of the capability of this method to deposit conformal thin films on three-dimensional structures.

<figure><img src="image_4.png" /><figcaption>**Figure 1.** (a) Photograph of the MEMS device. (b) SEM image of a silicon micro-cantilever coated by the [Fe(Htrz)₂(trz)](BF₄)@P(VDF-TrFE) composite film. (c,d) Larger scale SEM images showing the MEMS, the PCB and the bonding wires, which are also uniformly coated.</figcaption></figure>

MEMS can be operated in static or dynamic regime (or a combination of both) [23,24]. In the static mode, the MEMS is used to transform a static deformation induced by an external perturbation into a variation in electrical current or voltage (or vice-versa). In the dynamic mode, the cantilever is oscillated in resonance conditions and the MEMS is used to track changes induced by an external perturbation on the resonance frequency ($f_r$) and the damping of the oscillation. The damping is characterized by the quality factor (Q), which can be defined as the ratio of $f_r$ and the full width at half maximum (FWHM) of the resonance peak. Higher Q indicates a lower rate of energy loss (per oscillation cycle) and allows for higher sensitivity.

Figure 2a shows the change in the resonance curve of the micro-cantilevers due to the deposition of the 1.5 µm-thick composite film. One can depict a decrease in the resonance frequency by ca. 13 kHz and an important broadening of the peaks (i.e., an overall increase in Q). The resonance shift is mainly related to the addition of extra mass to the cantilever, whereas the change in the quality factor stems chiefly from a change in the elastic properties
of the resonator with the addition of the SCO@P(VDF-TrFE) coating. We also acquired resonance curves of the coated cantilever at 90 °C on heating as well as on cooling back from 130 °C. Since 90 °C falls in the middle of the SCO hysteresis (vide infra), the heating (and cooling) spectrum corresponds to the LS (resp. HS) state of the particles in the composite. As displayed in Figure 2b, the resonance is substantially downshifted in the HS state (by ca. 0.7 kHz). This finding is in agreement with previous results reported for the SCO-MEMS systems with [Fe(H₂B(pz)₂)₂(phen)] [11] and [Fe(HB(tz)₃)₂] [12] spin crossover coatings. In the next step, we carried out more than 100 thermal cycles with simultaneous recording of the resonance frequency of the coated cantilever. Following a small drift during the first few cycles, the thermomechanical behavior of the MEMS devices turned out remarkably robust and well-reproducible (Figure 2c). When the temperature increases (resp. decreases) between 40 °C and 130 °C the resonance frequency reversibly decreases (resp. increases) between 45,700 ± 100 Hz and 49,600 ± 100 Hz. It shall be noted that this reversibility of dynamical MEMS properties upon cycling refers mostly to the robustness of the thermomechanical properties of the composite. On the other hand, the reversibility of the SCO properties is more clearly inferred from the static MEMS data (vide infra).

<figure><img src="image_5.png" /><figcaption>**Figure 2.** (a) Resonance peak of a cantilever before and after coating with a 1.5 µm-thick [Fe(Htrz)₂(trz)](BF₄)@P(VDF-TrFE) composite film (T = 25 °C). (b) Resonance peak of the coated cantilever acquired in the LS and HS states within the hysteresis region (T = 90 °C). (c) Dynamic tracking of the MEMS resonance frequency for the coated cantilever upon 105 successive thermal cycles between 40–130 °C.</figcaption></figure>

The temperature-dependent static deflection of the cantilever coated with the nanocom-
posite material is compared with the behavior of the uncoated and P(VDF-TrFE) coated
cantilevers in Figure 3. The uncoated cantilever displays no substantial deflection in the
investigated temperature range (Figure 3a). As can be expected, the bilayer cantilever,
coated with the pure P(VDF-TrFE) film, exhibits a downward bending—corresponding to
an increase in the piezoresistance, *R* (Figure 3b). On heating, the cantilever coated with
the nanocomposite film exhibits first a similar, continuous downward bending due to the thermal expansion mismatch of the different layers (Figure 3c). However, one can depict a very clear hysteresis associated with the spin transition with an abrupt upward (downward) bending of the cantilever around 114 °C (83 °C) corresponding to the LS-to-HS (HS-to-LS) transitions. Overall, ten thermal cycles were carried out to verify the reproducibility of the SCO actuation in the static regime. After the first heating (“run-in”), the amplitude of actuation becomes stable and perfectly reversible. The spin transition temperatures (extracted from the derivative curves $dR/dT$) in these experiments are slightly upshifted when compared to the transition temperatures recorded by variable temperature UV absorption measurements on the same film (118 °C and 72 °C, see Figure 3d), which occurs most likely due to the different thermalization conditions in the MEMS and optical absorbance measurements (freestanding cantilever in 5 mbar vacuum vs. film on substrate in 1 atm ambient air).

<figure><img src="image_6.png" /><figcaption>Figure 3. Actuating response upon heating and cooling of (a) an uncoated silicon micro-cantilever, (b) a microcantilever coated with a pure P(VDF-TrFE) film, and (c) a microcantilever coated with a [Fe(Htrz)₂(trz)](BF₄)@P(VDF-TrFE) composite film. In (b,c) both the measured static piezoresistance change (left axis) and the calculated tip deflection (right axis) are given and the derivative curves ($dR/dT$ vs. $T$) are also plotted. (d) UV absorbance of a [Fe(Htrz)₂(trz)](BF₄)@P(VDF-TrFE) composite film recorded at 350 nm both for heating and cooling. The derivative curve ($dA/dT$ vs. $T$) is also shown.</figcaption></figure>

### 3. Discussion

We used the MEMS dynamical data primarily to evaluate the Young’s modulus of the spray-coated films. Indeed, the resonant angular frequency ($\omega$) is determined by the geometry of the resonator and the mechanical properties of the materials used. For our (idealized) bilayer, single-clamped cantilever, the ratio of the resonance frequencies of the coated ($\omega_{f+s}$) and uncoated ($\omega_s$) cantilevers is given by [25]:

$$ \left(\frac{\omega_{f+s}}{\omega_s}\right)^2 = \frac{(E_r t_r^3 + 1)(E_r t_r + 1) + 3E_r t_r (t_r + 1)^2}{(\rho_r t_r + 1)(E_r t_r + 1)} \quad (1) $$
with $E$, $t$ and $\rho$ the Young’s modulus, thickness and mass density, respectively. The subscript $r$ refers to the ratio between the film and substrate properties. As shown by Equation (1), the resonance shift (Figure 2a) is related primarily to changes in the Young’s modulus, mass density and thickness of the resonator. In these calculations, we considered that the mechanical properties of the silicon cantilever are constant with temperature, which is a reasonable approximation for this small temperature range [26]. The model parameters for the silicon cantilever were $t = 2 \ \mu m$, $E = 169 \ GPa$ and $\rho = 2330 \ kg/m^3$. As shown in Table 1 the extracted Young’s modulus ($E$) of the SCO@P(VDF-TrFE) composite film is ca. 3.2 GPa at 40 °C and 0.5 GPa at 95 °C, respectively. This marked softening shows the same tendency as reported in [27] for the pure P(VDF$_{70}$-TrFE$_{30}$). We shall note here that, in principle, the resonance frequency shift upon the SCO (Figure 2b) could be used to assess the SCO-induced changes in mechanical properties. In our composite films, however, the spin transition in the particles and the Curie transition in the matrix occurs in nearly the same temperature range ($T_{Curie} = 110 \ °C$ and $70 \ °C$ on heating and cooling, respectively [20]). As a result, the MEMS dynamical behavior, which is governed chiefly by the matrix properties, did not allow us to discriminate the SCO phenomenon properly.

On the other hand, it turned out that the Curie transition has a much smaller effect on the static MEMS properties (Figure 3b), which were thus used to assess the actuating properties associated with the SCO. The two key equations, which are widely used for the analysis of the (quasi-static) flexural bending of thin film—substrate, bilayer beam systems have been developed by Timoshenko [28] and Stoney [29]. Timoshenko’s theory was originally developed to model metal bimorphs that bend upon heating due to different coefficients of thermal expansion. However, it can be applied for other kinds of eigenstrains and can thus provide connection between the curvature of the bilayer cantilever ($k$) and the actuation strain $\varepsilon_{SCO}$, associated with the SCO [8]:

$$k = \frac{6(1 + m)^2}{h \left[ 3(1 + m^2) + (1 + mn) \left( m^2 + \frac{1}{mn} \right) \right]} \varepsilon_{SCO} \quad (2)$$

where $m = t_f/t_s$, $n = E_f/E_s$ and $h = t_f + t_s$. The subscripts $s$ and $f$ stand for substrate and film, whereas $t$ and $E$ denote thickness and Young’s modulus, respectively. For small tip deflections ($\delta$), the beam curvature can be simply estimated as:

$$k = \frac{2\delta}{l^2} \quad (3)$$

where $l$ is the beam length. Equations (2) and (3), together with the experimentally measured beam deflection (Figure 3c), beam geometry and Young’s modulus data (extracted from the dynamical data using Equation (1) can then be used to determine the actuation strain.

On the other hand, Stoney’s equation can be used to calculate the film stress $\sigma_f$ causing the beam curvature $k$ [29]:

$$\sigma_f = \frac{-E_s t_s^2}{6(1 - v_s) t_f} k \quad (4)$$

where $v$ is the Poisson’s coefficient ($v_{si} = 0.22$). This equation provides good numerical results for thin films ($t_f \ll t_s$) deposited on rectangular plates ($z \ll x, y$).

At this point, we must note that the cantilevers are actually formed by three layers. The silicon base is electrically insulated by an $SiO_2$ layer on top of which the SCO@P(VDF-TrFE) composite films are spray-coated. Metallic paths are also deposited on the $SiO_2$ layer. In addition, preferential orientation of the rod-shaped particles can be expected, leading to anisotropic properties. Hence, the calculation results reported below, which refer to a bilayer beam with a thin, isotropic deposit, should be considered only as a relatively crude estimation of the actual material properties.
For the spray-coated [Fe(Htrz)₂(trz)]BF₄)@P(VDF-TrFE) composite film, the inputs and outputs of the Timoshenko and Stoney equations are shown in Table 1. The spin transition gives rise to a reversible tip displacement of 0.88 µm, which corresponds to a curvature change of ca. 80 m⁻¹. This allows us to estimate the actuation strain, εₔ₅CO = 1.3%, the film stress associated with the SCO, σₔ₅CO = 7.7 MPa and the volumetric work density, w/v = Eεₔ₅CO²/2 = 40 mJ/cm³. It is instructive to compare these results with the precedent work of Manrique et al. [13], where the same SCO complex was used to spray-coat an SU8 epoxy-based composite on a MEMS cantilever with similar geometry. As shown in Table 1, we find similar actuation strain and lower values for the actuation stress and work density. The latter finding is obviously related to the fact that the Young's modulus of the epoxy composite is significantly higher (3.2 GPa) than that of the P(VDF-TrFE) (0.5 GPa) at the working temperature of the actuator.

**4. Materials and Methods**

[Fe(Htrz)₂(trz)](BF₄) particles were synthesized as described in [21]. Pellets of P(VDF-TrFE) were provided by Piezotech. A dispersion of the SCO particles (33 wt%) and P(VDF-TrFE) in 2-butanone was spray-coated on the MEMS cantilevers using a SUSS MicroTec AltaSpray 8 manual spray-coater. Before the coating process, the cantilevers were cleaned by means of a Jelight UVO-Cleaner 42 UV-Ozone cleaning instrument. Details of the MEMS fabrication and operation principles are given in refs. [23,24], whereas the cantilever geometries are shown in Table 1. Briefly, for the dynamic measurements, the device was actuated at its resonance frequency by a magnetic field and the mechanical stress induced by the vibration of the cantilever was detected by integrated piezoresistors. Static measurements were conducted either separately or simultaneously with the dynamic measurements using the same piezoresistances. To quantify the amplitude of the actuation, the piezoresistance signal was calibrated by an external optical vibrometer [13]. To detect the actuation associated with the SCO phenomenon, the MEMS were cycled between 40 °C and 130 °C using a HFS350V heating-cooling stage (Linkam Scientific). All experiments were conducted in a controlled vacuum pressure of 5 mbar.

**5. Conclusions**

A primary interest of our spray-coating approach is that the resulting films can be
easily incorporated into complex micro- and nanosystems. Indeed, the spray-coating
method allows depositing the films on practically any kind of surface, even if it displays
high topography features. This advantage was clearly demonstrated in the present paper
by spray-coating MEMS silicon micro-cantilevers with [Fe(Htrz)₂(trz)](BF₄)@P(VDF-TrFE)
composite films. The obtained bilayer cantilevers were then investigated for their dynamic
(vibration) and static (deflection) thermomechanical properties, which were found to
be highly reversible over more than 100 thermal cycles. From the dynamical data, we
extracted the Young's modulus of the spray-coated composite, which was comparable with
literature data. Most importantly, we evidenced a reproducible actuation around the spin
transition temperatures in the composite with an actuation strain of ~1.3%, an actuation
stress of ~7.7 MPa and a work density of ~40 mJ/cm³. These results provide proof for
the conservation of the attractive thermal actuating properties of the particles in the thin,
spray-coated films. Further work is planned to exploit this actuation in conjunction with
the piezoelectric properties of the matrix for thermal energy harvesting purposes.

**Author Contributions:** Conceptualization: A.B., methodology: T.L. and L.N., investigation: J.E.A.-C., F.M., M.P.-B., data analysis: J.E.A.-C., F.M., L.S., G.M., funding acquisition: L.S., A.B., supervision: G.M., L.N., T.L., A.B., writing—original draft: J.E.A.-C., G.M. All authors have read and agreed to the published version of the manuscript.

**Funding:** This work was supported by the CONACYT—Mexico (PhD grant 471690 of J.E.A.-C.) and the Occitanie Region/Federal University of Toulouse (PhD grant of MPB), the ANR-France (19-CE09-0008-01), and by the LAAS-CNRS technology platform, a member of Renatech network.
**Institutional Review Board Statement:** Not applicable.

**Informed Consent Statement:** Not applicable.

**Data Availability Statement:** Data is available from the authors on reasonable request.

**Conflicts of Interest:** The authors declare no conflict of interest.

## References

1. Gütlich, P.; Hauser, A.; Spiering, H. Thermal and Optical Switching of Iron(II) Complexes. *Angew. Chem. Int. Ed.* **1994**, *33*, 2024–2054. [CrossRef]

2. Gütlich, P.; Goodwin, H.A. *Spin Crossover in Transition Metal Compounds I–III*; Springer: Berlin/Heidelberg, Germany, 2004; Volume 233–235.

3. Bousseksou, A.; Molnar, G.; Salmon, L.; Nicolazzi, W. Molecular spin crossover phenomenon: Recent achievements and prospects. *Chem. Soc. Rev.* **2011**, *40*, 3313–3335. [CrossRef] [PubMed]

4. Halcrow, M.A. (Ed.) *Spin-Crossover Materials: Properties and Applications*; John Wiley & Sons: Oxford, UK, 2013.

5. Collet, E.; Guionneau, P. Structural analysis of spin-crossover materials: From molecules to materials. *Comptes Rendus Chim.* **2018**, *21*, 1133–1151. [CrossRef]

6. Guionneau, P. Crystallography and spin-crossover. A view of breathing materials. *Dalton Trans.* **2014**, *43*, 382–393. [CrossRef] [PubMed]

7. Halcrow, M. Structure: function relationships in molecular spin-crossover complexes. *Chem. Soc. Rev.* **2011**, *40*, 4119–4142. [CrossRef]

8. Shepherd, H.J.; Gural'Skiy, I.A.; Quintero, C.M.; Tricard, S.; Salmon, L.; Molnar, G.; Bousseksou, A. Molecular actuators driven by cooperative spin-state switching. *Nat. Commun.* **2013**, *4*, 2607. [CrossRef]

9. Gural'Skiy, I.A.; Quintero, C.M.; Costa, J.S.; Demont, P.; Molnar, G.; Salmon, L.; Shepherd, H.; Bousseksou, A. Spin crossover composite materials for electrothermomechanical actuators. *J. Mater. Chem. C* **2014**, *2*, 2949–2955. [CrossRef]

10. Chen, Y.-C.; Meng, Y.; Ni, Z.-P.; Tong, M.-L. Synergistic electrical bistability in a conductive spin crossover heterostructure. *J. Mater. Chem. C* **2014**, *3*, 945–949. [CrossRef]

11. Manrique-Juarez, M.D.; Rat, S.; Mathieu, F.; Saya, D.; Séguy, I.; Leïchlé, T.; Nicu, L.; Salmon, L.; Molnar, G.; Bousseksou, A. Microelectromechanical systems integrating molecular spin crossover actuators. *Appl. Phys. Lett.* **2016**, *109*, 061903. [CrossRef]

12. Manrique-Juarez, M.D.; Mathieu, F.; Shalabaeva, V.; Cacheux, J.; Rat, S.; Nicu, L.; Leïchlé, T.; Salmon, L.; Molnár, G.; Bousseksou, A. A Bistable Microelectromechanical System Actuated by Spin-Crossover Molecules. *Angew. Chem. Int. Ed.* **2017**, *56*, 8074–8078. [CrossRef]

13. Manrique-Juárez, M.D.; Mathieu, F.; Laborde, A.; Rat, S.; Shalabaeva, V.; Demont, P.; Thomas, O.; Salmon, L.; Leïchlé, T.; Nicu, L.; et al. Micromachining-Compatible, Facile Fabrication of Polymer Nanocomposite Spin Crossover Actuators. *Adv. Funct. Mater.* **2018**, *28*, 1801970. [CrossRef]

14. Mikolasek, M.; Manrique-Juarez, M.D.; Shepherd, H.J.; Ridier, K.; Rat, S.; Shalabaeva, V.; Bas, A.-C.; Collings, I.E.; Mathieu, F.; Cacheux, J.; et al. Complete Set of Elastic Moduli of a Spin-Crossover Solid: Spin-State Dependence and Mechanical Actuation. *J. Am. Chem. Soc.* **2018**, *140*, 8970–8979. [CrossRef]

15. Urdampilleta, M.; Ayela, C.; Ducrot, P.-H.; Rosario-Amorin, D.; Mondal, A.; Rouzières, M.; Dechambenoit, P.; Mathonière, C.; Mathieu, F.; Dufour, I.; et al. Molecule-based microelectromechanical sensors. *Sci. Rep.* **2018**, *8*, 1–6. [CrossRef]

16. Dugay, J.; Giménez-Marqués, M.; Venstra, W.J.; Torres-Cavanillas, R.; Sheombarsing, U.N.; Manca, N.; Coronado, E.; Van Der Zant, H.S.J. Sensing of the Molecular Spin in Spin-Crossover Nanoparticles with Micromechanical Resonators. *J. Phys. Chem. C* **2019**, *123*, 6778–6786. [CrossRef]

17. Piedrahita-Bello, M.; Angulo-Cervera, J.E.; Courson, R.; Molnár, G.; Malaquin, L.; Thibault, C.; Tondu, B.; Salmon, L.; Bousseksou, A. 4D printing with spin-crossover polymer composites. *J. Mater. Chem. C* **2020**, *8*, 6001–6005. [CrossRef]

18. Molnár, G.; Rat, S.; Salmon, L.; Nicolazzi, W.; Bousseksou, A. Spin Crossover Nanomaterials: From Fundamental Concepts to Devices. *Adv. Mater.* **2018**, *30*, 1703862. [CrossRef]

19. Rat, S.; Piedrahita-Bello, M.; Salmon, L.; Molnár, G.; Demont, P.; Bousseksou, A. Coupling Mechanical and Electrical Properties in Spin Crossover Polymer Composites—Supporting information. *Adv. Mater.* **2018**, *30*, 1705275. [CrossRef]

20. Piedrahita-Bello, M.; Martin, B.; Salmon, L.; Molnár, G.; Demont, P.; Bousseksou, A. Mechano–Electric coupling in P(VDF–TrFE)/spin crossover composites. *J. Mater. Chem. C* **2020**, *8*, 6042–6051. [CrossRef]

21. Suleimanov, I.; Costa, J.S.; Molnár, G.; Salmon, L.; Fritsky, I.; Bousseksou, A. Effect of ligand substitution in [Fe(H-trz)2(trz)]BF4 spin crossover nanoparticles. *Fr. Ukr. J. Chem.* **2015**, *3*, 66–72. [CrossRef]

22. Létard, J.-F.; Guionneau, P.; Goux-Capes, L. Towards Spin Crossover Applications. *Top. Curr. Chem.* **2006**, *235*, 221–249. [CrossRef]

23. Alava, T.; Mathieu, F.; Mazenq, L.; Soyer, C.; Remiens, D.; Nicu, L. Silicon–Based micromembranes with piezoelectric actuation and piezoresistive detection for sensing purposes in liquid media. *J. Micromech. Microeng.* **2010**, *20*, 075014. [CrossRef]

24. Dufour, I.; Maali, A.; Amarouchene, Y.; Ayela, C.; Caillard, B.; Darwiche, A.; Guirardel, M.; Kellay, H.; Lemaire, E.; Mathieu, F.; et al. The Microcantilever: A Versatile Tool for Measuring the Rheological Properties of Complex Fluids. *J. Sens.* **2012**, *2012*, 719898. [CrossRef]
25. Whiting, R.; Angadi, M.; Tripathi, S. Evaluation of elastic moduli in thin-film/substrate systems by the two-layer vibrating reed method. *Mater. Sci. Eng. B* **1995**, *30*, 35–38. [CrossRef]

26. Vanhellemont, J.; Swarnakar, A.K.; Van Der Biest, O. Temperature Dependent Young's Modulus of Si and Ge. *ECS Trans.* **2014**, *64*, 283–292. [CrossRef]

27. Hafner, J.; Teuschel, M.; Schneider, M.; Schmid, U. Origin of the strong temperature effect on the piezoelectric response of the ferroelectric (co-)polymer P(VDF70-TrFE30). *Polymer* **2019**, *170*, 1–6. [CrossRef]

28. Timoshenko, S. Analysis of Bi-Metal Thermostats. *J. Opt. Soc. Am.* **1925**, *11*, 233–255. [CrossRef]

29. Stoney, G.G. The tension of metallic films deposited by electrolysis. *Proc. R. Soc. London. Ser. A Math. Phys. Sci.* **1909**, *82*, 172–175. [CrossRef]