High-performance piezoelectric composites via β phase programming

Polymer-ceramic piezoelectric composites, combining high piezoelectricity and mechanical flexibility, have attracted increasing interest in both academia and industry. However, their piezoelectric activity is largely limited by intrinsically low crystallinity and weak spontaneous polarization. Here, we propose a Ti3C2Tx MXene anchoring method to manipulate the intermolecular interactions within the all-trans conformation of a polymer matrix. Employing phase-field simulation and molecular dynamics calculations, we show that OH surface terminations on the Ti3C2Tx nanosheets offer hydrogen bonding with the fluoropolymer matrix, leading to dipole alignment and enhanced net spontaneous polarization of the polymer-ceramic composites. We then translated this interfacial bonding strategy into electrospinning to boost the piezoelectric response of samarium doped Pb (Mg1/3Nb2/3)O3-PbTiO3/polyvinylidene fluoride composite nanofibers by 160% via Ti3C2Tx nanosheets inclusion. With excellent piezoelectric and mechanical attributes, the as-electrospun piezoelectric nanofibers can be easily integrated into the conventional shoe insoles to form a foot sensor network for all-around gait patterns monitoring, walking habits identification and Metatarsalgi prognosis. This work utilizes the interfacial coupling mechanism of intermolecular anchoring as a strategy to develop high-performance piezoelectric composites for wearable electronics.

properties 14 . Accordingly, obtaining the β phase is essential for improving the piezoelectricity of fluoropolymers. To this end, various processing methods were utilized to induce dipole alignment for increasing the electroactive β phase content including mechanical stretching 15 , electric poling 16 and thermal annealing 17 . Meanwhile, the inclusion of TrFE monomers into PVDF chains can also modify the crystallization kinetics of PVDF films and lead to morphotropic phase boundary (MPB)-like behaviors (e.g., the high piezoelectric activity) in piezoelectric polymers 8,11,[18][19][20] . By synergizing the merits of both highperformance piezoelectric ceramics and flexible fluoropolymers, the piezoelectric composites are the choice of materials for nextgeneration on-body bioelectronics toward widespread scenarios, such as soft robotics 21 , biomonitoring [22][23][24][25][26] , and human-machine interface 5,13,27 . Nevertheless, the random distribution and complex connectivity of ceramic fillers within the fluoropolymer matrix cause inhomogeneity and discontinuity of inorganic-organic interfaces, which hampers the long-range alignment of dipole moments in the polymeric chains and thus the all-trans conformation (i.e., polar βphase) [28][29][30] . Meanwhile, the large discrepancy in dielectric permittivity between ceramic fillers and polymer matrix dramatically weakens the applied electric field in the ceramics during electric poling and thereby restricts the domain evolution and permanent polarization [31][32][33] . Consequently, it is highly desired to develop high-performance piezoelectric composites with long-range all-trans conformation and strong polarization.
In order to promote and stabilize the ratio of the polar β-phase, versatile fillers such as carbon nanotubes 34 , silver nanowires (AgNWs) 6,35 , ferrites [36][37][38][39] , graphene oxide (GO) and reduced graphene oxide (rGO) salts [40][41][42] have been introduced as a nucleating agent to maximize and maintain the zigzag alignment of methylene (−CH 2 ) and difluoromethylene (−CF 2 ) groups within the polymer chains. Note that functional groups of the GOs and their derivates are conducive to forming dipole moment reinforcements to sustain the residual orientation of β-phase nanocrystals 40 . Meanwhile, the essential factor for the nucleation of the β-phase in the PVDF nanocomposites lies in the static electric interaction between the fillers with a negative zeta potential and the −CH 2 groups having a positive charge density [35][36][37][38] . Despite the formation of β-phase could be facilitated partially by doping nanofillers, the sparse functional groups, discontinuous interface, low compatibility and inability to align domain on a large-scale of conductive fillers make it rather challenging to attain a homogeneous and long-range molecular interaction with PVDF molecular chains 43,44 , which significantly hinders the dipole polarization and thus the alltrans conformation of the as-synthesized polymer composites. Meanwhile, the mechanism of domain alignment enabled by intermolecular interaction at the inorganic-organic interfaces in piezoelectric composites is still unclear and underexplored.
Herein, we demonstrate a simple and efficient strategy for tailoring the local dipole moment and β phase content of piezoelectric polymer composites by introducing Ti 3 C 2 T x MXene nanosheets. Employing phase-field simulation and molecular dynamics (MD) calculations, we show that OH surface terminations on the Ti 3 C 2 T x nanosheets offer hydrogen bonding with the fluoropolymer matrix, leading to dipole alignment and enhanced net spontaneous polarization of the polymer-ceramic composites. We then translated this interfacial bonding method into electrospinning to rationalize the piezoelectric properties of samarium doped PMN-PT/PVDF (Sm-PMN-PT/PVDF) nanofibers via Ti 3 C 2 T x nanosheets inclusion. It was found that an appropriate loading (2.5 wt%) of Ti 3 C 2 T x flakes below the percolation threshold effectively strengthen the polarization efficiency and interfacial coupling between the inorganic nanofillers and the organic polymer matrix, promoting the piezoelectricity by 160 % in comparison with the undoped version. Toward practical application, a soft piezoelectric textile (PT) sensor was developed by electrospinning the MXene-enabled piezoelectric composite (MPC), which can be integrated into shoe insole for continuous gait patterns monitoring, walking habits identification, and Metatarsalgia prognosis.

