Emerging haptic technology based on piezoelectric actuators enables to realize innovative tactile human–machine interface. The standard solution is based on stand-alone bulk ceramics glued directly on the haptic device. Thin-film actuators with metal–insulator–metal structure have been developed to directly integrate actuators on haptic plates. The thickness of thin films is limited to 2 μm, leading to large capacitance and, thus, too high-power consumption. To solve this issue, we developed haptic devices based on a 10 μm-thick PZT film deposited on a 0.65 mm-thick platinized silicon substrate. These thick films are made of a PZT composite slurry associated with sol-gel sol infiltration. They are dense and exhibit a permittivity of 1000 and dielectric loss lower than 0.05. Our fabricated haptic device containing three actuators connected in series exhibits an antisymmetric Lamb wave resonant mode at 62.0 kHz, in line with finite element modeling. At the limit of touch detection (1 μm out of plane deflection), the power consumption of the haptic device is 150 mW at 40 V. This represents a 15-fold consumption reduction with respect to the same haptic device made with 0.5 μm-thick PZT thin films.

Although human–machine interface (HMI) strongly depends on vision and hearing, it calls for adding the sense of touch.1,2 Hence, surface haptic based on piezoelectric actuators is developed to generate tactile effect on surfaces.3–5 Piezoelectric actuators can induce an ultrasonic wave at certain frequencies.6–8 A squeezed-air film appears between the vibrating surface and a sliding finger, which in turn decreases the coefficient of friction.9–11 This enables real-time control of the friction coefficient of someone's finger, thus emulating artificial textures.

Bulk piezo-ceramic haptic actuators stand for the nominal technological solution. They enable large deflections and are typically made of PZT ceramics with thickness in the mm-range,12,13 which require large voltage, typically beyond 100 V, and which prevents using collective fabrication means, inducing fabrication costs difficult to overcome.

One alternative enabling acceptable deflection and technology compatible with mass production is based on piezoelectric thin films.14,15 However, the main drawback of piezoelectric thin film actuators for haptic is large power consumption induced by their large capacitance. For example, Hua et al. recently demonstrated a transparent haptic device based on 0.5 μm-ITO/2 μm-PZT/0.5 μm-ITO metal–insulator–metal (MIM) stack on Corning Lotus™ NXT Glass (71 × 60 × 0.5 mm3).15 The haptic device showed an out-of-plane deflection of 1 μm at 10 Vpp at 22.5 kHz. The power consumption of the whole haptic device reached 3.6 W.

Efforts have nonetheless been devoted to developing power-efficient devices. Glinsek et al. demonstrated a transparent haptic device based on a PZT thin film and ITO interdigitated (IDE) electrodes on fused silica glass.16 The device exhibited a Lamb mode at 73 kHz and an out-of-plane peak-to-peak displacement of 3 μm at 150 V unipolar. The benefit of the IDE structure is a smaller capacitance than MIM (120 pF for each capacitor in this case), which therefore leads to much lower power consumption, namely, 75 mW per actuator (at 150 V unipolar = 92 Vrms). Note though that the footprint of this device is very much reduced (15 × 3 × 0.5 mm3). A simple linear increase in the width of the actuators, in order to make a fair comparison with the previous work, would multiply by 20 the value of the capacitance, reaching then 1.5 W per actuator. These examples not only show that consumption is substantial in various haptic devices from the literature but also that it is not straightforward to come up with a fair comparison between designs.

To obtain low-power haptic devices, our idea is to develop piezoelectric thick films actuated with a MIM structure. This enables decreasing the capacitance value (C) while keeping a similar out of plane deflection at a given voltage, as explained hereafter. The 1D-lateral force F induced by voltage U applied to a piezoelectric layer on a plate is17 

F=e31,fUtfilmAfilm,
(1)

