# Interface combinatorial pulsed laser deposition to enhance heterostructures functional properties

Jérôme Wolfman, Béatrice Negulescu, Antoine Ruyter, Ndioba Niang, Nazir Jaber

► To cite this version:

Jérôme Wolfman, Béatrice Negulescu, Antoine Ruyter, Ndioba Niang, Nazir Jaber. Interface combinatorial pulsed laser deposition to enhance heterostructures functional properties. Dr. Dongfang Yang. Practical Applications of Laser Ablation, , 2020, 978-1-83968-304-6. 10.5772/intechopen.94415. hal-03038828v3

HAL Id: hal-03038828

https://univ-tours.hal.science/hal-03038828v3

Submitted on 14 Feb 2022

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

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

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

Distributed under a Creative Commons Attribution 4.0 International License
## Chapter

# Interface Combinatorial Pulsed Laser Deposition to Enhance Heterostructures Functional Properties

Jérôme Wolfman, Beatrice Negulescu, Antoine Ruyter, Ndioba Niang and Nazir Jaber

## Abstract

In this chapter we will describe a new development of combinatorial pulsed laser deposition (CPLD) which targets the exploration of interface libraries. The idea is to modulate continuously the composition of interfaces on a few atomic layers in order to alter their functional properties. This unique combinatorial synthesis of interfaces is possible due to very specific PLD characteristics. The first one is its well-known ability for complex oxide stoichiometry transfer from the target to the film. The second one is the layer by layer control of thin film growth at the atomic level using in-situ RHEED characterization. The third one relates to the directionality of the ablated plume which allows for selective area deposition on the substrate using a mobile shadow-mask. However PLD also has some limitations and important PLD aspects to be considered for reliable CPLD are reviewed. Multiple examples regarding the control of interface magnetism in magnetic tunnel junctions and energy band and Schottky barrier height tuning in ferroelectric tunable capacitors are presented.

**Keywords:** combinatorial synthesis, oxide thin films and multilayers, functional heterostructure

## 1. Introduction

The first report of laser ablation generated plasma to deposit a film dates back to 1965, using a continuous ruby laser.[1] As the obtained film quality was inferior to those made with other deposition techniques at that time, the laser deposition remained confidential for two decades. The discovery mid 80's of the high Tc superconductor YBa₂Cu₃O₇-χ (YBCO) [2-3] triggered a frantic search for room temperature superconducting cuprates materials, including strong effort for thin film synthesis. The complex cuprate cationic composition makes film growth difficult using conventional physical vapor deposition (PVD) techniques like magnetron sputtering or evaporation. Indeed, to grow films with the right composition it is mandatory to start with a stoichiometric plasma. The Ar⁺ sputtering rate of multi-cationic targets and the distribution in the plasma strongly depends on the cations mass, which precludes the use of a stoichiometric target to
obtain a stoichiometric plasma. The offset target composition leading to stoichiometric films is unpredictable and a trial-and-error method is usually employed that renders magnetron sputtering an impractical technique for multicationic new materials research. Evaporation temperatures depend on the nature of the cations, which make a thermal vaporization of stoichiometric target inappropriate too. The attempt to produce cuprates films using sputtering and evaporation led to poor quality films.

Ceramic target vaporization based on laser ablation does not depend on the nature of the atoms within the target and result in a stoichiometric plasma as long as the energy per surface unit (i.e. the laser fluence) is above the ablation threshold. Venkatesan and co-workers were the first to vaporize an YBCO target using a pulsed excimer laser (UV KrF). After fine tuning the ablation conditions by checking the spatial distribution of the cations they rapidly produced superconducting films having the best physical properties of that time. [4] This first successful synthesis of high Tc cuprates films using a pulsed laser drew the interest of the community and several laboratories started working in the field now known as pulsed laser deposition (PLD). In the following decade, beside cuprates superconductors, strong interest in colossal magneto-resistive manganites and multiferroic ferrites has arisen. This intense scientific activity on multi-cation oxides associated to PLD ease of use, versatility and reasonable cost contributed to its rapid expansion in the 90's

Since then PLD based thin film research thrived over a wide range of materials, physical properties and applications. PLD has benefited from in-situ real-time characterization tools developed for other deposition technique to mature into an epitaxial film growth method with control at the atomic scale and able to produce heterostructures with sharp interfaces. At the end of the 90's, material scientists considered transposing combinatorial synthesis, a high throughput synthesis method developed by the pharmaceutical industry, to the field of multi-cation oxides research. Combinatorial PLD (CPLD) rapidly emerged, the basic idea being to grow on a single substrate, in a single deposition run, a film with continuous in-plane chemical modulation. In this way, a wide range of chemical compositions are produced within the same sample and can be quickly scanned to identify compounds with optimum targeted properties. Since the new millennium, CPLD has been refined and its field extended to ternary phase diagram exploration. More recently CPLD research field has been extended to a new territory: the exploration of interface compositions in heterostructures with enhanced functional properties a.k.a. Interface Combinatorial Pulsed Laser Deposition (ICPLD)

## 2. Pulsed laser deposition for combinatorial synthesis

The aim of this chapter is not to thoroughly describe PLD but to emphasize its most important aspects, advantages and limitations with regard to combinatorial synthesis of oxide films and heterostructures. The reader interested in an exhaustive description of PLD is encouraged to consult PLD's introduction reference book [5] or the other chapters devoted to the subject in this book.

The first pre-requisite for CPLD is the growth of films with uniform thickness and homogeneous composition over the entire sample surface. Although PLD has a reputation for stoichiometric transfer from the target to the film, this is not however straightforward and several deposition parameters have to be fine-tuned. Starting from the formation of a stoichiometric plasma, target thermal vaporization occurring below the ablation threshold should be reduced as much as possible. This means that local fluence everywhere on the laser beam spot should be above the ablation threshold, implying very steep sidewalls of the laser beam energy distribution. Such a distribution, called a top-hat, necessitate laser beam
shaping with a beam homogenizer. The laser beam divergence depends on the discharge voltage and could affect the beam spot size and energy distribution depending on the beam shaping method, so the discharge voltage should be kept constant. The plasma expand from the target toward the substrate in the form of a plume which interacts with ambient gas molecules or atoms. Multiple collisions per atom or ions occur leading to a thermalized but still highly directional plasma reaching the substrate. It follows a radial distribution of thickness and composition at the substrate surface. In order to get a uniform and homogeneous film, one has to scan the plume with respect to the substrate, by either scanning the laser spot over the target surface at constant focus and/or moving the substrate (rotation, translation). Doing so will result in a film with much increased thickness uniformity and homogeneous composition. However the composition may not match the target stoichiometry. Indeed, it has been shown that cationic off-stoichiometry up to a few percent occurs depending on the fluence for a fixed spot size and target-substrate distance [6]. So cationic composition/uniformity should always be checked versus fluence with high sensitivity characterization techniques. Once the right fluence has been determined and stabilized over the laser spot on the target, it is crucial to ensure that it stays constant in time. Regarding this point, two factors have to be taken into account. The first one is the degradation of the excimer laser gas charge with time, leading to a drop of the output energy. As already specified, it is advisable to work at constant voltage discharge, so one should not compensate the energy drop with voltage. The solution is to start with a higher energy output than required and to modulate the beam energy with an external attenuator. The second factor is the laser attenuation at the chamber introduction viewport. Although the visible part of the plasma, i.e. the plume, is very forward peaked, some low energetic species are slowly deposited everywhere in the chamber, including on the introduction viewport. This leads to a time varying attenuation which will affect the fluence at the target surface. To compensate for this attenuation, it is mandatory to measure the fluence after the viewport and to increase the input energy accordingly. One can reorient the beam out of the chamber through a pollution preserved port with a translatable UV mirror to measure its energy without breaking the vacuum.

