Bipolar loop-like non-volatile strain in the (001)-oriented Pb(Mg1/3Nb2/3)O3-PbTiO3 single crystals.

Strain has been widely used to manipulate the properties of various kinds of materials, such as ferroelectrics, semiconductors, superconductors, magnetic materials, and "strain engineering" has become a very active field. For strain-based information storage, the non-volatile strain is very useful and highly desired. However, in most cases, the strain induced by converse piezoelectric effect is volatile. In this work, we report a non-volatile strain in the (001)-oriented Pb(Mg1/3Nb2/3)O3-PbTiO3 single crystals and demonstrate an approach to measure the non-volatile strain. A bipolar loop-like S-E curve is revealed and a mechanism involving 109° ferroelastic domain switching is proposed. The non-volatile high and low strain states should be significant for applications in information storage.

Strain has been widely used to manipulate the properties of various kinds of materials1234, such as ferroelectrics1, semiconductors2, superconductors3, magnetic materials4, etc., and “strain engineering” has become a very active field1. One popular way to introduce strain is via lattice misfit between the epitaxial film and the substrate, and dramatic variation of strain can be realized by selecting suitable film and substrate1. However, in this case, the strain can not be tuned continuously and reversibly, and different substrates may induce other effects besides the strain effect when comparing the properties of films grown on different substrates. In this regard, converse piezoelectric effect is very useful for introducing continuous and reversible strain by electric fields. Recently the strain-mediated multiferroic heterostructures, composed of magnetic thin films and ferroelectrics (FE), have attracted much attention due to the magnetoelectric coupling between magnetic thin films and FE567891011121314. In these multiferroic heterostructures, strain is produced via the converse piezoelectric effect of FE and can be modulated continuously and reversibly. This kind of strain becomes a primary pathway to control the electronic transport and magnetic properties of the magnetic films in the multiferroic heterostructures, and “straintronics” has been proposed15. Owing to the superior piezoelectric properties16, relaxor ferroelectrics such as (1-x)Pb(Mg1/3Nb2/3)O3-xPbTiO3 (PMN-PT) single crystals, which are binary solid solution system (x represents the nominal composition of PbTiO3), have been widely used in the strain-mediated multiferroic heterostructures, besides their applications in ultrasonic devices, telecommunication, medical imaging17. So far, most reported in-plane strain-electric field (S-E) behavior of the (001)-oriented PMN-PT is a symmetric butterfly-like curve without remanent strain18, i.e. the electric-field induced strain is volatile. However, non-volatile strain is more important and highly desired for the strain-based information storage in the multiferroic heterostructures19. For electric-field induced reversible non-volatile strain, there are only a few reports on FE materials with remnant strain states in the S-E curves, which can be classified into two types. For type I, the polarization-electric field (P-E) hysteresis loop is offset about the origin as shown in polycrystalline PbZr1-xTixO3 due to an internal field which is generated by different kinds of defects20. Because the correlation between strain and polarization, the S-E curves also deviate from the symmetric butterfly-like curves and show an asymmetric behavior. In this case, there is indeed remanent strain at zero electric field, but the asymmetric behavior of the butterfly-like curves are determined by the polarity of the space charges21, accompanied with fatigue2223. For type II, Wu et al.242526 recently obtained the remanent strain in the (011)-oriented PMN-PT single crystals. When a critical reversed electric field was applied, the polarizations reoriented to the four possible in-plane <111> directions, which acted as a metastable state during the polarization reversal process. After decreasing the electric field to zero, hysteresis occurred for both P-E and S-E curves and two different strain states arise. However, a critical reversed electric field is required to induce the effect in these systems. If the electric fields applied to PMN-PT are larger or smaller than the critical field, the remanent strain will decrease dramatically24. Since the critical reversed electrical field is a bit smaller than the coercive field, the cycle of the electric field applied to the (011)-oriented PMN-PT is actually unipolar as mentioned in Wu's paper24, while for the bipolar case, there is no remanent strain24. In this paper, we report a non-volatile strain in the (001)-oriented PMN-PT single crystals. A bipolar loop-like S-E curve is revealed and a mechanism involving 109° ferroelastic domain switching is proposed. We also put forward and demonstrate an approach for measuring the non-volatile strain.

Results

Methods for strain measurement and loop-like strain behaviors in PMN-PT single crystals

