# Domain wall motion in Pb(Zr\_0.20Ti\_0.80)O\_3 epitaxial thin films

Caroline Borderon, Alan Brunier, Kevin Nadaud, Raphaël Renoud, Marin

Alexe, Hartmut Wolfgang Gundel

► To cite this version:

Caroline Borderon, Alan Brunier, Kevin Nadaud, Raphaël Renoud, Marin Alexe, et al. Domain wall motion in Pb(Zr\_0.20Ti\_0.80)O\_3 epitaxial thin films. Scientific Reports, 2017, 7 (1), pp.3444. 10.1038/s41598-017-03757-y. hal-01762561

HAL Id: hal-01762561

https://hal.science/hal-01762561

Submitted on 27 Jun 2018

HAL is a multi-disciplinary open access
archive for the deposit and dissemination of sci-
entific research documents, whether they are pub-
lished or not. The documents may come from
teaching and research institutions in France or
abroad, or from public or private research centers.

L'archive ouverte pluridisciplinaire HAL, est
destinée au dépôt et à la diffusion de documents
scientifiques de niveau recherche, publiés ou non,
émanant des établissements d'enseignement et de
recherche français ou étrangers, des laboratoires
publics ou privés.

<figure><img src="image_1.png" /></figure>

Distributed under a Creative Commons Attribution 4.0 International License
OPEN

# Domain wall motion in Pb(Zr<sub>0.20</sub>Ti<sub>0.80</sub>)O<sub>3</sub> epitaxial thin films

C. Borderon<fnref n="1" />, A. E. Brunier², K. Nadaud<fnref n="1" />,³, R. Renoud<fnref n="1" />, M. Alexe² & H. W. Gundel<fnref n="1" />

Received: 1 February 2017

Accepted: 3 May 2017

Published online: 13 June 2017

Two Pb(Zr<sub>0.20</sub>Ti<sub>0.80</sub>)O<sub>3</sub> samples of different thickness and domain configuration have been studied. The c-domain sample was found to have a higher coercive field E<sub>c</sub> and higher dielectric losses than the other which presents approximately 60% of c-domains and 40% of a-domains as observed by piezo force microscopy (PFM) characterization. Hyperbolic law measurements reveal that the higher coercive field is due to domain wall pinning in deeper defects and hence a higher field E<sub>th</sub> is required for unpinning. The dissipation factors due to domain wall motion, however, are similar in both samples since the domain wall density is low and there is almost no interaction between domain walls. The higher dielectric losses in the c-domain oriented sample are a result of a greater contribution from the lattice and seem to be due to strain from the substrate, which is not relieved in a thin sample. PFM and dielectric characterization are complementary methods which provide a better understanding of the domain wall motion.

The dielectric response of domain wall motions have been already investigated for several years in polycrystalline Pb(Zr,Ti)O<sub>3</sub> (PZT) thin films<sup>1-4</sup> which present a great contribution of domain wall jump to the permittivity. In order to study those domain wall movements, the Rayleigh law has been employed at low field (E < E<sub>c</sub>/2) to maintain the domain wall density (Rayleigh region)<sup>5-8</sup>. This law is composed of two coefficients where one corresponds to the domain wall pinning/unpinning (jump) phenomena. Recently, however a hyperbolic law which is more complete has been developed<sup>9</sup>. This law permits the contribution from domain wall vibration and the lattice to be distinguished. The dielectric permittivity is then described as a function of the amplitude of the exciting electric field<sup>9</sup>:

$$ \epsilon_r = \epsilon_{rl} + \epsilon_{rw} = \epsilon_{rl} + \sqrt{\epsilon_{r-rev}^2 + (\alpha_r E_0)^2}. \quad (1) $$

