ORIGINAL RESEARCH article

Front. Chem., 17 September 2021
Sec. Physical Chemistry and Chemical Physics
https://doi.org/10.3389/fchem.2021.746229

The Influence of the Chemical Potential on Defects and Function of Perovskites in Catalysis

www.frontiersin.orgGregor Koch1, www.frontiersin.orgMichael Hävecker2, www.frontiersin.orgPierre Kube1, www.frontiersin.orgAndrey Tarasov1, www.frontiersin.orgRobert Schlögl1,2 and www.frontiersin.orgAnnette Trunschke1*
  • 1Department of Inorganic Chemistry, Fritz-Haber-Institut der Max-Planck-Gesellschaft, Berlin, Germany
  • 2Max Planck Institute for Chemical Energy Conversion, Heterogeneous Reactions, Max-Planck-Gesellschaft, Mühlheim, Germany

A Sm-deficient Sm0.96MnO3 perovskite was prepared on a gram scale to investigate the influence of the chemical potential of the gas phase on the defect concentration, the oxidation states of the metals and the nature of the oxygen species at the surface. The oxide was treated at 450°C in nitrogen, synthetic air, oxygen, water vapor or CO and investigated for its properties as a catalyst in the oxidative dehydrogenation of propane both before and after treatment. After treatment in water vapor, but especially after treatment with CO, increased selectivity to propene was observed, but only when water vapor was added to the reaction gas. As shown by XRD, SEM, EDX and XRF, the bulk structure of the oxide remained stable under all conditions. In contrast, the surface underwent strong changes. This was shown by AP-XPS and AP-NEXAFS measurements in the presence of the different gas atmospheres at elevated temperatures. The treatment with CO caused a partial reduction of the metals at the surface, leading to changes in the charge of the cations, which was compensated by an increased concentration of oxygen defects. Based on the present experiments, the influence of defects and concentration of electrophilic oxygen species at the catalyst surface on the selectivity in propane oxidation is discussed.

Introduction

The A or B positions in the crystal lattice of ABO3 compounds with perovskite structure can be occupied by almost all metallic elements of the periodic table according to Goldschmidt’s rule (Goldschmidt, 1926). In this way, a wide variety of compounds is obtained. Other materials with special properties are accessible by partial substitution of either A or B leading to AA′BB′O3 formulas (Kurbakov, 2010; Zhu et al., 2015; Yang et al., 2016; Bartel et al., 2019). Chemical substitution is generally also accompanied by changes in defect structure. For example, the replacement of an A3+ with A2+ leads to lattice distortion (Kurbakov, 2010) and the formation of oxygen vacancies. Alternatively, partial oxidation of the B element may occur if this is possible (Ciambelli et al., 2000; McFarland and Metiu, 2013; Toniolo and Schmal, 2016). Both phenomena can also occur simultaneously. Defect formation can be achieved not only by the partial or complete substitution of A or B. Thus, reduction of the available amount of A or B also provides access to defect-rich compounds A1-xB1-yO3 that differ significantly from the ideal stoichiometry ABO3 (Golikov et al., 2005; Neagu et al., 2013). Due to the dense packing of ions in the perovskite structure, the occupation of interstitial sites is very unlikely. Therefore, for defect formation, only vacancy formation or migration of smaller B elements to vacant A positions, which offer more space, can be considered (Wołcyrz et al., 2003; Ignatans et al., 2019).

Defect-rich structures are of particular importance with respect to the application of perovskites in heterogeneous catalysis (Royer et al., 2014; Najjar and Batis, 2016). In addition to the increasing interest in perovskites as catalysts for oxygen reduction and oxygen evolution reactions (Suntivich et al., 2011; Gupta et al., 2016; Zhu et al., 2019), this class of materials is also widely used in photocatalysis (Grabowska, 2016; Tasleem and Tahir, 2020) or total oxidation to remove volatile organic and hazardous compounds (Gil et al., 2004; Liu et al., 2018a; Liu et al., 2018b; Liu et al., 2019). However, perovskites are also suitable catalysts for selective oxidations at low temperatures (Polo-Garzon and Wu, 2018; Kamata, 2019), especially for the oxidative dehydrogenation of propane (Koch et al., 2020). The catalytic properties can be influenced by various factors, such as the surface composition (Fierro and Tejuca, 1987; Koch et al., 2020), the concentration of defects, or the partial reduction of metal ions (Evarestov et al., 2005; Mierwaldt et al., 2014; Ignatans et al., 2019; Koch et al., 2020). Oxygen species at the catalyst surface are responsible for the activation of the propane molecule. X-ray Photoelectron Spectroscopy (XPS) allows the identification of different oxygen species in the surface region of perovskites (Mierwaldt et al., 2014; Stoerzinger et al., 2014; Stoerzinger et al., 2015). Furthermore, the dynamics of these species can be studied as a function of the chemical potential of the gas phase under reaction conditions of the selective oxidation of propane using ambient pressure (AP) XPS (Koch et al., 2020).

Mn-based perovskites are generally catalyzing the total combustion of propane already at low reaction temperatures. However, the addition of water vapor to the feed of propane and oxygen leads to higher selectivity of the desired selective oxidation product propene. The improved selectivity is apparently related to the partial reduction of surface manganese ions and the formation of defects and hydroxides (Koch et al., 2020). To examine these relationships in more detail, a samarium-deficient Sm0.96MnO3 perovskite was prepared by combustion synthesis and subjected to various post-treatments at 450°C in inert gas and gases containing H2O and CO. The compound was chosen because it is structurally stable under relevant reaction conditions (Kamata et al., 1979; Hui et al., 1990; Golikov et al., 2005) and exhibits good catalytic activity at low temperatures as well as significant changes in catalytic performance upon variation of the feed composition. In contrast, LaMnO3 (Koch et al., 2020), which is often studied in oxidation reactions, is more sensitive to changes in oxygen partial pressure. This can ultimately lead to defect restructuring associated with the transformation of crystal symmetry from rhombohedral to orthorhombic and thus with the change of bulk structure (Miyoshi et al., 2003).

The catalyst was examined before and after the treatments in the different gas phases with respect to its catalytic performance, defect structure and surface composition. AP-XPS experiments and ambient pressure Near Edge X-Ray Absorption Fine Structure (AP-NEXAFS) measurements in the presence of the gas phase were performed to study in detail the changes in the O 1s spectra and the surface oxidation state of the Mn and Sm ions as a function of the pretreatment conditions and the feed composition. Based on the spectroscopic results, possible relationships between the surface structure of the perovskite and its properties as a catalyst in propane oxidation are discussed.

Materials and Methods

Synthesis

A self-igniting gel was first prepared by the Pecchini method (Pechini, 1967). Glycine was used as a fuel in a stoichiometric ratio to the metal nitrate concentration, so that combustion occurred after ignition of the resulting gel without the participation of oxygen from the gas phase (Cooper and Kurowski, 1996; Hwang et al., 2005; Koch et al., 2020). At first, 13.74 g of Sm(NO3)3·6H2O (Aldrich, 99.9% trace metal basis, LOT: MKCG 1993) were dissolved in 75 ml H2O. Afterward, 8.17 g of Mn(NO3)2·4H2O (ROTH, 98%, LOT: 46827460) were added, yielding a slightly yellow solution. While stirring 6.27 g of glycine (Lancester Synthesis, Batch 10052309) was added. After stirring for 1 h, the solution was transferred into an evaporation basin placed on a hot plate and the water was evaporated at 120°C. The formed resin was self-ignited by heating to 350°C. The dark powder (ID 29987) was obtained with a yield of 82%. On basis of thermogravimetry, the optimal calcination temperature was determined (Supplementary Figure S1). For calcination, 6.27 g of combustion product was put into a combustion boat which was placed in a tube furnace. While streaming in 200 ml min−1 a mixture consisting of 20% O2 in Ar, the powder was heated to 800°C at 5°C min−1, calcined at this temperature for 5 h, and then cooled to 25°C. The calcination yield was 91%, which corresponded to an overall yield of 76%. The catalyst ID is 30036.

X-Ray Diffraction

X-ray diffraction (XRD) patterns were measured on a Bruker D8 ADVANCE II theta/theta diffractometer, using Ni-filtered Cu Kα radiation and a position-sensitive LynxEye silicon strip detector. Structure parameters were calculated by a least-square fitting of the diffraction patterns by using corresponding structure models taken from the ICSD database. Refinement was performed utilizing the program package TOPAS (version 4.2, copy-right 1999–2009 Bruker AXS).

X-Ray Fluorescence

The metal content of the perovskite catalysts was determined by X-ray fluorescence spectroscopy using a Bruker S4 Pioneer wavelength dispersive X-ray fluorescence spectrometer. For sample preparation, a mixture of 0.05 g of the catalyst and 8.9 g of lithium tetraborate (>99.995%, Aldrich) was fused into a disk using an automated fusion machine (Vulcan 2 MA, Fluxana).

Thermogravimetry Coupled with Mass Spectrometry

The weight loss of the powder was measured using a NETZSCH TG 209F1 Libra thermo balance and the gas phase composition was determined using a Pfeiffer Omnistar mass spectrometer. 11.782 mg of the sample were put into an Al2O3 crucible, which was then placed on the sample holder. After equilibration time, the sample was heated to 900°C at 10°C min−1 while a gas flow of 70 ml min−1 streamed through the heating chamber. The gas was composed of 21% O2 in Ar. Previously measured correction without sample was considered in data analysis.