MD simulation of Ti 3 C 2 T x enabled intermolecular anchoring effect
Ti 3 C 2 T x nanosheets were exfoliated from the Ti 3 AlC 2 MAX phase by selectively etching Al layers using a HCl/LiF solution, during which the surface transition metal spontaneously reacts with water or fluoride ions to produce hydroxyl (−OH), oxygen (−O), and fluoro (−F) surface terminations. These abundant functional groups terminated on the Ti 3 C 2 T x lamellae arising from chemical etching are favorable for designing hydrogen, ionic and covalent linkage with surrounding dielectrics [45][46][47] . To this end, Ti 3 C 2 T x nanosheets were employed to establish hydrogen bonds and electron-dipole interaction with C−H and C−F moieties of the PVDF molecules, enabling strong intermolecular bindings with PVDF polymer chains as shown in Fig. 1a. The local anchoring of PVDF chains on the two-dimensional MXene platelets directs in-situ alignment and orientation of CH 2 and CF 2 moieties and therefore transition from initial randomly coiled conformations (left part of Fig. 1a) to long-range all-trans conformation (right part of Fig. 1a), enlarging the macroscopic out-of-plane polarization and piezoelectricity.
To understand the interaction between the Ti 3 C 2 T x nanosheets and fluoropolymer, MD calculations were implemented using the periodic lattice of Ti 3 C 2 T x with −OH surface terminations and 60 "mer" chains of PVDF (Fig. 1b). To better illustrate the functionality and mechanism of hydrogen bonds between −OH group and −CF 2 moieties of the PVDF, the corresponding MXene flakes functionalized with −O termination (Fig. 1c) were also constructed and calculated with the same number of "mer" chains. Apparently, after 2000 ps interaction, the PVDF chains anchored on the Ti 3 C 2 (OH) 2 lamellae accomplish an out-of-plane polarization of 371.42 D under electrical poling while the counterpart versions on the Ti 3 C 2 O 2 platelets trigger an out-of-plane polarization of 308.13 D ( Supplementary Fig. 1), suggesting the enhancement of self-assembly of highly aligned PVDF chains on the Ti 3 C 2 (OH) 2 flakes via the intermolecular interaction enabled by hydrogen bonding. Additionally, the pulling force needed to apply on each atom of the single PVDF chain is respectively 2.919 pN and 1.807 pN for Ti 3 C 2 (OH) 2 and Ti 3 C 2 O 2 to initialize the separation (Fig. 1d), implying a stronger interaction and tighter locking between Ti 3 C 2 (OH) 2 flakes and PVDF chains. Without hydroxyl groups to guide the alignment of moieties of the PVDF molecule, the dipoles within the PVDF matrix would present irregular distribution to minimize the internal electrostatic energy (Supplementary Figs. 2,3).
The phase transition could be clearly reflected from the variation of dihedral angles in bond conformation within the PVDF chains, where the most favorable torsional bond arrangement has substituents at 180°(trans or T) to each other instead of the ones at ±60°(gauche or G) 14 . Notably, the −OH functionalized Ti 3 C 2 T x platelets possess a higher proportion of bond conformation near 180° (Fig. 1f) but a lower fraction of bond conformation near −60°( Fig. 1e and Supplementary  Fig. 4) vice versa, indicating the transform from the gauche conformation to trans conformation as a consequence of hydrogen bonding induced anchoring effect. Interestingly, according to the MD calculation, Ti 3 C 2 (OH) 2 platelets give rise to a higher fraction of T bond fraction ( Supplementary Fig. 5) and TTTT (all-trans) conformation (Fig. 1g) in the PVDF chains in comparison with that of Ti 3 C 2 O 2 , which validates the role and functionality of hydrogen bonds in terms of arrangement and orientation of polymer chains toward all-trans conformation (i.e., polar β phase), promoting the net spontaneous polarization and the piezoelectricity of the fluoropolymer composites. Aside from the PVDF, the β-phase content and polarization of poly(vinylidene fluoride-co-trifluoroethylene) (P(VDF-TrFE)) copolymers can also be enhanced by this proposed MXene anchoring method according to the MD simulation (Supplementary Figs. 6, 7), validating the universality and feasibility.

