y of chemical reaction equilibria in pores of activated carbons †

Stanisław Staszic University of Applied S 64-620 Piła, Poland. E-mail: sfurmaniak@p Physicochemistry of Carbon Materials R Nicolaus Copernicus University in Toruń, G School of Engineering and Information Te 6150, WA, Australia Department for the Modelling of Physico-C Maria Curie Skłodowska University in Lub Poland † Electronic supplementary informa 10.1039/c7ra08992a Cite this: RSC Adv., 2017, 7, 53667


Introduction
In 1994, Johnson et al. 1 and, independently, Smith and Triska 2 published a description of methodology allowing the modelling of chemical reactions with the help of the Monte Carlo simulation method.This technique (called the reaction ensemble Monte Carlo, RxMC) predicts the locations of equilibrium states according to the laws of statistical physics.In order to use the technique, one needs to know the ideal-gas free energies (including the intramolecular contributions due to vibrational, rotational and electronic degrees of freedom) for the reacting species and the potentials describing the interactions in the system (adequately dened force eld).Since the chemical reactions are modelled by some additional "reaction" steps (a simultaneous deletion of the reactant molecules and insertion of the product molecules according to the reaction stoichiometry), the RxMC simulations can be performed under different conditions, like isobaric or isochoric, and in different statistical ensembles. 3The RxMC simulation method and its various applications have been very well reviewed in the papers by Turner et al. 4 and by Dubbeldam et al. 3 One should note that the RxMC simulation technique has been successfully used to study reactions at high temperatures and/or under high pressure, reactions in solutions, at the interfaces between different phases and, what is of particular interest here, the equilibrium states of reactions under nonideal conditions occurring in porous solids or near solid surfaces.6][7][8] One should note that in the case of reactions occurring in pores, the accurate experimental determinations of equilibrium constants may be extremely challenging (if possible at all).This is caused by, among other things, the difficulties in accurate measurements of concentrations of the reacting species in pores. 4,9,10orówko et al. 5 used the RxMC method to simulate the reaction process on the model surface (being either a noninteracting hard wall and an attractive wall).The molecules of the reagents were modelled as non-interacting hard-spheres.Despite the fact that this pioneering work used some rather articial systems (lacking the attractive interactions between the reacting molecules), it conrmed the usefulness of the method for modelling the inuence of external elds (the adsorbing wall) on the chemical equilibrium.In the next paper, Borówko and Zagórski 6 extended their considerations to ideal slit-like pores.They studied the equilibrium states of chemical reactions, taking into account both the interactions between reactant or product and the pore walls, as well as the interactions between the reacting molecules.Their results gave clear evidence of the signicant impact of an external eld generated by the pore walls on chemical equilibria, i.e., the changes (an increase or decrease) of the yield of reaction in the pore in comparison with the bulk system.The studies conducted over the last two decades have involved some more complex porous systems.For example, Turner et al. 11,12 studied the ammonia synthesis in the pores being replicas of some real carbons obtained from the reverse Monte Carlo method.The comparison between the ideal slit-like pores and more realistic models of carbon demonstrated that the corrugation of pore walls and the pore connectivity did not affect the yield of ammonia synthesis much.This conclusion was the consequence of a relatively high simulation temperature, because a lower reaction temperature could possibly magnify the effects of the geometric pore characteristics.7][8][11][12][13][14][15][16][17][18][19][20][21][22][23] The results of RxMC simulations of reactions in porous systems published hitherto may be grouped into the following classes involving different pore models: (a) ideal pores having a well-dened geometry (innite slit-like and cylindrical pores/nanotubes), [4][5][6][7][8][12][13][14][15][16][17]22,23 (b) inorganic ordered adsorbents, i.e. silica, 17,22 silicalite-1, 18,19 pillared clays, 20 MCM-41, 20 and zeolites (MFI, TON, LTL and FER) 19 or (c) replicas of real carbons obtained from the Reverse Monte Carlo method (RMC). 8,11,21 On he other hand, different processes were simulated: (a) simple model reactions involving LJ-type interaction potentials: isomerisation reaction A 4 B, 15 dimerisation/ association reaction A + A 4 B (with the formation of a spherical product) 5,16,24 or A + A 4 A 2 (with the formation of a nonspherical linear product) 5 and (b) real reactions: ammonia synthesis, 8,[11][12][13]20 nitric oxide dimerisation, 8,13,14,21 the hydrogen iodide decomposition, 7,8,12 the propene metathesis, 18,19 and CO 2 methanation, 17,22 and xylene isomerisation.23 It should be pointed out that the pore chemistry and the pore morphology can also be investigated.22 To summarize, the results of the above-described selected works (on the modelling of chemical reactions using RxMC method) allow a better understanding of the factors affecting the reaction equilibrium in pores.However, many issues remain unsolved. In pat, this is connected to the method of modelling the pore structure. On on hand, the isolated innite ideal slit-like pores are only a rough approximation of the real carbons. Th][26][27] On the other hand, the replicas of real carbons (although they are more realistic and reect the properties of such materials rather well) do not allow the modications to their properties, e.g.porosity. Hece, a systematic study aiming at the evaluation of the inuence of porosity changes on the reaction equilibrium may be impossible to conduct with the use of the "inexible" replicas.In order to perform such a study, other virtual porous carbons (VPCs) models seem to be more useful.28,29 The use of VPCs makes it possible to generate a series of materials with a gradually changing structure. However according to our knowledge, there are no reports based on the series of realistic model carbons with systematically modied porosities.Typically, such studies are based only on one or few replicas of particular real carbons.8,11,21 The study described in this paper is an attempt to ll the gap.
We have considered here different simple model reactions, such as the isomerisation (A 4 B), the dimerisation (A + A 4 B) and the synthesis (A + B 4 C), on the series of VPCs of systematically changing porosity.We have not strictly dened the chemical nature of the reacting molecules (they have been modelled as the single Lennard-Jones centres).This facilitates a systematic study allowing to determine qualitative relationships between the porous structure of carbons and the reaction equilibria, for various combinations of the energetic parameters of the reacting molecules and/or the different values of the equilibrium constant.