The same setup was used for the TG-MS experiment simulating the CO treatment. 11.873 mg of the calcined sample were loaded in an Al2O3 crucible. After equilibrating, the calcined sample was heated to 450°C at 10°C min−1 followed by a dwelling period of 20 min. The gas consisted of 21% O2 in Ar and streamed with 70 ml min−1 to remove all surface carbonates possibly formed. Subsequently, the atmosphere was changed to 2% CO in Ar at the same flow rate and the temperature program was repeated.

N2 Physisorption

N2 adsorption isotherms were measured at −196°C using the Autosorb-6B analyzer (Quantachrome) after outgassing the catalysts in the vacuum for 2 h at 150°C. Data analysis was done using the Quantachrome Autosorb software package. The specific surface area was calculated according to the multipoint Brunauer-Emmett-Teller method (BET) in the range 0.05 < p/p0 < 0.15 assuming the N2 cross-sectional area of 16.2 Å2.

Scanning Electron Microscopy and Energy Dispersive X-Ray Fluorescence Spectroscopy

The sample was fixed on carbon tape on the top of an Al-sample stage before the transfer into the SEM load chamber. SEM images were conducted using a Hitachi S-4800 microscope equipped with a cold field emission gun. For imaging, 1.5 kV acceleration voltage and 4 mm working distance was set to display the image by using both, upper and lower secondary electron detector. Additionally, EDX mapping was done setting the acceleration voltage to 15 kV and working distance to 10 mm. An energy-dispersive QUANTAX 800 EDX spectrometer working with a XFlash®6│30 detector was used to map the elemental distribution. Usually, mapping was done for 2 min when sufficient counts were gained (>10 kps). Recorded spectra were analyzed considering all found peaks and corresponding elemental ratios were calculated for each map/spectrum. Standard deviations were calculated on the results of all maps.

Catalytic Testing

The catalytic tests were carried out as described before (Koch et al., 2020) using a setup for partial oxidation (Integrated Lab Solutions) with 10 fixed bed tubular reactors in parallel. All reactors have an inner diameter of 2 mm and all reactors were equipped with a thermocouple placed inside the catalyst bed. The catalytic performance was determined at atmospheric pressure under steady-state conditions. The reactant feed consisted of 5% propane and 10% O2 in N2 (C3H8/O2/N2 = 5/10/85). Furthermore, 40% H2O was added to the feed using an evaporator (C3H8/O2/N2/H2O = 5/10/45/40). Before testing, the catalysts were pressed with a force of 2 t for 1 min, crushed, and sieved to particles, which sizes ranged from 100 to 200 µm. Generally, a total flow of 15 ml min−1 was set. Mass of catalysts of 52, 96 and 142 mg were loaded into single steel reactors to achieve contact times of W/F = 0.22 g∙s∙ml−1, W/F = 0.38 g∙s∙ml−1 and W/F = 0.57 g∙s∙ml−1.

An online gas chromatograph (Agilent 7890A) was used for gas analysis. A combination of Plot-Q (length 30 m, 0.53 mm internal diameter, 40 µm film thickness) and Plot-MoleSieve 5 A columns (30 m length, 0.53 mm internal diameter, 50 µm film thickness), connected to a thermal conductivity detector (TCD), was used to analyze the permanent gases CO, CO2, N2, O2, and CH4. A system of an FFAP (length 30 m, 0.53 mm internal diameter, 1 µm film thickness) and a Plot-Q column (length 30 m, 0.53 mm internal diameter, 40 µm film thickness) connected to a flame ionization detector (FID) was used to analyze the C2-C3 hydrocarbons.

Conversion of propane XC3H8, and selectivity Si of product i in percentage, were calculated based on the carbon number and the sum of all products formed:

XC3H8=i=1nNicii=1nNici+3cC3H8,out100
Si=Nicii=1nNici100,

where Ni is the number of carbon atoms in the product i, ci is the concentration of the product i in the reactor exit gas, and  cC3H8,out  is the concentration of propane in the exit gas. The conversion of oxygen in percentage was calculated based on the difference between the oxygen concentration in the reactant feed cO2,in and in the exit gas cO2,out:

XO2=cO2,incO2,outcO2,in100.

Integral reaction rates rpropane normalized to the specific surface area in µmol∙h−1∙m−2 for propane consumption were determined based on the first derivation of the linear function

XC3H8=f(WF),

resulting in:

rpropane=dXC3H8d(WF)S,

where XC3H8 is used not in percentage, but as a value between 0 and 1, W is the catalyst mass in g, F is the flow rate of propane in µmol∙h−1, and S is the specific surface area in m2∙g−1.

Initial formation rates rinitial of propene and CO2 were calculated through extrapolation. Corresponding formation rates at each contact time were fitted with a linear function. The function was used to calculate the initial rates at 0% conversion of propane.

Ambient Pressure Photoelectron Spectroscopy and Ambient Pressure Near-Edge X-Ray Absorption Spectroscopy

Spectra were collected according to the procedure described in (Koch et al., 2020). AP-XPS and AP-NEXAFS measurements were performed at the ISISS (Innovative Station for In Situ Spectroscopy) facility operated by the MPG at the synchrotron radiation source BESSY II in Berlin, Germany. 12.2 mg of sample was pressed into a pellet of 8 mm in diameter. The pellet was mounted on a steel plate and covered by a steel cover lid having a 4 mm central aperture. The sample was connected to a thermocouple and then the sample holder was loaded to the instrument as described previously (Salmeron and Schlögl, 2008; Knop-Gericke et al., 2009).

Firstly, spectra were collected at room temperature in a vacuum, then heated in a stream of 5 ml⋅min−1 atmosphere consisting of propane/O2/He/H2O of 0/20/80/0 (heating rate 5°C min−1) to 270°C. All gases and water vapor were dosed to the XPS cell via calibrated mass flow controllers (BRONKHORST). The gas composition was then switched consecutively to 0/10/90/0, and 5/10/45/40 before switching to 25% CO in He. This was followed by heating to 450°C at 5°C min−1. After spectra collection, the sample was cooled to 270°C at 5°C min−1 again and the atmosphere was set to 5/10/45/40, and then to 5/10/85/0. The total gas pressure at all reaction conditions was 0.25 mbar.

AP-XPS spectra were collected with an exit slit width of the beamline of 111 μm. The pass energy of the hemispherical electron analyzer was 20 eV for all core level spectra. The photon energy scale of the beamline has been adjusted by measuring the Fermi edge of a Pd foil at various photon energies. The measurements were carried out at a constant photoelectron kinetic energy of ≈170 eV for all core levels, examining the outermost surface layer with an approximate thickness corresponding to an approximate electron inelastic mean free path IMFP of λ = 0.6 nm. For this purpose, the photon energies were set to 255 eV (Mn 2s), 810 eV (Mn 2p), 700 eV (O 1s), 455 eV (C 1s), 1,260 eV (Sm 3d), 315 eV (Sm 4d) and 170 eV (valence band). Hence, the transmission of electrons through the gas phase and in the lens of the electron spectrometer was the same for all elements and does not need to be included in the quantitative calculation of element abundances. Additionally, the subsurface region of the material was studied by using the larger kinetic energy of ≈770 eV using photon energies of 1,300 eV (O 1s) and 1,055 eV (C 1s) that corresponds to an IMFP of λ = 1.6 nm. These values of λ were estimated using the predictive TPP-2M formula by Tanuma, Powell and Penn (Tanuma et al., 1991; Tanuma et al., 2003) (Powell and Jablonski, 2010) (NIST Electron Inelastic-Mean-Free-Path Database Version 1.2). For the quantitative analysis of the core level intensities, the atomic subshell photoionization cross sections and asymmetry parameters from numerical calculations by Yeh and Lindau (Yeh and Lindau, 1985; Yeh, 1993) were used taking the monochromatic photon energy-dependent photon flux into account. The core-level spectra were deconvoluted using Gaussian-Lorentzian product functions after subtracting a Shirley background with the CasaXPS software (Neal Fairley, Version 2.3.15, © 1999–2009 Casa Software Ltd.). Gas-phase peaks have been suppressed during the AP-XPS measurements by applying a bias of +90 V to the entrance aperture of the spectrometer.

AP-NEXAFS data processing included averaging two spectra of each edge followed by background subtraction and normalization. Total electron yield was used for Mn L-edge spectra and Sm M-edge spectra. Auger electron yield was used for O K-edge spectra with the analyzer set to a kinetic energy of Ekin = 385 eV and a pass energy of Epass = 50 eV.

Results

Phase-pure Sm0.96MnO3 (Supplementary Figure S2, Table 1) ICSD# 95491 was prepared by the self-combustion synthesis route (Mori et al., 2002). The catalyst exhibited pores and voids of different sizes due to the high amounts of gases formed during ignition synthesis. The pores were retained during calcination. Thermal treatment of the synthesized material in oxygen is generally required to remove remaining organic components and carbonates formed by storage in air (Supplementary Figure S1). During the calcination of the combustion product, small crystallites were formed that adhere to each other and creating cavity-rich macrostructures (Supplementary Figure S3). The predominant size fraction of the nanocrystals was about 50 nm according to SEM, which is in good agreement with the calculated crystallite size from XRD refinement (51 nm, Table 1). However, larger and smaller crystallites are also present, but only in small fractions (Supplementary Figure S3). Due to the irregular macrostructure and the voids, a specific surface area of 7.6 m2⋅g−1 was obtained (Table 2). The desired non-stoichiometry of Sm and Mn in the volume was reached by using a reduced mass of the samarium precursor (Table 2). However, the surface of the catalyst showed a 1:1 stoichiometry of Sm and Mn (Table 2). Analysis of the volume revealed a slight oxygen deficit of the compound (Sm0.96MnO2.93, Table 2). Therefore, the segregation of nano-sized binary oxides seems unlikely (Golikov et al., 2005; Neagu et al., 2013). Segregated phases were also not detected by XRD (Supplementary Figure S2) and EDX mapping (Supplementary Figure S4).