This relation is valid for real and imaginary part of the permittivity. E<sub>0</sub> is the amplitude of the electric field E(t) = E<sub>0</sub> sin ωt. ε<sub>rl</sub> is the intrinsic lattice contribution which has its origin from the deformation of the unit cell<sup>5</sup>,<sup>6</sup>, and ε<sub>r-rev</sub> corresponds to domain wall vibration which is a reversible process and can be a measure of the domain wall mobility<sup>5</sup>. The sum of these two parameters corresponds to the initial permittivity (when E<sub>0</sub> = 0). The parameter α<sub>r</sub> represents the linear dependence of the permittivity on the exciting electric field. It is associated with domain wall pinning and is related to an irreversible modification of the local polarization<sup>5, 10</sup>. All parameters depend on the crystal structure but α<sub>r</sub> is especially sensitive to the presence of impurities, dopants or defects<sup>11</sup>. The parameter α<sub>r</sub> is generally high in polycrystalline PZT thin films (~1 × 10<sup>-5</sup> m/V)<sup>12, 13</sup> because the jump distance is important due to the material's columnar structure and the low degree of domain wall pinning. As a consequence, the threshold field E<sub>th</sub> necessary for unpinning is low and the Rayleigh region (linear variation of the permittivity) is clearly visible. In the present study we address the contribution of the ferroelastic 90 degree domain walls by studying epitaxial PZT thin films, one is mainly c-domain oriented and the other consists of c and a-domains. These two distinct structures allow us to correlate the polarization mechanisms to the respective domain wall motions.

## Experimental Procedure

The Pb(Zr<sub>0.20</sub>Ti<sub>0.80</sub>)O<sub>3</sub> (PZT) films were synthesized by pulsed-laser deposition (PLD) onto (001)-oriented SrTiO<sub>3</sub> (STO) substrates coated with a SrRuO<sub>3</sub> (SRO) layer which is used as a bottom electrode. The PZT composition was chosen because it has a lattice parameter of 3.935 Å which is close to that of SRO (3.928 Å)<sup>14</sup>. Two (001) oriented samples of thickness of 60 nm and 100 nm have been fabricated using similar process as in ref. 14. Briefly,

<footnote n="1">¹IETR UMR CNRS 6164, Université de Nantes, 2 rue de la Houssinière, 44322, Nantes, France. ²University of Warwick, Department of Physics, Gibbet Hill Road, Coventry, CV4 7AL, United Kingdom. ³GREMAN, CNRS UMR 7347, 16 rue Pierre et Marie Curie, 37071, Tours Cedex 2, France. Correspondence and requests for materials should be addressed to C.B. (email: caroline.borderon@univ-nantes.fr)</footnote>
<table><thead><tr><th>Sample</th><th>c-domain</th><th>a/c-domain</th></tr></thead><tbody><tr><td>Topography</td><td></td><td></td></tr><tr><td>Out of plane amplitude</td><td></td><td></td></tr><tr><td>Out of plane phase</td><td></td><td></td></tr><tr><td>In plane amplitude</td><td></td><td></td></tr><tr><td>In plane phase</td><td></td><td></td></tr></tbody></table>

<figure><img src="image_2.png" /></figure>

<figure><img src="image_3.png" /></figure>

<figure><img src="image_4.png" /></figure>

<figure><img src="image_5.png" /></figure>

<figure><img src="image_6.png" /></figure>

<figure><img src="image_7.png" /></figure>

<figure><img src="image_8.png" /></figure>

<figure><img src="image_9.png" /></figure>

<figure><img src="image_10.png" /></figure>

<figure><img src="image_11.png" /><figcaption>**Figure 1.** Piezoelectric Force Microscopy (PFM) characterization of the epitaxial PZT thin films. The scan size is 2 µm × 1 µm for the c-domain sample and 1 µm × 0.5 µm for the a/c-domain sample.</figcaption></figure>

