
==== Front
J Am Chem Soc
J Am Chem Soc
ja
jacsat
Journal of the American Chemical Society
0002-7863
1520-5126
American Chemical Society

39178385
10.1021/jacs.4c06587
Article
N2O Decomposition on Singly and Doubly (K and Li)-Doped Co3O4 Nanocubes—Establishing Key Factors Governing Redox Behavior of Catalysts
https://orcid.org/0000-0002-0617-7355
Nowakowski Leszek †∥
Hudy Camillo †
https://orcid.org/0000-0003-2815-0644
Zasada Filip †
https://orcid.org/0000-0002-9777-0268
Gryboś Joanna †
https://orcid.org/0000-0002-3462-436X
Piskorz Witold †
https://orcid.org/0000-0003-3112-2759
Wach Anna ‡
Kayser Yves §⊥
https://orcid.org/0000-0002-7526-3097
Szlachetko Jakub ‡
https://orcid.org/0000-0001-7226-7626
Sojka Zbigniew *†
† Faculty of Chemistry Jagiellonian University, ul. Gronostajowa 2, Krakow 30-387, Poland
‡ National Synchrotron Radiation Centre SOLARIS Jagiellonian University, ul. Czerwone Maki 98, Kraków 30-392, Poland
§ Physikalisch-Technische Bundesanstalt (PTB), Abbestr. 2-12, Berlin 10587, Germany
∥ Doctoral School of Exact and Natural Sciences, Jagiellonian University, Prof. St. Łojasiewicza St 11, Krakow 30-348, Poland
* Email: sojka@chemia.uj.edu.pl.
23 08 2024
04 09 2024
146 35 2445024466
14 05 2024
09 08 2024
19 07 2024
© 2024 The Authors. Published by American Chemical Society
2024
The Authors
https://creativecommons.org/licenses/by/4.0/ Permits the broadest form of re-use including for commercial purposes, provided that author attribution and integrity are maintained (https://creativecommons.org/licenses/by/4.0/).

The intimate mechanism of N2O decomposition on bare and redox-tuned Co3O4 nanocubes (achieved by single (Li or K) and double (Li and K) doping) was elucidated. The catalysts synthesized by the hydrothermal method were characterized by X-ray electron absorption fine structure measurements, X-ray diffraction, Raman spectroscopy, scanning electron microscopy, transmission electron microscopy, and Kelvin Probe techniques. TPSR and steady-state isothermal catalytic tests reveal that the N2O turnover frequencies are critically sensitive to the work function of the catalysts, adjusted purposely by doping. For the catalysts obtained by one-pot hydrothermal synthesis, lithiation of the Co3O4 nanocubes leads to the formation of {Li’8a, Co·16d} species, decreasing steadily the work function and the activity, while for the catalysts prepared by postsynthesis impregnation, formation of {Li’8a, Co’16d, Co··16c} species leads to a volcano-type dependence of the catalytic activity and the work function in parallel. The beneficial effect of potassium was discussed in terms of mitigation of surface potential buildup due to the accumulation of ionosorbed oxygen intermediates (surface electrostatics), which hinders the interfacial electron transfer. Analysis of the catalytic activity response to the redox tuning of Co3O4, substantiated by DFT calculations, allowed for a straightforward conceptualization of the redox nature of the N2O decomposition in terms of the lineup of frontier orbitals of the N2O/N2O– and O2–/O2 reactants with the surface DOS structure and the resultant molecular orbital interactions. The positions of the virtual bonding 3πg0(N2O)−α-3dz2 and the occupied 2πg1(O2–)−α-3dz2 states relative to the Fermi energy level play a crucial role in the regulation of the forward and backward interfacial electron transfer events, which drive the redox process.

Narodowe Centrum Nauki 10.13039/501100004281 2017/27/B/ST4/01155 Polish Ministry of Higher Education NA 1/SOL/2021/2 Infrastruktura PL-Grid 10.13039/501100011089 465001049 Narodowe Centrum Nauki 10.13039/501100004281 2020/37/B/ST3/00555 document-id-old-9ja4c06587
document-id-new-14ja4c06587
ccc-price
==== Body
pmc1 Introduction

Catalytic decomposition of N2O has been widely studied over simple and mixed oxides,1−5 ABO3 perovskites,6,7 AB2O4 spinels,8−10 hydrotalcite-derived oxides,11,12 mesoporous silica materials,13 zeolites,14−16 and a variety of supported catalysts.17−20 Among them, cobalt spinel catalysts belong to one of the most promising systems for fundamental and applied investigations.21−23 Depending on the synthesis conditions, cobalt spinel (Co3O4) can be easily obtained in the form of well-shaped cubic, octahedral, or cuboctahedral nanocrystals,24 and its composition modified via substitution at A and B sites, providing excellent materials for rigorous fundamental investigations into the structure-redox reactivity relationships, using N2O decomposition as a convenient model process.1,2,25,26

Substantial improvement in the redox performance of the catalyst may be achieved by bulk structure modification (doping with alien cations)1,2,27−29 or by tuning the electronic properties with alkali promoters.26,30−33 The latter effect is usually significant and results in a lowering of the reaction temperature of N2O decomposition even by 200 °C.1,2,25 In our previous articles, we have shown that such beneficial action is mainly of electronic origin and is associated with the presence of alkali adspecies (K, Cs)25,26,33 but not with their diffusion into the catalyst bulk as previously suggested.27 The alkali promoters present on the catalyst surface, by lowering the work function of the cobalt spinel, facilitate the forth and back interfacial electron transfer processes that occur between the spinel active sites and the N2O reactant (dissociative reduction step) and the reaction of the O2–(ads)/O22–(ads) intermediates (oxidative oxygen evolution step). Although the activation energy of the decomposition of nitrous oxide over potassium or cesium-promoted cobalt spinel is fairly well correlated with changes in the work function of the catalyst,2,25,26,34 the intimate molecular nature of this process and its relationship with the redox properties of the catalyst have not been elucidated in the required detail so far.

Incorporation of the nonredox dopant ions into the spinel matrix provides another way of tuning the electronic properties of the catalyst on purpose.2,35 The framework positions 8a (tetrahedral) and 16c (octahedral) or interstitial 16d (octahedral) and 48f (tetrahedral) are the possible loci for such substitution (for details, see Electronic Supporting Information, Figure S1, and the associated text). This has allowed for straightforward identification of the cobalt active site due to the selective removal of Co2+ in the 8a positions by doping with redox inert Mg2+ cations or replacement of Co3+ in the 16c positions with Al3+.35 Valence pinning of cobalt cations can be achieved in turn by aliovalent substitution. In this context, lithium ions, due to their small radius, are particularly suitable species to be incorporated.36−38 The presence of hypovalent Li cations in the spinel matrix gives rise to an alternation of the cobalt oxidation state imposed by a charge compensation constraint, modifying the catalyst Fermi level energy (EF) and the associated catalytic properties.2,39−42

The inherent intricacy of the interfacial redox processes, resulting from the complex DOS structure of the semiconducting oxide catalyst, and a variety of alignment patterns between the energy bands of the active sites and frontier (HOMO–LUMO) orbitals of the reactants, makes a comprehensive, in-depth understanding of all involved molecular events still a challenging endeavor. As a functional approach, several redox reactivity descriptors have been proposed (particularly in accounting for the ORR and OER electrocatalytic activity), including eg level filling, antibonding orbital occupancy, the density of the states near the Fermi level, and the d-band center or position of the highest occupied d-state vs the Fermi level.43−46 In the field of thermal oxide catalysis, the establishment of analogous descriptors is still less advanced.47−49 However, by realizing the complex nature of the interfacial redox processes, the validity of such singular parameters focused on only the catalyst component is of rather limited value for establishing a more general account. In this context, the relative energies of the electronic levels of the reactants/intermediates and the catalyst active sites along with the corresponding molecular orbital coupling strengths play a pivotal role.

A concise epitomization of the heterogeneous redox process in terms of the electronic structure of both interacting moieties (catalyst surface DOS structure and the frontier molecular orbital (FMO) pattern of the reactant admolecules), taking into account the corresponding orbital interactions (p-d hybridization), the effects resulting from the position of the Fermi level, and variable surface electrostatics in particular, has not been achieved as yet.

The catalytic gas/solid redox processes are inherently associated with interfacial forth electron transfer (reduction of reactants/intermediates) and back electron transfer (oxidation of reactants/intermediates) events. Due to intimate contact between the reaction pair, allowing for sizable orbital interactions, they usually proceed along the inner-sphere, adiabatic routes featured by strong orbital coupling.50 In such electron transfer processes, bonds are often formed and broken, and semiclassical approaches based, e.g., on the candid Gerischer model,51 cannot be used directly for interpretation of the redox behavior of the involved catalyst/reactant couples.

Despite the apparent mechanistic simplicity of the N2O decomposition over oxide catalysts,1,2,52 the kinetics of this reaction can often be complex, or even elusive, due to a varying number of the active sites with temperature, the presence of the charged Onz–(ads) intermediates, and mixed reaction-surface diffusion kinetics.53−56 Therefore, it may exhibit unexpectedly intriguing characteristics regarding, e.g., a wide variation of the activation energies values (ΔEa), reported even when nominally the same catalysts (such as cobalt spinel) are used.1,2,57−59

In the present study, we address these points for catalytic N2O decomposition, which served as a useful redox probe reaction, choosing an euhedral cobalt spinel and its doped derivatives as a dedicated model catalytic system. The aim is to evaluate the redox behavior of cobalt spinel catalysts of a well-defined nanocubic shape and unravel the effects resulting from bulk and surface doping of Co3O4 with two different alkali cations (lithium and potassium), introduced alone or jointly. The resulting controlled electronic structure perturbations of the cobalt spinel allow for straightforward clarification of the observed changes in the catalytic activity in terms of the Fermi level and surface potential variation (due to the accumulation of the Onz–(ads) intermediates on the catalyst surface), playing along with the substantial interfacial orbital interactions a central function.

The constructed conceptual model of N2O decomposition on the cobalt spinel p-type semiconductor catalysts rationalizes the mechanism of the interfacial back and forth electron transfer events, which govern the course of this redox reaction. Due to the focus on the redox nature of the deN2O reaction only, the effect of typical contaminants (O2, H2O, and NO) was intentionally omitted, as it has been thoroughly elucidated for cobalt spinel-based catalysts in our previous paper.60

2 Experimental Section

2.1 Catalysts Synthesis

Cube-shaped Co3O4 nanoparticles were synthesized by the template-free hydrothermal method, described previously by us.24,61 For this purpose, Co(NO3)2·6H2O and NaOH precursors with the molar ratio 2:1 were dissolved in 10 mL of deionized H2O and transferred to a 20 mL Teflon-lined steel autoclave. The reaction mixture was heated at 180 °C for 5 h and then cooled down to room temperature. The obtained black precipitate was separated by centrifugation, washed several times with distilled water to remove sodium ions, and dried at 60 °C in air overnight (this sample is labeled as h-Co3O4). The final product was annealed at 500 °C in the air for 5 h to form monodispersed cobalt spinel nanocubes of high crystallinity and submicrometer size (c-Co3O4 samples). The lithium-doped h-LixCo3–xO4 series (0.040 ≤ x ≤ 0.099) was obtained via an analogous one-pot hydrothermal synthesis by adding NaOH along with the LiOH dopant in the 2:1, 1:1, and 1:2 molar ratios to the reaction mixture before the hydrothermal reaction. For the synthesis of the i-LixCo3O4 (0.033 ≤ x ≤ 0.060) and i-Ky/Li0.045Co3O4 (0.0013 ≤ y ≤ 0.012) series, the calcined c-Co3O4 nanocubes were next doped with lithium and potassium by incipient wetting impregnation using aqueous solutions of LiNO3 or KNO3, respectively. The impregnated spinels were finally calcined at 600 (for the Li doping step) and 450 °C (for the K doping step). Following the literature, after such treatment, the small Li cations are incorporated into the spinel matrix,62 whereas the larger potassium ions remain on the surface.2,25 An auxiliary series of the lithiated samples with an extended Li/Co content ranging from 3.29 × 10–3 to 2.31 × 10–1 was obtained by dry impregnation of Co3O4 (obtained by precipitation with NaOH33,35) using aqueous solutions of lithium nitrate of appropriate concentrations.

