White Rose Research Online URL for this paper: Version: Accepted Version

Similar documents
This is a repository copy of Conceptualizing a circular framework of supply chain resource sustainability.

This is a repository copy of An in vivo platform for identifying inhibitors of protein aggregation.

White Rose Research Online URL for this paper: Version: Accepted Version

This is a repository copy of The Uptake of HO ₂ on Meteoric Smoke Analogues.

This is a repository copy of Aggregation of growing crystals in suspension: III. Accounting for adhesion and repulsion.

Piezo materials. Actuators Sensors Generators Transducers. Piezoelectric materials may be used to produce e.g.: Piezo materials Ver1404

This is a repository copy of Observation of electrolytic capacitor ageing behaviour for the purpose of prognostics.

This is a repository copy of Myths about HIV and AIDS among serodiscordant couples in Malawi.

Non-linear Properties of Piezoelectric Ceramics

This is a repository copy of Microbiology: Mind the gaps in cellular evolution.

White Rose Research Online URL for this paper: Version: Accepted Version

1990. Temperature dependence of soft-doped / hard-doped PZT material properties under large signal excitation and impact on the design choice

Ferroelectric Ceramic Technology for Sensors. Jim McIntosh Ceramic Operations Manager PCB Piezotronics

White Rose Research Online URL for this paper: Version: Accepted Version

An equivalent dipole analysis of PZT ceramics and lead-free piezoelectric single crystals

This is a repository copy of Optimised PEI impregnation of activated carbons - Enhancement of CO2 capture under post-combustion conditions.

Addition 1. Shear Stack Piezoelectric Elements and Shear Effect Basics

This is a repository copy of All-acrylic film-forming colloidal polymer/silica nanocomposite particles prepared by aqueous emulsion polymerization.

Depolarization of a piezoelectric film under an alternating current field

This is a repository copy of The Use of Stream Power as an Indicator of Channel Sensitivity to Erosion and Deposition Processes.

Online publication date: 23 June 2010 PLEASE SCROLL DOWN FOR ARTICLE

This is a repository copy of Effective thermal conductivity of porous solder layers.

This is a repository copy of Improved analytical model for predicting the magnetic field distribution in brushless permanent-magnet machines.

Micromechanical modeling and simulation of piezoceramic materials

Mechanical stability of polar defects in ferroelectric perovskites

THERMAL DEGRADATION AND AGING OF HIGH TEMPERATURE PIEZOELECTRIC CERAMICS

This is a repository copy of Defining a 3-dimensional trishear parameter space to understand the temporal evolution of fault propagation folds.

White Rose Research Online URL for this paper: Version: Accepted Version

Supplementary Information

This is a repository copy of Measure for Measure. Semi-cylindrical illuminance: a semi-conceived measure?.

PIEZOELECTRIC ceramics for advanced actuator applications often

This is a repository copy of Bringing heaven down to earth : the purpose and place of religion in UK food aid.

I. INTRODUCTION II. SAMPLE PREPARATION JOURNAL OF APPLIED PHYSICS VOLUME 92, NUMBER 5 1 SEPTEMBER

The electric field induced strain behavior of single. PZT piezoelectric ceramic fiber

Determination of the lead titanate zirconate phase diagram by the measurements of the internal friction and Young s modulus

Title. Author(s)H. H. PAN; C.K. CHIANG; R.H. YANG; Y.H. WU; C.S. CHA. Issue Date Doc URL. Type. Note. File Information CONTAINING SLAG

Laser Interferometric Displacement Measurements of Multi-Layer Actuators and PZT Ceramics

Domain and Phase Change Contributions to Response in High Strain Piezoelectric Actuators

White Rose Research Online URL for this paper: Version: Supplemental Material

BIASED PZT MATERIALS FOR ACOUSTIC TRANSDUCERS

Double-peak switching current in soft ferroelectric lead zirconate titanate Kamel, T.M.; de With, G.

DEVELOPMENT OF HIGH-STRAIN LOW-HYSTERESIS ACTUATORS USING ELECTROSTRICTIVE LEAD MAGNESIUM NIOBATE (PMN)