The second pre-requisite for CPLD is a layer by layer growth mode for a control of the composition at the unit cell level. Layer by layer PLD growth of various perovskites has been demonstrated multiple times using in-situ real time grazing incidence electron diffraction (RHEED). This ability for 2d growth comes from the very peculiar and sequential PLD surface crystallization kinetic. Indeed, although average PLD deposition rate is quite slow, instantaneous deposition rate a few µs after laser pulses is extremely high, creating a supersaturation close to the surface which results in high nucleation rate. The large number of nucleus present after the first pulse favors 2d growth. Furthermore, as PLD relies on photon to vaporize the target, the ambient gas pressure can be varied over a very wide range from vacuum level to a few mbar. The deposition pressure allows for the control of the energy of the species reaching the substrate and ultimately their remaining kinetic energy to explore the surface and find the nucleation sites.

The last pre-requisite for CPLD is the ability to lower the average deposition rate per pulse to produce the smallest composition step when mixing materials from different targets. For a fixed fluence and target-substrate distance, the deposition rate per pulse strongly correlates with the laser spot size on the target which is easily adjustable. Deposition rate as low as a few hundred pulses per perovskite unit cell thick layer can be reached. One drawback of PLD is the lack of deposition rate stability. Indeed the structure of the irradiated target surface evolves with laser exposure; the deposition rate lowering with target aging is more severe for a brand new surface and decays with target exposure. Scanning over a
large area of the target reduces the number of laser spots per location and makes this problem less stringent. It is however mandatory to calibrate the deposition rate prior a deposition or to measure it in-situ using RHEED.

**3. Combinatorial Pulsed Laser Deposition**

The concept of combinatorial research, initially introduced in the pharmaceutical industry, aims at synthesizing a large number of compounds with complex and systematically varied composition in a single batch. A high speed characterization technique is then used to scan through this material library and identify the compounds presenting the targeted property.

This approach is very appealing regarding multi-cationic perovskite-related oxides. Indeed their physical properties, covering a very wide span (superconductivity, ferroelectricity, ferromagnetism, colossal magnetoresistance, tunable resistivity...), result from the subtle equilibrium between competing interactions involving charges, spins and orbitals. This results in an extreme sensitivity of the properties to the cationic composition and requires a fine, thorough and systematic scanning in order to optimize the performances.

The first thin film combinatorial synthesis attempts implied multiple successive room temperature depositions to vary the composition. Several annealing steps were then required for thermal diffusion and crystallization. This approach, where reaction products are dominated by thermodynamic like in ceramic sintering, is not appropriate for epitaxial thin films and heterostructures. On the contrary, epitaxial oxide thin films growth by PLD does not require post-deposition annealing. Moreover the directionality of the plume allows for localized deposition on the substrate through a shadow-mask. These multiple advantages were soon associated to combinatorial synthesis of oxide materials at the end of the 90's simultaneously in the United States and in Japan. The original approach is schematically described in **Figure 1**.

The single unit cell combinatorial cycle consists of two successive depositions from different targets. The material deposited from each target is distributed across the substrate surface using a moving shadow-mask in order to vary locally the layer completeness (from 0% to 100%). After one deposition cycle, the resulting 1µc thick layer has a composition which varies laterally from one target composition to the other (composition spread). Then this deposition cycle is repeated N times in order to achieve targeted final film thickness.

This idealized combinatorial PLD synthesis produces a continuous variation of compositions on a single sample, guarantying identical growth conditions for all compositions and eliminating the risk of sample variability.

<figure><img src="image_2.png" /><figcaption>**Figure 1.**
Schematic of the combinatorial single layer synthesis and of the resulting composition spread library (e.g. Bi₁-xGaxFeO₃ composition library)</figcaption></figure>
A consequence of CPLD is the necessity to use local probes to scan through the compositions library and assess performances. There is a trade-off between the minimum probe size required to measure the targeted physical property and the lateral composition gradient in order to characterize a “homogeneous” compound at the probe level.

The main difficulty regarding CPLD synthesis is to reach adequate control of both local composition and thickness using PLD. As explained in the previous chapter this requires long in-depth preparation work involving chemical characterization and exhaustive optimization of all the deposition parameters.

A prerequisite to reliable CPLD synthesis is to produce a uniform and smooth film with constant composition and thickness from each target over the surface of the future CPLD samples. One should not attempt CPLD before proving this achievement. Another important point is to keep a statistical approach to the characterization of the libraries. It is tempting to produce ternary phase diagrams using 3 different targets on a single substrate. However, in this case there is only one location on the sample per compound which is statistically insufficient. To our knowledge, any published work concerning CPLD ternary phase diagrams on a single sample does not present a proper characterization that demonstrates the control of the composition and thickness across the library. This is not very surprising knowing how difficult it already is to master a single gradient for a binary diagram across one sample. Unfortunately even for binary samples a large number of CPLD articles have been published without any evidence of control of composition/thickness which brought some discredit on CPLD among the scientific community. To change this perception and reinforce CPLD credit, we proposed a decade ago a statistical approach to the characterization of both the composition and the physical properties for binary phase diagram. [7] More recently we developed an alternative approach to explore ternary phase diagrams: instead of trying to produce the full ternary phase diagram on a single sample, we select lines of compositions cutting through the corresponding triangular diagram. Each synthesized sample has a composition gradient in one direction with the compositions defined for one line. [8] This way statistical characterization are possible along the direction orthogonal to the gradient. Scanning through the ternary phase diagram along multiple lines thus requires the synthesis of a few samples with only three targets.