Figure 1a is the schematic of the sample and the experimental configuration. (001)-oriented (1-x)PMN-xPT (x = 0.28, PMN-PT with 28% of PT) single crystals were used in our work. The single crystal orientation and microstructure of the polished surface were characterized as shown in Fig. S1 in Supplementary Information. Figure 1b shows the P-E hysteresis loop for PMN-PT measured by using a Radiant Technologies Precision Premier II system. The in-plane S-E curves were measured by using the continuous method as shown in Figure 1c, which has been commonly used in the literature1718. We measured the in-plane strains along the [110] and [1-10] directions of the (001) plane, whose converse piezoelectric responses are more sensitive27. The S-E curve along the [110] direction is shown in Figure 2a. In contrast to the normal symmetric butterfly-like curve18, it shows an asymmetric behavior with two remanent strain states at zero electric field. The in-plane S-E curve along the perpendicular direction [1-10] also shows an asymmetric butterfly-like behavior (Figure S2 in Supplementary Information) and is complementary with that of the [110] direction, providing a clue to the origin of the phenomenon as discussed later.

Figure 1

Experimental configuration, ferroelectric property and two different methods for measuring strain.

(a) Schematic of the sample and the experimental configuration. (b) P-E loop of PMN-PT. (c) and (d) Schematic of the continuous and pulsed measurement methods for S-E curves, respectively. The strain measurement starts at the end of interval as pointed by the solid blue triangle symbols.

Figure 2

Electric-field tuned strain behavior along [110] direction.

(a) and (b) S-E curves along the [110] direction measured by continuous and pulsed measurements, respectively. (c) The curve obtained by subtracting the strain values in the loop-like S-E curves (Figure 2b) from those in asymmetric butterfly curves (Figure 2a). (d) The stable and remarkable high/low remanent strain states achieved by switching the polarity of the pulsed electric field.

As shown in Fig. 1b, the hysteresis loop of our sample is symmetric about the origin, which means that the internal field induced by defects is negligible. This is different from the aforementioned type I asymmetric S-E curve, which is related to the asymmetry of the P-E hysteresis loop. Moreover, unlike the unipolar S-E curves for the aforementioned type II24, the S-E curve of our sample is bipolar. In addition, the electric fields applied in our work (±8 kV/cm) are much higher than the coercive field (1.2 kV/cm) characterized in the P-E loop (Fig. 1b). So the sample is in the saturated polarization states, which is different from the aforementioned type II case. Thus, a different mechanism for the asymmetry of S-E curves should be involved in our samples.

In order to get the non-volatile part in the asymmetric butterfly-like curves, we proposed an approach as shown in Fig. 1d, which is denoted as the pulsed measurements method. The detailed descriptions for the continuous and pulsed measurements are in the experimental section. In both measurements, we begin to measure the strain at the end of the interval. The difference is that electric voltage is always applied to the sample in the continuous measurement, while no electric voltage during the interval for the pulsed measurement. Thereby, in the pulsed measurement, the strain is measured at zero electric filed and reflects the remanent strain. Figure 2b shows the S-E curve obtained by the pulsed measurement. Interestingly, it exhibits a bipolar loop-like behavior, similar to the P-E hysteresis loop (Fig. 1b). The loop-like S-E curve suggests that one can easily control the strain states by electric fields. For the perpendicular [1-10] direction, we also got the asymmetric butterfly-like and loop-like S-E curves (Fig. S2), which exhibit a reverse behavior compared to those of the [110] direction. Figure 2c shows the S-E curve obtained by subtracting the strain values in the loop-like curve (Fig. 2b) from those in the asymmetric butterfly-like curve (Fig. 2a) at the same voltage point by point. The S-E curve in Fig. 2c is a symmetric butterfly-like curve, which is common in the previous reports18. Similar results were also obtained along the [1-10] direction (Fig. S2). Based on the above results and analyses, it can be deduced that the observed asymmetric butterfly-like S-E curves are composed of the symmetric butterfly-like and the loop-like curves. Stable and remarkable high/low remanent strain states can be achieved by switching the polarities of the pulsed electric fields as shown in Fig. 2d for the [110] direction and Fig. S2d for the [1-10] direction. It can be seen that the pulsed electric field polarity dependence of the high/low strain states for these two directions is converse, which is consistent with the complementary behavior in their S-E curves (Fig. 2 and Fig. S2). Instead of a specific value or range of electric field required in the type II case, any fields larger than the coercive field of the sample can induce the remarkable high/low non-volatile strain states. Moreover, as the sample is in the saturated polarization states, the non-volatile strain states were stable and didn't show an obvious relaxation with time, which is in contrast to the behavior of the (011)-oriented PMN-PT28. And the high/low strain states can be switched by pulsed electric fields for thousands of times without obvious fatigue effect (Fig. S3 in Supplementary Information).