the films were deposited using KrF excimer laser with an energy fluence of approximately 1 J/cm² on vicinal STO obtained by standard etching and annealing process of nominal 0.1° miscut crystals. Nominal 25 nm thick SRO epitaxial layers were firstly deposited at 700 °C in 150 mbar O₂ pressure followed by PZT deposited at 600 °C in 250 mbar O₂ pressure. The thickness of the PZT layer was chosen below and above the critical thickness at which misfit strain relaxes in order to obtain one fully strained PZT film showing only c-domains and one partly relaxed film showing mixed a/c-domains. Detailed microstructural investigations have been done in ref. 14.

The surface morphology of the two samples was investigated by piezoelectric force microscopy (PFM) using a Park XE-100 AFM with ac voltage of 2 V and frequency of 22.36 kHz. ATEC EFM 20 tips were used for the a/c-domain sample and Pt coated NSC14 tips for the c domain sample respectively. Stanford SR830 lock in amplifiers were used and out of plane and in plane PFM images were obtained simultaneously. The amplitude and phase characterization show the two samples (Fig. 1). The phase and amplitude images of the sample of 60 nm thickness have no variation in contrast suggesting that this sample seems to be mainly c-oriented. For the sample of 100 nm thickness, the PFM images show different regions and indicate that a-domains embedded in a c-domain matrix make up approximately 40% of the sample. In a thin layer the strain from the substrate is not relieved by the formation of ferroelastic 90 degree domain walls and the sample is mostly c-domain¹⁴. In the thicker PZT layer the total strain energy is relaxed by the formation of a polydomain structure in the epitaxial films¹⁴. Thus, we have two samples with similar crystallographic characteristics (both oriented (001)) but with different domain configurations.

Circular platinum electrodes of 160 µm radius were deposited by RF sputtering on both samples in order to realize a Metal-Insulator-Metal (MIM) capacitor. The P-E loops have been measured using ferroelectric analyzer (aixaCCt TF2000E) at 300 K and 25 kHz to avoid dielectric losses at low frequencies. The capacitance and the
<figure><img src="image_12.png" /></figure>

<figure><img src="image_13.png" /><figcaption>**Figure 2.** Ferroelectric polarizations and switching current hysteresis loops of (a) c-domain and (b) a/c-domain samples.</figcaption></figure>

dielectric losses tan δ have been determined with an Agilent 4294 A impedance meter at an AC-field amplitude of $E_0 = 5$ kV/cm. The real and imaginary parts of the permittivity have been calculated from the capacitance and the dielectric loss values. The variation of the permittivity has been measured at 10 kHz as a function the AC field amplitude (0.5 kV/cm to 100 kV/cm corresponding to an applied voltage from 5 mV to 1 V, respectively). All the measurements were performed at room temperature (20 °C).

## Results and Discussion

Ferroelectric polarization and switching current hysteresis loops are shown in Fig. 2 for the two samples. The c-domain sample has a square hysteresis loop with a high remnant polarization ($P_r \approx 95 \pm 5$ µC/cm²) which is almost equal to the saturation polarization $P_s$. This value is close to the value reported for PZT epitaxial films of the same composition¹⁵, ¹⁶. The polarization hysteresis curve of the a/c-domain sample has also a square shape, typical of hard ferroelectric. The remnant polarization $P_r$ is similar to the c-domain sample, the coercive field, however, is considerably lower ($E_c \approx 280 \pm 10$ kV/cm for the a/c-domain sample whereas $E_c \approx 780 \pm 10$ kV/cm in the case of the c-domain sample). As a consequence, the c-domain sample is less saturated than the a/c-domain sample which is easier to polarize due to the strain relieved with the presence of a-domains. The polydomain sample hence can be more easily polarized. The c-domain sample presents also a larger distribution of coercive field which indicates the presence of a higher dc conductivity of fully strained high quality PZT films. This high conductivity is due to a higher mobility related to a low density of extended defects (mainly threading dislocations) and to an influence of polarization on the Schottky barrier. Both effects lead to a higher conductance as a direct influence on the Schottky-Simmons conduction mechanism, as shown in previous work on similar films¹⁷.