To illustrate the effectiveness of CPLD we will discuss about the search for new lead-free piezoelectrics to replace (1-x)PbZrO₃-xPbTiO₃ (PZT), the most used material in microelectromechanical systems for sensing, actuating, or energy harvesting applications. PZT ferroelectric films present large piezoelectric coefficients and electromechanical coupling, enabling long range motions and high energy densities. [9] A unique characteristic of lead-based solid solutions presenting high piezoelectric coefficients is the strong enhancement of their piezoelectric response in the vicinity of a composition induced phase transition between ferroelectric phases with different crystalline symmetries, called a morphotropic phase boundary (MPB). In PZT the MPB lies between a rhombohedral ferroelectric phase and a tetragonal ferroelectric phase. [10] The microscopic origin of this enhanced piezoelectric activity is still being debated but usually involves easiness of polarization rotation at the MPBs. [11-13] PZT is lead-based and thus targeted by environmental regulations (e.g. RoHS EU Directive). Thus, alternative piezoelectric lead-free materials are required, and an obvious direction is to look for MPBs in other ferroelectric solid solutions. The rhombohedral perovskite BiFeO₃ (BFO), being a robust ferroelectric (Tc~1100K) with record polarization (100 µC·cm⁻²), is a good starting point. [14, 15] Solid solutions of BFO with tetragonal ferroelectric perovskite phase like PbTiO₃ and BaTiO₃ have been synthesized, and MPBs have been found in both cases. [16, 17].
As previously discussed, the first step is to produce uniform films of controlled composition before aiming to CPLD. The Bi element being volatile, the stabilization of pure BFO in thin films is not straightforward, as several parasitic phases can coexist. To compensate for Bi volatility we used an enriched Bi<sub>1.1</sub>FeO<sub>3</sub> target and found the range of temperature, oxygen pressure and fluence which lead to pure BFO films. The structure was studied by x-ray micro-diffraction (μXRD) and Bi/Fe ratio by Rutherford Back Scattering (RBS) versus laser fluence. A fluence of 1.72 J·cm<sup>-2</sup> was identified as giving Bi/Fe = 1 with a deposition temperature of 700 °C, an oxygen pressure of 0.2 mbar, a laser repetition rate of 6 Hz with a target-substrate distance of 4.5 cm. [18] RBS is an averaging technique (spot size of 2 x 2 mm<sup>2</sup>), so to assess local uniformity of thickness and composition we turned to an Electron Probe Micro-Analyzer (EPMA) equipped for Wavelength Dispersive X-Ray Spectroscopy (WDS) (see [18] for details). A specific thin film analysis program has been used to determine BFO's composition (TFA/WDS layerf, Cameca). BFO film's density × thickness product (ρ·t) and composition were simultaneously computed. Thirty measurements were realized with a beam diameter of 20 µm (20 keV, 100 nA) every 300 nm along the film. **Figure 2a** represents the ρ·t product and the weight percentages of Bi, Fe and O after self-consistent analysis of the raw data. The average value of ρ·t is 242 µg/cm<sup>2</sup> with a standard deviation of 5.4 µg/cm<sup>2</sup> equivalent to a relative variation of 2.2%. Considering the bulk BFO density (d = 8.38 g/cm<sup>3</sup>) we find a thickness t = 289 nm ± 6.5 nm (1σ). The stability of composition is even greater along the sample. The average weight percentages and standard deviations for the Bi and Fe, transposed into atomic percentages are respectively 20.02% ± 0.08% and 19.98% ± 0.08% for an expected value of 20%. With these statistical analysis, we find that the composition dispersion is Bi<sub>1.001</sub>±<sub>0.004</sub>Fe<sub>0.999</sub>±<sub>0.004</sub>O<sub>3</sub> (1σ) and thickness standard deviation σ(t)≤2.2% i.e. a very good thickness and composition uniformity along the sample surface. As the aim was to measure the piezoelectric coefficient of BFO-based solid solution, we deposited BFO onto La<sub>0.8</sub>Sr<sub>0.2</sub>MnO<sub>3</sub> (LSMO) an epitaxial conductive oxide electrode grown on SrTiO<sub>3</sub> (001) substrate. A zoom of the θ-2θ X-ray diffraction pattern around (001) diffraction peak is shown in **Figure 2b**. Only (00l)pc (pseudo-cubic notation) oriented diffraction peaks are visible on this pattern while no parasitic phase could be detected [19]. The thickness (Pendellösung) fringes observed around both LSMO and BFO (001) reflections demonstrate the crystalline quality and the smoothness of surface and interfaces.

We choose GaFeO₃ (GFO) as the second member of the solid solution to be explored in order to find a MPB. GFO does not have a perovskite structure but crystallizes in a much more complex orthorhombic structure (SG Pc21n). [20] In

<figure><img src="image_3.png" /><figcaption>**Figure 2.**</figcaption></figure>

(a) Evolutions of density × thickness product (ρ·t) and weight percentages of Bi, Fe and O as a function of the position and (b) X-Ray diffraction pattern of BFO/LSMO/STO heterostructure
the (1-x)BFO - (x)GFO system (BGFO) designed here, Ga<sup>3+</sup> cations have incentives to substitute for Bi<sup>3+</sup>. Small ions like Ga<sup>3+</sup> do not occupy perovskite A-site in conventionally synthesized ceramics. However, perovskites with small A-site cations like Sc<sup>3+</sup> and Mg<sup>2+</sup> have recently been stabilized using high-pressure high-temperature synthesis, and Ga<sup>3+</sup> is envisaged as A-site cation in this emerging field.[21] Epitaxial strain during PLD growth has longtime proven to be an alternative to high pressure synthesis for metastable phase stabilization. [22] As GFO is not a perovskite, we do not expect to obtain a solid solution at high x values, so we limited our range of investigation to 0 ≤ x ≤ 0.12. In this range, the Goldsmith tolerance factor t is greater than 0.87, not too far from non-substituted BFO's tolerance factor (t = 0.89), and compatible with a distorted perovskite structure. So, it is plausible that some Ga<sup>3+</sup> ions occupy the perovskite A-site in our films, although it is probable too that part of the Ga<sup>3+</sup> shares the B-site together with Fe<sup>3+</sup>. The Ga<sup>3+</sup> substitution for Bi<sup>3+</sup> being limited to 0%-12%, the BFO deposition conditions were used for GFO. WDS analysis were done along and across the composition gradient (1 point each 300 µm). Assuming the Bi<sub>w</sub>Ga<sub>x</sub>Fe<sub>y</sub>O<sub>z</sub> formula with z=3, cationic contents w, x, and y have been extracted from these measurements, the error being estimated to 0.005. The extracted Ga content is plotted versus position in **Figure 3a**, showing a linear increase from 0% up to 12%, in good agreement with the nominal concentrations. We note that both Bi and Fe contents decrease from 1 in pure BiFeO<sub>3</sub> to about 0.97 for 6% Ga doping. As x increases from 6% to 12%, Fe content goes back to about 1.0 while Bi content continue to decrease. [19] From these values, one could suspect that Ga is substituted partly for Bi and partly for Fe.

X-ray Reciprocal Space Mapping (RSM) around the (103)pc reflection have shown that BGFO and LSMO are epitaxial on STO (not shown here, see [19]). The lattice parameters evolutions confirm that Ga effectively enters into the BFO structure. Furthermore a characteristic splitting of the (103)pc reflection of BFO strained by cubic STO disappears for a Ga content 5% < x < 7%. This indicates a change of symmetry of the film and might be the signature of a MPB. Piezoelectric characterizations were made on 30 x 30 µm² top Pt electrodes (dc-sputter deposited via a lift-off process) using a laser scanning vibrometer (LSV model MSA-500, Polytec, V<sub>AC</sub>=1V) (see schematic of the heterostructure **Figure 3b** top). A typical mapping of the extracted effective piezoelectric coefficient $d_{33}^{eff}$ across one electrode, using a 3µm laser spot size, is presented in **Figure 3b** (bottom), showing a uniform displacement. The $d_{33}^{eff}$ coefficients were extracted from fifteen LSV measurements on each electrode and three different lines across the composition gradient were measured. **Figure 3c** shows the variation of the $d_{33}^{eff}$ as a function of Ga content with standard deviation represented as error bars. After a slow increase of $d_{33}^{eff}$ at low Ga content, a sharp peak centered at about 6.5% is