Simulation boxes
We have used the series of nine VPCs described in detail in the earlier work. 30The model carbons (shown schematically in Fig. 1a) have been generated using the simple Metropolis Monte Carlo method, and one of the most sophisticated carbon force elds, i.e. the environment-dependent interaction potential for carbon (EDIP) proposed by Marks. 31,32All the structures have been placed in the cubic simulation box (4.5 Â 4.5 Â 4.5 nm) with periodic boundary conditions in all three directions.The subsequent VPCs in the series differ by the density of carbon atoms.The following values of the density have been considered 0.5, 0.6, 0.7, 0.8, 0.9, 1.0, 1.1, 1.2, and 1.3 g cm À3 .Each structure is denoted as dx.xwhere x.x is the carbon density.
The porosity of the carbonaceous adsorbents have been characterized by a geometrical method proposed by Bhattacharya and Gubbins (BG). 33The implementation of the method was described in details elsewhere. 30,34The BG method makes it possible to determine the histograms of effective pore diameters (d eff ).The data may be also used to calculate the average diameter of pores accessible for the molecules having a given diameter (i.e. the pores wider than the assumed minimum diameter (d min )): where P(d eff ) is a probability of the pores having the effective diameter equal to d eff .Besides, the combination of the BG method and Monte Carlo integration 30,35 has been used to determine the accessible volume within the pores (V acc (d min )).