TABLE 1
www.frontiersin.org

TABLE 1. Lattice parameters obtained form XRD analysis of fresh and spent Sm0.96MnO3 catalysts.

TABLE 2
www.frontiersin.org

TABLE 2. Surface area and elemental composition of Sm0.96MnO3 obtained from bulk analysis (XRF and O-content) and surface analysis (XPS and NEXAFS).

The Sm0.96MnO3 catalyst is oxidizing propane already at temperatures below 300°C. (Merino et al., 2005; Koch et al., 2020). In the present work, the reaction was performed in either dry feed, i.e. 5% propane and 10% O2 in nitrogen (5/10/85), or wet feed, i.e. 5% propane, 10% O2 and 40% H2O in nitrogen (5/10/45/40) (Koch et al., 2020). The dry feed at 300°C was chosen as a reference point to compare the influence of the different pretreatments on activity and selectivity (Supplementary Figure S5). The selectivity to propene measured at the reference points as a function of propane conversion is shown in Figure 1A for different contact times. Regardless of the different pretreatment conditions, represented by different colors, and regardless of the contact time, represented by full, half-filled and empty symbols, all the measurement points lie on the same S-X curve. Thus, in dry feed, it can be seen that the different pretreatments only affected the integral rates. Compared to the conversion after calcination in air at 450°C, the conversion decreased slightly after treatment in pure N2 at the same temperature. This loss of activity could not be reversed by subsequent treatment in pure oxygen and further treatment in synthetic air at 450°C. The decreased conversion could be caused by either slight sintering or surface modifications. An increase in the apparent activation energy from 61 ± 6 kJ⋅mol−1 to 76 ± 4 kJ⋅mol−1 after N2 pretreatment suggests more for the latter (Supplementary Table S1). After the treatments in 10% water vapor and 2% CO, the activity decreases more significantly.

FIGURE 1
www.frontiersin.org

FIGURE 1. Selectivity to propene in (A) dry feed (C3H8/O2/inert = 5/10/85) and (B) in wet feed (C3H8/O2/inert/H2O = 5/10/45/40) after specific pretreatments at 450°C. Additionally, in (B) selectivity at different conversions in dry feed in the temperature range from 250 to 310°C is depicted for comparison (red data points). F/W = 0.29 L⋅min−1∙g−1 (full symbols); F/W = 0.16 L⋅min−1∙g−1 (half-filled symbols); F/W = 0.11 L⋅min−1∙g−1 (empty symbols).

Using wet feed, the selectivity to propene remarkably increased significantly (by about 18–83% over the entire range of measured conversions) at comparable conversion (Figure 1B) indicating a change in the properties of the catalyst probably related to the addition of water vapor. An increase in selectivity in the wet feed has been observed previously for Mn-based perovskites that had a nearly balanced Mn/A element ratio or an excess of Mn at the surface (Koch et al., 2020). In contrast, an Mn deficiency on the surface resulted in reduced selectivity (Koch et al., 2020). This effect was particularly pronounced after the pretreatments in steam and CO at 450°C, and the selectivity improvement was more evident at low propane conversion (Figures 2A,B). Apparently, both steam treatment and reducing conditions during pretreatment promote the formation of a more selective catalyst under reaction conditions, but only when propane oxidation is carried out in the presence of H2O in the feed. When switched to dry feed, the catalyst returns to its non-selective mode (Figure 1A and Figures 2A,B).

FIGURE 2
www.frontiersin.org

FIGURE 2. Selectivity to propene shown at the same conversion of propane: 4% (A,C) and 8% (B,D). Comparison between dry and wet feed (A,B) and initial wet feed (after calcination in air at 450°C) and wet feed after the different treatments (C,D). C3H8/O2/inert = 5/10/85 (dry feed); C3H8/O2/inert/H2O = 5/10/45/40 (wet feed); F/W = 0.29 L min−1∙g−1; T = 300°C. The numbers above the columns indicate the respective increase in selectivity.

The used catalyst showed the same bulk structural properties as the fresh catalyst (Table 1, Supplementary Figure S2). Only minor changes in lattice parameters were observed, which did not affect the total unit cell volume. Therefore, an extensively increased formation of oxygen vacancies in the bulk could be excluded, because this should have led to lattice expansion and increased cell volume (Miyoshi et al., 2003). Moreover, the excess of Mn in the bulk remained after catalytic testing (data not shown). Recrystallization or crystallite growth did not occur at these low temperatures applied because the crystallite sizes before and after catalytic testing were similar.

In contrast to the stable bulk structure, the surface structure changed significantly when the chemical potential was varied, which was achieved both by changing the feed and varying the pretreatment atmosphere. To understand the significant increase in selectivity to propene in the wet feed after pretreatment with carbon monoxide (Figure 1B and Figures 2C,D), the surface was analyzed under reaction conditions.

Fitting the Mn L3,2-edges using a linear combination of relevant reference spectra for Mn2+, Mn3+ and Mn4+ species in the perovskite structure (Koch et al., 2020) showed that the oxidation state of manganese in the surface region changed under the different conditions (Figure 3). The dominant oxidation state of manganese in perovskites is Mn3+, accompanied by smaller proportions of Mn2+ and Mn4+ (Figure 4). (Levasseur and Kaliaguine, 2008; Liu et al., 2018b; Liu et al., 2019) Mn2+ species were already present in the fresh Sm0.96MnO3 catalyst (Koch et al., 2020) balanced by about 5% Mn4+ resulting in an average oxidation state of 3.03. After heating in 20% O2, the average oxidation state decreased to 3.00. Switching to wet feed resulted in a further reduction. All Mn4+ vanished, and the fraction of Mn2+ increased to 4% balanced by Mn3+. While treating the sample in CO at 450°C, the Mn2+ fraction was growing to 5% leading to an average oxidation state of 2.95. Having switched back to wet feed, the Mn2+ fraction decreased to the level reached in the wet feed before the CO treatment at higher temperatures and dropped even further after switching back again to dry feed. Hence, the partial reduction of Mnn+ species at higher temperatures was caused by the gas phase composition and is a reversible phenomenon of the Sm0.96MnO3 perovskite catalyst (Koch et al., 2020).

FIGURE 3
www.frontiersin.org

FIGURE 3. Measured and normalized TEY Mn L3-edges (black) and corresponding fitted spectra (red) using a linear combination of reference spectra of Mn2+ (pink), Mn3+ (purple) and Mn4+ (brown) (Koch et al., 2020). The conditions during the measurements at the synchrotron facility are indicated in each diagram in terms of temperature and feed compositions C3H8/O2/He or C3H8/O2/He/H2O.

FIGURE 4
www.frontiersin.org

FIGURE 4. Surface characteristics of Sm0.96MnO3 at 0.25 mbar at different temperatures and different feeds (C3H8/O2/He/H2O) as given on the x-axis. Bottom: Fraction of Mn2+ and Mn4+ of all Mn ions measured with NEXAFS. Values obtained from NEXAFS Mn L2,3-edge fitting. Mn3+ balances the rest. Middle: Corresponding average oxidation state n of Mnn+. Top: Mn fraction obtained from Mn 2p and Sm 3d XPS peak area analysis.

Under the more reducing conditions in CO, also Sm3+ was slightly reduced indicated by a small peak at 1,072 eV in the Sm 3d5/2 spectrum (Figure 5) (Krill et al., 1980; Brunckova et al., 2019). This peak was observed after treatment in CO in wet feed and also still in the following dry feed (red spectrum in Figure 5). Its intensity decreased in dry feed compared to wet feed indicating that the process was reversible. But the reoxidation of the partially reduced A, i.e., Sm, was slow (Figure 5, inset) compared to the fast reoxidation of Mn2+. However, the reduction of Sm3+ was only detected by AP-XPS and not by AP-NEXAFS at the Sm M5,4-edges (Supplementary Figure S6) showing that only the outermost surface fraction of Sm3+ was reduced. Since the electron kinetic energy in the AP-XPS experiments was with 170 eV distinctly smaller compared to the primary kinetic energy of the electrons created in the XAS process at an excitation energy in the NEXAFS experiment (about 800 eV), the resulting depth of information was much smaller in the Sm 3d AP-XPS measurements compared to the AP-NEXAFS measurement at the Sm M-edges. In the NEXAFS spectra at the Sm M-edges, the small surface fraction of reduced Sm2+ species was apparently hidden by the larger fraction of Sm3+ species in the deeper layers.

FIGURE 5
www.frontiersin.org

FIGURE 5. Normalized Sm 3d5/2 spectra of Sm0.96MnO3 collected at different conditions given as C3H8/O2/He/H2O (initial–black, 10% O2/He–green, wet feed–blue, dry feed–red). Inset depicts selected spectra ranges to show the Sm2+ peak.