where e31,f is the effective piezoelectric coefficient, tfilm is the film thickness, Afilm is the cross-sectional area of the film, which equals film thickness times film width (i.e., Afilm = tfilm × wfilm). Equation (1), therefore, simply reads F = e31,f × U × wfilm. This indicates that the induced force does not depend on PZT thickness. As will be shown in the finite element modeling (FEM) section, the increase in the structure's stiffness remains negligible if the film's thickness remains within 10% of the substrate's thickness. In addition, increasing the film's thickness proportionally decreases C, the capacitance of the actuator. As the power consumption is proportional to the capacitance value, this in turn linearly decreases the overall consumption of the haptic device. Indeed, power consumption P equals C × f × Urms2, with Urms the applied voltage at frequency f.

In this work, we first investigated the impact of PZT film thickness on the performance of a simple haptic device using FEM. It reveals that the device deflection is hardly ever impacted by PZT thickness. A thicker film enables decreasing power consumption as discussed before. Based on this analysis, we fabricated two haptic transducers with identical dimensions but actuated, respectively, with thin (0.5 μm) and thick (10 μm) Pb(Zr0.53Ti0.47)O3 (PZT) films. The thin and thick films were based, on the one hand, on a standard sol-gel approach and, on the other hand, on a composite slurry coupled with an infiltration sol-gel process. The fabricated haptic devices moved by three actuators each exhibited standing Lamb waves at 62.0 kHz (thick films) and 71.5 kHz (thin films). The capacitance of one thin film actuator is 160 nF whereas that of one thick film actuator is only 12 nF, both measured at 1 kHz. When 13 Vpp (peak to peak) were applied to each actuator, the thin and thick film haptic devices exhibited, respectively, a peak-to-peak deflection of 0.4 and 1 μm (limit of detection with a finger) and a consumption of 2.3 W (thin film) and 150 mW (thick film). This experiment, therefore, proves that increasing the thickness of PZT films in haptic devices decreases drastically the power consumption and not at the expense of reducing the haptic effect.

A 2D FEM was performed using COMSOL 5.5 to study the relationship between PZT thickness and deflection of a simple haptic plate. A scheme of the haptic device (34 × 10 × 0.65 mm3) is shown in Fig. 1(a). Three PZT capacitors (2 × 10 mm2) are placed symmetrically on the silicon substrate with spacing of 11.3 mm between each other. Electrodes, passivation, and barrier layers were not modeled because they are too thin to influence the device performance. A voltage of 40 V was applied to each PZT actuator. The piezoelectric coefficient e31 and relative permittivity εr were set to −7 C m−2 and 1000, respectively. The damping loss factor calculated as Δf/fr, where fr and Δf are resonant frequency and resonance width at full width at half maximum, respectively, is 0.0063 for 0.65 mm-thick platinized silicon. Figure 1(b) shows the modeled maximum deflection of the haptic device from Fig. 1(a) and each actuator's capacitance as a function of PZT thickness. A linear decrease in the capacitance is observed with PZT thickness, meaning that 1, 10, and 100 μm-thick films, respectively, exhibit capacitances of 180, 18, and 1.8 nF. As expected, the corresponding peak displacements vs PZT thickness of the haptic device spread in a much smaller range, varying from 1.6 to 1.8 μm only. The displacement is slightly improved when PZT increases because adding PZT has a softening effect, visible on the acoustic resonance. Indeed, the frequency of the latter decreases when PZT thickness increases, as displayed in Fig. 1(c). It makes sense because PZT is softer and heavier than Si. Figure 1(d) represents the mode shape of a 10 μm-PZT haptic device at its A0-Lamb wave mode resonance (66.1 kHz) where the maximum deflection appears.

FIG. 1.

Finite element modeling. (a) Geometry of the model built with three PZT capacitors on a 0.65 mm-thick silicon substrate. (b) The displacement of the haptic device (black symbol) and the capacitance of each actuator (red symbol) as a function of PZT thickness. (c) Resonant frequency plotted as a function of PZT thickness. (d) Modeled shape of haptic device based on a 10 μm-thick PZT film.

FIG. 1.