2.2 Methods

X-ray diffraction patterns were recorded with a Rigaku Miniflex X-ray diffractometer in the 2θ range of 15–85° with a resolution of 0.02° and a time interval of 1s per step. The chemical composition of lithium-doped samples was determined by inductively coupled plasma mass spectrometry (ICP-MS, PerkinElmer with a FIAS system). For this purpose, the spinel sample was dissolved in concentrated HCl (∼35%).

The loading of potassium deposited on the surface of the spinel catalysts was assessed by means of energy-dispersive X-ray spectrometry (XRF, Thermo Scientific, ARL QUANT’X), equipped with the Rh anode operating at an acceleration voltage of 4–50 kV (1 kV steps) and beam size of 1 mm.

The μ-Raman spectra were recorded in the range of 100–900 cm–1 with 1 cm–1 resolution using a Renishaw InVia spectrometer, equipped with the CCD detector and a confocal Leica DMLM microscope, with an excitation wavelength of 514 nm. Nine scans were accumulated for each measurement to achieve a satisfactory signal-to-noise ratio. Scanning electron microscopy (SEM) imaging of the spinel nanocrystals’ morphology was performed on a Tescan VEGA 3 equipped with a LaB6 cathode. The samples were gold-coated before the microscopic observations. Scanning transmission electron microscopy (STEM) imaging of the samples was performed with a FEI Tecnai/Osiris electron microscope equipped with an X-FEG Schottky emitter (200 kV). For energy dispersive X-ray (EDX) mapping, a 4-sector silicon drift windowless detector and the Bruker Esprit software were used. Prior to the microscopic analyses, the samples were deposited on a lacey carbon film supported on a copper grid (Agar Scientific, London, UK, 300 mesh).

Contact potential difference (CPD) measurements were carried out by the Kelvin dynamic condenser method with a KP6500 probe (McAllister Technical Services). To standardize the catalyst surface, the measurements were carried out in vacuum (10–7 mbar) at 150 °C after annealing the samples at 400 °C. The work function values, Φ, were calculated from the relation: CPD = Φreference – Φsample, using a standardized stainless steel plate as a reference (d = 3 mm, Φreference = 4.3 eV).

X-ray absorption fine structure (XAFS) measurements in the soft X-ray energy range were performed at the plane grating monochromator beamline of Physikalisch-Technische Bundesanstalt at the electron storage ring BESSY II. The XAFS spectra were collected in the fluorescence mode after installing the samples on a goniometer in an ultrahigh vacuum chamber equipped with a windowless silicon drift detector, which was oriented at 90° and in the polarization plane of the synchrotron radiation.63 The samples were oriented at 45° to the incidence direction and with the normal vector to the sample surface plane contained within the polarization plane of synchrotron radiation. The X-ray beam spot size on the sample was around 0.14 mm2. The Co L2,3-edge XAFS spectra were acquired in the 750–830 eV energy range with energy increments of 500 meV between data points. The OK-edge spectra were collected in the energy range of 500 to 600 eV with the same energy increments. The XAFS spectra were normalized to the 0 to 1 range for subsequent analysis. The X-ray absorption spectra in the hard X-ray energy range were collected on a laboratory spectrometer. The setup was composed of an X-ray source (XOS X-Beam Superflux PF with a Mo anode) and a von Hámos geometry-based X-ray spectrometer to detect the X-ray radiation transmitted through the sample. The spectrometer was adjusted to diffract radiation by the Si(110) crystal in the diffraction order of n = 4 and detect the signal on the two-dimensional X-ray camera (Andor Newton DO920P). The size of the X-ray beam spot on the sample was 100 μm × 100 μm. The energy calibration was performed numerically by finding the best fit of the reference spectrum (from the XAS database) to the measured spectrum of a cobalt foil, assuming a linear pixel-to-energy conversion function. The recorded spectra were normalized in the 0–1 range.

2.3 Catalytic Measurements

The isothermal and temperature-programmed (TPSR) catalytic N2O decomposition tests were performed in the 25–600 °C range using a homemade setup operating in the LabView environment and equipped with a Hiden Analytical HPR20 QMS detector, Brooks mass controllers, and automatic Valco switching valves. The experiments were carried out using a quasi-CSTR-type quartz flow reactor and 150 mg of the catalyst (sieve fraction of 0.2–0.3 mm). The feed flow of 42–102 mL·min–1 and the heating rate of 10 °C·min–1 were used. The N2O conversion was calculated based on the QMS signals normalized against the helium gas balance as1

whereas the reaction turnover frequencies, TOFN2O, were calculated using the equation2

where F = 30 mL/min is the flow rate, Vm = 22,400 cm3 is the gas molar volume (T = 298 K, p = 1 bar), mcat = 0.15 g is the catalyst mass, SBET is its specific surface area, and nsurf = 6 × 1018 1/m2 is the surface concentration of Co3+ cations on the (100) facet.64 The reproducibility of the catalytic tests was checked by performing 3 consecutive measurements of N2O decomposition on bare Co3O4, Li-, K-, and (Li,K)-doped cobalt spinel catalysts in the TPSR mode. The obtained N2O conversion profiles were virtually identical (see Figure S10 in the ESI Section). The accuracy of activation energy determination was assessed using standard error analysis of the slope in the Arrhenius equation. The apparent reaction order (m = ∂ ln rN2O/∂ ln pN2O) was determined from the pressure dependence of the reaction rate (rN2O = kpN2Om).

2.4 Molecular Modeling

DFT calculations were carried out using the VASP code and the projector-augmented plane wave method (PAW).65 The preliminary geometry optimization was performed at the PB/DFT + U level of theory with the Hubbard parameter U = 3.5 eV applied to the Co-3d states,66,67 whereas the electronic structure calculations were performed with the hybrid HSE06 functional. Typically, we used the standard Monkhorst–Pack68 grid with the 5 × 5 × 5 sampling mesh for the bulk and the 3 × 3 × 1 mesh for slab calculations; the cutoff energy was set to 500 eV, and the SCF convergence criterion was set to 1 × 10–5 eV. For the atomic position relaxation, the conjugate-gradient method, improved by Brent’s steps corrector algorithm, was employed,69 within the convergence criterion of 1 × 10–4 eV·Å–1. The population analysis was carried out by means of the Bader method. The bulk cobalt spinel unit cell was obtained by fully optimizing all internal degrees of freedom of the cubic (1 × 1 × 1) cell containing 56 ions (Co24O32). For the modeling of the localized electronic states, the occupation matrix control method (OCC-MAT) was used.70 The work function of the parent and doped catalysts was calculated using slab models of different sizes (1 × 1, 2 × 1, and 3 × 1 for modeling the dopant loading effect), with 11 Å of the oxide thickness (for details, see Figure S2 in the ESI Section). The atomic positions in the four top and four bottom layers of the slab models were relaxed within the criterion of 1 × 10–4 eV·Å–1. The applied vacuum spacing was sufficiently large for the successful convergence of the planar averaged electrostatic potential to a constant value far from the surface, taken as the reference vacuum level (Figure S3 in ESI). Projected crystal orbital Hamilton population (pCOHP) was analyzed using the Lobster v5.0 package.71 This method allows for partitioning the band-structure energy into bonding, nonbonding, and antibonding contributions. The integrated projected crystal orbital hamilton population (IpCOHP) calculated up to the Fermi level serves to quantify the strength of the chemical bonds between the involved atomic moieties. Following the Newns-Anderson type treatment, the lifetime of the activated adspecies was estimated from the halfwidth of the chemisorption (hybridization) function approximated by the Lorentzian shape as τ = (η/2)Δ–1.72

3 Results and Discussion

3.1 Characterization of the Catalyst

To determine the critical value of Li doping of Co3O4 that preserves the spinel phase homogeneity (lack of undesired, spurious LiCoO2), we examined at first an auxiliary series of the lithiated samples with the extended Li/Co content ranging from 3.29 × 10–3 to 2.31 × 10–1 (see Table S1 in ESI). As revealed by the X-ray diffraction (XRD) and Raman spectra (Figure S4), the spinel phase is preserved for Li/Co < 4 × 10–2. Above this threshold value, gradual formation of the surface LiCoO2 phase was observed, as revealed by the development of the (006) and (104) diffraction peaks (see Figure S4b1,b2, respectively). The formation of the LiCoO2 phase was also confirmed by the appearance of bands at 490 and 595 cm–1, assigned to the Eg and A1g vibrations, respectively (Figure S4c). As a result, for further investigations, we restricted the Li/Co ratio below the ascertained critical limit since the LiCoO2 phase appeared completely inactive in the N2O decomposition in comparison to the Co3O4 phase, as revealed in a reference experiment (see Figure S4d).

A list of the spinel pure phase samples doped by Li and K (the prime series), used for further investigation, and their compositions determined by ICP-MS (Li) and XRF (K) techniques are shown in Tables S2 and S3 in the ESI Section. In the case of potassium, noting its accommodation on the surface (see below), its areal atomic concentration, nK, was calculated based on the BET surface area value.

The phase purity of the parent cobalt spinel and the Li and K-doped catalysts was verified by XRD and Raman spectroscopic measurements (Figures 1, S5 and S6). The diffraction patterns (Figures 1a and S5a, S6a) show that the spinel catalysts obtained are highly crystalline (all the observed XRD lines can be indexed within the Fd3̅m space group), 24210-ICSD, and no segregated phases comprising Li or K were identified in the applied doping range. Their possible presence in high surface dispersions (preventing XRD detection) was also excluded by Raman measurements (Figures 1b, S5b, and S6b). The Raman spectra of the Co3O4, i-LixCo3O4, h-LixCo3–xO4, i-Ky/Co3O4, and i-Ky/Li0.045CoO4 samples are all characterized by 5 bands located at 194 (F2g), 480 (Eg), 520 (F2g), 620 (F2g), and 690 (A1g) cm–1,73 confirming good purity and high crystallinity of the singly and doubly doped spinel catalysts.

Figure 1 Selected XRD patterns (a) of the parent Co3O4 (a1—green line), i-Li-2-Co (a2—blue), h-Li-2-Co (a3—purple), and i-K-5/i-Li-2-Co samples (a4—yellow), and the corresponding Raman spectra (b1–b4) together with SEM pictures (c–f). Exemplary STEM/EDX map of the potassium distribution in the i-K-5/i-Li-2-Co sample (g). The results of the characterization of the remaining samples are shown in Figures S5 and S6 in ESI.

The selected SEM pictures (Figure 1c–f), representative for the given series of the catalyst, reveal that the synthesized spinel nanocrystals exhibit clear-cut cubic morphology (see reconstructed shapes in the inset in Figure 1c), with a size of 780 ± 60 nm (gauged by the body diagonal), regardless of the presence or absence of the dopants and the way along which they were introduced.

STEM/EDX mapping confirmed accumulation of potassium essentially on the surface of the spinel nanocubes for both singly K- and doubly (K,Li)-promoted catalysts, as expected for the postsynthesis impregnation applied for potassium deposition and the fact that the large ionic radius prevents its effective incorporation into the spinel matrix. Exemplary map of K-distribution for the most active i-K-5/i-Li-2-Co catalyst is shown in Figure 1g.

For lithium doping via one-pot hydrothermal synthesis of the catalysts, the Li cations are directly incorporated into the spinel matrix in the course of the reaction. Whatever is the framework location of the lithium ions (in the tetrahedral 8a or octahedral 16d sites), in this case, the parent cobalt cations always become oxidized to maintain the charge balance. Since disclosing the actual Li locus by the applied spectroscopic methods appeared elusive (due to the detection limitations), we established the preferred tetrahedral 8a site of lithium localization by means of DFT modeling, applying the occupational matrix (OCC-MAT) scheme (see below).

