Constructing A Rational Kinetic Model of the Selective Propane Oxidation Over A Mixed Metal Oxide Catalyst

This research presents a kinetic investigation of the selective oxidation of propane to acrylic acid over a MoVTeNb oxide (M1 phase) catalyst. The paper contains both an overview of the related literature, and original results with a focus on kinetic aspects. Two types of kinetic experiments were performed in a plug flow reactor, observing (i) steady-state conditions (partial pressure variations) and (ii) the catalyst evolution as a function of time-on-stream. For this, the catalyst was treated in reducing atmosphere, before re-oxidising it. These observations in long term behaviour were used to distinguish different catalytic routes, namely for the formation of propene, acetic acid, acrylic acid, carbon monoxide and carbon dioxide. A partial carbon balance was introduced, which is a ‘kinetic fingerprint’, that distinguishes one type of active site from another. Furthermore, an ‘active site’ was found to consist of one or more ‘active centres’. A rational mechanism was developed based on the theory of graphs and includes two time scales belonging to (i) the catalytic cycle and (ii) the catalyst evolution. Several different types of active sites exist, at least as many, as kinetically independent product molecules are formed over a catalyst surface.


Introduction
Light alkenes are largely produced from crude oil feedstock by steam-cracking.These hydrocarbons are transformed into valuable intermediates, such as ethylene oxide, acrylic acid, and maleic anhydride.Current research efforts are undertaken to establish alkanes instead of alkenes as feedstock for the latter processes, hence, facilitating the change from a crude oil to natural gas-based economy in this respect.The synthesis of maleic anhydride from butane has already been realised in larger scale, whereas the analogous selective oxidation of propane and ethane still faces challenges in their performance to compete with the established alkene-based processes [1,2].
For the complex reaction of propane oxidation on oxide catalysts, particularly the MoVTeNbO x in the M1 structure, has been studied by many researchers under different conditions, including atmospheric conditions and high vacuum (cf.Table 1).
In this article, we will analyse this process as an example of the selective oxidation process over an MoVTeNb oxide catalyst.a PFR: Plug Flow Reactor; TAP: Temporal Analysis of Products; b balance to 100% is inert gas; c relative to total composition, balance nitrogen; d CO 2 0-50%; e relative to C 3 H 8 ; f this work.
As many catalytic heterogeneous reactions, the selective oxidation of propane over MoVTeNbO x phase, is a complex physico-chemical process.Three interwoven aspects have to be stressed: (a) it is a multi step process, (b) it occurs over the active catalytic site (centre) of complex structure, and (c) it is a process with multi-scale temporal behaviour.A categorisation is presented in Figure 1 to visualise different levels of complexity.(a) Multi step process: The main steps of this process are: (i) interactions between components of the reaction mixture and the catalytic surface, i.e., the impact reaction between the molecule of gaseous propane and catalyst oxygen, adsorption of water etc; (ii) mass transfer processes, i.e., oxygen exchange between the catalyst surface and bulk; (iii) transport of charge (electrons and/or ions) within the catalytic unit.