Orientation dependence of electromechanical properties of relaxor based ferroelectric single crystals

This is a repository copy of Analysis and design optimization of an improved axially magnetized tubular permanent-magnet machine.

Effect of Crystal Orientation on Piezoelectric Response of Single-Crystal Piezoelectric Layers

This is a repository copy of An application of Kozeny Carman flow resistivity model to predict the acoustical properties of polyester fibre.

Temperature Dependence of the Dielectric, Elastic and Piezoelectric Material Constants of Lead Zirconate Titanate Ceramics

This is a repository copy of Experimentation and modelling of near field explosions.

Ferroelectric ceramics for high-power applications,

This is a repository copy of Evaluation of output frequency responses of nonlinear systems under multiple inputs.

Modern Piezoelectrics

This is a repository copy of Incorporating intraspecific trait variation into functional diversity: Impacts of selective logging on birds in Borneo.

Domain-size dependence of piezoelectric properties of ferroelectrics

Dae Won Ji and Sang-Joo Kim

Piezoelectric materials for MEMS applications Hiroshi Funakubo Tokyo Institute of Technology

This is a repository copy of Study of wind profile prediction with a combination of signal processing and computational fluid dynamics.

Effects of niobium doping on the piezoelectric properties of sol gel-derived lead zirconate titanate films

crystals ISSN

This is a repository copy of A three-dimensional numerical model of borehole heat exchanger heat transfer and fluid flow.

An evaluation of switching criteria for ferroelectrics under stress and electric field

Piezoelectric Composites as Bender Actuators

Stress-Induced Phase Transition in Pb(Zr 1/2 Ti 1/2 )O 3

White Rose Research Online URL for this paper:

ACOUSTIC EMISSION MEASUREMENTS ON PIEZOELECTRIC/ FERROELECTRIC MATERIALS

PIEZOELECTRIC TECHNOLOGY PRIMER

Lecture 7: Extrinsic semiconductors - Fermi level

Transduction Based on Changes in the Energy Stored in an Electrical Field

This is a repository copy of The effect of the dynamical state of clusters on gas expulsion and infant mortality.

Electric field dependent sound velocity change in Ba 1 x Ca x TiO 3 ferroelectric perovskites

MODELING AND SIMULATION OF PIEZOCERAMIC MATERIALS USING MICROMECHANICAL APPROACH

PHYSICAL REVIEW B 81,

Chapter 6 ELECTRICAL CONDUCTIVITY ANALYSIS

This is a repository copy of TACTIC : The TRIUMF Annular Chamber for Tracking and Identification of Charged particles.

This is a repository copy of Improving the associative rule chaining architecture.

Analysis of the state of poling of lead zirconate titanate (PZT) particles in a Zn-ionomer composite

Sensing behavior of ferromagnetic shape memory Ni-Mn-Ga

Piezoelectric nonlinearity and frequency dispersion of the direct piezoelectric response. of BiFeO 3 ceramics. Switzerland

Classification of Dielectrics & Applications

Piezo Theory: Chapter 1 - Physics & Design

Modifying the Electrical Properties of Ba 0 85 Ca 0 15 Zr 0 1 Ti 0 9 O 3 Ceramics by the Nanocrystals-Induced Method

ANALYSIS AND NUMERICAL MODELLING OF CERAMIC PIEZOELECTRIC BEAM BEHAVIOR UNDER THE EFFECT OF EXTERNAL SOLICITATIONS

Domain switching and electromechanical properties of pulse poled Pb Zn 1Õ3 Nb 2Õ3 O 3 PbTiO 3 crystals

This is a repository copy of Modelling Outliers and Missing Values in traffic Count Data Using the ARIMA Model..

Dielectric Materials: Properties and Applications

STABILITY OF MECHANICAL AND DIELECTRIC PARAMETERS IN PZT BASED CERAMICS JAN ILCZUK, JUSTYNA BLUSZCZ, RADOSŁAW ZACHARIASZ