The partial reduction of Mn and Sm requires a charge compensation and the redox processes should be reflected in both the O 1s spectra and the O K-edge spectra. First, the O 1s spectra are discussed. In Figure 6 the measured O 1s and the corresponding fitted O 1s spectra are shown for the wet feeds and the CO treatment as surface scans (A), (B) and (C) as well as deep scans (D) and (E). Remarkably, the fraction of each oxygen species, which are introduced and discussed below, changes when comparing the outermost surface region (Figures 6A–C) with the deeper scans (Figures 6D,E). To describe all O 1s spectra measured under different conditions, different types of oxygen species should be considered according to the work published so far (Koch et al., 2020). The assignment and the fit model are summarized in Table 3. Olattice represents the oxygen positioned in regular coordination in the crystal lattice. Odefect at 530.1 eV is oxygen, which is located in the vicinity of a defect (vacancy, off-centered position, hydroxyl group, etc.) (Fan and Goodenough, 1977; Sunding et al., 2011; Mierwaldt et al., 2014). Here it is noteworthy to mention that the oxygen atoms occupy two different positions in the orthorhombic lattice. Both oxygen atoms have each built up a distorted octahedral polyhedron consisting of two diametrically placed Mn atoms and four additional Sm atoms (Figure 7). These two octahedra differ slightly in their distortion. The distances between each oxygen position and its metal neighbors vary slightly, but the average distances between oxygen atoms and Mn and Sm atoms are similar (Supplementary Table S2). Since the two different oxygen species each have a different local environment, they should be spectroscopically distinguishable. However, it is not clear whether these differences can be resolved experimentally. This requires further experiments and the support of theory. From the available experiments, there is no clear evidence to further distinguish between an Olattice1 and an Olattice2, which lay below the signal of Olattice and could also overlap with the signal of Odefect. Therefore, according to the literature, only the Olattice and Odefect species are considered to describe the volume of the perovskite in the following discussion. The species Osurface at 530.9 eV is attributed to an electrophilic oxygen species adsorbed on the surface (Merkle and Maier, 2008; Mierwaldt et al., 2014; Stoerzinger et al., 2014). The binding energy 531.2 eV was uniquely assigned to Ocarbonate because the peak intensity of this signal correlates with the peak intensity of Ccarbonate from C 1s spectra at a binding energy of 288–289 eV (Stoerzinger et al., 2014; Koch et al., 2020). In addition, OOH species are assigned to the binding energy at 531.8 eV (Nesbitt and Banerjee, 1998; Banger et al., 2011). Adsorbed H2O and adsorbed oxygenates at binding energies of 532.2 eV and 534–535 eV, respectively, are present at higher temperatures only with a fraction of less than 1%.

FIGURE 6
www.frontiersin.org

FIGURE 6. Measured O 1s spectrum (black), fitted spectrum (red) and background (black dots) measured in wet feed (C3H8/O2/He/H2O = 5/10/45/40) at 270°C (A,C–E) and in 25% CO balanced by He at 450°C (B). Fitted peaks of O species are depicted as follows: Olattice (olive) Odefect (green), Osurface (magenta), Ocarbonate (orange), and OOH (blue). (A–C) were measured with 700 eV excitation energy representing the surface region (surface). (D,E) were measured with 1,300 eV excitation energy representing information of deeper regions (deep).

TABLE 3
www.frontiersin.org

TABLE 3. Energy positions (BE) of oxygen species in O 1s spectra and applied peak width (FWHM) for fitting all O 1s spectra (Koch et al., 2020).

FIGURE 7
www.frontiersin.org

FIGURE 7. Representation of the two different lattice oxygen (red balls) O1 (A) and O2 (B) centered between next neighboring Mn atoms (purple balls) and Sm atoms (yellow balls). Numbers next to the illustrated connection between oxygen and metal atom (orange) give the distance between the atoms. Corresponding octahedral polyhedra are depicted on the right.

Figure 8 and Figure 9 show the evolution of the fraction of oxygen species as a function of the conditions applied. Two photon energies were used for measuring the O 1s spectra. At 700 eV (inelastic mean free path IMFP = 0.6 nm), mainly the outermost surface is analyzed. At 1,300 eV (IMFP = 1.6 nm), deeper surface layers are also detected. In Figure 8 the depth profile for Olattice and Odefect is shown. In general, the amount of Olattice was larger in deeper layers than at the outermost surface. With increasing temperature, the fraction of Odefect increased, especially at the outermost surface (Figure 8), which supports the plausibility of the applied fit model.

FIGURE 8
www.frontiersin.org

FIGURE 8. Evolution of the fraction of Olattice (top) and Odefect (bottom) at different temperatures and feeds (C3H8/O2/He/H2O) given below. Deep, filled columns: Extracted from O 1s peak fitting of spectra conducted using 1,300 eV excitation energy. Surface, hatched columns: Extracted from O 1s peak fitting of spectra conducted using 700 eV excitation energy. Dark red, surface, hatched columns: Extracted from O 1s peak fitting of spectra conducted at 450°C using 700 eV excitation energy.

FIGURE 9
www.frontiersin.org

FIGURE 9. Evolution of the fraction of Osurface (magenta) and OOH (blue) extracted from O 1s peak fitting at different temperatures and feeds (C3H8/O2/He/H2O) given below. Hatched columns represent fractions at 450°C, filled columns represent fractions at 270°C.

After heating up in synthetic air to 270°C and subsequently reducing the oxygen content from 20 to 10%, the surface of the catalyst was healed and exhibited the largest fraction of Olattice. Minor fractions of Odefect, Osurface, and OOH were also observed indicating complex surface structures. Carbonates were already largely decomposed at the reaction temperature, so their contribution to the O 1s spectrum could be neglected. Remarkably, the increase in Odefect fraction after heating indicates that defect-rich surface structures have already been generated by the increased temperature and/or reduced oxygen partial pressure. This is consistent with model calculations that have shown that defect-rich surfaces of LaMnO3 are energetically favored (Evarestov et al., 2003; Evarestov et al., 2004). The change in chemical potential by switching to wet feed caused an increase in the concentration of Odefect, OOH and Osurface at the expense of Olattice (Figure 8 and Figure 9). Hence, either the H2O in the feed or the alkane and its derivatives or both actively modified the surface of the catalyst. After treatment with CO at 450°C and decreasing the temperature in CO to 270°C, a further little decrease of Olattice and a distinct decrease of Osurface and OOH was observed. At the same time, the fraction of Odefect increased. This is consistent with the partial reduction of Mn3+ and Sm3+ observed using AP-NEXAFS and AP-XPS, respectively (Figure 4 and Figure 5). The surface oxygen vacancies could have been formed by the oxidation of CO to CO2 using surface oxygen as an oxidant (Hui et al., 1993). This was confirmed by repeating the CO treatment in a thermogravimetry experiment using a mass spectrometer for gas analysis. Only the formation of CO2 was observed while heating up to 450°C in 2% CO balanced by Ar (Supplementary Figure S7). Switching to the wet reaction feed caused an increase in Osurface and OOH and a simultaneous slight decrease in Odefect. The fraction of Olattice at the outermost surface remained unchanged. Comparing the wet feed before CO treatment with the wet feed after CO treatment, the amounts of Osurface and Olattice in the surface region were slightly lower after CO treatment, but the Odefect concentration was clearly higher, while the amount of OOH remained similar (Figure 8 and Figure 9). Thus, the distribution of oxygen species in the surface region resembled the state before CO treatment, except for Odefect species, which was present in higher concentrations. The higher concentration of Odefect is supported by the thermogravimetry in CO-containing atmosphere. The found mass loss of about 0.6% corresponds theoretically to the removal of all oxygen atoms linked to a Mn or a Sm of the outermost surface layer (Supplementary Figure S8, right). This would imply the reduction of the surface to the metallic state by removal of all oxygen atoms of the outermost surface layer (Supplementary Figure S8, left), which is not found spectroscopically. Only 0.2% mass loss is estimated assuming a 001 surface showing a balanced A/B ratio as found for this catalyst (Table 2) and assuming fully saturated with oxygen atoms (Supplementary Figure S8, left). This means also oxygen from the bulk or at least from the near-surface range reacts with CO. Hence, more oxygen vacancies are likely formed and thus the number of Odefect increases, which is oxygen in the vicinity of a defective site. By removing H2O from the feed, the Olattice fraction in the deep and surface areas increased slightly again (Figure 8). At the same time, the proportion of Odefect on the surface decreased slightly. This decrease was more evident in the deep region, indicating an ongoing healing process in the bulk. Moreover, the amount of Osurface increased and the amount of OOH slightly decreased (Figure 9). According to data measured in this study, Osurface appears to be required to start the healing process, reduce the amount of Odefect and increase the amount of Olattice again.