(b)
Complex structure of the active site/centre: The 'structure of the catalyst' tackles the nature of the catalyst surface; i.e., the architecture of catalytically active site/centre, how many active centres are considered per site, and furthermore, how many sites are distinguishable etc.An empirical concept of site isolation proposed by Grasselli [19][20][21] can be used as one of basic structure pre-assumptions for the rational kinetic model.This concept is presented as follows: "The catalytic reaction sequence should occur on an ensemble of atoms ('the active site') that holds together with the educt molecule all other reactants and exchanges electrons between the reactants.For these operations neither the influx of active species during the unit conversion time (as long as the educt is adsorbed) nor the participation of electrons from outside of the active site are required . . .The catalytic cycle is only finished when the site is regenerated into its initial active form, being here fully oxidised and holding all oxygen atoms belonging to its initial structure" [22].
In this concept, many additional processes are not considered, such as processes of the formation of active sites, processes of spill-over, rather to say, all processes external regarding the isolated sites.From the kinetic point of view, the process over the isolated site will be characterised by the specific relationship between the rate of reactants and products which participate in this process.However, these relationships are still not recognised.
In the literature, there are many considerations on the complex functional structure of the active sites which is responsible for generation of different products.It makes sense to distinguish: a 'whole active site' and 'active centre' (or 'centres') within this site, which generates (or generate) the specific product(s).One 'active site' consists of at least one 'active centre'.There are speculations about dimer-and trimer structures within the whole active site.These different centres can be considered as ingredients over which different catalytic routes are performed and different products are generated.In particular, it was assumed that an active site for oxidation of propane may consist of dimer-and trimer-ingredients.Dimers and more complex structures (trimers and even 'dimers of dimers') can store redox-equivalents as a part of the process of oxygen activation.
It is clear that structures of active sites and/or active centres mean surface structures.The surface of a catalyst material (e.g., a metal or metal oxide) may significantly differ from its bulk structure and composition, and may additionally be affected strongly by the gas phase surrounding it [12,[23][24][25][26][27].

(c)
Multi-scale temporal behaviour: As for multi-scale temporal behaviour of the catalytic process, at least three characteristic times have to be distinguished; 'fast' and 'moderate' intrinsic times of the catalytic cycle, and the 'slow' time, which is external.The latter is caused by the modification and evolution of the catalytic centre/site under the influence of reactive medium.In the catalytic literature, this phenomenon was described and conceptualised by Boreskov [28], see also Yablonskii et al. [29].A variety of physical methods for the investigation of the surface and bulk of catalysts and different methods of kinetic analysis, mostly steady state methods, have been used to investigate this process.
1.1.Investigations of the Selective Oxidation of Propane.Methods and Results

Physical Methods
Several methods were applied to characterise the catalyst before and after reaction to report their change.Micro calorimetric experiments showed a lower adsorption enthalpy of propane after reaction, suggesting a significant change in the nature of the catalyst surface [25].Furthermore, ESR (electron spin resonance) and IR (infra-red) investigations found different sites for propane activation and nucleophilic oxygen addition on MoO 3 catalysts, as well as different types of lattice oxygen and their role during the reaction are discussed [30].A combined LEIS (low-energy ion scattering), SEM (scanning electron microscopy) and TEM (transmission electron microscopy) investigation of M1 catalyst found that the (001) plane of the M1 crystals being not solely responsible for the selectivity during propane (amm)oxidation [31].
In situ experiments are very valuable for the characterisation of catalysts under working conditions.Most important, and still missing, is their link to kinetic interpretations.In situ XPS (X-ray photoelectron spectroscopy) investigation of MoVTeNbO x M1 catalysts at 623 K and dry and wet C 3 H 8 /O 2 feed found the production of acrylic acid correlated with the Te/V ratio and inversely correlated with the abundance of Mo [25,32].
Bulk and surface V content correlated with the selectivity to acrylic acid, determined by in situ XPS [5].
Experiments were performed in propane/oxygen feed with and without steam.The addition of steam was beneficial for the selectivity to acrylic acid, as it was found before, which was observed with an enrichment in V 5+ and Te 6+ , and depletion in Mo 6+ on the surface of the catalyst.Simultaneously, the conductivity decreased together with the work function and a modification of the valence band structure.A sub-surface space charge region was formed, which was depleted in electrons under reaction conditions.The extension of band bending was not affected by the addition of steam, hence no affect on the number of charge carriers in this space charge region.Near-edge X-ray absorption fine structure (NEXAFS) investigations point to a significantly modified electronic structure of the surface, especially the vanadyl bond motif (suggesting an enrichment of covalently bonded V 5+ ) [33].

Kinetic Methods. Model-Free Kinetic Data
For the investigation of the selective oxidation of propane; mostly the plug flow reactor (PFR) has been applied.Examples for these investigations are summarised in Table 1.
Typically, there was investigated an influence of the three reactants, i.e., propane, oxygen, and water.
The essential role of water in activating the catalyst is widely accepted.Without the presence of its large excess (in comparison to stoichiometric requirements), neither acetic acid nor acrylic acid would form on the catalyst.Centi et al. demonstrated the necessity of surface modification by water and/or ammonia prior to observable selective (amm)oxidation products [34].
Novakova et al. [6] determined apparent kinetic orders for propane disappearance, first order with respect to propane, 0.26 with respect to oxygen.The formation rate of propene is first order with respect to propane and 0.26 order with respect to oxygen.Formation of acrylic acid is close to first order in propane and one half with respect to oxygen.Water increased propane conversion and selectivity to oxygenated products.Space time variations show non-zero intercept for acrylic acid and acetic acid.Propene is considered a major intermediate, for the formation of acetic acid and acrylic acid.
The kinetic order was determined to be 0.88 and 0.06 with respect to propane and oxygen partial pressure, respectively [10].The independence of oxygen partial pressure was considered for all rates, if its concentration is above 1 vol-%.Below that threshold the re-oxidation becomes very slow and the oxygen partial pressure becomes kinetically relevant [10].
Several O 2 /C 3 H 8 and H 2 O/C 3 H 8 ratios have been investigated for testing the predictability of power law, Eley-Rideal (ER) and Mars-van-Krevelen (MvK) models.Mazloom et al. [11] analysed these models and found that the oxygen dependency was described most accurately by the MvK model.In addition, they concluded that under water free conditions the propane activation is rate determining, whereas under wet conditions it is the oxygen activation.Generally, most of researchers support the MvK mechanism for the selective oxidation reactions over mixed metal oxide catalysts [1,2,36,37].
Naumann d'Alnoncourt et al. investigated the propane oxidation over MoVTeNbO x in the M1 structure.As found before by Novakova et al. [6] steam is beneficial for acrylic acid production.The apparent activation energy was found to depend on the feed composition [12].
Summing up, one can say, that most of the kinetic data are likely obtained in the small conversion domain.The apparent kinetic order of propane was determined close to first order with respect to propane and zero or small with respect to oxygen.As for the influence of water, its beneficial character can be considered as a reasonable consideration.However, literature data regarding the kinetic order of water are still uncertain.

Mechanisms
Since the 1970s, in the literature there was discussion regarding the details of selective propane oxidation focused on main questions: • Is the mechanism parallel or consecutive, or parallel-consecutive, regarding the different gas products, in particular propene?• Does this complex reaction occur on a single site or multiple sites?
• What is the nature of an active site?
In early isotopic labelling experiments for the reaction of propene with H 18  2 O over a MoWSnTeO x catalyst, the incorporation of 18 O into the acrylic acid molecule was observed [15].
A comprehensive isotopic labelling (D and 13 C) investigation over Fe-Sb oxide catalysts was performed.Starting from propene; no significant H/D scrambling in the propene precursor could be observed, neither a double bond isomerisation.The slow step is the allyl formation from adsorbed propene [38].
Two different active sites were proposed, i.e., an acidic and an oxidizing site (cf. Figure 2), in investigation on BiMo-oxide catalysts [39,40].Either site may contribute to scrambling of the deuterium labelled carbon atom being oxidised.The rate determining step is the abstraction of the allyl hydrogen to form the π-allyl complex.Ammoxidation and oxidation of propene over molybdate and antimonate catalysts were investigated.For molybdates and oxidation, the allyl formation was considered as slow compared to other reaction steps, and the first oxygen insertion reversible [41].
Investigation of the role of water was performed by Saleh-Alhamed et al. on the partial oxidation of propene [42,43].Water suppresses the CO 2 formation, by blocking the most active sites and also keeps the catalyst at a high oxidation state.Water competes with acrylic acid for adsorption sites.
According to Barrault et al. [3] the addition of water changed selectivity, not the activity of the catalyst.Mo 5+ species are responsible for the formation of acetic acid, where as Mo 6+ species support the formation of acrylic acid, via acrolein.Both reaction pathways branch off from propene as a central intermediate product.
Furthermore, Barrault et al. [3] found a parallel reaction pathway for acetic and acrylic acid.Allylic oxidation is considered to occur over coordinatively unsaturated Mo sites in combination with lattice oxygen.If a site is highly reduced and slowly reoxidised with molecular oxygen, the lattice may not be able to provide enough lattice oxygen to achieve high selectivity for the allylic oxidation.
For the similar ODH of propane process, Chen et al. [13] proposed a first order reaction network (propane → propene, propane → CO x , propene → CO x ), and detailed mechanism for ODH over MoO x /ZrO 2 and VO x /ZrO 2 catalysts, and these mechanisms were similar.The initial C-H bond cleavage was to be rate determining.It was suggested, that the relative adsorption enthalpies of propene and propane influence the relative reaction rates of these molecules, as well as the Lewis acidity of the cations involved in the adsorption.Propane ODH and propane combustion show similar activation energies, hence, the fate of reaction is decided after the kinetically relevant step.Activation energies are higher for propene combustion.Additionally, a large number of isotopic labelling experiments was performed.
Lin et al. [44] investigated the selective oxidation of propane on several catalyst systems, mainly the mixed metal oxide MoVTeNbO x in the M1 structure, and proposed a first order reaction network [45].Also Lin et al. [44] proposed propene as the central intermediate from which the major pathway to acrylic acid via acrolein branches off, together with the minor pathway to acetic acid via acetone.
Luo et al. [46] studied the partial oxidation of propane, emphasising the differences in reaction pathways over zeolite VO-H-beta and MoVTeNb-oxide.The schemes in Figure 3 compare both catalysts suggesting the reaction in zeolite H-beta being driven by its acidity, whereas in MoVTeNb-oxide the oxidation of a C 3 -chain rather than the cracking and total oxidation was favoured.Furthermore, two sites for MoVTeNb-oxide are suggested; the first site catalyses the reaction propane-propene-acrylic acid, whereas the second supports the pathway propene-acetic acid.Centi et al. [34] proposed a reaction mechanisms for the selective (amm)oxidation reaction.It was pointed out the role of surface dynamics, i.e., the modification of the surface reactivity due to chemisorbed species (reactants, intermediates, and products), hence, the dynamics of catalyst reconstruction.Multiple pathways exist even for the same reaction product.
Vitry et al. [7] assumed at least three catalytic oxidation sites to be present on the catalyst surface (i.e., for propane activation, allylic oxidation, and acrolein oxidation) and supported this assumption by data on different kinetic orders.
A more comprehensive mechanism including acetic acid and propionic acid was presented by Novakova et al. [6].Their investigations were performed over a Mo-V-Sb-Nb oxide catalyst.The scheme is presented in Figure 4a.The order of CO 2 and CO formation with respect to oxygen is 1.1 and 0.8, respectively.Acrylic acid was found to be first order in propane and close to one half in oxygen [6].
Balcells et al. [47] employed a first-order power rate law network, where propene is considered to be the only primary intermediate, form which acrolein, acrylic acid and by-products are formed in parallel pathways.
Grabowski et al. compared two kinetic models, namely related to the Langmuir-Hinshelwood (LH) mechanism and Eley-Ridel (ER) mechanism.The latter one achieved a better description of the experimental data [48].Kinetics of the oxidative dehydrogenation of C 2 -C 3 alkanes was summarised in a comprehensive Grabowski article [49].
Based on the network by Lin et al. [45], Naraschewski et al. [10] modified the scheme of the reaction by introducing two additional reaction pathways (cf. Figure 4b).The rate equations were fitted to experimental data [10].It is assumed that only one active site associated with vanadium exists, and furthermore, that not all vanadium present is part of the catalytically active cycle.Hence, there are active and inactive sites containing vanadium.The initial C-H bond cleavage is considered to be rate determining, however, the subtle change in V concentration, hence, substitution for Mo in the M1 structure, affects the activity and selectivity of the catalyst.Not all V sites are catalytically active.Sites are distinguished for propane adsorption and oxygen insertion, and C-C bond cleavage.The formation of acetic acid via acetone can be lowered by decreasing the Brønsted acid sites [10].
Temporal analysis of products (TAP) pulse-response investigations were performed over Mo-V-Te-Nb oxide catalysts [9].The proposed reaction scheme is presented in Figure 5.The main features are the competition of water, propene and oxygen for the same vacancy.Acrylic acid may originate from propane directly (no re-adsorption of propene).Acetic acid, however, originates from re-adsorbed propene.
Selective oxidation of propane was investigated by Naumann et al. [12,50] leading to acrylic acid.An interesting peculiarity was found, i.e., a maximum in the selectivity of acrylic acid as a function of propane conversion.Furthermore, the dependence of apparent activation energies was determined as a function of reaction conditions.
Ammoxidation of propane provides insight into the reducibility of the MoVNbTeO x catalyst.The catalyst could be reduced with C 3 H 8 /NH 3 pulses (analogy to the results of reduction in C 3 H 8 /H 2 O feed of this work) up to a depth of 70 layers.The catalyst could be re-oxidised without structural collapse using O 2 pulses [14].The measurement of kinetic isotope effects (KIE) in propane oxidation revealed essential differences in the activation of the propane molecule over MoVTeNb oxide with M1 structure and silica-supported vanadium oxide.Whereas the abstraction of the first hydrogen atom at the methylene group of the C 3 H 8 molecule is rate limiting over silica-supported vanadium oxide, activation of propane over M1 clearly involves both methyl and methylene groups simultaneously in rate-limiting process(es) [51].The formation pathways of carbon oxides during propane oxidation were studied by means of pulse experiments with 13 C-labeled/unlabeled propane/oxygen mixtures.Over M1, oxidation of acrolein to acrylic acid via oxygen insertion into the C-H bond of the CHO group is faster than oxidation of the vinyl group followed by decarbonylation of the resulting surface species and acetaldehyde formation.In contrast, decarbonylation of acrolein is preferred over silica-supported vanadium oxide showing that surface structure and electronic properties of the catalyst have an impact on the reaction mechanism.

Catalyst Structure and Nature of the Active Site
Generally it is agreed, that vanadium sites are required for the activation of alkanes.However in their investigation, Holmberg et al. [52] found the importance of the chemical surrounding of such V-sites, thus, only the orthorhombic (not the hexagonal) phase of a MoVNbTeO x catalyst was active for paraffin activation.Furthermore, they found that both orthorhombic and hexagonal phases are selective for the ammoxidation of propene.
Vitry et al. [53] found that the addition of niobium to a MoVTeO x catalyst increased the selectivity to acrylic acid.
Grasselli et al. [14] investigated MoVNbTeO x catalysts in M1 and M2 structure, namely their physical mixture as particles and powders.The latter exhibit a much higher selectivity for the selective oxidation of propane.Furthermore, the redox behaviour was investigated.The catalyst changed its selectivity depending on its oxidation-reduction state.This reduction-oxidation process was found to be reversible.
Oxygen exchange between the catalyst bulk and surface is traditionally taken into account for the explanation of kinetic behaviour of selective oxidation of propane.As a development of this concept there was proposed, that a reactive site acts as a collective electron storage unit.[18] The first consideration is: this 'exchange' idea contradicts to the 'site isolation' concept.However, it maybe considered as an apparent contradiction: it is a subject of the further analysis.
The parameter for this exchange, which can be found by conductivity analysis, may be used as a coefficient in our kinetic model for oxygen exchange [12].

Aim
In the present work, we are aiming to develop a methodology for designing rational kinetic models employing the selective oxidation of propane over MoVTeNbO x catalysts as an example.Also we attempt to link features of temporal kinetic behaviour to mechanistic and structural aspects, proposing a kinetic fingerprint for the justification for the isolated site concept.

Characterisation
The catalyst has been characterised by several methods.One of the most important information is presented in Figure 6.Both the 'fresh' and 'used' catalyst (i.e., after at least the here presented kinetic investigations) show typical reflexes for the M1 structure.The reflex positions and width are, in comparison to the 'fresh' sample, preserved for the 'used' sample.The structural features of M1 phase remain unchanged, despite repeated reducing/oxidising treatment to the catalyst.The BET surface area was determined to be 7.31 m 2 g −1 .X-ray fluorescence (XRF) investigations gave an element distribution of Mo/V/Te/Nb/O = 46.27%/6.74%/6.82%/7.01%/33.2%.The catalyst appeared in highly crystalline rods of several 100 nm in length with an homogeneously distributed elemental composition, revealed by TEM and SEM characterisation, respectively.

Justification of Kinetic Data
The first step of kinetic analysis is testing the validity of obtained data using elemental balances.In this work it was done via checking the elemental balances for carbon and oxygen with an accuracy of ±3%.7 and 8.All product formation rates are proportional to propane partial pressure, most pronounced for the production of acrylic acid and propene.CO, acetic acid and CO 2 form a group of similar dependence (cf. Figure 7a).Only acrylic acid showed a slight positive dependence on the partial pressure of oxygen, whereas all other products were observed to be independent (cf. Figure 7b).The dependence on water partial pressure is presented in Figure 8. Acrylic acid, acetic acid and CO are clearly proportional with respect to water partial pressure.From the dependence of the latter two, it was immediately concluded, that the acetic acid production was accompanied by the production of CO, which is also shown in Equation ( 2).We distinguish CO, which originates solely from non-selective oxidation (cf.Equation ( 4)) and CO related to the formation of acetic acid, as presented in Equation (2).To evaluate the CO production separately from the production of acetic acid, the latter was subtracted from CO, determined according to Equation ( 6). [CO Thus, the turquoise bullets in Figures 7 and 8 represent the previously described formation rate of CO, independent of the CO production due to acetic acid formation.Table 2 summarises the apparent kinetic orders of product formation rates with respect to all three reactants, i.e., propane, oxygen, and water.The apparent kinetic order n app of a product formation rate R j was determined according to Equation (7).The apparent kinetic order of CO 2 , propene, and acrylic acid with respect to propane are similar, in the range from 0.6-0.8,however, their kinetic coefficients vary within one order of magnitude.The trends of CO 2 and CO-aca are rather similar, certainly with respect to oxygen and water, in which both kinetic coefficients and apparent kinetic orders are much alike.Propene and both carboxylic acids have similar apparent orders with respect to propane and oxygen, also their kinetic coefficients with respect to propane are rather similar.A significant difference is observed in their kinetics with respect to water.Propene clearly does not depend, whereas both acetic and acrylic acid show an apparent order of one half.Water is, thus, necessary for the formation of acetic and acrylic acid, whereas it is at the same time a reaction product, hence, to some extend we can call it auto-catalysis regarding water.
Overall, the sum of apparent kinetic orders is smaller or about one for CO 2 , CO, and propene.Acetic and acrylic acid apparent kinetic orders sum up to a value of 1.5.

Kinetic Experiments during Catalyst Evolution ('Slow Dynamics' Scale of Min and Hours)
The experimental data in Figures 9 and 10a was recorded after a treatment of the catalyst in reducing conditions (i.e., 12 h in C 3 H 8 /O 2 /H 2 O/N 2 = 3/0/10/87 at 573 K and 10 4 mL g −1 h −1 ).Subsequently, the oxygen feed was added (i.e., C 3 H 8 /O 2 /H 2 O/N 2 = 3/12/10/75 at 573 K and 10 4 mL g −1 h −1 ).The gas analysis was started in this moment, representing time zero.The chemical composition of gases from the PFR exit was measured every 20 min to observe the 'slow' process of the kinetic steady-state (rather say pseudo-steady-state) of the catalyst.Before the observed different trends are described, the non-destructive treatment in reducing conditions needs to be explained.Reducing conditions were applied several times according to the same procedure described before.The reproducibility of sequential treatment in reducing and oxidising conditions is presented in Figure 9.The same catalyst loading was treated for 2 h in reducing (i.e., propane/water/inert) feed for 2 h, before applying an oxidising feed (i.e., propane/oxygen/water/inert).This procedure was repeated after the catalyst has reached stable performance in terms of product evolution (i.e., not changing with time).Figure 9 confirms the reproducibility of such sequential treatment in reducing and oxidising conditions.Both experiments show essentially the same product evolution as a function of time-on-stream.Certainly, the initial state of the catalyst (at time zero) was affected by the length of pre-treatment time (compare Figures 9 and 10a  Clearly, the catalyst was affected by the pretreatment, since the catalyst regains, thus, recovers its original performance to produce acrylic acid with a high selectivity.In Figure A1 the course of propane conversion and product evolution during oxidising and reducing feed is shown.In reducing feed, the propane conversion drops rapidly (cf. Figure A1a).After 10 min the conversion approaches practically zero.The concentration of all products except acetone decreases, the concentration of the acids decreases faster than the concentration of CO and CO 2 (cf. Figure A1b).
Figure 10a shows the evolution of all product species as a function of time, after a treatment in reducing conditions for 12 h.Significant trends are observed.CO and CO 2 decreased in a similar fashion, with an simultaneous increase of acrylic acid.Both acetic acid and propene respond much faster to the suddenly added oxygen in the feed, in comparison to the formerly mentioned products.
A rapid increase was observed for acetic acid from the first to the second data point and decreasing subsequently.The concentration of propene increases slightly with time-on-stream.Furthermore, the ratio of (CO-aca)/CO 2 is 0.8 and constant as a function of time-on-stream.All concentrations achieved a plateau, eventually.

Data Interpretation
In our interpretation and discussion, we are going to combine data of both types of experiments with 'slow' and 'fast' times.First, we would like to stress the obvious fact, that we do have this 'slow' trend in kinetic dependences, and this tendency is exhibited for special conditions of the catalyst, i.e., for a reduced catalyst.As experiments show, this trend is very reproducible, and the catalyst, analysed after all experimental treatments, is not damaged, i.e., reduced vs. oxidised.The fact whether this 'slow' trend is observed or not, will only depend on reduction and oxidation degree of the catalyst.The existence of this trend, reveals the fact of some 'slow' process, which exist inside the whole catalytic process.Regarding the whole catalytic cycle (cycles), it is likely external, e.g., it may combine different catalytic cycles.
The plateau of all observed dependencies can be identified with the stage when the catalyst is already re-oxidised.Obviously all the data presented in Figure 10a can be attributed to the stage of catalyst evolution from its reduced state to the oxidised one.
Figure 10a contains information about a significant difference in the temporal behaviour of all products, regarding the TOS-scale.Three cases can be distinguished: • no temporal trend for propene concentration; • a strongly decreasing trend in CO and CO 2 concentration; and a strongly increasing trend in acrylic acid concentration; • some trend in acetic acid concentration, which can be considered between 'small' and 'significant', depending on the degree of reduction, which is governed by the pretreatment time.This trend was small if the pretreatment time was up to 2-3 h (cf. Figure A2).
Because of this difference in trends, we can hypothesise that in our process at least three different active cites are present.
The cite for CO, CO 2 , and acrylic acid formation was affected by oxidation, whereas the oxidation to propene was not affected at all.
In Figure 10b the mass balances are plotted as a function of time-on-stream.The propane conversion (named as all) is the sum of all components, and it is slightly decreasing.
The interesting feature observed during TOS-investigations was the following: at any time the relationship 3 see Figure 10b.It can be termed as a 'partial carbon balance'.This partial carbon balance may be employed as a fingerprint for the existence of the special 'isolated site'.(Generally speaking, this site can have a complex structure, combining different centres.)How to interpret this balance?It is not a total carbon balance, since this equation does not include concentrations of other carbon-containing substances, i.e., propane, propene and acetic acid.
It is possible that Equation (8) relates to the corresponding overall reaction in which acrylic acid is oxidised into CO, CO 2 , and water.
The partial balance in Equation ( 8) consists of the products of Equations ( 3)-( 5).Let's present Equation (3) with a factor of two, and Equations ( 4) and ( 5) together as follows.
Subtracting Equation ( 9) from Equation (10) gives In accordance with this stoichiometric equation, the following relationships can be written: Hence, This linear relationship in Equation ( 16) is a consequence of the stoichiometric Equation ( 11).However, a hypothesis about the oxidation of acrylic acid into CO and CO 2 via reaction (11), or even the formation of acrylic acid from CO and CO 2 (the reverse of reaction ( 11)), can be eliminated easily.It makes sense to say, that Reaction (11) has no physical meaning.Equation ( 11) reflects a hidden process of surface evolution, a bridge mechanism, as it will be explained later.
From our point of view, this intriguing invariant [CO] + [CO 2 ] + [aa]-balance (in Figure 10b) (the CO employed in this balances is only produced other than from the reaction forming acetic acid), which does not depend on the propane conversion, gives a hint on a 'closed' structure, within which the complete transformation of admitted substances (propane) occurs, and substances acrylic acid, CO, and CO 2 are produced.This structure can be identified with an isolated site.Generally, this conclusion relates to Grasselli's concept, that a reaction may occur on an isolated site, and phase cooperation plays a significant role in the process as well [19,20,39].Several catalysts are composed of different phases, which are proven to outperform as a mixture compared to the single phase catalysts.Site isolation is the distribution and isolation of active sites/components/centres from each other.However, since these products (acrylic acid from one side and CO/CO 2 from another ones) are characterised by different kinetic dependences, we hypothesize that this site has two centres, which are responsible for two different routes.
One of these centres is responsible for the production of CO and CO 2 , and the second one performs production of acrylic acid.Our hypothesis about these two processes (i.e., 'two processes'-'two centres') is justified by the significant difference in long term kinetic dependencies for acrylic acid, from one site, and CO and CO 2 , from another one.
The mentioned centres, say X and X , may be linked by a slow oxidation reaction.Therefore, such reaction will be responsible for slow kinetic behaviour presented in Figure 10a.The detailed transformation over the 'isolated site' is depicted schematically in Figure 11.The transformation X X is fed by [O (b/g) ], which represents oxygen from either the bulk catalyst phase or the gas phase.The possible sources of oxygen shall not be limited to one or the other phase here.
Since the slow kinetic evolution results in some steady-state with the 'non-zeroth' concentration, process of site transformation has to be considered reversible, and the schematic drawing will be modified.The X X transformation is 'slow' in comparison with the 'fast' cycles of product generation.Therefore, cycle A (CO and CO 2 production) and cycle B (acrylic acid production) can be presented and analysed separately.Slow dynamics of the separate 'α − β' transformation is qualitatively presented in Equation ( 17): where α and β are dimensionless concentrations of the two centres, respectively, and the balance equation Figure 11.Scheme of a bridge mechanism.The 'active site' consists of two 'active centres' X and X'.On centre X the total oxidation (cycle A) takes place, and the selective oxidation (cycle B) takes place on centre X'.X(TO) corresponds to the sum of all intermediates, which are part of cycle A (for total oxidation).X'(SO) corresponds to the sum of all intermediates, which are part of cycle B (for selective oxidation).The bridge mechanism is fed by [O (b/g) ], i.e., oxygen from the catalyst bulk or the gas phase.
Has to be fulfilled, with α = θ X (t) + θ X(TO) (t) and β = θ X (t) + θ X (SO) (t) and t is in the scale of time-on-stream.θ X (t) and θ X(TO) (t) are surface/bulk intermediates/phases which belong to the total oxidation (TO) cycle and θ X (t) and θ X (SO) (t) are intermediates/phases which belong to the selective oxidation (SO) cycle, respectively.

Rational Mechanism
Comparing data on both types of experiments, it is evident that long term-measurements (TOS experiments) exhibit dramatic differences between temporal behaviour of products, in comparison with short time experiments (partial pressure variation experiments).
In fact, partial pressure experiments do not show significant differences in the kinetic trends of the reaction products with respect to oxygen, i.e., all kinetic orders are zero.All product formation rates depend positively on propane partial pressure.The only exception is found for the influence of water partial pressure, both acetic and acrylic acids are positively depending on water partial pressure (auto-catalysis), whereas all other products (i.e., propene, CO, and CO 2 ) are independent.
Oppositely, in long term experiments, as mentioned before, the big difference in time-dependences for all product formation rates was observed, shown in Figures 9 and 10a the preliminary analysis above.We did not find a stoichiometric relationships, which relate these dependences to each other.There was clearly a strong coupling observed between the formation rate of acrylic acid and the formation rates of CO and CO 2 .At the same time, there is no kinetic correlation between other product formation rates: i.e., between (i) propene and acetic acid, (ii) propene and CO, CO 2 , (iii) propene and acrylic acid (iv) acetic acid and CO 2 .
Difference in long term behaviour for some products reflects a difference in specific catalytic routes of these product formation.Using this hypothesis, we are proposing the following mechanism.It has such routes: (a) route of propene formation, (b) route of acetic acid formation, and (c) and (d) routes of CO 2 and CO formation bridged by a slow oxidation step to the route of acrylic acid formation.
Route of CO and CO 2 formation have one joined step, in which propane is consumed and some specific steps, in which CO and CO 2 are produced, respectively.CO produced in route (c) and also in route (b), together with acetic acid.Therefore, every specific product is characterised by a specific route over a specific centre (principle 'specific product-specific centre').There is an obvious advantage of such a principle, i.e., simplicity in writing the steady-state (or pseudo-steady-state) rate equation.
From Figures 7-10a it can be easily seen, that the ratio of (CO-aca)/CO 2 is roughly around one, at least not strictly bound to each other.For general purpose, CO and CO 2 are produced by different centres, but our experimental data is not enough to distinguish these.It can, however, be concluded, that CO and CO 2 do not originate from the same process.Otherwise, a constant 2:1 or 1:2 ratio of (CO-aca)/CO 2 would be expected, since propane has three carbon atoms.Thus, propane forms solely CO via one route, and solely CO 2 via another route.
Characterising our mechanisms, we would like to stress, that all mechanisms are rational.Mechanisms for propene and CO and CO 2 formation, are two step mechanisms, i.e., they are minimal catalytic cycles.As for mechanisms for the production of acetic and acrylic acid, they are four step cycles, containing one step of interaction with propane, one step for the reaction with water and two steps for the reaction with oxygen.In our opinion, it is a reasonable level of simplification.Also, these two mechanisms for the formation of acrylic and acetic acid, have a buffer step (step out of the cycle), in which water is a participant.As indicated in the introduction, the presence of water is necessary to activate the cycles for the production of selective oxidation products.Additionally, based on theoretical kinetic considerations [54] we need to have a reversible buffer step of the interaction with water to obtain a sum of kinetic orders greater than one, which is in accord with experimental findings (cf.Table 2).
All these mechanisms relate to 'fast' processes/reactions, and their rates were monitored on the level of seconds (Figures 7 and 8).The 'slow' process of oxidation/reduction (X X ), is presented in Figure 10a.The slow dynamics of this process controls the evolution of pseudo-steady-state kinetic dependences for the formation of CO and CO 2 , and acrylic acid.
The mechanism presented in Figure 12 can be modified to take into account the peculiarities of acetic acid kinetic behaviour, i.e., the concentration goes through a maximum and then a slow decrease occurred.Presently, we do not have enough data to secure such interpretation.

Kinetic Equations
In this section, the kinetic equations related to our rational mechanism together with their corresponding apparent kinetic orders are presented route-by-route.As mentioned above, many routes are uncoupled, and can be easily derived based on theory of graphs [54].For all cases R is the rate of product formation and k i is the kinetic coefficient of corresponding reaction.[C 3 H 8 ] and [H 2 O] are the partial pressures of propane and water, respectively, and [O (b) ] is the concentration of oxygen in the bulk catalyst phase, which participates in the reaction.Thus, the foundation for our mechanism and the resulting kinetics is the well accepted Mars-van-Krevelen mechanism.n A,B is the apparent kinetic order of rate of substance A with respect to component B.
Rate of propene formation: Equation ( 19) becomes: The apparent kinetic orders for this cycle are: where θ Z and θ ZO are dimensionless reduced and oxidised surface intermediates in active site Z in the route of propene formation.Rate of acetic acid formation: with If terms containing [O (b) ] 2 will dominate over the other terms, and the rate expression reduces to where The apparent kinetic orders for this cycle are: The physical meaning of these characteristics is understandable.Typically, the apparent order regarding some substance, is equal to the dimensionless concentration (or sum of them), which react with this substance [54].As found for this cycle, the apparent order regarding oxygen, is equal to the sum of concentration of intermediates θ Y + θ YC3H6 , the apparent order regarding water is equal to θ Y O + θ YC2H4 , as for the apparent order regarding propane, it is equal to show, that it is equal to the sum of concentrations of θ YO + θ Y O .
Rate of acrylic acid formation: with If terms containing [O (b) ] 2 will dominate over the other terms, and the rate expression reduces to The apparent kinetic orders for this cycle are: It can be noted, that the kinetic expressions and apparent kinetic orders for both acids are similar, reflecting the similarity of both mechanisms on different centres.
Rate of CO and CO 2 formation: CO formation (only total oxidation) in Figure 12: CO 2 formation in Figure 12: with If terms containing [O (b) ] will dominate over the other terms, the rate expressions reduce to and The apparent kinetic orders for this cycle are: Table 3 summarises all apparent kinetic orders of the proposed mechanism in Figure 12.The obtained apparent kinetic orders allow to make a comparison with experimental values, presented previously.The non-trivial experimental fact is, that the sum of experimental kinetic orders is larger than one, 1.64 for acetic acid and 1.32 for acrylic acid.It is confirmed by our theoretical analysis.The theoretical values of these sums are for acetic and acrylic acid are 1 + θ Y O and 1 + θ XO , respectively.The introduction of a buffer step, increased the sum of kinetic orders on the value of concentration of buffer intermediate, which reacts with water.12.
Rate of 'slow' dynamics: Slow dynamics is represented by the differential equation: where The starting point is α = 1, when all the catalyst is reduced.At the final equilibrated state The solution of the differential equation is: where t is in the scale of time-on-stream.Or and Hypothetically, starting from the completely oxidised catalyst (β = 1), we will have and It is easy to show, that dividing Equation (51) by Equation ( 54) which is constant for a given temperature and partial pressure of oxygen.This apparent equilibrium constant for the oxidation-reduction process can be estimated as follows: For estimating α eq and β eq , let's present Equation (38) as where For a fixed temperature and gas composition R(CO 2 ) is constant.For a completely reduced catalyst at (t = 0; α = 1) Similarly, for the catalyst state at which we assume an oxidation-reduction equilibrium Roughly considering proportionality between measured CO 2 concentration and the rate of [CO 2 ]-change (inlet concentration of CO 2 = 0), we have where [CO 2 ] 0 and [CO 2 ] eq are the CO 2 concentrations at t = 0 and at equilibrium.Therefore, and Thus, β eq = 1 − α eq = 0.71, and 3. Discussion

Kinetics and Mechanism
The catalyst has been treated in reducing and oxidising reactant medium.Figures 9 and 10a report the catalyst evolution under oxidising atmosphere after a reductive treatment.The data suggest, that such reducing treatment of the M1 catalyst was reversible.These observations are in line with experimental findings on reduction and oxidation behaviour of MoVO x from Ueda et al. [55,56].In their contribution, two different types of lattice oxygen were distinguished.In early stages of treatments in reducing conditions type-1 oxygen evolves from the structure irreversibly, whereas type-2 oxygen leaves the structure continuously and returns reversibly in oxidising atmosphere [55].The structure of the M1 catalyst, however, is stable, even under extreme conditions [12].
We correspond to Ueda et al.'s [55] concept, reflecting it kinetically by our 'slow dynamics' in Figure 11.This mechanism visualises this reversible reduction-oxidation behaviour of the catalyst from a kinetic point of view, expressed through the kinetic equation for this slow dynamics (cf.Equation ( 17)).The apparent equilibrium constant K app in Equations ( 55) and ( 56) may therefore be correlated with the degree of reduction (parameter δ in [55]).The alkane conversion was significantly affected by the degree of reduction [55], which confirms the importance of oxygen (and water) in providing a certain chemical potential in the gas phase to stabilise the catalyst surface structure [2,12].According to Barrault et al. [3], and supporting the statements above, a high selectivity towards selective oxidation is only provided by an oxidised site (not a reduced site).This is expressed in our proposed bridge mechanism in Figure 11.
The Mars-van-Krevelen (MvK) mechanism is widely accepted for selective oxidation reactions over mixed metal oxide catalysts [1,2,36,37].Further significant inspiration to our mechanistic development was provided by different concepts: (i) the site isolation concept by Grasselli [19][20][21]; (ii) the concept of an ensemble of atoms forming an active site by Schlögl [22]; (iii) the proposed multi-route parallel reaction pathways by Schlögl et al. [8] and Barrault et al. [3]; and (iv) the change in product distribution upon surface acidity by Barrault et al. [3] and Luo et al. [46].
Based on these concepts, the experimental observations in Figures 7, 8, and 10a were interpreted.Our mechanism was presented in Figure 12 and consists of four independent catalytic cycles for the formation of propene, acetic acid, acrylic acid, and total oxidation, respectively.All products (except total oxidation) originate from a specific catalytic cycle, a specific active centre, which is in line with Schlögl [22] and Grasselli [19][20][21].Acetic acid and acrylic acid are formed via different reaction pathways [3].The total oxidation products CO and CO 2 are proposed to as well originate each from an individual specific active centre.This was also earlier proposed by Naumann d'Alnoncourt et al. [12].The data presented in this work is, however, not sufficient to support this hypothesis.Therefore both products are summarised in one cycle.Such a separate reaction pathway from propane to total oxidation products was earlier proposed [8,46], together with the stability of other reaction products against further oxidation under the given reaction condition [8].In most mechanisms, however, the total oxidation is considered to be the final product of all species involved (e.g., [6,9,10]).
Water has a special role in many mechanisms presented in the past.In our mechanism, water is both a reactant in the selective oxidation of propane, and at the same time an activator for both catalytic cycles forming acetic acid and acrylic acid.For example Centi et al. [34] and Novakova et al. [6] confirm the necessity of water for the activation of the surface and to increase the selectivity to selective oxidation products, respectively.Schlögl et al. confirm these observations by in situ XPS investigation [25,32].
Propene has been chosen as a product during the selective oxidation of propane, which re-adsorbs for further oxidation, e.g., [3,6,9,10,12,46].Certainly, if the data in Figure 8 would be plotted as selectivity against conversion or feed concentration, the propene-selectivity would drop (since the selectivity is affected by the concentration of all products).This would easily lead to the conclusion that propene re-adsorbs to be further converted to acrylic acid.Due to the observations in Figures 8  and 10a, i.e., a constant propene concentration during increasing acrylic acid concentration (and, thus, increasing propane conversion), propene is in our mechanism (cf. Figure 12) kinetically independent.The tools for the mechanistic development were provided by the theory of graphs [54].
Our apparent kinetic orders were presented in Table 2 are in good agreement with other investigations.Apparent kinetic orders are high with respect to propane for propene and intermediate high for acetic and acrylic acid.All products show an apparent kinetic order close to zero with respect to oxygen (which could be expected for a MvK-type reaction) [6][7][8].Kinetic orders with respect to water are hardly reported.

Structure-Kinetic Interrelationship
Based on the results, interpretations and discussions in this work, an overall model of our kinetic understanding for catalytic surface reactions is provided in Figure 13.It incorporates the concepts of isolated active sites presented by Grasselli [19][20][21] and the concept of an ensemble of atoms forming an active site presented by Schlögl [22].
Figure 13 depicts our understanding of the combination of different 'active sites' on a catalyst surface.An active site can be characterised by a complex structural architecture.An 'active site' consists of one or more 'active centres'.There are possibly as many active centres as there are kinetically independent products of a reaction.Thus, if a reactant splits on one active site into two product molecules, they are not kinetically independent.Therefore, an active site is in our understanding characterised by a 'kinetic fingerprint' of the product(s) evolving from it.
Furthermore, in Figure 13 are two levels of isolation expressed: (i) Two active sites are separated on a slow and on a fast time scale.(ii) Two active centres (on the same site) are separated on a fast time scale, and joint together by a slow process.
The final hypothesis in our understanding is: Each kinetically independent product originates from a specific active centre.The chemical nature of an active centre determines uniquely how and to which product the reactant molecule(s) are converted.Manifestation of the isolated sites showed in the previous analysis, using the partial balance of carbon atoms [CO] + [CO 2 ] + [acrylic acid], which was found experimentally.The complex nature of the isolated active site, was shown based on the difference long-term time-on-stream dependences for CO and CO 2 from one side and acrylic acid from another side, respectively.We assumed, that the 'fast' processes of product formation are linked by the 'slow' process of centre oxidation-reduction.For acetic acid and propene formation, we also assume different centres, considering the differences in long-term kinetic trends.
As for distinguishing multi-route and single-route scenarios, there is a potential to answer this question, based on selectivity analysis.
An analysis of the selectivity, rather to say, the differential selectivity, in particular the ratio of rates of different products or the ratio of product rate over the reactant rate can be very fruitful.As known, the typical kinetic rate expression has the form of Langmuir-Hinshelwood-Hougen-Watson (LHHW) expression, and the form of this expression is 'a numerator over denominator'-equation.For the multi-route reaction on the same centre, kinetic equations for all routes have the same denominator.In this case, the differential selectivity will be presented by the ratio 'of numerators' which can be a quite simple function.For the multi-route reactions on different centres, denominators of kinetic expressions (LHHW-equations) will be not the same, and selectivities will be functionally not so simple.Therefore, the observed partial balances can be considered as a potential fingerprint of the multi-route process over the isolated site.
Comparing rate of CO and CO 2 formation, we conclude: the corresponding route rates (according to Temkin-Horiuti terminology [29]) are equal, because the routes include the same step of interaction between the centre X and propane.Rate of oxygen consumption and water production for these routes will be different, according to the stoichiometric equations.
Analysing kinetic equations for different products, and assuming oxygen excess, we can consider that steps, which do not contain oxygen are limiting in the steady-state process.These steps may be assumed slow in comparison with oxygen-containing steps of the same route.However, they are certainly 'fast' in comparison with the 'slow' oxidation-reduction step.
Generally, if the reaction is considered to occur on different sites via a single route, the apparent kinetic orders and apparent energy of activation can be used as a traditional tool for estimating surface intermediates and energetic parameters of elementary steps.We presented analytic expressions for the apparent orders of all five products with respect to all three reactants, i.e., for propane, oxygen, and water (cf.Table 3).

Synthesis and Characterisation
A phase pure M1 MoVTeNbO x catalyst (internal identification number 17817) was synthesised according to Kolen'ko et al. [57], following a slurry-spray-drying protocol.The as-prepared catalyst was characterised by powder X-ray diffractometry (XRD), nitrogen adsorption (Brunauer-Emmett-Teller method-BET), Raman, and electron microscopy (i.e., SEM and TEM) methods.A second series of TEM and XRD analysis was carried out for the used catalyst, after (at least) all experiments described in the present work.

Reactor Set-Up and Reaction Conditions
The kinetic experiments were carried out in one reactor of a 10-fold parallel reactor unit, which has been presented in more detail before [50].A gas chromatograph (7890A, Agilent, Santa Clara, California, United States) was used to perform online gas analysis.It was equipped with a combination of PoraplotQ and Plot Molesieve columns connected to a thermal conductivity detector, and a combination of HP-FFAP and PoraplotQ columns connected to a flame ionisation detector.All reactant and product gases were analysed (including water).100 mg of catalyst were placed over 50 mg of SiC, both in a sieve fraction of 250-355 µm.The catalyst was pressed (2.5 t) and sieved to receive the desired sieve fraction.
All experiments were carried out at 573 K, 1 bar pressure and a weight hourly space velocity (WHSV) of 10 4 mL g −1 h −1 .A typical reactant composition was C 3 H 8 /O 2 /H 2 O/N 2 = 3/12/10/75.At partial pressure variation experiments, one reactant partial pressure was changed.The flow of nitrogen was adjusted accordingly to keep the WHSV constant.The propane conversion was kept below 10%.
The formation rate of products is determined as follows: where R j is the formation rate of product j in mol g −1 s −1 , v the molar total flow rate in mol s −1 , m cat the catalyst mass in g, and [vol − %] j the determined outlet concentration of product j by the gas chromatograph divided by 100%.
• Both experiments were performed in a tube reactor (the plug flow reactor, PFR).In both experiments, the steady-state kinetic data have been produced.Obviously, there is a big difference between the two experiments.
The first one corresponds to the 'fast' and 'moderate' processes of gas transformation over the solid catalyst at the given catalyst state (catalyst composition).The second one relates to the 'slow' process within the 'gas-solid' catalytic system.Such 'slow' kinetic behaviour was observed in studies over the preliminary reduced catalyst.In experiments over oxide catalysts, which have been described in the literature previously, such 'slow' trends were not observed.Kinetic evolution demonstrated in our experiments can be explained by gradual oxidation of the catalyst under the influence of the reactive medium.
All experiments, both in with a scale of seconds and with a scale of minute/hours, were reproducible and undermine its validity.

Conclusions
The selective oxidation of propane to acrylic acid was kinetically investigated over a phase pure M1 MoVTeNbO x catalyst at 573 K and 1 bar pressure in a plug flow reactor.
The investigation was performed on two time scales, i.e., 'slow' (minutes and hours) and 'fast' (seconds).Characterisation by XRD showed no structural changes of the catalyst, despite repeated oxidation and reduction treatments.Observed differences in long term behaviour, have been used for distinguishing catalytic routes for the formation of propene, acetic acid, acrylic acid, CO, and CO 2 .Based on the partial carbon balance ([CO] + [CO 2 ] + [aa] = const.),the existence of an isolated site was assumed for the production of these substances.Two different centres are assumed within this site for the production of CO and CO 2 and for the production of acrylic acid, respectively.The relationship between these routes is linked by a slow step of catalyst oxidation.
The partial carbon balance, established here, can be considered as a 'kinetic fingerprint', which is a special combination of products originating from one active site, for the consideration of the known concept of isolated sites.A 'kinetic fingerprint' distinguishes types of active sites.An active site consists of one or more active centres.Therefore, each kinetically independent product originates from a specific active centre.The chemical nature of an active centre determines uniquely how and to which product the reactant molecule(s) are converted.
In the development of the rational mechanism, the principle was used, 'every specific route is distinguished by a specific centre'.Kinetic rate expressions have been found for all routes and for a slow oxidation-reduction process.
The main message of our article is: the kinetic investigations of 'slow dynamics' in the scale of min and hours reveal long term trends in catalyst evolution.In combination with investigations over a catalyst at steady-state, such experiments give an insight in mechanistic aspects.The performed research can be considered as a development of methodology for revealing a rational mechanism of complex catalytic reactions.
It is interesting to mention, at fixed external parameters (i.e., flow rate, concentration, and temperature) the selectivity of highly molecular substances (i.e., acrylic acid) is growing during the slow re-oxidation of the catalyst, exhibiting a simple example of chemical evolution.

Figure 1 .
Figure 1.Scheme of three levels of kinetic complexity.

Figure 3 .
Figure 3. Reaction schemes taken from Ref. [46].Scheme 4 (left): "Reaction network for propane oxidation over Mo 1 V 0.3 Te 0.23 Nb 0.12 O x .Shaded arrow refers to main reaction pathway; compounds with underlines refer to the main products involved in propane and propene oxidation reaction."and Scheme 5 (right): "Reaction network for propane oxidation over VO-H beta.Shaded arrow refers to main reaction pathway; compounds with underlines refer to the main products involved in propane and propene oxidation reaction".

Figure 6 .
Figure 6.XRD diffractogram of the fresh and used sample -typical reflexes for M1 structure appear, the bulk structure of the catalyst is preserved.Each diffractogram was normalised with respect to its corresponding maximum value.

2. 4 .
Kinetic Experiments under Catalyst Steady-State Conditions Partial pressure variation experiments were performed under steady-state conditions (C 3 H 8 /O 2 /H 2 O/N 2 = 3/12/10/75, 573 K and a weight hourly space velocity (WHSV) of 10 4 mL g −1 h −1 ).The feed of one reactant (i.e., C 3 H 8 , O 2 , and H 2 O, respectively) was changed, maintaining the feed of the respective other reactants and adjusting the flow of nitrogen in order to keep the WHSV constant.The observations are shown in Figures

Figure 12 .
Figure 12.Summary of the proposed reaction mechanism of selective oxidation of propane.

Figure 13 .
Figure 13.The data interpretation lead to the following model: (a) A catalyst has different types of 'active sites' (here, 1 and 2).A specific product (mixture) is formed over each active site.Active sites have one or more 'active centres'.(b) The evolution of the active centres on site 1 is linked together (one appears, the other disappears), thus, the sum of product evolution as a function of time over site 1 is constant.This phenomenon is called a 'kinetic fingerprint'.
(a) after 1 h in reducing condition (b) after 3 h in reducing condition

Table 1 .
Experiments: catalysts and working conditions.

Table 2 .
Apparent kinetic orders of product formation rates presented in Figures7 and 8as a function of the three reactants.

Table 3 .
Apparent kinetic orders with respect to the three reactants ([C 3 H 8 ], [O 2 ], and [H 2 O], respectively) of the mechanism proposed in Figure

Kinetic experiments during catalyst evolution ('slow dynamics' scale of min and hours):
Observation of catalyst evolution with time-on-stream (TOS) under C 3 H 8 /O 2 /H 2 O/N 2 = 3/12/10/75 feed, after treatment in reducing conditions for 12 h in C 3 H 8 /O 2 /H 2 O/N 2 = 3/0/10/87.Both parts of the experiment were carried out at 573 K and a constant WHSV of 10 4 mL g −1 h −1 .