<figure><img src="image_4.png" /><figcaption>**Figure 3.**
(a) Ga concentration measured by WDS along the composition gradient. (b) Heterostructure schematic (top) and $d_{33}^{eff}$ mapping of a top electrode measured by laser vibrometry (bottom). (c) BGFO $d_{33}^{eff}$ as a function of Ga content. From [18] with permission</figcaption></figure>
<figure><img src="image_5.png" /><figcaption>**Figure 4.**
(a) Ferroelectric cycles for various Ga contents measured at 77K. The corresponding coercive fields, remanent and saturation polarizations and electric field offsets are plotted in (b), (c) and (d) respectively</figcaption></figure>

observed. The maximum $d_{33}^{eff}$ value, about 53 pm/V, is twice larger than the value obtained for undoped BFO. Finally, above x=8%, $d_{33}^{eff}$ falls to 15 pm/V. The sharp enhancement of piezoelectric properties occurring around 6.5% of Ga doping is correlated with the symmetry change observed around the same composition by RSM. [19] To confirm the presence of a MPB we looked for a change in the ferroelectric properties of BGFO with Ga content x. We used a ferroelectric tester (Radiant LC II) to measure polarization hysteresis P(E). The hysteresis cycles presented **Figure 4** were acquired at liquid nitrogen temperature.

A clear transition in the ferroelectric cycle shapes is visible **Figure 4a** as Ga content increases. As the Ga content gets over 6%, a strong increase of the coercive field (**Figure 4b**) associated to a strong decrease of both saturation and remnant polarization (**Figure 4c**) are observed, together with a change of sign of the electric field offset (hysteresis imprint **Figure 4d**). This demonstrates that a change of ferroelectric phase is occurring, correlated to the $d_{33}^{eff}$ peak, and implies that a MPB is present in BGFO at about 6.4%. It is important to note that the $d_{33}^{eff}$ peak is very sharp in composition and that a conventional ceramic approach with 1% doping steps would have miss it, emphasizing the power of the continuous composition spread in CPLD.

# 4. Interface Combinatorial Pulsed Laser Deposition

## 4.1 Interfaces of oxide heterostructures: the new territory

In strongly correlated complex oxides, charge, spin, orbital and lattice degrees of freedom co-exist and interplay cooperatively. In particular the complex balance between these degrees of freedom and related interactions generates a rich spectrum of competing phases in perovskites or perovskite-derived materials (e.g. high Tc superconductor, metal-insulator transitions, magnetism, ferroelectricity, piezoelectricity...). The recent progress in deposition techniques allowed the production of complex perovskites heterostructures with atomically sharp
inter interfaces, which expanded material researcher's horizon. Fascinating phenomena and novel states of matter at complex oxide heterointerfaces have been reported. One can cite for instance the existence of high mobility 2d electron gas at LaAlO3 / SrTiO3 interface, even becoming superconducting at low temperature, while both materials taken separately are insulating. [23, 24] Another striking example is the transition of CaTiO3 from its usual non-polar state into a high-temperature polar oxide thanks to interfacial tilt epitaxy [25]. The isomorphism of the ABO3 perovskite oxide structure allows for a wide range of chemically modulated interfaces.

Some of the phenomenon occurring at perovskite interfaces are reported **Figure 5**. Rumping, polar discontinuity, interfacial B-site cation environment asymmetry, BO₆ octahedral rotations are all potential levers to modulate interface properties. Their complex interplay is strongly affected by cationic substitutions and a complete and fine exploration of the possible interface compositions is required in order to identify new physics phenomena or enhanced properties. Interface CPLD (i.e. ICPLD) is a powerful tool in that respect.

One application where oxide interfaces plays a crucial role is the ferroelectric (FE) voltage tunable capacitor envisaged for future RF communication technologies (5G and Near Field Communication NFC). [27] The relative dielectric permittivity ε<sub>r</sub>(E) of a FE has a large electric field dependence. [27] The perovskite solid solution Ba<sub>1-x</sub>Sr<sub>x</sub>TiO<sub>3</sub> (BST) is the most widely used FE in current 4G thin film parallel plates varactors because of its excellent tunability/losses compromise. New specifications for 5G and NFC (higher frequency and reduced driving voltage) call for improved varactor properties. Reducing the FE film thickness from 240nm (4G) down to 50-100 nm range is one option toward meeting these new specifications. However, in such thin FE films the metal / FE (M/FE) interface influence is reinforced in a damaging way, due to two interfacial phenomena. The first one is the existence of FE “dead-layers” with degraded ε<sub>r</sub> and spontaneous polarization close to the electrodes, producing an effective non-tunable interfacial capacitance. [28, 29] The second one is the increased leakage current due to insufficient Schottky Barrier Height (SBH).

Interface engineering can be used to tailor band alignment and interface polarizability. The insertion of a thin layer with different atomic element(s) at the

<figure><img src="image_6.png" /><figcaption>**Figure 5.**
Schematic of competing interactions and phenomenon at perovskite ABO₃ interfaces.
From [26] with permission.</figcaption></figure>
interface allows to manipulate the chemical bonding and promotes atomic rearrangement. Let's consider for instance the anti-displacement of anions and cations predicted at $Ba^{2+}O^{2-}/M$ and $Sr^{2+}O^{2-}/M$ interfaces and quantified by a rumpling parameter R. [30] R depends on the chemical bonding and is responsible for an interface dipole, which in turn modulate the SBH. Interestingly, the insertion of e.g. a single Al atomic plane at the $BaO/M$ interface strongly affects R and SBH. Indeed, for M = Pd the SBH goes from 1.4eV to 2.6eV. [30] Significant rumpling has been experimentally shown for $SrTiO_3$ (STO) in contact with $La_{2/3}Sr_{1/3}MnO_3$ (LSMO), a metallic perovskite electrode, inducing a polarization in the non-ferroelectric STO. [31] The continuity of the perovskite structure through the LSMO/STO interface and its ionic character offer new ways to control electronic properties. In $La_{1-x}Sr_xMnO_3$ (LSMOx), the B-site cation ratio $Mn^{3+}/Mn^{4+}$ is determined by the A-site ratio $La^{3+}/Sr^{2+}$. Along [100], successive AO and $BO_2$ planes are polar for LSMOx and charge neutral for BST. Interfacing LSMOx with BST leads to tunable interfacial polar discontinuity which can induce lattice polar distortion and result in SBH modulation. [32-34]

