Binod Paudelad,
Igor Vasilieva,
Mahmoud Hammourib,
Dmitry Karpovc,
Aiping Chend,
Valeria Lautere and
Edwin Fohtung*af
aDepartment of Physics, New Mexico State University, Las Cruces, NM 88001, USA. E-mail: efohtung@nmsu.edu
bDepartment of Physics, California State University, Los Angeles, CA 90032, USA
cSwiss Light Source, Paul Scherrer Institute, Villigen, 5232, Switzerland
dCenter for Integrated Nanotechnologies (CINT), Los Alamos National Laboratory, Los Alamos, NM 87544, USA
eNeutron Scattering Division, Oak Ridge National Laboratory, Oak Ridge, TN 37831, USA
fLos Alamos National Laboratory, Los Alamos, NM 87544, USA
First published on 29th April 2019
We utilize polarized neutron reflectometry (PNR) in consort with ab initio based density functional theory (DFT) calculations to study magnetoelectric coupling at the interface of a ferroelectric PbZr0.2Ti0.8O3 (PZT) and magnetic La0.67Sr0.33MnO3 (LSMO) heterostructure grown on a Nb-doped SrTiO3 (001) substrate. Functional device working conditions are mimicked by gating the heterostructure with a Pt top electrode to apply an external electric field, which alters the magnitude and switches the direction of the ferroelectric (FE) polarization, across the PZT layer. PNR results show that the gated PZT/LSMO exhibits interfacial magnetic phase modulation attributed to ferromagnetic (FM) to A-antiferromagnetic (A-AF) phase transitions resulting from hole accumulation. When the net FE polarization points towards the interface (positive), the interface doesn't undergo a magnetic phase transition and retains its global FM ordered state. In addition to changes in the interfacial magnetic ordering, the global magnetization of LSMO increases while switching the polarization from positive to negative and decreases vice versa. DFT calculations indicate that this enhanced magnetization also correlates with an out of plane tensile strain, whereas the suppressed magnetization for positive polarization is attributed to out of plane compressive strain. These calculations also show the coexistence of FM and A-AF phases at zero out of plane strain. Charge modulations throughout the LSMO layer appear to be unaffected by strain, suggesting that these charge mediated effects do not significantly change the global magnetization. Our PNR results and DFT calculations are in consort to verify that the interfacial magnetic modulations are due to co-action of strain and charge mediated effects with the strain and charge effects dominant at different length scale.
In this article, we report ME coupling at the interface of La0.67Sr0.33MnO3 (LSMO) and PbZr0.2Ti0.8O3 (PZT) layers grown on a Nb-doped STO (001) substrate and investigate interfacial magnetism modulations by switching the ferroelectric polarization orientations in the PZT layer. Hole doped lanthanum manganites such as La1−xAx MnO3 (A = Sr, Ca, Ba) well recognized for their colossal magnetoresistance property exhibit varied phases as a function of hole dopant concentration and temperature, as they carry strong interplay between charge, lattice, orbital and spin degrees of freedom.11–15 La1−xSrx MnO3 with x = 0.33 is the optimum doped manganite owing to the half-metallic behavior and optimized spin polarization with relatively high Curie temperature (TC = 360 K); which is manifested through dominant double-exchange mechanisms that delocalize eg electrons.16–18 PbZrxTi1−xO3 the solid solution of PbZrO3 and PbTiO3 with x = 0.2 was chosen in our experiments because of its high Curie temperature and large saturation polarization.19,20 For this reason, the LSMO/PZT heterostructure is considered as one of the best building blocks to fabricate ultra-low power spintronic devices that can be controlled with electric fields.
The methods of inducing and controlling magnetism across the interface of ferroelectric and magnetic layers in heterostructures have been profoundly studied in recent years since they are promising to replace energy consuming current controlled spin transfer torque (STT) memory devices.21–23 Recent studies have shown that the interface magnetoelectric coupling persists from a few unit cells to tenths of nanometers. However, understanding the mechanism of magnetoelectric coupling across such length scales under device working conditions requires a combination of theoretical and non-invasive experimental probes. These probes can access both the local distribution and changes in the magnetic polarization with atomic resolution across tens of unit cells. The electronic reconstructions at atomic scales via electric field induced strain coupled magnetism tuning in manganites, due to changes in electron hopping rates and associated changes in the atomic bond angles and lengths, have been widely studied.5,24–29 Electric field induced charge (electron–hole pairs) mediated coupling across FM/FE interfaces also alters the charge carrier concentrations at the FM oxide layers via an accumulation or depletion process depending on the direction of FE polarization.30 The interfacial magnetism and thus ME coupling can be manipulated through changing the magnetic moments due to the flipping of spin ordering or orbital reconstruction, changing the exchange interactions (resulting in competition between different magnetic phases), and by changing global magnetization (a result of changes in magnetic anisotropy).23 Charge mediated magnetism and ME coupling in LSMO/PZT multiferroic heterostructures has been extensively reported31–34 and some studies have ME coupling associated with the strain-mediated effects across the FE/FM interface. However, to understand the ME coupling that is persistent from Angstroms to tens of nm length-scales, we need probes that are capable of disentangling coupled charge and strain effects with local and atomic scale sensitivity. In this article, we combine PNR and DFT to elucidate the mechanisms of charge vs. strain mediated ME coupling across the PZT/LSMO interface and we explore contributions from an external electric field driven strain and charge mediated ME coupling.
Extensive experimental and theoretical work has been performed to understand magnetic modulations at LSMO/PZT interfaces.18,31,33–40 Vaz et al.31 studied ME coupling effects across the LSMO/PZT interface and attributed the observed ME coupling to interface charge mediations. Using X-ray absorption spectroscopy (XAS), they demonstrated that the modulation of Mn-ion valency is responsible for polarization induced magnetism changes at the LSMO interface region. They also argued that the observed suppression of magnetization at the interfacial layer associated with charge accumulation was attributed to strongly depopulated 3d z2 orbitals. This, in turn, favours the idea that super exchange interactions lead to the A-AF coupling at this interface.23 The authors concluded that the origin of the observed ME coupling was completely charge mediated i.e. not strain – driven. Spurgeon et al.34 studied the polarization surface charge screening, which induced magnetic phase gradients at LSMO/PZT interfaces. The authors utilized electron energy loss spectroscopy (EELS), PNR and DFT calculations to study the nature of magnetic phase transition across the interfaces.34,35 They also observed an asymmetrical magnetic phase gradient accompanied by a 30% difference in local FE polarization at the top and bottom interfaces. Interfacial spin state control of LSMO/PZT interfaces as a function of polarization orientations was studied by X. Ma et al.36 They observed that the spin states of a one-unit cell containing Mn ions changed from FM to A-AF due to an exchange interaction without a net change in the bulk magnetization. Zhou et al.37 studied the thickness dependence of LSMO/PZT heterostructures on interfacial ME coupling by measuring valence states of Mn atoms from XAS and observed an increase in Mn valency with increasing thickness of LSMO/PZT. Recently, Mayer et al. investigated the evolution of magnetic moments at the LSMO/PZT interfaces using PNR.38 The authors noticed that the surface and interface magnetization of LSMO was enhanced with the introduction of an FE capping layer. They observed enhanced interfacial magnetization for the polarization direction oriented towards interface and suppressed for the opposite direction. The thickness of the layer at the interface where the magnetization changed was found to be 2 nm, which agrees with the field screening.
Burton and Tsymbal39 studied the electric field induced magnetic reconstruction on La1−xAxMnO3 (LAMO)/BTO (A = Sr, Ca or Ba) interfaces for LAMO residing at FM and A-AF phase boundaries. This comprehensive study confirmed the occurrence of an A-AF phase at the interfaces for every phase boundary of all hole doped manganites. Chen and Ismail-Beigi18 studied ME coupling across LSMO/PTO interfaces using DFT. They found that the interfacial spin reconstruction of Mn atoms, resulted from carrier modulation due to polarization, are accounted for the changes in magnetization. Recently our group (M. Hammouri et al.)40 utilized DFT to study ME coupling across the LSMO/PZT interface. We observed the modulation of interfacial magnetization due to the presence of antiferromagnetic (A-AF) and FM phases at the LSMO/PZT interface. We also noticed that the balance of the magnetic (A-AF and FM) modulation could be altered by an electric field induced FE polarization.
However, all the charge mediated results reported above are not convergent and consistent with the bulk phase diagram of La1−xSrxMnO3.15 Vaz et al. and Meyer et al. found enhanced magnetization at LSMO (x = 0.2)/PZT interface even though hole accumulation brings LSMO towards the A-AF insulating phase and suppressed magnetization for hole accumulation that brings LSMO towards the FM metallic phase.31,38 Moreover, magnetic suppression at the LSMO/PZT interface is frequently encountered due to polar discontinuity and interface mixing.41,42 Spurgeon et al. found local polarization asymmetry at two LSMO interfaces with PZT and reported hole depletion for both interfaces instead of accumulation on bottom and depletion on top, which is not in agreement with the phase diagram of LSMO.34 Along with charge screening effects, many effects and phenomena that can contribute to magnetic modulations at the LSMO/PZT interfaces co-exist simultaneously. Issues that have not been addressed include (i) whether strain and charge mediation co-exist, (ii) and if so at what length scales do charge screening and the strain effect persists to contribute to the reported ME coupling?
To address these issues, we utilized PNR in consort with DFT calculations to understand the role of strain and charge mediated interfacial ME coupling in a heterostructure of LSMO/PZT. We utilized the Magnetism Reflectometer at the Spallation Neutron Source, Oak Ridge National Laboratory to provide polarized neutron beams to the LSMP/PZT sample.43 PNR is a probe that is particularly well suited for determining the nanostructures of magnetic thin films and multilayers.44,45
Reciprocal space maps, obtained with the X-ray diffraction (XRD), of the system, show the films are epitaxial and of good quality. Fig. 1(a) shows the peaks corresponding to PZT, STO and LSMO layers in the absence of an external electric field, whereas Fig. 1(b) displays the reciprocal space map of LSMO/PZT heterostructure on STO for the applied field of +2 V, showing the layers are still epitaxial with change in shape of PZT peak, which shows the polarization in PZT layer. 2Θ–ω (001) peaks for PZT for different electric fields as shown in Fig. 1(c). Negative and remnant fields have increased the out of plane tensile strain while positive fields decreased it. Applied electric fields on PZT/LSMO heterostructures cause lattice distortion along [001] and [100] directions, which is verified by shifting of all peaks towards smaller Qz values and larger Qx values respectively.
The crystal structure of bulk La0.67Sr0.33MnO3 was modeled using a 1 × 1 × 6 tetragonal supercell containing 30 atoms. The Brillouin zone was sampled with an 8 × 8 × 3 k-point mesh. To improve the accuracy of the theoretical description of the magnetic properties of LSMO, we used the Hubbard-corrected GGA+U exchange–correlation functional63–65 for Mn atoms. The Hubbard parameters for the Coulomb and exchange interactions were selected to be U = 2 eV and J = 0 eV, respectively. These parameters have been used in previous DFT studies of the electronic and magnetic properties of LSMO.40 The crystal structure of bulk PbZr0.2Ti0.8O3 was modeled using a 1 × 1 × 5 tetragonal supercell containing 25 atoms. The Brillouin zone was sampled with a 5 × 5 × 4 k-point mesh. The La0.66Sr0.33MnO3/PbZr0.2Ti0.8O3 (001) interface was modeled using a tetragonal supercell containing 6 unit cells of LSMO and 5 unit cells of PZT. The supercell was constructed using the LSMO layer containing Sr atoms in the 2-nd and 4-th unit cells and the PZT layer containing a Zr atom in the middle cell. The termination of the interfacial layers of LSMO and PZT in our theoretical model was consistent with the experimental studies, in which PZT was connected to the predominantly La-terminated surface of LSMO. We selected left-right symmetric positions for the substitutional Sr atoms in the modeled layer of LSMO to allow a more meaningful comparison of the distribution of Mn magnetic moments at the left and right sides of the LSMO layer. The Brillouin zone was sampled with a 6 × 6 × 3 k-point mesh. The total thickness of the modeled LSMO/PZT interfacial region was approximately 40 Å.
The optimized structures of tetragonal supercells for bulk LSMO, bulk PZT, and the LSMO/PZT (001) interface were calculated by DFT.40 A positive and negative uniaxial strain was applied to the modeled structure by expanding and contracting the size of the supercell in the direction orthogonal to the LSMO/PZT interface by ±1.5% and ±3%. After that, the structure of the LSMO/PZT interface was optimized by simultaneously relaxing the positions of all atoms and the lattice parameters in the plane parallel to the interface, while keeping the lattice parameter in the direction orthogonal to the interface constant.
To understand the ferroelectric polarization induced magnetic modulations at the PZT/LSMO interface, PNR magnetic depth profiles of LSMO for different polarization states are shown in Fig. 2(d). Magnetic depth profiles in the figure are not uniform and a few unit cells of magnetic dead layers are distinctly visible at left interfaces, which is attributed to the polarization of the ferroelectric PZT layer. Such changes in magnetization profile at the interface due to electric polarization are caused by magnetoelectric coupling. Polarization surface charges are produced at both surfaces of PZT, with polarization pointing either towards or away from the interfaces. These PZT surface charges screen the charges of equal and opposite sign at the LSMO interface to modulate the carrier concentrations either by accumulation or depletion of holes. Carrier modulation is restricted to only a few unit cells because of Thomas-Fermi screening lengths.66 Our LSMO is optimum doped half metal with a concentration of 0.33, which has larger screening depths as compared to metallic conductors. Therefore, the electric field induced interfacial magnetic modulations on LSMO are strongly noticed as shown in Fig. 2(d). For positive polarization with a field of +2 V (i.e. polarization pointing towards the interface), the magnetic modulation is almost zero because the depletion of holes cannot switch FM into A-AF phase as the screening charges at the LSMO surface would not be enough to achieve that. The hole depletion drives the magnetic phase towards the left of ferromagnetic metal (FMM) phase (x < 0.33) that may bring LSMO interface towards the ferromagnetic insulator (FMI) phase, thereby stabilizing the ferromagnetic phase at the expense of Curie temperature.67 If the polarization is negative, i.e. it points away from the interface due to the field of −2 V, accumulation of holes takes place at the LSMO interface due to positive screening charges. Few unit cells of LSMO at the interface are modulated as it undergoes a phase transition from FMM to AFM close to x = 0.5, thereby suppressing the FM phase. This phase transition favors the super-exchange interactions flipping the magnetic moments of Mn atoms close to the interface and hence magnetic reconstruction may arise within a few unit cells (2 nm) of LSMO. The PZT layer will be in a self-polarized state at zero field with polarization pointing away from the interface and will still accumulate holes at the LSMO interface, modulating magnetization due to the presence of almost same AF layer that was for −2 V. This supports the strong interplay between charge and spin degrees of freedom in manganites. Such interfacial modulation is accompanied by an increase in bulk magnetization for two accumulation states (i.e. for −2 V and 0 V) as compared to depletion (i.e. for +2 V). Within the FM-metallic phase, the magnetization of the LSMO layer becomes proportional to doping concentration as a higher hole concentration provides strong double exchange interactions to hop more carriers between Mn3+ and Mn4+. This is consistent with the increase in bulk magnetization of LSMO while going from +2 V to −2 V through 0 V as a result of tensile and compressive strains along with charge modulations due to applied fields. The strain coupling on magnetostrictive LSMO from the PZT layer modulates the magnetism and such modulations vary from interface to the surface.35 Therefore, strain effects are mostly contributing to bulk magnetic modulations whereas charge effects are only limited to interfaces modulations as shown in PNR. Fig. 2(e) shows the spin asymmetry ratio as a function of Qz. This gives the orientation of magnetic moments perpendicular to the interface.68
To investigate the influence of PZT polarization and applied uniaxial strain on the magnetic properties of the LSMO/PZT interface, we calculated the total energies and magnetic moments of the LSMO/PZT interfacial region in different magnetic configurations. According to the phase diagram of tetragonal manganites reported by Fang et al.,69 bulk LSMO has ferromagnetic (FM) order of Mn magnetic moments at the Sr doping level of x = 0.33 and the lattice ratio c/a = 1. An A-type antiferromagnetic (A-AF) order is favored when the c/a lattice ratio is reduced to approximately 0.97. Other types of magnetic order, such as C-type AF and G-type AF, are expected to be significantly higher in energy if the c/a lattice ratio is close to 1.0. In agreement with the phase diagram of tetragonal manganites, our calculations predicted that the FM configuration was lower in energy than the A-AF configuration for bulk La0.66Sr0.33MnO3. In contrast, we found that the LSMO/PZT interfacial region generally favored a mixed magnetic configuration with FM layers on the right side where the polarization of PZT pointed towards the surface of LSMO and A-AF layers on the left side where the polarization of PZT was directed away from the surface of LSMO. The dependence of the total energy of the LSMO/PZT interfacial region on the magnetic configuration of LSMO for different values of applied uniaxial strain is shown in Fig. 3(b).
In the absence of applied strain, the magnetic configuration with a single A-AF layer on the left side of the LSMO slab had the lowest total energy.40 A positive strain stabilizes the FM magnetic configuration of LSMO and decreases the number of A-AF layers on the left side of the LSMO slab, whereas a negative strain stabilizes the A-AF magnetic configuration of LSMO and increases the number of A-AF layers on the left side of the LSMO slab. This result is consistent with the increase in bulk magnetization while switching from +2 V to −2 V in experiment. The lowest-energy magnetic state of the LSMO slab was found to be completely FM for ε ≥ 3%, and completely A-AF for ε ≤ −3%.
Fig. 3(c) presents a schematic illustration of the influence of applied uniaxial strain on the distribution of magnetization across the LSMO/PZT interfacial region. The number of A-AF layers on the left side of the LSMO slab in the lowest-energy magnetic configuration increased with a negative applied strain (compression) and decreased with a positive applied strain (expansion). A simulated magnetic profile of the LSMO layer for different values of applied strain is displayed in Fig. 3(b). The simulated profile was obtained as a statistical distribution of the FM, A-AF, and mixed FM/A-AF LSMO/PZT magnetic configurations at T = 300 K and B = 0.1 T. Our study demonstrated the possibility of the existence of a layer of A-AF LSMO at the left LSMO/PZT interface, where the polarization of PZT was directed away from the surface of LSMO. The presence of the A-AF phase of LSMO in the LSMO/PZT interfacial region could be responsible for the presence of a nonmagnetic layer near the LSMO/PZT interface shown in Fig. 2(c). We found that a positive strain increased, and a negative strain decreased the magnetization of LSMO. The results of our calculations demonstrated that the spatial distribution of magnetization inside the LSMO slab is strongly influenced by the direction of the PZT polarization and by the applied strain.
When an electric field is applied to the ferroelectric layers, the modulation of carriers at the interfaces and distortion of lattices in response to the direction of the field take place simultaneously. This implies that charge and strain mediated effects are contributing collectively to modulate magnetic properties across LSMO as shown from PNR and DFT. The lattice distortion results as the strain is transferred to the magnetic layers from the FE layers, thereby changing the magnetic properties of the interfacial region due to magnetostrictive effects. To explore the interplay between electric field induced strains and bulk magnetic moments, we simulated magnetic phases on LSMO for different values of applied strain because the magnetic properties of LSMO strongly depend on the lattice parameters. The number of A-AF layers at the left interface shown in Fig. 3(b) and (c) corresponds to the magnitudes of bulk magnetization/magnetic moments of LSMO but does not represent the length of interfacial ME coupling as charge effects are less noticeable, as explained by the DFT results. Fig. 3(c) shows that out of plane compressive strain decreases magnetization of LSMO and out of plane tensile strain increases magnetization of LSMO regardless of the polarization direction of PZT i.e. for both accumulation and depletion. This is because the coupling of strain and magnetism is directly related to the distortion of MnO6 octahedra,12 where altering Mn–O bond angles triggers the Jahn–Teller distortion that directly affects the double-exchange mechanism. Even at the zero strain the co-existence of FM and A-AF phases is exhibited. This is attributed to the competition between two phases as LSMO shifts towards the A-AF phase for the negative polarization. The hopping of eg electrons is influenced by the strain states as the out of plane compressive strain reduces the hopping probability by increasing the bond lengths and decreasing the kinetic energy of the hopped electrons.14 Therefore, the A-AF phase starts to dominate the FM phase for larger compressive strains. For tensile strains, the electron hopping probability is very high, as there is a reduction in the Mn–O bond lengths and 3-D hopping starts to work, which favors double exchange interactions and increases the FM phase for larger tensile strains as illustrated in Fig. 3(c). The phase evolution of LSMO can also be explained in terms of orbital reconstructions. Since the eg electron orbital is doubly degenerated, such degeneracy can be lifted by either compression or elongation of Mn–O octahedra called Jahn–Teller distortions; which implies either favoring d3Z2−r2 or dx2−y2, i.e. either imposing tensile or compressive strains.18,70 Ferromagnetic ordering is isotropic along x-, y- and z-direction, but A-type antiferromagnetic ordering is intralayer ferromagnetic and out of plane interlayer antiferromagnetic.18 Under a compressive strain, the d3Z2−r2 occupancy decreases resulting in a lowering of the fraction of the FM phase. In contrast, a tensile strain increases the d3Z2−r2 occupancy that increases the fraction of the FM phase thereby increasing hopping along the c-direction.
Since the applied electric field not only distorts the lattices but also polarizes the PZT layer, these results can be explained with the strain states simulated in DFT and the charge modulation on LSMO due to charge screening from PZT polarization. According to the phase diagram of tetragonal manganites reported by Fang et al.,69 the magnetic phase modulations at interface and bulk regions are jointly influenced by strain states from tetragonal distortion of LSMO and hole concentration.64 The entire region above the phase boundary line is FM and below is A-AF phase. For DFT calculations, the magnetic phases are only studied residing on the nominal hole doping (x = 0.33), shown by vertical red line in Fig. 4. The obtained results from both studies are tentatively consistent with the phase diagram from Fang et al. that LSMO shows the mixed phase of FM and A-AF phases at c/a = 0.97 and x = 0.33. When accumulation tends towards x > 0.33, LSMO undergoes transition to the A-AF phase, but the depletion brings towards x < 0.33 that leads to FM-phase as shown in Fig. 4. Moreover, for c/a ≥ 0.97 with hole depletion, LSMO shifts towards the FM phase and for c/a ≤ 0.97 with hole accumulation LSMO shifts towards the A-AF phase as shown in the phase diagram.71 A larger c/a ratio results in the FM phase, whereas a smaller c/a ratio results in the A-AF phase. Interfacial A-AF layers appear when the PZT polarization points away from the interface and the FM phase appears when the PZT polarization points towards the interface. Our DFT calculations demonstrate that coupling between the PZT and LSMO layers significantly enhances the sensitivity of the magnetic configuration of LSMO to the applied strain. The high sensitivity of the magnetic configuration of LSMO to the applied strain can be explained by the coexistence of FM and AFM phases in the PZT/LSMO interfacial region. The simultaneous presence of FM and AFM phases near the PZT/LSMO interface can explain the experimentally observed variations of the magnetic profiles within the LSMO layer. Experimental results can be compared with the phase diagram residing at three LSMO doping concentrations corresponding to −2 V, 0 V and +2 V. The interfacial magnetization is suppressed with the formation of a magnetically dead layer up to ∼2 nm accompanied by the increase in bulk magnetization for 0 V and −2 V. However, the width of the interfacial magnetically dead layer for a negative polarization decreases while going from 0 V to −2 V, i.e. increasing the out of plane tensile strains, which is consistent with the DFT predictions shown in Fig. 3(c). The bulk magnetization has increased while going from 0 V to −2 V, as more Mn3+ ions find their Mn4+ counterparts to have double exchange coupling due to the increase in hole concentration from x = 0.33. The LSMO layer has suppressed magnetization up to larger distances for +2 V, but still remains in the FM phase. This result can be explained by the depletion of holes which decrease the total magnetization. The bulk and interfacial magnetic modulations for the three applied voltages are consistent with the results of our calculations, as the magnetization is increasing with increasing the values of the c/a ratio while going from +2 V to −2 V with 0 V in between. However, the length of A-AF layer has not increased in PNR while applying −2 V from 0 V, this might be a reason that the presence of Pt capping layer generally reduces the nominal hole concentration.71 Therefore, a tensile strain increases the magnetization and introduces the interfacial A-AF layer, whereas a compressive strain reduces the magnetization and retains its ferromagnetic phase intact at the interface.
Fig. 4 Phase diagram of LSMO as a function of Sr concentration (x) and tetragonality c/a.69 LSMO experiences a compressive strain when grown on STO. The phase transitions due to accumulation and depletion of holes in LSMO and varied tetragonal ratios of the sample are shown.69,71 |
This journal is © The Royal Society of Chemistry 2019 |