Surface modifications and the increase of Odefect species were confirmed by near-edge X-ray absorption spectroscopy at the O K-edge (Figures 10A,B). In transition metal oxides containing 3d B metals, the first absorption above the threshold at ∼530 eV corresponds to the excitation from the O 1s state to the hybridized O 2p-B 3d state. The spectrum is difficult to interpret because the intensity and shape of this feature are related in a complex way to a multitude of overlapping and interconnected phenomena, such as oxidation state, spin state, degree of hybridization, strain and defects in the oxide structure. The 3d band of Sm0.96MnO3 is split into two peaks at 529.5 and 530.8 eV, which are in the simplified picture of the ligand field theory assigned to Mn eg and t2g states, respectively (Mierwaldt et al., 2014; Frati et al., 2020). In the spectra shown here (Figures 10A,B) the positions match well to those reported. The eg peak is invariant at 529.5 eV while the t2g peak maximum shifts from 530.5 eV in dry feed to 530.7 eV in wet feed after CO treatment. When the chemical potential of the atmosphere was changed, a change in the onset (during CO treatment at 450°C, grey spectrum) and intensity of the pre-edge peak was clearly observed (Figure 10B left), indicating that charge transfer had occurred. Also, the shape of the peak changed drastically. The overall area of the pre-peak was integrated taking into account a step function for the ionization to the continuum (Figure 10B right). The area below the pre-peak ranges from 2.79 during CO treatment to 3.45 in wet feed after CO treatment. The O K-edge after the CO treatment in wet feed differs strongly from those measured in the wet feed before the CO treatment. Since the pre-edge peak represents holes in the band, the number of holes is consequently reduced when the catalyst is reduced during CO treatment, resulting in a small area of 2.79. Switching back to wet feed yields in a significant increase in the area to 3.45 which is clearly visible from the increased t2g peak at 530.7 eV (Figures 10A,B arrow). This number is also higher compared to 2.96 measured in the wet feed before CO treatment. Hence, the CO treatment tremendously modified the surface and the sub-surface of the catalysts in the way that the CO reduced partially the surface metals Mn and Sm (Figure 4 and Figure 5) and generated more Odefects (Figure 8), which is supported by thermogravimetry measurement. Finally, these surface modifications yielded in an increased number of holes in the t2g band when the water of the wet feed or one of the other feed ingredients or reaction intermediates come into play. After switching to the dry feed, the area below the pre-peak is with 3.17 still increased compared to 2.96 in the wet feed before, but it already indicates the slow healing.

FIGURE 10
www.frontiersin.org

FIGURE 10. (A) Normalized O K-edge spectra measured in different atmospheres and temperatures. (B) Corresponding pre-edge range of O K-edges. Left: enlarged overlay at pre-peak onset. Right: Pre-peak including the step function in grey dots representing the background. Resulting areas of pre-edge peak are given on the right. Arrows in (A,B) indicate the same position of most changing range in the pre-edge region. (C) Corresponding normalized valence band spectra measured in different atmospheres and temperatures. In addition, difference spectra are plotted in dots at the bottom. Enlarged difference spectra are depicted in the inset at the top left. Difference spectra show the difference to previous conditions. Conditions are given as C3H8/O2/He/H2O (initial–black, 10% O2/He–green, wet feed–blue, CO/He–grey, dry feed–red).

The indication of partial reduction of Mn3+ to Mn2+ in CO at 450°C was complementarily confirmed by Mn L2,3-edge spectra (Figure 4) and valence band spectroscopy (Figure 10C). In addition, adsorption of propane and propene and the formation of OH groups were evidenced by the appearance of peaks at 8.8 and 10.2 eV (Wernbacher et al., 2019). Thus, the reactants are not only converted on the catalyst surface but also modify the catalyst surface, which results in a reactant-dependent catalytic behavior.

Discussion

The optimized combustion method is suitable to achieve batches in gram scale of phase pure Sm0.96MnO3 powders having low bulk density due to the agglomerated crystallites forming caves, hollow spheres, and large networks (Supplementary Figure S3). The in this way prepared Sm0.96MnO3 can be operated as a catalyst in the oxidative dehydrogenation of propane (Koch et al., 2020). The Sm0.96MnO3 is characterized by a stable bulk crystal structure under the applied reaction conditions. Neither crystallite growth (Table 1) nor lattice expansion, which were observed for similar LaMnO3 perovskites (Miyoshi et al., 2003; Koch et al., 2020) occurred. Both show the stability of the lattice and imply that extensive re-organization of the lattice can be neglected.

However, the functionality of the catalyst changes in terms of both activity and selectivity when different pretreatments and wet feed are employed showing that the chemical potential affects the surface of the catalyst and influences, thus, the catalytic properties. The chemical potential reversibly modifies the oxidation state of metals in the surface region. Even heating to 270°C in an oxygen-containing atmosphere and then switching to wet feed leads to full reduction of all initially measured Mn4+ species to Mn3+ and also to an increased formation of Mn2+. This illustrates the importance of operando techniques, which can reveal the dynamic redox behavior of the Sm0.96MnO3 catalyst under various conditions. The partial reduction of Mnn+ in perovskites, e.g. in the oxygen evolution reaction (Mierwaldt et al., 2014) or the oxidative dehydrogenation of propane (Koch et al., 2020) has been related to H2O in the atmosphere and is feasible when Mnn+ occupies an A position (Ignatans et al., 2019; Koch et al., 2020) in the ABO3 perovskite structure. In the A-deficient Sm0.96MnO3 catalyst prepared here, some small Mnn+ ions probably sit on spacious A positions. The changes in oxidation state are balanced by types and amount of oxygen in the surface region. When the Mn is partially reduced, also the amount of Olattice/(Sm + Mn) is reduced showing that the partial reduction is accompanied by loss of oxygen atoms sitting on a perfect lattice position (Figure 11 grey and purple circles). Thus, the increase in Odefect is expected, which supports the applied fit model. In general, more oxygen (Oall is the sum of all species) was found using AP-XPS (1.66 < Oall/(Sm + Mn) < 2.98) than with the chemical analysis of the bulk (Obulk/(Sm + Mn) = 1.49), showing a highly defective surface. The chemical analysis cannot distinguish between oxygen in different configurations (e.g., Olattice, Odefect), but gives just the sum Obulk (oxygen determined by bulk chemical analysis). The Olattice content determined by surface analysis varies as a function of conditions as follows: 0.72 < Olattice/(Sm + Mn) < 1.56. It decreases to a minimum in the wet feed after treatment in CO at 450°C and slowly increases in the subsequent dry feed again (Figure 11). Hence, the amount of oxygen in the surface region in relation to A and B does not only depend on the type of B as previously discussed (Fierro and Tejuca, 1987) but also on the feed and the temperature. Interestingly, when only oxygen was present in the atmosphere at 270°C, bulk and surface analysis yielded the most similar figures for Oall/(Sm + Mn), namely 1.49 and 1.66, respectively, with Olattice having the largest fraction in the surface analysis. This was larger here than in all other conditions. However, the crystal lattice was not perfect, as the Olattice/(A + B) ratio of 1.18 in the outmost surface layer was lower than the expected ratio of 1.49 for a perfect Sm0.96MnO3 lattice (Figure 11, dotted line). This is because there are always unsaturated atoms on the surface. Therefore, only a perfect, defect-free lattice would result likely in a consistent bulk and surface composition. Conversely, it can be deduced that the over-stoichiometry of oxygen in the surface region already indicates additional oxygen species, with the interaction of this oxygen with the metals being lower than the average M-O interaction in the bulk. Noteworthy to mention, the additional oxygen is not coming simply from adsorbed H2O wetting the surface because the corresponding water signal at 533.2 eV (Table 3) decreases after heating to 270°C in all atmospheres due to high temperature (Figure 8). The oxygen excess on the surface is based on OOH, Osurface, a little amount of Ocarbonate and Odefect. Especially, the Odefect is the major non-lattice oxygen species. The amount of Odefect is increased in wet feed (Mierwaldt et al., 2014; Koch et al., 2020) and further increased by reductive treatment in CO. It slowly decreases in the wet feed again and decreases, even more, when dry feed is employed (Figure 8). Simultaneously, Mn and Sm are partially reduced and reoxidized, respectively. Moreover, the Odefect increases apparently at the expense of Olattice when the reaction feed streams over the catalyst (Figure 8). The oxygen can undergo hydrogenation through propane and its derivatives creating thus an OOH, also indicated in the difference plots of the valence bands (Figure 10C), when the oxygen is bound to an Mn (Li et al., 2019). In vicinity to this in situ formed hydroxyl group also Odefects are created causing thus the increase in surface concentration of Odefect species in the feed. These findings suggest that the creation of the Odefect is related to both hydroxylation and oxygen vacancy formation. Moreover, the peak increase in the pre-peak region of the O K-edge NEXAFS spectra at 530.8 eV in wet feed underlines the formation of defective structures. This is also reflected in the valence band measurements, which provide information on the occupied states (Figure 10C).

FIGURE 11
www.frontiersin.org

FIGURE 11. The ratio of the sum of all fitted oxygen (Oall) species and the sum of Mn + Sm obtained from XPS analysis (empty circles) in dependence on conditions applied and the ratio of the amount of Olattice and the sum of Mn + Sm obtained from XPS analysis (full circles) in dependence on conditions applied. Purple full circles depict the Mn oxidation state obtained from Mn L-edge fitting. C3H8/O2/inert = 5/10/85 (dry); C3H8/O2/inert/H2O = 5/10/45/40 (wet). The dotted line represents the number from bulk analysis for comparison (XRF and O content).

In conclusion, the assignment of the O 1s signal at 530.1 eV to oxygen, which is part of irregular surface structures (Odefect), is strongly supported by simultaneously observed charge transfer in the O K-edge spectra and the valence band spectra. Hence, Odefect represents the defective state of the surface and this defective perovskite surface catalyzes the oxidative dehydrogenation of propane more to propene compared to a less defective perovskite surface. Once the healing is finished, the initial catalytic properties of Sm0.96MnO3 are recovered, which is reflected in the similar performance in the dry feed after the different pretreatments (Figure 1A). Moreover, remembering that Mn3+ is partially reduced in wet feed but simultaneously more holes are observed indicating reactant-induced charge distribution in the outermost surface region. Therefore, the healing cannot be easily described as the filling of formed oxygen vacancies. Many processes cause charge transfer and may thus contribute to the catalytic properties of the Sm0.96MnO3 catalyst in oxidative dehydrogenation of propane.