Finite element modeling. (a) Geometry of the model built with three PZT capacitors on a 0.65 mm-thick silicon substrate. (b) The displacement of the haptic device (black symbol) and the capacitance of each actuator (red symbol) as a function of PZT thickness. (c) Resonant frequency plotted as a function of PZT thickness. (d) Modeled shape of haptic device based on a 10 μm-thick PZT film.

Close modal

In this part, only the thick films fabrication process is detailed. The thin film process is very well established, for instance, in Ref. 16 and the details are given in the supplementary material. The comparison between thick and thin film haptic devices will though be performed in the main paper.

PZT thick films were fabricated by the deposition of a composite slurry (1 g/mL), which was synthesized by dispersing crystalline Nb-PZT [Pb(Zr0.53Ti0.47)0.98Nb0.02O3] powder into PZT [Pb(Zr0.53Ti0.47)O3] chemical solution deposition (CSD) solution. More details on processing and characterization are provided in the supplementary material. The slurry was deposited on platinized silicon via spin coating and was annealed at 700 °C to sinter the powders and crystallize the amorphous deposit. The deposition-heating-crystallization process was repeated up to seven times to obtain the desired thickness. The cross section scanning electron microscopy (SEM) micrograph of a four-layer-thick film reveals a porous microstructure [Fig. 2(a)] in which individual deposited layers can be observed. Its thickness is 6 μm.

FIG. 2.

Cross-sectional SEM micrographs of (a) a PZT composite thick film after four depositions without infiltration. Sol infiltrated PZT composite thick films with four (b) and seven (c) depositions. (d) The XRD pattern of the 10 μm-thick PZT film shown in (c). Peaks are denoted with pseudo-cubic Miller indices corresponding to the PZT perovskite phase.19 

FIG. 2.

Cross-sectional SEM micrographs of (a) a PZT composite thick film after four depositions without infiltration. Sol infiltrated PZT composite thick films with four (b) and seven (c) depositions. (d) The XRD pattern of the 10 μm-thick PZT film shown in (c). Peaks are denoted with pseudo-cubic Miller indices corresponding to the PZT perovskite phase.19 

Close modal

To improve the density of the films, a so-called “x(1C + 4S)” procedure18 was deployed, in which one composite slurry deposition (1C) is followed by four CSD solution depositions (4S) to enable infiltration of the sol into the composite. x denotes the number of deposition-infiltration cycles. Details about this process can be found in the supplementary material (Fig. S2). The cross section of 4(1C + 4S) films displays a much denser microstructure [Fig. 2(b)] than the sample without infiltration [Fig. 2(a)]. 10 μm-thick films were achieved with a 7(1 C + 4S) procedure [Fig. 2(c)].

The x-ray diffraction (XRD) pattern of the 10 μm-thick film with sol-infiltration [7(1C + 4S)] is shown in Fig. 2(d). A pure perovskite phase is identified using the powder diffraction file (PDF) No. 01–070-4264.19 The presence and intensity of {100}, {110}, and {200} peaks indicate random orientation, which is a typical feature of nanoparticles-based thick films.20 

Figure 3(a) shows the polarization-electric field [P(E)] loops of 6 μm-thick PZT film deposited without infiltration. The remanent polarization (Pr) is 25 μC cm−1, and the coercive field (Ec) is 125 kV cm−1. Ec is rather large and linearly shifts as voltage increases. It suggests substantial leakage current as observed in the current density j(E) loops [Fig. 3(b)]. εr and tanδ are, respectively, 400 and 0.2, as shown in Fig. 3(c).

FIG. 3.

Electrical characterizations of PZT thick films (I) deposition of PZT slurry without infiltration (6 μm-thick) and (II) with infiltration (10 μm-thick). (a) and (d) polarization-electric field P(E) hysteresis loops, (b) and (e) current density-electric field j(E) curves, and (c) and (f) permittivity εr (black curve) and loss tangent tanδ (red curves) measurements. Electrode area is 3.14 mm2.

FIG. 3.