LSMO is a ferromagnetic (FM) half-metal, i.e. having a 100% spin-polarization at the Fermi level. For the latter reason it has been intensively studied as a spin-polarized electrode in LSMO/STO/LSMO magnetic tunnel junction (MTJ). MTJs are used e.g. as memory bits in magnetic MRAMs. The tunnel resistance depends on electrode spin-polarization and on the relative orientation of the electrode magnetic moments, with high resistance $R_{AP}$ (resp. low resistance $R_P$) for antiparallel (resp. parallel) states. A 100% spin polarized electrode leads to a theoretical infinite $R_{AP}$ which is ideal for the cited application. In LSMO/STO/LSMO, a record tunnel magneto-resistance ($TMR= (R_{AP}-R_P)/R_P$) of about 2000% was reported, but unfortunately for temperature far below the Curie temperature $T_C$. [35] The vast majority of the electrons tunnel from the interfaces, their spin-polarization being affected by the nature of the chemical bonding. FM correlations at manganite interfaces are known to be weaker than in bulk, causing a magnetic “dead layer” which probably explains the diminution of TMR close to $T_C$. [36-38] Attempts have been reported at creating a doping profile at the interfaces by inserting a 2 uc thick $LaMnO_3$ layer [39, 40] or a single uc thick $La_{0.33}Sr_{0.67}MnO_3$ [41] layer to overcome this problem with some improvement of interface magnetism but still not a full recovery of bulk properties. As for SBH and interface polarizability, multiple factors might participate to interface magnetism weakening, like charge discontinuity driven intermixing, octahedral tilt induced in the first LSMO layers by octahedral connectivity at the interface, substrate strain and so on. A combinatorial heuristic approach to the definition of interface composition is a powerful tool to help understanding all these factors interplay and to enhance the interface magnetism, SBH or interface polarization.

## 4.2 The LSMO/STO interfaces

Incorporating a few uc of combinatorial $LSMO_x$ ($0 \le x \le 1$) at the LSMO/STO interface to modulate the chemical bonding, the carrier density and the polar discontinuity could potentially induce STO lattice polar distortion, SBH modulation, as well as restoring interface ferromagnetism.

### 4.2.1 Ferromagnetism at STO/LSMO interface

Before producing the described ICPLD heterostructures, we first optimized the LSMO physical properties, composition and thickness uniformity. The magnetic properties of the film were then characterized versus temperature using
a commercial Kerr magnetometer equipped with a cryostat (NanoMoke II, Durham Magneto Optics). As the magnetism at the LSMO/STO interface is weakened, the Curie temperature will depend on LSMO thickness for very thin films. To avoid this regime, we worked with 30nm thick LSMO films (~80uc). Several films were deposited on TiO₂ terminated (100) STO substrates with high-pressure RHEED monitoring (Staib/TSST) at various fluence, temperature and oxygen pressure. The optimized deposition conditions leading to a Tc = 341K were P₀₂=0.2 mbar, Tsub = 850°C, f = 5Hz and a fluence of 0.83 J/cm². RHEED oscillations were clearly visible during all the deposition process implying a layer by layer growth. X-ray diffraction patterns (Θ−2Θ) showed only (00l)ₚc reflections with thickness fringes attesting for the crystalline quality and the surface and interface smoothness. RSM confirmed epitaxial “cube on cube” growth of LSMO on STO. The homogeneity of the films in term of composition and magnetic properties over a 1cm² STO substrate was verified for thinner films, in the range where Tc is thickness dependent. A 20uc thick sample was deposited with vertical and horizontal scanning of the laser, staying in focus at the target position, and of the substrate respectively.

The Kerr magnetometer laser spot (diameter < 5µm) was scanned on the sample surface at fixed temperatures to measure magnetic hysteresis curves in 311 points spread across the sample surface. This (x,y) scan was repeated every 2.5K from room temperature to 350K after thermal stabilization. Each hysteresis curve was processed in order to extract saturation and remnant magnetization (Msat and Mr respectively). M(T) curves can then be reconstructed for each point on the sample surface, allowing to assess for Tc in each location. Maps reporting FM and paramagnetic (PM) areas of the sample are reproduced in **Figure 6a** (top) for various temperatures, the measurement points being indicated with black dots. The distribution of Tc is reported as an FM area percentage in **Figure 6a** bottom. Over 91% of the surface transit from FM to PM states on a temperature range less than 5K wide (325K < T < 330K) and 100% inside a 10K range. As LSMO's Tc is very thickness and composition sensitive, the tight Tc distribution indicates a good composition and thickness uniformity. We confirmed this uniformity with a WDS characterization over a 9x9 mm² area (25 x 25 = 625 points) of the same sample for La, Sr and Mn (JEOL 8530F). The small film thickness conjugated to the presence of Sr in the substrate did not allow to compute the composition with cationic ratios of the film. However, the WDS sensitivity is high enough to provide maps of relative variations for each element (see **Figure 6b**). The Sr map, with signal originating mostly from STO substrate, illustrates the electron beam stability (σSr~1%) which is crucial for point to point comparison. Note that the drop in the corner, corresponding to silver paste contact to evacuate the charges, was excluded from the statistical analysis. On La and Mn maps a similar slight slope is visible with a corresponding standard deviation of 4.6% and 4.9% respectively. WDS signal is strongly correlated to the thickness, therefore we can conclude that the thickness distribution is relatively tight.

The interface contribution to the overall magnetic signal increases as the LSMO thickness decreases. It is however difficult to predict the optimal LSMO thickness as the overall also decreases with thickness. A powerful aspect of CPLD is the possibility to deposit wedge-shaped layers with continuous thickness variation using shadow-masking (see **Figure 7a**). Before inserting the ICPLD layer, we checked the thickness control on two LSMO wedges, spanning from 8 uc to 76 uc, by measuring Tc versus (x,y) and temperature. The obtained Tc are represented **Figure 7b** with standard deviation represented as error bars. Tc noticeably decreases below 30 uc with an acceleration below 20 uc. In the inset of **Figure 7b** is represented a Tc map of wedge#2, with the measured points represented as black dots. Constant nominal thickness levels are vertical, with thickness variation along x. A 10K color increment is used and one can see that the
<figure><img src="image_7.png" /><figcaption>**Figure 6.**
(a) Magnetic state maps of LSMO 20 uc film at various temperature (top) and Tc distribution at the sample surface (bottom). (b) WDS signal for La, Sr and Mn over 9x9 mm² of the same sample.</figcaption></figure>

<figure><img src="image_8.png" /><figcaption>**Figure 7.**
(a) Schematic of LSMO wedge #2. (b) Tc versus thickness and (inset) Tc map for wedge #2</figcaption></figure>

lines separating the adjacent areas are almost vertical, attesting the good control of thickness variation in the wedge.

To synthesize the ICPLD LSMOx layer we used LaMnO₃ (LMO) and SrMnO₃ (SMO) targets with the deposition parameters identified for LSMO, including laser and substrate stage scans. Deposition rate was evaluated using RHEED oscillations. A 3 uc thick LSMOx layer (0 ≤ x ≤ 1) was deposited onto TiO₂-terminated STO substrate, followed by a LSMO wedge with thickness variation direction perpendicular to LSMOx composition gradient. A schematic representation of this sample is represented **Figure 8a**.