SURFACE BARKHAUSEN NOISE INVESTIGATIONS OF STRESS AND LEAKAGE FLUX

Striking similarity of ferroelectric aging effect in tetragonal, orthorhombic and rhombohedral crystal structures

Aging effect evolution during ferroelectricferroelectric phase transition: A mechanism study

Electro-shape-memory effect in Mn-doped BaTiO 3 single crystals and in situ observation of the reversible domain switching

Effect of Domain Wall Motion and Phase Transformations on Nonlinear Hysteretic Constitutive Behavior in Ferroelectric Materials

Hardening-softening transition in Fe-doped Pb(Zr,Ti)O 3 ceramics and evolution of the third harmonic of the polarization response

Bridging to the Continuum Scale for Ferroelectric Applications

Transducer Control Algorithms

Final Project Report. Constitutive Behavior of Relaxor Single Crystals. Submitted to W. Smith, ONR. August 21, Christopher S.

First-principles calculations of piezoelectricity and polarization rotation in Pb Zr 0.5 Ti 0.5 O 3

White Rose Research Online URL for this paper:

Low Voltage Ionic Wind Generation using Piezoelectric Transformers

Ferroelectricity. Phase transition. Material properties. 4/12/2011 Physics 403 Spring

SUPPLEMENTARY INFORMATION

APPENDIX A Landau Free-Energy Coefficients

Transcription:

This is a repository copy of Temperature Dependence of Domain Contributions as a Function of Ageing in Soft and Hard Lead Zirconate Titanate Piezoelectric Ceramics. White Rose Research Online URL for this paper: http://eprints.whiterose.ac.uk/117448/ Version: Accepted Version Article: Qaisar, SA, Comyn, TP and Bell, AJ (2017) Temperature Dependence of Domain Contributions as a Function of Ageing in Soft and Hard Lead Zirconate Titanate Piezoelectric Ceramics. IEEE Transactions on Ultrasonics, Ferroelectrics, and Frequency Control, 64 (6). pp. 1023-1028. ISSN 0885-3010 https://doi.org/10.1109/tuffc.2017.2686490 2017, IEEE. This is an author produced version of a paper published in IEEE Transactions on Ultrasonics, Ferroelectrics and Frequency Control. Personal use of this material is permitted. Permission from IEEE must be obtained for all other users, including reprinting/ republishing this material for advertising or promotional purposes, creating new collective works for resale or redistribution to servers or lists, or reuse of any copyrighted components of this work in other works. Uploaded in accordance with the publisher s self-archiving policy. Reuse Unless indicated otherwise, fulltext items are protected by copyright with all rights reserved. The copyright exception in section 29 of the Copyright, Designs and Patents Act 1988 allows the making of a single copy solely for the purpose of non-commercial research or private study within the limits of fair dealing. The publisher or other rights-holder may allow further reproduction and re-use of this version - refer to the White Rose Research Online record for this item. Where records identify the publisher as the copyright holder, users can verify any specific terms of use on the publisher s website. Takedown If you consider content in White Rose Research Online to be in breach of UK law, please notify us by emailing eprints@whiterose.ac.uk including the URL of the record and the reason for the withdrawal request. eprints@whiterose.ac.uk https://eprints.whiterose.ac.uk/

