Electromagnetic field controlled domain wall displacement for induced strain tailoring in BaTiO3-epoxy nanocomposite

Failure in an epoxy polymer composite material is prone to initiate by the coalescence of microcracks in its polymer matrix. As such, matrix toughening via addition of a second phase as rigid or/and rubber nano/micro-particles is one of the most popular approaches to improve the fracture toughness across multiple scales in a polymer composite, which dissipates fracture energy via deformation mechanisms and microcracks arrest. Few studies have focused on tailorable and variable toughening, so-called ‘active toughening’, mainly suggesting thermally induced strains which offer slow and irreversible toughening due to polymer’s poor thermal conductivity. The research presented in the current article has developed an instantaneous, reversible extrinsic strain field via remote electromagnetic radiation. Quantification of the extrinsic strain evolving in the composite with the microwave energy has been conducted using in-situ real-time fibre optic sensing. A theoretical constitutive equation correlating the exposure energy to micro-strains has been developed, with its solution validating the experimental data and describing their underlying physics. The research has utilised functionalised dielectric ferroelectric nanomaterials, barium titanate (BaTiO3), as a second phase dispersed in an epoxy matrix, able to introduce microscopic electro-strains to their surrounding rigid epoxy subjected to an external electric field (microwaves, herein), as result of their domain walls dipole displacements. Epoxy Araldite LY1564, a diglycidyl ether of bisphenol A associated with the curing agent Aradur 3487 were embedded with the BaTiO3 nanoparticles. The silane coupling agent for the nanoparticles’ surface functionalisation was 3-glycidoxypropyl trimethoxysilane (3-GPS). Hydrogen peroxide (H2O2, 30%) and acetic acid (C2H4O2, 99.9%) used as functionalisation aids, and the ethanol (C2H6O, 99.9%) used for BaTiO3 dispersion. Firstly, the crystal microstructure of the functionalised nanoparticles and the thermal and dielectric properties of the achieved epoxy composite materials have been characterised. It has been observed that the addition of the dielectric nanoparticles has a slight impact on the curing extent of the epoxy. Secondly, the surface-bonded fibre Bragg grating (FBG) sensors have been employed to investigate the real-time variation of strain and temperature in the epoxy composites exposed to microwaves at 2.45 GHz and at different exposure energy. The strains developed due to the in-situ exposure at composite, adhesive and their holding fixture material were evaluated using the FBG. The domain wall induced extrinsic strains were distinguished from the thermally induced strains, and found that the increasing exposure energy has an instantaneously increasing effect on the development of such strains. Post-exposure Raman spectra showed no residual field in the composite indicating no remnant strain field examined under microwave powers < 1000 W, thus suggesting a reversible strain introduction mechanism, i.e. the composite retaining its nominal properties post exposure. The dielectric composite development and quantifications presented in this article proposes a novel active toughening technology for high-performance composite applications in numerous sectors.


Research hypothesis
Dielectric nanomaterials exhibit electric field induced strain which is attributed to intrinsic mechanisms from lattice deformation and extrinsic mechanisms due to domain wall (DW) movement 31 , extensively used as actuators and transducers. The inclusion of such material within a rigid epoxy FRP materials can impose a compressive stress field in its surrounding epoxy matrix when activated its DW movements by external electric field stimulation. As result of the DW movements, a microwave stimulation at GHz frequencies induces effective dipolar displacement (leading to intrinsic strains) to the nanomaterial's molecules that, at the interface with their surrounding rigid polymer, is converted to compressive mechanical strain. The hypothesis of this research was based upon suggesting that microcrack propagation during dynamic and impact events would be suppressed under such microwave induce d compressive stress-strain field, i.e. higher strain energy would be required to create new fracture surfaces, however the current article presents attempts on the quantification of the field induced strains. Therefore, the study focused on: 1-development of a modified epoxy with uniformly distributed dielectric nanomaterials that exhibit electro-strain, 2-in-situ quantification of its strain response under an electromagnetic field (e.g. microwave), and 3-establishment of a theoretical constitutive equation underpinning the correlation between the induced mechanical field and the microwave field.
To investigate the microstructural response of such dielectric nanomaterial embedded nanocomposite under an electromagnetic field, firstly, we developed a modified epoxy system with uniformly distributed barium titanate (BaTiO 3 ) nanoparticles. Materials' properties, functionalisation and fabrication methods developed in the research are described in "Materials and fabrication" section. Constitutive equations underpinning the multi-physics field induced strain and temperature have been provided in "Theoretical multi-physics constitutive equations" section. With a second phase added to the epoxy system, the dielectric properties have been improved. Meanwhile, it also has effects on the curing extent of the epoxy system. Basic dielectric and cure state characterisation methods and their results are presented and described in "General material characterisations" section and "General material characterisations" section, respectively. To experimentally measure the perceivable strain effect induced by the stimulation of embedded BaTiO 3 particles' strains under microwave stimulation, www.nature.com/scientificreports/

Interaction of ferroelectric materials with the microwave field
The mechanism of the material and microwave interaction is particularly complex and it involves the orientation of dipoles, free-to-move electrons, domain wall movement and electron spins, which are activated by the electric and/or magnetic field components 43 . The essence of the interaction between the dielectric material and the electromagnetic field is energy transformation at the molecular level. Therefore, the dielectric properties of a material, which is measured as the complex permittivity quantifies the ability of a dielectric material to absorb and store electrical potential energy when exposed to a microwave field 44 . The complex permittivity that represents the microwave-material interactions are expressed by: where ε ′ is the real parts of the permittivity that indicating the ability to store microwave energy of the material while ε ′′ is the imaginary part of the permittivity that indicating the ability to dissipate the microwave energy 44 . The dielectric loss tangent is expressed by: The loss tangent is the measure of the material to convert absorbed energy into other forms of energy such as heat 45 . The dipoles in an alternating electric field will reorient themselves to align with the new field direction. However, this process cannot occur instantaneously due to inherent non-zero inertia, and thus some time is required to fulfil the reorientation of the dipoles. At microwave frequencies, dipole polarisation is assumed to be the most dominant one for energy transfer at the molecular level.