Design of MXene/Sm-PMN-PT/PVDF nonwoven piezoelectric composites
To bring the aforementioned proposed principle and strategy into reality, electrospinning was adopted to construct MXene doped fluoropolymer due to its inherent merit of integrating in-situ stretching with local poling in one step. Figure 2a elucidates the synthesis procedure and composition of the MXene/Sm-PMN-PT/PVDF nonwoven piezoelectric composites, in which a suspension of MXene powder, Sm-PMN-PT particles and PVDF was used as a precursor for subsequent electrospinning (Supplementary Fig. 8). Sm-PMN-PT was selected as the ceramic fillers because of its large piezoelectric coefficient (~1500 pC/N) 2 in comparison with the conventional piezoelectric ceramics such as BTO (~190 pC/N) 48 , PZT (560 pC/N) 49 and PMN-PT (~620 pC/N) 3 (bottom inset of Fig. 2a). Meanwhile, in addition to intermolecular linkage bound with PVDF polymer chain (upper inset of Fig. 2a), the introduction of MXene lamellae efficiently enhances the conductivity and mechanical ductility of the precursor and thus promotes the subsequent local electric poling and in-situ mechanical stretching of polymer composites during electrospinning, which is beneficial to achieving the zigzag alignment of −CH 2 and −CF 2 and thus the all-trans conformation within PVDF polymer (Supplementary Fig. 9). Figure 2b reveals the morphology and element structure of the assynthesized piezoelectric nonwoven textiles consisting of randomly stacked molten PVDF/Sm-PMN-PT/MXene nanofibers through the electrospinning method. The incorporation of MXene nanosheets exhibits a negligible impact on the morphology of the fibrous configuration ( Supplementary Fig. 10). Energy-dispersive spectrometer (EDS) mapping spectra of Pb, F, and Ti elements further confirms the uniform distribution of Sm-PMN-PT, PVDF and MXene in the composite films ( Fig. 2b and Supplementary Fig. 11). Note that the formation of TiO 2 originates from the oxidation of MXene on the surface (Supplementary Figs. 12,13). According to the high-resolution transmission electron microscope images, PVDF polymeric chains aligning along the inorganic-organic interface validate the zigzag alignment of −CH 2 and −CF 2 (Fig. 2c). In addition, X-ray photoelectron spectrometer (XPS) spectra ( Fig. 2d and Supplementary Fig. 14) evidently demonstrates that the oxygen functionalities and hydrogen bonds were anchored on the surface of the as-synthesized nanofibers, where the emerging peaks of C−O and C=O groups resulting from the MXene incorporation suggest that strong interactions such as hydrogen bonds are formed between the −OH and −O x functional groups in the MXene flakes and the −CF 2 groups within PVDF, evidencing the aforementioned hypothesis. Higher binding energy peaks observed in O 1s region ( Supplementary Fig. 15) for C−Ti−O and C−Ti−OH demonstrate that hydroxy termination is dominant on the surface of the flakes. Figure 2e elucidates the X-ray diffraction (XRD) pattern of the assynthesized nanofibers and pristine MXene and Sm-PMN-PT powders. The characteristic 2θ peak of the Ti 3 C 2 T x lamellae is 6.48°corresponding to (002) lattice plane with a d-spacing of 13.66 Å, indicative of successful removal of the Al atom layer from the multilayer Ti 3 AlC 2 MAX phase during the chemical etching process ( Supplementary  Fig. 16) 50 . Apparently, the peak intensity (002) of the composites is lower than that of the pristine MXene as the addition of PVDF molecules enhances the distance between MXene flakes. perovskite structure 51 . In addition, the typical diffraction peak at 20.5°f or the crystal plane of (110) maximizes at the Ti 3 C 2 T x mass fraction of 2.5 wt%, indicating a larger ratio of local all-trans conformation (polar β phase) when compared with other versions.
The frequency dependencies of dielectric relative permittivity and dielectric loss of as-prepared MPC textiles and undoped textiles were respectively presented in Fig. 2f, g as a function of MXene inclusion amount at room temperature. The dielectric permittivity for all the samples slightly declines with increasing frequency as a consequence of relaxation (Fig. 2f). It is found that the dielectric permittivity pronouncedly undulates with MXene doping amount and reaches a maximum at 2.5 wt% over a frequency range of 1 kHz-1 MHz, while the dielectric loss maximizes at 5 wt%. This is because the proper inclusion of conductive MXene nanosheets is favorable for interfacial polarization between fillers and polymer matrix, whereas the excessive addition of MXene into polymer composites gives rise to an augmented leaking current and thereby a larger dielectric loss.

