Poling field versus piezoelectric property for [001] c oriented 91%Pb(Zn 1/3 Nb 2/3 )O 3 9%PbTiO 3 single crystals

Similar documents
Orientation dependence of electromechanical properties of relaxor based ferroelectric single crystals

Intermediate ferroelectric orthorhombic and monoclinic M B phases in [110] electric-field-cooled Pb Mg 1/3 Nb 2/3 O 3 30%PbTiO 3 crystals

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

Temperature-dependent phase transitions in Pb(Zn1/3Nb2/3)0.93Ti0.07O3 crystal

Micro-Brilouin scattering study of field cooling effects on ferroelectric relaxor PZN-9%PT single crystals

Phase diagram of the ferroelectric relaxor (1-x)PbMg1/3Nb2/3O3-xPbTiO3 Noheda, B.; Cox, D.E.; Shirane, G.; Gao, J.; Ye, Z.-G.

Electric Field- and Temperature-Induced Phase Transitions in High-Strain Relaxor- Based Ferroelectric Pb(Mg1 /3Nb2/3)1 - xtixo3 Single Crystals

I. Introduction. Feiming Bai, Naigang Wang, Jiefang Li, and D. Viehland 1, P.M. Gehring 2, Guangyong Xu and G. Shirane 3. (Received Feb.

Physics Department, Brookhaven National Laboratory, Upton, NY 11973, USA

The Monoclinic Phase in PZT: New Light on Morphotropic Phase Boundaries

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

Additional information on J. Appl. Phys. J. Peräntie, J. Hagberg, A. Uusimäki, J. Tian, and P. Han

Dynamic origin of the morphotropic phase boundary by Hu Cao et al.

Piezoelectric materials for MEMS applications Hiroshi Funakubo Tokyo Institute of Technology

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

Supporting Information

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

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

Domain-size dependence of piezoelectric properties of ferroelectrics

Supplementary Information

PLEASE SCROLL DOWN FOR ARTICLE

Mechanical confinement: An effective way of tuning properties of piezoelectric crystals

T d T C. Rhombohedral Tetragonal Cubic (%) 0.1 (222) Δa/a 292K 0.0 (022) (002) Temperature (K)

Brillouin scattering and molecular dynamics study of the elastic properties of Pb(Mg 1/3 Nb 2/3 )O 3

Influence of Ceramic Particle Sizes on Electrical Properties of Lead Zirconate Titanate (PZT)/Nylon57 Composites

Adaptive ferroelectric states in systems with low domain wall energy: Tetragonal microdomains

Relaxor characteristics of ferroelectric BaZr 0.2 Ti 0.8 O 3 ceramics

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

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

Waste heat energy harvesting using the

Dielectric Properties and Lattice Distortion in Rhombohedral Phase Region and Phase Coexistence Region of PZT Ceramics

Ferroelectric materials contain one or more polar axes along which a spontaneous

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

First-principles study of BiScO 3 1Àx - PbTiO 3 x piezoelectric alloys

Effect of domains configuration on crystal structure in ferroelectric ceramics as revealed by XRD and dielectric spectrum

Structural Phase Transition and Dielectric Relaxation in Pb(Zn 1/3 Nb 2/3 )O 3 Single Crystals

Electromechanical-induced antiferroelectric ferroelectric phase transition in PbLa(Zr,Sn,Ti)O 3 ceramic

Singularities in Landau-Devonshire potentials for Ferroelectrics Phase-Transitions

Effect of Zr/Ti Ratio Content on Some Physical Properties of Low Temperature Sintering PZT PZN PMnN Ceramics

Probing local polar structures in PZN-xPT and PMN-xPT relaxor ferroelectrics with neutron and x- ray scattering

Microstructure, phase transition, and electrical properties of K 0.5 Na x Li x Nb 1 y Ta y O 3 lead-free piezoelectric ceramics

Effect of substrate-induced strains on the spontaneous polarization of epitaxial BiFeO 3 thin films

Phase transition behavior and high piezoelectric properties in lead-free BaTiO 3 CaTiO 3 BaHfO 3 ceramics

Physical and Dielectric Properties of Silver Lithium Niobate Mixed Ceramic System

Dielectric Properties of Two-Stage Sintered PMN-PT Ceramics Prepared by Corundum Route

CMP Seminar April 3 rd, 2008

Seung-Ho PARK, Seyit URAL, Cheol-Woo AHN 1, Sahn NAHM 1 and Kenji UCHINO

Depolarization of a piezoelectric film under an alternating current field

Theerachai Bongkarn 1, Naratip Vittayakorn 2 and Gobwute Rujijanagul 3. Phitsanulok 65000, Thailand. Chiang Mai Thailand

Modern Piezoelectrics

Advanced Piezoelectric Single Crystal Based Actuators

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

crystals ISSN

Pressure as a Probe of the Physics of ABO 3 Relaxor Ferroelectrics

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

High tunable dielectric response of Pb 0.87 Ba 0.1 La 0.02 (Zr 0.6 Sn 0.33 Ti 0.07 ) O 3 thin film

Single-phase driven ultrasonic motor using two orthogonal bending modes of sandwiching. piezo-ceramic plates

INVESTIGATIONS OF POLARIZATION SWITCHING OVER BROAD TIME AND FIELD DOMAINS IN VARIOUS FERROELECTRICS

Aging effect evolution during ferroelectricferroelectric phase transition: A mechanism study

Generalized continuum theory for ferroelectric thin films

Advanced. piezoelectric. materials. Science and technology. Edited by. Kenji Uchino WOODHEAD PUBLISHING. Oxford Cambridge Philadelphia New Delhi

Ferroelectricity. Phase transition. Material properties

arxiv:cond-mat/ v1 [cond-mat.mtrl-sci] 14 Feb 2001

PIEZOELECTRIC ceramics for advanced actuator applications often

We are IntechOpen, the world s leading publisher of Open Access books Built by scientists, for scientists. International authors and editors

Experimental search for high Curie temperature piezoelectric ceramics with combinatorial approaches

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

Temperature dependence of piezoelectric properties of high-tc Bi(Mg1/2Ti1/ 2)O3 PbTiO3

BIASED PZT MATERIALS FOR ACOUSTIC TRANSDUCERS

Growth and Characterization of Lead Zirconate- Titanate (PbZr 1-x Ti x O 3 )-Based Novel Piezo- /Ferroelectric Single Crystals

Electric field-induced phase transitions in (111)-, (110)-, and (100)-oriented Pb(Mg1 3Nb2 3)O3 single crystals

Piezoelectric Response from Rotating Polarization

G. Ravichandran Aeronautics & Mechanical Engineering Graduate Aeronautical Laboratories California Institute of Technology

Electromechanical stability of a Mooney Rivlin-type dielectric elastomer with nonlinear variable permittivity

Monte Carlo Simulation of Ferroelectric Domain Structure: Electrostatic and Elastic Strain Energy Contributions

Phase diagram and piezoelectric response of (Ba 1 x Ca x )(Zr 0.1 Ti 0.9 )O 3 solid solution

Coexistence of monoclinic and tetragonal phases in PMN-PT single. crystal

Supporting information

First principle based phase stability in PMN-xPT. M. Sepliarsky (1) and R. E. Cohen (2)

Review on high temperature piezoelectric ceramics and actuators based on BiScO 3 PbTiO 3 solid solutions

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

Monoclinic phases arising across thermal inter-ferroelectric phase transitions

Effect of grain size on the electrical properties of Ba,Ca Zr,Ti O 3 relaxor ferroelectric ceramics

Epitaxial piezoelectric heterostructures for ultrasound micro-transducers

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

Micromechanical modeling and simulation of piezoceramic materials

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

Impedance Analysis and Low-Frequency Dispersion Behavior of Bi 4 Ti 3 O 12 Glass

Dielectric and micromechanical studies of barium titanate substituted (1-y)Pb (Zn 1/3 Nb 2/3 )O 3 -ypt ferroelectric ceramics

In situ Transmission Electron Microscopy Study on the Phase Transitionsin Lead-Free (1 x)(bi1/ 2Na1/2)TiO3 xbatio3 Ceramics

Constructing Landau-Ginzburg-Devonshire Type Models for Ferroelectric Systems Based on Symmetry

OPTIMAL DESIGN OF FERROELECTRIC CERAMICS MICROSTRUCTURE

Dielectric behaviors of Pb Fe 2/3 W 1/3 -PbTiO 3 relaxors: Models comparison and numerical calculations

Barrier Layer; PTCR; Diffuse Ferroelectric Transition.

lead-free perovskite piezoelectric ceramics Cheuk W. Tai * and Y. Lereah Department of Physical Electronics, School of Electrical Engineering,

2 ( º ) Intensity (a.u.) Supplementary Figure 1. Crystal structure for composition Bi0.75Pb0.25Fe0.7Mn0.05Ti0.25O3. Highresolution

Characteristics of Lead Free Ferroelectric Thin Films Consisted of (Na 0.5 Bi 0.5 )TiO 3 and Bi 4 Ti 3 O 12

ABSTRACT. strontium calcium titanate-barium zirconate titanate xbsct-(1-x)bzt

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

Effect of interfacial dislocations on ferroelectric phase stability and domain morphology in a thin film a phase-field model

Transcription:

J Mater Sci (2011) 46:1839 1843 DOI 10.1007/s10853-010-5009-z Poling field versus piezoelectric property for [001] c oriented 91%Pb(Zn 1/3 Nb 2/3 )O 3 9%PbTiO 3 single crystals Yang Xiang Rui Zhang Wenwu Cao Received: 9 August 2010 / Accepted: 14 October 2010 / Published online: 30 October 2010 Ó Springer Science+Business Media, LLC 2010 Abstract Electromechanical property measurements and microstructure observations using optical microscopy were performed on a [001] c oriented k 33 resonator made of 91%Pb(Zn 1/3 Nb 2/3 )O 3 9%PbTiO 3 single crystal, which was polarized under different electric fields. At room temperature, when the poling field is 1100 V/mm, the electromechanical coupling factor k 33 is 0.90 and piezoelectric coefficient d 33 is 1665 pc/n. Such superior electromechanical properties could be attributed to the formation of monoclinic multi-domain structure, which transforms to tetragonal phase at 46 C. While at a higher poling field of 1200 V/mm, the crystal becomes singledomain tetragonal state and its k 33 and d 33 are only about 0.69 and 850 pc/n, respectively. The critical poling field to transform the monoclinic phase to the tetragonal phase is found to be *1120 V/mm. Introduction As one of the most promising next-generation piezoelectric materials for high-sensitivity small size sensors, transducers and actuators, relaxor-based ferroelectric single crystals (1-x)Pb(Zn 1/3 Nb 2/3 )O 3 xpbtio 3 (PZN xpt) have been under extensive investigations for their extraordinary high piezoelectric behaviors [1 3]. Due to special structural stability, the crystal with composition near the Y. Xiang R. Zhang (&) W. Cao Department of Physics, Harbin Institute of Technology, Harbin 150080, Heilongjiang, China e-mail: ruizhang_ccmst@hit.edu.cn W. Cao Materials Research Institute, The Pennsylvania State University, University Park, PA 16802, USA morphotropic phase boundary (MPB, x & 9%) exhibit unique physical properties. There are numerous reports on unpolarized PZN 9%PT single crystals. Uesu et al. [4] observed both orthorhombic and tetragonal phases in unpoled PZN 9%PT single crystal at room temperature. Kiat et al. [5] found monoclinic and tetragonal phases in unpoled PZN 9%PT powder after being annealed at 800 C for 1 h. There were also some investigations about the coexistence of rhombohedral and tetragonal phases in PZN 9%PT single crystals [6 8]. However, up to date, there were little reports on poling influence to the properties of PZN 9%PT single crystals. In general, the macroscopic properties of domain engineered (1-x)Pb(Zn 1/3 Nb 2/3 )O 3 xpbtio 3 single crystals are strongly related to their domain structures. The aim of the present work is to experimentally establish a relationship between the macroscopic properties and domain structures of [001] c poled PZN 9%PT single crystals. Experimental procedure PZN 9%PT single crystals used in this work were supplied by the Microfine Materials Technologies P/L (Singapore). The sample was oriented, cut and polished into a k 33 resonator with the dimensions of 1 mm//[100] L 9 1 mm// [010] W 9 6 mm//[001] T. Thin layer gold electrodes were sputtered on both ends (1 9 1 mm). The bar sample was poled progressively in silicone oil at room temperature from 100 to 1400 V/mm. Before each repoling at different poling field, the sample was heat treated at 400 C for 3 h to release the internal stresses and completely depole the sample. After each poling step, the dielectric permittivity e T 33 at 1 khz, the electromechanical coupling factors k 33 and elastic compliance coefficient s E 33 were determined by

1840 J Mater Sci (2011) 46:1839 1843 the resonance and antiresonance frequencies measured by an impedance analyzer (Agilent 4294A). The piezoelectric constant d 33 was measured by a Berlincourt-type d 33 meter. Results and discussions The measured properties of [001] c oriented PZN 9%PT single crystal as a function of the poling field are given in Fig. 1. If the poling field is less than 600 V/mm, the corresponding piezoelectric coefficient d 33 is less than 400 pc/n. When the poling field is increased to 700 V/mm, d 33 increases to 1380 pc/n and the coefficients s E 33, k 33, and e T 33 are 133 9 10-12 m 2 /N, 0.87 and 2090e 0, respectively. When the poling field is further increased to 1110 V/mm, d 33 s E 33, e T 33, and k 33 reach 1670 pc/n, 136 9 10-12 m 2 /N, 2807e 0 and 0.90, respectively. But when the poling field is increased to 1120 V/mm, all the electromechanical properties start to decrease abruptly. When the field level E reached 1200 V/mm, the elastic compliance coefficient s 33 becomes 120 9 10-12 m 2 /N, and the dielectric permittivity e 33 (1426e 0 ), electromechanical coupling factors k 33 T (0.69), and piezoelectric constant d 33 (849 pc/n) of the sample become much smaller than the corresponding Fig. 1 Dielectric permittivity e T 33, elastic compliance s E 33, electromechanical coupling factor k 33, and piezoelectric coefficient d 33 of PZN 9%PT single crystal as functions of poling field values of the sample poled under a field of 1100 V/mm (e T 33 = 2807e 0, k 33 = 0.90, and d 33 = 1665 pc/n, respectively). Therefore, the critical poling field for this sample is 1120 V/mm. In order to explain this phenomenon, we have carried out a detailed study on the electromechanical property and domain structures of PZN 9%PT single crystal poled under 1100 and 1200 V/mm. Upon heating the sample from room temperature with a heating/cooling stage (Linkam TEMS600), the dielectric permittivity e T 33 at 1 khz, electromechanical coupling factors k 33 and elastic compliance coefficient s E 33 were measured by the resonance technique using an impedance analyzer and the piezoelectric constant d 33 was calculated from the following equation: qffiffiffiffiffiffiffiffiffiffiffi d 33 ¼ k 33 : ð1þ e T 33 se 33 Ferroelectric domain configurations were observed by polarized light microscopy (ZEISS AXIOSKOP40). To be consistent for all measurements, the same k 33 resonator was used for all studies, including measurements of dielectric, piezoelectric, elastic compliance coefficients, and optical microscopy studies of domain structures. Figure 2 shows the temperature dependence of e T 33, s E 33, k 33, and d 33 for the bar-shaped sample poled under a field of 1100 V/mm. During heating, all electromechanical properties exhibit two peaks at about 46 and 177 C. The lower temperature abnormality at 46 C corresponds to the phase transition from monoclinic to the tetragonal phase, while the peak at 177 C indicates tetragonal to cubic paraelectric phase transition (T c ). The maximum values of e T 33 (50614e 0 ) and d 33 (4209 pc/n) appear at T c * 177 C. Monoclinic phase M B -type and M A -type with space group C m, and M C -type with space group P m were observed in these type of piezo crystals [9 11]. For the case of PZN 9%PT poled along [001] c, the M B phase can be excluded, which was observed in [110] c oriented sample under electric bias [5]. Domain structure observation were performed on the PZN 9%PT single crystal after being poled under a field of 1100 V/mm at room temperature. We found that there were no complete extinctions observed on the three pairs of surfaces in any direction under crossed polarizers due to the overlapping of different domains and the reflection of light from dense domain walls [12]. When the polarizer is along \001[, the sample becomes darker than when the crossed polarizers are in any other directions. This indicates that the domain structure belongs to the M C -type monoclinic phase [13]. Figure 3 shows domain configurations at room temperature on (100), (010), and (001) surfaces under crossed polarizers. The domain configurations on (100) shown in Fig. 3a and on (010) shown in Fig. 3b are similar, exhibit mutually perpendicular dark stripes. When focusing on the

J Mater Sci (2011) 46:1839 1843 1841 Fig. 2 Temperature dependence of dielectric permittivity e T 33, elastic compliance s E 33, electromechanical coupling factor k 33, and piezoelectric coefficient d 33 of the sample poled under a field of 1100 V/mm same area but different depth along [100] c or [010] c, the dark stripes will shift but retaining their directions along \011[ or \101[, which indicates that the domain walls are not perfectly perpendicular to the surfaces but are slightly inclined. Therefore, the different incline angles make the domain walls appear to have different thickness [12]. Crossing domain structure can be observed on the (001) surface, as shown in Fig. 3c. When the poling field is along [001] c, four different polarization orientations remained for the M C phase, which are denoted by PI, PII, PIII, and PIV as shown in Fig. 4a. According to the polarization orientation in the M C phase provided by Devonshire theory [14] as well as domain configurations observed in PZN 9%PT single crystal shown in Fig. 3, typical domain patterns are depicted in Fig. 4b. All domain walls in Fig. 4b are permissible by the mechanical compatibility argument in ferroelectric crystals with M C monoclinic symmetry [12]. From the observation and analysis above, it can be concluded that M C phase exists in [001] c oriented PZN 9%PT crystal when it is poled under a field level of 1100 V/mm at room temperature. It should be noted that the domain structures described here are different from that reported in Ref. [15] for the poled PZN 9%PT crystals. The authors of Ref. [15] observed broad, straight domains separated by the walls along [010] c, which may belong to monoclinic C m space group [12]. However, [001] c electric field could not induce monoclinic M A and M B phases in PZN 9%PT single crystal [16], so the crystal used in Ref. [15] might have PT content lower than 9%. For the PZN 9%PT single crystal sample poled under an electric field of 1200 V/mm, the dielectric permittivity e T 33, elastic compliance coefficient s E 33, electromechanical coupling factors k 33 and piezoelectric coefficient d 33 vs. temperature curves are shown in Fig. 5. Only one peak appears at T c * 177 C in each curve, which corresponds to the tetragonal to cubic phase transition. This means that the sample poled under an electric field of 1200 V/mm is in tetragonal phase at room temperature. At room temperature, when the microscope is focused on the (001) surface of the sample poled under the field of 1200 V/mm, the crystal appears in complete extinction under crossed polarizers. By focusing the microscope on (100) and (010) surfaces of the k 33 vibrator, complete extinction only occurs when the direction of one of the polarizers is along [001] c. These phenomena indicate that the crystal poled under the field of 1200 V/mm is in singledomain tetragonal state with the polarization vector along [001] c at room temperature. Fig. 3 Domain structures at room temperature on a (100), b (010), and c (001) surfaces of a k 33 resonator made of PZN 9%PT single crystal

1842 J Mater Sci (2011) 46:1839 1843 Fig. 4 a Four different polarization orientations in [001] c poled PZN 9%PT single crystal. b Illustration of domain patterns formed in a [001] c poled PZN 9%PT single crystal differently polarized domains will be completely switched to the field direction [001] c and the crystal becomes singledomain tetragonal phase. This field induced phase transition leads to the drastic degradation of the electromechanical properties. Summary and conclusions Fig. 5 Temperature dependence of dielectric permittivity e T 33, elastic compliance s E 33 electromechanical coupling factor k 33, and piezoelectric coefficient d 33 of the sample poled under a field of 1200 V/mm Moreover, it is clear that 1110 V/mm is the optimum poling field for the PZN 9%PT single crystal, which produces monoclinic multi-domain structures (Fig. 3) and gives the best electromechanical properties. If the field increases further, some monoclinic domains will be switched to the tetragonal phase, resulting in the degradation of crystal performance. This can be verified by the decrease of electromechanical property values of the crystal when the poling field is increased beyond 1120 V/mm, as shown in Fig. 1. When the field is larger than 1150 V/mm, all four In summary, electromechanical property measurements and optical microscopy observations of domain microstructures were performed on a [001] c oriented k 33 resonator made of 91%Pb(Zn 1/3 Nb 2/3 )O 3 9%PbTiO 3 single crystal under different poling electric fields, and as a function of temperature. The optimum poling field for PZN 9%PT single crystal is determined to be 1110 V/mm, which produced M C -type monoclinic multi-domain structures, the electromechanical coupling factor k 33 reaches 0.90 and the piezoelectric coefficient d 33 is 1670 pc/n. Upon heating, the M C phase transformed into tetragonal phase at about 46 C, then to cubic phase at T c * 177 C. At room temperature, the critical poling field value was found to be 1120 V/mm, beyond which some domains will be switched to tetragonal phase so that the macroscopic electromechanical properties will degrade. When the poling field is beyond 1150 V/mm, all domains will be completely switched to the field direction [001] c so that the crystal becomes single-domain tetragonal phase, resulting a drastic decrease of the piezoelectric coefficient d 33 (*850 pc/n) and dielectric permittivity e T 33 (*1400e 0 ) at room temperature. Upon heating, these single-domain resonators underwent only one phase transition from tetragonal to cubic phase at T c * 177 C. Acknowledgements This research was supported in part by the NSFC under Grant no. 50602009, 50972034, Program of the Ministry of Education of China for New Century Excellent Talents in University under Grant no. NCET-06-0345, Postdoctoral Science Research Developmental Foundation of Heilongjiang Province under

J Mater Sci (2011) 46:1839 1843 1843 Grant no. LBH-Q06068, SRF for ROCS, SEM, and NIH under Grant no. P41-EB21820. References 1. Park SE, Shrout TR (1997) J Appl Phys 82:1804 2. Kuwara J, Uchino K, Nomura S (1981) Ferroelectrics 37:579 3. Kuwata J, Uchino K, Nomura S (1982) Jpn J Appl Phys 21:1298 4. Uesu Y, Matsuda M, Yamada Y, Fujishiro K, Cox DE, Noheda B, Shirane G (2002) J Phys Soc Jpn 71:960 5. Kiat M, Uese Y, Dkhil B, Matsuda M, Malibert C, Calvarin G (2002) Phys Rev B 65:064106 6. Uesu Y, Yamada Y, Fujishiro K, Tazawa H, Enokido S, Kiat JM, Dikhil B (1998) Ferroelectrics 217:319 7. Zhang L, Dong M, Ye ZG (2000) Mater Sci Eng B 78:96 8. Singh G, Bhaumik I, Ganesamoorthy S, Karnal AK, Tiwari VS (2007) Cryst Ref Technol 42:378 9. Cao H, Bai FM, Wang NG, Li JF, Viehland D, Xu GY, Shirane G (2005) Phys Rev B 72:064104 10. Ye ZG, Noheda B, Dong M, Cox D, Shirane G (2001) Phys Rev B 64:184114 11. Noheda B, Cox DE, Shirane G, Gao J, Ye ZG (2002) Phys Rev B 66:054104 12. Bokov AA, Ye ZG (2004) J Appl Phys 95:6347 13. Bao P, Yan F, Liu XM, Zhu JS, Shen HM, Wang YN, Luo HS (2006) Appl Phys Lett 88:092905 14. Vanderbilt D, Cohen MH (2001) Phys Rev B 63:094108 15. Renault AE, Dammak H, Calvarin G, Thi MP, Gaucher P (2002) Jpn J Appl Phys 41:3846 16. Noheda B, Zhong Z, Cox DE, Shirane G, Park SE, Rehrig P (2002) Phys Rev B 65:224101