The impregnation of the Co3O4 nanocubes with the LiNO3 solution and subsequent calcination were monitored by using QMS and Raman techniques. In the temperature range of 300–450 °C, an evolvement of NO and O2 (with a small amount of NO2) was observed due to the decomposition of the nitrate anions (Figure S7). The incorporation of the Li+ cations into the spinel lattice was confirmed by a systematic hypsochromic shift of the Raman A1g band with an increasing lithium content (Figure S8). There are two possible scenarios of this process depending on the mechanistic subtleties of the adsorbed nitrate decomposition. When this step is concomitant with the substitution of Li+ for the tetrahedral Co2+ cations, the latter are shifted from the framework 8a into the empty interstitial octahedral 16c positions.74 This is accompanied by a reduction of the Co3+ cations located at the 16d sites and a commensurate evolution of NO and O2 in equal amounts3

as observed in Figure S7, indeed.

In an alternative mechanism of lithium incorporation, the nitrate precursor is first decomposed into Li2O adspecies, with an equal formation of the NO, O2, and NO2 gas products.4a

The resultant xLi2O entities are next integrated with the spinel host, and the Li+ ions entering the 8a sites oxidize the Co3+ cations in the 16d positions to maintain the charge balance.4b

In the first case, Li doping is tantamount to the formation of two Li’8a and Co’16d species, which are charge-balancing the interstitial Co··i centers. In the second case, the Li’8a defects are balanced by the hole Co·16d centers, as in the case of the one-pot synthesis. The evolution of small amounts of NO2 (Figure S7) speaks in favor of the Co3+ reduction described in eq 3. Nevertheless, to ascertain which of these two possibilities actually takes place, we performed X-ray absorption spectroscopy (XAFS) investigations into the alterations of the cobalt oxidation state upon the lithiation. It should be emphasized, however, that in both cases Li+ cations are located in the tetrahedral (8a) sites in the spinel matrix, but the cobalt valence pinning effects imposed by the disparate routes of their incorporation are different.

The redox and structural changes induced by Li doping were examined using the oxygen K, and cobalt L and K-edges for the bare Co3O4, and the most active (among the lithiated spinels) catalyst i-Li-2-Co (i-Li0.045Co3O4), see Figure 2a,b. The collected Co L-edge spectra show two main regions, lower energy L3-edge and higher energy L2-edge, which arise from transitions of 2p3/2 → 3d and 2p1/2 → 3d, respectively. The L2,3-edge spectra provide rich electronic structure information, being highly sensitive to valency, spin states, and coordination geometries.75 Generally, the recorded XAFS Co L-spectrum of the Li-impregnated spinel exhibits clear features (marked by orange arrows) indicative of a partial reduction of Co3+ to Co2+. This observation is further confirmed by the analysis presented in Figure 2b, where the differences between the O K-edge spectra of the CoO, Co3O4, and i-LixCo3O4 samples are directly compared. The peak position of the O K-edge at 531 eV, assigned to O 1s → O 2p transitions,76 remains essentially intact for the c-Co3O4, i-LixCo3O4, and CoO samples, as expected for the unchanged oxidation state of the O2– anions, but its intensity is significantly reduced upon Li doping. Together with the associated more dramatic changes in the postedge region, these variations speak in favor of the incorporation of lithium into the cobalt spinel according to eq 3 when introduced via impregnation. Inspection of the Co L2,3-edges properly confirmed this finding. The γ and δ features in the L3-edge, assigned to Co3+ in the octahedral low spin coordination, decrease in line with the fact that the content of the Co3+16d (S = 0) cations is reduced upon Li doping (feature α corresponds to high spin octahedral Co2+, characteristic of CoO),77 [and references therein]. The concomitant decline and broadening of the β features, due to high spin tetrahedral Co2+, are consistent with the presence of Li cations in the 8a sites, resolving the incorporation mechanism definitely.

Figure 2 Co L-edge (a), O K-edge (b), and Co K-edge XAFS spectra (c) together with the cobalt <AOS> correlation diagram (d) for i-Li0.045Co3O4 (a1–d1), c-Co3O4 (a2–d2), h-Co3O4 (c3,d3) h-Li0.086Co2.94O4 (c4,d4), and the reference CoO samples.

The cobalt K-edge spectra of Co3O4 and CoO show a significant energy shift between the valence states 2+ and 3+, allowing reliable elucidation of the cobalt average oxidation state in the Li-doped spinel catalysts, based on the position of the inflection point in the recorded K-edge spectra (Figure 2c). The resulting average oxidation state values of cobalt, <AOS>, for the bare h-Co3O4 (<AOS> = 2.76) and c-Co3O4 (<AOS> = 2.67) samples, and the lithiated i-Li0.045Co3O4 (<AOS> = 2.55) and h-Li0.086Co2.94O4 (<AOS> = 2.89) catalysts are shown in Figure 2d.

Enhancement of the <AOS> values of cobalt upon calcination of the h-Co3O4 catalyst obtained by the hydrothermal method results from the annihilation of the cationic vacancies (the missing surface Co2+ cations), which is accompanied by oxygen release, as previously shown by us.78 The calcination process, Co1–xCo2O4 →(1 – x)Co3O4 + 2xO2, leads to a decrease in the <AOS> value from 8/(3 – x) = 2.76 to 8/3, as observed in the XAFS spectra (Figure 2c,d). Thus, the composition of the hydrothermally obtained cobalt spinel nanocubes (h-Co3O4) corresponds actually to , whereas the calcined c-Co3O4 nanocubes assume the stoichiometric composition with <AOS> = 2.67. The determined cobalt <AOS> values for the Li-doped samples definitely confirm that in the case of the impregnated lithium catalysts (i-LixCo3O4), part of the octahedral Co3+ cations is reduced to Co2+, which is accounted for by eq 3. In contrast, for the lithiated cobalt spinel obtained via one-pot hydrothermal synthesis (h-LixCo3–xO4), a fraction of the Co3+ cations is oxidized into Co4+, following eq 4a. As a result, both cobalt spinels that are lithiated in different ways exhibit dissimilar redox features, which can be easily controlled by the way Li is introduced.

The work function, Φ, measurements in a vacuum (Figure 3a) and in oxygen (reaction product of N2O decomposition) at 295 and 350 °C (Figure 3b1,b2) were next performed to characterize the redox properties of the bare (c-Co3O4) and selected lithiated (h-Li-2-Co and i-Li-2-Co) samples, along with the doubly promoted catalyst (i-K-5/Li-2-Co) of the highest deN2O activity among all the investigated samples (see below). Variation in the work function with the K and Li loading is shown in Figure S9 (ESI Section). In both cases, a volcano-type behavior of Φ as a function of the lithium or potassium content was observed. The nonmonotonous changes in the work function for surface K doping are typically observed, and can well be accounted for by the Topping model2 [and references therein]. The origin of the minimum in the work function for the lithiated samples (Figure S9b) can be associated with the reduction of cobalt in the low concentration range (see Figure 2d), followed by the gradual development of LiCoO2 surface entities as the Li content increases above the critical value (see Figure S4).

Figure 3 Work function measurements for bare c-Co3O4, singly doped h-Li-2-Co and i-Li-2-Co, and doubly doped i-K-5/Li-2-Co spinel samples measured in vacuum (a) and in the presence of O2 at 295 °C (b1) and 350 °C (b2). Comparison of the experimental (blue line) and DFT-calculated (orange line) work function variations for the investigated catalysts (c).

For the impregnated catalysts, the presence of Co2+ in octahedral positions is well reflected by the lowering of the work function from Φ = 4.98 eV (c-Co3O4) down to 4.79 eV (i–Li-Co-2) and to 4.8 eV for the i-K-5/Li-2-Co catalysts when measured in a vacuum. For h-Li-2-Co, a slight work function enhancement to 5.04 eV remains in line with the formation of octahedral Co4+ holes. The Φ values for the bare c-Co3O4 and the lithiated spinels distinctly increase (ΔΦ = 0.19–0.1 eV) when measured at 295 °C in the presence of dioxygen, i.e., the reactant molecule of the highest electron affinity (EA = 0.45 eV79). On the contrary, for the doubly doped catalyst, the corresponding change was significantly smaller (ΔΦ = 0.03), the reason for which is explained below. When passing to a higher temperature (350 °C), the ΔΦ values augment further to ΔΦ = 0.33–0.2 and 0.04 eV for the singly and doubly doped catalysts, respectively.

In analogy to the previous treatment,80 the observed variation of the work function can be discussed in terms of the following factors: Φ = −EF + eΔV0S + ΔΦdip. The last two terms gauge the electrostatic barrier for electron removal to the vacuum level of zero energy, caused by the inherent surface potential of the bare surface (ΔV0S), and the Helmholtz surface potential created by the charged (ionosorbed) oxygen adspecies (ΔΦdip = e θads μdip/εε0, where μdip is the magnitude of the surface dipole, θads is the areal density of anionic oxygen adspecies, and ε and ε0 is the electric permittivity of the environment and vacuum, respectively81).

The formation, electronic structure, and stability of the negatively charged reactive oxygen adspecies, such as O2– and O22– on the (100) surface of cobalt spinel below 295 °C, were previously described by us in detail.61 An electron acceptor character of the adsorption gives rise to the formation of surface dipoles (ΔΦdip) with the protruding negative charge, leading to the observed enhancement of the work function, as reported elsewhere.82,83 The observed increase in the Φ values upon rising the temperature to 350 °C, in turn, results from partial dissociation of the diatomic O2z– into the monatomic O– adspecies of a larger dipole moment (see below).82,84 The much smaller sensitivity of the work function changes to the oxygen presence in the case of the doubly doped i-K-5/Li-2-Co catalyst can easily be explained by the formation of O–Kδ+ surface dipoles of an inversed polarity, which counterbalance the surface potential produced by Onz– adspecies.

As a result, by proper selection of the lithium doping technique, the redox state of the cobalt cations can intentionally be tuned owing to the valence pinning effect. The changes in the electronic structure and electrostatic state of the spinel surface upon doping and adsorption of the electron acceptor molecules, ascertained by XAFS and work function measurements, were next accounted for by DFT modeling. The experimental work function values are reasonably well reproduced by the calculation results (see Figure 3c), confirming definitely the proposed interpretation.

3.2 Catalytic Performance

For probing the redox behavior of the bare and Li/K-doped spinels, a catalytic N2O decomposition was used as a suitable model reaction. Neglecting, for the sake of simplicity, its involved mechanism (see Chapter S4 in ESI for more details), the reaction proceeds in two essential steps of prime dissociative reduction, epitomized as 2N2O + 2e– → 2O–ads + 2N2(g), and successive oxidative oxygen evolution reaction, epitomized as 2O–ads → O2(g) + 2e–,2,25,31,52 probing directly the redox properties of the catalyst (e– stands for the electron of the Fermi level energy). Such account remains in accordance with the apparent reaction order of 0.95–0.83 (see Figure S12c, and the associated discussion), where the leading reaction consists of the kinetically most demanding N2O dissociation. Its course is, however, significantly influenced by the subsequent oxygen evolution process, deviating the apparent reaction order from unity. Noteworthily, though both molecular events are sensitive to the applied doping, the overall reaction mechanism is still maintained, and the first step remains the rate-determining.

The N2O conversion curves for the lithiated Co3O4 samples measured in the TPSR mode as a function of the Li doping level are shown in Figure 4a. In the case of the i-LixCo3O4 series, a nonmonotonous dependence of the N2O conversion on the lithium content was found, with the optimal value for the i-Li-2-Co sample. All of the catalysts in this series were more active compared to the parent c-Co3O4, in contrast to the h-LixCo3–xO4 catalysts, which are noticeably less active and also less sensitive to the doping level. The corresponding Arrhenius plots are collated in Figure S11, along with the variation of the activation energy with an increasing concentration of lithium. While for the h-LixCo3–xO4 catalysts, the ΔEa values increase steadily with the Li doping level, in the case of the i-LixCo3O4 series, it passes through a well-pronounced minimum for the i-Li-2-Co sample, revealing a dramatically opposite catalytic behavior in response to the spinel redox tuning, which is governed by the route the lithium dopant was introduced and its concentration (Figure S11a’).