Chemical equilibria modelling
The majority of the studies 5-8,11-16,18-21 using the RxMC method to model reaction equilibria in pores described in literature may be divided into two groups.Here, the simulation system is composed of two subsystems.One of them is the box containing the adsorbent, while the other box contains only the bulk gaseous phase.Since the transfer of molecules between the subsystems is allowed and the rst box has constant volume, the size of the bulk system is modied during the simulations in order to keep the constant pressure.This method is closer to the real experimental situation, in which the reactants are usually adsorbed from the gaseous phase.The use of Grand Canonical Monte Carlo technique (GCMC) may be an alternative for the direct modelling of the bulk reservoir.Hansen et al. 18 found that the GCMC simulations are superior to the Gibbs ensemble simulations for reactions where the bulk-phase equilibrium can be calculated in advance, and does not have to be simulated simultaneously with the molecules inside the pore.In the case of such simulations, the assumption of chemical equilibrium in the bulk phase determines the values of chemical potential of the reacting species.In the GCMC simulations these chemical potentials are the same also in the pores.So, the conditions of chemical equilibrium are fullled also in the adsorbed phase.This state may be achieved only with the use of GCMC trial moves (creation and annihilation of the molecules of all the reacting species).The additionally performed attempts of reaction steps (especially when they are occurring with a low probability) may not give an important contribution to the determination of equilibrium composition.
The same conclusions may be derived also from thermodynamic consideration.Chemical potentials of the mixtures components in two different phases are equal if the phases are in equilibrium.One can assume that there is chemical equilibrium in one of the phases (the sum of the products of the stoichiometric coefficients of all the reacting species and their chemical potentials are equal to zero, i.e.Sm i n i ¼ 0). 1,2,4In such a situation (equality of m i in both phases) the condition of chemical equilibrium is fullled also in the other phase even if the reaction in this phase does not occur.The fact that reaction steps are not necessary in such simulations is indirectly also conrmed by the results obtained by Lísal et al. 21They showed that the composition of the mixture reacting in pores may be predicted with the use of single component adsorption isotherms and ideal adsorbed solution model.Since our aim is to directly simulate reactions in pores, including adsorption equilibrium with bulk phase, we have modied the scheme of GCMC simulations with the chemical reaction proposed by Hansen et al. 18 In order to make the reaction steps responsible for determination of the adsorbed mixture composition (mainly the amount of the product molecules), we have assumed that the reaction can occur only in the pores (i.e. the adsorbent may be formally treated as a catalyst).Thus, the bulk phase contains only the reactants (the desorption of the product molecules has not been directly considered).In this case, the modelling of gaseous reservoir is not necessary.The adsorption of the reactants may be modelled with the use of the trial moves utilized during the Grand Canonical Monte Carlo simulations (GCMC), i.e. the creation of molecules at randomly selected positions and the annihilation of randomly chosen molecules in the simulation box. 36The other trial moves that have been used to modify the state of the system involved the random displacements of the also randomly chosen molecules and the (forward and backward) reactions steps, which are responsible for the creation and destruction of the product molecules in the pores.
The equation representing the chemical reaction may be written in the general form as: where S is the number of the reacting species (reactants and products), n i refers to the stoichiometric coefficients (negative for the reactants and positive for the products), and R i refers to the formulas of the species.The simplest reaction step may be realized by the destruction of a certain numbers (dened by the corresponding values of n i 's) of the randomly chosen reactant molecules (forward step) or product molecules (backward step) and the random creation of the equivalent numbers of product or reactants.The acceptation probability of such a trial move is dened as: 1,2,4 where V is volume of the simulation box, N i is the number of molecules of the i-th species, b ¼ 1/k B T (k B -Boltzmann constant, Ttemperature), DU is the change in the congurational energy, and nis the sum of stoichiometric coefficients: The parameter G (the reaction quotient) is related to the bulk equilibrium constant, and is dened as: where C i,bulk,eq refers to the equilibrium concentrations (in number of molecules per volume) of all the species.The value of reaction quotient is equal to the equilibrium constant in the case of isomerisation reaction.In the case of other two kinds of considered reactions (i.e.dimerisation and synthesis) these two quantities are mutually proportional. 4One should also note that the constant G may have some units (if ns 0).The optional sign (i.e.'AE') in the eqn (3) is positive ('+') for the forward reaction step and negative ('À') for the backward one.
We have used the above-described approach to model chemical equilibria in the pores of all the VPCs studied (cf.Fig. 1 and Table 1).Three model reactions have been considered: (a) isomerisation: (b) dimerisation: and (c) synthesis: In order to simplify the calculations, we have assumed spherical symmetry for all the reacting molecules (A, B and C).The physical interaction between all the species involved has been assumed as described by the Lennard-Jones (LJ) potential.The cross-interaction parameters have been calculated using the usual Lorentz-Berthelot mixing rules. 36It should be noted that for simplicity, during the discussion in paragraph 3, we have given only the values of the potential well depth of interactions between the molecules of the same type (i.e. 3 ii ).However, the analysed changes are also the consequence of differences in energy of interactions with other species including adsorbent atoms and here the assumed mixing rules are also crucial.For the carbon atoms, we have used the following values of LJ parameters: s SS ¼ 0.34 nm and 3 SS /k B ¼ 28 K. 37 The interactions between all the pairs have been cut at the distance equal to 5 Â s ij .We have not applied the long range tail correction since the assumed conditions are far from any points of phase transitions.Our study has had qualitative character and this numerically costly procedure was not necessary.We have also xed the LJ parameters for the A molecules by taking s AA ¼ 0.35 nm and 3 AA /k B ¼ 100 K (in all the cases).Besides, in the case of synthesis reaction, we have assumed that the collision diameter for molecules of reactant B is the same as for A, i.e. s BB ¼ s AA ¼ 0.35 nm.The collision parameters of the product molecules have been calculated assuming that the volume of the product molecule is equal to the sum of the volumes of reactant molecules.Only in the isomerisation reaction the volumes of the reactants and products are the same.This arbitrary assumption reduces the number of parameters whose inuence should be tested.Analogical simplications were also applied by others. 6hus, in the case of the isomerisation reaction, s BB ¼ s AA , for the dimerisation and the synthesis reactions, we have taken ) 1/3 , respectively.In consequence, the product collision diameter for the isomerisation reaction was equal to 0.3500 nm and for the dimerisation and the synthesis reactions it was equal to 0.4410 nm.We are aware, of course, that in the case of dimerisation and synthesis reactions in the narrow pores, the assumption of spherical symmetry for the products of reactions is a rather crude approximation.The non-sphericity of the products may affect the results considerably, and this problem will be addressed in the future papers.
All the reaction equilibria have been studied at the xed temperature of T ¼ 298 K.Each simulation (for the chosen set of parameters: the equilibrium constant (G), the LJ parameters, and the pressure of the reactant(s) in the gaseous phase) has consisted of 1 Â 10 7 cycles.During one cycle, 5000 iterations have been preformed.Each iteration has been a single attempt to change the system state via a randomly chosen trial move (the creation or the annihilation of the reactant molecule and displacement of randomly chosen molecule (both the reactant as well as the product)).In addition, during randomly chosen iterations the (forward or backward) reaction step has been performed (with probability equal to 1/500).During the reaction steps, the position of created molecule (one of the molecules) in half of the randomly chosen cases was the same as the position of removed molecule (one of the molecules).In other cases, new molecules have been created in the randomly chosen position in the box.Table S1 in the ESI † collects the acceptation criteria of different trial moves.The rst 4 Â 10 6 cycles (2 Â 10 10 GCMC trial moves and ca. 4 Â 10 7 attempts of the reaction step) have been treated as equilibration and the data collected during the next 6 Â 10 6 cycles (3 Â 10 10 GCMC trial moves and ca.6 Â 10 7 attempts of the reaction step) have been used to calculate averages.In the case of isomerisation and dimerisation reactions, we have assumed the pressure of the reactant (A) in the gaseous phase to be equal to 0.1 MPa.In the case of the synthesis reaction, we have assumed the total pressure of the reactants in the gaseous phase (i.e.p A + p B ) to be equal to 0.1 MPa.For all the considered reactions and the combination of parameters, we have also simulated the reaction equilibria in the bulk phase.The methodology of these computations has been exactly the same as described above (i.e. the reactant(s) molecules have been created and destructed according to the GCMC trial moves and the product molecules have been introduced/removed into/from the simulation box via the reaction steps in the absence of adsorbent), but we have used the box of the size 20 Â 20 Â 20 nm without the adsorbent.The results for the bulk phase have been treated as reference.