M(H) cycles were acquired versus position (512 sites) and temperature (120 temperatures with 90K < T < 340K) automatically during a few days. Then M<sub>Sat</sub> and M<sub>r</sub> were extracted for each loop, M(T) curves reconstructed and T<sub>c</sub> estimated for each (composition x, thickness t) doublet. **Figure 8b** presents the T<sub>c</sub> curves plotted versus t<sub>LSMO</sub> for various Sr content x. One can see that the variation of T<sub>c</sub> versus t<sub>LSMO</sub> depends on x, and in particular for t<sub>LSMO</sub> > 7 uc, the less Sr the more rapid is the T<sub>c</sub> decrease. Going from LMO to SMO at t<sub>LSMO</sub> = 10 uc, T<sub>c</sub> is increased by 60K (blue arrows **Figure 8b**). Furthermore, to reach a given Curie temperature
<figure><img src="image_9.png" /><figcaption>**Figure 8.**
(a) Schematic representation of the ICPLD LSMOₓ / wedge LSMO sample. (b) Curie temperature curves versus tLSMO for various Sr content x in LSMOₓ layer</figcaption></figure>

of 240K, one needs 8 uc of LSMO on top of SMO and more than 11 uc of LSMO on
top of LMO (red arrows in **Figure 8b**). One can compare these results to the one
obtained for x=0.29 (pink curve in **Figure 8**) where the heterostructure is similar
to a simple LSMO/STO interface. Inserting a 3 uc SMO layer at the LSMO/STO
interface proves to be beneficial in terms of Tc for tLSMO >7 uc. However we observe
a cross-over for tLSMO ≤ 7 uc. The Tc decrease with tLSMO accelerates for Sr rich
compositions, and no magnetism could be detected at tLSMO = 5 uc for 0.29 ≤ x ≤ 1.
On the contrary, the lower the Sr content the higher the Tc for 0 ≤ x ≤ 0.21 at tLSMO
= 5 uc. This reinforcement of FM for LMO coincides with an important increase of
the coercive field to values higher than usually observed for LSMO (Hc > 300 Oe
at T = 100K). This is compatible with a second FM phase, harder than LSMO and
in contact with it. LMO is antiferromagnetic (AFM) in bulk form. However, several
studies reported FM LMO films on STO substrate down to 6 uc. [e.g. 42] In this
article, the transition from AFM to FM has been attributed to an electronic
reconstruction at the interface originating from the polar nature of the LMO. In
our case the LMO layer is topped by LSMO, and it is quite possible that by
proximity effect and/or stress LMO becomes FM at 3 uc thick.

4.2.2 Band alignment at LSMO/STO interface

We now turn to interface issues arising in tunable capacitors with thinned FE
film i.e. the increased influence of dead ferroelectric layer on tunability and the
increased leakage current. As discussed previously the insertion of a LSMOx
ICPLD layer at LSMO/BST interface may increase interface polarizability and
modulate SBH. In order to easily disentangle spontaneous and chemically induced
polarizations we choose to work with a non-polar composition of BST i.e. STO. We
deposited onto TiO₂-terminated STO substrate 38 uc of LSMO followed by 3 uc of
LSMOₓ (0 ≤ x ≤ 1) and in the direction perpendicular to the gradient a STO wedge
(3-15 uc) keeping an access to both LSMO and LSMOₓ with the deposition
parameters described above.

A schematic of the sample structure is represented **Figure 9a**. The sample was transferred into an ultra-high vacuum atomic force microscope chamber (UHV-AFM Omicron) without breaking the vacuum. The AFM image presented in **Figure 9b** was taken about the red dot in **Figure 9a** with a total thickness of 56 uc. Terraces separated by steps of about 4Å, i.e. one perovskite cell parameter, are clearly visible (see profile in **Figure 9c**) attesting of the layer by layer growth up to 56 uc. There exists however some 2Å height features on the terraces indicating the probable existence of two terminations at the surface (SrO and TiO₂).
The sample was then air-exposed and inserted into a UV photoelectron spectroscopy chamber (UPS ESCALAB 250Xi Thermo Fisher) to evaluate the work function $W_F$ as a function of position. Although the surface contamination due to air exposure prevented to extract absolute $W_F$ values, the relative variations of $W_F$ with Sr content and STO thickness could be determined assuming a “uniform” surface contamination. UPS spectra were taken at various x content for STO thickness ranging from 0 to 9 uc. A zoom around the emission threshold of the He II UPS spectra (He II energy 40.8V, bias 4V) is shown in **Figure 10a** for LSMOx ($t_{STO}$ = 0 uc). From the thresholds one can estimated the $W_F$ reported in **Figure 10b**. A clear continuous decrease of the work function is observed as the Sr content

<figure><img src="image_10.png" /><figcaption>**Figure 9.**
(a) Schematic of the STO (sub.)/LSMO/LSMOx/STO heterostructure. (b) UHV-AFM image. (c) Profile from the white arrow in b).</figcaption></figure>

<figure><img src="image_11.png" /><figcaption>**Figure 10.**
(a) HeII UPS threshold for various x of air-exposed LSMOx $U_{bias}$ = 4V. (b) Corresponding extracted $W_F$ versus Sr content x</figcaption></figure>

<figure><img src="image_12.png" /><figcaption>**Figure 11.**
(a) HeII UPS threshold for various $t_{STO}$ at $x$=0.2 $U_{bias}$ = 4V. (b) Corresponding extracted $W_F$ versus $t_{STO}$ for various Sr content x</figcaption></figure>
increases. This trend is opposite to the downward Fermi level shift inferred from core-level XPS shift as a function of x reported in the literature [43] and seen by us (not shown). The counter-intuitive decrease of W<sub>F</sub> while E<sub>F</sub> decreases too is due to the LSMO<sub>x</sub> induced charge discontinuity variation at the surface. Going from LMO to SMO, the LSMO<sub>x</sub> terminal plane changes from Mn<sup>3+</sup>O<sub>2</sub><sup>2-</sup> to Mn<sup>4+</sup>O<sub>2</sub><sup>2-</sup>, i.e. with a surface charge per unit cell going from -1 to 0. The more negatively charged a surface is, the harder for an electron to escape from the surface, the higher the W<sub>F</sub>. [44]

The electrical nature of the contact between a metal and a semiconductor directly depends on the relative values of the metal W<sub>F</sub> and semiconductor electronic affinity E<sub>a</sub> for ionic semiconductor [45]. For E<sub>a</sub> > W<sub>F</sub> an ohmic contact forms, while for E<sub>a</sub> < W<sub>F</sub> a Schottky barrier is created. STO is generally considered an n-type ionic semiconductor with a Fermi level very close to the conduction band (i.e. E<sub>a</sub>~W<sub>F</sub>). As the LSMO<sub>x</sub> W<sub>F</sub> varies the LSMO<sub>x</sub>/STO contact nature might be affected. UPS spectra were acquired for various LSMO<sub>x</sub> Sr content every 200 µm along the STO wedge. A zoom of the corresponding UPS emission thresholds obtained for x=0.2 and 0 ≤ t<sub>STO</sub> ≤ 9 uc is shown in **Figure 11a**. The threshold position varies rapidly with t<sub>STO</sub> for thin STO layers then stabilizes. The W<sub>F</sub> estimated from the UPS thresholds for various (x, t<sub>STO</sub>) doublet are reported in **Figure 11b**. One can see the curves folding together toward a W<sub>F</sub> value of about 3.58 eV (relative) corresponding to the intrinsic STO work function.