Figure 4 TPSR conversion curves of the N2O decomposition for the lithiated cobalt spinel catalysts obtained via impregnation (i-LixCo3O4 series) and one-pot hydrothermal synthesis (h-LixCo3-xO4 series) (a), together with the profiles of the potassium impregnated (i-Ky/Co3O4 series) (b) and the Li and K impregnated catalysts (i-Ky/Li0.045Co3O4 series) (c). The dotted lines correspond to the reference bare c-Co3O4 catalyst.

For the most active catalyst i-Li-2-Co and the reference c-Co3O4, the conversion dependence on the N2O pressure, pN2O, and on the contact time, τ, was additionally measured under the steady-state conditions (Figure 5a,b) at various temperatures (300 °C–500 °C). These results allow for accurate determination of the rate constants, k, and their temperature dependence at various surface coverages controlled by pN2O. According to our previous studies,56,60 within the mean-field approximation, a simple relation, k(T)·τ = X(T)[1 – X(T)], can be used for this purpose, for the type of the applied quasi-CSRT reactor. Next, a more detailed analysis of the temperature (and the surface coverage) dependence of the rate constants, k(T)θ, was performed for the two limiting cases corresponding to pN2O = 50 hPa (θ = max in the investigated pressure range) and pN2O → 0 hPa (θ = 0). The latter values were determined by extrapolating the pressure dependence of the conversion at constant contact time and temperature to 0 hPa (see Figure S12a,b). This approach allows for the determination of the intrinsic conversion value at zero surface coverage (X0, θ = 0), and the N2O conversion value, Xθ, at the surface coverage corresponding to the selected pressure, pN2O = 50 hPa. These values may then easily be converted to the rate constants k(T)θ, as described above. The corresponding plots are shown in Figure 5a1,a2,b1,b2 for bare and lithiated Co3O4 (i-Li-2-Co), respectively.

Figure 5 N2O conversion as a function of p(N2O) and the contact time for c-Co3O4 (a) and for the most active impregnated i-Li-2-Co (i-Li0.045Co3O4) catalyst (b). The corresponding X/(X – 1) plots are versus the contact time at various temperatures at pN2O = 0 Pa (a1,b1) and 50 hPa (a2,b2). The associated Arrhenius plots are shown in the insets.

The activation energies, ΔEa(pN2O), determined from the Arrhenius plots (see the corresponding insets in Figure 5a1,a2,b1,b2) for pN2O = 50 hPa and pN2O → 0 hPa are equal to ΔEa = 1.03 ± 0.05 and 0.96 ± 0.04 eV for c-Co3O4 and ΔEa = 1.16 ± 0.07 and 0.66 ± 0.07 eV for the i-Li-2-Co catalyst, respectively. Thus, in the case of the most active i-Li-2-Co catalyst, there is a striking sensitivity of the activation energy value to the N2O pressure, which is almost doubled when passing from pN2O → 0 to 50 hPa. The corresponding changes in the activation energy for the reference bare cobalt spinel are apparently less pronounced. Such large variations in ΔEa for N2O decomposition, repeatedly encountered in the literature,1,2,57,58 can be associated with the formation of negatively charged oxygen intermediates. The resultant buildup of a surface electrostatic field strongly influences the dynamics of the interfacial electron transfer processes. The latter constitute the very nature of the N2O decomposition redox, and this issue is addressed in more detail in the next section.

The effect of potassium doping on cobalt spinel activity is shown in Figure 4b. The observed nonmonotonous behavior of the catalyst performance as a function of the potassium loading is in agreement with our previous finding.2,25 It is associated with similar changes in the work function (see Figure S9a), which are well accounted for by the Topping equation.2 The influence of potassium doping on the catalytic performance of the best lithiated sample (i-Li-2-Co) is shown in Figure 4c. There is a pronounced synergy resulting from double doping with the most active catalyst, i-K-5/i-Li-2-Co, where the full conversion of N2O is already reached at T < 350 °C, while for the reference bare Co3O4 at this temperature, the conversion still remains below 20%.

The clear-cut nanocubic morphology of the investigated spinel catalysts (see Figure 1) allows for a reliable translation of the N2O conversion at the representative temperature of 350 °C into the turnover frequencies (eq 2). The resultant TOF values were plotted as a function of the doping level for the i-LixCo3O4, h-LixCo3–xO4, i-Ky/Co3O4, and i-Ky/Li0.045Co3O4 series (Figure 6a–d). The i-LixCo3O4 catalysts are characterized by a distinct volcano-type dependence of TOF on the Li/Co ratio (Figure 6a), in contrast to the h-LixCo3–xO4 series, where the TOF values are steadily decreasing in a rather flat way (Figure 6b).

Figure 6 TOF values of the N2O decomposition at 350 °C for i-LixCo3O4 (a), h-LixCo3–xO4 (b), i-Ky/Co3O4 (c), and i-Ky/LixCo3O4 (d) series of the cobalt spinel catalysts. The upper profile (marked green) in the panel (c) indicates the changes in the work function with surface potassium concentration (it represents a section of the complete plot of Φ, shown in Figure S9 in ESI).

The corresponding TOF profiles for the i-Ky/Co3O4 and i-Ky/Li0.045Co3O4 catalysts are plotted in Figure 6c,d. Additionally, the top panel in Figure 6c shows the changes in the work function, Φ, with the increasing surface potassium concentration. Like in the case of i-LixCo3O4, the TOF values for i-Ky/Co3O4 are also changing in a volcano fashion with increasing potassium loading. These changes are directly correlated with the work function in such a way that the lower the work function Φ (easier the electron transfer), the higher the activity (TOF value).

The performance of doubly promoted i-Ky/LixCo3O4 catalysts significantly outperforms those of the singly doped i-LixCo3O4 and i-Ky/Co3O4 counterparts, revealing a clear synergy of the lithium and potassium doping when both are introduced by the impregnation. The origin of this effect and the appearance of the TOF maximum are discussed in Chapter 3.4.

3.3 DFT Investigations into Structure-Redox Relationships

The optimized geometry of the parent bulk Co3O4 is shown in Figure S1, and the selected structural parameters are collated in Table S4. For modeling the incorporation of lithium into the cobalt spinel matrix, we examined four models of different lithium loci (Figure 7) to account for the two ways of Li doping, described in eqs 3 and 4a, 4b. In the case of the h–Li-Co3O4 catalysts, two possibilities of Li incorporation in the framework, tetrahedral or octahedral sites, were considered (see Figure 7a1,a2). The HSE06/DFT modeling confirmed the Li+ state of the lithium dopant, with the tetrahedral 8a (Li’8a) localization being more stable than the octahedral 16d (Li”16d) one, in accordance with the experimental findings.85 The corroborative OCC-MAT calculations indicated, in turn, that Li’8a generates an electron hole, which is preferentially localized on the 16d cobalt, producing Co4+ (Co·16d) species. The formation of Co4+ is in accordance with the increased <AOS> value obtained from XAFS measurements (Figure 2) and the work function enhancement described above (Figure 3). The alternative location of the hole on an adjacent tetrahedral Co2+ cation (Co·8a) is disfavored by 1.25 eV. In the case of the less stable incorporation of Li in the octahedral sites, the created Li”16d defect (Figure 7a2) is charge compensated by two holes, which are preferentially localized on the adjacent octahedral cations (formation of 2Co·16d species). Other possible configurations, such as Co·8a/Co·16d and 2Co·8a, are higher in energy (by 0.18 and 0.86 eV, respectively).

Figure 7 Location possibilities of the Li ions in the cobalt spinel host and the generated electron and hole-bearing defect centers in the case of the lithiated spinel catalyst obtained via one-pot synthesis (a1,a2) and by postsynthetic impregnation (b1,b2).

The two scenarios of Li incorporation through the impregnation method (i-LixCo3O4 series) are illustrated in Figure 7b. The calculations reveal that the Li location in the framework 8a position (Li’8a, marked yellow in Figure 7b1) with the concomitant displacement of Co2+ from the 8a to 16c sites (Co··16c, marked in purple) is energetically more favorable than the location of Li at the interstitial 48f position (formation of Li·48f species, marked green in Figure 7b2) by 0.77 and 0.95 eV for the charge-balancing electrons accommodated at the 16d (Co’16d) or 8a (Co’8a) sites, respectively. For the preferred Li’8a center, the reduction of the octahedral framework Co3+ cations (Co’16d) is preferred over the interstitial octahedral Co2+ (Co·16c) and the framework tetrahedral Co2+ (Co’8a) cations by 1.69 and 2.65 eV, respectively (Figure 7b1).

The particular mechanisms of Li incorporation, ascertained by the OCC-MAT/DFT calculations, lead to a good agreement (within the error <5 ± 1%) between the experimental <AOS> values of cobalt derived from the XAFS plot (Figure 2d) and the <AOS> values calculated based on the established defect structure of the corresponding spinel catalysts [2.89 (XAFS) vs 2.72 (DFT) for h-Li0.086Co2.94O4 and 2.55 (XAFS) vs 2.65 (DFT) for i-Li0.045Co3O4]. This finding remains also in accordance with the observed and calculated work function changes as well (see Figure 3c,b).

In conclusion, the OCC-MAT/DFT computational results resolve the locus of Li dopants and the resulting valence pinning of the cobalt cations. In the case of the catalyst obtained via one-pot hydrothermal synthesis, the formation of pair {Li’8a, Co·16d} defects takes place upon the incorporation of lithium into the spinel matrix. The lithiation of spinel via postsynthesis impregnation leads, in turn, to the formation of triplet {Li’8a, Co’16d, Co··16c} entities, upon the substitution of Li+ for Co2+. This is accompanied by the displacement of the tetrahedral Co2+ cations into the 16c interstitial octahedral sites.86

For the lithiated catalyst obtained by one-pot hydrothermal synthesis (h-Lix(Co3-xO4)), the formation of the {Li’8a, Co·16d} defects decreases the work function due to the EF energy lowering, while for the Li-doped spinel obtained via postsynthesis impregnation, the formation of {Li’8a, Co’16d, Co··16c} leads to an increase in the EF level and the resultant lowering of the work function.

In the case of the (100)Co3O4 surface covered by the charged monatomic and diatomic oxygen intermediates, the observed increase in the ΔΦ value (Figure 3b,c) results from the formation of the surface dipoles, (ΔΦdip = e θads μdip/εε0), as already mentioned above. The drop in the ΔΦ value after potassium addition is caused by the inversion of the dipole moment (μ = 0.44 D) with respect to that produced by the anionic oxygen adspecies. The dipole moment resulting from the ionosorbed oxygen intermediates (calculated from the plane-averaged charge density difference along the z-direction87) is weaker (μ = 1.28 D) for the diatomic O2– adspecies with the charge spread over both atoms, increasing significantly for the monatomic O– intermediates (μ = 2.02 D). This explains properly the corresponding changes in the work function upon the temperature increase in terms of thermal dissociation of superoxide adspecies, which almost doubles the surface dipole moment. Notably, the sign and size of the surface dipoles produced are fully consistent with the observed and calculated work function changes (see Figure 3 and the associated discussion).

3.4 Conceptual Account of the Catalyst Redox Behavior

The experimental results obtained, complemented by the calculated DOS structure of the spinel (100) surface and the energies of the FMO levels of N2O and O2, allowed us to construct a conceptual model of the redox interactions between the catalyst and the reacting molecules, build upon the tenets of molecular structure-based approaches to adsorption and catalysis.41,43,48 It provides a rational background for the elucidation of the effect of doping and charge buildup on the catalyst surface in the course of N2O decomposition. In the proposed treatment, the conceivable redox steps of nitrous oxide decomposition5