1 Temperature Dependence of Domain Contributions as a Function of Ageing in Soft and Hard Lead Zirconate Titanate Piezoelectric Ceramics S. Adam Qaisar, Tim P. Comyn, and Andrew J. Bell, Senior Member, IEEE Abstract The stress and temperature dependence of three different PZT ceramics have been measured at 48 hours and 1 month after poling. The data were fitted to a modified Rayleigh law in which a negative quadratic stress term accounts for a saturation of the domain wall contributions at higher stress. The relative changes to the fitting parameters can be explained by the materials possessing differences in the concentration of orientable defect dipoles, but with similar distributions of relaxation time constants. P Index Terms piezoelectricity I. INTRODUCTION IEZOELECTRIC materials are important contributors to electronic technology in actuator, sensor and signal processing applications. For the most widely used materials such as lead zirconate titanate ceramics, Pb(Zr xti 1-x)O 3 or PZT, it is well established that the presence of ferroelectric domains contributes significantly to the piezoelectric effect [1 6]. Specifically, applied electric field or stress can induce local reorientation of the crystal lattice, observed as the displacement of domains walls. Non-180 domain wall motion augments the field-induced strain or stress induced charges and therefore increases the value of the most widely used measure of piezoelectric activity, the charge coefficient, d ij. Piezoelectric materials are generally categorized as either soft or hard. Soft materials are engineered to maximize the charge coefficient at the expense of increasing the dielectric and mechanical losses, whilst hard materials are designed to have low losses and a high electromechanical Q factor, at the expense of a reduced charge coefficient. The observed differences between hard and soft piezoelectric materials are strongly correlated with the defect chemistry of the materials which is thought to control domain wall mobility [7]. In soft, donor doped PZT, there is evidence to suggest that up to 50% of the value of the d coefficient can be accounted for by non-180 domain wall motion. It is also widely accepted that the observed non-linearity of the piezoelectric effect with respect to applied field and stress is principally due to the influence of domain walls [7]. However, the form of the non-linearity is dependent on a number of issues. Whilst many materials, including soft PZTs, exhibit a linear increase in the d coefficient with increasing field or stress, known as the Rayleigh law, others, notably hard materials, exhibit higher order stress/field dependence, but with significant evolution over time. The non-linear behavior of piezoelectric materials facilitates an understanding of how the material will perform in real world scenarios where the material may be subject to a variety of static or dynamic mechanical loads and fields over time. As many applications call for wide operating temperature range, it is also important to understand how non-linear behavior is modified by temperature. Hence the work presented here investigates the non-linearity of the direct piezoelectric effect of 3 types of commercial PZT ceramics as a function of temperature and ageing. These measurements were conducted up to a maximum of 150 C for all samples in a single experimental step, and after two different ageing periods. II. EXPERIMENTAL The material selection consisted of a soft (Navy Type VI - Noliac A/S, Kvistgaard, Denmark), a hard (Navy Type I - PI Ceramic GmbH, Lederhose, Germany) and an intermediate (Navy Type II - PI Ceramic GmbH, Lederhose, Germany) PZT compositions, as shown in Table I. All materials were obtained off the shelf and cut to ensure an aspect ratio of at least 1:1 (height width) to minimize Financial support was provided by the UK Engineering and Physical Sciences Research Council, and Aero Engine Controls plc. The authors would like to thank Physik Instrumente GmbH and Noliac A/S for assistance with sample supply and preparation. At the time the work was carried out all authors were with the School of Chemical and Process Engineering, University of Leeds, Leeds, LS2 9JT, UK. Corresponding author: Andrew J. Bell, School of Chemical and Process Engineering, University of Leeds, Leeds, LS2 9JT, UK. (a.j.bell@leeds.ac.uk). S. Adam Qaisar is now with the National Nuclear Laboratory, Sellafield, Cumbria, CA20 1PG, UK (adam.qaisar@nnl.co.uk) Tim P. Comyn is now with Ionix Advanced Technologies Ltd., 3M Buckley Innovation Centre, Firth St, Huddersfield, HD1 3BD, UK (tim.comyn@ionix.at).