Characteristics of VPCs
In Table 1 we have summarized the parameters characterizing the VPCs considered, and the histograms of their pore diameters have been shown in Fig. 1b.It is evident that all the structures are strictly microporous, i.e. the pore sizes do not exceed 2 nm.The lack of wider pores results from the relatively small sizes of the model structures.27]38 One should note that the synthesis methods leading to highly microporous carbons have also been reported, and these materials seem to be very promising adsorbents. 25,27rom the results given in Fig. 1b it follows that the porosity of the VPCs considered here changes systematically, i.e. passing from the d0.5 carbon to the d1.3 one, the contributions vanish gradually due to the wider pores.The last two VPCs (i.e.d1.2 and d1.3) have only quite small pores, about 1 nm in diameter.In consequence, the average pore diameter decreases as well (Table 1).The only exception is the d1.1 structure, which contains a small amount of wider pores (ca.1.6 nm) and its average pore diameter deviates slightly from the general trend.

Isomerisation reaction
We have begun the discussion of the chemical reactions in VPCs by considering the isomerisation reaction.Fig. S1 in the ESI † compares the densities of the reacting molecules in pores and the mole fractions of the product in different VPCs, obtained for different values of the equilibrium constant (G).The assumed values of G correspond to the following ratios of average numbers of A and B molecules for reaction in the bulk phase: 4 : 1, 2 : 1, 1 : 1, 1 : 2 and 1 : 4. The simulations have been carried out assuming that the potential well depth for the product molecules (3 BB /k B ¼ 120 K) is slightly higher than for the reactant molecules (3 AA /k B ¼ 100 K).The concentrations of the reactant and product molecules increase with passing from the carbon d0.5 to d1.3, regardless of the value of equilibrium constant, since a gradual decrease of the pore diameters causes the rise in adsorption.In general, an increase of the equilibrium constant is expected to cause a rise in the amount of the product in the system, and this is illustrated by the results given in Fig. S1 in the ESI.† The carbon porosity affects the densities of The pores accessible to the reactants and also to the products of isomerisation reaction.
c The pores accessible to the products of dimerisation and synthesis reactions.
reacting molecules, but it does not affect appreciably the reaction yield.The analysis of the product mole fractions reveals that in all the systems the reacting mixture in pores is richer in B molecules, than for the same reaction taking place in the bulk phase.Moreover, the effect of connement on the reaction yield is gradually enhanced when the average pore size decreases.
The observed behaviour results from the differences in the energetic parameters of the reacting molecules (higher 3 value for the product).Thus, the rise in the carbon density and/or the reduction in the pore diameters enhances the adsorption of the product molecules (B).Besides, the connement and the resulting adsorption of the species result in their higher concentration, as compared with the bulk.
In order to check the effects of the differences in the adsorption energetics of reacting molecules, we have performed the simulation for different values of the potential well depth of the product, assuming that G ¼ 1, i.e. for the A to B ratio equal to 1 : 1, when the reaction takes place in the bulk phase.The obtained results are shown in Fig. S2 in the ESI.† As expected, when 3 BB ¼ 3 AA (3 AA /k B ¼ 100 K) the average concentrations of the reactant and the product molecules are the same.The differences in adsorption energetics of the reacting species cause that the preferentially adsorbed compound (A or B) dominates (it should be noted that the number of (A) molecules is mainly determined by its pressure in the bulk phase which is the same for all the systems).Of course, these energetic parameters hardly affect the composition of the bulk reacting mixture.Since the densities of the reacting species are low, intermolecular interactions have a negligible impact.When the product molecules (B) are preferentially adsorbed, the reaction yield in pores is higher than in the bulk phase.On the other hand, when 3 BB < 3 AA, , an opposite effect appears, and the reaction yield in pores is lower than in the bulk.A comparison of the data for different VPCs clearly reveals some regularities.The reduction in the pore sizes and/or the rise in the carbon density cause the product mole fraction to increase (decrease), when the product adsorption energy becomes higher (lower).