are restricted to the two molecular events of an initial dissociative reduction of N2O (entrance into the catalytic cycle through an interfacial electron transfer, Figure 8), and a terminal evolution of O2 (exit from the catalytic cycle upon back electron transfer, Figure 9). The in-between molecular surface events of diffusive association of the O–(ads) intermediates into transient O22–(ads), and their subsequent oxidative transformation into O2–(ads) adspecies have been previously discussed by us,61 while a detailed molecular resolution of these surface steps will be elucidated in a forthcoming article.

Figure 8 Molecular cascade of events in the course of N2O dissociation induced by an interfacial electron transfer, and the alignment of the energy levels of the prime molecular orbitals of N2O and the cobalt active centers. DOS structure of the (100) surface (gray shadowing) with the superimposed p-DOS of 3dz2 (marked in blue) and the FMO energy levels of the linear (marked in red) and bent (orange) N2O molecule in a distal configuration (a1), the key interaction diagram of the dz2-10a′ molecular orbitals before electron transfer (a2), constructed based on the DOS structure in a proximal configuration (a3), and the corresponding -pCOHP profile (a4). DOS structure in the proximal configuration upon electron transfer (a5), together with DOS of the resulting O–ads intermediates produced upon N2O– dissociation (a6). The spin density contours before (a2) and after electron transfer (a5) are colored pale green and pale violet, whereas the partial charge densities are colored pale yellow.

Figure 9 DOS structure of the (100) surface (marked by gray shadowing with the 3dz2 p-DOS component in marked blue) with the O2– adspecies in the proximal configuration (marked in red) (a1). Molecular orbital interaction diagram for dz2-2π||*in the proximal configuration before (a2) and in the distal configuration after the back electron transfer (a3). The corresponding -pCOHP profiles are shown in (b1) and (b2), respectively. The spin density contour of SOMO (a1) is colored pale violet and green, whereas the partial charge densities are pale yellow.

Redox interactions at the N2O|(100)Co3O4 interface governing the catalyst activity depend critically on the surface DOS structure and the actual manifold of the FMO energy levels of the reactants and intermediates involved. To drive the N2O decomposition reaction, the surface cobalt active sites act as persistent interfacial electron shuttling centers and molecular templates, providing orbitals of proper energy and symmetry for N2O capture and efficient electron transfer. Because the neutral N2O is linear (C∞v), whereas the N2O– transient is bent (Cs), according to the Franck–Condon (FC) principle (the nuclear positions remain frozen during ET), the instantaneous geometries of both species must match to make the electron transfer feasible. This is associated with a significant internal reorganization energy of the N2O molecule [λ(N2O) = 2.3 eV]88, the distortion of which in a gas phase is used here as an explanatory example (see Figure S13). Upon the bending, resulting from the vibronic ν2 excitation (C∞v → Cs transformation), the acceptor 3π* LUMO is significantly stabilized and becomes a 10a′ orbital (Figure S13a1). This is reflected in a remarkable enhancement in the electron affinity of N2O [from −2.1 (VEA) to = 0.22 eV (AEA)],52 which favors electron transference. The bending of N2O also leads to its inherent preactivation (Figure S13a2), manifested by a notable accumulation of the negative charge in the oxygen moiety (qO = −0.2 |e|), and a convergence of the charge in both nitrogen atoms at the common value qN = 0.1 |e|. These changes are accompanied by a substantial lengthening of the NN–O bond (by ∼0.3 Å), defining the actual dissociation coordinate of nitrous oxide. Once the electron is captured in the 10a′ orbital of the bent N2O molecule, the potential energy surface becomes dramatically flattened (cf. Figure S13b1,b2), making the resulting N2O– transient species highly unstable. The FC locus for electron transfer is defined by the lowest energy crossing point between the N2O and N2O– potential energy surfaces (see Figure S13b3, blue circle). Owing to the shallow minimum of the N2O– potential energy surface, this transient species splits into N2 and O– moieties spontaneously since the FC point is situated well above the N2O– dissociation barrier. This means that the NN–O bond-breaking process is inherently induced by ET, and a significant part of the activation energy for the N–O bond-breaking results from the geometrical reorganization of the N2O molecule [λ(N2O)], which is imposed by the Franck–Condon constraint.

As reported by us, the sustainable active sites for the decomposition of N2O over the Co3O4 are constituted by the octahedral cobalt (Co16d) cations.35 The DOS structure of the surface (100) (see Figure S14) shows that the pentacoordinate Co3+ active sites (C4v) assume an open shell (dxy2, dyz2, dxz1, dz21) configuration. In the case of the dz2 orbitals protruding toward the reactants, a strong electron exchange locates the empty α-dz2 state 1.8 eV above the EF level, whereas the occupied β-dz2 counterpart situated below EF is significantly broadened, ranging from −0.5 to −8 eV, with the maximum around −6 eV. For N2O reduction, the requisite electron has to be accommodated in LUMO (3π*), which is located too high in energy with respect to the Fermi level to make ET thermodynamically feasible. The corresponding lineup of the (100) surface DOS with the frontier orbitals of the N2O reactant is shown in Figure 8.

In the distal configuration, the HOMO (2π) level of N2O is deeply (−6 eV) plunged in the spinel valence band (VB), whereas the antibonding LUMO (3π*) is situated in the conduction band (Figure 8a1, see DOS features marked in red). Both are then redox inactive, and since E(LUMO) ≫ EF prevents the electron transfer to occur, the N2O molecule remains intact. However, vibrational excitation (ν2 = 596.3 cm–1) leads to N2O bending (3π* → 10a′ transformation), which strongly stabilizes the energy of the in-plane LUMO (see DOS features marked in yellow), as discussed above. Within the constraints imposed by the MO symmetry and energy matching rules, the 10a′ orbital of the bent N2O molecule can overlap most efficiently with the close-lying empty α-dz2 orbital of the cobalt active site (marked in blue). Following the Hoffman treatment of bonding on surfaces89 in a proximal configuration, such interaction leads to the formation of virtual bonding (ϕb = 10a’ + dz2) and antibonding (ϕa = dz2 – 10a′) states. The resultant molecular orbital diagram is shown in Figure 8a2, where the positions and energy dispersion of the ϕb and ϕa levels were determined from the calculated DOS structure of the N2O|Co3O4 system in the proximal configuration (Figure 8a3). Their bonding and antibonding character was ascertained by using the corresponding pCOHP diagram (Figure 8a4). Hitherto, a small interaction between the occupied β-dz2 and β-10a′ states contributes to an incipient association of the N2O molecule, similarly to the weaker out-of-plane overlap between the 3a” and 3dyz orbitals, but these interactions were neglected for the sake of simplicity. The extent of ϕa stabilization is proportional to the S2/ΔχN2O term (where S is the overlap integral and ΔχN2O is the energy difference between the parent dz2 and 10a′ levels in the distal position). The ϕb state, although initially empty (see the corresponding spin density contour in Figure 8a2 with the spurious amount of ρspin on the N2O moiety only), plays a critical role in N2O activation. The position of this energy level relative to the Fermi energy gauges for the amount of energy needed to trigger the interfacial electron transfer. An increasing hybridization (∼S2/ΔχN2O) between the 10a′(N2O) and 3dz2 Co orbitals, induced by shrinking the Co-ON2 distance, can drive the virtual ϕb below EF (at d(Co-ON2) ∼ 1.7 Å), making ET thermodynamically feasible. This is accompanied by dramatic changes in the hybridization of the resultant (ϕ′b)1 orbital, caused by the lengthening of the O–NN bond (from 1.26 to 1.67 Å), and incipient severance of the dinitrogen molecule (see the developing in-plane π|| and out-of-plane π⊥ orbitals of the nascent N2 in Figure 8a5). At this stage of the reaction, the ϕ′b state, apart from the dominant parent 2pz (0.38) and 3dz2 (0.1) contributions, also acquires considerable 2px(0.14) and 3dxz (0.04) character. Such dramatic changes in the molecular structure of the N2O moiety, caused by the electron capture, are reflected in a pronounced enhancement of the total IpCOHP value of the Co–N2O–Co bond from −3.6 to −9.1 eV. The neutral N2O molecule is more strongly attached via the N atom (IpCOHP = −2.28 eV) than the O atom (IpCOHP = −1.33 eV), whereas in the case of the N2O– ligand, the O–Co bond becomes slightly stronger than the N–Co one (−4.84 eV vs −4.29 eV). In particular, the interfacial electron transfer results in a spectacular enhancement of the IpCOHP value for the key α-dz2-2pz(O–Co) interaction realizing this process, which increases from −0.26 to −1.99 eV. The lifetime of the N2O–(ads) transient, τ(N2O–) ∼ 400 fs, estimated from the hybridization function of LUMO using the Newns-Anderson model,72 is sufficiently large in comparison to the NN–O stretching time,τN–O ∼ 5 fs, which triggers successful dissociation of the N2O– molecule. Since within this approach, the N2O– lifetime is determined by the broadening of the (ϕ′b)1 state due to the local interactions with the 3d orbitals of the cobalt adsorption site, it is not directly sensitive to the distal/neighboring dopants. Notably, the estimated value of τ(N2O–) is of the same order in magnitude as that reported for N2O– trapped on a titania catalyst (100 fs).50

The corresponding spin density plot (Figure 8a5) confirms the pronounced accumulation of the unpaired electron on the oxygen atom of the N2O– transient, consistent with the SOMO character of ϕ′b. Upon final detachment of the N2 moiety, ϕ′b evolves into the occupied pz(O)–dz2 and px(O)–dxz (ϕ”b) states, revealing a strong electronic relaxation of the redox core once the electron transfer process is fully accomplished. After N2 release, the unpaired electron is shifted into the ϕ”b orbital of the ensuing O– adspecies (Figure 8a6), while the pz(O)–dz2 part of the original 10a’–3dz2 interaction, left after N2O dissociation, becomes deeply immersed in the VB band.

As mentioned above (see eq 5), the evolution of the final O2(g) is initiated by the recombination of the two O–(ads) intermediates into a transient O22–(ads) adspecies, which is next oxidized into O2(g) by a gradual release of the 2 electrons back to the cobalt active sites (closing the redox cycle). The orbital interaction diagram for the last step of this process, O2–(ads) O2(g), is shown in Figure 9, where the key role is played by the singly occupied bonding orbital, ψb, that is located just below the EF level, and resulting from an overlap between the α-2π||* orbital of the superoxide intermediate and the α-3dz2 orbital of cobalt (see Figure 9a1,a2).

The bonding character between Co–O2–, and the antibonding character between O–O is ascertained by analysis of the pCOHP profiles shown in Figure 9b1,b2, respectively. Since the ψb level is located below EF, the back electron transfer can be achieved by decreasing the overlap term (S2/ΔχO2), realized simply by lengthening the Co–O2– distance above 2.4 Å. The mutual DOS-MO alignment of the released O2 molecule (in distal configuration) and the α-3dz2 state of the restored cobalt active center is shown in Figure 9a3. The orthogonal β-2π⊥* overlapping with β-3dxz corresponds to SOMO (see Figure 9a1, and the corresponding spin density spatial distribution). Thus, the electronic configuration of the bound O2– moiety can be formulated as (β-2π⊥*)1(β-2π||*)1(α-2π||*)1(α-2π⊥*)0. As a result, the spin-polarized back electron transfer (removal of an electron from α-2π||*) leads directly to dioxygen release in its triplet ground state 3Σ–g.