Interestingly this value is inside the range of W<sub>F</sub> spanned within LSMO<sub>x</sub> (see pink part of **Figure 11b**). Looking at the evolution of W<sub>F</sub> vs t<sub>STO</sub> for thicknesses up to 3 uc, there is a clear transition from a downward to an upward bending as x increases. This reflects the band bending that occurs at the LSMO<sub>x</sub> / STO interface and implies that the contact is modified from Ohmic type to Schottky type.

This result is of importance regarding the optimization of the SBH in BST FE tunable capacitor in particular, but more generally for any metal/semiconductor contacts.

**5. Conclusion**

In this chapter we reviewed the qualities and limitations of PLD for the synthesis of oxides in general and for its use in combinatorial PLD synthesis (CPLD) in particular. We listed some counter-actions to mitigate the PLD limitations together with the mandatory steps to take before attempting reliable CPLD synthesis, i.e. demonstrating the control of both thickness and composition over the whole sample surface. We then detailed a statistical characterization approach to reliably interpret results from CPLD libraries of compounds. An example of this approach is presented, regarding the exploration of lead-free Ga-doped BiFeO<sub>3</sub> solid solution for MPB-related piezoelectric properties enhancement. Finally we described a new interface CPLD development (ICPLD) for the exploration of functional interface libraries. This combinatorial interface synthesis approach, with continuous lateral chemical modulation of a few atomic layers, is unique to the best of our knowledge. The effectiveness of ICPLD regarding the control of interface magnetism for magnetic tunnel junctions and energy band and Schottky barrier height tuning in ferroelectric tunable capacitors was demonstrated. This shows that ICPLD is a powerful tool to accelerate heterostructures functional properties enhancement.

**Acknowledgments**
Authors would like to thanks J-L Longuet from CEA Le Ripault (France) for the WDS characterizations presented here and Xavier Wallart from IEMN, UMR CNRS 8520, Villeneuve d'Ascq (France), together with Pascal Andreazza from ICMN, UMR CNRS 7374, for UPS measurements.

This work was funded by Région Centre within the projects INTIM-C 2015 and FLEXIGEN.

**Conflict of Interest**

No conflict of interest

**Author details**

Jérôme Wolfman*, Beatrice Negulescu, Antoine Ruyter, Ndioba Niang and Nazir Jaber†

GREMAN Laboratory, University of Tours, Tours, France

*Address all correspondence to: wolfman@univ-tours.fr

† Current address: LCM Laboratory, CEA Leti, Grenoble, France.

**IntechOpen**

© 2020 The Author(s). Licensee IntechOpen. This chapter is distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/3.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

<figure><img src="image_13.png" /></figure>
References

[1] Smith HM, Turner AF. Vacuum Deposited Thin Films Using a Ruby Laser. Applied Optics.1965;4:147-148

[2] Wu MK, Ashburn JR, Torng CJ, Hor PH, Meng RL, Gao L, Huang ZJ, Wang YQ, Chu CW. Superconductivity at 93 K in a New Mixed-Phase Y-Ba-Cu-O Compound System at Ambient Pressure. Physical Review Letters. 1987;58:908-910. DOI: 10.1103/PhysRevLett.58.908

[3] Chu CW. Cuprates—Superconductors with a Tc up to 164 K. In: Rogalla H, Kes PH, editors. 100 years of superconductivity. Boca Raton: CRC Press/Taylor & Francis Group. pp. 244–253. ISBN 9781439849484.

[4] Dijkkamp D, Venkatesan T, Wu XD, Shaheen SA, Jisrawi N, Min-Lee YH, McLean WL and Croft M. Preparation of Y-Ba-Cu oxide superconductor thin films using pulsed laser evaporation from high Tc bulk material. Appl. Phys. Lett. 1987;51:619. DOI: 10.1063/1.98366

[5] Chrisey DB, Hubler GK. Pulsed laser deposition of thin films. New York: Wiley; 1984. ISBN 0-471-59218-8

[6] Ohnishi T., Shibuya K., Yamamoto T. and Lippmaa M. Defects and transport in complex oxide thin films. J. Appl. Phys. 2008;103:103703. DOI: 10.1063/1.2921972

[7] Guozhen L, Wolfman J, Autret-Lambert C, Sakai J, Roger S, Gervais M, Gervais F. Microstructural and dielectric properties of Ba<sub>0.6</sub>Sr<sub>0.4</sub>Ti<sub>1-x</sub>Zr<sub>x</sub>O<sub>3</sub> based combinatorial thin film capacitors library. Journ. Appl. Phys. 2010;108:114108

[8] Daumont CJM, Simon Q et al. Tunability, dielectric, and piezoelectric properties of Ba<sub>1-x</sub>Ca<sub>x</sub>Ti<sub>1-y</sub>Zr<sub>y</sub>O<sub>3</sub> ferroelectric thin films. Journ. Appl. Phys. 2016;119:094107. DOI: 10.1063/1.4942924

[9] Trolier-McKinstry S, Muralt P, J. Thin Film Piezoelectrics for MEMS. J. Electroceram. 2004;12:7-17

[10] Jaffe B, Cook WR, Jaffe H. Piezoelectric Ceramics. London: Academic Press, 1971. p136. ISBN-0123795508

[11] Davis M. Picturing the elephant: Giant piezoelectric activity and the monoclinic phases of relaxor-ferroelectric single crystals. J. Electroceram. 2007;19:25

[12] Ahart M, Somayazulu M, Cohen RE, Ganesh P, Dera P, Mao H-K, Hemley RJ, Ren Y, Liermann P, Wu Z. Origin of morphotropic phase boundaries in ferroelectrics. Nature 2008;451:545-548 DOI: 10.1038/nature06459

[13] Noheda B, Cox DE, Shirane G, Gonzalo JA, Cross LE, Park S-E. A monoclinic ferroelectric phase in the Pb(Zr<sub>1-x</sub>Ti<sub>x</sub>)O<sub>3</sub> solid solution. Appl. Phys. Lett. 1999;74: 2059. DOI: 10.1063/1.123756

[14] Kaczmarek W and Pajak Z, Differential thermal analysis of phase transitions in (Bi<sub>1-x</sub>La<sub>x</sub>)FeO<sub>3</sub> solid solution. Solid State Commun. 1975;17: 807-810. DOI: 10.1016/0038-1098(75)90726-7

[15] Roginskaya YE, Tomoshpolskii YY, Venevstev YN, Petrov VM, Zhdanov GS. The Nature of the Dielectric and Magnetic Properties of BiFeO<sub>3</sub>. Sov. Phys.-JETP 1966;23:47