Remarkably, switching from dry to wet feed results in an increase in the selectivity to propene. This is more pronounced after additional treatment at higher temperatures and more reducing chemical potentials. Moreover, it is reversible after switching to dry feed. The products formed in propane oxidation are only propene and CO2. This implies that possibly formed intermediates are fully oxidized under the conditions applied. Changes in selectivity can be caused by either enhanced formation of propene or by reduced formation of CO2. In order to compare the intrinsic selectivity of the catalyst after various treatments, the initial formation rates of propene and CO2 are compared (rates extrapolated to 0% conversion). The ratio between the initial propene formation rate and the initial CO2 formation rate is depicted in Figure 12A. The ratio is almost independent of the pretreatment when the dry feed is set. It slightly increases after steaming and CO treatment (red columns). In contrast to the dry feed, the ratio distinctly increases when the wet feed is employed. In addition, steaming and even more CO treatment cause a further distinct increase of the ratio. Whether the distinct improvements in the wet feed are more related to the increase in propene formation rate or by reduced CO2 formation rate becomes apparent when the formation rate of each product, i.e., propene or CO2, after respective treatment are compared to the initially measured formation rate of each product before pretreatments (Figure 12B). The loss in activity after treatments (Figure 1) is clearly mirrored in the loss of the CO2 formation rate. The corresponding ratio of the initial CO2 formation rate after each treatment and the initially measured CO2 formation rate (Figure 12, “initial”) decreases accordingly. The decrease of the CO2 formation rate appears significantly after treatment in steam and even more in CO. In particular, the CO2 formation rate strongly decreases in the wet feed. Therefore, treating the catalyst with reducing atmospheres results mainly in the loss of formation rate of the total combustion product CO2. Presumably, the surface structure responsible for total oxidation is degraded by the treatments, because the feed was the same before and after treatments and thus their effects related to the reactants can be excluded. However, at higher conversions of propane, the performance of the catalyst approaches to those initially measured suggesting that the structures required for total oxidation are slowly recovered. It still cannot be deduced whether the CO2 which is less formed due to the treatment originates from direct oxidation of propane or oxidation of intermediately formed propene. However, the selectivity changes could indicate an important role of oxygen on perovskite catalysts in propane oxidation. Thus, Osurface species could be responsible for both the high activity and the high total oxidation yield. These oxygen species are present in all conditions applied. Although the selectivity to propene is improved by reductive treatments in wet feed, CO2 formation cannot be fully suppressed. However, the CO2 formation rate increases when switching from the wet feed to the dry feed, which is in agreement with the observed decreased selectivity to propene and thus to increased total oxidation. In this discussion, it is assumed that the changes in selectivity are governed by the changes in CO2 formation rate while the apparent formation rate of propene varies little (Figure 12B). The Osurface species might be more important than OOH and Odefect in terms of selectivity, which explains the lower selectivity to propene in dry feed, although more Odefect was measured here (Figure 8 and Figure 9). The Osurface species could also be required for the healing of the oxygen vacancies and thus reduction of Odefect species, which have been created by treatments at higher temperatures. This is shown in Figure 13. The selectivity to propene increases after CO treatment with increasing proportion of Odefect, but decreases again after switching to dry feed, because Odefect decreases and Osurface increases. The concentration of surface defects and the concentration of electrophilic oxygen on the surface are interrelated and also have opposite and simultaneous effects on the selectivity to propene, which makes direct correlations difficult.

FIGURE 12
www.frontiersin.org

