Advances in Materials Science and Engineering
Volume 2020 (2020), Article ID 6345691, 7 pages
https://doi.org/10.1155/2020/6345691

Research Article

A Comparative Analysis of Asymmetric (BaTiO₃) /(BaZrO₃) Superlattices via X-Ray Diffraction and Raman Spectroscopy

O. A. Maslova,1,2† Yu. I. Yuzyuk,³ and S. A. Barannikova1,2,4

¹National Research Tomsk State University, Tomsk, Russia
²Institute of Strength Physics and Materials Science, Siberian Branch, Russian Academy of Sciences, Tomsk, Russia
³Southern Federal University, Rostov-on-Don, Russia
⁴Tomsk State University of Architecture and Building, Tomsk, Russia
†Deceased

Correspondence should be addressed to O. A. Maslova; o_maslova@rambler.ru

Received 23 February 2020; Revised 17 May 2020; Accepted 28 May 2020; Published 24 June 2020

Academic Editor: Fernando Rubio-Marcos

Copyright © 2020 O. A. Maslova et al. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

Abstract.

This work is aimed at studying asymmetric ((BT) /(BZ)) superlattices, grown by pulsed laser deposition onto (001) MgO substrates. The thicknesses of BT (ferroelectric) and BZ (paraelectric) layers were varied so that х ranged from 0 to 1 at a modulation period of about 80 Å. The films were 400 nm thick. The out-of-plane lattice parameters of constituents were assessed using X-ray diffraction. The lattice dynamic peculiarities of superlattices were probed via Raman spectroscopy; special attention is paid to the analysis of E(1TO) and A1(2TO) ferroelectric soft modes. A comparative analysis of data acquired via both experimental techniques reveals the enhancement of stress between BT and BZ layers with a decrease in symmetry from the tetragonal to a monoclinic phase due to strains induced by the lattice parameter mismatch between the constituents.

1. Introduction

Over the last 50 years, materials based on ferroelectric perovskites, such as artificial superlattices, have attracted attention of researchers due to their outstanding properties, such as high dielectric constants, high polarization, low dielectric losses, and high Curie temperature, which find application in various optoelectronic and microelectronic devices: nonvolatile random access memory devices, integrated systems, and tunable microwave elements [1–10]. In contradistinction to single-composition ferroelectric thin films, artificial superlattices (SLs) ensure high permittivity and great polarization in the relevant systems [3, 9]. The structures, where one constituent is ferroelectric BaTiO₃ (ВТ) and another is dielectric SrTiO₃ (ST) or PbTiO₃ (PT) or paraelectric BaZrO₃ (BZ), are among the most extensively studied objects. This is due to the fact that their physical characteristics are determined by epitaxial strain that is caused by the lattice parameter mismatch between constituents. Another important moment is the possibility to tune strain in heterostructures by varying the composition (stoichiometry) or the thickness of one of constituents composing a superlattice. Strains induced by the lattice parameter mismatch of the adjacent layers alter the ion positions and lattice vibrations, especially in low-frequency ferroelectric soft modes (E(1TO)) that are highly sensitive to ion displacements in thin films. Thus, the soft modes may furnish detailed information on misfit strains and, consequently, on the phase transitions in ferroelectric superlattices, which is essential in terms of the fundamental science and the practical use of these structures.

Earlier, we have studied the PLD grown superlattices onto MgO substrates, where one constituent was BaTiO₃ and another was a solid solution with Ba atoms partly substituted with Sr, that is, (BT/BST) [11–13]. In those systems, E(1TO) ferroelectric soft modes were monitored via Raman spectroscopy that is a powerful nondestructive tool for this purpose [2, 4, 9, 14, 15]. The E(1TO) frequency was found to depend to a large extent on the Ba/Sr content ratio due to the increase in its frequency by three times, while going from a pure BT film towards a pure ST film [12, 13]. This drastic variation in the low-frequency soft mode component was referred to stresses arising between the epitaxial constituent layers in the superlattice.