Characterization and phase-field simulation
To further analyze the interfacial coupling effect between the fluoropolymer matrix and the inclusion, piezoelectric force microscopy (PFM) is employed to assess the ferroelectric and piezoelectric properties of as-electrospun nanofibers with or without MXene doping. The PFM amplitude, PFM phase, phase-field simulated domain structure, and MD calculated polarization were presented for a single undoped Sm-PMN-PT/PVDF fiber ( Fig. 3a-d) and a MPC fiber (Fig. 3e-h). The PFM results show a larger area of single ferroelectric domain observed in the MPC fiber in comparison with that in the undoped fiber, implying an enhanced piezoelectric response after MXene addition. This was both certified by the phase-field simulation ( Fig. 3c, g) of the ferroelectric domain structure and MD calculation results (Fig. 3d, h) that the out-of-plane domain size of polymer composite becomes larger in the vicinity of the ceramic particle by adding MXene flakes. The phase-field simulation reveals a multi-domain state with 180 degrees domains in the polymer matrix for both undoped and MPC fibers, where the MPC fiber possesses a larger domain fraction with polarization aligned along the poling direction (PD) of electrospinning compared with that of the undoped fiber. Moreover, the embedded Sm-PMN-PT particle maintains a polycrystalline structure as its coercive field is much higher than the in-situ poling field (1.2 kV/cm) of electrospinning 2 .
The modulated dipole structure of MPC fibers with increased polarization further gives rise to an increased piezoelectric response. PFM measurements show that all the as-synthesized electrospun fibers exhibit a butterfly-shape hysteresis of the amplitude versus bias, where the MPC fiber displays a higher amplitude (Fig. 3j) than the undoped one (Fig. 3i), implying a larger piezoelectric coefficient (d 33 ). As displayed in Fig. 3k, the d 33 derived from PFM test first increases and then decreases upon increasing the MXene doping amount, with a maximum d 33 achieved at 2.5 wt% of MXene. It is worth noting that both phase-field simulation and MD calculation follow a similar changing tendency of polarization as well as out-of-plane domain proportion to the PFM results (Fig. 3l), which verifies the accuracy and reliability of our proposed theoretical modeling.
To achieve a more intuitive understanding of the functionality and mechanism of MXene lamellae addition, a schematic was utilized to elaborate the β crystalline phase conformation and orientation during the electrospinning of the MXene/Sm-PMN-PT/PVDF composites. For bare PVDF nanofibers (Fig. 4a), β phase crystalline can be partially aligned by the in-situ poling while some localized amorphous microstructures still exist. The situation was changed once a small amount of MXene (1.25 wt% and 2.5 wt% in our case) is added into the PVDF-based composite, as shown in Fig. 4b. In addition to the hydrogen bonds formed between PVDF chains and MXene sheets, the applied electric field yields inductive charges on the surface of conductive MXene nanosheets and thereby provokes a greater Coulomb force during the electrospinning. This electrostatic force further attracts PVDF chains to crystallize on the MXenes surface in the zigzag arrangement, which expedites the transformation of local amorphous regions into the β crystalline phase. Consequently, the β phase is enhanced when compared with pure PVDF nanofibers (Fig. 4a). However, superfluous addition leads to a stack of MXene nanosheets and suppresses the orientation and alignment of polymeric chains in the way of forming β phase crystalline. In the meantime, the connection of these stacked MXene flakes boosts the leaking current, which hinders the effective dielectric permittivity (ε 33 ) and piezoelectric coefficient (d 33 ) of the composite. To verify this proposed modeling on the crystallinity and phase shift, the Fourier transform infrared (FTIR) spectrum of the undoped textiles and MPC textiles doped with various MXene mass fractions were revealed in Fig. 4d. The adsorption peaks at 763 cm −1 and 973 cm −1 are attributed to α-phase crystalline structure of PVDF and the adsorption peak at 841 cm −1 , 1276 cm −1 and 1429 cm −1 are ascribed to a βphase crystalline structure 52 . According to Beer-Lambert equation 5 , the fraction of the β phase (F β ) of as-synthesized nanofibers can be attained as a function of MXene loading amount ( Supplementary  Fig. 17), where the highest β phase of 83.4% occurs at the loading of 2.5 wt% MXene. Consequently, the phase evolution trend upon MXene doping revealed by the FTIR test agrees well with the proposed principle.
To further testify the above hypothesis and shine some light on the mechanism, phase-field simulation was carried out to numerically illustrate the impact of conductive fillers on the domain structure and polarization of the ferroelectric composites. As shown in Fig. 4e-g, incorporation of a proper content (2.5 wt%) of MXene lamellae notably improves the polarization of the as-electrospun composite, whereas further accretion jeopardizes the polarization and may even cause it to fall below the undoped counterpart (Fig. 4h), which is in good accordance with foregoing deduction and FTIR test. Furthermore, as for the simulated domain structure (Fig. 4i-k), the incorporation of MXene drives the polarization toward the out-of-plane direction parallel to the electrospinning poling field when compared with the undoped one, which is also consistent with the PFM results. A loading of 2.5 wt% MXene contributes to an enhancement of the volume fraction of PD out-of-plane domain from 67.8 to 72.6% (Fig. 4l and Supplementary Figs. 18, 19). It is believed that the net out-of-plane domain is responsible for the augmented β phase content and corresponding larger piezoelectric response. Given that the piezoelectric coefficient is proportional to spontaneous polarization and permittivity (d 33 = 2P s ε 33 Q 33 ), a proper doping of the conductive component is proven to be an efficient approach to upgrade the piezoelectric performance of polymer composites.
To understand the role of MXene doping on the percolation effect in the as-electrospun composites nanofibers, a phase-field simulation of effective properties was implemented to analyze the dependence of the MXene doping amount on the dielectric, piezoelectric and mechanical properties of the electrospun unmodified and MPC textiles. As displayed in Fig. 5a, a three-dimensional fibrous configuration was randomly computer-generated to imitate the structure of electrospun nonwoven textiles ( Supplementary Fig. 10). The stress distribution, electric field and electric potential of MPC textiles under constant applied stress of 1 × 10 6 Pa along the out-of-plane (z) direction are calculated. As displayed in Fig. 5b, the calculated piezoelectric coefficient d 33 and dielectric permittivity ε 33 undulate with increasing MXene mass fraction from 0 to 5.0 wt%, reaching a maximum value at 2.5 wt%. A similar trend is observed on both the d 33 meter measured d 33 value and the calculated strain-induced piezoelectric polarization, two of which fluctuate with the loading of MXene flakes and reaches a maximum at 2.5 wt% (Fig. 5c and Supplementary Fig. 20). On the other hand, the stiffness coefficient retains almost the same at low MXene doping amount but remarkably decreases at high MXene doping level, indicating a hardening of the composites with high contents of doped MXene flakes.