Electrical characterizations of PZT thick films (I) deposition of PZT slurry without infiltration (6 μm-thick) and (II) with infiltration (10 μm-thick). (a) and (d) polarization-electric field P(E) hysteresis loops, (b) and (e) current density-electric field j(E) curves, and (c) and (f) permittivity εr (black curve) and loss tangent tanδ (red curves) measurements. Electrode area is 3.14 mm2.

Close modal

Figure 3(d) shows the P(E) loops of 10 μm-thick PZT films deposited by the 7(1 C + 4S) sol-infiltration process. Pr is 27 μC cm−1 and Ec is 42 kV cm−1, values that are standard for PZT films. Ec is divided by three after sol-infiltration. It suggests lower leakage current, which is confirmed in the j(E) loops of Fig. 3(e). At high field, current is indeed much smaller. The j(E) loops also show two sharp peaks in current that are linked with the switching of PZT ferroelectric domains, as expected. εr is 1000 and tanδ is 0.05 [Fig. 3(f)], values comparable to the values of our sol-gel deposited PZT thin films (εr = 900 and tanδ = 0.08, as shown in Fig. S5).

Two slider haptic devices with three actuating areas each were fabricated according to the modeling [cf. Fig. 1(a)], one based on 0.5 μm-thick sol gel films (for comparison) and one with 10 μm-thick sol-infiltrated PZT films. Both devices have the same dimensions, namely, 34 × 10 × 0.65 mm3.

The three PZT actuators were obtained by depositing and patterning, respectively, 100 nm-thick and 500 nm-thick Pt top electrodes of 2 × 8 mm2 on PZT thin and thick films. In Fig. 4(b), the P(E) loop of the marked actuator confirms its ferroelectric nature with Pr = 15 μC cm−2. The corresponding current curve does not display excessive leakage at high field [Fig. 4(c)]. In addition, the C(E) loop exhibits an expected value of capacitance around 12 nF per actuator [cf. Fig. 4(d)]. The losses are around 5% at maximum, which is a rather low value. The respective curves for the thin film actuators are displayed in Fig. 4(d) [P(E) loop] and Fig. 4(e) [C(E) loop]. The polarization is very similar to the one obtained for thick films. However, the capacitance is much larger, as expected, reaching 160 nF at maximum, with associated dielectric losses around 8%.

FIG. 4.

(a) Top-view and (b) cross-section schematic diagram of fabricated haptic devices based on a 10 μm-thick or 0.5 μm-thick PZT film. (b) and (d) P(E) hysteresis loop and (c) and (f) C(E) and tanδ of one of the fabricated actuators based on thick and thin films, respectively.

FIG. 4.

(a) Top-view and (b) cross-section schematic diagram of fabricated haptic devices based on a 10 μm-thick or 0.5 μm-thick PZT film. (b) and (d) P(E) hysteresis loop and (c) and (f) C(E) and tanδ of one of the fabricated actuators based on thick and thin films, respectively.

Close modal

The observed resonant frequencies of the thin and thick film actuators are, respectively, 71.5 and 62.0 kHz, corresponding to the frequency at which the deflection is the largest. Note that these values are not exactly the same as the ones predicted by modeling (67 kHz for the thin films and 66 kHz for the thick films) because the actual final dimensions, and especially the lengths, was slightly different from initial design due to imperfect cut.

The out-of-plane displacement of the thick film haptic device was recorded by scanning across the area of the device, as plotted in Fig. 5(a). 120 V peak-to-peak and unipolar (which means 60 V DC + 60 V AC) were applied to the three actuators in series, which corresponds to 40 V applied to each actuator. A typical anti-symmetric Lamb wave A0 was observed, which agrees well with the mode shape predicted by FEM. The maximum displacement reaches 3.0 μm peak-to-peak. We also observed that only 13.3 V is needed to achieve the minimum deflection that can be detected by a finger, namely 1 μm. A 2D-scan over the entire device surface [Fig. 5(b)] of the deflection at resonance has also been performed. The same Lamb stationary wave with four nodes equally spaced along the length was observed. Note also that the deflection along the y-axis stays unchanged as expected. The thin film haptic device could only withstand 20 V peak-to-peak and unipolar (which means 10 V DC + 10 V AC) for the three capacitors in series, which means 6.7 V per each actuator. The maximum deflection at the Lamb wave resonance at 71.5 kHz only reached 0.2 μm peak-to-peak, as shown in Fig. S7. The deflection is proportional to the applied voltage, so peak-to-peak deflection at 13.3 V should be 0.4 μm, which is 2.5 times less than the thick film device.