Dimerisation reaction
Before we start the discussion of results for the dimerisation reaction, one important fact should be noted.In the case of interactions with complex structures (for example carbonaceous), the energy of adsorbed molecules depends not only on the potential well depth but also on the collision diameter as, for example, in the analytical Steele potential (interactions with ideal at walls). 37However, this dependence is not clearly dened in the case of disordered structure.For the dimerisation (and also synthesis reactions) we have assumed that the size (collision diameter) of the product molecules (B) is bigger than of the reactant molecules (A).Therefore, it is very hard to nd the value of 3 for (B) that would result in the same energetics of adsorption of the A and B species.However, we have found the base value (3 BB /k B ¼ 60 K) which results in the similar energetics of solid-uid interactions.
Fig. S3 in the ESI † presents the composition of the reacting phase in the dimerisation reaction occurring in our porous materials for different values of the equilibrium constant (in this case we have assumed that 3 BB is 20% higher than the base value).The assumed values of G correspond to the following ratios of average numbers of A and B molecules for reaction in the bulk phase: 4 : 1, 2 : 1, 1 : 1, 1 : 2 and 1 : 4. The observed regularities are very similar to those found for the isomerisation reaction (cf.Fig. S2 in the ESI †).However, the average densities of the product molecules (B) are lower, since the volume occupied by a single (B) molecule is now higher.The assumed value of the potential well depth for the product molecules results in their higher adsorption, and hence in the higher reaction yield than in the gaseous phase.Regardless of the value of the equilibrium constant, the mole fractions of the product increases up to the d1.0 carbon, while it is nearly constant, or decreases slightly when G is sufficiently high, for the VPCs of high concentration of micropores.A careful analysis of the results reveals also a small anomaly observed for the carbon d1.2, which points to a slightly lower reaction yield than the adjacent structures (d1.1 and d1.3).Fig. 2 presents the product mole fraction in different porous carbons obtained for the systems characterized by different values of potential well depth for B molecules.The middle value of 3 BB /k B ¼ 60 K has been chosen is such a way that the reaction yields for all the carbons have been similar to the values corresponding to the reaction in the bulk.This value of 3 BB corresponds to a similar adsorption energetics of (A and B) molecules.Similarly to the case of the isomerisation reaction, the values of 3 BB higher (lower) than 60 K, lead to the higher (lower) product mole fraction than in the bulk phase.The reduction of pore sizes and/or the rise of the carbon density signicantly changes the product mole fraction.This effect is strong for the carbons up to d1.0, and becomes rather small for the carbons d1.1, d1.2 and 1.3.A small negative anomaly for the carbon d1.2 is observed again (with the exception of the case with the lowest value of 3 BB /k B ¼ 36 K).