MPC textiles for energy harvesting and sensing
To explore the device performance of the fabricated MPC textiles in response to external stimulation, the output piezoelectric signals were systematically characterized under a diversity of applied forces varying from 1 to 9 N, as plotted in Fig. 5d, e. Notably, the electric outputs of the MPC textiles are proportional to applied stress and the MXene doped devices deliver larger voltage and current than the unmodified version. Furthermore, under a fixed external force of 3 N, the signal intensity varies with the changing doping amount of MXene and the 2.5 wt% MXene doped nanocomposite yields a gain of 160% when compared with the undoped sample ( Supplementary Figs. 21,   22). It is worth pointing out that the electric outputs follow an analogous changing tendency as the calculated piezoelectric coefficient upon conductive filler doping (Fig. 5e), which convincingly confirmed the validation and feasibility of our proposed theoretical modeling and corresponding calculation. Consequently, all the following measurements were carried out using an optimal MXene mass fraction of 2.5 vol% because of its favorable piezoelectric properties. Figure 5f shows the impedance response of the prepared MPC textile doped with 2.5 wt% MXene under various resistances from 10 KΩ to 100 GΩ, where the maximum output power was achieved at 80 MΩ. The mechanical stretchability and robustness of as-synthesized MPC textiles were testified by the long-term fatigue test as illustrated in Fig. 5g, where the output voltage retains a constant value when subjected to more than 7000 cycles of 5 N impact force application, demonstrating an excellent mechanical robustness and reliability for practical scenarios.
To verify the sensing capability of MPC textile for active biomonitoring, a self-powered gait monitoring system was developed. As schemed in Fig. 6a, by attaching as-synthesized devices onto five different positions (M1-M5) of a conventional insole (see Supplementary Note 1, Supplementary Figs. 23, 24 for details), the stress distribution of the foot can be spontaneously detected. It can be clearly seen that the incorporation of MXene triggers a much larger amplitude of output charge than the undoped one when the big toe pressing the device on site of M5 (Fig. 6b) regardless of walking postures including pigeontoed, normal and splayfooted (Fig. 6c). Furthermore, various gait patterns of walking, running, jumping, falling forward, and falling backward can be accurately recognized and distinguished with regards to the dynamic signal mapping among the five units ( Fig. 6d and Supplementary Movie 1). It is interesting that the falling forward and falling backward possess adverse signal mapping profiles due to the  Fig. 25). Note that even though the as-electrospun PTs might deliver smaller output signals than the triboelectric devices 53 , they do not need an extra mobile component like triboelectric devices to trigger electricity generation, allowing for allin-one configuration and biocompatibility 12 . To demonstrate the potential of MPC textiles in identifying walking habits, the posture abnormity like pigeon-toed or splayfooted was further distinguished using the as-prepared smart insole ( Fig. 6e and Supplementary Fig. 26). Apparently, the splayfoot renders a much higher pressure on the lateral side but lower pressure on the inner side of insole in comparison with the normal posture, which gives rise to a larger signal amplitude on position M1 and M3. In contrast, the pigeon-toed walking posture triggers a much higher pressure on the inner side but lower pressure on the lateral side of the insole when compared with the normal posture, leading to a larger amplitude on position M2, M4 and M5. Aside from the gait monitoring, the prepared smart insole based on MPC textiles also reveals capability in clinic diagnosis. As a widespread foot condition, Metatarsalgia brings about a shape deformation of the metatarsal region and subsequent sharp pain in the ball of the foot. Figure 6f unveils and compares the walking signal profiles of a healthy person and a Metatarsalgia patient on the metatarsal region (M3, M4, M5). Due to the classical symptom of the curled toe, a signal contour with larger intensity on M5 and smaller intensity on M3 was recognized for the Metatarsalgia patient when compared to that of the healthy tester ( Supplementary Fig. 27), which demonstrates the capacity of MPC textiles in clinic prognosis of Metatarsalgia.