Proportions of ferroelastic switching in PMN-PT single crystals with different S-E behaviors

It should be mentioned that not all the samples exhibit the asymmetric butterfly-like S-E curves. We also obtained the symmetric butterfly-like curve as shown in Fig. S4a in Supplementary Information, similar to that reported by Thiele18, in samples with the same nominal composition as those showing the asymmetric butterfly-like behavior. Figure S4b shows the S-E curves measured by the pulsed measurement. Obviously, there is no remanent strain at zero field in this kind of samples. Hereafter, we denote the samples with the asymmetric and symmetric S-E curves as sample A and sample S, respectively. To identify the differences between these two kinds of samples, we carried out XRD-RSM with in situ electric fields since it can provide straightforward information on the microstructure of the samples and the switching category of FE domains. PMN-PT samples are in the rhombohedral phase with eight spontaneous polarizations along the <111> directions, denoted as r1+/r1−, r2+/r2−, r3+/r3−, r4+/r4− as shown in Fig. 3a and 3b. When the polarity of electric field was changed, the spontaneous polarizations can be switched by 71°, 109° and 180°, respectively. Polarization variants changing from r1+/r3+ to r2−/r4− or r2+/r4+ to r1−/r3− (vice versa) correspond to the 109° switching, while polarization variants changing from r1+ to r3− or r3+ to r1− (the same case for r2/r4) correspond to the 71° switching, and polarization variants without any category change (e.g. r1+ to r1− and so on) correspond to the 180° switching. Figure 3c and 3d show the reflections around the (113) peak for sample A under −8 kV/cm and +8 kV/cm, respectively. The distortions corresponding to polarization variants r2 and r4 are degenerate in the (113) reflection peak, so there are three spots corresponding to r1, r3, and r2/r4, respectively. According to the lattice parameters in the rhombohedral structure29 and the analytical method30, we analyzed the positions and intensities of the spots (see SI for details), which can give the percentages of different polarization variants under positive and negative electric fields, respectively. Table 1 shows the results for sample A, suggesting that the proportion of r2/r4 under +8 kV/cm is about 20% less than that under −8 kV/cm. Since the change between r2/r4 and r1/r3 from positively poled cases to negatively poled cases (vice versa) corresponds to the 109° switching, the 20% reduction of r2/r4 proportion represents the net proportion of 109° switching. For comparison, we also measured the reflections around the (113) peak with in situ electric fields for sample S. The results are shown in Figure S4c and S4d. Similar to sample A, the (113) reflection peak also shows three spots. We analyzed the data using the same method as for sample A and the results are shown in Table S1. The percentages of different polarization variants are almost unchanged with the polarities of electric field. This means that there are 71/180° switching, but almost no net 109° switching since the percentage of r2/r4 is nearly unchanged under positive and negative electric fields.

Figure 3

The schematic diagram of the polarization variants and XRD-RSMs around (113) reflections.

(a) and (b) Configurations of the (113) crystal face (translucent light blue plane) and polarization variants of rhombohedral phase under negative and positive electric fields, respectively. (c) and (d) (113) reflection peak of XRD-RSMs for sample A with in situ electric fields of −8 kV/cm and +8 kV/cm, respectively. Ellipses with different colors in RSMs outline the contours of the deduced spots.

Table 1

Proportions of different polarization variants under positive and negative electric fields, respectively

E (kV/cm)	r1	r2/r4	r3	Correlation Coefficient
−8	41.7%	55.5%	2.8%	0.96
+8	17.1%	35.6%	47.3%	0.94

The crucial role of the net 109° switching for the non-volatile strain