FIG. 5.

(a) Line scan along the haptic device length and (b) 2D-map measured atop the PZT thick films haptic device when 120 V unipolar were applied to the three actuators in series, at its A0 Lamb wave mode resonance frequency (62.0 kHz).

FIG. 5.

(a) Line scan along the haptic device length and (b) 2D-map measured atop the PZT thick films haptic device when 120 V unipolar were applied to the three actuators in series, at its A0 Lamb wave mode resonance frequency (62.0 kHz).

Close modal

Note that increasing the film thickness could induce a change in piezoelectric properties due to stress relaxation or clamping effect, as described in Ref. 21. In our specific case, e31,f values obtained by fitting the model to experimental deflection were −7 and −5 C m−2 for thick and thin films, respectively. Therefore, there is a significant increase in piezoelectric properties in our thick films (+40%). The converse e31,f coefficient of the PZT thick film was also experimentally measured through a cantilever structure.22 Its value is −7.6 C m−2, which is in line with the value extracted from the haptic experiment. All the details are described in the supplementary material.

As exposed earlier, the purpose of developing PZT thick film-based haptic devices is to reduce the power consumption. The power consumption P to drive one capacitor can be estimated as P = CfUrms2, where C, Urms, and f are capacitance, rms voltage, and resonant frequency, respectively. In order to perform a fair comparison, we applied the same voltage to the thin and the thick film haptic devices, namely, 40 V unipolar (20 V DC + 20 V AC, at 24 Vrms). The results are displayed in Table I and show that one thin film actuator consumes 750 mW, whereas one thick film actuator consumes 50 mW. Therefore, the haptic device with thick films consumes 15 times less than the one with thin films. As shown before, the thick film haptic device can withstand 120 V with the actuators in series. In this case, the consumption of all actuators is 450 mW. When 40 V (unipolar, peak-to-peak) are applied to the actuator at the resonance, the directly measured power consumption (Urms × Irms) is 500 mW, which is in line with the value (450 mW) obtained through Cf(Urms)2. More measurement details are described in the supplementary material.

TABLE I.

Comparison of 10 μm-PZT thick film-based haptic device fabricated in this work with 0.5 μm-PZT thin film-based haptic device. For a fair comparison, the haptic devices have the same dimensions. The same voltage has been applied to both (40 V unipolar, at 24 Vrms) at their respective resonant frequency.

CharacteristicsThick film haptic deviceThin film haptic device
Permittivity of PZT measured at 1 kHz 1000 900 
Dielectric loss of PZT 0.05 0.08 
Remanent polarization of PZT 27 μC cm−2 15 μC cm−2 
Haptic device frequency 62.0 kHz 71.5 kHz 
Capacitance of one actuator 12 nF 160 nF 
Peak-to-peak deflection 1 μ0.4 μ
Power consumption of one actuator (of the three actuators) 50 mW (150 mW) 750 mW (2.3 W) 
CharacteristicsThick film haptic deviceThin film haptic device
Permittivity of PZT measured at 1 kHz 1000 900 
Dielectric loss of PZT 0.05 0.08 
Remanent polarization of PZT 27 μC cm−2 15 μC cm−2 
Haptic device frequency 62.0 kHz 71.5 kHz 
Capacitance of one actuator 12 nF 160 nF 
Peak-to-peak deflection 1 μ0.4 μ
Power consumption of one actuator (of the three actuators) 50 mW (150 mW) 750 mW (2.3 W) 