Discussion
A high-performance piezoelectric composite was fabricated via termination engineering using Ti 3  moieties in PVDF chains is responsible for the long-range growth of TTTT configuration and alignment of domain orientation, which efficiently boosts the β phase content and piezoelectric coefficient. We then translated the hydrogen bond anchoring strategy into reality by synthesizing the Ti 3 C 2 T x /Sm-PMN-PT/PVDF composites-based soft PTs via electrospinning. It is found that the incorporation of 2.5 wt% Ti 3 C 2 T x into Sm-PMN-PT/PVDF composites leads to an output gain of 160% when compared with the undoped version. Furthermore, the as-electrospun MPC could be easily integrated into the conventional insole for active gait patterns monitoring, walking habits identification and Metatarsalgi prognosis. This work theoretically and experimentally looks into the underlying mechanism of the interfacial coupling effect of piezoelectric nanocomposites, opening up a paradigm for the development of high-performance wearable electronics.

Synthesis of Sm-PMN-PT piezoelectric ceramic particles
The composition of the Sm-PMN-PT is Pb 0.97 Sm 0.02 [(Mg 1/3 Nb 2/3 ) 0.7 Ti 0.3 ]O 3 3 . MgNb 2 O 6 powder was initially prepared using a solid-state reaction at 1200°C. Pb 3 O 4 , MgNb 2 O 6 , TiO 2 , and Sm 2 O 3 powders were mixed by ball-milling overnight as wet dispersion. Calcination of the mixed powders was conducted in a non-contaminated furnace at 850°C for 2 h.

Synthesis of layered MXene nanosheets
Firstly, 10 ml ultrapure water was utilized to dilute 30 ml 12 M hydrochloric acid solution (Chron chemicals Co., Ltd, China) to obtain 40 ml 9 M HCl. Subsequently, 2 g LiF (Aladdin, China) was added into 40 ml 9M HCl followed by magnetic stirring for 5 min in order to obtain a hydrofluoric acid solution. Next, 1.2 g Ti 3 AlC 2 powder (MAX phase, purchased from 11 Technology Co., Ltd., China) was weighed and slowly added to the mixed solution for water bath stirring under a rotation speed of 400 rpm for 24 h at 35°C. Then, the reaction product was centrifugated at 1737 × g for 5 min several times and rinsed with deionized (DI) water until the pH value reached 6. After suction filtration and vacuum drying (45°C, 24 h) followed by an ice bath ultrasonic, the MXene nanosheets were attained for spare.