2 TABLE I: MANUFACTURER MATERIAL PARAMETERS. Material Type d 33 / pc N -1 33/ 0 s 33 / TPa -1 k 33 tan Q m T C / C PI PIC181 Hard 265 1200 14.2 0.66 0.003 2000 330 PI PIC255 Intermediate 400 1750 20.7 0.69 0.020 80 350 Noliac NCE55 Soft 670 5000 21.0 0.72 0.022 70 159 measured 48 hours after poling, from room temperature to 150 C in 25 C increments, using the apparatus described above, with a 10 MPa compressive pre-stress and an additional cyclic compressive stress, varying sinusoidally between zero and a maximum amplitude with a frequency of 1 Hz, as shown in Fig. 2. The measurements on each sample were conducted in one experimental run, to ensure consistency between measurements at different stresses and to minimize the influence of thermal cycling on the material behavior. The measurements were then repeated after 1 month ageing at room temperature. Fig. 1. Schematic of the test apparatus - comprising of band heater, stainless steel and zirconia stand-offs and thermal isolation. clamping effects [1], [8]. The soft and intermediate samples took the form of a cylinder (Ø 4.5mm and Ø 4 mm respectively), with the hard PZT as a rectangular prism (2 x 2 x 4 mm). For consistency, they were thermally de-poled by progressive heating to beyond the material-specific depolarization temperature. All three were then re-poled at 2.25 kv mm -1 in silicone oil for 10 minutes at 70 C. All materials were left at room temperature for 48 hours before further measurements. Direct longitudinal charge coefficient measurements were carried out as a function of stress amplitude using an in-house constructed system comprising a material test frame (3382, Instron) capable of providing a constant stress pre-load and a PZT-based actuator (P-225, Physik Instrumente UK Ltd, UK) to apply a cyclic stress at 1 Hz, measured using an in-line quartz load cell and charge amplifier (9641B and 5011B, Kistler Instruments Ltd, UK). The charge generated by the sample was measured using an identical charge amplifier connected directly to the material under test. The sample was positioned within a temperature controlled region within the test apparatus shown in Fig. 1. To mitigate non-parallel loading of the sample surfaces, the apparatus was designed with a ±5 accommodation in the orientation of the platens contacting the sample. The experimental temperature was measured using a thermocouple, connected directly to thermal acquisition hardware. All signal generation and data acquisition was synchronous to 1 ms, with data saved in a streaming time division multiplexed format. Fig. 2 shows a typical example of raw data obtained from one complete stress cycle, with a linear regression fit to calculate d 33. The direct piezoelectric coefficient of all samples was Fig. 2. Typical example of raw data collected from the experimental apparatus, showing one complete stress cycle charge response and linear regression fit (red) to calculate d 33. III. RESULTS The experimental data are presented in Fig. 3 as plots of d 33 vs stress amplitude X 0 as a function of temperature at both 48 hours and 1 month after poling. We attribute the differences between the manufacturers own values of d 33 and our measured values to the probable differences in poling regimes. After the initial measurements at 48 hours, the samples were re-measured using the Berlincourt meter at room temperature and found to have no more than a 3% change in charge coefficient compared to the value before the temperature excursion. On initial analysis, the hard composition exhibits substantially linear behavior between d 33 and X 0, that is it conforms to the Rayleigh law:, (1) where d i is the piezoelectric charge coefficient at zero stress and is the Rayleigh coefficient. However, the intermediate

3 (a) (b) (c) (d) (e) (f) Fig. 3. d 33 as a function of stress amplitude at 1 Hz and temperature for (a) Hard PZT at 48 hours; (b) Hard PZT at 1 month; (c) Intermediate PZT at 48 hours; (d) Intermediate PZT at 1 month; (e) Soft PZT at 48 hours; (f) Soft PZT at 1 month ageing. and soft compounds are initially non-linear, with higher order stress dependence, hence all the data have been analyzed according to:. (2) The coefficients d i, and are compared in Fig. 4 for all three materials as a function of temperature. It can be observed that at 48 hours the hard material is substantially linear in terms of the stress dependence of the charge coefficient, hence the value of is insignificant. There is also little temperature dependence of both d i and the first order coefficient. After 1