[16] Wei Y, Wang X, Zhu J, Wang X, Jia J. Dielectric, ferroelectric, and piezoelectric properties of BiFeO<sub>3</sub>–BaTiO<sub>3</sub> ceramics. J. Am. Ceram. Soc. 2013;96: 3163-3168. DOI: 10.1111/jace.12475

[17] Woodward DI, Reaney IM, Eitel RE, Randall CA. Crystal and domain structure of the BiFeO<sub>3</sub>–PbTiO<sub>3</sub> solid
solution. J. Appl. Phys. 2003;94:3313.
DOI: 10.1063/1.1595726

[18] Jaber N, Wolfman J, Daumont C, Négulescu B, Ruyter A et al. Laser fluence and spot size effect on compositional and structural properties of BiFeO₃ thin films grown by Pulsed Laser Deposition. Thin Solid Films 2017;634:107–111. DOI: 10.1016/j.tsf.2017.05.003

[19] Jaber N, Wolfman J, Daumont C, Négulescu B et al. Enhancement of piezoelectric response in Ga doped BiFeO₃ epitaxial thin films. J. Appl. Phys. 2015;117:244107; doi: 10.1063/1.4923217

[20] Abrahams SC, Reddy JM, Bernstein JL. Crystal structure of piezoelectric ferromagnetic gallium iron oxide. J. Chem. Phys. 1965;42:3957. DOI: 10.1063/1.1695868

[21] Belik A, Yi W. High-pressure synthesis, crystal chemistry and physics of perovskites with small cations at the A site. J. Phys.: Condens. Matter 2014;26:163201. DOI: 10.1088/0953-8984/26/16/163201

[22] Norton DP, Chakoumakos BC, Budai J, Lowndes DH, Sales BC, Thompson JR, Christen D. Superconductivity in SrCuO₂-BaCuO₂ superlattices: formation of artificially layered superconducting materials. Science 1994;265:2074-2077. DOI: 10.1126/science.265.5181.2074

[23] Ohtomo A, Hwang HY, A high-mobility electron gas at the LaAlO₃/SrTiO₃ heterointerface. Nature 2004;427:423–426. DOI: 10.1038/nature02308

[24] Reyren R, Thiel S et al. Superconducting interfaces between insulating oxides. Science 2007;317:1196-1199. DOI: 10.1126/science.1146006

[25] Yuan Y et al. Three-dimensional atomic scale electron density reconstruction of octahedral tilt epitaxy

in functional perovskites. Nature Com. 2018;9:5220. DOI: 10.1038/s41467-018-07665-1

[26] Huang W, Yin Y, Xiaoguang L. Atomic-scale mapping of interface reconstructions in multiferroic heterostructures. Appl. Phys. Rev. 2018;5:041110. DOI: 10.1063/1.5053426

[27] Tiggelman M et al. On the trade-off between quality factor and tuning ratio in tunable high-frequency capacitors. IEEE trans. Elec. Dev. 2009;56:2128-2136. Doi: 10.1109/TED.2009.2026391.

[28] York R, Tunable Dielectric Materials in: M. Steer and W.D. Palmer editors, Multifunctional Adaptive Microwave Circuits and Systems, New York: Scitech: 2009, chap. 4

[29] Lichtensteiger C et al. Ferroelectricity in ultrathin capacitors in: Pacchioni G, Valeri S editors, Oxide Ultrathin Films: Science and Technology, Weinheim: WYLEY;2012 Chap.12 ISBN 978-3-527-33016-4

[30] Nuñez M and Nardelli MB. First-principles theory of metal-alkaline earth oxide interfaces. Phys. Rev. B 2006;73:235422. DOI: 10.1103/PhysRevB.73.235422

[31] Wang Z et al. Anomalously deep polarization in (001) interfaced with an epitaxial ultrathin manganite film. Phys. Rev. B 2016;94:155307. DOI: 10.1103/PhysRevB.94.155307

[32] Yu P, Chu Y-H, Ramesh R. Oxide interfaces: pathways to novel phenomena. Mat. Today 2012;15:320-327. ISSN:1369 7021

[33] Koohfar S et al. Structural distortions at polar manganite interfaces. Phys. Rev. B 2017;96:024108. DOI: 10.1103/PhysRevB.96.024108

[34] Hikita Y et al. Termination control of the interface dipole in La₀.₇Sr₀.₃MnO₃/Nb:SrTiO₃ (001) Schottky junctions. Phys. Rev. B
2009;79:073101 (2009) DOI:
10.1103/PhysRevB.79.073101

[35] Werner R, Petrov A Yu, Alvarez Miño L, Kleiner R, Koelle D, Davidson B A. Improved tunneling magnetoresistance at low temperature in manganite junctions grown by molecular beam epitaxy. Appl. Phys. Lett. 2011;98:162505. DOI: 10.1063/1.3581885

[36] Sun J Z, Abraham D W, Rao R A, Eom C B. Thickness-dependent magnetotransport in ultrathin manganite films. Appl. Phys. Lett. 1999;74:3017. DOI: 10.1063/1.124050

[37] Bibes M, Balcells L, Valencia S, Fontcuberta J, Wojcik M, Jedryka E, Nadolski S. Nanoscale Multiphase Separation at La₂/₃Ca₁/₃MnO₃/SrTiO₃ Interfaces. Phys. Rev. Lett. 2001;87:067210. DOI: 10.1103/PhysRevLett.87.067210

[38] Tebano A, et al. Evidence of orbital reconstruction at interfaces in ultrathin La₀.₆₇Sr₀.₃₃MnO₃ films. Phys. Rev. Lett. 100, 137401 (2008).

[39] Yamada H, Ogawa Y, Ishii Y, Sato H, Kawasaki M, Akoh H, Tokura Y. Engineered interface of magnetic oxides. Science 2004;305:64-648. DOI: 10.1126/science.1098867

[40] Kavich JJ et al. Nanoscale suppression of magnetization at atomically assembled manganite interfaces: XMCD and XRMS measurements. Phys. Rev. B 2007;76:014410. DOI: 10.1103/PhysRevB.76.014410

[41] Boschker H et al. Preventing the reconstruction of the polar discontinuity at oxide heterointerfaces. Adv. Funct. Mater. 2012;22:2235-2240. DOI: 10.1002/adfm.201102763

[42] Wang R et al. Imaging and control of ferromagnetism in LaMnO₃/SrTiO₃ heterostructures. Science 2015;349:716-719. DOI: 10.1126/science.aaa5198

[43] Saitoh T, Bocquet AE et al. Electronic structure of La₁-xSrxMnO₃ studied by photoemission and x-ray-absorption spectroscopy. Phys. Rev. B 1995;51:13942. DOI: 10.1103/PhysRevB.51.13942

[44] Ryan Jacobs, John Booske, and Dane Morgan. Understanding and Controlling the Work Function of Perovskite Oxides Using Density Functional Theory. Adv. Funct. Mater. 2016, 26, 5471-5482

[45] R. Tung, Appl. Phys. Rev. 1, 011304 (2014)
https://doi.org/10.1063/1.4858400