FIGURE 12. (A) Ratio of initial propene formation rate and CO2 formation rate in dry feed (red) and wet feed (blue). (B) Ratio of initial propene (C3H6) formation rates (rinitial (C3H6)) after the different treatments and initial propene formation rate at the beginning of the experiment (condition: initial in dry feed (red) and in wet feed (blue) on the top. The ratio of initial CO2 formation rates after the different treatments and initial CO2 formation rate (rinitial (CO2)) at the beginning of the experiment (condition: initial in dry feed (light grey) and wet feed (dark grey) on the bottom. Data obtained from contact time variation at 300°C. Contact time (W/F) ranges from 0.00347 g·min·L−1 to 0.00947 g·min·L−1 in dry feed (C3H8/O2/inert = 5/10/85) and wet feed (C3H8/O2/inert/H2O = 5/10/45/40).

FIGURE 13
www.frontiersin.org

FIGURE 13. Selectivity to propene at 4% conversion of propane in wet feed (5/10/45/40), in wet feed after CO treatment and in dry feed (5/10/85/0) in the top panel. Corresponding fractions of Odefect and Osurface are plotted in the bottom panel.

In summary, it can be concluded that the catalyst is stable in its bulk structure but the surface changes significantly under propane oxidation conditions as a function of the feed applied and thus as a function of the chemical potential. The local geometry of the terminating atoms changes and is limited to the outermost surface of the catalyst. These changes are clearly reflected in the core level and valence band spectra, as well as in the X-ray absorption spectra. However, in total these structural surface transformations are of minor extent. Therefore, the positive effect on the selectivity in propane oxidation is only slight and only significant at low conversion rates. But these results clearly show that for the optimization of a catalyst, not only its solid state and surface chemistry must be known, but also the kinetics of the chemical reactions in the solid state under different process conditions. In the search for more selective catalysts, other methods of catalyst synthesis should be used that more effectively increase the number of defects and thus the dynamic transformation of the catalytically active surface structures.

Data Availability Statement

The raw data supporting the conclusion of this article will be made available by the authors, without undue reservation.

Author Contributions

The manuscript was written by GK in direct iterative discussion with MH and ATr. GK prepared the catalysts and characterized them using standard methods and the assistance of colleagues named in the acknowledgements. PK carried out and evaluated the catalyst tests. The measurements at the synchrotron were carried out by GK and MH with the support of colleagues mentioned in the acknowledgements. The measurements were evaluated by MH and GK. ATa carried out the thermoanalytical investigation. RS contributed to the general discussion and interpretation.

Funding

This work was conducted in the framework of the BasCat collaboration among BASF SE, Technical University Berlin, Fritz-Haber-Institut der Max-Planck-Gesellschaft, and the cluster of excellence “Unified Concepts in Catalysis” (UniCat http://www.unicat.tu-berlin.de).

Conflict of Interest

The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.

Publisher’s Note

All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.

Acknowledgments

The Helmholtz-Zentrum Berlin is acknowledged for providing beamtime (proposal number 191-08014) and for the continuous support of the ambient-pressure activities of the MPG at BESSY II. Spencer J. Carey, Yuanqing Wang, Giulia Bellini and Detre Teschner are acknowledged for assistance in collecting spectra at BESSY. Jasmin Allan and Frank Girgsdies are acknowledged for TG and XRD measurements, respectively. Maike Hashagen is acknowledged for N2 physisorption measurement. Olaf Timpe and Gudrun Auffermann are acknowledged for elemental analysis. Frank Rosowski, BASF and BasCat, is kindly acknowledged for discussions during the initial phase of the experiments.

Supplementary Material

The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fchem.2021.746229/full#supplementary-material

References

Banger, K. K., Yamashita, Y., Mori, K., Peterson, R. L., Leedham, T., Rickard, J., et al. (2011). Low-temperature, High-Performance Solution-Processed Metal Oxide Thin-Film Transistors Formed by a ‘sol-Gel on Chip' Process. Nat. Mater. 10, 45–50. doi:10.1038/nmat2914

PubMed Abstract | CrossRef Full Text | Google Scholar

Bartel, C. J., Sutton, C., Goldsmith, B. R., Ouyang, R., Musgrave, C. B., Ghiringhelli, L. M., et al. (2019). New Tolerance Factor to Predict the Stability of Perovskite Oxides and Halides. Sci. Adv. 5, eaav0693. doi:10.1126/sciadv.aav0693

PubMed Abstract | CrossRef Full Text | Google Scholar

Brunckova, H., Kanuchova, M., Kolev, H., Mudra, E., and Medvecky, L. (2019). XPS Characterization of SmNbO4 and SmTaO4 Precursors Prepared by Sol-Gel Method. Appl. Surf. Sci. 473, 1–5. doi:10.1016/j.apsusc.2018.12.143

CrossRef Full Text | Google Scholar

Ciambelli, P., Cimino, S., De Rossi, S., Faticanti, M., Lisi, L., Minelli, G., et al. (2000). AMnO3 (A=La, Nd, Sm) and Sm1−xSrxMnO3 Perovskites as Combustion Catalysts: Structural, Redox and Catalytic Properties. Appl. Catal. B: Environ. 24, 243–253. doi:10.1016/s0926-3373(99)00110-1

CrossRef Full Text | Google Scholar

Cooper, P. W., and Kurowski, S. R. (1996). Technology of Explosives. New York, NY: Wiley.

Evarestov, R. A., Kotomin, E. A., Fuks, D., Felsteiner, J., and Maier, J. (2004). Ab Initio calculations of the LaMnO3 Surface Properties. Appl. Surf. Sci. 238, 457–463. doi:10.1016/j.apsusc.2004.05.171

CrossRef Full Text | Google Scholar

Evarestov, R. A., Kotomin, E. A., Heifets, E., Maier, J., and Borstel, G. (2003). Ab Initio Hartree-Fock Calculations of LaMnO3 (110) Surfaces. Solid State. Commun. 127, 367–371. doi:10.1016/s0038-1098(03)00440-x

CrossRef Full Text | Google Scholar

Evarestov, R. A., Kotomin, E. A., Mastrikov, Y. A., Gryaznov, D., Heifets, E., and Maier, J. (2005). Comparative Density-Functional LCAO and Plane-Wave Calculations of LaMnO3 Surfaces. Phys. Rev. B. 72, 214411. doi:10.1103/physrevb.72.214411

CrossRef Full Text | Google Scholar

Fan, J. C. C., and Goodenough, J. B. (1977). X‐ray Photoemission Spectroscopy Studies of Sn‐doped Indium‐oxide Films. J. Appl. Phys. 48, 3524–3531. doi:10.1063/1.324149

CrossRef Full Text | Google Scholar

Fierro, J. L. G., and Tejuca, L. G. (1987). Non-stoichiometric Surface Behaviour of LaMO3 Oxides as Evidenced by XPS. Appl. Surf. Sci. 27, 453–457. doi:10.1016/0169-4332(87)90154-1

CrossRef Full Text | Google Scholar

Frati, F., Hunault, M. O. J. Y., and De Groot, F. M. F. (2020). Oxygen K-Edge X-ray Absorption Spectra. Chem. Rev. 120, 4056–4110. doi:10.1021/acs.chemrev.9b00439

PubMed Abstract | CrossRef Full Text | Google Scholar

Galakhov, V. R., Demeter, M., Bartkowski, S., Neumann, M., Ovechkina, N. A., Kurmaev, E. Z., et al. (2002). Mn3s Exchange Splitting in Mixed-Valence Manganites. Phys. Rev. B. 65. doi:10.1103/physrevb.65.113102

CrossRef Full Text | Google Scholar

Gil, A., Gandía, L. M., and Korili, S. A. (2004). Effect of the Temperature of Calcination on the Catalytic Performance of Manganese- and Samarium-Manganese-Based Oxides in the Complete Oxidation of Acetone. Appl. Catal. A: Gen. 274, 229–235. doi:10.1016/j.apcata.2004.07.004

CrossRef Full Text | Google Scholar

Goldschmidt, V. M. (1926). Die Gesetze der Krystallochemie. Naturwissenschaften 14, 477–485. doi:10.1007/bf01507527

CrossRef Full Text | Google Scholar

Golikov, Y. V., Balakirev, V. F., and Fedorova, O. M. (2005). Homogeneity Range of Samarium Manganite in Air. Inorg. Mater. 41, 1010–1012. doi:10.1007/s10789-005-0253-3

CrossRef Full Text | Google Scholar

Grabowska, E. (2016). Selected Perovskite Oxides: Characterization, Preparation and Photocatalytic Properties-A Review. Appl. Catal. B: Environ. 186, 97–126. doi:10.1016/j.apcatb.2015.12.035

CrossRef Full Text | Google Scholar

Gupta, S., Kellogg, W., Xu, H., Liu, X., Cho, J., and Wu, G. (2016). Bifunctional Perovskite Oxide Catalysts for Oxygen Reduction and Evolution in Alkaline Media. Chem. Asian J. 11, 10–21. doi:10.1002/asia.201500640

PubMed Abstract | CrossRef Full Text | Google Scholar

Hui, L., Futai, M., and Yonghua, C. (1990). Study of the Preparation Process of Perovskite Oxides SmMnO3 by XRD and TPR Technique. React. Kinet Catal. Lett. 42, 151–156. doi:10.1007/bf02137632

CrossRef Full Text | Google Scholar

Hui, L., Futai, M., Zuzhu, M., Haiyang, L., and Yabo, X. (1993). Study of the Reactivity of Surface Oxygen Species on SmMnO3 and SmFeO3 Catalysts. J. Alloys Comp. 193, 68–69. doi:10.1016/0925-8388(93)90312-b

CrossRef Full Text | Google Scholar

Hwang, C.-C., Tsai, J.-S., Huang, T.-H., Peng, C.-H., and Chen, S.-Y. (2005). Combustion Synthesis of Ni-Zn Ferrite Powder-Influence of Oxygen Balance Value. J. Solid State. Chem. 178, 382–389. doi:10.1016/j.jssc.2004.10.045

CrossRef Full Text | Google Scholar

Ignatans, R., Mallia, G., Ahmad, E. A., Spillane, L., Stoerzinger, K. A., Shao-Horn, Y., et al. (2019). The Effect of Surface Reconstruction on the Oxygen Reduction Reaction Properties of LaMnO3. J. Phys. Chem. C. 123, 11621–11627. doi:10.1021/acs.jpcc.9b00458

CrossRef Full Text | Google Scholar

Kamata, K., Nakajima, T., and Nakamura, T. (1979). Thermogravimetric Study of Rare Earth Manganites AMnO3 (A=Sm,Dy,Y,Er,Yb) at 1200°C. Mater. Res. Bull. 14, 1007–1012. doi:10.1016/0025-5408(79)90065-5

CrossRef Full Text | Google Scholar

Kamata, K. (2019). Perovskite Oxide Catalysts for Liquid-phase Organic Reactions. Bcsj 92, 133–151. doi:10.1246/bcsj.20180260

CrossRef Full Text | Google Scholar

Knop‐Gericke, A., Kleimenov, E., Hävecker, M., Blume, R., Teschner, D., Zafeiratos, S., et al. (2009). “Chapter 4 X‐Ray Photoelectron Spectroscopy for Investigation of Heterogeneous Catalytic Processes,” in Advances in Catalysis (Academic Press), 213–272. doi:10.1016/S0360-0564(08)00004-7

CrossRef Full Text | Google Scholar

Koch, G., Hävecker, M., Teschner, D., Carey, S. J., Wang, Y., Kube, P., et al. (2020). Surface Conditions that Constrain Alkane Oxidation on Perovskites. ACS Catal. 10, 7007–7020. doi:10.1021/acscatal.0c01289

CrossRef Full Text | Google Scholar

Krill, G., Amamou, A., and Senateur, J. P. (1980). Valence Changes of Samarium Ions in Mixed SmS1-xPxcompounds Studied by Photoemission (XPS and UPS). J. Phys. F: Met. Phys. 10, 1889–1897. doi:10.1088/0305-4608/10/8/024

CrossRef Full Text | Google Scholar

Kurbakov, A. I. (2010). Electronic, Structural and Magnetic Phase Diagram of Sm1-xSrxMnO3 Manganites. J. Magnetism Magn. Mater. 322, 967–972. doi:10.1016/j.jmmm.2009.11.034

CrossRef Full Text | Google Scholar

Levasseur, B., and Kaliaguine, S. (2008). Effect of the Rare Earth in the Perovskite-type Mixed Oxides AMnO3 (A=Y, La, Pr, Sm, Dy) as Catalysts in Methanol Oxidation. J. Solid State. Chem. 181, 2953–2963. doi:10.1016/j.jssc.2008.07.029

CrossRef Full Text | Google Scholar

Li, X., Teschner, D., Streibel, V., Lunkenbein, T., Masliuk, L., Fu, T., et al. (2019). How to Control Selectivity in Alkane Oxidation?. Chem. Sci. 10, 2429–2443. doi:10.1039/c8sc04641g

PubMed Abstract | CrossRef Full Text | Google Scholar

Liu, L., Jia, J., Sun, T., and Zhang, H. (2018a). A Facile Method for Scalable Preparation of Mesoporous Structured SmMnO3 Perovskites Sheets for Efficient Catalytic Oxidation of Toluene. Mater. Lett. 212, 107–110. doi:10.1016/j.matlet.2017.10.048

CrossRef Full Text | Google Scholar

Liu, L., Sun, J., Ding, J., Zhang, Y., Jia, J., and Sun, T. (2019). Catalytic Oxidation of VOCs over SmMnO3 Perovskites: Catalyst Synthesis, Change Mechanism of Active Species, and Degradation Path of Toluene. Inorg. Chem. 58, 14275–14283. doi:10.1021/acs.inorgchem.9b02518

PubMed Abstract | CrossRef Full Text | Google Scholar

Liu, L., Zhang, H., Jia, J., Sun, T., and Sun, M. (2018b). Direct Molten Polymerization Synthesis of Highly Active Samarium Manganese Perovskites with Different Morphologies for VOC Removal. Inorg. Chem. 57, 8451–8457. doi:10.1021/acs.inorgchem.8b01125

PubMed Abstract | CrossRef Full Text | Google Scholar

Mcfarland, E. W., and Metiu, H. (2013). Catalysis by Doped Oxides. Chem. Rev. 113, 4391–4427. doi:10.1021/cr300418s

PubMed Abstract | CrossRef Full Text | Google Scholar

Merino, N., Barbero, B., Grange, P., and Cadús, L. (2005). LaCaCoO Perovskite-type Oxides: Preparation, Characterisation, Stability, and Catalytic Potentiality for the Total Oxidation of Propane. J. Catal. 231, 232–244. doi:10.1016/j.jcat.2005.01.003

CrossRef Full Text | Google Scholar

Merkle, R., and Maier, J. (2008). How Is Oxygen Incorporated into Oxides? A Comprehensive Kinetic Study of a Simple Solid-State Reaction with SrTiO3as a Model Material. Angew. Chem. Int. Ed. 47, 3874–3894. doi:10.1002/anie.200700987

PubMed Abstract | CrossRef Full Text | Google Scholar

Mierwaldt, D., Mildner, S., Arrigo, R., Knop-Gericke, A., Franke, E., Blumenstein, A., et al. (2014). In Situ XANES/XPS Investigation of Doped Manganese Perovskite Catalysts. Catalysts 4, 129–145. doi:10.3390/catal4020129

CrossRef Full Text | Google Scholar

Miyoshi, S., Hong, J.-O., Yashiro, K., Kaimai, A., Nigara, Y., Kawamura, K., et al. (2003). Lattice Expansion upon Reduction of Perovskite-type LaMnO3 with Oxygen-Deficit Nonstoichiometry. Solid State Ionics 161, 209–217. doi:10.1016/s0167-2738(03)00281-9

CrossRef Full Text | Google Scholar

Mori, T., Kamegashira, N., Aoki, K., Shishido, T., and Fukuda, T. (2002). Crystal Growth and crystal Structures of the LnMnO3 Perovskites: Ln=Nd, Sm, Eu and Gd. Mater. Lett. 54, 238–243. doi:10.1016/s0167-577x(01)00569-9

CrossRef Full Text | Google Scholar

Najjar, H., and Batis, H. (2016). Development of Mn-Based Perovskite Materials: Chemical Structure and Applications. Catal. Rev. 58, 371–438. doi:10.1080/01614940.2016.1198203

CrossRef Full Text | Google Scholar

Neagu, D., Tsekouras, G., Miller, D. N., Ménard, H., and Irvine, J. T. S. (2013). In Situ growth of Nanoparticles through Control of Non-stoichiometry. Nat. Chem. 5, 916–923. doi:10.1038/nchem.1773

PubMed Abstract | CrossRef Full Text | Google Scholar

Nesbitt, H. W., and Banerjee, D. (1998). Interpretation of XPS Mn(2p) Spectra of Mn Oxyhydroxides and Constraints on the Mechanism of MnO2precipitation. Am. Mineral. 83, 305–315. doi:10.2138/am-1998-3-414

CrossRef Full Text | Google Scholar

Pechini, P. M. (1967). Method of Preparing lead and Alkaline Earth Titanates and Niobates and Coating Method Using the Same to Form a Capacitor. US Patent No. 3330697.

Polo-Garzon, F., and Wu, Z. (2018). Acid-base Catalysis over Perovskites: a Review. J. Mater. Chem. A. 6, 2877–2894. doi:10.1039/c7ta10591f

CrossRef Full Text | Google Scholar

Powell, C. J., and Jablonski, A. (2010). NIST Electron Inelastic-Mean-Free-Path Database. Gaithersburg, MD: National Institute of Standards and Technology, U.S. Department of Commerce.

Royer, S., Duprez, D., Can, F., Courtois, X., Batiot-Dupeyrat, C., Laassiri, S., et al. (2014). Perovskites as Substitutes of noble Metals for Heterogeneous Catalysis: Dream or Reality. Chem. Rev. 114, 10292–10368. doi:10.1021/cr500032a

PubMed Abstract | CrossRef Full Text | Google Scholar

Salmeron, M., and Schlögl, R. (2008). Ambient Pressure Photoelectron Spectroscopy: A New Tool for Surface Science and Nanotechnology. Surf. Sci. Rep. 63, 169–199. doi:10.1016/j.surfrep.2008.01.001

CrossRef Full Text | Google Scholar

Stoerzinger, K. A., Hong, W. T., Azimi, G., Giordano, L., Lee, Y.-L., Crumlin, E. J., et al. (2015). Reactivity of Perovskites with Water: Role of Hydroxylation in Wetting and Implications for Oxygen Electrocatalysis. J. Phys. Chem. C. 119, 18504–18512. doi:10.1021/acs.jpcc.5b06621

CrossRef Full Text | Google Scholar

Stoerzinger, K. A., Hong, W. T., Crumlin, E. J., Bluhm, H., Biegalski, M. D., and Shao-Horn, Y. (2014). Water Reactivity on the LaCoO3 (001) Surface: An Ambient Pressure X-ray Photoelectron Spectroscopy Study. J. Phys. Chem. C. 118, 19733–19741. doi:10.1021/jp502970r

CrossRef Full Text | Google Scholar

Sunding, M. F., Hadidi, K., Diplas, S., Løvvik, O. M., Norby, T. E., and Gunnæs, A. E. (2011). XPS Characterisation of In Situ Treated Lanthanum Oxide and Hydroxide Using Tailored Charge Referencing and Peak Fitting Procedures. J. Electron Spectrosc. Relat. Phenomena 184, 399–409. doi:10.1016/j.elspec.2011.04.002

CrossRef Full Text | Google Scholar

Suntivich, J., Gasteiger, H. A., Yabuuchi, N., Nakanishi, H., Goodenough, J. B., and Shao-Horn, Y. (2011). Design Principles for Oxygen-Reduction Activity on Perovskite Oxide Catalysts for Fuel Cells and Metal-Air Batteries. Nat. Chem. 3, 546–550. doi:10.1038/nchem.1069

PubMed Abstract | CrossRef Full Text | Google Scholar

Tanuma, S., Powell, C. J., and Penn, D. R. (2003). Calculation of Electron Inelastic Mean Free Paths (IMFPs) VII. Reliability of the TPP-2M IMFP Predictive Equation. Surf. Interf. Anal. 35, 268–275. doi:10.1002/sia.1526

CrossRef Full Text | Google Scholar

Tanuma, S., Powell, C. J., and Penn, D. R. (1991). Calculations of Electorn Inelastic Mean Free Paths. II. Data for 27 Elements over the 50-2000 eV Range. Surf. Interf. Anal. 17, 911–926. doi:10.1002/sia.740171304

CrossRef Full Text | Google Scholar

Tasleem, S., and Tahir, M. (2020). Recent Progress in Structural Development and Band Engineering of Perovskites Materials for Photocatalytic Solar Hydrogen Production: A Review. Int. J. Hydrogen Energ. 45, 19078–19111. doi:10.1016/j.ijhydene.2020.05.090

CrossRef Full Text | Google Scholar

Toniolo, F. S., and Schmal, M. (2016). “Improvement of Catalytic Performance of Perovskites by Partial Substitution of Cations and Supporting on High Surface Area Materials,” in Perovskite Materials - Synthesis, Characterisation, Properties, and Applications. Editors L. Pan, and G. Zhu (IntechOpen). doi:10.5772/61279

CrossRef Full Text | Google Scholar

Wernbacher, A. M., Kube, P., Hävecker, M., Schlögl, R., and Trunschke, A. (2019). Electronic and Dielectric Properties of MoV-Oxide (M1 Phase) under Alkane Oxidation Conditions. J. Phys. Chem. C. 123, 13269–13282. doi:10.1021/acs.jpcc.9b01273

CrossRef Full Text | Google Scholar

Wołcyrz, M., Horyń, R., Bourée, F., and Bukowska, E. (2003). Structural Defects in LaMnO3 Phase Studied by Neutron Diffraction. J. Alloys Comp. 353, 170–174. doi:10.1016/S0925-8388(02)01313-0

CrossRef Full Text | Google Scholar

Yang, J. B., Kim, M. S., Creel, T. F., Zhao, H., Chen, X. G., Yelon, W. B., et al. (2016). “Structural, Magnetic and Transport Properties of B-Site Substituted Perovskite La0.7Sr0.3MnO3,” in Perovskite Materials - Synthesis, Characterisation, Properties, and Applications. Editors L. Pan, and G. Zhu (IntechOpen). doi:10.5772/61770

CrossRef Full Text | Google Scholar

Yeh, J. J. (1993). Atomic Calculation of Photoionization Cross-Sections and Asymmetry Parameters. Langhorn PE (USA): Gordon & Breach Science Publishers.

Yeh, J. J., and Lindau, I. (1985). Atomic Subshell Photoionization Cross Sections and Asymmetry Parameters: 1 ⩽ Z ⩽ 103. At. Data Nucl. Data Tables. 32, 1–155. doi:10.1016/0092-640x(85)90016-6

CrossRef Full Text | Google Scholar

Zhu, H., Zhang, P., and Dai, S. (2015). Recent Advances of Lanthanum-Based Perovskite Oxides for Catalysis. ACS Catal. 5, 6370–6385. doi:10.1021/acscatal.5b01667

CrossRef Full Text | Google Scholar

Zhu, Y., Liu, X., Jin, S., Chen, H., Lee, W., Liu, M., et al. (2019). Anionic Defect Engineering of Transition Metal Oxides for Oxygen Reduction and Evolution Reactions. J. Mater. Chem. A. 7, 5875–5897. doi:10.1039/c8ta12477a

CrossRef Full Text | Google Scholar

Keywords: SmMnO3, perovskite, surface modification, defects, heterogeneous catalysis, oxidative dehydrogenation of propane, AP-XPS, AP-NEXAFS

Citation: Koch G, Hävecker M, Kube P, Tarasov A, Schlögl R and Trunschke A (2021) The Influence of the Chemical Potential on Defects and Function of Perovskites in Catalysis. Front. Chem. 9:746229. doi: 10.3389/fchem.2021.746229

Received: 23 July 2021; Accepted: 06 September 2021;
Published: 17 September 2021.

Edited by:

Guowei Li, Max Planck Institute for Chemical Physics of Solids, Germany

Reviewed by:

Yingfen Wei, École Polytechnique Fédérale de Lausanne, Switzerland
Ravindra Pandey, Indian Institute of Technology Roorkee, India

Copyright © 2021 Koch, Hävecker, Kube, Tarasov, Schlögl and Trunschke. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.

*Correspondence: Annette Trunschke, trunschke@fhi-berlin.mpg.de

Download