4 month the value of d i has increased slightly but with somewhat less temperature dependence. If it is assumed that d i represents the reversible contribution to piezoelectricity, whilst the value of indicates the size of the extrinsic, or irreversible contribution, often associated with domain wall motion, then both contributions are relatively stable as a function of time and temperature. At 48 hours the intermediate material exhibits non-linear behavior with low temperature dependence of d i but significant temperature dependence of both and. The former has a positive value and decreases with increasing temperature, whilst the latter is negative and increases with increasing temperature. However, after 1 month, there is a substantial decrease in d i and both the amplitude and temperature dependence of both and has decreased markedly. The soft material exhibits the largest magnitudes and the greatest temperature dependence of all three parameters. The initial value of charge coefficient, d i, decreases with increasing temperature, with little change between the 48 hour and 1 month values. Ageing for 1 month sees an approximate threefold reduction in the value of and a reduction in its temperature dependence, whilst the magnitude of decreases with increasing temperature, but is relatively unchanged during the ageing period. due to the ferroelastic domain wall motion. A positive value of represents increasing strain contributions (i.e. greater domain wall motion) as a function of increasing stress amplitude. We interpret, which in the current analysis is always negative, as being the extent to which this effect diminishes or saturates at di (pc N -1 ) (a) (b) IV. DISCUSSION The most widely employed interpretation of the difference in properties between hard and soft piezoelectric materials relies on the model of the ferroelastic domain wall motion contribution to the charge coefficient being pinned to a greater or lesser extent by defect dipoles formed between oxygen vacancies and acceptor dopants or impurities. Under this model, following poling, defect dipoles align themselves with the direction of local polarization. As the time constant for reorientation of such dipoles is usually much longer than the period of the signal used to measure the charge coefficient, the defect dipoles inhibit the domain wall motion. In hard materials, it is assumed that there is a high concentration of such dipoles, whilst soft materials have a lower defect dipole concentration. The ageing process in piezoelectric ceramics, by which the piezoelectric activity decreases as a function of time is therefore associated with the gradual alignment of defect dipoles with the local polarization. Whilst the Rayleigh law (Eqn. (1)) is cited by a number of authors as indicative of domain wall processes, it is clear that the majority of data here do not adhere to a simple linear dependence of charge coefficient upon stress, hence a higher order term has been introduced (Eqn. (2)). This term is introduced as the simplest means of characterizing the departure from Rayleigh linearity; it is not assigned any mechanistic significance. It is assumed that d i is a purely reversible contribution to the piezoelectric coefficient and is dominated by the lattice contributions, which, in addition to extensions of the polarization vector along the polar axis, include polarization rotation. The coefficients and represent the extrinsic contributions (c) Fig. 4. Coefficients of the second order Rayleigh fit (a) d i, (b) and (c) as a function temperature.