The conceptual diagram for the redox interaction of a N2O molecule with the (100) surface of the parent, singly, and doubly doped Co3O4 is shown in Figure 10. For the determination of the local vacuum level (Evac) for the bare and covered (100) surface of Co3O4, we calculated the plane-averaged electrostatic potential profiles along the 100 direction, Vxy(z).87 The distance between the Fermi level (EF) and Evac, shown in Figure 10a, corresponds to the work function for the bare surface, ΦDFT = 4.8 eV, which agrees quite well with the experimental value of 4.98 eV (Figure 3c). The alignment of the LUMO levels for the linear (3π*||) and bent (10a′) N2O reactant, and the SOMO of the O2 product (2π*) with respect to the valence and conduction bands of the (100) Co3O4 surface (their composition is represented by the 3D partial charge densities, averaged in the range of 0÷0.5 eV below VBM and above the CBM), is shown in Figure 10b. At distal configuration, the MO energy levels of N2O(g) and O2(g) are referenced to a common local Evac level. The distance of LUMO (N2O) and SOMO (O2–) to the catalyst Fermi energy position controls the feasibility of forth/back electron transfer that turns the catalytic cycle. According to the Fermi–Dirac distribution, it is also influenced by the temperature of the reaction. As a result, the EF energy, which is equivalent to the chemical potential of the electrons in the spinel catalyst, plays a key threshold function in these processes.

Figure 10 Conceptual redox diagram for N2O decomposition on cobalt spinel catalysts. Plane-averaged electrostatic potential profile, Vxy(z), defining the local vacuum level (Evac), the work function for the bare Co3O4 (a), the surface DOS structure with the energy levels of the redox orbitals of N2O, O2, and O2– (b), changes in the surface potential and the Fermi level, caused by the accumulation of charged oxygen intermediates on the surface and the doping of spinel by Li and K (c).

As explained above, the LUMO level of N2O is lowered upon bending and hybridization with the cobalt 3dz2 orbital (∼S2/ΔχN2O) to reach the threshold level of EF, where ET can occur spontaneously. Above EF, the electron transfer may arise as a thermally activated process only. Thus, an apparent ΔEa value depends, i.a., on the amount of the 3dz2–10a′ overlap, controlled by the N2O approach to the surface cobalt active sites, and the reorganization energy needed for N2O bending [λ(N2O)].

For the evolution of dioxygen, which closes the catalytic cycle, the (3dz2-2π*||) level must be elevated above EF to transform the O2–(ads) intermediates into the O2(gas) molecules via back electron transfer. This can be achieved by lowering the hybridization extent with the catalyst surface (S2/ΔχO2), simply by increasing the Co–O2 distance above 1.6 Å. The uneven gap between the 3dz2-10a’ (blue arrow) and 3dz2-2π*|| levels (orange arrow) and the Fermi energy implies that the reduction of N2O should energetically be more demanding than the oxidation of the O2–(ads) intermediate into the final dioxygen product, in a straight accordance with the apparent reaction order close to one (see Figure S12c and the associated text in the ESI Section).

The described picture corresponds to N2O decomposition for pN2O → 0, where the surface coverage by the charged oxygen intermediates (O– and O2–) is negligible. In the case of pN2O = 50 hPa, the observed increase of the activation energy by 0.5 eV (from ΔEa = 0.66 to 1.16 eV, see Section 3.2) is associated with the buildup of the surface potential ΔVs,O ≈ 0.49 eV (reflected by an increase of the work function value) (ΔΦdip = ΔVs,O) due to partial coverage of the surface by the anionic oxygen intermediates. The resultant shift of the local Evac level moves the 10a′ energy away from EF (marked by the blue arrow in Figure 10c), making the interfacial electron transfer more demanding. This leads to the apposite augmentation of the activation energy for N2O dissociation by 0.5 eV, as experimentally observed.

The key role of the Fermi level in the redox behavior of the spinel catalyst is further substantiated by taking into account the electronic effects of the Li and K doping, described in Sections 3.1 and 3.3. In the case of the h-LixCo3-xO4 series obtained by hydrothermal synthesis, the generated {Li’8a, Co·16d} species (defects) lower the EF level (Φ increases). This augments the barrier for the ET, in accordance with the observed lower activity of these catalysts, which steadily varies in parallel with Φ (see TOF values in Figure 6b and Ea in Figure S11).

The opposite effect is observed for the i-LixCo3O4 series obtained by the impregnation, which leads to the generation of the {Li’8a, Co’16d, Co··16c} entities. The latter raises the EF level, decreasing the barrier of the fourth electron transfer, which is reflected in the significant enhancement of the TOF values in comparison to the bare spinel (see Figure 6a). However, the catalytic activity pattern is actually more involved, passing through the maximum with increasing Li/Co ratio. This phenomenon can be accounted for by a trade-off between two opposite effects. Enhancement of the EF level increases the deN2O activity, which favors the accumulation of the ionosorbed oxygen intermediates (by shifting the reaction to lower temperatures where the oxygen adspecies are more stable). However, the latter effect by increasing the ΔVS,O value hinders the ET activation of N2O, lowering the TOF values for higher Li/Co ratios (see Figure 6a). Furthermore, the presence of the interstitial cobalt (Co··16c) cations increases the energy of the dioxygen desorption from 0.65 to 0.91 eV for bare Co3O4 (Co16d-O2-Co16d) and i-LixCo3O4 (Co16d-O2-Co16c) catalysts, respectively (see Figure S15), favoring thereby an accumulation of the charged oxygen intermediates on the catalyst surface.

The resultant surface potential buildup of the oxygen intermediates can be mitigated by the addition of potassium. The created surface dipole (μ = 0.44 D) due to potassium adspecies is of the opposite direction to that produced by the anionic oxygen intermediates (μ = 1.28 and 2.02 D for O2– and O–, respectively), lowering the energy separation between the 10a′ LUMO level of N2O and EF (marked by the red arrow), see Figure 10c. Within this account, the nonmonotonous changes of the TOF values for the i-Ky/Co3O4 catalyst (Figure 6c) result from the parallel variation of the work function upon doping (see the green profile in the upper panel).

The volcano dependence on K loading is repeated in the case of the doubly doped i-Ky/LixCo3O4 catalysts of the top activity (see Figure 6d). The appearance of the TOF maximum can be traced back to the synergy between the beneficial enhancement of the EF level (Li doping) and the lowering of the surface potential built-up (K-doping). The left part in Figure 6d corresponds to the work function changes controlled by the Topping model, in full analogy to the singly K-doped spinel (Figure 6c). Yet, excessive loading of potassium leads to an apparent steric blocking of surface active sites, and the TOF values drop eventually below the value observed for the unpromoted spinel (see Figure 6d right part).

4 Conclusions

The redox properties of cobalt spinel can be tuned on purpose by aliovalent doping with Li ions. The one-pot hydrothermal lithiation leads to the generation of {Li’8a, Co·16d} species (defects) that lower the EF level and the deN2O activity in a commensurate way. Postsynthetic lithiation via impregnation leads to the generation of the {Li’8a, Co’16d, Co··16c} species and volcano-type dependence of the TOF values on the Li content. The key factors governing redox activity of the cobalt spinel catalyst include: lineup of the MO energy levels of the reactants and the surface DOS structure, reorganization energy associated with the N2O bending to satisfy the Franck–Condon constraint, strength of the mutual interactions of the LUMO (N2O) and SOMO (O2–ads) with the 3dz2 (Co) orbitals, and buildup of surface potential due to accumulation of anionic oxygen intermediates (surface electrostatics). The established molecular orbital pattern and the position of the Fermi level control the trade-off between the interfacial electron transfer that triggers N2O dissociation and the back electron transfer that drives the subsequent dioxygen evolution. The crucial role of overlap between the virtual 10a′-3dz2 orbitals in attaining the transference of electrons between the N2O reactant (acceptor center) and the octahedral cobalt active sites (donor center) was shown for the first time.

The constructed conceptual framework allows for a concise account of the multifaceted redox features of the N2O decomposition reaction on cobalt spinels. It extends the classic picture of catalysis on semiconductors, pioneered by Wolkenstein and Hauffe, primarily by including a key hybridization of the molecular orbitals of the reacting species with the catalysts DOS features, surface electrostatics, and the reorganization energy of reactants. The proposed model of redox behavior can be applied to any p-type semiconductor catalysts for a straightforward rationalization of the molecular structure-redox reactivity relationships.

Supporting Information Available

The Supporting Information is available free of charge at https://pubs.acs.org/doi/10.1021/jacs.4c06587.Structure of Co3O4, slab models for cobalt spinel molecular modeling, characterization of the catalysts, catalytic performance, basic structural, and electronic characterization of Co3O4, Walsh diagram and potential energy surfaces of N2O and N2O–, DOS structure of the (100) surface of Co3O4, pDOS profiles of the exposed truncated octahedral Cooct5C cations, and energetic profiles for oxygen evolution on bare and Li-doped Co3O4 (PDF)

Supplementary Material

ja4c06587_si_001.pdf

Author Present Address

⊥ Max Planck Insitute for Chemical Energy Conversion, Stiftstr. 34-36, 45470 Mülheim an der Ruhr, Germany

The authors declare no competing financial interest.

Acknowledgments

This work was supported by Project OPUS No. 2017/27/B/ST4/01155, financed by the Polish National Science Center. We acknowledge the Polish computing infrastructure PLGrid for access to the LUMI supercomputer owned by the EuroHPC JU, and hosted by CSC, Finland, and the LUMI consortium for the computing grant no. 465001229. A.W. and J.S. acknowledge the funding from the Polish Ministry and Higher Education project: “Support for research and development with the use of research infrastructure of the National Synchrotron Radiation Centre SOLARIS” under contract nr 1/SOL/2021/2. J.S. acknowledges the National Science Centre in Poland under grant number 2020/37/B/ST3/00555.
==== Refs
References