Fabrication of MPC textile-based sensor
First, the as-synthesized MXene powders of 0, 29.5, 59, 88.5, 118 mg, corresponding to the mass fraction of 0, 1.25, 2.5, 3.75, 5 wt% respectively, were dispersed in a beaker of 6 mL dimethylformamide (DMF) and 4 mL acetone mixed solvent and ultrasonically treated for 20 min to evenly disperse nanosheets in DMF solvent. Then, 118 mg Sm-PMN-PT nanoparticles together with 2.25 g PVDF powder were added into the mixed solution via magnetic stirring at 50°C for 2.5 h in a water bath to obtain a stable and uniform electrospinning precursor solution. After solution agitation, ultrasonic treatment for 30 min allows nanoparticles to disperse evenly. Subsequently, the as-received solution was added into a 10 mL plastic syringe with a needle orifice of 25 gauge for electrospinning. The electrospinning was carried out for 2 h with a gap distance of 15 cm, a pushing speed of 0.6 mL/h and an electric field of 1.2 KV/cm at room temperature and 42% humidity (Zhiyan Technology Co., Sichuan). Then, as-electrospun textile was adhered by two pieces of aluminum tapes as electrodes and tailored into small pieces with a size of 3.0 cm × 3.0 cm. Finally, the whole device was packaged using PET substrates and transparent medical tape to build up a MPC textile-based sensor.

Characterization and measurement
The d 33 values were characterized using a ZJ quasi-static curvature radius of 35 nm and a tip height of 14-16 µm, where an alternating-current voltage is applied to the metallic disk (=150 kHz, U RMS = 0.5 V) and two Stanford Research 830R lock-in amplifiers are used to monitor the amplitude and phase of the deflection and torsion of the cantilever. Conductive Platinum-Iridium silicon cantilevers (NSG01/Pt) were used for the PFM characterization with detailed parameters listed in Supplementary Table 1. The phase purity and crystal properties of the samples were analyzed by an X-ray diffractometer (XRD: D8 Advance, Bruker-AXS, Germany) equipped with Cu-Kα radiation (λ = 1.5418 Å). An XPS (ESCALAB 250Xi) was used with an Al Kα excitation (1486.8 eV). The chemical composition and valence states were characterized by X-ray photoelectron spectroscopy (XPS, model ESCALAB 250, Al Kα, hυ = 1486.6 eV). The FTIR spectra were performed on Fourier infrared spectrometer (Smartlab 9, Rigaku, Japan). A linear motor with tunable frequency, amplitude, velocity and force was adopted to apply a diversity of stress and strain on the samples.

Phase-field simulation of the ferroelectric domain structures
In the phase-field simulation, the time-dependent Ginzburg-Landau equation was employed to describe the evolution of the polarization 54 [Eq. 1], where L is the kinetic coefficient, t is the time, F is the total free energy of the system, r = (x, y, z) is the spatial position vector, and Pi (r, t) is the polarization field. δF ∕ δ Pi (r, t) is the thermodynamic driving force for the spatial and temporal evolution of polarization. The total free energy of the system includes contributions from the bulk free energy f bulk , the elastic energy f elastic , the electrostatic energy f electric , and the gradient energy f grad : [Eq. 2] In particular, the bulk free energy density f bulk is expressed following the Landau theory. The bulk free energy density of Sm-PMN-PT is written as a 6-th order polynomial 55 f bulk = α 1 P 2 x + P 2 y + P 2 z + α 11 P 4 x + P 4 y + P 4 z + α 12 P 2 x P 2 y + P 2 x P 2 z + P 2 y P 2 z + α 111 P 6 x + P 6 y + P 6 z + α 112 P 2 x ðP 4 y + P 4 z Þ + P 2 y ðP 4 z + P 4 x Þ + P 2 z ðP 4 x + P 4 y Þ The expression is established following a uniaxial anisotropy with spontaneous polarization along the z axis (which is the PD during electrospinning), while polarization components along x and y axes assume paraelectric properties with higher order terms neglected.
Incorporation of MXene was considered through treating different background dielectric constants ε b by solving the electrostatic equilibrium equation 56 , where ε b correlated to MXene doping amount (Fig. 2f) to account for the variation of β phase fraction [Eq. 5]. The evolution of the domain structure was simulated for a single electrospun fiber under an applied electric field of 1.2 × 10 5 V/m, following the in-situ poling field (1.2 KV/cm) of electrospinning. The simulations were performed in a three-dimensional system with a total size of 512 nm × 512 nm × 512 nm containing a single Sm-PMN-PT particle (diameter = 100 nm) in a PVDF matrix, which was then discretized into a three-dimensional array of 128 × 128 × 128 grids with a grid size of Δx = Δy = Δz = 4 nm. Periodic boundary conditions were employed for the polarization field, the electric field, and the mechanical displacement field. The material constants of Sm-PMN-PT and PVDF are listed in Supplementary Tables 2 and 3, respectively. All phase-field simulations in the present work were performed using the MuPRO Ferroelectric module for linux cluster software released from Mu-PRO LLC.

