P. S.
Miedema
*a,
N.
Thielemann-Kühn†
bc,
I. Alonso
Calafell‡
b,
C.
Schüßler-Langeheine
b and
M.
Beye
a
aDeutsches Elektronen Synchrotron DESY, Notkestraße 85, 22607 Hamburg, Germany. E-mail: p.s.miedema@gmail.com
bInstitute Methods and Instrumentation for Synchrotron Radiation Research (FG-ISRR), Helmholtz-Zentrum Berlin für Materialien und Energie GmbH, Albert-Einstein-Straße 15, 12489 Berlin, Germany
cInstitut für Physik und Astronomie, Universität Potsdam, Karl-Liebknecht-Straße 24/25, 14476 Potsdam, Germany
First published on 20th September 2019
Electronic structure modifications due to strain are an effective method for tailoring nano-scale functional materials. Demonstrated on nickel oxide (NiO) thin films, Resonant Inelastic X-ray Scattering (RIXS) at the transition-metal M2,3-edge is shown to be a powerful tool for measuring the electronic structure modification due to strain in the near-surface region. Analyses from the M2,3-edge RIXS in comparison with dedicated crystal field multiplet calculations show distortions in 40 nm NiO grown on a magnesium oxide (MgO) substrate (NiO/MgO) similar to those caused by surface relaxation of bulk NiO. The films of 20 and 10 nm NiO/MgO show slightly larger differences from bulk NiO. Quantitatively, the NiO/MgO samples all are distorted from perfect octahedral (Oh) symmetry with a tetragonal parameter Ds of about −0.1 eV, very close to the Ds distortion from octahedral (Oh) symmetry parameter of −0.11 eV obtained for the surface-near region from a bulk NiO crystal. Comparing the spectra of a 20 nm film of NiO grown on a 20 nm magnetite (Fe3O4) film on a MgO substrate (NiO/Fe3O4/MgO) with the calculated multiplet analyses, the distortion parameter Ds appears to be closer to zero, showing that the surface-near region of this templated film is less distorted from Oh symmetry than the surface-near region in bulk NiO. Finally, the potential of M2,3-edge RIXS for other investigations of strain on electronic structure is discussed.
Lattice strain is typically measured through the analysis of lattice peak positions measured with hard X-ray diffraction. On nano-sized materials the low amount of material together with small interaction cross sections for the hard X-ray radiation only provides small signals, averaged over the whole film thickness. Especially for the analysis of strain in surface-near regions, signals are largely shielded by the stronger bulk diffraction. In addition, it is not the lattice changes due to strain themselves that are important in strain engineering, but it is the associated change in electronic structure, e.g., the band gap and valence and conduction band properties, that ultimately leads to enhanced properties of the material. Already for simple binary correlated oxides, e.g., cobalt(II) oxide (CoO)8,11 and nickel oxide (NiO),12,13 the strain of the system has been shown to have an influence on properties such as the orbital polarization or the spin and corresponding relative magnetization, or on catalytic activity. For more complex materials strain engineering provides an attractive alternative method for tailoring materials’ properties beyond what is possible through doping or chemical substitution.10,14
As combination of X-ray absorption and emission spectroscopy, Resonant Inelastic X-ray Scattering (RIXS), is specific for determining electronic structure properties such as the band gap,15 covalent interactions,16 potential energy surfaces17 and excited states18 with element- and oxidation-state specificity. In this report it will be shown that RIXS is also sensitive to changes in electronic structure due to induced strain. Usually for 3d-transition metal materials, Metal L2,3-edge RIXS (schematically 2p63dn → 2p53dn+1 → 2p63dn) is performed for bulk electronic structure studies, partially because it turns out to provide the best possible energy resolution due to sharpest core levels and thus highest chemical and functional information.19 However, to provide information on the electronic structure near surfaces, M2,3-edge RIXS (3p63dn → 3p53dn+1 → 3p63dn) is a viable alternative. In particular, for heterogeneous solid catalyst systems, the (near-) surface region provides the chemical activity and a detailed understanding of the changes of the electronic structure in this region is crucial. The X-ray attenuation length in (bulk) NiO is much shorter when addressing the M2,3 resonances instead of the L2,3 levels: across the full L2,3-edge resonances an average attenuation length of about 50 nm is observed for NiO, with detuning from the L3-edge maximum yielding probing depths of around 300 nm. For the M2,3-edge the probing depth is roughly 16 nm20 and thus provides much higher sensitivity to the near-surface region, while L2,3-edge RIXS probes a more significant part of the bulk of the sample. Therefore the M2,3-edge RIXS of 3d-transition metal materials is particularly well-suited to study the electronic structure changes in thin films induced by strain from the substrate.
In this research, we will provide theoretical Ni M-edge RIXS multiplet calculations with tetragonal distortions from octahedral (Oh) symmetry. We will show that as a function of distortion, certain signals in the RIXS spectra are changing and in particular that ratios between particular peak areas can be used to quantify the distortion from Oh symmetry. As such, the theoretically obtained ratios can be matched onto ratios obtained for experimental spectra of different NiO systems.
Nickel M2,3-edge RIXS measurements were performed with a predecessor of the SolidFlexRIXS-setup21 attached to the UE-112/PGM1 beamline22 at the BESSY-II synchrotron in Berlin, Germany. The Nickel M2,3-edge absorption spectra were recorded via sample drain current measurements.23 The first M3 resonance maximum occurred at 67.2 eV and was selected as excitation energy for the RIXS experiments, since it is sufficiently below the threshold of the 3p core ionization. The RIXS spectra for the nano-layer NiO-films required acquisition times of typically three to four hours with the used setup, while for the NiO bulk sample a measurement time of only 90 minutes was considered sufficient. The incident X-rays were horizontally polarized in order to reduce the elastic scattering contribution to the RIXS spectra and the samples were rotated with the sample normal around 60° to the incoming beam, resulting in the normal being 30° off the direction of the spectrometer mounted at a right angle to the incoming beam in the horizontal polarization plane.
The Quanty package24–26 was used to calculate Ni M-edge horizontally polarized X-ray absorption spectra (XAS) and RIXS calculations with horizontal incoming polarization (X in) and the orthogonal outgoing polarizations (Y + Z out), thereby simulating the experimental geometry with the RIXS spectrometer in the plane of the incoming polarization. Simulations were performed for Ni2+ in octahedral (Oh) symmetry with a crystal field value of 10Dq = +1.05 eV and with a negligible magnetic field of H = 0.000001 T. Distortions from Oh symmetry were taken into account with two parameters called Δ1 and Δ2, where Δ1 splits the original t2g levels and Δ2 the original eg levels of the Oh symmetry as shown in Fig. 1A.
![]() | ||
Fig. 1 (A) 3d-Orbital splitting from Oh (left) in a D4h-type symmetry with e (dxz, dyz), b2 (dxy), a1 (dz2) and b1 (dx2−y2). One-electron energies due to the (distorted) crystal field splitting are shown as function of Dq, Ds and Dt27,28 and the splitting due to Δ1 and Δ2 is shown as used in the quantity multiplet simulations performed here. (B) Calculated Ni M-edge RIXS with horizontally polarized excitation and emission from the orthogonally polarized components. (C) Calculated Ni M-edge horizontally polarized XAS. For panels B and C the parameters in the multiplet calculations with Ni2+ were 10Dq = 1.05, Δ1 = Δ2 = 0 eV and 0.3 eV Gaussian broadening and the excitation energy scale was shifted by 66.8 eV. The vertical line in panel C indicates the excitation energy of the RIXS slices for all calculations with varying Δ1, Δ2 (or Ds) used in further analyses. |
The energetic splitting of the d-orbitals can be equivalently formulated with Δ1 and Δ2 or with the symmetry parameters Ds and Dt. The one-electron energies as function of Dq, Dt and Ds are shown in Fig. 1A. The Δ1 and Δ2 approach can be used to test the behavior of the t2g and eg splitting independently. In addition, this approach can be used in tests of symmetry distortions that are different from the D4h-type. While ignoring Dt (Dt = 0 eV) to account only for small distortions to octahedral symmetry, for example due to Jahn–Teller distortions, the relationship between Ds and Δ1 and Δ2 is: Δ1 = 3 × Ds, Δ2 = 4 × Ds. As for such simulations with Ds distortions, the Δ2 = 4/3 × Δ1 relationship was fixed.
The simulations for Ni2+ with 10Dq = 1.05 eV and various distortions Ds in eV yield linearly polarized XAS and full RIXS maps. One set is shown in Fig. 1B and C for comparison with the experimental data. RIXS slices were taken at the vertical line indicated in the XAS spectrum of Fig. 1C. We chose a single value for 10Dq, for both ground and X-ray excited state for sake of simplicity, however it is known that 10Dq may be different in the ground state (and RIXS final state) and the X-ray excited state.29 Overall, the first band in these simulations agreed with the chosen 67.2 eV excitation energy of the experimental data. In addition, the overall calculated XAS shape agreed well with the experimental XAS of bulk NiO.30
In order to have an analysis of the areas of peak 2 and peak 3 in the simulated spectra that can directly be related to experimental data, were fitted six Voigt-profiles to the calculated spectra for a range of Ds-values in between −0.15 and +0.15 eV. An example of this analysis for a M2,3-edge RIXS slice of Ni2+ with 10Dq = 1.05 and Ds = 0.0033333 eV is shown in Fig. 2B. Fig. 2C summarizes the obtained peak area fitting results for different values of Ds, displaying the trends in spectral weight of peak 2 and 3 as a function of Ds: the peak 2 weight increases, while the peak 3 decreases with increasing Ds. These trends are becoming less clear for Ds ≥ 0.07 eV as will be discussed shortly. The peak 2/peak 3 area ratio as a function of Ds is shown in Fig. 2D, which may be directly compared to the values found for experimental data shown below. The peak area ratio in panel D is roughly linear with the distortion Ds, up to Ds ≈ 0.07 eV. For larger Ds values, the fitting procedure using six peaks begins to deviate from the simulated spectra, since peak 3 starts to split into two components starting around this value of Ds = 0.07 eV (for example, the green solid line with Ds = 0.1 eV in Fig. 2A shows a side-band component in peak 3). For the remaining discussion we restrict ourselves to Ds values where this splitting is small enough to be ignored.
We further mention that the energetic positions of peak 2 and 3, which for both peaks increase slightly with Ds (see e.g., Fig. 2A), yields additional qualitative information about the value and sign of the distortion parameter Ds, especially when symmetry distortions due to strain become large, e.g., when the value of Ds is on the order of Dq (= 10Dq/10).
The experimental Ni M2,3-edge XAS of all the investigated samples appear similar to the NiO bulk spectrum measured by Chiuzbăian et al.30 The 67.2 eV feature was chosen as the excitation energy for our experimental M2,3-edge RIXS. The measured RIXS spectra at an excitation energy of 67.2 eV for the different film thicknesses (10, 20 and 40 nm) layers of NiO on MgO (NiO/MgO), as well as a 20 nm film of NiO on a 20 nm film of Fe3O4 on MgO (NiO/Fe3O4/MgO), and the spectrum from a NiO bulk crystal are shown in Fig. 3A as function of energy loss in eV. Compared to the RIXS spectrum of the NiO bulk (blue line) in Fig. 3A, the inelastic features at about 1 eV (peak 2) and 1.6–1.7 eV (peak 3) are enhanced for all the nano-layered NiO samples. The elastic line (peak 1) is about 100 to 1000 times higher than the energy loss features and has been cut off for clarity. Due to potential saturation effects in the elastic peak, its shape gets distorted. Therefore the peak analysis was performed as follows: the elastic line (peak 1) is assumed to be perfectly symmetric and the background found at negative energy loss is subtracted from the positive energy loss side before integrating the areas of peaks 2 or peak 3. The peak analysis results for the different samples are displayed in Fig. 3B and C. Both peak 2 and peak 3 areas in Fig. 3B are decreasing as function of increasing nano-layer film size on MgO and are approaching the peak areas obtained for NiO bulk.
The peak ratios of the experimental data in Fig. 3C for the 10, 20, 40 nm NiO/MgO and bulk samples are around 0.5–0.55. Comparing these ratios of the experimental data with the analyses of the Quanty calculations in Fig. 2D suggests negative values of Ds of about −0.10 ± 0.02 eV for all those samples, e.g., the ratios of 0.54, 0.52, 0.52 and 0.50 relate to rounded-off Ds values of −0.10, −0.11, −0.11 and −0.11 eV for the respective 10 nm, 20 nm, 40 nm NiO on MgO and bulk NiO samples. Observing that the sign of Ds is negative is logical, since the lattice constants of the MgO and Fe3O4 substrates are larger than the NiO lattice constant. That allows the Ni–O bonds in-plane to extend and in a crystal field description that leads to energy lowering of for example the 3dx2−y2 orbital (b1 in Fig. 1A) compared to the 3dz2 (a1) orbital. In order to simulate such an effect of a lower relative energy for the b1 orbital compared to the a1 orbital, Ds should be negative.
Since the signal in the M2,3-edge RIXS experiment is dominated by the topmost 15 nm, the peak area ratios show that the growth of NiO layers on MgO beyond a 20 nm thickness already becomes more bulk-like, e.g., the probed surface-near layers further away from the MgO substrate (e.g., the surface-near region of both the 20 and 40 nm NiO on MgO) are having a similar lattice as NiO bulk. Strain effects due to the substrate are expected to decrease with increasing layer thickness, ultimately leading to surfaces that resemble the electronic structure similar of the bulk material. This trend is indeed been found in our M-edge RIXS electronic structure analyses for NiO/MgO, e.g., the average peak area ratio of 10 nm NiO/MgO is different from the NiO bulk, while the average peak area ratio of 40 nm NiO/MgO is already very similar to the NiO bulk, showing the robustness of the M-edge RIXS method for electronic structure characterization due to strain. Taking into account the error bars all nanolayered NiO films on MgO can be considered similar to bulk NiO. (Note that the error bars in Fig. 3C are rather large due to the analysis procedure applied to the experimental data and derives from the noise level of the acquired data.)
Focusing on the experimental analysis obtained for the bulk NiO, the symmetry of bulk NiO is often considered to be Ni2+ in an Oh surrounding. However, it is known since the 1940s that crystalline bulk NiO has a small orthorhombic distortion.31 In addition, in the case of M-edge RIXS, the (near-)surface region of bulk NiO is probed which may be different from the overall bulk NiO as pointed out by Freitag et al.32 This validates our observation of a substantial distortion from Oh symmetry for bulk NiO.
In order to test the Ni M-edge RIXS procedure to a different induced strain, a 20 nm NiO/Fe3O4/MgO was grown and the experimental RIXS spectrum is shown as well in Fig. 3A (light blue line). From the experimental RIXS spectrum of 20 nm NiO/Fe3O4/MgO in Fig. 3A, it is already clear that the energy loss peaks 2 and 3 are more intense than for the 20 nm NiO/MgO (Fig. 3B): one easily observes that the peak areas for 20 nm NiO/Fe3O4/MgO (marked as 20 nm/Fe3O4 on the NiO thickness scale in Fig. 3B) are approximately 1.5–2 times higher than the 20 nm NiO/MgO. The peak 2/peak 3 area ratios in Fig. 3C for 20 nm NiO/MgO and NiO bulk samples compared to the 20 nm NiO/Fe3O4/MgO substantiate that the 20 nm NiO/Fe3O4/MgO has a different induced strain than the NiO/MgO samples. The 20 nm NiO/Fe3O4/MgO has a peak ratio of ∼0.75 relating to a Ds distortion in the quanty calculations of −0.025 eV, e.g., a value which is much closer to zero than the previously discussed NiO/MgO and bulk NiO samples. As mentioned before, the lattice parameter mismatch of Fe3O4 with NiO (0.48%) should be smaller than the mismatch of MgO with NiO (0.72%). Thus, surprisingly, we observed that the obtained Ds distortion of 20 nm NiO/Fe3O4/MgO differs more from the bulk NiO than the 20 nm NiO/MgO does. We should point out here that the Fe3O4 unit cell has to match two NiO unit cells so that could somewhat explain our unexpected observation. In addition, the fact that Fe3O4 grown on MgO also has a lattice mismatch, which might not be compensated yet after 20 nm of Fe3O4 layers, needs also to be taken into account. Ultimately this seems to lead to NiO nano-layer growth on top that is more Oh-symmetric than the near-surface crystal structure of a bulk NiO crystal.
For the theoretical analyses shown above, 10Dq = 1.05 eV was fixed agreeing well with the experimental data (peak 2 energy loss position) and published values for NiO.30 In order to analyze the robustness of our Ds calibration with respect to the choice of 10Dq, the dependence of the above analyses on different values of 10Dq is checked by varying the distortion Ds for fixed 10Dq = 0.8 and 10Dq = 1.2 eV. For 10Dq = 0.8 eV similar trends as above remain present, although absolute values for the peak ratios are somewhat different. Since 10Dq = 0.8 eV shows the energetic position of peak 2 at 0.8 eV and peak 3 at about 1.6 eV, this 10Dq is considered not to be realistic for matching the experimental data. For a 10Dq value of 1.2 eV, peak 3 generally shows a double structure, which is in disagreement with the experimental data shown in Fig. 3A, so the energetic position and shape of peaks 2 and 3 strongly depend on the choice of 10Dq and agreement with the experiment enforces values 1.05 ± 0.05 eV.
The provided multiplet analyses demonstrate the potential of obtaining lattice distortion related electronic structure information with M2,3-edge RIXS. Furthermore we performed M-edge RIXS calculations for other 3d-transition metal ions (e.g., Co2+, Fe2+) and observe clear relationships between distortions of the ligand field around the metal ion and features in the RIXS spectra. Thus, a very similar M-edge analysis procedure for electronic structure modifications could also be used for analysis of, for example, cobalt oxide nano-films.8,11
Although the SolidFlexRIXS setup that we used was not set to the optimal resolution in order to improve count rates, we want to stress here that for M-edge RIXS nowadays sub-100 meV resolutions are achieved with rather compact instruments, very much in contrast to L2,3-edge RIXS of 3d-transition metal materials, where such energy resolution can only be achieved with huge experimental effort. Also it should be noted that the energy range of the M2,3-edges of 3d-transition metal materials can be achieved with lab-high-harmonic-generation (HHG) sources. Although current HHG sources do not provide sufficient flux to perform the here presented RIXS experiments, we expect future developments increasing photon flux in combination with increased efficiency of RIXS optics,33,34 such that one could foresee that this type of strain engineering analyses with M-edge RIXS may be performed in the lab in the future.
Footnotes |
† Present address: Fachbereich Physik, Freie Universität Berlin, Arnimallee 14, 14195 Berlin, Germany. |
‡ Present address: Vienna Center for Quantum Science and Technology (VCQ), Faculty of Physics, University of Vienna, Boltzmanngasse 5, Vienna A-1090, Austria. |
This journal is © the Owner Societies 2019 |