The present work is aimed at investigating the influence of layer thickness ratio in superlattices on the soft mode behavior by the example of asymmetric (BaTiO₃)/(BaZrO₃) ((BT)/(BZ)) superlattices prepared via pulsed laser deposition onto (001) MgO substrates, where La0.5Sr0.5CoO₃ oxides served as buffer layers. The microstructural characteristics (phase composition and lattice parameters) and lattice dynamics of these systems were experimentally probed via X-ray diffraction and Raman spectroscopy. X-ray diffraction is known to be one of the most reliable methods for studying epitaxial strain in superlattices [16–18], because it ensures information on the degree of structural perfection of films and allows one to determine the lattice parameters and to establish the orientation between the film and the substrate. Raman spectroscopy is a nondestructive and powerful method for lattice dynamics characterization in perovskites [2, 4, 9, 14, 15, 19] and allows one to detect even the finest atomic displacement in a system under consideration by tracking transformations in ferroelectric soft modes that are the main indicators of misfit strain in such systems.

2. Experimental

2.1. Synthesis

Thin-film ((BT) /(BZ)) superlattices (with a total thickness of 400 nm) were synthesized onto (001) MgO substrates via pulsed laser deposition using a KrF excimer laser operating at 248 nm. The 4 nm thick La0.5Sr0.5CoO₃ (LSCO, aLSCO = 3.805 Å)) buffer layers were applied under an oxygen partial pressure of 0.2 mbar onto the substrates preheated to a temperature of 750°C. The oxygen pressure upon the deposition of BT and BZ constituents was kept at a value of 0.1 mbar. The BT and BZ layers were alternately grown onto MgO substrates with a constant periodicity Λ of about 80 Å, and their thickness ratio was varied so that x was 0, 0.15, 0.3, 0.5, 0.7, 0.85, and 1.

2.2. X-Ray Diffraction Measurements

X-ray diffraction measurements were made at room temperature using a Shimadzu XRD 6000 X-ray diffractometer with CuKα₁ radiation in the θ−2θ scanning mode. The scan step was 2θ = 0.01 deg.

2.3. Raman Spectroscopy Measurements

Raman spectra were excited at room temperature by an Ar⁺ laser (λ = 514.5 nm) and recorded in the backscattering geometry with an inVia Renishaw Reflex spectrometer equipped with a coupled-charge device (CCD detector). The laser beam was focused on the film surface in a ∼2 µm spot using a 100x objective (the Leica microscope). The spectra were acquired in “side-view” (, , and ) backscattering geometries (hereinafter referred to as ZX, XX, and ZZ spectra) in obedience to the crystallographic axes of the MgO cubic substrate: (for details, see [18]). The obtained spectrograms were then corrected for the Bose–Einstein temperature factor.

3. Results and Discussion

3.1. X-Ray Diffraction

Figure 1(a) displays the XRD spectrograms of ((BT) /(BZ)) SLs with x from 0 to 1; here x = 0 means a pure BT film and x = 1 refers to a BZ film; both are deposited using a similar route to that for SLs (see the Experimental section). According to one series of θ−2θ peaks observed within the detection limits of the instrument, all samples represent a single-phase well-oriented structure with no extraneous phases. Moreover, as seen in Figure 1(a), the (002) peaks of BZ ﬁlm and MgO substrate match each other. Using the Bragg law enabled us to evaluate the out-of-plane lattice parameters for parent BT and BZ films. For BT film, the lattice parameter enriches 4.014 Å, being intermediate between the a (3.992 Å) and c (4.036 Å) lattice parameters of the bulk BT crystal. In turn, the c-axis parameter of BZ thin film is 4.236 Å, exceeding the a parameter value of the bulk BZ (4.192 Å). Thus, the BZ thin film is apparently under a compressive tetragonal distortion due to the lattice parameter mismatch between BZ and MgO substrate.

(a)

(b)