Finally, we can also compare approximately these values with the ones reported in the literature in the introduction by normalizing the consumption with the actuators' width. Hence, if all actuators were 10 mm in width and if they were driven in order to reach 3–4 μm in deflection, Hua's15 and Glinsek's16 would, respectively, consume 3.5 W and 250 mW. In these conditions, our thick film haptic device consumes 450 mW, which is much less than Hua's because our capacitance is smaller. It is though twice as much as Glinsek's but our device enables working at lower voltage (40 Vpp per actuator vs 150 Vpp). Note also that we used a 0.65 mm-thick Si substrate as haptic plate, which is more difficult to deflect than the glass substrate used in Ref. 16 (0.5 mm thick). At constant voltage, modeling shows that the deflection of a haptic device with 0.5 mm-glass is four times larger than the same one with 0.65 mm-Si. Therefore, this indicates that our thick film actuators should only need 10 Vpp to actuate the same glass haptic plate with the same dimension to reach a deflection of 3 μm. Our solution clearly calls for even thicker films, as suggested by our modeling in Fig. 1(b).

As suggested by simple equations, finite element modeling demonstrated that using thick piezoelectric films for haptic devices enables decreasing power consumption. Dense 10 μm-thick PZT films were fabricated. They were based on a PZT composite slurry complemented with a sol-infiltration process. Thanks to their dense microstructure, these PZT thick films exhibit ferroelectric properties comparable with state-of-art sol-gel PZT thin films. Two haptic devices were fabricated, respectively, based on PZT thin and thick films. Both are functional and exhibit an A0 Lamb wave resonant mode at 71.5 and 62.0 kHz, respectively. The power consumption per thick-film actuator is 50 mW for generating a 1 μm out-of-plane displacement, which represents a 15-fold reduction compared with our 500 nm-PZT thin film haptic device. Our PZT deposition technique enables reaching 10 μm in thickness. Modeling clearly shows that reaching 100 μm in PZT thickness would allow the same deflection with the same voltage, but this would also reduce consumption by one order of magnitude, down to a few mW per actuator.

See the supplementary material for further details about processing of solutions, powders, slurries, and films. Microstructural and electrical characterization of the films and haptic devices is also further described.

The authors acknowledge the Fonds National de la Recherche (FNR) of Luxembourg for supporting this work through the Project No. PRIDE17/12246511/PACE. S.D. and B.M. acknowledge the support of Slovenian Research Agency (core funding P2-0105).

The authors have no conflicts to disclose.

Longfei Song: Data curation (equal); Investigation (equal); Writing – original draft (equal); Writing – review & editing (equal). Sebastjan Glinsek: Conceptualization (equal); Methodology (equal); Supervision (supporting); Validation (equal); Writing – review & editing (equal). Silvo Drnovsek: Data curation (supporting); Investigation (supporting). Veronika Kovacova: Investigation (supporting). Barbara Malic: Investigation (supporting); Writing – review & editing (supporting). Emmanuel Defay: Conceptualization (equal); Supervision (equal); Validation (equal); Writing – review & editing (lead).

The data that support the findings of this study are available from the corresponding authors upon reasonable request.