Konsolakis M. Recent Advances on Nitrous Oxide (N2O) Decomposition over Non-Noble-Metal Oxide Catalysts: Catalytic Performance, Mechanistic Considerations, and Surface Chemistry Aspects. ACS Catal. 2015, 5 , 6397–6421. 10.1021/acscatal.5b01605.
Wójcik S. ; Grzybek G. ; Stelmachowski P. ; Sojka Z. ; Kotarba A. Bulk, surface and interface promotion of Co3O4 for the low-temperature N2O decomposition Catalysis. Catalysts 2020, 10 , 41 10.3390/catal10010041.
Wu Y. ; Zheng Y. ; Granger P. Past, Present, and Future in the Development of Medium and High-Temperature Catalytic Processes for N2O Decomposition. Advances in Catalysts Research; Advances in Material Research and Technology; Springer: Cham, 2004; Vol. 7 , pp 221–257.
Zabilskiy M. ; Djinović P. ; Tchernychova E. ; Tkachenko O. P. ; Kustov L. M. ; Pintar A. Nanoshaped CuO/CeO2 materials: Effect of the exposed ceria surfaces on catalytic activity in N2O decomposition reaction. ACS Catal. 2015, 5 , 5357–5365. 10.1021/acscatal.5b01044.
Ruszak M. ; Inger M. ; Witkowski S. ; Wilk M. ; Kotarba A. ; Sojka Z. Selective N2O removal from the process gas of nitric acid plants over ceramic 12CaO • 7Al2O3 catalyst. Catal. Lett. 2008, 126 , 72–77. 10.1007/s10562-008-9619-x.
Dacquin J. P. ; Lancelot C. ; Dujardin C. ; Da Costa P. ; Djega-Mariadassou G. ; Beaunier P. ; Kaliaguine S. ; Vaudreuil S. ; Royer S. ; Granger P. Influence of preparation methods of LaCoO3 on the catalytic performances in the decomposition of N2O. Appl. Catal., B 2009, 91 , 596–604. 10.1016/j.apcatb.2009.06.032.
Richards N. ; Carter J. H. ; Parker L. A. ; Pattisson S. ; Hewes D. G. ; Morgan D. J. ; Davies T. E. ; Dummer N. F. ; Golunski S. ; Hutchings G. J. Lowering the operating temperature of perovskite catalysts for N2O decomposition through control of preparation methods. ACS Catal. 2020, 10 , 5430–5442. 10.1021/acscatal.0c00698.
Franken T. ; Palkovits R. Investigation of potassium doped mixed spinels CuxCo3-xO4 as catalysts for an efficient N2O decomposition in real reaction conditions. Appl. Catal., B 2015, 176–177 , 298–305. 10.1016/j.apcatb.2015.04.002.
Chellam U. ; Xu Z. P. ; Zeng H. C. Low-Temperature Synthesis of MgxCo1-xCo2O4 Spinel Catalysts for N2O Decomposition. Chem. Mater. 2000, 12 (3 ), 650–658. 10.1021/cm990355l.
Maniak G. ; Stelmachowski P. ; Stanek J. J. ; Kotarba A. ; Sojka Z. Catalytic properties in N2O decomposition of mixed cobalt-iron spinels. Catal. Commun. 2011, 15 , 127–131. 10.1016/j.catcom.2011.08.027.
Obalová L. ; Karásková K. ; Jirátová K. ; Kovanda F. Effect of potassium in calcined Co-Mn-Al layered double hydroxide on the catalytic decomposition of N2O. Appl. Catal., B 2009, 90 , 132–140. 10.1016/j.apcatb.2009.03.002.
Karásková K. ; Obalová L. ; Jirátová K. ; Kovanda F. Effect of promoters in Co-Mn-Al mixed oxide catalyst on N2O decomposition. Chem. Eng. J. 2010, 160 , 480–487. 10.1016/j.cej.2010.03.058.
Chmielarz L. ; Kuśtrowski P. ; Drozdek M. ; Rutkowska M. ; Dziembaj R. ; Michalik M. ; Cool P. ; Vansant E. F. SBA-15 mesoporous silica modified with rhodium by MDD method and its catalytic role for N2O decomposition reaction. J. Porous Mater. 2011, 18 , 483–491. 10.1007/s10934-010-9401-5.
Tabor E. ; Jíša K. ; Nováková J. ; Bastl Z. ; Vondrová A. ; Závěta K. ; Sobalík Z. Role of Pt(0) in bimetallic (Pt,Fe)-FER catalysts in the N2O decomposition. Micropor. Mesopor. Mater. 2013, 165 , 40–47. 10.1016/j.micromeso.2012.07.035.
Abu-Zied B. M. ; Schwieger W. ; Unger A. Nitrous oxide decomposition over transition metal exchanged ZSM-5 zeolites prepared by the solid-state ion-exchange method. Appl. Catal., B 2008, 84 , 277–288. 10.1016/j.apcatb.2008.04.004.
Li S. ; Zhang C. ; Li J. ; Li Y. ; Wang H. ; Li C. ; Song Y. Direct catalytic decomposition of N2O over Co(x)/RPSA Catalysts. Res. Chem. Intermed. 2019, 45 , 3601–3616. 10.1007/s11164-019-03811-6.
Dacquin J. P. ; Dujardin C. ; Granger P. Surface reconstruction of supported Pd on LaCoO3: Consequences on the catalytic properties in the decomposition of N2O. J. Catal. 2008, 253 , 37–49. 10.1016/j.jcat.2007.10.023.
Parres-Esclapez S. ; Illán-Gómez M. J. ; de Lecea C. S.-M. ; Bueno-López A. On the importance of the catalyst redox properties in the N2O decomposition over alumina and ceria supported Rh, Pd and Pt. Appl. Catal., B 2010, 96 , 370–378. 10.1016/j.apcatb.2010.02.034.
Zabilskiy M. ; Arčon I. ; Djinović P. ; Tchernychova E. ; Pintar A. In-situ XAS study of catalytic N2O decomposition over CuO/CeO2 catalysts. ChemCatChem 2021, 13 , 1814–1823. 10.1002/cctc.202001829.
Jing Y. ; Taketoshi K. ; Zhang N. ; He C. ; Toyao T. ; Maeno Z. ; Ohori T. ; Ishikawa N. ; Shimizu K. Catalytic decomposition of N2O in the presence of O2 through redox of Rh oxide in a RhOx/ZrO2 catalyst. ACS Catal. 2022, 12 , 6325–6333. 10.1021/acscatal.2c01321.
Xue L. ; Zhang C. ; He H. ; Teraoka Y. Catalytic decomposition of N2O over CeO2 promoted Co3O4 spinel catalyst. Appl. Catal., B 2007, 75 , 167–174. 10.1016/j.apcatb.2007.04.013.
Russo N. ; Fino D. ; Saracco G. ; Specchia V. N2O catalytic decomposition over various spinel-type oxides. Catal. Today 2007, 119 , 228–232. 10.1016/j.cattod.2006.08.012.
Li R. ; Li Y. ; Liu Z. Recent advances in the catalytic removal of NOx and N2O over spinel oxide-based catalyst. Fuel 2024, 355 , 129405 10.1016/j.fuel.2023.129405.
Gryboś J. ; Hudy C. ; Gryczynska A. ; Piskorz W. ; Sojka Z. Hydrothermal Synthesis of Euhedral Co3O4 Nanocrystals via Nutrient-Assisted Topotactic Transformation of the Layered Co(OH)2 Precursor under Anoxic Conditions: Insights into Intricate Routes Leading to Spinel Phase Development and Shape Perfection. Cryst. Growth Des. 2020, 20 (12 ), 7771–7787. 10.1021/acs.cgd.0c01031.
Zasada F. ; Stelmachowski P. ; Maniak G. ; Paul J.-F. ; Kotarba A. ; Sojka Z. Potassium promotion of cobalt spinel catalyst for N2O decomposition-accounted by work function measurements and DFT modelling. Catal. Lett. 2009, 127 , 126–131. 10.1007/s10562-008-9655-6.
Maniak G. ; Stelmachowski P. ; Zasada F. ; Piskorz W. ; Kotarba A. ; Sojka Z. Guidelines for optimization of catalytic activity of 3d transition metal oxide catalysts in N2O decomposition by potassium promotion. Catal. Today 2011, 176 , 369–372. 10.1016/j.cattod.2010.11.043.
Omata K. ; Takada T. ; Kasahara S. ; Yamada M. Active site of substituted cobalt spinel oxide for selective oxidation of CO/H2. Part II. Appl. Catal., A 1996, 146 , 255–267. 10.1016/S0926-860X(96)00151-2.
Yan L. ; Ren T. ; Wang X. ; Ji D. ; Suo J. Catalytic decomposition of N2O over MxCo1-xCo2O4 (M = Ni, Mg) spinel oxides. Appl. Catal., A 2003, 45 , 85–90. 10.1016/S0926-3373(03)00174-7.
Yan L. ; Ren T. ; Wang X. ; Gao Q. ; Ji D. ; Suo J. Excellent catalytic performance of Zn Co1–Co2O4 spinel catalysts for the decomposition of nitrous oxide. Catal. Commun. 2003, 4 , 505–509. 10.1016/S1566-7367(03)00131-6.
Ohnishi C. ; Asano K. ; Iwamoto S. ; Chikama K. ; Inoue M. Alkali-doped Co3O4 catalysts for direct decomposition of N2O in the presence of oxygen. Catal. Today 2007, 120 , 145–150. 10.1016/j.cattod.2006.07.042.
Grzybek G. ; Stelmachowski P. ; Gudyka S. ; Duch J. ; Ćmil K. ; Kotarba A. ; Sojka Z. Insights into the twofold role of Cs doping on deN2O activity of cobalt spinel catalyst-towards rational optimization of the precursor and loading. Appl. Catal., B 2015, 168–169 , 509–514. 10.1016/j.apcatb.2015.01.005.
Sun J. ; Song A. ; Tian Y. ; Zhan H. ; Deng J. ; Wang H. ; Ke M. Unravelling the Effect of Alkali Metal Depositionon Co3O4 for Catalytic Decomposition of N2O. ChemCatChem 2023, 15 , e20220144 10.1002/cctc.202201449.
Stelmachowski P. ; Maniak G. ; Kotarba A. ; Sojka Z. Strong electronic promotion of Co3O4 towards N2O decomposition by surface alkali dopants. Catal. Commun. 2009, 10 , 1062–1065. 10.1016/j.catcom.2008.12.057.
Obalová L. ; Maniak G. ; Karásková K. ; Kovanda F. ; Kotarba A. Electronic nature of potassium promotion effect in Co–Mn–Al mixed oxide on the catalytic decomposition of N2O. Catal. Commun. 2011, 12 , 1055–1058. 10.1016/j.catcom.2011.03.017.
Stelmachowski P. ; Maniak G. ; Kaczmarczyk J. ; Zasada F. ; Piskorz W. ; Kotarba A. ; Sojka Z. Mg and Al substituted cobalt spinels as catalysts for low temperature deN2O—Evidence for octahedral cobalt active sites. Appl. Catal., B 2014, 146 , 105–111. 10.1016/j.apcatb.2013.05.027.
Zhecheva E. ; Stoyanova R. ; Angelov S. Doping of Co3O4 with lithium by a solid-state reaction in air I. Oxidation degree and coordination of cations. Mater. Chem. Phys. 1990, 25 , 351–360. 10.1016/0254-0584(90)90124-s.
Bhattacharya J. ; Wolverton C. Relative stability of normal vs. inverse spinel for 3d transition metal oxides as lithium intercalation cathodes. Phys. Chem. Chem. Phys. 2013, 15 , 6486–6498. 10.1039/c3cp50910a.23529669
Wu X. ; Scott K. A Li-doped Co3O4 oxygen evolution catalyst for nonprecious metal alkaline anion exchange membrane water electrolysers. Int. J. Hydrogen Energy 2013, 38 , 3123–3129. 10.1016/j.ijhydene.2012.12.087.
Wolkenstein T. The electron theory of catalysis on semiconductors. Adv. Catal. 1960, 12 , 189–264. 10.1016/S0360-0564(08)60603-3.
Eichelbaum M. ; Hävecker M. ; Heine C. ; Wernbacher A. M. ; Rosowski F. ; Trunschke A. ; Schlogl R. The electronic factor in alkane oxidation catalysis. Angew. Chem., Int. Ed. 2015, 54 , 2922–2926. 10.1002/anie.201410632.
Haber J. ; Witko M. Oxidation catalysis—electronic theory revisited. J. Catal. 2003, 216 , 416–424. 10.1016/S0021-9517(02)00037-4.
Nandakumar N. K. ; Seebauer E. G. Relating catalytic activity of d0 semiconducting metal oxides to the Fermi level position. J. Phys. Chem. C 2014, 118 , 6873–6881. 10.1021/jp4127946.
García-Mota M. ; Vojvodic A. ; Abild-Pedersen F. ; Nørskov J. K. Electronic origin of the surface reactivity of transition-metal-doped TiO2(110). J. Phys. Chem. C 2013, 117 , 460–465. 10.1021/jp310667r.
Tao H. B. ; Fang L. ; Chen J. ; Yang H. B. ; Gao J. ; Miao J. ; Chen S. ; Liu B. Identification of surface reactivity descriptor for transition metal oxides in oxygen evolution reaction. J. Am. Chem. Soc. 2016, 138 , 9978–9985. 10.1021/jacs.6b05398.27441842
Xu D. ; Qu W. ; Gu X. ; Li M. ; Liu R. ; Dong Y. ; Ma Z. ; Chen Y. ; Tang X. An activity descriptor for perovskite oxides in catalysis. Chem Catal. 2022, 2 , 1163–1176. 10.1016/j.checat.2022.03.011.
Gao D. ; Deng P. ; Zhang J. ; Zhang L. ; Wang X. ; Yu H. ; Yu J. Reversing free-electron transfer of MoS2+x cocatalyst for optimizing antibonding-orbital occupancy enables high photocatalytic H2 evolution. Angew. Chem., Int. Ed. 2023, 62 , e202304559 10.1002/anie.202304559.
Simböck J. ; Ghiasi M. ; Schönebaum S. ; Simon U. ; de Groot F. M. F. ; Palkovits R. Electronic parameters in cobalt-based perovskite-type oxides as descriptors for chemocatalytic reactions. Nat. Commun. 2020, 11 , 652 10.1038/s41467-020-14305-0.32005805
Yang C. ; Grimaud A. Factors controlling the redox activity of oxygen in Perovskites: From theory to application for catalytic reactions. Catalysts 2017, 7 , 149 10.3390/catal7050149.
Yu L. ; Yan Q. ; Ruzsinszky A. Key role of antibonding electron transfer in bonding on solid surfaces. Phys. Rev. Mater. 2019, 3 , 092801(R) 10.1103/physrevmaterials.3.092801.
Cheng C. ; Prezhdo O. V. ; Long R. ; Fang W.-H. Photolysis versus photothermolysis of N2O on a semiconductor surface revealed by nonadiabatic molecular dynamics. J. Am. Chem. Soc. 2023, 145 , 476–486. 10.1021/jacs.2c10643.36541604
Gerischer H. Electron-transfer kinetics of redox reactions at the semiconductor/electrolyte contact. A new approach. J. Phys. Chem. 1991, 95 , 1356–1359. 10.1021/j100156a060.
Pietrzyk P. ; Zasada F. ; Piskorz W. ; Kotarba A. ; Sojka Z. Computational spectroscopy and DFT investigations into nitrogen and oxygen bond breaking and bond making processes in model deNOx and deN2O reactions. Catal. Today 2007, 119 , 219–227. 10.1016/j.cattod.2006.08.054.
Kobayashi M. ; Kobayashi H. Decomposition kinetics of nitrous oxide on manganese dioxide—application of transient response method. J. Chem. Soc. Japan 1974, 6 , 438–442. 10.1252/jcej.6.438.
Pérez-Ramırez J. ; Kapteijn F. ; Schöffel K. ; Moulijn J. A. Formation and control of N2O in nitric acid production: where do we stand today?. Appl. Catal., B 2003, 44 , 117–151. 10.1016/S0926-3373(03)00026-2.
Obalova L. ; Fila V. Kinetic analysis of N2O decomposition over calcined hydrotalcites. Appl. Catal., B 2007, 70 , 353–359. 10.1016/j.apcatb.2005.11.031.
Piskorz W. ; Zasada F. ; Stelmachowski P. ; Kotarba A. ; Sojka Z. DFT modeling of reaction mechanism and ab initio microkinetics of catalytic N2O decomposition over alkaline earth oxides: from molecular orbital picture account to simulation of transient and stationary rate profiles. J. Phys. Chem. C 2013, 117 , 18488–18501. 10.1021/jp405459g.
Yu H. ; Wang X. Apparent activation energies and reaction rates of N2O decomposition via different routes over Co3O4. Catal. Commun. 2018, 106 , 40–43. 10.1016/j.catcom.2017.12.004.
Wang Y. ; Wei X. ; Hu X. ; Zhou W. ; Zhao Y. Effect of formic acid treatment on the structure and catalytic activity of Co3O4 for N2O Decomposition. Catal. Lett. 2019, 149 , 1026–1036. 10.1007/s10562-019-02681-2.
Duan Y. ; Zhang Q. ; Song Z. ; Wang J. ; Tang X. ; Liu Q. ; Zhang T. Effect of preparation methods on the catalytic activity of Co3O4 for the decomposition of N2O. Res. Chem. Intermed. 2017, 43 , 7241–7255. 10.1007/s11164-017-3071-8.
Grzybek G. ; Gryboś J. ; Indyka P. ; Janas J. ; Ciura K. ; Leszczyński B. ; Zasada F. ; Kotarba A. ; Sojka Z. Evaluation of the inhibiting effect of H2O, O2, and NO on the performance of laboratory and pilot K-ZnxCo3-xO4 catalysts supported on α-Al2O3 for low-temperature N2O decomposition. Appl. Catal., B 2021, 297 , 120435 10.1016/j.apcatb.2021.120435.
Zasada F. ; Piskorz W. ; Janas J. ; Budiyanto E. ; Sojka Z. Dioxygen activation pathways over cobalt spinel nanocubes – from molecular mechanism into ab initio thermodynamics and 16O2/18O2 exchange microkinetics. J. Phys. Chem. C 2017, 121 , 24128–24143. 10.1021/acs.jpcc.7b09597.
Julien C. ; Camacho-Lopez M. A. ; Mohan T. ; Chitra S. ; Kalyani P. ; Gopukumar S. Combustion synthesis and characterization of substituted lithium cobalt oxides in lithium batteries. Solid State Ionics 2000, 135 , 241–248. 10.1016/S0167-2738(00)00370-2.
Lubeck J. ; Beckhoff B. ; Fliegauf R. ; Holfelder I. ; Hönicke P. ; Müller M. ; Pollakowski B. ; Reinhardt F. ; Weser J. A novel instrument for quantitative nanoanalytics involving complementary X-ray methodologies. Rev. Sci. Instrum. 2013, 84 , 045106 10.1063/1.4798299.23635231
Zasada F. ; Piskorz W. ; Cristol S. ; Paul J.-F. ; Kotarba A. ; Sojka Z. Periodic density functional theory and atomistic thermodynamic studies of cobalt spinel nanocrystals in wet environment: molecular interpretation of water adsorption equilibria. J. Phys. Chem. C 2010, 114 , 22245–22253. 10.1021/jp109264b.
Kresse G. ; Joubert D. From ultrasoft pseudopotentials to the projector augmented-wave method. Phys. Rev. B 1999, 59 , 1758–1775. 10.1103/physrevb.59.1758.
Selcuk S. ; Selloni A. DFT+U study of the surface structure and stability of Co3O4(110): Dependence on U. J. Phys. Chem. C 2015, 119 , 9973–9979. 10.1021/acs.jpcc.5b02298.
Hafner J. Ab-initio simulations of materials using VASP: Density-functional theory and beyond. J. Comput. Chem. 2008, 29 , 2044–2078. 10.1002/jcc.21057.18623101
Monkhorst H. J. ; Pack J. D. Special points for Brillouin-zone integrations. Phys. Rev. B 1976, 13 , 5188–5192. 10.1103/PhysRevB.13.5188.
Press W. H. ; Flannery B. P. ; Teukolsky S. A. ; Vetterling W. T. Numerical Recipes in Pascal. The Art of Scientific Computing; Cambridge University Press: Cambridge, NY, 1989 [ISBN 0–521–37516–9].
Allen J. P. ; Watson G. W. Occupation matrix control of d- and f-electron localisations using DFT + U. Phys. Chem. Chem. Phys. 2014, 16 , 21016–21031. 10.1039/c4cp01083c.24832683
Nelson R. ; Ertural C. ; George J. ; Deringer V. ; Hautier G. ; Dronskowski R. LOBSTER: Local orbital projections, atomic charges, and chemical-bonding analysis from projector-augmented-wave-based density-functional theory. J. Comput. Chem. 2020, 41 , 1931–1940. 10.1002/jcc.26353.32531113
Ghan S. ; Diesen E. ; Kunkel Ch. ; Reuter K. ; Oberhofer H. Interpreting ultrafast electron transfer on surfaces with a converged first-principles Newns–Anderson chemisorption function. J. Chem. Phys. 2023, 158 , 234103 10.1063/5.0151009.37318168
Hadjiev V. G. ; Iliev M. N. ; Vergilov I. V. The Raman spectra of Co3O4. J. Phys. C: Solid State Phys. 1988, 21 , 199–201. 10.1088/0022-3719/21/7/007.
Yao Z. ; Kim S. ; Aykol M. ; Li Q. ; Wu J. ; He J. ; Wolverton C. Revealing the conversion mechanism of transition metal oxide electrodes during lithiation from first-principles. Chem. Mater. 2017, 29 , 9011–9022. 10.1021/acs.chemmater.7b02058.
Baker M. L. ; Mara M. W. ; Yan J. J. ; Hodgson K. O. ; Hedman B. ; Solomon E. I. K- and L-edge X-ray absorption spectroscopy (XAS) and resonant inelastic X-ray scattering (RIXS) determination of differential orbital covalency (DOC) of transition metal sites. Coord. Chem. Rev. 2017, 345 , 182–208. 10.1016/j.ccr.2017.02.004.28970624
Van Oversteeg C. H. M. ; Doan H. Q. ; de Groot F. M. F. ; Cuk T. In situ X-ray absorption spectroscopy of transition metal based water oxidation catalysts. Chem. Soc. Rev. 2017, 46 , 102–125. 10.1039/C6CS00230G.27834973
Edla R. ; Braglia L. ; Bonanni V. ; Miotello A. ; Rossi G. ; Torelli P. Study of gaseous interactions on Co3O4 thin film coatings by ambient pressure soft X-ray Absorption Spectroscopy. J. Phys. Chem. C 2019, 123 , 24511–24519. 10.1021/acs.jpcc.9b05433.
Zasada F. ; Piskorz W. ; Sojka Z. Cobalt spinel at various redox conditions: DFT+U investigations into the structure and surface thermodynamics of the (100) facet. J. Phys. Chem. C 2015, 119 , 19180–19191. 10.1021/acs.jpcc.5b05136.
Chen E. S. ; Wentworth W. E. ; Chen E. C. M. The electron affinities of NO and O2. J. Mol. Str. 2002, 606 , 1–7. 10.1016/S0022-2860(01)00813-4.
Greiner M. T. ; Chai L. ; Helander M. G. ; Tang W.-M. ; Lu Z.-H. Transition metal oxide work functions: The influence of cation oxidation state and oxygen vacancies. Adv. Funct. Mater. 2012, 22 , 4557–4568. 10.1002/adfm.201200615.
Hofmann O. T. ; Deinert J.-C. ; Xu Y. ; Rinke P. ; Stahler J. ; Wolf M. ; Scheffler M. Large work function reduction by adsorption of a molecule with a negative electron affinity: Pyridine on ZnO(101̅0)(101̅). J. Chem. Phys. 2013, 139 , 174701 10.1063/1.4827017.24206316
Lin L. ; Jacobs R. ; Ma T. ; Chen D. ; Booske J. ; Morgan D. Work Function: fundamentals, measurement, calculation, engineering, and applications. Phys. Rev. Appl. 2023, 19 , 037001 10.1103/PhysRevApplied.19.037001.
Nowotny J. ; Bak T. Electron probe for surface science and surface defect engineering of oxide semiconductors for sustainable energy conversion. J. Phys. Chem. C 2020, 124 , 20617–20642. 10.1021/acs.jpcc.0c05816.
Zasada F. ; Piskorz W. ; Janas J. ; Gryboś J. ; Indyka P. ; Sojka Z. Reactive oxygen species on the (100) facet of cobalt spinel nanocatalyst and their relevance in 16O2/18O2 isotopic exchange, deN2O, and deCH4 processes—a theoretical and experimental account. ACS Catal. 2015, 5 , 6879–6892. 10.1021/acscatal.5b01900.
Stoyanova R. ; Zhecheva E. ; Angelov S. Doping of Co3O4 with lithium by a solid state reaction in air. III. EPR evidence of intrinsic disorder reactions in the tetrahedral interstitials. Mater. Chem. Phys. 1990, 26 , 239–244. 10.1016/0254-0584(90)90013-Z.
Casas-Cabanas M. ; Binotto G. ; Larcher D. ; Lecup A. ; Giordani V. ; Tarascon J.-M. Defect chemistry and catalytic activity of nanosized Co3O4. Chem. Mater. 2009, 21 , 1939–1947. 10.1021/cm900328g.
Liu D. ; Meng X. ; Qin G. Insights into unidirectional migration of photo-excited electrons at ZnFe2O4/graphene van der Waals interface. Comput. Mater. Sci. 2019, 157 , 60–66. 10.1016/j.commatsci.2018.10.039.
Sojka Z. ; Che M. EPR investigation of the activation of N2O on Mo/SiO2 catalysts via electron transfer: From N2O as a ligand to adsorbed O– ion. J. Phys. Chem. 1996, 100 , 14776–14785. 10.1021/jp960871o.
Hoffmann R. A chemical and theoretical way to look at bonding on surfaces. Rev. Mod. Phys. 1988, 60 , 601–628. 10.1103/RevModPhys.60.601.
