A Grand Canonical Monte Carlo Simulation for the Evaluation of Pore Size Distribution of Nuclear-Grade Graphite from Kr Adsorption Isotherms

: Characterizing materials with low surface areas or with very small sample sizes requires innovative approaches beyond traditional N 2 and Ar adsorption measurements. The measurement of Kr adsorption isotherms is often employed to serve this purpose, yet its potential remains limited by the lack of models for the interpretation of the experimental results in terms of pore size distribution. In this work, simulated adsorption isotherms of Kr onto graphite in slit-shaped pores are generated with a Grand Canonical Monte Carlo method. The pore size distributions of nuclear-grade graphite samples and activated carbon are modelled by ﬁtting simulated isotherms to the experimental data. The resulting distributions are favourably compared with those generated by commercially available modelling packages, based on the use of N 2 adsorption isotherms using GCMC and BJH methods. The new GCMC-Kr kernel developed in this study offers an alternative method for the evaluation of the distribution of pore sizes in nuclear graphite and other low surface area materials, which can be employed when N 2 and Ar adsorption measurements cannot be carried out.


Introduction
Low-temperature gas adsorption isotherms, in conjunction with the BET (Brunauer, Emmett and Teller) method, have long been used as a standard procedure for the evaluation of the surface area of solids [1,2].N 2 at its boiling point of 77.4 K is the standard adsorptive of choice for surface area characterization [3][4][5].Ar and CO 2 adsorption at 87.3 K and 273 K, respectively, can be used in the study of materials with pores smaller than 1 nm.In the case of very low surface area materials, such as nuclear grade graphites, Kr adsorption is measured at 77.4 K [6,7].Historically, it has been shown that the adsorption of Kr onto graphite surfaces at temperatures below its triple point results in the formation of isotherms that display several clearly defined steps [8][9][10].More recently, this stepwise adsorption isotherm, has been attributed to the phase transitions which occur during the formation of successive adsorption layers, and has been the object of extensive experimental and modelling studies [11][12][13][14][15].
Gas adsorption data can be used for the calculation of the pore size distribution (PSD) of porous materials.In the case of mesoporous solids, capillary condensation takes place in the pores; using a modified Kelvin equation to describe the liquid-vapour equilibrium in the pore volume, it is possible to perform PSD calculations, as in the Barrett, Joyner and Halenda (BJH) method [16].Modern methodologies, such as the Non-Localized Density Functional Theory (NLDFT) [17][18][19] and Grand Canonical Monte Carlo (GCMC) [14,[20][21][22], represent an alternative approach for the evaluation of PSD from gas adsorption data.Gas adsorption and phase equilibria in simulated porous materials with various pore geometries, such as cylinders, hexagonal pores, wedges and slits can be investigated using GCMC [11,14,23,24].

Materials and Methods
The adsorption isotherm of Kr onto a non-porous graphitized carbon black (GBC) was studied in order to validate the chosen model of the interactions between gas particles and the graphite's solid surface.The GCB was supplied by MicrotracBEL (Osaka, Japan) and is used as the standard calibration sample for the gas adsorption instrument.Its specific surface area is 49.52 ± 1.5 m 2 g −1 .
Two different nuclear graphite grades, IG110 and IG430, produced by Toyo Tanso Ltd., Osaka, Japan, were chosen as examples of low-surface area graphite materials of industrial interest.IG110 and IG430 are fine-grained, isotropic graphites which exhibit high thermal durability and strength and are consequently candidates for use as a moderating and structural material in the core of Generation-IV nuclear reactors.The density and porosity values of these nuclear graphites can be found in the Supplementary Information in Table S1.The characterization of the porous space of nuclear-grade graphites is a complex endeavour [35][36][37].The complex porous structure of these polycrystalline materials can be seen in the scanning electron micrographs provided in the Supplementary Information in Figures S1 and S2.These show porous features as large as tens of µm, too large to be characterized by gas adsorption, and previous studies [36] identified a pore size distribution centered around 2 µm.However, a small volume percentage of features in the meso-and micro-porous range are also present, and these are the focus of this work.Despite its difficulties, reaching a better quantitative understanding of nuclear graphite's voids, including the small percentage of pores and cracks smaller than 100 nm, and of their evolution as irradiation and oxidative phenomena take place, is highly important to correctly predict the ageing of graphite components in a nuclear reactor [38].It is often difficult to obtain large amounts of irradiated nuclear graphite samples, and the more traditional N 2 adsorption fails to produce reliable isotherms on small amounts of such low surface area materials.The use of Kr and the GCMC-generated kernel developed in this study offer a valid alternative for the characterization of these types of samples.
Finally, in order to evaluate the applicability of the method presented in this work to carbon materials with higher surface area, a sample of activated carbon (AC) Darco KB 100 mesh (Aldrich, Dorset, UK) was also characterized using the methodology described in the following section.

Experimental
All samples were dried under vacuum for a minimum of 3 h at 305 • C using the BELPREP II-vac (MicrotracBEL, Japan) in order to remove all traces of moisture from their pore volume.Adsorption and desorption isotherms of N 2 and Kr at 77.4 K were obtained using a BELSORP-Max instrument (Microtrac-BEL, Japan).The system is fitted with three pressure transducers with scale ranges of 0.0133 kPa, 1.33 kPa and 133 kPa, respectively.This allows for high-resolution measurements to be obtained over extremely low relative pressure regions.N 2 isotherms were measured for pressures up to 101 kPa.Kr adsorption data were measured for pressures up to around 220 Pa.This is a much lower value of maximum relative pressure compared to that reached during N 2 measurements.Beyond this pressure, at the temperature used in these experiments, which is below the triple point of Kr, condensation takes place in the porous structure of the samples, and adsorption measurements are no longer reliable.GCB was characterized using only Kr adsorption, while the graphite samples and the activated carbon were analyzed using both gases.In order to overcome the issues associated with the acquisition of a reliable N 2 adsorption measurement for low surface area materials, a large amount of sample, in excess of 3 g, was analyzed using a large sample chamber with a volume of 5 cm 3 .This was possible because the samples under investigation were unirradiated virgin graphite samples, so sample size and availability were not an experimental constraint.Standard sample chambers, employed in the case of materials with higher surface area, have a volume of 1.8 cm 3 .For each sample, 3 replicate measurements were carried out.
N 2 adsorption isotherms were used as input to the BELMaster software (Microtrac-BEL, Japan), which provides several models for the evaluation of the surface area and PSD, including BET, BJH modelling of cylindrical mesopores and GCMC simulation of N 2 onto slit-shaped graphite micro-to macropores.Kr adsorption isotherms were modelled using the GCMC simulation described in the following section.

Modelling
In this work, simulated adsorption isotherms of krypton confined into slit-shaped pores were generated using a Grand Canonical Monte Carlo method, described in greater detail elsewhere [39,40].While pores of different shapes could be present in carbonaceous materials [41,42], the pore shape chosen in this work is consistent with work by Dubinin and Plavnik [43], which shows that activated carbons obtained from vegetable feed consist of nanocrystalline carbon.Analysis of the results of N 2 adsorption/desorption hysteresis on graphite, presented in Section 3, supports the choice of pore shape adopted in this study [36].
A simulation box with periodic boundary conditions in the x and y directions and with planar solid walls delimiting its width in the z direction is in equilibrium with a reservoir of gas particles.The volume V, temperature T and chemical potential µ of the system are set.The system is in equilibrium with an infinite reservoir of particles, and the pressure P is calculated from the chemical potential using the bulk equation of state of the ideal gas [44].The GCMC algorithm generates a large number of system configurations with a Markov chain process of random insertion, removal or translation of fluid particles in the simulation box.The probability of each move being accepted is determined by the potential energy change of the system caused by the move.
Krypton is modelled using a Lennard-Jones 12-6 potential: the interaction between 2 fluid particles, i and j, is given by: where r ij is the distance between the 2 fluid particles, Kr−Kr is the potential well depth and σ Kr−Kr is the kinetic collision diameter.The cut-off distance for the LJ potential is set as 5 σ Kr−Kr .Sophisticated models of adsorbent-adsorbate interaction potentials, taking into account electrostatic and induction effects, have been developed [45].However, in this work, due to the non-polar nature of Kr, the interaction between each fluid particle i and a solid graphite plane wall is modelled with a Steele 10-4-3 potential [46]: where ρ C = 114 nm −3 is the density of carbon atoms per unit volume of a graphene layer, ∆ = 0.335 nm is the distance between adjacent graphitic layers making up the graphite wall and z i is the distance of a gas particle i from the solid wall.The potential parameters Kr−C and σ Kr−C of the solid-fluid interaction are calculated from the parameters of the Kr-Kr and C-C interactions using Lorentz-Berthelot mixing rules.The values of the potential parameters used in the GCMC calculations are shown in Table 1.In the case of a slit-shaped pore, there are two walls and the total potential due to the interaction between a gas particles and the graphite solid walls is given by: where H is the distance between the carbon centres across the pore.The effective pore width accounts for the reduction in accessible pore volume caused by the physical size of the carbon atoms in the pore walls: where z 0 is the distance of closest approach of the adsorbate to the solid wall, calculated as the root of the adsorbate-solid wall potential equation [31,47].Following a similar approach to that of Miyahara et al. [48], the simulation box has dimensions 70 σ Kr−C , 10 σ Kr−C and H in the x, y and z direction respectively.Initial attempts to run less time-consuming simulations with a smaller number of equilibration and sampling steps (between 10 4 and 10 5 each) resulted in poorly simulated adsorption isotherms, as evidenced by the decrease of the adsorptive capacity with an increase in pressure at some pressure points, and by poorly defined adsorption step-wise behaviour.This is caused by the large number of rejected insertion/deletion attempts during equilibration and sampling.In this work, the simulation of each data point of each adsorption isotherm was run for 2 × 10 6 equilibration steps, followed by 5 × 10 6 sampling steps.
Fitting the experimental adsorption isotherm to the simulated kernel in order to determine the PSD is an ill-posed discrete optimization problem: the adsorption isotherm integral equation is [33] b where b(P) is the experimental adsorption isotherm, K(H, P) is the kernel of simulated isotherms for different pore sizes H and x(H) is the pore size distribution.The work by Herdes et al. [49] details the numerical procedure needed in order to achieve a physically meaningful solution to this inversion problem by employing the Tikhonov regularization method.
In this work, the optimal PSD was obtained using CVXPY [50], a Python library for the resolution of convex optimization problems.

Interaction between Kr and Graphite
The results of GCMC adsorption simulations are strongly dependent on the potential chosen to describe the interaction between adsorbent and adsorbate, and on the values of the potential parameters.Therefore, an initial validation of the chosen potential and parameters (described in Section 2.2) was carried out.A GCMC simulation of Kr adsorption on a single graphite surface, modelled using the Steele potential, produced an adsorption isotherm that was compared to historical data on exfoliated graphite [9,51], shown in Figure 1.The modelling results were also compared to data obtained experimentally on non-porous graphitized carbon black, Figure 2. Figure 1 shows an excellent agreement between the simulated and experimental data at very low pressure, as the modelled data correctly predict the formation of the first adsorption layer.As the pressure increases, the modelled system shows the formation of a second and third adsorption layers, closely matching those of the experimental data.A similar situation can be observed in Figure 2; in this case, a minor discrepancy in the pressure onset of the first adsorption layer, while the model's predictions of the second and third adsorption steps are in very good agreement with the measurements carried out on GCB.
The differences between the experimental and modelled adsorption isotherms can be explained by taking into account factors not included in the simplified structureless Steele potential, such as the presence of nano-scale void features between the graphite crystallites [11,13] and the apparent porosity created by the packing of the GCB.The Steele potential does not directly account for the interactions between each C atom on the graphite surface and each Kr atom in the fluid phase; however, the chosen potential is able to quantitatively model the step-wise behaviour of the adsorption of Kr on a non-porous graphite material, successfully simulating the 2-D transitions occurring as the successive adsorption layers are formed.It is also computationally faster than an explicit atomistic representation of the solid phase, which also supports its use in this work.

Slit Pores and GCMC-Generated Kernel
A set of adsorption isotherms for slit-pores was constructed.This comprises 95 isotherms covering pore sizes from 0.87 to 147 nm, each simulating adsorption at 75 pressure values and a temperature of 77.4 K, matching the temperature of the experimental adsorption isotherms.Figure 3a shows Kr adsorption within a simulated slit pore of width 3.7 nm at a pressure of 260 Pa.In order to clearly visualize the Kr particles, the pore walls are not shown in the Figure .It is possible to observe how Kr formed highly ordered layers on top of the pore walls, with localized structures that can be clearly recognized as close hexagonal packing in the xy plane, shown in Figure 3b.As well as the initial 2-D transitions corresponding to the multi-layer adsorption of Kr onto the graphite walls, Figure 3a shows the formation of a liquid bridge with 2 menisci.This is in agreement with the simulations shown in the work of Miyahara and Tanaka [24], who observed the formation of a simulated liquid bridge of N 2 in graphite pores.In their work on phase equilibria in confined spaces, they adopted a potential field which linearly decays to zero at the extremities of the simulation cell.This results in a liquid bridge centred in the middle of the simulation volume, where the full potential field acts, and it allows for the precise evaluation of liquid-vapour equilibrium in the pore.The liquid bridge in Figure 3a is not perfectly centred in the middle of the pore along the x axis, as in this work, the full potential field acts over the entire simulation cell; however, the focus of this work is on the evaluation of PSD distribution from Kr adsorption data and not the study of phase equilibria, and extending the full field potential to the entire simulation volume was considered to be a suitable approach.The local density of Kr particles as a function of their distance in the z direction from the centre of the pore is shown in Figure 4.The density distribution follows the trend shown by other authors [22,24,32], with clearly identified peaks where the successive adsorption layers form.The distance between the peaks is close to krypton's kinetic collision diameter σ Kr−Kr .
The 95 isotherms forming the GCMC-generated kernel constructed for this work are shown in Figure 5.The vertical axis shows the number of gas particles adsorbed per unit volume of the pore, while the horizontal axis shows the gas pressure.As the narrowest pores can only accommodate one layer of Kr particles within their void volume, they display one single sharp adsorption step, occurring at increasing pressure values for increasing pore width.In particular, the pore with a width of 0.87 nm, the narrowest in this kernel, is entirely filled with Kr even at very low pressure, and its adsorption isotherm is a nearly horizontal straight line.This suggests that, for narrow micropores smaller than 0.87 nm, the simulation method is not sensitive enough, and such pores should not be included in the kernel for the evaluation of PSD of porous materials.Intermediate pore sizes show a reduction in the magnitude of the first adsorption step.However, a second and, as the pore width increases, a third adsorption step become clearly visible, showing the formation of successive adsorption layers.At larger pore sizes, the interactions with the pore walls only affect the particles closer to the solid surfaces, while the particles in the central part of the void volume are virtually unaffected.This results in a further reduction in the relative magnitude of the first adsorption step.Following the formation of the first adsorption layer, the adsorption isotherms of larger pores follow a nearly linear trend with pressure.Local density of particles in a simulated 3.7 nm wide slit pore as a function of their distance from the pore centre.The particle density is expressed as a percentage of the total particles in the pore.

Pore Size Distribution of Nuclear-Grade Graphite and Activated Carbon
Kr and N 2 isotherms for IG110 and IG430 graphites and for activated carbon were obtained as described in Section 2.1.Kr adsorption curves are shown below in Figures 6-8, while N 2 adsorption and desorption isotherms studied in this work can be found in the Supplementary Information in Figures S3-S5.Replicate measurements were taken for each sample, but as no discernible differences could be identified between replicates, only one of the replicates is shown.The desorption hysteresis loop shown by the graphite samples is classified as type H3, which is usually associated with plate-like particles and slit-shaped pores [6].This supports the choice of adopting a slit geometry in our simulation of the void features of these samples.Activated carbon displays large initial adsorption of Kr at low pressure.It also shows an H3 desorption hysteresis.The presence of hysteresis on the desorption branch of the isotherm is associated with the presence of mesoporosity and demonstrates adsorption of Kr at high relative pressures had no limits [52].This confirms that, albeit as a small percentage of the overall void volume for the graphite samples, mesopores are present in all these samples, while it is more likely for the activated carbon pores to be in the microporous range.Table 2 shows the BET-specific surface area of the samples.Closely matching values of specific surface area were derived from the use of both gases in the case of the low surface area graphite samples.The BET surface area measured for activated carbon is 4 orders of magnitude larger than that of graphite; however, the BET method is not designed for microporous adsorbents such as activated carbon, and the results should be interpreted as an apparent surface area [6].The kernel of simulated Kr adsorption onto slit pores of varying sizes obtained using the GCMC method was used to generate modelled adsorption isotherms and PSDs of the three samples.This was accomplished by fitting a linear combination of the kernel simulated isotherms to the experimental gas adsorption curve, as described in Section 2.2.The PSD produced from the Kr GCMC-generated kernel is also compared to PSDs generated with a commercial GCMC algorithm and BJH model from N 2 adsorption data.
Figure 6a shows the results of the PSD simulations of IG110 graphite carried out using the GCMC-Kr adsorption simulations developed in this study and the PSD generated by a GCMC-N 2 and BJH N 2 from a commercial software package.The GCMC-Kr, in agreement with the GCMC-N 2 and BJH-N 2 models, shows no pores smaller than 1 nm.The N 2derived PSDs also suggest the presence of pore volume for pores 1.5 to 2 nm wide, which is also found in the modelled PSD from simulated Kr adsorption.In the 1-30 nm range, the trends shown by all models, and in particular by the two GCMC simulations, are similar, with a peak around 3.5-4 nm and 10-12 nm; however, the volume contribution of this peak is much larger, in relative terms, in the case of the GCMC-Kr simulation than in the results from both GCMC-N 2 and BJH-N 2 .While the GCMC-Kr PSD simulation shows very little pore volume for pores larger than 30 nm, both the PSDs generated by GCMC-N 2 and BJH-N 2 suggest the existence of pores larger than 30 nm.The simulated and experimental Kr adsorption isotherms are shown in Figure 6b.The match between the model and experiment is excellent up to a pressure of 200 Pa, where the simulation fails to predict the onset of the third step on the adsorption isotherm, that can be observed experimentally.
The simulated PSDs of sample IG430 are shown in Figure 7a: similarly to what was observed for IG110, the PSD from the GCMC-Kr method matches very closely those generated by the other models in the 1-30 nm range.However, GCMC-N 2 PSD shows pores of sub-nanometre sizes, which are not present in the GCMC-Kr PSD.The GCMC-Kr results show a very small amount of porosity for pore sizes larger than 30 nm, when compared with the results from other models.Figure 7b shows the simulated and experimental Kr adsorption isotherms.The agreement is excellent for pressures up to 170 Pa, including the first and second step of the adsorption isotherm.As in the IG110 case, the third step in the adsorption isotherm develops at a pressure of around 170 Pa, while the simulated isotherm shows a linear increase in adsorption volume for pressures up to 220 Pa.
A comparison between the results of the simulation developed in this work and those produced by modelling of N 2 adsorption data was also carried out for an activated carbon with high surface area (Figure 8a).The GCMC-N 2 method produces a PSD with a significant fraction of micropores smaller than 0.87 nm, which cannot be simulated by the model developed in this study, for the reason explained in Section 3.2.The GCMC-Kr distribution shows peaks between 1.3 and 20 nm, which matches very closely that in the GCMC-N 2 distribution.The match between the results of different simulations for the remaining pore sizes is excellent, with very few pores larger than 20 nm.The experimental and simulated Kr adsorption isotherms are also closely matching.However, at low pressures, the simulated isotherm has a stepwise behaviour, while the experimental adsorption isotherm is smooth over the entire range of pressures.
Table 3 shows the total void volume modelled by the two GCMC methods and by BJH-N 2 , with values quite closely matching for each sample.In the case of the graphite samples, the GCMC-Kr modelling results in the smallest value of total void volume.All models show IG110 to have a larger total void volume than IG430.The total void volume resulting from these models is very small, and, as to be expected in the case of nuclear-grade graphites, it accounts only for a very small fraction of the overall porosity of the samples.The activated carbon total void volume values estimated by the two GCMC methods are closely matching and 3 orders of magnitude larger than those of the graphite samples.The BJH model, when applied to activated carbon, results in a lower void volume.Before discussing in detail the results of the PSD simulations, it is important to point out that the direct comparison of results produced by the two GCMC models is complicated by the fact that the fitting methodology used in the commercial software used for the interpretation of N 2 adsorption data is not known.
Another issue to be noted is the resolution of the experimental measurements [49]: the minimum pressure at which gas adsorption was measured experimentally is around 0.1 Pa. Figure 5b shows how, according to the GCMC-Kr simulations, pores smaller than 1 nm would be completely filled at pressures below 0.1 Pa.In some of the larger pores, a mono-layer has already formed at this pressure.By obtaining higher resolution experimental Kr adsorption isotherms, it would be possible to gain better predictions of the PSD of porous materials.The graph of GCMC-N 2 kernel published in the work of Miyahara et al. [48] suggests that the same issue could also affect, to a lesser extent, the interpretation of N 2 adsorption data, and consequently, the PSD generated by such methods at the sub-nanometre scale.This explains the discrepancies between the PSDs obtained from the two GCMC methods for narrow pores (H < 1 nm) observed in the case of IG430 and activated carbon.The resolution of the experimental data causes the GCMC-Kr model to fail in modelling such pores; consequently, the simulation attributes their volume contribution to larger pores.However, it should also be noted that there is a level of uncertainty in the quantification of the volume contribution from the narrowest pores as determined by GCMC-N 2 modelling, especially in the case of IG430.The highest magnification SEMs shown in the Supplementary Information (Figures S1 and S2) show a randomly arranged structure of relatively large graphite crystals.While it is not possible directly to observe or quantify the presence of porous features smaller than 1 nm, it is likely that these features are mostly due to gaps between the graphite crystals, exposed crystal edges or, to a lesser extent, very small features on the surface of the crystals.So, while the presence of void features smaller than 1 nm cannot be visually verified from image analysis, their contribution to the overall porevolume of the samples is likely to be very small.It is, therefore, possible to postulate that the commercial GCMC-N 2 method is overestimating the presence of the narrowest pores.The shortcomings of both N 2 and Kr experimental adsorption measurements in probing the narrowest voids could cause the differences in the distributions of the narrowest pore sizes, resulting from the two corresponding simulation methods, to be exacerbated.
Figures 6a and 7a show good agreement between the two GCMC models, with similar PSDs for pore widths up to 50 nm in IG110 and IG430 samples.For pores wider than 30 nm, the GCMC-Kr seems to underestimate the pore volume when compared to the results produced by modelling based on N 2 data.It could be postulated that Kr adsorption isotherms do not contain sufficient information about wider pores.One of the contributing factors could be capillary condensation occurring in the voids between the crystallites; the presence of heterogeneous sites in larger pores may induce condensation.The experimental adsorption isotherms of IG110 and IG430 show sharp steps in the amount of adsorbed Kr at pressures of around 170 Pa, which become completely vertical for pressures greater than 220 Pa.The condensation taking place in the complex graphite void space cannot be predicted by the GCMC simulations, as shown in Figures 6b and 7b, which generate linear isotherms for pressures greater than 100 Pa.The inability of the GCMC-Kr model developed in this study to simulate the adsorption isotherms of graphite for the entire pressure range also explains why its estimate of the total void volume of graphite samples is smaller than the total pore volume obtained by GCMC and BJH methods from N 2 data.
The PSDs of activated carbon generated by the two models show the best agreement, especially for pores larger than 1 nm.The overall void volume predicted by the two methods is also closely matching, as, in this case, the GCMC-Kr model can closely simulate the adsorption isotherm over the entire range of pressures.However, there is a noticeable difference between the simulated (stepwise) and experimental (smooth) Kr adsorption isotherms at low pressure.This can be explained by the lack of preferential adsorption sites on the smooth walls of the simulated slit pores compared to the non-uniform surfaces and exposed layer edges found in activated carbon [11,51].These non-uniformities lead to a smooth adsorption isotherm that cannot be reproduced by the GCMC model in its current formulation, characterized by the sharp step of Kr uptake in smaller pores.In order to reduce the stepwise behaviour of adsorption in the case of activated carbon, it would be necessary to modify the model to take into account the roughness of the pore surface.

Conclusions
We presented the development of a kernel of simulated Kr adsorption isotherms generated using a GCMC adsorption model of Kr in slit-shaped carbon pores.Simulated adsorption isotherms and pore size distributions of porous graphite and activated carbon samples were generated by fitting the simulated kernel to experimental Kr adsorption data.The pore size distributions were compared with those generated by other modelling methods requiring N 2 adsorption isotherms.The agreement between the method developed in this work and other methodologies, which is good in the pore size range from 1 to 50 nm, confirms that Kr adsorption and the kernel we developed offer an alternative when the traditional characterization by N 2 adsorption and the standard PSD models available in commercial software fails to achieve meaningful or reproducible results for low-surface area samples, such as nuclear graphite samples.

Figure 1 .
Figure 1.Comparison of GCMC modelling results for adsorption of Kr on a single graphite wall with experimental adsorption isotherm of Kr on exfoliated graphite, as measured by Thomy et al. [9] (also shown in Diao et al. [13]).

Figure 2 .
Figure 2. Comparison of GCMC modelling results for adsorption of Kr on a single graphite wall with experimental adsorption isotherm of Kr on graphitized carbon black.

Figure 3 .
(a) 3-D representation of a simulated 3.7 nm wide slit-shaped pore with Kr adsorption taking place onto the solid graphite walls.The pore length is 25.9 nm, and its depth is 3.7 nm.(b) Detail of the first adsorption layer, showing close hexagonal packing of the particles.

Figure 4 .
Figure 4. Local density of particles in a simulated 3.7 nm wide slit pore as a function of their distance from the pore centre.The particle density is expressed as a percentage of the total particles in the pore.

Figure 5 .
GCMC-generated kernel of Kr adsorption isotherms, expressed as the number of adsorbed gas particles per nm 3 of slit-shaped pore volume.Some of the widths of the pores are identified on the graph.(a) Linear pressure axis.(b) Logarithmic pressure axis.

Figure 6 .Figure 7 .Table 3 .Figure 8 .
IG110 graphite.(a) Comparison between PSDs generated by GCMC and BJH models from N 2 isotherms and the results of the algorithm developed for the interpretation of Kr isotherms.(b) Experimental and simulated Kr adsorption isotherms.IG430 graphite.(a) Comparison between PSDs generated by GCMC and BJH models from N 2 isotherms and the results of the algorithm developed for the interpretation of Kr isotherms.(b) Experimental and simulated Kr adsorption isotherms.Simulated void volume of IG110, IG430 and activated carbon derived from different modelling approaches.Activated carbon.(a) Comparison between PSDs generated by GCMC and BJH models from N 2 isotherms and the results of the algorithm developed for the interpretation of Kr isotherms.(b) Experimental and simulated Kr adsorption isotherms.

Table 1 .
Potential [13]meters used in GCMC simulations.kB is the Boltzmann constant.The values of the Kr-Kr and C-C parameters are those found in Diao et al.[13]The Kr-C parameters are obtained by Lorentz-Berthelot mixing rules.

Table 2 .
BET specific surface area of IG110, IG430 and activated carbon derived from Kr and N 2 adsorption data.