Theoretical multi-physics constitutive equations
In this section, the theoretical analysis of the strain field induced in BaTiO 3 nanoparticles embedded epoxy under microwave exposure is established. Firstly, the theory of dynamic elastic response of tetragonal BaTiO 3 nanoparticles under alternating electric field at the microwave exposure is presented and discussed 46 . The displacement of the DWs is utilised as the eigenstrain for the estimation of the strain field in the BaTiO 3 -epoxy nanocomposite. A novel approach based on the principle of virtual work of the equivalent eigenstrain is incorporated in the calculation of the inhomogeneous inclusion 47 . Based on this approach, a theoretical framework of microwave activated dynamic elastic response of the nanocomposite has been established. Finally, the non-uniform distribution of the microwave field and the interaction with the samples have been incorporated in the framework.
Deformation in a representative crystallite of BaTiO 3 under the microwave field. The 90° DW movement is considered to have the major contribution to the extrinsic field-induced strain. Therefore, the tetragonal BaTiO 3 nanoparticle is assumed to be formed by several multi-domain crystallites that contain a laminar 90° domain structure. A conceptual schematic of a crystallite is shown in Fig. 2. l is the displacement www.nature.com/scientificreports/ of domain walls under microwave radiation; d is domain size and g (identified in the figure) is the particulate size, i.e. in our case and grain size in polycrystal. The domain wall movement in such crystallite is assumed to be not affected by that in other crystallites. Before the application of the microwave field, it is assumed that the DWs are in an equilibrium state, and all the domains are of an identical size, dielectric properties and elastic properties. Therefore, the interactions between adjacent domains are not included as it is assumed negligible compared with other factors as suggested in a previous study 49 . Moreover, the internal electric and mechanical fields are assumed to be negligible, hence no forces acting on the DWs prior to the microwave exposure 50 . When the microwave field is applied, the interactions with the magnetic component are negligible since BaTiO 3 is non-ferrite. The alternating electric component of the microwave field is assumed to be homogenous across the whole nanoparticle in this research. It is assumed that the spontaneous polarisations in the 90° domains re-orient themselves with the direction of an applied alternating electric field. Consequently, a driving force f A (t) as the function of time is introduced on the DWs that enables the DW motion due to polarisation change 51 . The displacement of the DW movement is denoted as l in Fig. 2. Due to the BaTiO 3 ferroelectricity, such movement introduces an internal electric field E in the nanoparticle due to the change in the polarisation direction. As a result, the internal stress field is also no longer negligible due to the newly developed internal electric field and DW movement from its equilibrium position.
An additional force f I on the DWs is developed subsequently due to the change in the internal electric and stress fields from zero state. Due to the displacements in DWs, the electric dipole moment p induced by the internal electric field, and the elastic dipole moments v induced by the internal mechanical stress field are given by 52 : where P 0 is the spontaneous polarisation, and S 0 is the spontaneous deformation. The correlation between l and the dipole moments is assumed to be linear, and their internal energy W electric under the internal electric field E , and W elastic under mechanical stress field T 13 are: where E is the internal electric field, and T 13 is the internal mechanical stress T ij ( i , j = 1, 2, 3) as illustrated in a cartesian coordinate system ( The expression for the DW motion is 51 : By the derivation of Eqs. (5) and (6), the force f I due to the internal fields per unit area k l is obtained from where k is the force constant per unit DW area, and A DW is the area of the DW. The displacements for the DWs from their equilibrium state under the microwave field are assumed to have the same magnitude, |�l(t)| . The force f I per unit DW area at a random point r , and time t can then be expressed by Eq. (8): The DWs motion is then described by the Newton's second law of motion: where m is the mass of a DW per unit area, f I (t) is the mean value of the DWs driven force f I from the internal fields, f A (t) is the aforementioned exerted force on the DWs by the applied microwave field. The f I (t) is then calculated by averaging the f I (r, t) over the summation of the DWs surfaces ( d ) in a nanoparticle, as presented below 46 .
The DW displacements induced internal electric field E(r, t) and internal mechanical field T 13 (r, t) in Eq. (11) must firstly be calculated prior to the f I (t) calculation. The internal electric field was evaluated in a quasi-static approximation for a given microwave field, and E(r, t) is considered to be proportional to the DW displacement �l(t) . The interal mechanical stress field T ij was considered as a dynamic problem as the DW motions lag behind the 2.45 GHz microwave frequency 46 . In this study, only the extrinsic effect, i.e. the DW motions, is considered as they are commonly accepted to be the primary contributor to the electro-strain compared with the intrinsic effect such as piezoelectric effect and electrostriction.
To obtain the value of the internal mechanical stress field, the change of the spontaneous deformation S 0 introduced by the DW displacements in Eq. (11) was calculated, thus the �S 0 13 (r, t) and �S 0 31 (r, t) were approximated by 46 : www.nature.com/scientificreports/ where w is the thickness of a 90° DW, V DW is the volume of DWs without an applied microwave field, and δ(V DW ) is the Dirac delta function for V DW . The thickness w is approx. 100 Å, much smaller compared with the domain width d and the particle size 48 , therefore, Eq. (12) was simplified to Eq. (13) 46 .
where represents the midplanes of DWs when l = 0 , δ 3 (�) is the Dirac delta function of the with a unit normal directed along the x 3 axis 53 .
The electric component of the oscillating microwave field is a sinusoidal function of time and the DW displacement �l(t) is written as where l m is the maximum displacement of the DWs under the microwave field. The nanoparticle is assumed to be a homogeneous and isotropic dielectric medium, and the internal electric field E(r, t) is proportional to �l(t), where ε 0 is the permittivity of the free space ( ε 0 = 8.8542 × 10 −12 F/m), d is the width of the DW, and h e (r) are approximated as 46 : where ǫ 11 is first component of the dielectric permittivity tensor ǫ ij ( i, j = 1, 2, 3 ) that represents the dielectric permittivity of the crystallite as illustrated in Fig. 2, and ε * is the dielectric permittivity of the nanoparticle. The force f I was then calculated with the internal mechanical stress field T 13 (r, t) , and the internal electric field E(r, t) from Eqs. (15) and (16). The field-induced force f A is due to the effective electric field in the x 1 direction of the applied microwave field as illustrated in Fig. 2. The electric field component E A1 is a function of time as follows: where E m is the amplitude of the electric component of the microwave field, and ∅ is the phase lag of the DW motion relative to the oscillatory microwave field. Equation (17)  where h s , h c , and h e are the average of h s (r), h c (r)and h e (r) over DW surface , evaluated using the expressions outlined in 46 . From Eq. (19), it could be seen that the DW displacement l m is correlated with the properties and strength of the electric component of the microwave field E m and the dielectric permittivity of the nanoparticle. The value of the DW displacement �l m (ω) could be obtained by numerical computations of the h s (ω) and h c (ω).
To get an approximate estimate, a further simplification was made that the nanoparticles are sufficiently small to be considered as a multi-domain single crystal, namely a crystallite. Therefore, the crystallite as illustrated in Fig. 2 could be considered as a single BaTiO 3 nanoparticle without any adjacent crystallite, and hence no corresponding internal electric and stress field. Consequently, Eq. (10) is simplified to: It is assumed that the DW movement under the oscillatory microwave field is a simple harmonic motion with a driven force of f A (t) . Neglecting the lag in the dipoles reorientation, the DW motion becomes in phase with the oscillating field. The driving force at the maximum electric field strength f m is then obtained by: and the amplitude l m : where Y is the young's modulus of the BaTiO 3 nanoparticle and d is the domain width where the area of the domain is approximated to be d 2 . where e ij is the elastic strain, and ε * ij is the eigenstrain. The compatibility condition of total strain ε ij is expressed in terms of displacement u i , i.e.
The stress σ ij is then expressed by the Hooke's Law, where C ijkl is the effective stiffness tensor, C ijkl = C ijlk and C ijkl u k,l = C ijkl u l,k . In the region where ε * kl = 0 , Eq. (25) is re-arranged to: The inversed Eq. (25) is then expressed by: where S ijkl is the compliance tensor that has the same symmetries as the stiffness tensor. The equation of equilibrium is then expressed as follows: where f i is the body force, and the boundary condition is: where n j is the outward unit normal vector to the boundary. The corresponding stress and strain field of a homogeneous inclusion in an infinite medium could be obtained using Green's function method with the profile of the inclusion 56 : The Eshelby inclusion's problem is applied as illustrated in Fig. 3. It is assumed that the inclusion has a dynamic deformation of the eigenstrain ε(t) as a function of time under the microwave field. A fictitious load T is applied to the inclusion to maintain the original shape. The inclusion is then placed back into the hole within the infinite epoxy matrix with the original shape and size. After removing the applied T , the inclusion exerts a traction force F on the matrix, which equals −T.
It is assumed that the inclusion is in ellipsoids shape and have identical elastic properties with the matrix. Furthermore, the eigenstrain in the inclusion is assumed uniform. Therefore, the inhomogeneous inclusion of embedded BaTiO 3 nanoparticles has been solved by transforming the inhomogeneous inclusion into a homogenous inclusion. The nanoparticle with residual strain embedded within the epoxy was transformed into a homogenous inclusion problem with a distribution of equivalent eigenstrain. where δu i is the virtual displacement, δε ij is the virtual strain, f i is the body force, T i is the traction force on the boundary, and σ ij is the stress field in the nanoparticle. Equation (31) minus Eq. (30) gives: As the stress state in both conditions are identical, therefore, f i = f * i , and T i = T * i . The δε ij can take any arbitrary value, and Eq. (32) is then expressed as follows: From the stress expression in Eq. (25), the stress state is obtained as Eq. (34) over Ω, and Eq. (35) over Ω 0 : where C ijkl and C * ijkl are the elastic constant tensors of the homogeneous and inhomogeneous inclusion, respectively, ε kl is the total strain, ε * kl is the real residual strain in the inhomogeneous Ω, and ε 0 kl is the equivalent eigenstrain in the homogeneous Ω 0 . Since f i = f * i and T i = T * i , the deformation of the nanoparticle in epoxy in Ω is equivalent to that of the virtual homogeneous inclusion in Ω 0 by incorporating an equivalent eigenstrain ε 0 kl . The substitution of Eq. (33) with Eqs. (34) and (35) gives: The equivalent eigenstrain in Ω 0 is then expressed by: Therefore, with the equivalent eigenstrain ε 0 ij , the stress and strain field in the nanoparticle over Ω could be calculated using Green's function method as follows 47,56 , accounting for the displacement field for a 3D inclusion: where G ij,l X − X ′ is the Green's function, and www.nature.com/scientificreports/ The strain field is then obtained from Eq. (38), i.e.
The total strain in Eq. (37) is then expressed equivalent to the eigenstrain ε 0 ij given by Eq. (40), which gives: Therefore, the corresponding deformation field in the inclusion and the matrix can be calculated using Green's function method with the equivalent eigenstrain ε 0 ij . As illustrated in Fig. 5, a BaTiO 3 nanoparticle is embedded in the epoxy matrix, assumed to be in full bonding with the epoxy (i.e. no disbond or interfacial defect) and the epoxy is assumed to be infinite. The BaTiO 3 nanoparticle and the epoxy are assumed to be isotropic. A nominal eigenstrain ε * ij , due to the DW movement is prescribed as shown in the figure. It is assumed that the eigenstrain is the collective displacements of all the domains that may exhibit different polarisation in all possible directions.
The nominal eigenstrain is denoted as ε * ij : The anisotropic elastic constants of isotropic materials are C ijkl , and C * ijkl for Epoxy and BaTiO 3 respectively: where ( * ) and µ(µ * ) are Lame constants 57 , v is the Poisson ratio, and δ ij is the Kronecker delta. Inserting Eqs. (43) and (44) into Eq. (37): The equivalent eigenstrain ε 0 kl and total strain can then be re-written as ε kl :  www.nature.com/scientificreports/ where ε is a function of ε 0 56 : Therefore, Eq. (47) is: And the equivalent eigenstrain ε 0 is solved by Eq. (50): Therefore, the strain ε kl in the epoxy is expressed by: Interaction theory of BaTiO 3 -epoxy nanocomposite and microwave field. The major interaction between the dielectric material and the microwave field is the polarisation of the dipoles, i.e. under an oscillatory microwave field, the dipoles re-orient themselves with the external electric field component of the microwave field. Therefore, the BaTiO 3 nanoparticles interact with the electric field component of the microwave field, and the magnetic field effect is considered negligible. The main contribution of microwave power loss in the multidomain BaTiO 3 nanoparticles originates from the piezoelectric domain wall movement 58 . Compared with the BaTiO 3 nanoparticles, the dielectric loss tangent of epoxy represents the ability to dissipate the absorbed energy is relatively smaller (in the range of 10 −3 to 10 −4 ) 59 . The mobility of the polar segments of polymer decreased significantly when the frequency reaches the GHz range 60 . Therefore, the dielectric loss in BaTiO 3 -Epoxy nanocomposites observed at a higher frequency (e.g. microwave at 2.45 GHz) is mainly attributed to the BaTiO 3 nanoparticles.
The behaviour of the electric and magnetic components of the microwave field is governed by Maxwell's equations coupled with boundary conditions 61 .
where H and E are the varying magnetic (A/m) and electric field strength (V/m), respectively, ω is the angular frequency (rad/s), ε is the complex relative permittivity, µ is the complex relative permeability, and ρ is the charge density. The first two equations, so-called Ampere's law and Faraday's law, correlate the changes in the electric or magnetic field. The last two equations are Gauss's laws, which represents the net magnetic flux out of a region. It must be zero while the net electric flux out of a region is related to the charge density within the region 62 . The power absorbed by the material in the microwave field is dependent on many factors such as shape, size, and location of the material. For a relatively thin material that fully penetrated by the microwave field, the power dissipated P (W/m 3 ) as heat in the material under microwave field is expressed based on Maxwell's Equation as follows 63 : where the first term µ 0 µ ′′ H 2 rms describes the dissipation of power due to magnetic field and can be neglected for a non-ferrite material, e.g. BaTiO 3 , ω = 2πf , f is the frequency of the applied microwave field (= 2.45 GHz), ε 0 is the dielectric permittivity of the free space, ε ′′ represents the imaginary part of the complex dielectric permittivity, and E rms is the root-mean-square value of the local electric field strength (V/m). Therefore, Eq. (57) is simplified to: However, the explicit evaluation of the absorbed power by the material inside the oven cavity is difficult to obtain due to the non-uniformity of the microwave distribution, and hence a varying electric field strength. The power absorbed by the material P absorb , is expressed as follows: where P input is the input power from the magnetron, P reflected is the power reflected and attenuated, and P oven is the power absorbed by the cavity wall, waveguide, and magnetron. The increased temperature due to microwave heating also has a significant impact on the dielectric properties of the material, which in turn affect the degree www.nature.com/scientificreports/ of interaction of the nanocomposite with the microwave field. In previous studies, an empirical approach has been proposed to estimate the absorbed power. It is assumed that the amount of microwave energy absorbed by the material/field interaction is converted into heat during the exposure, and is expressed as follows 64,65 : where m is the mass of the material, C is the specific heat capacity of the material, T is the temperature change, V is the volume of the material, and t is the microwave exposure time.
The energy transferred into the material is also significantly determined by the depth up to which the microwave field travels into it. The degree of microwave penetration in the material is defined by a factor, Penetration depth ( d p ) 66,67 , which is the distance between the surface of the material and the place inside the material where the magnitude of the field strength drops to e −1 (36.8%) of that at the surface. The penetration depth d p neglecting the effects from the magnetic field for a non-ferrite material is expressed as follows 65 : where c is the velocity of light, ε ′ is the real part of the complex dielectric permittivity, and tanδ is the tangent loss of the material, tan δ = ε′′/ε′ . Generally speaking, the microwave power decays more quickly in lossy materials with higher tanδ due to energy dissipation hence a lower penetration depth. From Maxwell's equation, it is noted that the energy of the microwave field is proportional to the square of the amplitude, which represents the maximum strength of its electric and magnetic components. For a continuous sinusoidal microwave field in the air of an empty cavity, the average intensity I is the power carried by the microwave field per unit area, and is expressed as follows 68 : where E m is the maximum electric field strength, and the root-mean-square value of the electric field is: Equations 62 and 63 above could only be used to calculate the electric field strength in an empty oven. When a sample is placed in the multimode microwave oven (our case), the electromagnetic field is considered to be a superposition of multiple plane waves penetrate the material from different directions. Regardless of the differences in each mode, the superposition of all the waves could be assumed to have an approximately uniform and constant distribution pattern in the material within the localised region 69 . Assumptions are made that the absorbed microwave energy is converted to heating in the material with neglected surface heat loss, thermal diffusion, and energy absorbed by the oven (e.g. magnetron and cavity wall). Therefore, the energy absorbed by the material is corresponding to the temperature rise of the material, and it is expressed by re-arranging Eq. (60): It should be noted that the internal electric field is only localised in the material and is different from that in the surrounding cavity. Under the same assumptions that the material is surrounded by a uniform microwave field, incident microwaves are mostly transmitted and absorbed by the material while some are reflected. It is also assumed no reflected wave at the interface. An expression is presented as follows 70 : where V is the volume of the material, ρ is the density of the materials, c is the phase velocity of the microwave ( 3 × 10 8 m/s) 71 , A is the surface area, ξ is the absorption coefficient that represents the fraction of the absorbed power that generates heat, and Ŵ is the transmission coefficient. Solving the equation, the time-averaged strength of the external electric field E rms,external obtained by: The power absorption coefficient ξ is expressed as follows 72 : where P 0 is the incident power at the surface, and a is the thickness of the BaTiO 3 -epoxy sample. The transmission coefficient, Ŵ , is the energy fluxes per unit time and unit area at the interface. It is assumed that the microwave incident is normal to the material surface as the energy flux comes from all directions within the cavity, and the transmission coefficient Ŵ is simplified as 71 : www.nature.com/scientificreports/ where n 1 and n 2 are the refraction index of air and the material, respectively. The above equations are employed to theoretically obtain an estimate of the average and peak electric field strengths of the incident microwave with a given cavity size. Subsequently, the absorbed microwave power by the nanocomposite has been estimated.

Materials and experiments
Materials and fabrication. The epoxy used in this study was Araldite LY1564, a diglycidyl ether of bisphenol A (DGEBA) and the curing agent was Aradur 3487, an amine hardener, supplied by Huntsman, UK. This epoxy resin system has relatively low viscosity and high flexibility mainly for aerospace and industrial structural composites parts. The coupling agent for surface functionalisation selected in this study was 3-glycidoxypropyl trimethoxysilane ( Table 1, below 73 . All the chemicals except BaTiO 3 powders were used as received without further treatment. The phase composition of BaTiO 3 powders has been confirmed by X-Ray diffraction (XRD) using a Siemens D5005 X-Ray Diffractometer with Cu-Kα radiation (wavelength λ of 1.54 Å, at 40 kV and 30 mA). The scanning speed was 0.5°/s in the range of 20° to 90°. The peaks analysis was performed by data analysis software Origin using a Lorentzian curve fit function. Figure 6 shows a well-defined perovskite structure of BaTiO 3 without a noticeable second phase. The split peaks between 44°and 46° of the 2θ angle correspond to reflections (002) and (200) confirms the tetragonal phase of BaTiO 3 74 .
Surface functionalisation of BaTiO 3 and characterisation. The BaTiO 3 powders were prepared using combustion method, therefore, firstly they were pre-treated in H 2 O 2 for hydroxylation process to add hydroxyl group (-OH) to the surface 75   www.nature.com/scientificreports/ oven at 80 °C for 24 h. The reflux and particle retrieving processes were similar as in the surface functionalisation with 3-GPS as shown in Fig. 7. 3-GPS was then applied to BaTiO 3 nanoparticles to improve the processability and filler dispersion in nanocomposites; the solution of 1 wt.% of 3-GPS with respect to BaTiO 3 was prepared. 150 mL aqueous solution of ethanol and deionized water (9:1) was firstly mixed in a beaker. Adding the acetic acid drops using a pipette and stir vigorously after each drop until the pH value of 3.5-4 measured by a MET-TLER TOLEDO pH meter was reached, stirred vigorously again for 3 min to form a clear solution. The low pH values of the solution facilitate the silane functionalisation process 77 . After the addition of the 0.1 g 3-GPS solution to the acidified solution using a pipette, the mixture was left in an ultrasonic bath for 30 min to form a homogenous solution. 10 g hydroxylated BaTiO 3 powders was then added to the silane solution, and mixed under ultrasonic bath for 10 min for better filler wetting. Finally, the mixture was refluxed at the boiling temperature of ethanol, 78 °C 78 , at 100 rpm using a mechanical stirrer for 6 h using a silicone oil bath over a hotplate. After refluxing, the BaTiO 3 was washed three times with deionized water and retrieved using centrifugation at 4500 rpm. The silane treated BaTiO 3 (Si-BaTiO 3 ) powders were dried at 110 °C for 24 h to avoid any condensation of silanol groups at the surface. In the end, the powders were crushed in a mortar and pestle for the nanocomposites preparation. The size distribution of Si-BaTiO 3 particles was analysed in its epoxy nanocomposite form. The surface functionalisation didn't significantly affect the particle size as shown in the SEM images and it is approximately 200 nm (discussed in "Scanning electron microscope (SEM)." section). The silane-treated powders were characterised by Thermogravimetric analysis (TGA) and Fourier transform infrared (FTIR) spectrometry to confirm the presence of the functional group in the following sections.
Fourier transform infrared (FTIR) spectrometry. FTIR is performed to detect the silane functional groups on the modified BaTiO 3 powders in transmission mode using a Jasco FT/IR-6200 in the range of 400-4000 cm −1 with a resolution of 2 cm −1 at room temperature. The FTIR spectra of Si-BaTiO 3 and untreated powders are presented in Fig. 8. The peaks at 1437 cm −1 and 1630 cm −1 are due to the small trace of BaCO 3, and physically absorbed water on BaTiO 3 powders from the combustion fabrication process according to the manufacturer 79  Thermogravimetric analysis (TGA). The amount of organic silane group grafted onto the surface of BaTiO 3 powder was obtained by TGA using a TA instruments Q-500 at a rate of 20 °C/min under a dry nitrogen gas flow rate of 100 ml/min to 700 °C, shown in Fig. 9 for untreated and Si-BaTiO 3 powders. The weight loss percentage indicated as the lowest point of the green line is 0.35% of untreated BaTiO 3 and 1.2% for treated. Therefore, the amount of 3-GPS grafted onto the BaTiO 3 surface was calculated to be approx. 0.85%, which is close to the previously designed 1% in the surface functionalisation process. The weight curve of 3-GPS treated BaTiO 3 powders showed three well-defined degradations as labelled in the figure. Based on the FTIR results, the absorbed water is removed between 50 to 150 °C, (1) and (2) in the figure. The degradation step from 240 to 450 °C represents the free hydroxyl groups (-OH) and/or silane molecules that remained on the surface, labelled as (3). The pyrolysis of surface silane compounds was observed at around 565 °C, labelled as (4) 79,81 . The increase in weight between 650 and 700 °C is only observed in the untreated powder and it is attributed to the impurities in the as-received powder (supplied from US Research Nanomaterials Inc) reacting with nitrogen. www.nature.com/scientificreports/   www.nature.com/scientificreports/ Fabrication of BaTiO 3 nanocomposites. The epoxy resin nanosuspension with 1, 5, 10, 15 wt.% Si-BaTiO 3 nanoparticles and 5, 10, 15 wt.% untreated BaTiO 3 nanoparticles were prepared as follows: Firstly, the weighed amount of BaTiO 3 powders was added to ethanol and sonicated with an ice water bath for 2 min with a 10 s pulse to form a homogenous solution. Then a weighed amount of epoxy and the previous mixed solution were blended in a beaker using a mechanical stirrer at 300 rpm and 80 °C overnight under the fume hood to gradually remove the ethanol without precipitation of the particles. The mixture was weighed before and after the previous step to ensure the complete removal of ethanol. The curing agent was then added to the mixture with a weight ratio recommended by the company and stirred for a further 3 min. Finally, the mixture was placed in the vacuum oven at 30 °C for 1 h to remove bubbles at 29 inHg and achieve complete removal of ethanol. The whole mixture was poured into a mould made of two pieces of glass clamped with a 3 mm silicone gasket in-between, as shown in Fig. 10. A uniform thickness of each sample was achieved with the assistance of this type of glass mould. They were then cured in the oven for 8 h at 80 °C as prescribed by the manufacturer, then cooled down to room temperature. The final samples were of size 160 × 140 × 3 mm 3 , and cut to different sizes using a precision cut-off machine BRILLANT 220.
Field induced experimental methods. General material characterisations. The morphological characterisation of the nanocomposites, post cure, with different weight loading between 0 wt.% and 15 wt.% were carried out using a field emission scanning electron microscope TESCAN Vega 3. A thin layer of gold is coated to each sample before testing to enable better conductivity. The thermal characterisation of the cured nanocomposites was studied using a differential scanning calorimeter (DSC Q200, TA Instruments). The DSC measurements were performed in a nitrogen flow (50 mL/min). Samples 5 × 3 × 1 mm 3 were prepared using a Buehler ISOMAT low speed saw. The samples were first heated from 40 to 220 °C at 20 °C/min. The measured DSC curves of all epoxy nanocomposites presented a distinct exotherm peak. Thereafter, the second run at a heating rate of 20 °C/ min was conducted up to 220 °C after cooling down from the first scan. The glass-transition temperature ( T g ) was determined using the inflection-point method following the IS/DIS 11,357-2. A tangent was drawn at the point of inflection in the DSC spectra from the second run, and the T g was the midpoint between onset and end-set of the drawn tangent. The complex relative permittivity and permeability of the neat epoxy and the nanocomposites at different weight loading were measured with increasing frequency from 0.25 to 4.50 GHz at a 0.25 step increment using a vector network analyser (VNA) with microstrip methods 82 . Flat plates samples of size 50 × 50 × 3 mm 3 were prepared as substrates for conductive lines as shown in Fig. 11. They were then attached to a ground metal plate at the bottom. The dielectric constants of the samples were then measured using the VNA. The resonant approach can provide a more accurate values of the dielectric properties compared with the transmission line method. However, the resonant methods could only measure the dielectric properties at a single frequency instead of measuring in the frequency range 83 . www.nature.com/scientificreports/ Field-induced strain measurements. The microwave field was selected to be the external stimulator as BaTiO 3 exhibit a peak in a dielectric loss at the microwave frequency range. Most importantly, the design of the experiment becomes more feasible with remote stimulation from the microwave field owing to the spacious cavity. The mechanical strain evolution in the nanocomposite was investigated under a microwave field within a temperature-controlled microwave cavity Panasonic NN-SF464MBPQ running at 2.45 GHz and cavity size of 354 × 338 × 230 mm 3 . In contrast with the conventional one, this oven, equipped with inverter technology, has a circuit board that replaces the transformer, hence the output power can be adjusted linearly by varying the pulse width to ensure a more precise and continuous microwave exposure 84 . The unique flat-bed design of this model is equipped with a stationary ceramic plate that allowed more space to place the sample and its holder. The strain field introduced by the BaTiO 3 nanoparticles to the surrounding epoxy, activated by microwave field, was studied by real-time strain measurements on the surface of the nanocomposite samples with the incorporation of Fibre optic sensors utilising fibre Bragg gratings (FBG) technique. The sensors were placed apart in equal distance from one another within a 90 mm straight line. This is theoretically aligned with the microwave's half wavelength (~ 60 mm) as also was experimentally observed during real-time temperature measurements. The sensors were located to ensure overlapping with at least three nodes and antinodes of the microwave cycle. The FBG sensors were fabricated by a periodic intense laser light applied onto the core of an optic fibre. The laser light exposure introduced a permanent increase in the refractive index of the fibre's core and a fixed modulation was created subsequently. Each FBG was approx. five mm long with a grating period of one micron. Two arrays with three FBGs for each sample were fabricated. FBG arrays were then adhesively bonded to the surface of two geometrically identical 135 × 10 × 3 mm 3 samples made of Si-BaTiO 3 epoxy nanocomposites. The strain array was bonded to the surface by adhesives at the FBG regions while the temperature array is firstly packaged into a capillary glass tube, and then bonded parallel to the strain array. The arrays were placed at the same location with the same distance in-between. The strain experienced by the samples was transferred to the FBG sensors, and the measurements were also affected by the temperature. Therefore, the temperature array only measured the temperature change, and compensated the strain measurements accordingly.
A preliminary test via a FLIR One Pro LT Thermal Camera was carried out to inspect the temperature change in different samples under different microwave power levels. Samples were placed at a designed location as shown in Fig. 12. The exposure time was then carefully selected based on the heating profiles of each sample to control the temperature of samples well below T g during the exposure to avoid any interference in strain measurements from the post-cure shrinkage. The sample was clamped on one end by a designed sample holder made of polytetrafluoroethylene (PTFE) to ensure minimum interaction with the microwave field owing to its extremely low dielectric loss. The size and location of the load inside the cavity were two primary factors that affect the microwave field distribution 84 . To control the microwave field distribution, the sample and the holder were placed at the designed location shown in the figure, for all tests. The power level of 100 W and 440 W were initially selected to avoid temperature surges in samples within a short period of time. Exposure time was set to be 650 s for 100 W and 148 s for 440 W to limit the temperature below T g (80 °C) based on the preliminary results obtained by the thermal camera.
Neat epoxy, BaTiO 3 nanoparticles, and adhesive used for bonding FBGs possess different thermal expansion coefficients (CTE) and microwave heating patterns. Microwave field interaction with the neat epoxy and adhesive have been investigated, thereafter. Two arrays with five FBGs of five mm long and one micron grating period were fabricated for each test. First, the adhesive response under microwave radiation was studied by two arrays of FBGs adhesively bonded to the surface of a PTFE block. It is assumed that measured strain and temperature change are solely due to the microwave heating of adhesive as PTFE has neglectable thermal response (temperature rise) under the microwave. The PTFE block is placed in the microwave oven at the designed location to locate the first FBG sensor on the left end at the same location as the one in the previous exposure on the nanocomposites. The tests have been performed under 100 W for 650 s and 440 W for 110 s. Data was recorded on the interrogator 10 s prior to the microwave exposure for each run to ensure no data is missing after the microwave starts. The effect of neat epoxy with the FBG sensors bonded to its surface was investigated as a controlled group to the test of 15 wt.% nanocomposites at 10 W for 600 s and 440 W for 60 s. Neat epoxy sample geometrically identical with 15 wt.% nanocomposites in the nanocomposite's exposure test was placed in the designed location with surfacebonded arrays of FBGs. The whole measurements from adhesive and the neat epoxy were compared with the www.nature.com/scientificreports/ results from the previously performed test on the nanocomposites for better distinguishing the BaTiO 3 -epoxy nanocomposite's response to the microwave exposure.
Raman spectroscopy. Raman characterisation was carried out to investigate residual stress, any indication of remnant strain field or crystal structure change in BaTiO 3 nanoparticles embedded epoxy post microwave exposure. Neat epoxy samples and 15 wt.% BaTiO 3 epoxy nanocomposites with the identical geometric dimension of 135 × 10 × 3 mm 3 were characterised before and after 440 W microwave radiation of 110 s. The Raman spectra were obtained using a Raman spectroscope Horiba Scientific LabRAM HR with a 514 nm excitation wavelength laser for all samples.

General material characterisations. Scanning electron microscope (SEM). A silane coupling agent
3-GPS was introduced to modify the BaTiO 3 particles in this work to achieve a finer dispersion in epoxy. To examine its effectiveness, the dispersion and distribution of the nanoparticles at epoxy with different loadings have been characterised using the SEM. Different levels of aggregation events due to increased weight fraction of particles and enhancement of distribution due to surface functionalisation were present: Figure 13a-f presents the SEM images of the fractured cross-section of the silane functionalised and nonfunctionalised BaTiO 3 epoxy composites at various wt.%. It is shown that the nanoparticles are uniformly dispersed with the aid of surface functionalisation. As shown in Fig. 13a, the silane-treated particles are of the similar size of 200 nm as pristine BaTiO 3 particles. There is no apparent aggregation of Si-BaTiO 3 particles until small-scale clusters become visible in Fig. 13c,d when the loading of BaTiO 3 reaches 10 wt.% and over, due to nanofillers agglomeration as circled in Fig. 13c,d. Overall, the majority of the Si-BaTiO 3 particles were finely dispersed compared with the untreated ones. Figure 13e,f show images of the composites with 15 wt.% nonfunctionalised BaTiO 3 particles at different magnification. In contrast with the functionalised BaTiO 3 at the same weight loading 15 wt.% shown in Fig. 13d, the 15 wt.% untreated BaTiO 3 exhibit severe agglomeration to form large clusters that have left large-scaled lamellar structures. The efficacy of the silane functional group is in good agreement with works conducted by others [85][86][87] .
DSC. The DSC was conducted to examine the effect on curing extent and T g of cured epoxy nanocomposites with varying weight fractions of Si-treated and untreated BaTiO 3 nanoparticles. Figure 14 shows the DSC spectra of the neat epoxy and nanocomposites at various BaTiO 3 content as an attempt to provide a comparative analysis between the variation of Tg in the different multi-material systems. www.nature.com/scientificreports/ As observed from the insets, the DSC spectra of epoxy nanocomposites exhibit insignificant bumps denoted as exothermic peaks representing the post-curing process due to the impeded curing by the existence of BaTiO 3 nanoparticles, i.e. indicating the nanocomposite material is well cured with the current weight loadings of the BaTiO 3 particles. The T g transition for the coloured lines that represent nanocomposites are shifted slightly to lower temperatures compared with the black line of neat epoxy which is evidence of reduced T g as the nanoparticles hindered the curing process of the epoxy resin 88 . Furthermore, the weak endothermic peaks that emerge on the spectra represent the relief of the stress introduced during processing and handling or the thermal history. The thermal history and residual stress were removed from the first run of heating and the samples were postcured during the first DSC analysis up to 200 °C. The true T g of each cured nanocomposite sample was calculated based on the DSC spectra of a second run under the same condition. Figure 15a presents the DSC spectra of a second heating scan for each epoxy nanocomposite with functionalised BaTiO 3 after cooling down from the 220 °C previous runs. The glass temperatures per nanocomposite category have been identified in Fig. 16. It can be observed that epoxy nanocomposites are well cured without any perceptible exothermal peak. The red and black lines are drawn as a demonstration of T g variation along with   www.nature.com/scientificreports/ the weight fraction of the nanofiller. The red line shows a gradual increase of T g with the increasing Si-BaTiO 3 content. The T g has increased up to 90 °C compared with that of the neat epoxy, 87 °C, as shown in Fig. 17. The increased T g represents that more energy is needed from higher temperatures to complete the glass transition process, which indicates that the movement of polymer chains is inhibited by the functionalised BaTiO 3 nanoparticles due to improved interfacial bonding. However, by adding the non-functionalised BaTiO 3 , the T g shows a remarkable drop at 5 wt.% at the beginning as marked by a black line in Fig. 15b. When the weight loading reaches 10 wt.%, the T g shows an increment surge from 81 to 84 °C. The T g of the nanocomposites reaches 89 °C at 15 wt.%. The untreated particles at 5 wt.% have a negative effect on the Tg by obstructing the formation of the cross-linking structure of epoxy due to particle agglomeration as shown previously in the SEM images. However, the addition of the nanoparticles can enhance the mobility of the polymer chain by introducing increased free volume at the filler-matrix interface 89 . For the uncured pure epoxy resin, the mobility of the polymer segments decreases and Tg increases with the temperature rise during the curing process. When the Tg finally approaches the curing temperature and becomes higher, the curing process ceases due to the lack of mobility of the segments. The free volume created by introducing more untreated nanoparticles increases the mobility of the segments during curing, and hence favours the curing process and higher Tg. and 15 wt.%. As discussed formerly, the ferroelectricity character of tetragonal BaTiO 3 enables a spontaneous polarisation that can be reversed by an applied field 94 . Therefore the ε ′ of epoxy is enhanced due to the addition of BaTiO 3 nanoparticles, which agrees with the findings from other researches 90,91,93,95 . Generally speaking, a substantial increase in the dielectric properties of modified epoxy is only achieved by a generous amount of BaTiO 3 content as high as 90 wt.% 96 . The blue line that represents the dielectric response of 5 wt.% silane-treated BaTiO 3 epoxy nanocomposites, exhibiting clearly a significantly lower value of real permittivity ε ′ below the neat epoxy's. This drastic decline in the value of ε ′ is attributed to the testing errors from airgaps that were identified due to the uncertainties associated with handling the samples, between the sample and conductive line 97 while was not observed in the other samples. Figure 17b illustrates that the loss tangent in all samples increases with increasing frequency. On the contrary with the real permittivity, the loss tangent has slightly decreased with adding BaTiO 3 powders except from 1 wt.% content. The loss tangent of a composite system has contributions from dipole orientation, conduction loss and interfacial polarisation 98 . Silane treatment of the BaTiO 3 powders introduce a reduction of the concentration of ionizable hydroxyl groups (-OH) on the BaTiO 3 surface and hinders the mobility of the charge carriers on the surface. The lowered conduction loss could be the primary reason for lowered loss tangent 80 . are due to Bragg wavelengths moving into adjacent spectral windows that were set up on the sensor interrogator. The strain and temperature data exhibit a general increasing trend with the exposure. They drop gradually, immediately, after the microwave stopped. The temperature appeared to be higher at both ends of the specimen (sensor 1 and 3, located left and right respectively) compared to the middle (sensor 2), which is in accordance  www.nature.com/scientificreports/ with the 'hot spot' theory due to microwave nonuniformity as illustrated in Fig. 21. Sensors 1 and 3 data also follows more similar temperature increasing trend (magnitude and rate). The microwave oven cavity is designed to be a multimode resonant cavity, similar to a closed waveguide. The microwaves generated by the magnetron are reflected back and forth by the metal walls inside the cavity and eventually forms standing waves throughout the cavity's volume as illustrated in Fig. 21 99 . The illustration is accompanied by a thermal image captured from the sample under the microwave, showing the approximate locations of the nodes and anti-nodes. As seen, the antinodes (hot spots) are located at or surrounding the sensors. Under microwave exposure, the distance between the area with higher temperature (antinodes) and the cold spot with lower temperature (nodes) is a quarter of the wavelength, which is approximately 31 mm 100 . The sensors have been positioned equal distances apart, and accordingly, they fell in different temperature zones under the standing waves, identified schematically in the figure (assuming that the far-left edge of the sample represents a node). The strain measured from a FBG sensor is derived from the change in the wavelength of reflected ultraviolet light due to grating period change. The change with strain and temperature is expressed by 101 : where is wavelength shift, 0 is the base wavelength, k = 1 − p , p is the photo-elastic coefficient, p =0.22, ε is the strain due to both mechanical and thermal factors, T is the temperature change, and α δ is the change of the refraction index. As expressed in Eq. (69), the measured data from the strain sensor are a combination of two factors. The first part is the actual strain of the material including any mechanical strain and strain due to www.nature.com/scientificreports/ temperature changes. The second part is the wavelength shift due to a change in the glass-refraction index of optical fibre with temperature. According to the FBG sensors manufacturer, 1 °C temperature variation induced strain of optical fibre itself is approximately equivalent to eight micro-strain. Under the 100 W and 440 W, the strain measured at the beginning of the microwave radiation in strain sensors 1 and 3 have similar trends of increasing as predicted when the temperature sensors 1 and 3 measurements have close gradients, having higher temperatures than that measured by sensor 2. A sudden drop in the strain measurements of the middle sensors (strain 2) can be observed in both cases soon after the initial surge, labelled as A-B. It instantaneously alters the strains by nearly identical 1162 and 1008 micro-strains (difference between A and B) under 100 W and 440 W, respectively. Such decline in the strain measurements is attributed to an immediate development of a compressive strain in response to the microwave exposure, as hypothesised in Sect. 2. Note that sensor 2 temperature is the lowest amongst the three sensors, in which the thermally induced strain is negligible. Such phenomenon could not be associated with the temperature rise since the temperature change is approximately 12 °C and 5.7 °C from room temperature 19 °C when the sudden drop occurs at point A, which is remarkably lower than 170 °C when post-curing shrinkage is introduced as indicated by the exothermic peak in the DSC spectra in Fig. 14, and lower than the T g to have any detrimental effect. Moreover, it is observed that unloading the specimen from the 440 W exposure results in a residual compressive strain in all sensors' locations. This occurs after the high non-linear variation of strains in the case of 440 W, unlike the linear strain behaviour under 100 W. This is analogous to mechanical field introduction in which unloading an elastic-plastic specimen beyond its elastic regime under tensile loading may introduce a compressive residual strain distribution depending on the hardening behaviour (i.e. kinematic or isotropic). Further investigation was conducted on the phenomena observed, described as follows: Different FBG arrays have been adhesively bonded at the same location on the surface of two specimens made of 15 wt.% with identical geometric dimensions and the location of the two specimens are placed at the exact location inside the microwave cavity. Accordingly, the distribution of microwaves in both cases is considered to be identical, also verified using thermal camera. The strain and temperature change measurements of sensor 2 under 100 W and 440 W before and after the 'sudden drop' is presented in Fig. 22a-d. It can be observed from Fig. 22a,b that the measured strains demonstrated a steady step-wise increase with the increasing temperature that fluctuated within a constant range of approximately 8 °C. The fluctuations in temperature change exhibit a periodic pattern that is approximately in phase with every step of the increment in strain. Under 440 W as shown in Fig. 22c,d, the increase in the strain and temperature data follow a similar trend as that under 100 W but at a higher rate.
In contrast with the middle sensors, the strain and temperature increments from sensor 1 under 100 W and 440 W are predominately synchronized as presented in Fig. 23. It could be observed from Fig. 23a,c that the cycle time for each step is approximate 20 s and each step is approximately 100 micro-strain. A similar trend is also observed in the first 50 s of the microwave exposure under 440 W from Fig. 23b. However, the cycle time for each step shown in Fig. 23d falls sharply to 2 s accompanying a higher increment of 300 micro-strain compared with the 100 W data which are attributed to the increased heating rate at higher microwave power.
Adhesive FBG data: www.nature.com/scientificreports/ Despite the widely known utilisation of FBG strain sensors, there exists a layer of adhesive and protective coating that may affect the data from the nanocomposite 102 . The FBG sensors are adhesively bonded to the surface of specimens by PERMABOND 920 Cyanoacrylate. The dielectric and thermal properties of the adhesive and neat epoxy are presented in Table 2.
Five FBG sensor arrays, 20 mm apart, were utilised to quantify the adhesive's response distribution. Figure 24 presents the strain and temperature change measured via sensors bonded to a PTFE block with adhesive under 100 W and 440 W. It can be observed from Fig. 24 that temperature increases in the adhesive with microwave exposure time at both power levels as predicted while the temperature rise has a higher rate under 400 W. The maximum temperature changes are around 7 °C and 4.5 °C for 100 W and 440 W, respectively. The insignificant amount of temperature rise is the direct evidence of relative less interaction of adhesive and microwave. The strain values steadily increase with the temperature rise in both cases due to the thermal expansion of PTFE and adhesive. The calculated thermal strain of PTFE based on Fig. 24 and Table 2 under 100 W is approx. 868 micro-strain at 7 °C temperature change. The calculated value is almost consistent with the measured value of strain sensor 1 which is approx. 960 micro-strain at similar temperature variation range. Unlike 100 W, the measured strain of sensor 1 in 440 W shown as the blue line (when microwaves are removed) is approximately 10 micro-strain, which is much lower than 458 microstrain, the calculated thermal strain of PTFE at the same temperature change (3.7 °C).
At the microwave exposure, the adhesive solely interacts with the microwave, and is heated up owing to its higher dielectric loss compared with the PTFE. However, the PTFE have a higher CTE than the adhesive. Consequently, the measured strain under 100 W is mainly attributed to the thermal expansion of PTFE from the conduction of heat. The strain sensors under 100 W showed strain increase initiated after 200 s, indicating that a substantial amount of time is required to trigger the heating conduction between the adhesive and PTFE. On the contrary, 440 W has a higher rate of temperature rise within 110 s exposure. Based on the analysis of the calculated and measured strain, the strain under 100 W is mainly attributed to PTFE thermal expansion while the strain under 440 W is dominated by the thermal expansion of adhesive. Under 440 W, a similar trend in temperature rise is noted between sensors 1 and 3, as well as sensor 2 and 4, which is described by the 'hot spot' theory due to non-uniform microwave field.
Neat Epoxy FBG data: www.nature.com/scientificreports/ Strain and temperature evolution of the neat epoxy with adhesively bonded arrays together under microwave are presented in Fig. 25. The microwave running time has been estimated to maintain the temperature well below the T g of the BaTiO 3 -epoxy nanocomposite. The temperature increases as a function of the exposure duration (equivalently, microwave energy). Sensor 5 data present a higher rate of temperature variation compared with the other sensors conforming to the 'hot spot' theory. The region in which sensor 5 is located within the zone with a higher rate of temperature rise that is denoted as a 'hot spot' zone presented in Fig. 26, hence the next hot spot shown in Fig. 26 (red square symbol) is half (60 mm) of the microwave wavelength (120 mm) is located between sensors 3 and 2. As schematically demonstrated, the rest of the FBG sensors are located out of the hot spot region, and therefore they develop similar rates of temperature under microwave.
Strain and temperature response of the neat epoxy with bonded sensor arrays is investigated under microwave (Fig. 25). The temperature increases as a function of microwave exposure time and starts to drop as soon as the microwave stops as predicted. No sudden drop such as that observed in the nanocomposite is observed for the neat epoxy which indicates the phenomenon linked with the embedded BaTiO 3 . Sensor 5 presents a higher rate of temperature change compared with the other sensors, corresponding to the 'hot spot' theory discussed formerly, and illustrated in Fig. 26. As schematically shown in Fig. 26, the other FBG sensors are located out of the hot spot region, and therefore develop similar rates of temperature variation under microwave.  www.nature.com/scientificreports/ The temperature change in both cases reached the highest value of 30 °C compared with that measured for the adhesive only. The measured strain follows a increasing trend with the temperature change, which is a resultant of the thermal expansion from the adhesive and neat epoxy. Maximum strain and temperature values measured by left-side sensor for adhesive on PTFE (sensor 5), adhesive on epoxy (sensor 5), and adhesive on nanocomposite (sensor 1) are tabulated in Table 3, summarised from the data shown in Figs. 19 and 20, under 100 W and 440 W exposure.
The measured strain of the first two categories (adhesive on PTFE and adhesive on the neat epoxy) are temperature compensated which means the effect of thermal expansion from the optical fibre itself is eliminated. Approximately 8 micro-strain is achieved by a 1 °C change in temperature according to the manufacture. The true strain of the epoxy nanocomposites samples is then lower than the measured strain from sensor 1: For the case of 100 W exposure, the measured strain of adhesive on PTFE is mainly attributed to the thermal expansion of the PTFE, as formerly discussed. When the FBG arrays are adhesively bonded to a neat epoxy sample, the measured strain drops by approx. 52% compared to that (790 micro-strain) from the sensor adhesively bonded on PTFE at the same temperature variation of 5 °C. It is due to the lower CTEs of epoxy and adhesive compared to that of the PTFE. Furthermore, at this temperature change (5 °C), the measured strain without temperature compensation of epoxy nanocomposites loaded with BaTiO 3 drops 35% compared with the neat epoxy. However, as temperature change rises to 30 °C, the measured strain in both cases tends to be similar. Under 440 W, the value of strain is mostly dominated by adhesive as described earlier. When the temperature change firstly reaches 5 °C at the beginning of the microwave exposure, the neat epoxy and nanocomposites have similar strain (320 and 330 micro-strain, respectively).

Raman characterisation.
Raman spectroscopy was utilised to examine any residual/remnant strain induced by the particles to the surrounding epoxy matrix, or any crystal structure transformation of BaTiO 3 particles that might be introduced by the microwave exposure. The middle region of the samples where the 'sudden drop' in strain occurred was of particular interest, hence it was the region where the laser was focused to obtain the Raman spectra. Another run at microwave power of 1000 W for 50 s was performed for such effect.
Raman characterisations have initially been performed on the neat epoxy and BaTiO 3 -epoxy nanocomposite samples with 15 wt.% prior to microwave to establish baseline curves for identification of possible change in the position of Raman peaks for the samples post microwave exposure. Raman spectra of the epoxy are shown in Fig. 27a. Broader peaks at 600-800 cm −1 and 800-1000 cm −1 are presented in Fig. 27a for the neat epoxy   www.nature.com/scientificreports/ run at 1000 W for 50 s was performed with Raman performed at the same location as that at 440 W, presented in Fig. 27d. No conspicuous differences are observed, and it could be noticed that the 304 cm −1 still emerges without recognizable changing in peak height, which implies no changes in the crystal structure of BaTiO 3 after the microwave exposure, and no or slight remnant strain field.  Table 3 have not been corrected to account for the temperature effect of fibre refraction index, while the results  It could be observed from the figures that the measured strain of the nanocomposite (blue triangles) are generally smaller than that of the neat epoxy (red circles). Notably, the difference increases with increasing temperatures under 440 W.

Discussion
It is well known that the mechanical properties such as elastic modulus, Poisson's ratio, and shear modulus [110][111][112] and the CTE [113][114][115][116] are strongly dependent upon the microcracks in the matrix. Camahort et al. 117 correlated a reduction of the CTEs of thermal cycled composite laminates in the specimens with microcracks identified under the microscope. A theoretical study of the CTEs in cracked laminates demonstrated that CTEs change as a function of the internal stress field 115 . This discrepancy in CTEs under different internal stress fields is due to the fact that the internal stresses induced by thermal cycling could be either tensile or compressive. When the stress field is tensile, the microcracks further opened up and cause a reduction in the value of CTE. Therefore, this correlation between the microcracks and CTEs might provide an explanation of the reduced strain measurements presented previously. Under microwave, domain wall movements are triggered in the BaTiO 3 nanoparticles, hence stress fields around the microcracks might be changed and lead to a different CTE of the nanocomposite as illustrated in Fig. 30.
It is found by many studies that the CTE of composite material is highly dependent upon its constituents 118,119 . A decrease in the CTE of the nanocomposite material is observed and is attributed to a relative lower CTE of the fillers compared with that of the epoxy matrix  (2)]. The curing of BTO/Epoxy samples starts simultaneously with the temperature rising due to the microwave heating effect. However, under relatively lower temperatures, the curing effect is neglectable (exothermic peak can be seen from the DSC graph as curing is an exothermic reaction). The temperature of microwave heating is controlled below 50 °C which is below the exothermic peak around 170 °C. Therefore, the curing shrinkage effect is not significant during microwave radiation.
The average crystallite size of the purchased BaTiO 3 nanoparticles are determined by XRD results using Scherrer's equation 120 : where k is the dimensionless shape factor of the crystallite, is the X-ray wavelength, 2θ is the peak position, B(2θ) is the half-width of the maximum intensity of the peak and θ is the Bragg diffraction angle. A round-up value of 1.0 is selected in this study. The mean crystallite size is 35.9 nm, and it indicates that the corresponding grain size is 35.9 nm in the 200 nm BaTiO 3 nanoparticles as shown in Fig. 31. Note that the strain measured by the sensors is averaged over a 5 mm × 0.5 mm region (i.e. sensor size). Such sensing area covers approximately 890,000 200 nm-particles, or alternatively 70 particles over a 20 × 20 µm 2 area. www.nature.com/scientificreports/ The black square marks in Figs. 28 and 29 represent the theoretically calculated values of linear thermal expansion of the adhesive used for bonding the FBG arrays. It can be observed that the strains in the neat epoxy are slightly lower than those measured by the thermal expansion of the adhesive at an identical temperature. A schematic diagram of the cross-section of an FBG array bonded to the sample is presented in Fig. 32. As discussed for Figs. 24 and 25, the adhesive exhibits lower interaction with microwave, hence a poorer dielectric heating effect. Thus, when the exposure initiated, the neat epoxy is firstly heated up due to the microwave selective heating based upon its different dielectric loss. The energy is then transferred to the adhesive via heating conduction from the epoxy, which initiates a thermal expansion in the adhesive. However, the epoxy possesses a relatively lower CTE in comparison with that of the adhesive. Consequently, the thermal expanding rates are different due to the mismatch in CTEs. As such, the measured strain in the FBGs is a result of thermal expansion of the adhesive constrained (via bonding) by the epoxy, causing the FBG arrays to bend with increasing temperature and bend back when the temperature drops. This behaviour might also be an explanation of the 'zig-zag' step-wise increment observed in Fig. 23. When the temperature ramping rate is slower under 100 W, the heat transfer between the epoxy and the adhesive is steadier compared with that under 440 W. Subsequently, the strain measurements exhibit a clearer step-wise increment under 100 W.
It is also noteworthy that the unanticipated finding of a 'sudden drop' during real-time strain measurement of the middle FBG sensor (presented in Figs. 19 and 20) occurs at the beginning of the microwave exposure at 1045 micro-strain under 100 W and 1261 micro-strain under 440 W, which is a clear implication of a phenomenologically induced compressive strain in the multi-material system around the middle FBG sensor. The true reason for the occurrence of this phenomenon remains unclear, however since the temperature effects are excluded and the phenomenon is only seen in the nanocomposite (and not in the neat epoxy), such high extrinsic strain variation is strongly attributed to the collective domain wall movements due to the reorientation of dipoles at high frequency in each domain.
Effect of non-uniformity of the microwave cavity. As the microwave travels through the material, the microwave power decays due to the absorption by the material as it penetrates. As a result, the distribution of the microwave field within the material is non-uniform. For the 15 wt.% BaTiO 3 -epoxy nanocomposites, the penetration depth calculated by Eq. (61) is approximately 0.86 m, based on the measured dielectric properties  www.nature.com/scientificreports/ at 2.5 GHz. This value indicates that approximately 63% of the energy from the incident microwave is dissipated at this distance under the sample surface. Therefore, for a nanocomposite with 3-mm thickness, it would be a reasonable assumption to make that the microwave fully penetrates the sample. Apart from the non-uniformity caused by the penetration depth of the microwave, there is also another factor that affects the non-uniform electric field distribution in the material, which is the standing waves formed in the microwave oven, extensively studied [121][122][123] . When the microwave field transmits and reflects within the cavity of the microwave oven, the standing waves form, and the 'hot spot' phenomenon manifests. The distance between two 'hot spots' is half of the wavelength of the incident microwave, which is approximately 60 mm. Therefore, the temperature change measured by FBG sensors exhibit different heating rates along with the 130 mm sample.
Effect of absorbed microwave power. The power absorbed by the material is evaluated by the method of the energy balance of the microwave heating effect. At power 100 W, the average power per unit area in the 354 × 338 mm 2 microwave oven is calculated as follows: From Eqs. (62) and (63), the peak electric field strength E m and root-mean-square E rms of the incident microwave field is: Similarly, for the 440 W microwave power level, the peak electric field strength E m is 1.67 × 10 3 V/m, and the E rms is 1.18 × 10 3 V/m. The estimated electric field is the power absorbed by the air. However, it is assumed that when the microwave field generation is initiated by the magnetron, the electric field strength within the microwave oven is assumed to be equivalent to the value of the E rms as the dimension of the sample (4 cm 3 ) is relatively smaller than the cavity size (27,520 cm 3 ). The energy absorbed as a function of time calculated based on the E rms calculated above is estimated as follows: The strain and temperature measurements as a function of energy absorbed in Joules are presented below, under 100 W and 440 W.
It could be observed from Figs. 33 and 34 that, the energy absorption rate is close to the rate of temperature rise in the FBG sensors. At the beginning of the microwave exposure, the temperature rise exhibits a similar rate compared with the energy absorption rate. As the exposure time increases, the temperature rise becomes significantly slower. The 'hot spot' theory due to the non-uniformity of microwave fields is proposed to be the possible explanation of this phenomenon as described formerly.
The internal electric field inside the sample could only be estimated based on the assumption that 100% energy transfers into heat. Therefore, the external electric field of the incident microwave at the surface of the www.nature.com/scientificreports/ BaTiO 3 -Epoxy sample, E rms,external is obtained from Eq. (66), and employed to be the E rms to solving the Eq. (22) to estimate the eigenstrain of the BaTiO 3 .
where calculated ρ by given weight and volume is 1.3 g/cm 3 , C calculated from the DSC spectra is 1.238 J.g -1 .°C −1 , a is sample thickness 3 mm, c is the speed of microwave (3 × 10 10 cm/s), ε 0 is the permittivity of the free space (8.85 × 10 −14 F/cm), and ξ and Ŵ is calculated by Eq. (67), Eq. (68), as follows: where d p is calculated to be 850 mm, n 1 is 1 for air and n 2 = √ ε ′ = √ 4.42. At the beginning of the exposure when the microwave incident the sample at t = 10s under 100 W, the rate of temperature rise T t is obtained using the linear fit of the temperature change measurements from FBG sensor 1, which is 0.05.
The domain width d is highly dependent on the grain/crystallite size, therefore, it is estimated to be 4 × 10 −5 cm with the crystallite size assumed to be equal to the nanoparticle size 200 nm based on the measurements in previous studies 124 , and the strain ε 440 at t = 7s under 440 W with the rate of temperature rise of 0.3, is 144 × 10 −6 . This value could only represent the approximate amplitude of the induced DW displacement under such electric field that impinges the sample surface.
Post-microwave residual effect in the BaTiO 3 /epoxy nanocomposites. As discussed in Sect. 3, the electric field-induced strain is predominately attributed to the irreversible domain wall movement in multidomain ferroelectric crystals. Therefore, when the microwave field is removed, a residual stress field is hypothesised, and Raman spectroscopy was employed in this study to probe any possible residual strain effects that originate from microwave field activated BaTiO 3 nanoparticles. The Raman peaks shift indicated the changes in the atomic spacing of the material, namely the strain information 125 . In the meantime, it was also discovered that the altered stress state surrounding the BaTiO 3 could introduce a crystal structure change due to the piezoelectricity effect exhibit in ferroelectric materials. Raman data did not indicate a significant difference between the Raman spectra of 15 wt.% nanocomposite samples with and without microwave treatment. www.nature.com/scientificreports/

Conclusions
The quantitative research on microwave field-nanocomposite interaction was conducted to study the ferroelectric materials' response in high-performance (rigid) epoxy composite that offers reversible microwave activated electro-strains introduced by a second dielectric phase (BaTiO 3 ). The data presented a pioneering investigation as no investigations have been reported, thus far, on the micromechanical (extrinsic) strain response of rigidly constrained ferroelectric materials in polymer under external electromagnetic fields. The research has characterised the thermal and dielectric properties of ferroelectric nanoparticles (BaTiO 3 ) modified polymer (epoxy) at weight loading up to 15 wt.%, and subsequently investigated the extrinsic strain state induced by the BaTiO 3 nanoparticles' dipole displacements onto their surrounding epoxy, however observed at the phenomenological level in the multi-material system. Silane surface functionalisation of the nanoparticles was utilised to uniformly disperse and improve interfacial particle-matrix bonding quality. FBG sensors-based technique was employed, in-situ with the microwave exposure, for real-time monitoring of the strain and temperature response of the nanocomposite subjected to a 2.45 GHz microwave field at controlled exposure power and energy, followed by theoretical constitutive equations underpinning field-material interactions, and the remnant extrinsic strain investigations using Raman spectroscopy. The key remarks on the investigations conducted in this article are as follows: 1. A time-dependent strain field is introduced, proportional to the exposure time/energy, in the multi-domain BaTiO 3 nanoparticles under the stimulation of a microwave field, linearly under 100 W and non-linearly under 440 W. The strain is theoretically described via development of multi-physics constitutive expressions for the nanoparticle possessing a multi-domain crystallite based on an estimated electric field strength exhibiting at the geometric boundary of the particle. 2. The principle of virtual work was utilised, phenomenologically for the theoretical development, to estimate the strain within the epoxy matrix embedded with BaTiO 3 subjected to the microwave exposure, able to underpin the in-situ induced micro-straining performance of the BaTiO 3 -epoxy nanocomposite. 3. A sudden drop in micro-strain data by the value of approx. 1000 compressive micro-strains was observed in the BaTiO 3 -epoxy composites at the beginning phase of the exposure under the different exposure powers examined. Such phenomenon has not been reported thus far in the existing literature, however investigated herein using different characterisation and theoretical development at micromechanical level. The drop indicated an instantaneously induced compressive strain achieved within a small fraction of a second, which is recovered afterwards with the increasing exposure energy, followed by a gradual increase during the exposure. No arcing in the microwave was observed, and no damage in any part of the FBG, adhesive, composite and PTFE; The temperature variation due to the exposure was kept well below the T g , and the DSC investigations presented a fully cured nanocomposite. The theoretical-experimental study of the response led to the fact that such effect is caused by the DW movement response of the BaTiO 3 -epoxy composite under microwave exposure, which is attributed to a mechanical compressive field to its surrounding polymer, with its effect phenomenologically observed by the in-situ FBG technique. 4. The nanocomposite exhibited linear and highly non-linear strain variations under the 100 W and 440 W exposure, respectively. Residual compressive strains were developed post 440 W exposure, offering an analogous response to the mechanical field pre-straining which takes the strain hardening materials (e.g. polycrystalline) beyond its elastic regime. Such strain hardening response in the case of 440 W is attributed to the semi-crystalline phase that BaTiO 3 had introduced to the nanocomposite. 5. To achieve the full benefit from the electro-strain behaviour of the BaTiO 3 fillers, uniform dispersion of the nanoparticles in the epoxy matrix was considered to be the primary determinant factor. Large agglomerates were observed in non-functionalised BaTiO 3 -epoxy nanocomposites due to the attraction between the nanoscale particles, and their intrinsic incompatibility with hydro-carbon epoxy matrix. Uniform dispersion of the surface functionalised BaTiO 3 nanoparticles ensures a uniform application of the domain wall movements triggered stress field to their surrounding epoxy matrix. Furthermore, it also avoids the possible negative attributes from large agglomerations to the mechanical performances, as the aggregation of nanoparticles might act as a weak spot compared with the uniform dispersed nanoparticles. The efficiency of silane functionalisation on the uniform dispersion of BaTiO 3 nanoparticles was demonstrated.