In sample A, there are not only 71/180° switching but also about 20% of net 109° switching. While in sample S, there are only 71/180° switching and almost no net 109° switching. Combined with the complementary property in the loop-like curves in sample A (Fig. 2b and S2b) and the symmetric butterfly-like curves in sample S (Fig. S4a), it can be concluded that the loop-like S-E curve is related to the net 109° switching. In contrast, without net 109° switching, the S-E curve exhibits the butterfly-like behavior. Figure 4 is the schematic diagram of the 109° switching, 71/180° switching and corresponding changes of distortion/strain. For polarization switching from P+ to P−, we take r1+ and r3+ as an example to show the 109° switching (right) and 71/180° switching (left) as shown in the upper panel of Fig. 4. It can be seen that the projections of polarizations to the (001) plane before and after the 109° switching are perpendicular, while they are parallel before and after the 71/180° switching. As a result, the in-plane tensile distortions related to the polarizations are perpendicular before and after the 109° switching and parallel before and after the 71/180° switching (lower panel of Fig. 4). It is also noted that change of the tensile distortions along the [110] and [1-10] directions are complementary for the 109° switching. This is the origin of the complementary loop-like curves of sample A along the two perpendicular directions. In sample S, there are only 71/180° switching and the S-E curves show a symmetric butterfly-like behavior without remanent strain.

Figure 4

Correlation between domain switching and distortion.

Schematic diagrams of the 109° switching, 71/180° switching (upper panel) and corresponding changes of distortion (lower panel).

Discussion

An interesting question arises as why sample A and sample S, with the same nominal composition, show such a big difference in the S-E curves, related to the presence and absence of net 109° switching in sample A and sample S, respectively. It is likely that defects/defect configurations in the crystals may play an important role. As reported by Fu et al.31, polarization rotation along the different paths needs to overcome different energy barriers. And it has been shown that defects can strongly affect the barrier for polarization rotation, e. g. defects make the 180° polarization flip become a feasible process32. It was also shown that defect configurations play an important role in the polarization switching33. These may shed light on our case as follows. Some defects/defect configurations in sample A lowers the energy barrier between the two polarization configurations with the polar axis rotating by 109°, which is crucial for the loop-like S-E curves. In contrast, it is not the case for sample S and 109° domain switching is not favored due to the required high energy. The presence of key defects/defect configurations may be due to details during the crystal growth process. Further study on this issue is deserved.

In conclusion, we report a non-volatile strain in the (001)-oriented PMN-PT single crystals. A bipolar loop-like S-E curve is revealed and a mechanism involving 109° ferroelastic domain switching is proposed. We also put forward and demonstrate an approach for measuring the non-volatile strain. As a comparison, the non-volatile strain in the (011)-oriented PMN-PT reported previously is dependent on a critical reversed electric field and decreases dramatically if the electric field deviates from the critical value. While the non-volatile strain in the (001)-oriented PMN-PT in our work can be obtained by any fields larger than the coercive field and maintains a stable state. The non-volatile high and low strain states should be significant for applications in information storage.

Methods

Our samples are one-side-polished (001)-oriented PMN-PT with the edges cut along the pseudocubic [110] and [1-10] directions, and dimensions of 5 mm × 5 mm × 0.5 mm. Au layers of 100 nm thick were deposited on both sides of the samples as electrodes. The strain gauge was pasted on the sample surface with glue M-Bond 610, and kept at 120°C for 2 h to strengthen the paste effect.

Two kinds of measurements for S-E curves are shown in Fig. 1c and 1d, which are named as continuous and pulsed measurements, respectively. In the continuous measurement, the voltage applied to the sample is increased with a step of 20 V. After a period of delay, an interval is followed. At the end of the interval, the measurement begins, and then the voltage is changed for the next cycle. So the voltage is always applied. In the pulsed measurement, the voltage is also increased with a step of 20 V. After a period of delay, the voltage is removed for an interval. At the end of the interval, the strain measurement starts. After that the voltage is changed and the measurement goes into the next cycle. We carried out the XRD-RSM experiments with in situ electric fields of ±8 kV/cm using the synchrotron radiation x-ray source. All the measurements were carried out at room temperature.

Author Contributions

L.Y. and Y.Z. planned the experiments. L.Y. made the samples and performed all the strain measurements with the assistance from S.Z., P.L., P.M., Y.G. and C.N. The data analysis was performed by L.Y., Y.L. and A.C. The XRD-RSM experiment and analysis were performed by Y.Y., H.H. and C.G. The paper was written by L.Y. and Y.Z. All authors reviewed and commented on the manuscript.