Effective properties calculation of the piezoelectric composites
These simulations were performed using Mu-PRO Effective Properties module for linux cluster released from Mu-PRO LLC to systematically study and compare the stress transfer capability and electromechanical coupling effect of the electrospun nanofiber architecture with and without the MXene incorporation. The composite system consists of randomly oriented Sm-PMN-PT filler fibers with a diameter d = 200 nm at a total content of 3.23 vol% in a PVDF matrix. The total system size was taken as 12.8 µm × 12.8 µm × 12.8 µm, which was discretized into a three-dimensional array of 128 × 128 × 128 grids. A threedimensional periodic boundary condition was employed.
The effective dielectric constant ε ef f r and piezoelectric coefficient d ef f r of the composite system was calculated by modeling the coupled response of the electric displacement field D(r) and the strain ε(r) of the composites to a small applied testing electric field E app under a stress-free condition. The material response was simulated by solving the electrostatic and elastic equilibrium equations coupled with linear constitutive relations 5 , i.e., where ε 0 is the vacuum permittivity, ε r (r) is the relative dielectric constant of the local phase, c(r) is the elastic stiffness, and d(r) is the piezoelectric coefficient [Eqs. 6,7]. A Fourier-spectral iterativeperturbation method 57 is employed for obtaining the numerical solution. The effective material properties of the composites were then calculated following ε 0 ε ef f r E app = ΔD h iand d ef f r T E app = Δε h i, where ΔD h iand Δε h iare the spatial average of the displacement response and the strain response with the whole system, respectively. The materials constants of the Sm-PMN-PT fillers and PVDF matrix are listed in Supplementary Note 2.

Molecular dynamics (MD) calculation
In the MD simulations, 60 PVDF molecules with each 30 VDF monomers were constructed at 1 nm from the Ti 3 C 2 T x nanosheets having an initial film density of 1.3 g/cm 3 . The periodic lattice of Ti 3 C 2 T x nanosheets with OH or O surface terminations was directly taken from the literature 58 and was fixed at their lattice positions, interacting with the PVDF chains only via van der Waals and electrostatic interactions. The Universal Force Field (UFF) 59 was utilized to describe the intra-and intermolecular interactions for PVDF chains as well as the interaction between the PVDF and Ti 3 C 2 T x nanosheets. The partial charges of the atoms of PVDF chains were assigned using the bond increments method, whereas the charges for the atoms of the Ti 3 C 2 T x substrate were adopted from the first-principle calculations 60 . For equilibration, the PVDF-Ti 3 C 2 T x composite model was initially energy minimized with the energy and force convergences of 1 × 10 −4 kcal mol −1 and 1 × 10 −6 kcal mol −1 Å −1 , and then followed by a 2.5 ns run in the NVT ensemble. After that, an electrical field of 1.0 V/nm was applied to the PVDF chains normal and opposite to the Ti 3 C 2 T x substrate (i.e., z direction) for another 2.0 ns, to investigate the polarization performance of the PVDF. The induced polarization was then calculated by summing the dipole moments of all the atoms of the PVDF film within the film volume V as [Eq. 8] where q i and r i are the particle charge and Cartesian coordinate vector of the ith atom of the PVDF. In addition, the interaction strength between one single PVDF chain and the Ti 3 C 2 T x substrate was estimated by applying a constant force along the z direction to each atom of the PVDF chain. A series of MD simulations were performed with varying forces using an incremental step of 0.01 or 0.002 kcal mol −1 Å −1 to produce the minimum pulling force required for fully desorbing the PVDF chain away from the Ti 3 C 2 T x substrate. For all the simulations, the temperature was maintained at 300 K using the Nośe-Hoover thermostat with a relaxation time of 1.0 ps. Newton's equations of motion were integrated using the velocity-Verlet algorithm with a time step of 1.0 fs and periodic boundary conditions were applied in all three dimensions. The non-bonded van der Waals interactions were cut off at 10 Å and the long-range electrostatic interactions were calculated using the particle-particle particle-mesh solver 61 with 5 × 10 −6 accuracy tolerance. All the MD simulations were carried out using the Large-scale Atomic/Molecular Massively Parallel Simulator software 62 .

Data availability
The data that supports the findings of this study are available from the authors on reasonable request.
Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons license, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons license and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this license, visit http://creativecommons.org/ licenses/by/4.0/.