5 higher stresses. Whilst in soft materials, the low density of defect dipoles allows relative ease of domain wall translation, occasional deep traps can completely prevent further motion of individual walls, thereby leading to saturation at high stress. The soft and intermediate materials have significant values of, reflecting a sparse population of defects but with sufficient deep traps to diminish the extrinsic contributions at higher field. The hard material exhibits a vanishingly small value, characteristic of a material in which a high concentration of homogeneously distributed defects effectively inhibits domain wall motion at all stress levels; hence there is no diminution at higher stress. Furthermore, in the hard material, the value of does not change as a function of time or temperature. This suggests that the time constant for orientation of the defect dipoles is very short compared to the timescale of these experiments (<< 48 hours). Also, any extrinsic contribution that may be present is not thermally activated, as may be expected for domain wall motion. Indeed, it may be that the Rayleigh behavior exhibited in the hard material may not be domain related, but is due to intrinsic non-linearity in the elastic constants. This effect may be masked in softer materials by the domain contributions. In the soft material, and d i have similar negative temperature dependence. According to simple phenomenological arguments, their values both depend upon the spontaneous polarization, P s. As the Curie temperature of the soft materials is approximately 150 C less than the intermediate and hard samples, the decrease of P s as a function temperature dominates their behavior more than in the higher T C materials. The small change in d i as a function of ageing is not regarded as significant in this case. The reduction in both the amplitude and temperature dependence of after 1 month may indicate some stabilization of the orientation of the dipole population. However, stays approximately constant indicating that the sparse distribution of deep traps still results in diminution of the extrinsic contribution at higher stress. As demonstrated mainly by the values of, there appears to be a marked change in the material around 75 C. This has not been investigated further, but may be indicative of a shift in the morphotropic phase boundary as a function of temperature. The intermediate material has many similarities to the soft at 48 hours and certainly fulfils its categorization as being intermediate to the hard and soft behaviors. However, after 1 month it is almost indistinguishable from the hard material, other than having slightly higher values of both d i and. This suggests that the defect dipoles take longer to orient with the local polarization than in the hard material; i.e. after 48 hours the defects are still randomly oriented and so do not act to inhibit extrinsic contributions, whereas after 1 month they are completely oriented. Hence, we might express the differences in the three materials simply in terms of the time constant for defect dipoles to align with the local polarization ( ); the hard material has << 1.7 x 10 5 s, for the soft > 2.5 x 10 6 s and for the intermediate, 1.7 x 10 5 < < 2.5 x 10 6 s. It is widely understood that the origin of the difference between hard and soft materials Fig. 5 Conceptual plot of distribution of dipole orientation time constants for the three materials. is in concentration of acceptors and oxygen vacancies. So how is the behavior described here manifest itself as a difference in alignment time constants? We can assume that each material has a different total concentration of defect dipoles, but with a similar distribution of alignment time constants and that a threshold concentration of aligned dipoles is required to completely eliminate domain wall contributions (Fig. 5). For a large total concentration of dipoles, as in the hard material, a sufficient number of dipoles have time constants below 48 hours and therefore exhibit the hard state within that period. However, in the intermediate material the threshold concentration of aligned dipoles is not achieved until some time between 48 hours and 1 month. In the soft material, the concentration is so low that the threshold is never achieved, although some hardening is observed before 1 month. V. CONCLUSIONS The stress and temperature dependence of three different PZT ceramics have been measured at two different post-poling intervals. The data were fitted to a modified Rayleigh law in which a negative quadratic stress term accounts for a moderation in the domain wall contributions at higher stress. The relative changes to the fitting parameters can be explained by the materials possessing differences in the concentration of orientable defect dipoles, but with similar distributions of relaxation time constants. REFERENCES [1] D. Damjanovic and M. Demartin, The Rayleigh law in piezoelectric ceramics, Ferroelectrics, vol. 29, pp. 2057 2060, 1996. [2] T. P. Comyn and A. J. Bell, Study of intrinsic / extrinsic piezoelectric contributions in La-doped BiFeO 3- PbTiO 3 ceramics using the Rayleigh method, Proc. IEEE 16 th ISAF, Nara, 2007 pp 402-403. IEEE, 2007. [3] K. Kobayashi, K. Hatano, Y. Mizuno, and C. a. Randall, Rayleigh Behavior in the Lead Free Piezoelectric Li x(na 0.5K 0.5) 1-xNbO 3 Ceramic, Appl. Phys. Express, vol. 5, no. 3, p. 031501, Feb. 2012. [4] F. Li, S. Zhang, Z. Xu, X. Wei, J. Luo, and T. R. Shrout, Composition and phase dependence of the intrinsic and extrinsic piezoelectric activity of domain engineered (1-x)Pb(Mg 1/3Nb 2/3)O 3- xpbtio 3 crystals., J. Appl. Phys., vol. 108, no. 3, Aug. 2010.

[5] M. Davis, D. Damjanovic, and N. Setter, Direct piezoelectric effect in relaxor-ferroelectric single crystals, J. Appl. Phys., vol. 95, no. 10, p. 5679-5684, 2004. [6] A. Pramanick, D. Damjanovic, J. E. Daniels, J. C. Nino, and J. L. Jones, Origins of Electro-Mechanical Coupling in Polycrystalline Ferroelectrics During Subcoercive Electrical Loading, J. Am. Ceram. Soc., vol. 94, no. 2, pp. 293 309, Feb. 2011. [7] H. Jaffe, W. R. Cook and B. Jaffe, Piezoelectric Ceramics, Academic Press, London 1971. [8] A. Barzegar, D. Damjanovic, and N. Setter, The effect of boundary conditions and sample aspect ratio on apparent d 33 piezoelectric coefficient determined by direct quasistatic method., IEEE Trans. Ultrason. Ferroelectr. Freq. Control, vol. 51, no. 3, pp. 262 70, Mar. 2004. 6