Abstract
Revelation of unequivocal structural information at the atomic level for complex systems is uniquely important for deeper and generic understanding of the structure property connections and a key challenge in materials science. Here we report an experimental study of the local structure by applying total elastic scattering and Raman scattering analyses to an important non-relaxor ferroelectric solid solution exhibiting the so-called composition-induced morphotropic phase boundary (MPB), where concomitant enhancement of physical properties have been detected. The powerful combination of static and dynamic structural probes enabled us to derive direct correspondence between the atomic-level structural correlations and reported properties. The atomic pair distribution functions obtained from the neutron total scattering experiments were analysed through big-box atom-modelling implementing reverse Monte Carlo method, from which distributions of magnitudes and directions of off-centred cationic displacements were extracted. We found that an enhanced randomness of the displacement-directions for all ferroelectrically active cations combined with a strong dynamical coupling between the A- and B-site cations of the perovskite structure, can explain the abrupt amplification of piezoelectric response of the system near MPB. Altogether this provides a more fundamental basis in inferring structure-property connections in similar systems including important implications in designing novel and bespoke materials.
Similar content being viewed by others
Introduction
Ferroelectric materials with a composition-driven structural crossover, commonly known as morphotropic phase boundary (MPB), have become an indispensable part of many modern devices, particularly used as sensors, actuators and memories, utilizing their superior properties at the MPB. The term MPB, which literally refers to a boundary between two forms, was first coined to describe the chemically induced change in the ferroelectric long-range order of the famous PbZr x Ti1−x O3 (PZT). Since then, it is now a well established fact that invoking a structural instability by tweaking the composition may result in anomalous characteristics similar to PZT. The current understanding of an MPB and the associated enhancement of certain physical properties of a ferroelectric material primarily relies on the fact that the system acquires a state where the rotation of the unit-cell polarization vector becomes easier due to the development of additional degrees of freedom either in a single low-symmetry phase or in several coexisting phases. The concept of a bridging low-symmetry phase became noted after the discovery of a monoclinic Cm phase at the MPB of PZT1,2,3 and manifested a renewed interest in studying as well as in designing bespoke complex ferroelectric materials. However in recent years it has also been demonstrated that the analysis of the gross average structure broadly simplifies the ubiquitous and varied complexity of the mesoscopic-scale structural features, which are crucial to understand the occurrence of anomalous properties at the MPB4,5,6,7,8,9,10. Especially, with many competing interactions and structural frustrations at the atomic level, multi-component ferroelectric materials in general provide a unique challenge in developing precise structural models that would correspond to the physical properties11. Consequently, it is still a riveting topic to investigate aspects of structural behaviour even for immensely studied systems like PZT, in search of a more rigorous model than the existing concept of easy rotation of the net polarization12,13,14. It is also highly anticipated that the understanding would not be complete unless the static structural models are equally complemented by dynamical information, typically obtainable from inelastic scattering processes15,16,17,18.
Driven by the motivation of finding Pb-free or reduced-Pb alternatives of PZT, there has been a strong interest in Bi-containing ferroelectric solid solutions with the general formula xBiMeO3-(1 − x)PbTiO3, Me = Sc, Fe, Mg1/2Ti1/2, Ni1/2Ti1/2, Ni1/2Zr1/2 etc following the revelation of the MPB features in those systems19, 20. For example, xBiScO3-(1 − x)PbTiO3 (BS-PT) exhibited even superior physical properties with a higher operational range of temperature than PZT21, 22. The successful partial substitution of Pb by Bi without compromising cherished propertied of pure-Pb containing systems was considered as an important milestone for developing eco-friendly materials. Therefore BS-PT, in the form of ceramics, thin films as well as single crystals with desirable properties, gained rapid interest as a potential replacement of PZT. However the high price of Sc2O3 is making this particular system less accessible. Although not as attractive as BS-PT, xBiMg1/2Ti1/2O3-(1 − x)PbTiO3 (xBMT-PT), first reported in 2004, provides a reasonable alternative to the expensive Sc with promising MPB-properties at x ≈ 0.6323. It is also considered as an interesting system as BMT was reported to be a structural analogue of anti-ferroelectric PbZrO3 24. In addition, xBMT-PT has also been uniquely shown to possess zero thermal expansion coefficient in the range 0.2 ≤ x ≤ 0.4, and highly stable piezoelectric properties at non-ambient temperatures25, 26. In terms of the structural phase transition driven by the composition, very recently Upadhyay et al.27 proposed that there is a tetragonal (P4 mm) to monoclinic (P m) phase transition through a mixed phase (P4 mm + P m) region that exists in the range 0.60 ≤ x ≤ 0.67, based on the Rietveld refinements of the powder XRD pattern.
Morphotropic phase boundary in ferroelectric solid solutions has so far been detected through the typical average structural investigation, mostly applying standard powder diffraction technique for perovskite based oxide systems. Although there have been rigorous attempts to study particularly the local structural correlations, such as diffuse scattering studies on Pb-based complex systems revealing strong evidence for large deviations from the average structure8,9,10, 28,29,30, there is still lack of an experimentally conceived model at the atomic level to identify and correlate the properties to the different facets of the structure, which can then serve as a more fundamental basis in finding and designing superior as well as eco-friendly materials. Hence, there is a pressing need to elucidate mesoscopic-scale atomic correlations in ferroelectric solid solutions to better understand the physics of MPB.
With the present availability of high energy x-ray synchrotron facilities and spallation neutron sources, it is nowadays possible to obtain data for a wide range of reciprocal-lattice vectors Q for powder samples, and thereby to consider Fourier transformation of the data, taking both Bragg diffraction and diffuse scattering on equal weight. This is known as total scattering method which manifests pair distribution functions (PDFs). PDFs essentially describe the whole structure in terms of atom-atom distances weighted by their scattering power. Hence PDFs are critically sensitive to variations of the local correlations and are considered categorically as a powerful local probe31. The total scattering method has already been applied to a number of popular ferroelectric systems including both Pb-based and Pb-free compounds, and revealed hitherto unseen structural characteristics, such as large and persistent static displacements of cations from their crystallographic sites, distinct local and average polarisation, chemical ordering, formation of polar nano-regions as well as their developments with composition and temperature4, 5, 11, 13, 17, 32,33,34,35,36,37.
In this report we show the composition-driven evolution of the local cation-environment in xBMT-PT as well as their dynamical behaviour through a combined analysis of neutron PDFs and Raman scattering data at ambient conditions for the compounds covering the whole composition-range of stability across the phase diagram, which has not been reported so far. Our experimental results provide deeper insights into the structural phenomena occurring at the MPB of a perovskite-type ferroelectric solid solution and help to establish a comprehensive structure-property relationship for a broad range of systems.
Results and Discussion
Pair distribution function analysis
Figure 1(a and b) shows the development of the PDFs for xBMT-PT as a function of composition derived from the neutron total scattering data along with the {001} pc (pc refers to the pseudocubic setting) Bragg peaks extracted from in-house XRD data (Fig. 1c). The long-range correlations in the PDFs as shown in Fig. 1b in the range 40–46 Å exhibit a distinct change at the MPB composition, which can be easily linked to the observed trend in the Bragg peaks with x. However, the short-range correlations up to 8 Å, where peaks can be uniquely associated with the different first-neighbour atom-atom distances in the perovskite-type structure, do not imply for an abrupt change to single out the MPB composition (Fig. 1a). Nevertheless, it appears that the peaks of the PDF at x = 0.10 are most pronounced and then it gradually broadens with increasing x. The evident gradual changes continues up to the MPB composition and become almost negligible for the compositions x = 0.63, 0.65 and 0.70 in the distances longer than approximately 5.1 Å. The changes seen in the peak shape with increasing x, suggest enhanced structural disorder, which indicates the composition-induced structural phase transition is more of a order-disorder type than a displacive type. This distinction from a pure displacive-type structural phase transition is an important tag for a ferroelectric system and has impact on the properties especially under external stimuli and non-ambient conditions.
In order to extract more specific and quantitative information about the local structural changes across the MPB, we have carried out typical big-box modelling based on the reverse Monte Carlo (RMC) method against the experimental PDF. The analyses were restricted to the atomic-distance range 1–20 Å to ensure that the resulting model solely reveals the local correlations38, 39. Very recently, there have been a few reports on popular ferroelectric systems containing Pb and/or Bi using similar big-box modelling applying RMC technique, which have shown how the local structural diversity could be related to the anomalous macroscopic properties32, 34, 36, 37, 40. A similar method, first described by Keeble et al.35, was adopted here to map the distribution of the cationic displacement directions with respect to the centre of their corresponding oxygen polyhedra on a stereograph. These graphs essentially help to evaluate the behaviour of the individual cations with x. The magnitude and the direction of the polar displacements of the cations are both crucial for a ferroelectric material since it affects the polarization as well as the structure of the system, and consequently asserts the ensuing macroscopic properties6. Figure 2 shows such stereographs for different cations as a function of x. For both Pb and Bi, it is apparent that at low values of x the directions of off-centre displacements are consistently along the [001] pc direction (tetragonal distortion) with a gradual dispersion of the high-density region as x increases. In addition, an abrupt enlargement of the dispersion can be seen particularly for the A-site cations from x = 0.60 to 0.63 (see Fig. 2d). The strong coupled behaviour of the A-site cations could be justified with their similar electronic properties, however it should be noted that neutron scattering cannot distinguish well between Pb and Bi as they have very similar scattering lengths. Although Ti seems to follow the A-site cations with the composition, the change in the direction-distribution is however smoother, with a relatively higher level of dispersion on the stereographs for all x compared to that of the A-site cations. Mg seems to behave uniquely as it tends to scatter initially on the {001} pc planes. However for x ≥ 0.50 that preference becomes weaker, and at x = 0.70, the Mg shifts are predominantly along the [001] pc .
Apart from these direct observations, the maps can also help to derive further characteristics of the local structure. On the basis of the density distribution of the points in Fig. 2, it is apparent that the inclusion of Bi3+ and Mg2+ in PbTiO 3 triggers a structural instability at the local level, and the system is evolving from a strongly anisotropic to a rather isotropic in terms of the displacement directions of the individual cations, which could equally be seen as a reduction of the energy difference between the different orientation-states with the increase of x. In addition, this statistical variation of directions with x suggests the system gradually becomes more pliable to accommodate non-collinear polar shifts of the cations while approaching MPB from the tetragonal side. This feature in general has direct consequences on the properties of a system, where polarization is coupled to strain and the strain is considered as a primary constraint to the cooperative switching of the dipoles under external fields41. Altogether this can be seen as a direct observation of a thermodynamical picture where the individual potential functions for cations are flattened, which in turn provides greater flexibility for the local polarization6, 42.
Considering the densities of the favoured directions for x ≥ x MPB, the maps cannot refer a specific average symmetry for the system. Although the x-ray powder diffraction results suggest a single monoclinic Pm phase above the MPB (x = 0.70), where the polarization vector should remain preferably within the {100} pc planes (see Fig. 2c), it is inconclusive from the gross trend in Fig. 2d. However, the zoomed-in central regions of the stereographs around [001] pc for Pb and Ti (see Fig. 2e, Bi shows the same trend as Pb) further reveal that there is a subtle yet distinct difference in the direction-trends on a local scale. In particular, Ti exhibits an inclination to the [111] pc -type polarization even at x = 0.10 in comparison to the A-site polar displacements. The A-site cation polar shifts exhibits a deviation in tendency from [001] pc to [111] pc in the range 0.30 ≤ x ≤ 0.60, and beyond that it is highly degenerate. However at x = 0.70, there is a weak preference of staying on the {110} pc plane, whereas Ti polar displacements shows a tendency to be on the {100} pc plane. This as a whole, can be related to all types of monoclinic distortions (M A , M B and M C ), commonly considered for ferroelectric systems based on the perovskite structure. Nevertheless, this emphasizes the fact that the local structural disorder is greatly enhanced with large values of x for the family of xBiMeO3-(1 − x)PbTiO3, where the Bi-containing end members are mostly not stable at ambient conditions. As a result, conventional analyses of the powder diffraction data through the Rietveld method often face challenges in providing unambiguous average structural description particularly for x ≥ x MPB 36.
In order to quantify the observed distribution of the points on the stereographs, we have used a common orientation order parameter S = (1.5〈cos 2 θ〉 − 0.5)43, where θ represents the angle between the cation displacement direction and the [001] pc direction. The cosine is averaged over all displacement directions, and weighted by the corresponding density from the stereographs. The variation of values of S in Fig. 3 depicts the sequence of transition with x by a measure of ordering of the cations with respect to our chosen direction [001] pc , and intriguingly, it reaches its minimum value of around 0.4 at x MPB for all ferroelectrically active cations. This provides an unequivocal evidence that the system acquires a state with maximum structural instability at the MPB, where the piezoelectric coefficient d 33 reaches its maximum value, while the Curie temperature T c abruptly decreases.
Figure 4 demonstrates the mean values of the polar shifts of the cations with x as determined from the refined structural model, which can be directly related to the values of intrinsic polarization of the system. The shifts are comparable with the reported values based on experiments as well as theoretical calculations on several similar Pb and Bi containing solid solutions44,45,46. Apart from Mg, cations show a gradual increase in the average shifts with increasing x. However, this increase cannot account for the abrupt rise of the d 33 coefficient at the MPB. This is highly significant, because it indirectly underlines the critical influence of enhanced randomness of the shift-directions on the response functions driven by an external stimulus. Moreover, it reinstates the fact that the ferroelectric Curie temperature is heavily coupled to the microscopic strain of the system, and hardly influenced by the magnitudes of polar shifts of the cations. It is further interesting to note that the mean shifts become almost equal in magnitude for the A-site and the Ti cations on the onset of the MPB, which characterise competing A- and B-site driven ferroelectricity47, 48. The magnitudes of the displacements and the standard deviations for the Mg are unexpectedly large considering the chemical behaviour of the cation24, 46, 49. For low values of x however, bigger uncertainty can be expected in Mg-O distances, especially because the Mg-O peak sits on the shoulder between the negative Ti-O and the positive Pb/Bi-O peak (see Fig. 1a). However as the x increases, the respective mean shifts seem to become smaller gradually and therefore the values should be more reliable.
Raman scattering analysis
Figure 5 shows the composition-dependent Raman spectra for the xBMT-PT system with the corresponding assignment of different modes in complex perovskite-type ferroelectrics50, 51. It is evident that there is a gradual change in the spectra as a function of x and therefore, it is difficult to pin-point the MPB intuitively. However the detailed analysis of the individual spectra reveals distinct anomalies for several phonon modes. The lowest energy modes near 44 and 88 cm−1 (peaks 1 and 2), which are commonly dominated by A-site cationic vibrations50,51,52, soften around x = 0.4 − 0.5 < x MPB (see Fig. 6a). In particular peak 2, which is observed in undoped PbTiO3, shows a distinct minimum in the wavenumber vs x together with a maximum in FWHM vs x. This provides a clear evidence for a local-scale structural transformation driven by a rearrangement of heavy A-site cations18.
It is worth noting that in single crystals of PbTiO3, the Raman scattering near 88 cm−1 corresponds to the E(TO) soft mode. Polycrystalline ferroelectric solid solutions however possess oblique phonon states of mixed character. Hence, a direct ‘one-to-one’ assignment to a certain irreducible representation is not appropriate. Nevertheless, the x value at which peak 2 softens, matches the composition exactly at which σ seems to reach a level of saturation for all ferroelectrically-active cations (see Fig. 4), on the onset of an enhanced orientation disorder in polar shifts of the cations (see Fig. 3).
The Raman scattering at 44 cm−1 (peak 1) appears only in the presence of BMT. However, it cannot result exclusively from Bi vibrations developing a two-mode behaviour of the peak around 88 cm−1, because if that would be the case, then the intensities of both peaks would vary gradually with x 53. This is evidently not the case here (see Fig. 6a), and indicates that the new peak should be related to a different structural state of A-site cations rather than to a different chemical species. In fact, the intensity of peak 1 abruptly increases between x = 0.5 and x MPB, along with a sharp hike in the peak width. Altogether this affirms that the addition of BMT induces a distinct energetically favourable state in the A-site cation subsystem, which becomes dominant at x ≥ x MPB. Besides, larger peak widths at compositions x ≥ x MPB for both peak 1 and 2 indicate strong structural disorder, which could be equally conceived as a frustration in the coupling of coherent local distortions. This conforms wonderfully to the features of the orientation order parameter S (Fig. 3) and σ (Fig. 4b), revealed from the PDF analysis. The occurrence of the MPB is actually mirrored by the A-cation dynamics following the fact that the squared wavenumber difference reaches its minimum value (see the inset in Fig. 6a), i.e. the difference between the two energy states becomes minimal precisely at the MPB.
The phonon modes around 230 (peak 6) and 280 cm−1 (peak 7) are related to the B-site cationic vibrations50,51,52. These two modes evidently merge at x = 0.4, which could be related to the suppression of the prominent tetragonal distortion within the BO6 octahedra. It is important to note that at x = 0.4 the mean polar shifts of the A-site cations get to its maximum value but the Ti polar displacements continue to increase and become almost equal in magnitude with those of the A-site cations exactly at the MPB. This holds a favourable condition to constitute a strong coupling between the off-centred A- and B-site cations. Evidently, the phonon mode involving the BO3 vibrations against the A-site cation vibrations has a well pronounced minimum at x MPB (see Fig. 6c), which apparently drives the system to a phase transition.
The Raman peak near 350 cm−1 increases steadily in intensity with the increase in x, and it has even higher intensity than the B-cation vibration mode near 280 cm−1 for compositions above the MPB (see Fig. 4). Lattice-dynamics calculations for Pb-based complex perovskite-type oxides reveal that the Raman scattering near 350 cm−1 arises from a point phonon mode, which is silent (T2u symmetry) in the aristotype structure \(Pm\bar{3}m\) but may generate Raman intensity in a distorted double-perovskite structure54. The doubling of the unit cell may be induced by NaCl-type local chemical order at the B site and/or antiferrodistortive structural order, but only the latter may generate Raman activity of the phonon states near 350 cm−1. The inspection of the atomic vector displacements indicates that this mode is comprised of oxygen vibrations along the A-O bonds in the {111} pc planes, but it can be also thought as a rotation of the BO6 octahedra about the 〈111〉 pc directions (see inset in Fig. 6d). Consequently, the intensity of the Raman peak near 350 cm−1 is expected to increase with the development of the BO6 tilts. This was in fact detected by combined Raman scattering and neutron/synchrotron x-ray diffraction in a few relaxor ferroelectrics under high pressure55,56,57. Therefore, the gradual increase in the intensity of the phonon mode near 350 cm−1 with increasing x suggests that the dynamic BO6 tilting i.e. the antiferrodistortive ordering becomes significant for large values of x. A further crosscheck of this fact was found in the analysis of the RMC refined structural models where a gradual decrease of the B-O-B bond angle, which is typically used as a measure of the static tilts of the BO6 octahedra58, was recorded with the increase in x (see Supplementary Figure S5). Hence this ensuing development of local dynamic antiferrodistortive order upon doping can effectively interfere with the cationic polar shifts through suppressing the flexibility and the affinity to reorient under external field.
Summary and Implications
On the whole, our results provide an atomistic view of the development of composition-driven structural phase transition of a ferroelectric solid solution based on the perovskite structure with chemical disorder on both A and B sites, based on a combined analyses of the PDF and Raman scattering data. It is noted that the apparent change in the average structure with increasing content of BMT can be envisaged essentially as concurrent increase in local structural disorder in terms of off-centre static displacements of cations and randomness in their directions. However, it is the degree of stochasticity of the polar-shift directions (parameter S) which evolves systematically with composition and describes the development of structural instability, leading to the morphotropic phase boundary. But this parameter alone cannot explain the structure-property connections. We found that at MPB, all ferroelectrically active cations (A-site Pb2+ and Bi3+, B-site Ti4+) acquire very similar off-centre displacements, which together with the enhanced flexibility lead to a strong dynamic coupling between the polar shifts of the A- and B-site cations. This suggested that the combination of such structural instability and dynamic coupling is pivotal and an absence of any one of these factors may not bring about the expected boom in the properties, following the fact that the abrupt fall in the piezoelectric coefficient perfectly correlates with the diminished coupling factor for x > x MPB. These specific aspects of atomic-level structural correlations have not been conceived so far for a ferroelectric system, and should be however distinguished from the random orientation of the so-called polar nanoregions or local electric fields, commonly attributed to relaxor ferroelectrics. Our refined structural models did not suggest any notable clustering or chemical ordering at the A- or B-sites in support of that (See Supplementary Figure S7). The present random-direction model does conforms to the concept of polarization rotation or polarization extension but further put forward the necessity of the collaborative dynamic coupling effect, which should be an important part of the atomistic driving force for the enhanced properties often seen around an MPB.
Generally speaking, the proposed model featuring distinct static and dynamic characteristics might be applicable for a broad range of perovskite-based ferroelectric solid solutions especially with Pb and Bi where similar behaviour of the average structure have already been reported. However it entreats on a number of colligated issues, such as the distinction of the composition-driven phase boundaries from the temperature or pressure-induced phase boundaries, nature of the B-site chemical disorder (here Mg acts as a modifier of the coupling processes between the ferroelectrically active elements), combination of ferroelectric-antiferroelectric order (as in PZT), and especially the MPBs of Pb-free systems in order to develop efficient design rules.
Methods
Samples
Ceramic samples of xBMT-PT with x = 0.10, 0.20, 0.30, 0.40, 0.50, 0.63, 0.65, and 0.70 were prepared following the conventional solid-state synthesis route, details of which can be found elsewhere25. Room-temperature powder x-ray diffraction data (Stoe Stadi-MP powder X-ray diffractometer) were collected to verify the formation of a single perovskite phase. Electron microprobe analyses (wavelength-dispersive Cameca Microbeam SX100 SEM-system) averaging over 50 points from each compound were conducted to confirm the expected chemical compositions (see Supplementary Material Fig. S1). Commercially available powder of PbTO3 (Sigma Aldrich, purity ~99.9%) was used as a reference sample in Raman spectroscopic analyses.
Total neutron scattering and RMC modelling
Room-temperature neutron total scattering data were collected at the Nanoscale Ordered Materials Diffractometer (NOMAD) at the Spallation Neutron Source of Oak Ridge National Laboratory. NOMAD is a dedicated instrument for the total scattering experiments and allows to collect data for a wide range of reciprocal-space vector Q (=4π sinθ/λ), which is a necessary condition to produce reliable PDFs. For our measurements the Fourier transformations were done with Q max = 31.4 Å−1, which provided a real-space resolution of around 0.1 Å.
RMC modelling of the structure against the PDF data was performed using the RMCprofile package39. The starting models for different compositions were built using the structural parameters determined from prior Rietveld refinements of the neutron powder diffraction pattern. The modeling box size was approximately 54 × 54 × 54 Å3 and consisted of ~13000 atoms. There were 20 independent runs for each composition in order to have good statistics of the structural parameters. The DISCUS software59 was used to extract the various structural parameters from the refined models.
Raman scattering
Raman spectra were collected with a Horiba T64000 triple-grating spectrometer equipped with an Olympus B41 confocal microscope (50x objective), and a Symphony liquid-N2-cooled CCD detector. The spectra were recorded with a laser wavelength of 514.5 nm, on plate-shaped pellets in backscattering geometry, with a spectral resolution of ~2 cm−1, and peak-position precision of 0.35 cm−1. No polarization, orientation, and spatial dependence of the Raman spectra were detected. The measured spectra were temperature-reduced to account for the Bose-Einstein phonon occupation factor and fitted with pseudo-Voigt (PV) functions (PV = q * Lorentz + (1 − q) * Gauss) to determine the peak positions, full widths at half maximum (FWHMs), and integrated intensities. The criterion for the maximum number of fitted peaks was dI/I < 0.5 for all peaks, where I and dI are the calculated integrated intensity and the corresponding uncertainty, respectively60. In fact, for all compounds the achieved ratios dI/I were less than 0.25.
References
Guo, R. et al. Origin of the high piezoelectric response in PbZr(1−x)Ti x O0. Phys. Rev. Lett. 84, 5423–5426 (2000).
Noheda, B. et al. A monoclinic ferroelectric phase in the Pb(Zr x Ti1−x )O solid solution. Appl. Phys. Lett. 74, 2059–2061 (1999).
Ahart, M. et al. Origin of morphotropic phase boundaries in ferroelectrics. Nature 451, 545–549 (2008).
Egami, T., Teslic, S., Dmowski, W., Viehland, D. & Vakhrushev, S. Local atomic structure of relaxor ferroelectric solids determined by pulsed neutron and x-ray scattering. Ferroelectrics 199, 103–113 (1997).
Egami, T., Dmowski, W., Akbas, M. & Davies, P. K. Local structure and polarization in pb containing ferroelectric oxides. AIP Conference Proceedings 436, 1–10 (1998).
Grinberg, I., Cooper, V. & Rappe, A. Relationship between local structure and phase transitions of a disordered solid solution. Nature 419, 909–911 (2002).
George, A. M., Iñiguez, J. & Bellaiche, L. Effects of atomic short-range order on the properties of perovskite alloys in their morphotropic phase boundary. Phys. Rev. Lett. 91(1), 045504 (2003).
Schoenau, K. A. et al. Nanodomain structure of PbZr1−x Ti x O3 at its morphotropic phase boundary: Investigations from local to average structure. Phys. Rev. B 75(1), 184117 (2007).
Burkovsky, R. G. et al. Structural Heterogeneity and Diffuse Scattering in Morphotropic Lead Zirconate-Titanate Single Crystals. Phys. Rev. Lett. 109(1), 097603 (2012).
Goossens, D. J. Diffuse Scattering from Lead-Containing Ferroelectric Perovskite Oxides. ISRN Mater. Sci 2013, 1–17 (2013).
Egami, T. Local Structure of Ferroelectric Materials. Annu. Rev. Mater. Res. 37, 297–315 (2007).
Noheda, B. et al. Polarization rotation via a monoclinic phase in the piezoelectric 92% PbZn1/3Nb2/3O3-8% PbTiO3. Phys. Rev. Lett. 86, 3891–3894 (2001).
Jeong, I.-K., Lee, J. K. & Heffner, R. H. Local structural view on polarization rotation in relaxor ferroelectric (1-x)Pb(Zn1/3Nb2/3)O3-xPbTiO3. Appl. Phys. Lett. 92(1), 172911–172913 (2008).
Fu, H. & Cohen, R. Polarization rotation mechanism for ultrahigh electromechanical response in single-crystal piezoelectrics. Nature 403(1), 281–283 (2000).
Egami, T. Local Structure and Dynamics of Ferroelectric Solids. Struct Bond 124, 69–88 (2007).
Dmowski, W. et al. Local Lattice Dynamics and the Origin of the Relaxor Ferroelectric Behavior. Phys. Rev. Lett. 100(1), 137602 (2008).
Takenaka, H., Grinberg, I. & Rappe, A. Anisotropic Local Correlations and Dynamics in a Relaxor Ferroelectric. Phys. Rev. Lett. 110(1), 147602 (2013).
Husson, E. Raman spectroscopy applied to the study of phase transitions. Key Eng. Mater 155–156, 1–40 (1998).
Suchomel, M. R. & Davies, P. K. Predicting the position of the morphotropic phase boundary in high temperature PbTiO3-Bi(B′B′′)O3 based dielectric ceramics. J. Appl. Phys. 96, 4405–4410 (2004).
Grinberg, I., Suchomel, M. R., Davies, P. K. & Rappe, A. M. Predicting morphotropic phase boundary locations and transition temperatures in Pb- and Bi-based perovskite solid solutions from crystal chemical data and first-principles calculations. J. Appl. Phys. 98(1), 094111 (2005).
Eitel, R. E. et al. New High Temperature Morphotropic Phase Boundary Piezoelectrics Based on Bi(Me)O3–PbTiO3 Ceramics. Jpn. J. Appl. Phys. 40, 5999–6002 (2001).
Eitel, R. E., Randall, C. A., Shrout, T. R. & Park, S.-E. Preparation and Characterization of High Temperature Perovskite Ferroelectrics in the Solid-Solution (1−x)BiScO3–xPbTiO3. Jpn. J. Appl. Phys. 41, 2099–2104 (2002).
Randall, C. A. et al. Investigation of a high T c piezoelectric system: (1−x)Bi(Mg1/2Ti1/2)O3–(x)PbTiO3. J. Appl. Phys. 95, 3633–3639 (2004).
Khalyavin, D. D. et al. Crystal structure of metastable perovskite Bi(Mg1/2Ti1/2)O3: Bi-based structural analogue of antiferroelectric PbZrO3. Chem. Mater. 18, 5104–5110 (2006).
Chen, J., Tan, X., Jo, W. & Rödel, J. Temperature dependence of piezoelectric properties of high-T C Bi(Mg1/2Ti1/2)O3–PbTiO3. J. Appl. Phys. 106(1), 034109 (2009).
Zhang, Q., Li, Z., Li, F., Xu, Z. & Yao, X. Temperature Dependence of Dielectric/Piezoelectric Properties of (1−x)Bi(Mg1/2Ti1/2)O3–xPbTiO3 Ceramics with an MPB Composition. J. Am. Ceram. Soc 93, 3330–3334 (2010).
Upadhyay, A. & Singh, A. K. Grain size dependent phase stabilities and presence of a monoclinic (Pm) phase in the morphotropic phase boundary region of (1−x)Bi(Mg1/2Ti1/2)O3–xPbTiO3 piezoceramics. J. Appl. Phys. 117(1), 144102 (2015).
Baba-Kishi, K. Z., Welberry, T. R. & Withers, R. L. An electron diffraction and Monte Carlo simulation study of diffuse scattering in Pb(Zr,Ti)O3. Journal of Applied Crystallography 41, 930–938 (2008).
Ge, W. et al. Lead-free and lead-based abo3 perovskite relaxors with mixed-valence a-site and b-site disorder: Comparative neutron scattering structural study of (Na1/2Bi1/2)TiO3 and Pb(Mg1/3Nb2/3)O1/2. Phys. Rev. B 88(1), 174115 (2013).
Kopecký, M., Kub, J., Fábry, J. & Hlinka, J. Nanometer-range atomic order directly recovered from resonant diffuse scattering. Phys. Rev. B 93(1), 054202 (2016).
Egami, T. & Billinge, S. Underneath the Bragg Peaks Structural Analysis of Complex Materials (Pergamon Materials Series, 2012).
Jeong, I.-K. et al. Direct observation of the formation of polar nanoregions in Pb(Mg1/3Nb2/3)O3 using neutron pair distribution function analysis. Phys. Rev. Lett. 94(1), 147602 (2005).
Jeong, I.-K. Ferroelectric-relaxor crossover in Ba(Ti1−x Zr x )O3 studied using neutron total scattering measurements and reverse Monte Carlo modeling. Phys. Rev. B 81(1), 214119 (2010).
Levin, I., Krayzman, V. & Woicik, J. C. Local-structure origins of the sustained curie temperature in (Ba,Ca)TiO3 ferroelectrics. Appl. Phys. Lett. 102(1), 162906 (2013).
Keeble, D. S. et al. Bifurcated Polarization Rotation in Bismuth-Based Piezoelectrics. Adv. Funct. Mater. 23, 185–190 (2013).
Zhang, N. et al. The missing boundary in the phase diagram of PbZr1−x Ti x O3. Nat. Commun. 5(1), 5231–5239 (2014).
Datta, K., Richter, A., Göbbels, M., Keen, D. A. & Neder, R. B. Direct mapping of microscopic polarization in ferroelectric x(BiscO3) − (1 − x)(PbTiO3) throughout its morphotropic phase boundary. Phys. Rev. B 93(1), 064102 (2016).
McGreevy, R. L. Reverse monte carlo modelling. J. Phys. Cond. Matter 13, R877–R913 (2001).
Tucker, M. G., Dove, M. T. & Keen, D. A. Application of the reverse monte carlo method to crystalline materials. J Appl. Cryst. 34, 630–638 (2001).
Aksel, E. et al. Local atomic structure deviation from average structure of Na0.5Bi0.5TiO3: Combined x-ray and neutron total scattering study. Phys. Rev. B 87(1), 104113 (2013).
Li, J. Y., Rogan, R. C., Ustündag, E. & Bhattacharya, K. Domain switching in polycrystalline ferroelectric ceramics. Nature materials 4, 776–781 (2005).
Heitmann, A. A. & Rossetti, G. A. Thermodynamics of ferroelectric solid solutions with morphotropic phase boundaries. J. Am. Ceram. Soc. 97, 1661–1685 (2014).
Gramsbergen, E. F., Longa, L. & de Jeu, W. H. Landau Theory of the Nematic Isotropic Phase Transition. Physics Reports 135, 195–257 (1986).
Grinberg, I. & Rappe, A. M. Nonmonotonic T C trends in Bi-based ferroelectric perovskite solid solutions. Phys. Rev. Lett. 98(1), 037603 (2007).
Qi, T., Grinberg, I. & Rappe, A. M. Correlations between tetragonality, polarization, and ionic displacement in pbtio3-derived ferroelectric perovskite solid solutions. Phys. Rev. B 82, 134113 (2010).
Suewattana, M., Singh, D. J. & Limpijumnong, S. Crystal structure and cation off-centering in Bi(Mg1/2Ti1/2)O3. Phys. Rev. B 86(1), 064105 (2012).
Ghita, M., Fornari, M., Singh, D. & Halilov, S. Interplay between A-site and B-site driven instabilities in perovskites. Phys. Rev. B 72(1), 054114 (2005).
Shi, J., Grinberg, I., Wang, X. & Rappe, A. M. Atomic sublattice decomposition of piezoelectric response in tetragonal PbTiO3, BaTiO3, and KNbO3. Phys. Rev. B 89(1), 094105 (2014).
Miura, K., Kubota, M., Azuma, M. & Funakubo, H. Electronic and Structural Properties of BiZn0.5Ti0.5O3. Jpn. J. Appl. Phys. 48(1), 09KF05 (2009).
Welsch, A.-M. et al. Transformation processes in relaxor ferroelectric PbSc0.5Ta0.5O3 heavily doped with nb and sn. Z. Kristallogr. 226, 126–137 (2011).
Datta, K., Richter, A., Göbbels, M., Neder, R. B. & Mihailova, B. Mesoscopic-scale structure and dynamics near the morphotropic phase boundary of (1 − x)PbTiO3 − xBiScO3. Phys. Rev. B 92(1), 024107 (2015).
Jiang, K. et al. Relationship between negative thermal expansion and lattice dynamics in a tetragonal PbTiO3-Bi(Mg1/2Ti1/2)O3 perovskite single crystal. RSC Adv 6, 3159–3164 (2016).
Chang, I. & Mitra, S. Long wavelength optical phonons in mixed crystals. Adv. Phys 20, 359–404 (1971).
Mihailova, B., Bismayer, U., Guettler, B., Gospodinov, M. & Konstantinov, L. Local structure and dynamics in relaxor-ferroelectric PbSc1/2Nb1/2O3 and PbSc1/2Ta1/2O3 single crystals. J. Phys. Cond. Matter 14, 1091–1105 (2002).
Sani, A. et al. High-pressure phases in highly piezoelectric PbZr0.52Ti0.48O3. Phys. Rev. B 69(1), 020105 (2004).
Maier, B. J. et al. Structural state of relaxor ferroelectrics PbSc0.5Ta0.5O3 and PbSc0.5Nb0.5O3 at high pressures up to 30 GPa. Phys. Rev. B 84(1), 174104 (2011).
Waeselmann, N. et al. Pressure-induced structural transformations in pure and ru-doped 0.9PbZn1/3Nb2/3O3-0.1PbTiO3 near the morphotropic phase boundary. Phys. Rev. B 85(1), 014106 (2012).
Chong, S. Y. et al. Local Structure of a Pure Bi A Site Polar Perovskite Revealed by Pair Distribution Function Analysis and Reverse Monte Carlo Modeling: Correlated Off-Axis Displacements in a Rhombohedral Material. J. Am. Chem. Soc. 134, 5836–5849 (2012).
Neder, R. B. & Proffen, T. Diffuse scattering and defect structure simulations - A cook book using the program DISCUS (Oxford University Press, UK, 2007).
Datta, K., Richter, A., Göbbels, M., Neder, R. B. & Mihailova, B. Atomistic origin of huge response functions at the morphotropic phase boundary of (1 − x)Na0.5Bi0.5TiO3 − xBaTiO3. Phys. Rev. B 90(1), 064112 (2014).
Vanderbilt, D. & Cohen, M. H. Monoclinic and triclinic phases in higher-order devonshire theory. Phys. Rev. B 63(1), 094108 (2001).
Acknowledgements
Financial support by the Deutsche Forschungsgemeinschaft (MI 1127/8-1) is gratefully acknowledged. K.D. would like to thank Matt Tucker for several useful discussions on using the RMCprofile package. The macros created for the analysis of the structure files using the various available tool in the DISCUS package can be obtained from K.D. and R.B.N. A portion of this research at ORNL’s Spallation Neutron Source was sponsored by the Scientific User Facilities Division, Office of Basic Energy Sciences, US Department of Energy.
Author information
Authors and Affiliations
Contributions
K.D. and B.M. developed and planned the original project. J.C. provided all the ceramic samples. K.D. carried out the Raman measurements with the help from B.M. B.M. analysed the Raman data. J.N. conducted the neutron scattering experiments. K.D. analysed the PDF data in consultation with R.B.N. K.D. and B.M. wrote the manuscript. All authors reviewed the manuscript.
Corresponding author
Ethics declarations
Competing Interests
The authors declare that they have no competing interests.
Additional information
Publisher's note: Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Electronic supplementary material
Rights and permissions
This work is licensed under a Creative Commons Attribution 4.0 International License. The images or other third party material in this article are included in the article’s Creative Commons license, unless indicated otherwise in the credit line; if the material is not included under the Creative Commons license, users will need to obtain permission from the license holder to reproduce the material. To view a copy of this license, visit http://creativecommons.org/licenses/by/4.0/
About this article
Cite this article
Datta, K., Neder, R.B., Chen, J. et al. Atomic-level structural correlations across the morphotropic phase boundary of a ferroelectric solid solution: xBiMg1/2Ti1/2O3-(1 − x)PbTiO3 . Sci Rep 7, 471 (2017). https://doi.org/10.1038/s41598-017-00530-z
Received:
Accepted:
Published:
DOI: https://doi.org/10.1038/s41598-017-00530-z
This article is cited by
-
Structural, dielectric and Raman spectroscopy investigations of Ba1–xBi2x/3(Fe0.5Nb0.5)0.025Ti0.975O3 ceramics
Applied Physics A (2023)
-
Phase transformation in wet chemically synthesized Y2NiFeO6, and its magnetic and energy storage properties
Applied Physics A (2020)
-
Ethylene glycol/citric acid stabilized wet chemically synthesized Y2CoNiO6, and its structural, dielectric, magnetic and electrochemical behavior
Journal of Materials Science: Materials in Electronics (2020)
Comments
By submitting a comment you agree to abide by our Terms and Community Guidelines. If you find something abusive or that does not comply with our terms or guidelines please flag it as inappropriate.