The dielectric constant $\epsilon'_r$ and losses tan $\delta$ were measured as a function of frequency (Fig. 3). The permittivity of the c-domain sample is smaller ($\epsilon'_r \sim 58$ at 10 kHz) but of the same order than this of the c and a-domains ($\epsilon'_r \sim 75$ at 10 kHz). In ferroelectric materials, the permittivity of a polydomain structure is greater than that of an oriented one¹⁴, ¹⁸. This can be seen, for example, from the C-V curve where the permittivity is higher for a field near the coercive field due to the domain walls interaction¹⁸. The dielectric permittivity in the c-domain sample is the intrinsic value of dielectric permittivity, as there is no domain wall which may extrinsically enhance it ref. 19. The dielectric losses observed in the a/c-domain sample are lower than in the c-domain sample. This does not correspond to what has been reported before where the polydomain structure exhibits greater dielectric losses due to domain wall interactions¹⁸. The higher dielectric losses in our case may be due to the strain from the substrate which is not (or less) relieved in the case of the c-domain sample.
<figure><img src="image_14.png" /><figcaption>**Figure 3.** Dielectric properties (a) $\epsilon_r'$ and (b) tan$\delta$ as a function of frequency ($E_0 = 5$ kV/cm) for the two samples.</figcaption></figure>

In order to get a better understanding of the mechanisms responsible of the higher dielectric losses in the c-domain sample, the hyperbolic law has been used for decomposing the different components of the real and imaginary parts of the permittivity. The two parts of the permittivity (Fig. 4) can be perfectly fitted with the hyperbolic law (equation (1)). The evolution of the c-domain sample is less marked which shows that this sample has a lower contribution of domain wall pinning/unpinning. The calculated coefficients of the respective contributions are shown in Table 1. For the real part of the permittivity, the lattice parameter $\epsilon_r'$ is higher in the case of the a/c-domain sample due to its polydomain structure. The parameter $\epsilon_{r-rev}'$ due to reversible domain wall vibration is more than ten times larger in this sample than in the c-domain sample. According to Boser$^{20}$, this parameter is proportional to the domain wall density. In the a/c-domain sample, many domains have been observed by the PFM characterizations (Fig. 1) and the dielectric measurements and the hyperbolic law confirms the higher domain wall density. For the pinning/unpinning parameter $\alpha_r'$, the difference between the two samples is still more marked. $\alpha_r'$ in the c and a-domains sample is seventy times larger than in the c-domain sample. Again, according to Boser$^{19}$, this parameter is proportional to $\epsilon_{r-rev}'$ and can be expressed as:

$$ \alpha_r' \propto \frac{f_0 P_s}{\sqrt{NF_1}} \epsilon_{r-rev}' \quad (2) $$

where $f_0$ is a geometrical factor, $P_s$ is the polarization at saturation, $N$ is the number of obstacles or defects and $F_1$ is a function which is proportional to the pinning depth of the defect. As the two samples are strongly oriented, the polarization will have the same direction and the factor $f_0 P_s$ will be very similar. This can be seen from the P-E loops (Fig. 2) where the remnant polarization is almost identical for both samples. The difference of the parameter $\alpha_r'$ hence will be due to the product $NF_1$. This indicates that the a/c-domain sample, which has a higher parameter $\alpha_r'$, has a lower defect density and/or less deep defects. In the c-domain sample, defects are more important as the strain is not relieved due to the lower thickness. Moreover, the few domain walls present in this sample are
<figure><img src="image_15.png" /><figcaption>**Figure 4.** (a) Real and (b) imaginary part of the permittivity as a function of AC field amplitude (f = 10 kHz) for the two samples.</figcaption></figure>

<table><caption>**Table 1.** Coefficient of the lattice, domain wall pinning and domain vibration contributions for the real and the imaginary part of the total permittivity and the dissipation factors of the different contributions.</caption><thead><tr><th></th><th>ε'<sub>rl</sub></th><th>α'<sub>r</sub> (μm/V)</th><th>ε'<sub>r-rev</sub></th><th>ε''<sub>rl</sub></th><th>α''<sub>r</sub> (μm/V)</th><th>ε''<sub>r-rev</sub></th><th>m<sub>εrl</sub></th><th>m<sub>α</sub></th><th>m<sub>rev</sub></th></tr></thead><tbody><tr><td>c-domain sample</td><td>58.5</td><td>0.026</td><td>0.11</td><td>3.32</td><td>0.008</td><td>0.04</td><td>0.056</td><td>0.31</td><td>0.35</td></tr><tr><td>a/c-domain sample</td><td>73.9</td><td>1.832</td><td>1.49</td><td>1.72</td><td>0.635</td><td>0.51</td><td>0.023</td><td>0.35</td><td>0.34</td></tr></tbody></table>

blocked in deep defects. This can be seen from the threshold field for domain wall pinning/unpinning which is defined by ref. 9:

$$E_{th} = \frac{\epsilon'_{r-rev}}{\alpha'_{r}}, \qquad (3)$$

and is five times higher in the case of the c-domain sample. When the threshold field is low much domain walls pinning/unpinning occurs. This means that the domain walls are more free to move in the material. In the case of the c-domain sample a very high threshold field is observed ($E_{th} \sim 42.3$ kV/cm) and the domain walls are hence blocked. Switching the polarization becomes more difficult and consequently the coercive field is higher. The more important value of the threshold field in the c-domain sample is in agreement with the higher value of the coercive field observed in Fig. 2b.

Concerning the imaginary part of the domain wall motion coefficients ($\epsilon_{r-rev}$ and $\alpha_r$), the ratio between the two samples is quasi-identical. $\epsilon_{r-rev}''$ and $\alpha_r''$ are more than ten times higher in the a/c-domain sample. Consequently, the dissipation factors of the domain wall vibration $m_{rev}$ and of the pinning/unpinning $m_\alpha$ of the two samples are of the same order. These factors correspond to the ratio between the imaginary and the real parts of the respective contribution and are reported in Table 1. The $m_{rev}$ dissipation factors are equal to 0.34 and 0.35
and are similar to those that reported in the literature for BST thin films<sup>11, 18</sup>, indicating that the vibration of domain walls is a very dissipative phenomenon<sup>11</sup>. No information on the $m_{rev}$ values for PZT has been found in the literature. The dissipation parameters $m_\alpha$ are lower than those reported for PZT thin films ($m_\alpha \approx 0.42$)<sup>10, 21, 22</sup>. This might be due to the lower domain wall density in the two samples studied, as a higher density favors domain wall interaction and hence a higher dissipation factor<sup>18</sup>. The similar dissipation factor retrieved by our study indicates that the dissipation due to domain wall motion is not sensible to the domain wall density. According to Boser<sup>20</sup>, the ratio between the coercive field and the threshold field is equal to:

$$ \frac{E_C}{E_{th}} = \frac{2}{\sqrt{\pi}} \sqrt{\ln \frac{L_3}{2L_0}}, \qquad (4) $$

where $L_3$ is the distance between two domain walls and $L_0$ is the distance between two zero positions of the function which describes the distance between two defects. In our study the ratio is equal to 18 and 34 for the c-domain sample and the a/c-domain sample respectively, which is high compared to what is obtained by Boser ($E_C/E_{th} \approx 2.2$)<sup>20</sup>. The distance between domain walls is apparently so large that there is no interaction between them and the domain walls motions are not influenced by the domain wall density. As a consequence, the dissipation factors are identical for the two samples.

The dissipation of the lattice parameters $m_{crl}$ is two times higher in the c-domain sample than in the a/c-domain sample. This means that the higher dielectric losses in the c-domain sample is due to its structure, probably caused by the strain which is not relieved in the thinner film. Consequently, the a/c-domain sample has lower global dielectric losses.

## Conclusion

Two (001) oriented PZT samples of a different thickness have been studied. PFM characterization shows a majority of c-domains for the 60 nm thick sample whereas the sample of 100 nm thickness presents approximately 60% c-domains and 40% a-domains. The higher domain wall density in the c and a-domain sample is confirmed by the hyperbolic law measurements which show that the coefficients $\alpha'_r$ and $\epsilon'_{r-rev}$ describing domain wall motion are higher. In the case of the c-domain sample the threshold field $E_{th}$ is higher suggesting that the domain walls needs a higher energy to be unpinned. The mobility of the domain walls in this sample is reduced and a higher coercive field on the P-E loop is then observed. This sample presents also higher dielectric losses but the dissipation factors due to domain wall motions are similar to those of the c and a-domain sample. The higher dielectric losses seem to be due to the strain which increases the lattice losses. The combination of PFM and dielectric characterization in the present work has provided a more complete understanding of the dynamical behavior of domain wall motion in this system than has been previously possible. This study presents one aspect of the switching mechanism, which is the domain wall mobility under high field. Understanding this domain walls mobility allow us to better understand and potentially quantify the field dependence of ferroelectric switching.

## References

1. Garcia, J. E. et al. Evaluation of domain wall motion in lead zirconate titanate ceramics by nonlinear response measurements. *J. Appl. Phys.* **103**(5), 054108 (2008).

2. Bolten, D., Bottger, U. & Waser, R. Reversible and irreversible piezoelectric and ferroelectric response in ferroelectric ceramics and thin films. *Journal of the European Ceramic Society* **24**(5), 725–732 (2004).

3. Eitel, R. E., Shrout, T. R. & Randall, C. A. Nonlinear contributions to the dielectric permittivity and converse piezoelectric coefficient in piezoelectric ceramics. *J. Appl. Phys.* **99**(12), 124110 (2006).

4. Trolier-McKinstry, S., Gharb, N. B. & Damjanovic, D. Piezoelectric nonlinearity due to motion of 180° domain walls in ferroelectric materials at subcoercive fields: a dynamic poling model. *Appl. Phys. Lett.* **88**(20), 202901 (2006).

5. Taylor, D. V. & Damjanovic, D. Evidence of domain wall contribution to the dielectric permittivity in PZT thin films at sub-switching fields. *J. Appl. Phys.* **82**(4), 1973–1975 (1997).

6. Taylor, D. V., Damjanovic, D., Colla, E. & Setter, N. Fatigue and nonlinear dielectric response in sol-gel derived lead zirconate titanate thin films. *Ferroelectrics* **225**(1–4), 91–97 (1999).

7. Bharadwaja, S. S. N., Damjanovic, D. & Setter, N. Analysis of the nonlinear domain wall response in ferroelectric thin films. *Ferroelectrics* **303**, 59–63 (2004).

8. Damjanovic, D. Ferroelectric, dielectric and piezoelectric properties of ferroelectric thin films and ceramics. *Rep. Progr. Phys.* **61**(9), 1267–1324 (1998).

9. Borderon, C., Renoud, R., Ragheb, M. & Gundel, H. W. Description of the low field nonlinear dielectric properties of ferroelectric and multiferroic materials. *Appl. Phys. Lett.* **98**(1), 112903 (2011).

10. Garcia, J. E., Perez, R. & Albareda, A. Contribution of reversible processes to the non-linear dielectric response in hard lead zirconate titanate ceramics. *J. Phys. Condens. Matter.* **17**(44), 7143–7150 (2005).

11. Nadaud, K., Borderon, C., Renoud, R. & Gundel, H. W. Effect of Manganese Doping of BaSrTiO$_3$ on Diffusion and Domain Wall Pinning. *J. Appl. Phys.* **117**(8), 084104 (2015).

12. Kartawidjaja, F. C., Sim, C. H. & Wang, J. Heterolayered PZT thin films of different thicknesses and stacking sequence. *Journal of Materials Science* **44**(19), 5375–5382 (2009).

13. Lou, X. J. et al. Effect of polarization fatigue on the Rayleigh coefficients of ferroelectric lead zirconate titanate thin films: Experimental evidence and implications. *Applied Physics Letters* **105**(10), 102907 (2014).

14. Vrejoiu, I. et al. Intrinsic Ferroelectric Properties of Strained Tetragonal PbZr$_{0.2}$Ti$_{0.8}$O$_3$ Obtained on Layer-by-Layer Grown, Defect-Free Single-Crystalline Films. *Advanced materials* **18**(13), 1657–1661 (2006).

15. Morioka, H. S., Yokoyama, S., Oikawa, T., Funakubo, H. & Saito, K. Spontaneous polarization change with Zr/(Zr + Ti) ratios in perfectly polar-axis-orientated epitaxial tetragonal Pb(Zr,Ti)O$_3$ films. *Appl. Phys. Lett.* **86**(16), 3516–3518 (2004).

16. Nagarajan, V. et al. Size effects in ultrathin epitaxial ferroelectric heterostructures. *Appl. Phys. Lett.* **84**(25), 5225–5228 (2004).

17. Pintilie, L., Vrejoiu, I., Hesse, D., LeRhun, G. & Alexe, M. Ferroelectric polarization-leakage current relation in high quality epitaxial Pb(Zr,Ti)O$_3$ films. *Physical Review B* **75**(10), 104103 (2007).

18. Nadaud, K., Borderon, C., Renoud, R. & Gundel, H. W. Decomposition of the different contributions to permittivity, losses, and tunability in BaSrTiO$_3$ thin films using the hyperbolic law. *J. Appl. Phys.* **119**, 114101 (2016).
## Acknowledgements

M.A. acknowledges the Wolfson Research Merit and Theo Murphy Blue-sky Awards of Royal Society. The work was partly supported by the EPSRC (UK) through grant no. EP/M022706/1.

## Author Contributions

C. Borderon wrote the main manuscript text and prepared all the figures. She has also done all the fit for the hyperbolic law and the calculation of all the coefficients. A.E. Brunier has done the PFM measurements and P-E loop measurements. K. Nadaud has done the dielectric measurements (permittivity as a function of frequency and of the alternative field). R. Renoud has done the theoretical description of the domain walls. M. Alexe fabricated the samples and H.W. Gundel has supervised all the study. "All authors reviewed the manuscript".

## Additional Information

**Competing Interests:** The authors declare that they have no competing interests.

**Publisher's note:** Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

<figure><img src="image_16.png" /></figure>

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

© The Author(s) 2017

<footnote n="19">19. Xu, F. et al. Domain wall motion and its contribution to the dielectric and piezoelectric properties of lead zirconate titanate films. *Journal of Applied Physics* **89**(2), 1336–1348 (2001).</footnote>

<footnote n="20">20. Boser, O. Statistical theory of hysteresis in ferroelectric materials. *J. Appl. Phys.* **62**(4), 1344–1348 (1987).</footnote>

<footnote n="21">21. Garcia, J. E., Perez, R., Albareda, A. & Eiras, J. A. Non-linear dielectric and piezoelectric response in undoped and Nb<sup>5+</sup> or Fe<sup>3+</sup> doped PZT ceramic system. *Journal of the European Ceramic Society* **27**(13–15), 4029–4032 (2007).</footnote>

<footnote n="22">22. Eitel, R. E., Shrout, T. R. & Randall, C. A. Nonlinear contributions to the dielectric permittivity and converse piezoelectric coefficient in piezoelectric ceramics. *Journal of Applied Physics* **99**(12), 124110 (2006).</footnote>