(1−x)Pb(Mg1/3Nb2/3)O3-xPbTiO3 (PMN-PT) perovskite-like solid solutions are recognized for their outstanding electromechanical properties, which are of technological importance. However, some significant aspects of the crystal structures and domain assemblages in this system and the role of these characteristics in defining the functional performance of PMN-PT remain uncertain. Here, we used synchrotron x-ray diffraction to investigate the phase transition linking the paraelectric (cubic) and ferroelectric (tetragonal) phases in a single crystal of 0.65PMN-0.35PT. We analyzed the evolution of reciprocal-space maps across this transition. These maps were collected using small temperature step (1 K) and a high reciprocal-space resolution to reveal changes in the splitting of Bragg peaks caused by the formation of ferroelastic domains in the low-symmetry phase. Our results uncovered a two-phase state, cubic plus tetragonal phases, which exists over a narrow temperature range of only ≈4 K and exhibits a thermal hysteresis of ≈1.8 K. Remarkably, within this state, the lattice parameter of the cubic phase, , matches the orientational average of the lattice parameters for the tetragonal polymorph, . We discuss the implications of this matching, highlighting the possibility of it being realized by the formation of an assemblage of tetragonal twin domains separated from the cubic phase by a strain-free {110} boundary, as in the “adaptive phase” but without domain miniaturization.
(1−x)Pb(Mg1/3Nb2/3)O3-xPbTiO3 (PMN-PT) solid solutions attract interest because of their exceptional electromechanical properties. Single crystals of PMN-PT, with the composition near the morphotropic phase boundary (MPB), exhibit effective piezoelectric coefficients of ≈2500 pC/N,1,2 surpassing those of PbZr1−xTixO3 and α-SiO2 (quartz) by one and three orders of magnitude, respectively.3–6 In addition to its technological relevance, PMN-PT has been the subject of extensive academic research7–9 trying to understand the complex origins of the electromechanical response in this system.
The phase diagram of PMN-PT is featured by a transition between the high-temperature cubic paraelectric and the low-temperature tetragonal (x > 0.35) or rhombohedral (x < 0.28) phases (Fig. 1).10–14 These low-temperature phases are separated by a monoclinic phase (approximately, 0.28 < x < 0.35)15 that encompasses compositions displaying the enhanced properties and thus is of most practical and fundamental interest.16 All phase transitions in PMN-PT are associated with cooperative polar cation shifts,17 which remain partly disordered in the low-temperature phases with intricate local correlations.18,19 The nature of this correlated disorder resulting in complex assemblages of ferroelastic and ferroelectric domains remains debatable.
Here, we employ variable-temperature high-resolution reciprocal space (RS) mapping to probe the transition between the cubic (C) and tetragonal (T) polymorphs in a single crystal of 0.65PMN-0.35PT. We identify a narrow temperature range of ≈4 K, where these two phases coexist. This two-phase state exhibits a 1.8 K hysteresis between heating and cooling, supporting the first-order nature of this transition. We establish that the resulting assemblage of the cubic and tetragonal phases is characterized by a remarkable matching of the cubic and orientationally averaged tetragonal pseudocubic lattice parameters. The implications of this matching and its possible mechanism are discussed.
A single crystal of 0.65PMN-0.35PT was produced using a flux-growth technique. This method exploits spontaneous nucleation through the gradual cooling of a supersaturated solution from high temperature.20 Initial high-resolution RS mapping was performed on an “as-grown” unpoled sample using the four-circle x-ray diffractometer in Tel Aviv University.21 These experiments identified an “intermediate state” along the transition path between the C and T phases: Fig. 1 highlights the relevant phase-diagram region. Detailed studies were then carried out at the ID28 beamline22 of the European Synchrotron Radiation Facility. For these measurements, the same sample was etched to a cylinder-like shape with a diameter of ≈50 and a length of ≈200 . This sample preparation procedure was followed by two heating–cooling cycles, erasing any occasional poling history. The wavelength of the incident x-ray beam was , which yields an average x-ray penetration depth of . X-ray diffraction intensities were acquired with the PILATUS 1M detector by continuously rotating the crystal around the cylinder axis and saving the integrated detector images every . The measurements were performed at high-scattering angles, , ranging between 110 and 130 to enable highest possible RS resolution.
The data were collected at 25 temperatures between 403 and 458 K, with the sample temperature being controlled using the ID28-based heat blower. The anticipated phase-transition range (424–442 K) was scanned with the temperature step of 1 K. The heating/cooling rate was 0.3 K/min, but the RS measurements were taken only after achieving each target temperature. At each temperature, the data were used both to determine the crystal-orientation matrix using CrysAlis Pro software23 and to reconstruct high-resolution diffraction intensity distributions around 32 distinct Bragg peaks. The reconstruction of RS volumes (RSV) was accomplished using the custom MATLAB scripts, following procedures detailed elsewhere.21,24,25
Figure 2 shows three-dimensional diffraction intensity distribution around one of the collected Bragg peaks ( in this case), in the Cartesian coordinate system , with set nearly parallel to the scattering vector . The individual panels contain three projections of the intensity distribution along , and axes, defined as, e.g., . The individual sub-peaks in the RS maps are diffracted from a specific set of twin domains. The quantitative analysis of the separation between the sub-peaks can be done as described in Ref. 24. Specifically, the vectors between the individual components can be used as an input data for the assignment of the sub-peaks to the individual domain variants. The simplest version of this analysis relies on the RS separation along the Cartesian -axis because the component is nearly equal to the length of the scattering vector . The presence of one-to-one correspondence between and the scattering angle ( means that this separation can be observed using high-resolution powder diffraction experiment26–28 and allows determination of the free lattice parameters ( and for the case of T-domains).
The supplementary material contains an animated version of the same figure, showing the temperature dependence of the RSV around all the collected Bragg peaks. The transition to the C-phase is evidenced by the merging of all the sub-peaks. Figure 3(a) illustrates this temperature dependence using a one-dimensional intensity distribution . Panels (b)–(d) in Fig. 3 add three intensity projections at 440, 437, and 426 K, which correspond to the nominally cubic, “intermediate,” and tetragonal states, respectively. Figure 3(c) reveals a different topology of the sub-peaks, compared to that in the nominally tetragonal region [Fig. 3(d)]. The “intermediate state,” sandwiched between the C- and T-phases, exists over a narrow temperature interval of 4 K and is therefore difficult to detect. The C↔T transition temperature determined by our measurements is ≈10 K higher than that reported previously.10–14 Possible reasons for this discrepancy include differences between the ceramic samples used to determine the published phase diagrams and the single crystal studied here, and the coarser temperature steps employed in previous works.
Figure 3(c) confirms that one of the subpeaks (labeled “C”) ascribed to the intermediate state can be treated as a continuation of the corresponding cubic peak—an observation common to all the collected Bragg-peak maps (see the supplementary material). These subpeaks can be used to estimate the corresponding volume fraction of the C-phase as . Here, the sums in the numerator and the denominator run over the sub-peaks attributed to the C-phase and all the sub-peaks, respectively. The peak intensities, , were determined by integrating the intensity distributions over the radius around selected sub-peaks. Figure 4 displays the resulting dependence, highlighting the 4 K wide phase-coexistence region and the 1.8 K-wide hysteresis observed during the heating-cooling cycle. The presence of this hysteresis suggests that the observed phase coexistence is intrinsic rather than caused by a compositional segregation and consistent with the first-order type of the C↔T transition in PMN-PT.
Figure 5 reveals a remarkable matching of the “cubic” and “averaged tetragonal” lengths of the lattice basis vector in the phase-coexistence range. The difference is less than , which is below the relative standard uncertainty of , , and .
Thus, HP between the C phase and ATP can reside on the (011) plane (Fig. 6). The basis vectors, forming the corresponding 2D unit cell in this plane, would be and (these directions are explicitly displayed in Fig. 6). For the C-phase, the lengths of these vectors are and , respectively. For the ATP, these lengths are and . Neglecting the second power of , we can state that the (101)-micro-twinned tetragonal phase can grow mismatch-free on the (011) plane of the C-phase.
The original “adaptive phase” concept30–33 involves the miniaturization of tetragonal domains along the domain-wall normal, which, on average, results in the effective monoclinic symmetry. Our experiments show that the tetragonal domains in the two-phase state are large enough to produce Bragg peaks, which can be clearly assigned to the tetragonal, rather than a monoclinic phase. Assuming that the domain size is the only contribution to the width of the observed Bragg peaks (0.002–0.005 Å−1), we estimate this size as several hundred nanometers. Similar domain sizes, for 0.65PMN-0.35PT and other compositions near the MPB, were observed in Refs. 34–36. We propose that the lattice matching of the C- and T-phases in the intermediate, two-phase state observed here is enabled by the assembly of the T-domains in a fashion similar to ATP but without the domain miniaturization. We used the splitting of reflections in our 3D RS reconstructions to confirm the corresponding {110}-type HP orientations; the details of this analysis will be reported separately. Real space characterization of domain arrangements can be continued using, e.g., dark-field x-ray microscopy.37–39
In conclusion, we performed in situ high-resolution x-ray diffraction on a single crystal of 0.65PMN-0.35PT using a fine-step sampling of temperatures across the cubic-tetragonal phase transition. We discovered a two-phase state, only 4 K wide, characterized by the lattice parameter ( ) of the C-phase matching the orientationally averaged ) lattice metrics of the T-phase. This matching effectively corresponds to the unit-cell volume conservation. In principle, such conditions enable {110}-type habit planes connecting the C-phase and a micro-twinned assembly of T-domains.
SUPPLEMENTARY MATERIAL
See the supplementary material for the figures, similar to Figs. 2 and 3 but for all the measured Bragg peaks and temperatures.
We acknowledge Youli Li (University of California, Santa Barbara) for his help with and support of the Tel Aviv University-based x-ray diffractometer and the temperature controller. We acknowledge the following funding: Israel Science Foundation (Grant Nos. 1561/18, 3455/21, 1365/23 to Semën Gorfman); United States – Israel Binational Science Foundation (Award No. 2018161 to Semën Gorfman and Igor Levin); the U.S. Office of Naval Research (ONR Grant No. N00014-21-1-2085); and the Natural Science and Engineering Research Council of Canada (NSERC, DG, RGPIN-2023-04416).
Certain equipment, instruments, software, or materials, commercial or noncommercial, are identified in this paper in order to specify the experimental procedure adequately. Such identification is not intended to imply recommendation or endorsement of any product or service by NIST, nor is it intended to imply that the materials or equipment identified are necessarily the best available for the purpose.
AUTHOR DECLARATIONS
Conflict of Interest
The authors have no conflicts to disclose.
Author Contributions
Ido Biran: Data curation (equal); Investigation (equal); Methodology (equal); Validation (equal); Visualization (equal); Writing – original draft (equal); Writing – review & editing (equal). Alexei Bosak: Funding acquisition (equal); Investigation (equal); Methodology (equal); Resources (equal); Validation (equal); Writing – original draft (equal); Writing – review & editing (equal). Zuo-Guang Ye: Funding acquisition (equal); Investigation (supporting); Methodology (supporting); Resources (equal); Validation (equal); Writing – original draft (equal); Writing – review & editing (equal). Igor Levin: Conceptualization (lead); Funding acquisition (supporting); Investigation (supporting); Methodology (supporting); Project administration (lead); Resources (equal); Validation (equal); Writing – original draft (supporting); Writing – review & editing (equal). Semën Gorfman: Conceptualization (lead); Funding acquisition (lead); Investigation (equal); Methodology (equal); Project administration (lead); Resources (equal); Validation (equal); Writing – original draft (lead); Writing – review & editing (equal).
DATA AVAILABILITY
The data that support the findings of this study are available from the corresponding author upon reasonable request.