1.
S.
Standring
,
Gray's Anatomy E-Book: The Anatomical Basis of Clinical Practice
(
Elsevier Health Sciences
,
Philadelphia, PA
,
2015
).
2.
I.
Poupyrev
,
S.
Maruyama
, and
J.
Rekimoto
, in
Proceedings of the 15th Annual ACM Symposium on User Interface Software and Technology
,
2002
.
3.
T. H.
Yang
,
J. R.
Kim
,
H.
Jin
,
H.
Gil
,
H.
Koo
, and
H. J.
Kim
,
Adv. Funct. Mater.
31
,
2008831
(
2021
).
4.
C.
Basdogan
,
F.
Giraud
,
V.
Levesque
, and
S.
Choi
,
IEEE Trans. Haptics
13
,
450
(
2020
).
5.
Y.
Huang
,
K.
Yao
,
J.
Li
,
D.
Li
,
H.
Jia
,
Y.
Liu
,
C. K.
Yiu
,
W.
Park
, and
X.
Yu
,
Mater. Today Phys.
22
,
100602
(
2022
).
6.
F.
Bernard
,
F.
Casset
,
J.
Danel
,
C.
Chappaz
, and
S.
Basrour
,
J. Micromech. Microeng.
26
,
084007
(
2016
).
7.
M.
Wiertlewski
and
Colgate
,
J. E.
,
IEEE Trans. Haptics
8
,
43
(
2015
).
8.
F.
Casset
,
J.-S.
Danel
,
P.
Renaux
,
C.
Chappaz
,
F.
Bernard
,
T.
Sednaoui
,
S.
Basrour
,
B.
Desloges
, and
S.
Fanget
,
Mechatronics
40
,
264
(
2016
).
9.
M.
Biet
,
F.
Giraud
, and
B.
Lemaire-Semail
,
IEEE Trans. Ultrason, Ferroelect, Freq. Contr.
54
,
2678
(
2007
).
10.
M.
Wiertlewski
,
R. F.
Friesen
, and
J. E.
Colgate
,
Proc. Natl. Acad. Sci.
113
,
6907
(
2016
).
11.
T.
Watanabe
and
S.
Fukui
, in
IEEE International Conference on Robotics and Automation
(
IEEE
,
Piscataway, NJ
,
1995
), pp.
1134
1139
.
12.
See https://www.hap2u.net/haptic-technology/ for further information on piezoelectric haptics.
13.
See https://product.tdk.com/en/products/sw_piezo/haptic/piezohapt/index.html for further information on piezoelectric haptics.
14.
S.
Glinsek
,
L.
Song
,
V.
Kovacova
,
M. A.
Mahjoub
,
N.
Godard
,
S.
Girod
,
J.-L.
Biagi
,
R.
Quintana
,
T.
Schleeh
,
M.
Guedra
,
M.
Rupin
, and
E.
Defay
,
Adv. Mater. Technol.
7
,
2200147
(
2022
).
15.
H.
Hua
,
Y.
Chen
,
Y.
Tao
,
D.
Qi
, and
Y.
Li
,
Sens. Actuator, B
335
,
113396
(
2022
).
16.
S.
Glinsek
,
M. A.
Mahjoub
,
M.
Rupin
,
T.
Schenk
,
N.
Godard
,
S.
Girod
,
J. B.
Chemin
,
R.
Leturcq
,
N.
Valle
,
S.
Klein
,
C.
Chappaz
, and
E.
Defay
,
Adv. Funct. Mater.
30
,
2003539
(
2020
).
17.
K. R.
Oldham
,
J. S.
Pulskamp
,
R. G.
Polcawich
, and
M.
Dubey
,
J. Microelectromech. Syst.
17
,
890
(
2008
).
18.
R. A.
Dorey
,
S. B.
Stringfellow
, and
R. W.
Whatmore
,
J. Eur. Ceram. Soc.
22
,
2921
(
2002
).
19.
S.
Gates-Rector
and
T.
Blanton
,
Powder Diffr.
34
(
4
),
352
360
(
2019
).
20.
T.
Bakarič
,
B.
Malič
, and
D.
Kuscer
,
J. Eur. Ceram. Soc.
36
,
4031
(
2016
).
21.
G.
Han
,
J.
Ryu
,
W. H.
Yoon
,
J. J.
Choi
,
B. D.
Hahn
, and
D.-S.
Park
,
J. Am. Ceram. Soc.
94
,
1509
1513
(
2011
).
22.
A.
Mazzalai
,
D.
Balma
,
N.
Chidambaram
,
R.
Matloub
, and
P.
Muralt
,
J. Microelectromech. Syst.
24
,
831
838
(
2015
).

Supplementary Material