Synthesis reaction
In the case of synthesis reaction, the situation becomes more complex.The system contains two kinds of reactants (A and B) molecules in the bulk phase, and three kinds of particles (the reactants A and B, and the product C), and there are more parameters that have to be taken into account.The rst series of the calculations performed has concerned the effects due to the value of equilibrium constant on the reaction yield.Therefore, we have kept all other parameters xed.We have also assumed equimolar composition of the bulk phase, i.e. p A ¼ p B .However, the reactants A and B have been assumed to be characterized by different values of the potential well depth (3 AA /k B ¼ 100 K and 3 BB /k B ¼ 110 K).In consequence, the component B adsorbs more strongly and its density in the system is higher.For the assumed value of 3 CC /k B ¼ 72 K for the product C, its concentrations in porous materials are larger than in the gaseous phase (the used values of G correspond to the following ratios of the numbers of molecules A : B : C: 4 : 4 : 1, 2 : 2 : 1, 1 : 1 : 1, 1 : 1 : 2 and 1 : 1 : 4 in the bulk).The results obtained have been shown in Fig. S4 in the ESI.† Of course, the rise in the value of G increases the reaction yield.The higher concentration of C molecules hardly affects the concentrations of the reactants (A and B) for the carbons from d0.5 up to the d0.8, but leads to a decrease of their concentrations in highly microporous VPCs.The differences between the product mole fractions in pores and in the gaseous phase for all the values of G are always positive, and again a small anomaly for the carbon d1.2 is observed, i.e. the product (the reactants) mole fraction of product is slightly lower (are slightly higher) than in the case of adjacent VPCs (d1.1 and d1.3).
We have also studied the effects due to changes in the potential well depths of the product molecules (3 CC ) on the reaction equilibrium, and the results have been presented in Fig. 3.The assumed value of G corresponds to the A : B : C ratio in the bulk phase equal to 1 : 1 : 1 and the values of the other parameters have been kept the same as before.The observed changes are qualitatively consistent with the previously discussed effects of the changes of the energetic parameter of the product molecules for the dimerisation reaction (Fig. 2).Thus,  the middle value of 3 CC (60 K) results in similar mole fractions of the product in pores and in the bulk phase.The higher values of 3 CC lead to higher reaction yields than those obtained for the reaction in the gas phase, and the lower values of 3 CC result in the lower values of x C .The increase as well as the decrease of the reaction yield becomes higher when average pore size decreases up to the carbon d1.0.The last three VPCs in the series (d1.1-d1.3)have been found to have similar composition of the reacting mixture.One should note, however, that the anomaly for the d1.2 structure appears again (it is well seen for the higher values of 3 CC ).
The reaction yield depends not only on the energy adsorption of the product molecules (3 CC ), but it is also affected by the adsorption of the reactant molecules (A and B).We have considered the changes in 3 BB only, while 3 AA is xed and equal to 100 K. Fig. 4 shows the changes of the reacting phase composition in different carbons obtained for various choices of the potential well depth for the reactant B. As expected, the rise in 3 BB causes a gradual increase of the density and mole fraction of the adsorbed B molecules.One should note that the increase of 3 BB does not affect the adsorption of the other reactant A and of the product C in the carbons from d0.5 up to d0.9.In the case of highly microporous VPCs (d1.0-d1.3), a slight decrease of the numbers of A and C molecules has been found.It may be interpreted as the effect of blocking the highenergetic adsorption centres by the adsorbed B molecules.The effect of the rise of the number of B molecules in pores (connected with the increase of 3 BB ) on the number of product molecules is rather surprising.According to the Le Chatelier's principle, an increase of the reactant concentration should lead to the increase of the product concentration.This has not been observed.However, as already mentioned, the changes of the concentration of B in the carbons d0.5 to d0.9 do not affect the concentration of C for the carbons d0.5 d0.9, and for the most microporous VPCs (d1.0-d1.3)even cause a small decrease of the C density.In order to explain this observation, a simple thermodynamic argument can be used.One can assume that the reaction occurs in the pores as well as in the bulk phase, but the transport of the product molecules between gaseous phase and pores is forbidden.The partial pressure of C in the bulk for given values of the reactants partial pressures is determined by the equilibrium state.Under given conditions, the composition of the reacting gaseous phase is almost independent on the energetic parameters of the B molecules (see the horizontal lines in the lower panels of Fig. 4), the partial pressure of C is the same in all the systems.There are the following equilibria in the system: A bulk + B bulk 4 C bulk , A bulk 4 A pores , B bulk 4 B pores , and A pores + B pores 4 C pores .Of course, the chemical potential of every component in all the phases must be the same, since the molecules in pores and in the bulk phase are in equilibrium.Despite the fact that C molecules are not transferred from the bulk phase to the pores (and vice versa), the chemical potential of this compound in both phases is determined by the reaction equilibrium and it is the same.Therefore, the amount of the product in pores is mainly determined by its hypothetical partial pressure in the bulk (almost the same for all the Fig. 4 As in Fig. 3, but for the fixed values of the equilibrium constant (G ¼ 82.3 nm 3 ) and of the potential well depth of the product (3 CC /k B ¼ 72 K), and for different values of the potential well depth for reactant B (3 BB ).The values of the reactants pressure in the bulk phase have been fixed and equal to p A ¼ p B ¼ 0.05 MPa.The horizontal lines in lower panels mark the mole fractions for the reaction taking place in the bulk phase (it should be noted that these lines are overlapped).
systems), and hence similar product densities may be expected for the given carbon.The observed small decrease of the product density in the most microporous VPCs may be interpreted as the effect of blocking the adsorption places by preferentially adsorbed B molecules.
Finally, we have also studied the effects of the bulk phase composition on reaction equilibria in pores of the selected carbons.We have xed the energetic parameters of all the components, the equilibrium constant, and assumed a constant value of the total pressure of both reactants in the gaseous phase (i.e.p A + p B ¼ 0.1 MPa).Then, we have changed the partial pressures of the components.The simulation results have been shown in Fig. 5. Apart from the obvious changes of the reactants mole fractions, we have observed a clear maximum of the product mole fraction when the partial pressures of A and B are (the same in the bulk) nearly the same in porous systems.This, of course, results from the reaction stoichiometry, and the small asymmetry of x C in porous systems results from different values of 3 AA /k B ¼ 100 K and 3 BB /k B ¼ 110 K.It is noteworthy that the product mole fraction in porous systems is considerably higher than in the bulk.