Figure 1: (a) θ−2θ XRD profiles of asymmetric (BaTiO₃)/(BaZrO₃) superlattices at different thicknesses of BT and BZ constituents. (b) XRD profile calculated for (BaTiO₃)/(BaZrO₃) SL with x = 0.7 ((BT)0.3Λ/(BZ)0.7Λ (black) in comparison with a measured XRD spectrogram (red). The variance between experiment and theory can be due to imperfect layer/layer interfaces that are neglected in the simulation model.

All SLs in this study exhibit satellite peaks being characteristic of modulated heterostructures [17, 20], which evolve with constituent thickness x. Since the angular distance between the peaks is a function of stacking periodicity , the latter was found for each SL to be ∼80 Å in obedience to the Bragg formula: sinθn = nλx/2Λ, where n is an integer and λx = 1.5406 Å is the radiation wavelength used in this work.

To calculate out-of-plane lattice parameters c₁ and c₂ of BT and BZ layers, respectively, the relationship Λ = n₁с₁ + n₂с₂ was used, where n₁ and n₂ are the numbers of unit cells in BT and BZ layers, respectively. For this, the simulations underlying the SL perfect structure have been implemented using Matlab software according to the algorithm described in detail in [17, 20]. The diffraction amplitude was determined in a specular geometry based on the interplanar spacings of the BT and BZ layers along with their individual thicknesses as adjustable parameters. One of the examples of the best simulation of the XRD patterns of SLs is given in in Figure 1(b) for (BT)0.3Λ/(BZ)0.7Λ superlattice. It is important to mention that some theoretical satellites deviate from relevant experimental peaks, which is owing to the influence of imperfect layer/layer interfaces that have not been considered in the simulations [17]. In this respect, Figure 2 displays the out-of-plane parameters calculated for BT and BZ constituents in comparison to the bulk BT and BZ and their single thin ﬁlm values. As found, the out-of-plane parameter of BT layers decreases from 4.014 Å (x = 0) to 3.820 Å (x = 0.85), evidencing that the in-plane tensile stress in the BT layers in the SL is more pronounced than in the BT single ﬁlm. This implies that in such a system the a-axis of BT layers becomes directed along the normal to the substrate [20].

Figure 2: Out-of-plane parameters of BT and BZ layers as a function of BZ thickness.

In turn, the out-of-plane parameter of BZ layers changes slightly and its average value is about 4.23 Å. However, the parameter value for BZ layers is significantly larger than the a-axis bulk parameter, indicating that BZ layers are under in-plane compressive stress induced by the lattice parameter mismatch between BT and BZ layers. This in-plane compression results in enhanced tetragonal distortion of the BZ layers due to the discrepancy from BT layers.

3.2. Raman Spectroscopy Data

Figure 3 displays the room-temperature ZX, ZZ, and XX polarized Raman spectra of BT/BZ SLs at various constituent thicknesses. The Raman spectra of a pure BT film coincide with those obtained earlier in the same backscattering geometries in works [12, 13], evidencing a broad-wing low-frequency band in ZX spectrum, which is associated with E(1TO) ferroelectric soft mode. By analogy with [12, 13], this band is overdamped, revealing a maximum at about 30–40 cm and a width of approximately 200 cm⁻¹. However, adding the BZ constituents with increasing their thickness leads to a drastic transformation in the E(1TO) soft mode, as shown by the arrows in Figure 3 for ZX spectrum.

(a)

(b)

(c)

Figure 3: Room-temperature Raman spectra of (BaTiO₃)/(BaZrO₃) superlattices, recorded in “side-view” backscattering geometries at various BZ layer thicknesses. The arrows show the transformation of E(1TO) soft mode component at various x values and the dotted line outlines the position of another soft mode (full-symmetrical) component denoted as A1(2TO).

Noticeable changes in the position, intensity, and width of this Raman component, that is, a decrease in its intensity and width along with a dramatic shift towards the higher frequencies, are observed starting with a small thickness of BZ layers (x = 0.15), becoming even more pronounced with a further increase in BZ thickness. Thus, being initially overdamped, the E(1TO) component transforms into an underdamped mode.

Unlike the BT crystal, XX and ZZ spectra of the BT film are identical, exhibiting the same interference gap at ∼180 cm⁻¹. Meanwhile, introducing BZ layers leads to pronounced differences between these two types of spectra, which are observed in SLs with the periodicity up to x = 0.85 ((BT)0.15Λ(BZ)0.85Λ). However, arriving at x = 1, which corresponds to a pure BZ film, extinguishes the Raman lines in all three spectra, revealing the tetragonal paraelectric phase D4h¹ that allows no polar Raman-active modes [21].

An interesting observation is the appearance of a doublet below 200 cm⁻¹ in ZX Raman spectra of samples with concentrations of 0.3 ≤ x ≤ 0.7. This doublet even manifests in ZZ spectra (seeming to leak from ZX spectra), along with a pronounced pair of peaks in a range of 300–400 cm⁻¹, but is absent in XX spectra. The computer simulation of ZX spectra using Lorentzians in OriginLab processing medium (see Figure 4) reveals that the components of the doublet are centered at ∼136 and ∼168 cm⁻¹. Similar spectral features within the same wavenumber range were observed by Chemarin et al. [22] in polydomain BaCexZrO₃ solid solutions, exposed to hydrostatic pressure. They were assigned to E modes originating from Ba–O stretching and O–Ba–O bending force constants [22]. In BaCexZrO₃ specimens with x = 0.8 and x = 1 this doublet was observed even under ambient conditions but disappeared at pressures above 21 GPa. As mentioned in [22], the doublet seems to be due to nanodomains and nanophases, which then tend to vanish with increasing pressure and finally merge into a homogeneous phase. This leads to the complete connection of the frontier atoms between nanophases and, respectively, to the disappearance of a low-frequency doublet. Based on this fact, it appears that our samples with the periodicity asymmetry of 0.3 ≤ x ≤ 0.85 may undergo a pressure-like effect. Firth of all, there is the superposition of the signals from BT and BZ constituents due to the strong lattice parameter mismatch between the adjacent layers. According to XRD results (Figure 2), BZ layers are being under in-plane compression at simultaneous tension of BT constituents, and their lattice parameters in SLs slightly exceed the appropriate value in unconstrained BZ. Such a lattice distortion makes triply degenerated modes of cubic BZ split into the Raman-active components, as seen in Figure 3, thus playing the identical role of hydrostatic pressure, described by Chemarin et al. [22], and indicating the structural phase transitions to the lower symmetry. The appearance of the identical doublet was also observed in symmetrical (BT)Λ/2/(BZ)Λ/2 SLs in study [23], whose stacking periodicity varied between 32 Å and 256 Å. Furthermore, a visual comparison of ZX spectra presented in [23] with our ZX spectra reveals that the spectrum of (BT)Λ/2/(BZ)Λ/2 SL with Λ = 256 Å is similar to that of the (BT) Λ/(BZ)xΛ SL with x = 0.3; the spectral profile of the (BT)Λ/2/(BZ)Λ/2 SL with Λ = 62 Å almost matches that of the (BT) Λ/(BZ)xΛ SL with x = 0.5, and the signal of the (BT)Λ/2/(BZ)Λ/2 SL with Λ = 32 Å is quasi-identical to that of the (BT) /(BZ)xΛ SL with x = 0.7. So, the stacking periodicity also plays a significant role in the emergence of the above doublet, and, as seen from our analysis, even maintaining the total periodicity at a constant value (in our case about 80 Å) at varying the layer thickness of constituents allows one to achieve the pressure-induced effect in asymmetric superlattices by analogy with works [22, 23]. The absence of this doublet in SLs with x = 0.15 and x = 0.85 is apparently due to a quasi-homogeneous phase in specimens. In the first case, there is the dominance of the thicker BT layers within the SL structure and introducing the thinner BZ layers causes a drastic upshift in the E(1TO) soft mode component due to the strain induced by the lattice parameter mismatch. In the second one, the presence of the much thinner BT layers alternating with the thicker BZ layers reduces the symmetry of the BZ film and alters its Raman spectrum, exciting the E(1TO) soft mode. On the other hand, the disappearance of a doublet below 200 cm⁻¹ in the asymmetric SL with x = 0.85 seems to be a sign of the phase transition owing to the above described pressure-like effect of constituents which undergo noticeable misfit strain due to a pronounced lattice parameter mismatch between them, as seen from XRD data (see Figure 2) for a relevant SL.

Figure 4: Example of computer simulation of ZX Raman spectrum using Lorentzians for a (BaTiO₃)0.7Λ/(BaZrO₃)0.3Λ SL. E1 and E2 are presumably the doublet components described by Chemarin et al. [22]. E(1TO) and A1(2TO) are the soft modes.

As for soft modes, the computer simulation of ZX spectra in obedience to the procedure [12] reveals that the E(1TO) frequency of (BaTiO₃)0.15Λ/(BaZrO₃)0.85Λ superlattice increases by almost seven times in comparison with a pure BT film, achieving 218 cm⁻¹. A too high low-frequency soft-mode position (above 200 cm⁻¹) was frequently found in BT crystals in the ordered low-temperature rhombohedral phase [24–26]. Meanwhile, the symmetry of epitaxial film does not allow the emergence of the rhombohedral phase in them; thus one can assume the monoclinic r-phase that is the closest to the rhombohedral one [23]. The E(1TO) parameters (full width at a half-maximum (FWHM) and frequency (υ)) are plotted in Figure 5 as the functions of BZ layer thickness. As seen in Figure 5, the most drastic transformation of the E(1TO) component occurs upon the transition from a pure BT film to a SL with x = 0.30. As seen, the thicker BZ layers are, the smaller the variances in E(1TO) frequency and width are. Another soft mode component denoted as A1(2TO) in Figure 3, which is a full-symmetry component thoroughly described in works [24–26], is found to remain unchanged with increasing the BZ layer thickness (Figure 5). It is apparently may be interpreted by the fact that strains caused by the lattice parameter mismatch between BT and BZ constituents exert strong influence on the positions of Ti⁴⁺ ions in the plane of the substrate, reducing considerably their disorder. Meanwhile, the asymmetric thicknesses of constituents seem to amplify strains between the layers, leading to reorientation of polarization vectors in constituents so that the polarization vector in BT layers, being initially along the с-axis, lies in the plane of the substrate, whereas the polarization in BZ layers is oriented along the normal to the substrate plane, which coincides with the above-described XRD results.

(a)

(b)

Figure 5: Frequency (υ, cm⁻¹) and width (FWHM, cm⁻¹) of Е(1ТО) and A1(2TO) ferroelectric soft modes versus BZ layer thickness x.

4. Conclusions

A series of ferroelectric/paraelectric (BaTiO₃)/(BaZrO₃) ((BT) /(BZ)) artificial superlattices with varying constituent thicknesses and a constant stacking periodicity was thoroughly probed using room-temperature XRD and Raman spectroscopy techniques. In obedience to the X-ray patterns (both experimental and simulated), the BT layers in the SLs undergo the in-plane tensile stress, at which the a-axis becomes directed along the normal to the substrate. In turn, the BZ layers overcome the in-plane compressive stress because of misfit strains induced by the lattice parameter mismatch between BT and BZ layers. This observation agrees with polarized Raman spectroscopy data. Moreover, according to the Raman spectroscopy, the strain in BT constituents seems to induce the monoclinic phase in them, decreasing to a large extent the disorder of Ti⁴⁺ ions, which suggests the appearance of a polar phase in BZ layers, as follows from the drastically upshifted (by seven times) low-frequency E(1TO) soft mode and the unchanged A1(2TO) soft-mode component.

Therefore, in comparison with epitaxial structures with varying stoichiometry in one of the constituent layers [11–13], the design of superlattices with asymmetric constituent thickness at a constant modulation period allows one to achieve a much more drastic transformation of E(1TO) soft modes (by seven times in this study) over against the parent BT film. In this respect, the fabrication of superlattices with a predefined composition enables one to achieve the structures with tunable properties, which ensures their proper application.

Data Availability

The data used to support the findings of this study are available from the corresponding author upon request.

Conflicts of Interest

The authors declare that there are no conflicts of interest associated with this work.

Acknowledgments

O. A. Maslova is grateful to Yu. I. Yuzyuk for help in formulation of task and experimental characterization of samples, as well as for the useful discussion of results. This work was supported in the framework of the Program of Competitiveness by the National Research Tomsk State University for the period of 2018–2020 and the Program of Fundamental Scientific Research of Russian State Academies of Sciences for 2013–2020, direction III.23.

References

F. Ahmadi and H. Araghi, “Coexistence of magnetism and ferroelectricity in (PbTiO₃)m∕(BaTiO₃)n superlattices,” Superlattices and Microstructures, vol. 140, Article ID 106427, 2020.

R. R. Das, Y. I. Yuzyuk, P. Bhattacharya, V. Gupta, and R. S. Katiyar, “Folded acoustic phonons and soft mode dynamics in BaTiO₃/SrTiO₃ superlattices,” Physical Review B, vol. 69, p. 132302, 2004.

H. Tabata, H. Tanaka, and T. Kawai, “Formation of artificial BaTiO₃/SrTiO₃superlattices using pulsed laser deposition and their dielectric properties,” Applied Physics Letters, vol. 65, no. 15, pp. 1970–1972, 1994.

N. K. Karan, R. S. Katiyar, T. Maiti, R. Guo, and A. S. Bhalla, “Raman spectral studies of Zr4⁺rich BaZrxTi1−xO₃ phase diagram,” Journal of Raman Spectroscopy, vol. 40, no. 4, pp. 370–375, 2009.

K. Wang, W. Jiang, J.-N. Chen, and J.-Q. Huang, “Study of the electronic structure and half-metallicity of CaMnO₃/BaTiO3 superlattice,” Superlattices and Microstructures, vol. 97, pp. 116–124, 2016.

E. Bousquet, M. Dawber, N. Stucki et al., “Improper ferroelectricity in perovskite oxide artificial superlattices,” Nature, vol. 452, no. 7188, pp. 732–736, 2008.

P. R. Choudhury and S. B. Krupanidhi, “Dielectric response of BaZrO₃/BaTiO₃ and SrTiO₃/BaZrO₃ superlattices,” Journal of Applied Physics, vol. 104, no. 11, p. 114105, 2008.

N. A. Pertsev, P. E. Janolin, J. M. Kiat, and Y. Uesu, “Enhancing permittivity of ferroelectric superlattices via composition tuning,” Physics Review B, vol. 81, Article ID 144118, 2010.

B. D. Qu, M. Evstigneev, D. J. Johnson, and R. H. Prince, “Dielectric properties of BaTiO₃/SrTiO₃ multilayered thin films prepared by pulsed laser deposition,” Applied Physics Letters, vol. 72, no. 11, pp. 1394–1396, 1998.

A. I. Lebedev, “Ab initio studies of dielectric, piezoelectric, and elastic properties of BaTiO₃/SrTiO₃ ferroelectric superlattices,” Physics of the Solid State, vol. 51, no. 11, pp. 2324–2333, 2009.

O. A. Maslova, Y. I. Yuzyuk, N. Ortega, A. Kumar, R. Katiyar, and S. A. Barannikova, “A comparative analysis of BaTiO₃/(Ba, Sr)TiO³ and BaTiO₃/(Ba, Sr)TiO₃/SrTiO₃ artificial superlattices via Raman spectroscopy,” Materials Research, vol. 22, Article ID e20180389, pp. 1–5, 2019.

O. A. Maslova, Y. I. Yuzyuk, N. Ortega, A. Kumar, S. A. Barannikova, and R. Katiyar, “Phase transition peculiarities in BaTiO₃-based perovskite superlattices,” AIP Conf. Proceed, vol. 2051, no. 1-4, Article ID 020190, 2018.

N. Ortega, A. Kumar, O. A. Maslova, Y. I. Yuzyuk, J. F. Scott, and R. S. Katiyar, “Effect of periodicity and composition in artificial BaTiO₃/(Ba, Sr)TiO₃ superlattices,” Physics Review B, vol. 83, p. 144108, 2011.

T. Shimuta, O. Nakagawara, T. Makino, S. Arai, H. Tabata, and T. Kawai, “Enhancement of remanent polarization in epitaxial BaTiO₃/SrTiO₃ superlattices with asymmetric structure,” Journal of Applied Physics, vol. 91, no. 4, pp. 2290–2294, 2002.

O. A. Maslova, Yu.I. Yuzyuk, and S. A. Barannikova, “Lattice dynamics of barium titanate: single crystal, ceramic, and polycrystalline film,” physica Status Solidi, 2020.

T. Harigai, S.-M. Nam, H. Kakemoto, S. Wada, K. Saito, and T. Tsurumi, “Structural and dielectric properties of perovskite-type artificial superlattices,” Thin Solid Films, vol. 509, no. 1-2, pp. 13–17, 2006.

O. Nakagawara, T. Shimuta, T. Makino, S. Arai, H. Tabata, and T. Kawai, “Dependence of dielectric and ferroelectric behaviors on growth orientation in epitaxial BaTiO₃/SrTiO₃ superlattices,” Vacuum, vol. 66, no. 3-4, pp. 397–401, 2002.

Yu.I. Yuzyuk, A. Almeida, M. R. Chaves, V. A. Alyoshin, I. N. Zakharchenko, and E. V. Sviridov, “Soft mode in heteroepitaxial (Ba, Sr)TiO₃MgO thin films,” Physica Status Solidi, vol. 222, p. 235, 2000.

J.-P. Locquet, D. Neerinck, L. Stockman, Y. Bruynseraede, and I. K. Schuller, “Discrete and continuous disorder in superlattices,” Physical Review B, vol. 39, no. 18, pp. 13338–13342, 1989.

F. Le Marrec, R. Farhi, M. El Marssi, J. L. Dellis, M. G. Karkut, and D. Ariosa, “Ferroelectric PbTiO3/BaTiO3superlattices: growth anomalies and confined modes,” Physical Review B, vol. 61, no. 10, pp. R6447–R6450, 2000.

O. Diéguez, K. M. Rabe, and D. Vanderbilt, “First-principles study of epitaxial strain in perovskites,” Physics Review B, vol. 72, Article ID 144101, 2005.

C. Chemarin, N. Rosman, T. Pagnier, and G. Lucazeau, “A high-pressure Raman study of mixed perovskites BaCexZr1−xO₃,” Journal of Solid State Chemistry, vol. 149, no. 2, pp. 298–307, 2000.

M. El Marssi, Y. Gagou, J. Belhadi, F. De Guerville, Y. I. Yuzyuk, and I. P. Raevski, “Ferroelectric BaTiO₃/BaZrO₃ superlattices: X-ray diffraction, raman spectroscopy, and polarization hysteresis loops,” Journal of Applied Physics, vol. 108, Article ID 084104, 2010.

A. Scalabrin, A. S. Chaves, D. S. Shim, and S. P. S. Porto, “Temperature dependence of the A1 and E optical phonons in BaTiO₃,” Physica Status Solidi (B), vol. 79, no. 2, pp. 731–742, 1977.

C. H. Perry and D. B. Hall, “Temperature dependence of the Raman spectrum of BaTiO₃,” Physics Review Letter, vol. 15, p. 17, 1965.

K. Laabidi, M. D. Fontana, and B. Jannot, “Underdamped soft phonon in orthorhombic BaTiO₃,” Solid State Communications, vol. 76, no. 6, pp. 765–768, 1990.