General relationships
The above discussed qualitative relationships between the reaction equilibria in pores of activated carbons with systematic changes in their porosity have shown that the adsorption energetics plays an important role.Here, we try to determine some quantitative dependencies.One of the most obvious parameters characterizing porous structures is the average pore diameter.This parameter is also closely related to the energetics of adsorption.Fig. 6-8 show the correlations between the product mole faction and the reciprocal of the average diameter of pores accessible for the product molecules (see Table 1) for different reactions.The lowest value of 1/d eff,acc,av is related to the initial VPC in the series (i.e.d0.5) and the highest values are related to the nal VPCs (i.e.d1.2 and d1.3).In the case of isomerisation reaction (Fig. 6), we have found practically linear changes of x B with 1/d eff,acc,av .In this case, the reactant (A) and the product molecules (B) are of the same size, so that they demonstrate the same geometrical constrains within the porous structure.By increasing the equilibrium constant, in the systems with all the energetic parameters xed, the reaction yield increases when the average diameter of accessible pores becomes lower.On the other hand, for a given value of the equilibrium constant, the relative energy of adsorption of the reactant and of the product determine whether the product mole fraction increases (3 AA < 3 BB ) or decreases (3 AA > 3 BB ).This reects the fact that when the reactant is more strongly adsorbed than the product, the reaction is hindered, while in the case of stronger product adsorption, the reaction yield becomes higher.
In the case of the dimerisation (Fig. 7) and synthesis (Fig. 8) reactions, the situation is different, since the product molecules are larger than the reactant molecules.In the case of systems with large contribution of wider pores, the connement effects are rather limited, which causes the linear changes of the product mole fraction versus 1/d eff,acc,av as well as versus the adsorption energy of the product molecules to appear again.However, in the case of highly microporous structures, the deviations from linearity are observed, since the accessibility of very narrow pores to the product molecules is considerably hindered.This causes that the reaction yields are rather similar for those VPCs.Fig. 9 shows the ratio of the volume accessible only to the reactant molecules (0.3500 nm # d eff < 0.4410 nm) to the volume accessible to all the molecules (d eff $ 0.4410 nm).This ratio increases when the contribution due to narrow pores becomes higher.This fact (the presence of pores inaccessible to the product) explains why the yield of dimerisation and synthesis reactions does not increase proportionally to the reduction in the pore diameters for highly microporous carbons.The highest value of the analyzed parameter is observed for the structure d1.2, and it is responsible for the already mentioned anomalies in the behaviour of this carbon.
To summarize, the chemical equilibria in activated porous carbons have been found to be inuenced by the energetic parameters of the reacting molecules as well as by the geometrical connement.When the adsorption energy of the product is higher than that of the reactant or the reactants, the reaction yield is higher than in the bulk.In this case, the reduction of pore sizes leads to higher reaction yields.In the opposite situation (the preferential adsorption of the reactant or reactants), the yield of reaction under connement is lower than in the bulk phase.In such cases, the adsorbents with wider pores are a better medium for the reactions.Besides, when the   Fig. 9 The comparison of the ratio of the volume accessible only for the reactants molecules to the volume accessible for both the reactants and the product molecules in the dimerisation and synthesis reactions.
product molecules are larger than the reactant, the presence of pores inaccessible to the product causes a reduction of the reaction yield, regardless of the energetic parameters of the reacting species.
Finally, one may raise a question about the relation between the systems with chemical reaction in pores in equilibrium with the bulk containing non-reacting reactant/reactants (such systems have been considered in this paper) and the systems with chemical equilibria in both the pores and the gaseous phase.We have performed some additional simulations, taking the isomerisation and dimerisation reactions as examples, and using the selected carbons of different microporosity (d0.5, d0.9 and d1.3).The concentration of the reactant in the bulk phase has been controlled by its pressure (we have assumed a xed value).Assuming that the reaction occurs also in the gaseous phase (cf. the discussion of the results given in Fig. 4), one can expect the concentration and the partial pressure of the product to be determined by the value of the reaction constant.These values for the bulk can be calculated easily.We have considered only a low pressure, implying that the bulk density is also low.Under such conditions, the intermolecular interactions hardly affect the equilibrium state.We have considered the systems characterized by different values of G, and performed the simulations using two different schemes.In the rst scheme, the reactant molecules (A) have been created and annihilated in the pores according to the GCMC trial moves (adsorption equilibria with the bulk phase), and the product molecules (B) have been formed or destructed only via the reaction step.In the other scheme, adsorption equilibria concerned B molecules, and A molecules were created and destructed only during the reaction.The comparison of the results obtained with the help of both schemes has been shown in Fig. S5 and S6 in the ESI.† As can be observed, the densities of the reacting molecules in pores are the same in both cases.These results suggest that the methodology of simulations used here (assuming that the reaction takes place in the pores only) is equivalent to the systems in which the reaction occurs in the pores and in the bulk phase.

Comparison with the results of other studies
As it was mentioned in the Introduction, experimental studies of chemical equilibria are rather difficult to carry out.Hence, the reports on the measurements of compositions of mixtures reacting in carbonaceous pores are also rare.0][41][42][43] There is lack of reports based on the series of materials signicantly differing in porosity.If the series of materials were considered, the main idea was to study the effects of surface functionalities or other chemical modications. 10,44,45However, in the current study we have neglected the inuence of chemical heterogeneity (the raw carbons have been considered) and we plan to study the effects due to the presence of different surface groups on chemical equilibria in the future.So, the inuence of changing in the porous structure of activated carbons has not been investigated sufficiently.According to our knowledge, there is also a lack of theoretical studies designed to evaluate directly such effects in realistic models of heterogeneous carbon pores.[14]20,21 Several reports conrmed the positive effect of connement on the reaction yield. 6,8,11,12,21owever, such conclusions are the consequence of energetics of the reacting molecules adsorption (higher energy for product molecules).In our investigation we have tested systematically different combinations of energetic parameters of the species and we have also shown that the opposite situation (i.e.lowering of the reaction yield under connement) may occur.As we have discussed above, such behaviour has been observed if the energy of the product molecules adsorption is lower than for the reactant(s).Similar observations were reported by Domínguez. 15However, his study was restricted only to the isomerisation reaction in ideal slit-like pores.It should be also noted that, despite the lower reaction yield in the discussed case, the conned reacting mixture has higher density than he corresponding bulk phase.
12][13][14][15][16][17]22,23 In the case of energetic preferences of the product adsorption, the reduction of pore diameters causes the rise in the product mole fraction.Such behaviour was also observed for the NO dimerisation in slit-like and/or cylindrical pores, 8,13,14 the NH 3 synthesis 8,12,13 and the model dimerisation reaction (A + A 4 B) in slit-like pores. 6In the latest case, the potential well depth for the reactant and the product molecules was assumed to be the same.However, the collision diameter for the B molecules was higher and, as we have discussed in the Section 3.3, it resulted in the higher interaction energy of the B molecules with the walls.Clearly, the model of ideal pores is only the rough approximation.It does not reect the complex nature of the porosity of real adsorbents.7][28][29] Such simplications can also affect the results of the modelled reactions under connement.This was demonstrated, among others, by Lísal et al. 21They compared the results of RxMC simulation of the NO dimerisation in ideal slit-like pores and in three models/replicas of disordered nanoporous carbons obtained from the reverse Monte Carlo method.The use of such realistic replicas led to signicantly different behaviour of reacting mixtures in comparison to the ideal pores.
In the context of the above-mentioned experimental difficulties and insufficiency of other results of theoretical research, the current study is the rst systematic report on the relationships between the chemical equilibria and the porous structure of carbons and the energetic parameters of reacting molecules.Despite the fact that qualitative regularities described by us are based on the simple model reactions, we believe that our conclusions are universal and they can be adapted also to the real reactions.

Conclusions
We have presented a systematic Monte Carlo study concerning the effects of connement in pores of activated carbons of different distribution of pore diameter on the reaction equilibria.We have examined how the porosity and the values of the equilibrium constant and the energetic parameters of the reacting molecules affect the reaction yield.The performed simulations have conrmed that the proposed modication of the scheme of reactive Monte Carlo simulations in the grand canonical ensemble is useful to model equilibria of reactions occurring under connement, in porous materials.This has resulted in the determination of general qualitative relationships between the porous structure and chemical equilibria, despite the fact that our simulations have been based on simplied assumptions (the simplest models of the reacting molecules).
The obtained results have shown that the reaction equilibria in pores are considerably affected by the energetic parameters of the reacting species.However, when the product molecules are larger than the reactant molecules, the geometrical connement becomes an important factor as well.
We plan to study more complex systems, involving explicitly non-sphericity of the products of dimerisation and synthesis reactions.In such cases, the geometrical connement may lead to the phenomena and orderings that can not be treated using the simplest spherical models of molecules.

Fig. 1
Fig. 1 (a) The schematic representation of the considered VPCs (the frames reflect the size of the simulation box).(b) The Histograms of the pores diameters obtained with the BG method (subsequent histograms are shifted by 0.15 in comparison with the previous ones).

Fig. 2
Fig. 2 The comparison of the composition of the reacting phase (T ¼ 298 K) in pores of all the considered VPCs for the dimerisation reaction (A + A 4 B).The data corresponding to different values of the potential well depth for product molecules (3 BB ) and the same values of 3 AA and equilibrium constant (G ¼ 41.1 nm 3).The subsequent panels present the average densities of the reactant (r A ) and the product (r B ) in pores (the densities have been calculated per the volume of accessible pores) and the mole fraction of the product in the mixture (x B ).The x B mole fractions for the reactions in the bulk phase are shown as horizontal solid lines.It should be noted that these lines are overlapped.

Fig. 3
Fig. 3 As in Fig. 2 but for the synthesis reaction (A + B 4 C).The effects of the changes in potential well depth for product molecules (3 CC ) for the fixed values of the equilibrium constant (G ¼ 82.3 nm 3 ) and the potential well depth of the reactant B molecules (3 BB /k B ¼ 110 K).The upper panels present the average densities of the reacting species (r A , r B and r C ).The lower panels show their mole fractions in the mixture (x A , x B and x C ).The mole fractions for the reaction taking place in the bulk phase are marked by horizontal lines.

Fig. 5
Fig.5Effects of the bulk phase composition (total pressure of the reactants in the bulk phase p A + p B ¼ 0.1 MPa) on the composition of reacting phase in pores of selected VPCs (d0.5, d0.9 and d1.3) for the synthesis reaction (A + B 4 C).The applied values of parameters: G ¼ 82.3 nm 3 , 3 BB /k B ¼ 110 K, and 3 CC /k B ¼ 72 K.The average densities and mole fractions of the reacting species in pores have been shown in upper and lower panel, respectively.The mole fractions for the reaction in the bulk phase are presented as lines.

Fig. 6
Fig.6The changes of the product mole fraction (x B ) for the isomerisation reaction versus the reciprocal of the average diameter of the pores accessible for the product molecules (1/d eff,acc,av ), obtained for the systems presented in Fig.S1 † (a) and in Fig.S2 † (b).The meaning of symbols as in Fig.S1and S2 in the ESI.† The dashed lines have been drawn to guide the eye.

Fig. 7
Fig.7The changes of the product mole fraction (x B ) for the dimerisation reaction versus the reciprocal of the average diameter of the pores accessible for the product molecules (1/d eff,acc,av ) for all the systems presented in Fig.S3in the ESI † (a) and in Fig.2 (b).The meaning of symbols as in Fig.S3† and 2.

Fig. 8
Fig.8The changes of the product mole fraction (x C ) for the synthesis reaction versus the reciprocal of the average diameter of the pores accessible for the product molecules (1/d eff,acc,av ) for all the systems presented in Fig.S4in the ESI † (a) and in Fig.3 (b).The meaning of symbols as in Fig.S4† and 3.

Table 1
The characteristics of the studied VPCs VPC N S a d eff $ 0.3500 b nm d eff $ 0.4410 c nm V acc ðd eff $ 0:3500 nmÞ À V acc ðd eff $ 0:4410 nmÞ V acc ðd eff $ 0:4410 nmÞ The number of carbon atoms in the simulation box.b a