Biogas Cleaning via Vacuum Swing Adsorption Using a Calcium Metal–Organic Framework Adsorbent: A Multiscale Simulation Study

: Purifying biogas can enhance the performance of distributed smart grid systems while potentially yielding clean feedstock for downstream usage such as steam reforming. Recently, a novel anion-pillared metal–organic framework (MOF) was reported in the literature that shows good capacity to separate acetylene from carbon dioxide. The present study assesses the usefulness of this adsorbent for separating a typical biogas mixture (consisting of methane, nitrogen, oxygen, hydrogen, carbon dioxide, and hydrogen sulphide) using a multiscale approach. This approach couples atomistic Monte Carlo simulations in the grand canonical ensemble with the batch equilibrium modelling of a pressure swing adsorption system. The metal–organic framework displays selectivity at low pressures for carbon dioxide and especially hydrogen sulphide. An analysis of adsorption isotherm models coupled with statistical distributions of surface–gas interaction energies determined that both CH 4 and CO 2 exhibited Langmuir-type adsorption, while H 2 S displayed Langmuir-type behaviour at low pressures, with increasing adsorption site heterogeneity at high pressures. Batch equilibrium modelling of a vacuum swing adsorption system to purify a CH 4 /CO 2 feedstock demonstrated that such a system can be incorporated into a solar biogas reforming process since the target purity of 93–94 mol% methane for incorporation into the process was readily achievable.


Introduction
An attractive route for future distributed energy systems as part of a so-called "smart grid" is the adoption of micro-generation or micro-cogeneration systems at the scale of small to medium-sized companies, often within the context of cogeneration whereby a firm may generate both electricity and steam on-site and employ the steam for heating purposes.The use of a biodiesel-biogas fuel compression ignition engine was shown to have superior performance in terms of emissions [1], and such engines can be employed in the generation of electricity or cogeneration processes.The cogeneration of steam and electricity using biogas derived from pig manure and corn silage blends is an area of recent research [2] for which lifecycle aspects such as ozone depletion were examined for a variety of feedstocks, with a view to implement such processes on-site at pig farms.The micro-cogeneration of a compression ignition internal combustion engine using a lean homogeneous air/biogas mixture has been shown to be practicable for a company operating a mid-sized farm [3].Biomethane fuel obtained from biogas using a solid matrix of wood ash has been shown to be feasible for the purposes of injecting biomethane into the gas grid [4].The feasibility of biogas as an alternative fuel to generate electricity via a Stirling engine, within the context of distribution electricity production, has been shown in the literature [5].Fuel for some of the aforementioned systems could often be in the form of biogas, such as from agricultural wastes, manure, and domestic waste.The composition of biogas is highly variable, as outlined in Table 1, and there are often significant quantities of environmentally harmful compounds such as carbon dioxide and hydrogen sulphide.
Table 1.Typical biogas composition [5].In addition to the obvious threat to the environment posed by carbon dioxide, hydrogen sulphide is not only toxic, but, in the presence of water, it can cause sulphide stress cracking as it produces corrosive aqueous solutions.Natural gas naturally contains some hydrogen sulphide, with some wells yielding gas with a H 2 S content of over 10% [5,6].Biogas, on the other hand, typically contains up to 0.01% hydrogen sulphide, although it may contain up to 40% carbon dioxide [6], greatly lowering the calorific value of the gas.Moreover, biogas combustion in micro-generation activities will, of course, result in the production of carbon dioxide; hence, removing CO 2 pre-combustion will not only improve the calorific value of the biogas fuel, but will prevent downstream carbon dioxide emissions.

Component
Established biogas upgrading technologies include water scrubbing, physical organic scrubbing, chemical absorption, pressure swing adsorption, membrane separation, and cryogenic separation.These are extensively reviewed in the literature including good illustrations of their respective advantages and disadvantages along with techno-economic analyses examining features such as energy requirements, methane purity, and recovery and cost economics [7][8][9][10].The reader is directed to these thorough review articles for an exhaustive assessment of recent developments in biogas upgrading.Water scrubbing is the simplest and easiest process to operate.However, water contained in biogas often leads to pipe corrosion and clogging due to oxygen presence and bacteria growth, respectively.Furthermore, additional costs arise from the need to dry the final methane-rich stream.While adsorption does not require any heat or harmful chemicals, it exhibits higher methane losses than water scrubbing.The need to address this limitation by developing new adsorbents with improved CO 2 or H 2 S-to-methane selectivity motivated the investigation of a metal-organic framework, namely, CaMOF, in the present study.
Metal-organic frameworks (MOFs) represent a class of hybrid organic-inorganic porous materials in which an organic ligand (a nitrogen or oxygen donor) and a metal ion (alkaline metal, transition metal, or a lanthanide) are coordinated to form a tridimensional structure [11,12].Having emerged nearly three decades ago, MOFs are also referred to as porous coordination polymers, porous coordination networks, microporous coordination polymers, zeolite-like metal-organic frameworks, metal peptide frameworks, metal biomolecule frameworks, or mesoporous metal-organic frameworks [13].Owing to their advantageous properties such as their porous structure, high specific surface area, good structural stability, high tunability, the existence of active sites, and being readily available, these materials attracted researchers' attention in numerous fields.To date, the known applications of MOFs include drug delivery [14] (Lawson et al., 2021), optoelectronic materials [15], energy storage devices [16], catalysis [17], carbon capture and carbon dioxide conversion to valuable chemicals [18], wastewater treatment [19], and the separation of liquid [20] as well as gas mixtures [21].In relation to gas separation using MOFs, both membrane and adsorptive processes have been investigated [13].Selective adsorption, more especially for the separation of light gases, stands out as one of the most attractive areas of research related to MOFs.This is due to the fact that these gases (oxygen, nitrogen, carbon dioxide, methane) are relevant to environmental pollution remediation scenarios as well as renewable energy systems.These two research themes draw much attention from researchers and the public in the context of worldwide environmental concerns.Examples of extensively investigated systems include CO 2 /CH 4 and CO 2 /N 2 [22], CO 2 /H 2 [23], C 2 H 4 /C 2 H 6 [24], C 3 H 6 /C 3 H 8 [25], and n-C 4 H 10 /i-C 4 H 10 [26].Li and coworkers (2012) [13] and Lin and co-workers (2020) [27] provided further detail on numerous other separation problems examined in the literature.From both experimental and modelling studies, it emerged that structure and pore properties were identified through these studies as the predominant factors determining separation performance, including the selectivity of MOFs.Furthermore, MOFs were generally found to perform better than zeolites.This stems from challenges in the control of pore structure in zeolites [27].
Previous work has examined the upgrading of biogas using activated carbon [28] and metal-organic frameworks (MOFs) [29] via pressure swing adsorption to enhance the calorific value.With regard to MOFs in particular, core-shell type ionic liquid/MOF composites have been shown to possess exceptionally high selectivity for carbon dioxide over methane [30], while Mg-MOF-74 can be used at low pressure and is widely employed for carbon dioxide adsorption [31].MOF/polymer membrane composites have been shown to be effective for biogas separation [32][33][34][35].
This study's primary focus is on removing undesirable components from biogas using a recently developed Ca 2+ metal-organic framework (CaMOF) with a honeycomb framework structure [30], as shown in Figure 1 (the crystallographic information file is available in the Supplementary Materials).This MOF possesses high-density organic hydrogen-bonding acceptors and open-metal sites and is already known to be highly effective for separating acetylene from carbon dioxide; hence, it was selected for this exploratory investigation into its potential application as an adsorbent for cleaning biogas.Biogas cleaned of undesirable components can be beneficial in several ways.Firstly, with non-calorific species such as carbon dioxide and nitrogen removed, the heating value of biogas as a fuel can be increased.Secondly, the removal of so-called "sour" components (i.e., carbon dioxide and hydrogen sulphide) can improve the utility of biogas as a feedstock for steam reforming processes to produce synthesis gas and downstream chemicals using renewable biogas as a feedstock.A potential problem with respect to H 2 S removal, however, is its tendency to form corrosive solutions in the presence of moisture.This can, of course, limit the lifespan of materials, and so would remain a challenge in terms of the long-term operation of gas sweetening systems that employ the CaMOF or other MOFs.
The MOF in question has a unit cell with parameters a = 3.18746 nm, b = 3.18746 nm, c = 0.89890 nm, α = β = 90.0• , and γ = 120 • and the chemical formula C 162 H 180 N 90 O 54 Ca 18 .Coordinatively unsaturated metal sites, known as open-metal sites, are a common feature of MOF materials and are a key characteristic in terms of enhanced [31][32][33][34][35] and selective [36][37][38][39] gas adsorption at low pressures or fugacities.These sites may be unsaturated and made accessible via framework defects such as missing organic linkers [40][41][42] or, commonly, through chemical [43,44] or thermal [45][46][47][48][49] activation protocols.The CaMOF considered in this study has a honeycomb porous framework with regular 1D channels.The Ca 2+ cations have a pentagonal bipyramid coordinated geometry and are connected to one oxygen atom in DMA, three nitrogen atoms in the dtztp ligand, and three carboxylate oxygen atoms in the dtztp ligand.The aforementioned channels are lined with exposed Ca 2+ centres along with uncoordinated nitrogen atoms such that the effective window diameter is 1.5 nm and there is an accessible voidage of 68.9% [30].
The aim of this study is to assess the feasibility of the newly developed CaMOF to separate undesirable sour components from biogas, specifically CO 2 and H 2 S. Multicomponent molecular simulations of biogas blended with this adsorbent have not been published in the open literature.The batch equilibrium modelling of a pressure swing adsorption system using this adsorbent, another topic of analysis in this study, has also not been published in the open literature.This study also presents the first investigation into the incorporation of such a pressure swing adsorption system using the CaMOF into a biogas reforming process.[31][32][33][34][35] and selective [36][37][38][39] gas adsorption at low pressures or fugacities.These sites may be unsaturated and made accessible via framework defects such as missing organic linkers [40][41][42] or, commonly, through chemical [43,44] or thermal [45][46][47][48][49] activation protocols.The CaMOF considered in this study has a honeycomb porous framework with regular 1D channels.The Ca 2+ cations have a pentagonal bipyramid coordinated geometry and are connected to one oxygen atom in DMA, three nitrogen atoms in the dtztp ligand, and three carboxylate oxygen atoms in the dtztp ligand.The aforementioned channels are lined with exposed Ca 2+ centres along with uncoordinated nitrogen atoms such that the effective window diameter is 1.5 nm and there is an accessible voidage of 68.9% [30].
The aim of this study is to assess the feasibility of the newly developed CaMOF to separate undesirable sour components from biogas, specifically CO2 and H2S.Multicomponent molecular simulations of biogas blended with this adsorbent have not been published in the open literature.The batch equilibrium modelling of a pressure swing adsorption system using this adsorbent, another topic of analysis in this study, has also not been published in the open literature.This study also presents the first investigation into the incorporation of such a pressure swing adsorption system using the CaMOF into a biogas reforming process.

Methods
To mimic adsorption, Monte Carlo simulations following the Metropolis scheme [50] were employed in the grand canonical ensemble wherein the chemical potential of the fluid reservoir, adsorption cell volume, and system temperature (i.e., T = ambient temperature = 298 K) are fixed.The adsorbent atoms are fixed within the simulation cell while

Methods
To mimic adsorption, Monte Carlo simulations following the Metropolis scheme [50] were employed in the grand canonical ensemble wherein the chemical potential of the fluid reservoir, adsorption cell volume, and system temperature (i.e., T = ambient temperature = 298 K) are fixed.The adsorbent atoms are fixed within the simulation cell while the configuration of adsorbate molecules changes by the application of random moves that mimic real-world behaviour (with their probabilities of occurrence in parentheses): creation (23%), deletion (23%), rotation (24%), translation (24%), and regrowth (6%).The first two moves mimic adsorption and desorption, respectively, while the other three moves mimic the thermal motion of the adsorbed molecules within the simulation cell.Conversion between the chemical potential and the pressure of the gas was undertaken using the Peng-Robinson cubic equation of state [51], which is widely used in natural gas processing and is one of the leading two-constant cubic equations of state [52].The details of the Metropolis scheme approach to constructing Markov chains from which equilibrium properties can be determined are described in detail in the literature [53].To ensure equilibration of the system, no less than 2 × 10 6 moves were used to equilibrate, followed by a further 2 × 10 7 moves to generate results.Averaging across three independent simulations was undertaken for each data point in the adsorption isotherms, with the standard deviation being used for the estimated uncertainty.Molecular simulations were performed using Dassault Systèmes BIOVIA software 2020: in particular, BIOVIA Materials Studio 2020 [54] was used in this study.
The fully atomistic Condensed-phase Optimized Molecular Potentials for Atomistic Simulation Studies (COMPASS) forcefield [55] was used to describe all intra-and intermolecular interactions.The COMPASS forcefield has the form: wherein E b , E θ , E φ , and E χ are, respectively, contributions from bond (b) stretching, angle (θ) bending, torsional angle (φ) twisting, and out-of-plane angle (χ) distortion.The terms E b,b' , E b,θ , E b,φ , and E θ,φ are cross-coupling terms between internal coordinates.E elec and E LJ are non-bonded interaction terms accounting for electrostatic and dispersive van der Waals forces, respectively.The equilibrium bond length, bond angle, and torsional angle are denoted by b 0 , θ 0 , and φ, respectively.The terms k, k 1 , k 2 , and k 3 are coefficients specific to the COMPASS forcefield.For the non-bonded interactions, q is the magnitude of the electric charge of an atom, rij is the separation between sites i and j, and ε and σ are the Lennard-Jones parameters.It may be noted that the MOF adsorbent structure was rigid in the MC simulations; hence, all terms to do with bond bending, etc., only apply to the adsorbate molecules.The COMPASS forcefield contains numerous terms to account for all of the energetic contributions described in Equation ( 1), and the reader is referred to the MCCCS Towhee project [56], which contains a significant collection of the publicly available terms for this proprietary forcefield.An analytical tail correction was used for the van der Waals interactions beyond a cut-off radius of 1.5 nm.To accurately account for the electrostatic interactions of the particles in the system, the Ewald summation technique was employed [57].In the context of the COMPASS force field, van der Waals interactions are described using a 9-6 Lennard-Jones potential [58] coupled with the Waldman-Hagler combining rules [59] for unlike interactions.The COMPASS force field has been successfully used to describe a range of metal-organic frameworks (including those with open-metal sites) [60-64] as well as the gas components of interest [65][66][67][68]; hence, it was selected for this study.It may be noted that the limitations of the molecular model include simplifying assumptions used to develop functional forms for the bonded and non-bonded interactions.
To avoid artificial periodicity [69], it may be necessary to ensure a minimum size of the simulation cell such that the length of each side is at least twice the cut-off radius.This is a general problem potentially affecting any system being simulated at the molecular level.However, it can be noted that, for some systems, periodicity artefacts may be negligible if the simulation cell is sufficiently large, such as for clathrate hydrates [70], cement hydrate [66], and montmorillonite [71], wherein the aforementioned studies demonstrated that convergence in terms of adsorption could occur for these systems despite the length of each side of a simulation cell being less than twice the cut-off radius.Hence, tests were performed using carbon dioxide as the adsorbate to assess the adsorption at the temperature of interest (298 K) as a function of fugacity for simulation cells consisting of 1 × 1 × 1, 1 × 1 × 2, and 1 × 1 × 4 unit cells.It was found that a single unit cell was sufficient to describe the CaMOF adsorbent.Additional tests were conducted to assess the effect of the cut-off radius, wherein the cut-off radius was both decreased and increased, and it was observed that any effects in terms of adsorption converged for a cut-off radius of 1.5 nm.It may be noted here that the metal-organic framework has a unit cell with parameters a = 3.18746 nm, b = 3.18746 nm, c = 0.89890 nm, α = β = 90.0• , and γ = 120 • .The results of these tests can be found in the Supplementary Materials.
A useful property in terms of assessing the practicalities of implementing the proposed adsorption system is the isosteric heat of adsorption.The isosteric heat can be considered as a function of the loading of the adsorbent to determine the regenerability of the adsorbent [72].The higher the magnitude of h ad , the higher the energy costs for regeneration of the adsorbent.The isosteric heat of adsorption was computed from the adsorption isotherm results generated via Monte Carlo molecular simulation using the Clapeyron equation: As with the adsorption isotherms, the isosteric heats of adsorption were averaged across three independent simulations at ambient temperature (i.e., T = 298 K).
In addition to Monte Carlo molecular simulations, the diffusivity of all pure gas species in CaMOF was estimated using molecular dynamics, following a procedure in the literature for MOFs [73].Dynamical simulations in the canonical NVT ensemble at 298 K were performed with periodic boundary conditions using a 1 × 1 × 4 supercell with dimensions of 3.18746 × 3.18746 × 3.55956 nm, which was packed with a number of molecules of each gas species equivalent to their absolute loading at 100 kPa to compute pure species diffusivity.This supercell had the chemical formula C 648 H 720 N 360 O 216 Ca 72 .The Nosé-Hoover thermostat [74,75] was used to maintain the temperature of the system, and the velocity Verlet algorithm [76] with a timestep of 1 fs was utilised.As per the procedure described in the literature [73], 1 ns was used to equilibrate the system, after which 5 ns was used to generate data with configurational sampling every 100 fs.The self-diffusion coefficient D was then obtained from a linear fit of the mean squared displacement (MSD) according to Einstein's relation [77] using an interval of no less than 1 ns: All pure species adsorption isotherms were assessed in terms of the most widely used two-parameter adsorption isotherm models, namely, those of Langmuir [78] and Freundlich [79].The Langmuir adsorption isotherm model is possibly the simplest means of describing gas or liquid adsorption onto a substrate and is based on the following assumptions: monolayer adsorption, no lateral interaction or steric hindrance between adsorbed molecules, and identical and equivalent adsorption sites that are finite and localised.Mathematically, the Langmuir model is typically given as q = Qb Langmuir P/(1 + b Langmuir P) (4) in which q is the quantity of gas adsorbed, P is the gas pressure, Q is the saturation adsorption capacity, and b Langmuir is the so-called Langmuir equilibrium constant.The Freundlich model accounts for reversible and non-ideal adsorption, while also factoring in the non-uniform distribution of adsorption affinities and heats across heterogeneous surfaces.It is presented as wherein K is the Freundlich constant describing the adsorption capacity and 1/n represents the degree of heterogeneity of the surface.Another two-parameter model that was considered was the Temkin adsorption isotherm model [80]: in which b Temkin accounts for the total quantity adsorbed, and A describes the pressuresensitivity of adsorption.This model ignores extreme lows and highs of pressure and assumes that the heat of adsorption of all molecules decreases linearly rather than logarithmically with respect to adsorption coverage.The Temkin model is characterised by a uniform distribution of binding energies.
In addition to the two-parameter Langmuir, Freundlich, and Temkin models of adsorption, the Redlich-Peterson [81] and Sips [82] isotherms were also considered.The Redlich-Peterson isotherm is a hybrid isotherm that incorporates features of both the Langmuir and Freundlich approaches; at low pressures, it resembles the Langmuir model, while at high pressures, it approaches the Freundlich isotherm, as given by: q = kP/(1 + aP b' ) (7) in which k/a indicates the adsorption capacity, while b' is a constant indicating deviation from the simplistic Langmuir description (for which b' = 1).
The Sips adsorption isotherm is also a combination of the Langmuir and Freundlich isotherms, which was developed to describe heterogeneous adsorption systems while circumventing the limitation of increasing adsorbate concentrations associated with the Freundlich isotherm.At low pressures, it reduces to the Freundlich isotherm, while at high pressures, it results in monolayer adsorption characteristic of the Langmuir model.The Sips model is given by: q = kP bSips /(1 + aP bSips ) (8) in which k/a indicates the adsorption capacity, while b Sips is a constant indicating deviation from the simplistic Langmuir description (for which b Sips = 1).Note that fitting was undertaken using nonlinear least squares regression, as implemented in GNU Octave version 5.2.0 [83].
Assuming ideal mixture behaviour and non-competitive sorption, a batch equilibrium model [84] of the biogas upgrading system can be applied to this context-the reader is referred to the literature [85] for a detailed discussion regarding the computational aspects of this approach.The main limitation of the batch equilibrium approach is, of course, that it assumes equilibrium to have been reached.Hence, it illustrates the best-case performance of an adsorption system.Essentially, a material balance is solved across both the charging and discharging phases of a single adsorption bed.For vacuum swing adsorption (VSA), the following material balances need to be solved simultaneously: n i ad + Wq i ad = n i de + Wq i de + n i out (10) in which n i in is the number of moles of species i entering the system, n i ad is the number of moles in the gas phase in equilibrium with the adsorbed phase Wq i ad during the pressurisation step, n i de is the number of moles in the gas phase in equilibrium with the adsorbed phase Wq i de during desorption, and n i out is the number of moles in the discharge stream.

Pure Species Adsorption Behaviour
Figure 2 shows the absolute adsorption isotherms for pure carbon dioxide, nitrogen, methane, oxygen, hydrogen sulphide, and hydrogen in the CaMOF at ambient temperature (i.e., 298 K) generated by Monte Carlo molecular simulations in the grand canonical ensemble.The results of the fitted adsorption isotherms are presented in Tables 2-6.Pure species experimental data for carbon dioxide and methane from the literature [30] are also included, although measurements for the other species are lacking.It is apparent that the adsorbent possesses a large adsorption capacity for H 2 S and CO 2 over the conditions of interest, with potentially high selectivity for H 2 S over CO 2 .This feature suggests that the CaMOF may have reduced effectiveness in terms of reducing carbon emissions for biogas containing substantial quantities of hydrogen sulphide, which may impact process design work.It is also clear that CaMOF does not adsorb H 2 significantly over the range of conditions of interest to this study, suggesting that this material may be useful for hydrogen-processing applications, at least when H 2 occurs in conjunction with biogas constituents.Regarding the present study, the low uptake of hydrogen combined with the high uptake of CO 2 and H 2 S suggests that using CaMOF to pre-treat biogas may result in a cleaner biogas product with lower carbon and sulphur emissions.It is apparent that the COMPASS forcefield consistently overestimates the CO 2 uptake by a factor of about 2.27, whereas it underestimates the adsorption of methane by a factor of about 6.03.With these empirical factors applied to the adsorption isotherms (functionally, this would have the form q' = Fq where F is the appropriate factor), the computed adsorption isotherms match the available experimental data favourably.It can be noted that despite the deviations, the shapes of the adsorption isotherm curves for both methane and carbon dioxide are qualitatively similar to the experimental data.Hence, while the assumptions behind the self-consistent COMPASS force field may have some limitations for biogas constituent interactions with this particular CaMOF, qualitative observations and, thus, physical behaviour should be consistent with experimental observations, especially given that COMPASS has been used successfully for various MOFs and biogas components in other contexts, as outlined in Section 2. The deviation present may indicate that in this specific system, corrections such as, for example, steric adjustments may need to be undertaken for this forcefield on account of the approximately constant ratio between the experimental and simulation results in terms of sorption.
whereas it underestimates the adsorption of methane by a factor of about 6.03.With these empirical factors applied to the adsorption isotherms (functionally, this would have the form q' = Fq where F is the appropriate factor), the computed adsorption isotherms match the available experimental data favourably.It can be noted that despite the deviations, the shapes of the adsorption isotherm curves for both methane and carbon dioxide are qualitatively similar to the experimental data.Hence, while the assumptions behind the selfconsistent COMPASS force field may have some limitations for biogas constituent interactions with this particular CaMOF, qualitative observations and, thus, physical behaviour should be consistent with experimental observations, especially given that COMPASS has been used successfully for various MOFs and biogas components in other contexts, as outlined in Section 2. The deviation present may indicate that in this specific system, corrections such as, for example, steric adjustments may need to be undertaken for this forcefield on account of the approximately constant ratio between the experimental and simulation results in terms of sorption.With respect to carbon dioxide, a comparison between the adsorption isotherm models with the strongest correlation to experimental data (i.e., the Langmuir and Sips models) may be undertaken using the Akaike information criterion, AIC [86], adjusted for small sample sizes [87] (AICc).In this case, the most likely model to describe the results of the simulations would be the Langmuir model, given that it has an AICc of 14.356 compared to 36.764 for the Sips model.These values suggest that the Sips model is approximately 1.4 × 10 −5 times as probable as the Langmuir model.For methane, it is clear that the best fit among the adsorption isotherm models is the Langmuir model, since its correlation coefficient is practically unity.Experimental measurements [30] for methane and carbon dioxide over much narrower pressure ranges yielded behaviour incorporating features of both the Langmuir and Freundlich adsorption isotherm models.It can be noted that in this study, strictly employing the correlation coefficient as a measure of the goodness of fit would result in the Sips model being selected over the Langmuir model on account of its higher R 2 value; however, the Akaike information criterion suggests that it may, in fact, be overfitting the data, as the Sips model is only about 0.013 times as probable as the Langmuir model.An outcome of these analyses for carbon dioxide and methane is that the adsorption sites are likely homogeneous, since this is a key assumption underlying the Langmuir adsorption isotherm model.
In terms of fitting adsorption isotherm models to the data for H 2 S, it may be preferable to select either the Redlich-Peterson or Sips models.Purely on the basis of the R 2 as the value for comparison, the better choice would be the Sips model.However, computing the AICc for these adsorption isotherm models yields the following values for the Redlich-Peterson and Sips models, respectively: 40.5, and 41.7.Hence, the optimal model to select in terms of both agreeing with the results and minimising overfitting is the Redlich-Peterson model.On the basis of AICc, the Sips model is 0.54 times as probable as the Redlich-Peterson model.This suggests that the adsorption of hydrogen sulphide on the CaMOF adsorbent has features of both the Langmuir and Freundlich isotherms, tending to Langmuir-type behaviour at low pressure with greater site heterogeneity at high pressures as it shifts towards Freundlich-type behaviour.In terms of diffusion behaviour, it was difficult to discern trends in terms of the molecular properties of the gas species and the estimated diffusivity.When considering the molecular mass of the gas species, it was observed that the diffusivity generally decreased with increasing gas molecule mass, which is largely expected based on the kinetic behaviour of fluid species.This trend was not linear, however, and there may be steric effects related to the size and shape of gas molecules as well as the size and shape of channels within the CaMOF.Comparing the magnitudes of the diffusivities yields the following ranking in terms of rates of diffusion in the MOF: H 2 > N 2 > O 2 > CO 2 > CH 4 > H 2 S. The ranking of affinities in terms of adsorption (considering pressures from 10 kPa upwards) is H 2 S > CO 2 > O 2 > N 2 > CH 4 > H 2 .Comparing these two rankings as well as the results in Figure 3 suggests that it may not be straightforward to predict the diffusion of biogas constituents in CaMOF based solely on the quantity of gas adsorbed or a simplistic measure such as a single molecular property.

Isosteric Heat of Adsorption
Figure 4 shows the heats of adsorption for all species computed from the results of the pure species Monte Carlo molecular simulations.It is apparent that the CaMOF adsorbent is not well suited to remove CO2 as the energy costs of regeneration tend to increase as the quantity of carbon dioxide adsorbed increases, whereas the heat of adsorp-

Isosteric Heat of Adsorption
Figure 4 shows the heats of adsorption for all species computed from the results of the pure species Monte Carlo molecular simulations.It is apparent that the CaMOF adsorbent is not well suited to remove CO 2 as the energy costs of regeneration tend to increase as the quantity of carbon dioxide adsorbed increases, whereas the heat of adsorption of nitrogen does not vary with respect to the amount of gas adsorbed.While the trend for hydrogen sulphide is that of a rapid decrease in terms of the heat of adsorption with respect to gas uptake, the magnitude of h ad is rather high at 57 to 68 kJ/mol over the range of conditions considered in this study.The energy costs of regeneration in this study refer to the amount of energy needed to release the adsorbed gas molecules from CaMOF adsorbent regenerating the adsorbent for reuse.One of the key factors influencing these energy costs is the heat of adsorption.Higher h ad leads to stronger interactions between the adsorbent and the adsorbed gas molecules.This means that more energy will be required to release the gas molecules during regeneration processes.For comparison, the heat of adsorption of CO 2 varied from 29 to 40 kJ/mol in this study, while a previous report [72] on a favourable adsorbent in terms of regenerability yielded an isosteric heat of adsorption of approximately 20 kJ/mol for carbon dioxide on a different MOF-based adsorbent.The increase in h ad with the increase in the amount of CO 2 adsorbed implies that as the CO 2 is adsorbed, it becomes difficult and energetically expensive to desorb it.This also applies to the adsorption of H 2 S, because the h ad of H 2 S is high, which implies that regenerating adsorbent for H 2 S removal will be energy-intensive.The adsorption heat effect impacts not only the energy required for regenerating the adsorbent, but also the initial adsorption capacity and efficiency.As seen in the case of H 2 S, the heat of adsorption decreases with the increase in gas uptake, which may cause the adsorbent to become less efficient at higher gas concentrations.Hence, it can be expected that regenerability issues may be encountered when employing CaMOF in industrial operations, and overcoming this problem may increase the energetic costs of the gas separation process.

Mixture Sorption
Figure 5 shows the uptake of all biogas constituents on CaMOF for a typical biogas composed of 30 mol-% CO 2 , 14.99 mol-% N 2 , 50 mol-% CH 4 , 3 mol-% O 2 , 0.01 mol-% H 2 S, and 2 mol-% H 2 [5] in terms of absolute adsorption.It is apparent that for the enhanced selectivity of hydrogen sulphide and carbon dioxide, the vacuum pressure operation of pressure swing adsorption systems would be preferred.In terms of non-calorific species selectivity-defined as S = (Σx non-calorific /Σy non-calorific )/(Σx calorific /Σy calorific ), where x is the mole fraction of the relevant species in the adsorbed phase-optimal selectivity (S = 121) to remove the non-calorific species was observed at 1000 kPa, with a decreasing trend both below and above this pressure.However, the selectivity in favour of the adsorption of non-calorific species was higher at 10 kPa (where S = 98) than at 10 4 kPa.It can be noted that, in terms of the selectivity of a specific species i with respect to CH 4 , the expression for S would be defined as S = (x i /y i )/(x CH4 /y CH4 ).An additional biogas feedstock was also considered, as shown in Figure 6.This feed gas had a composition of 30 mol-% CO2, 14.9 mol-% N2, 50 mol-% CH4, 3 mol-% O2, 0.1 mol-% H2S, and 2 mol-% H2 [5], therefore possessing a higher proportion of sour gas constituents.The selectivity trends were similar to the previous biogas feedstock, with a peak selectivity in favour of non-calorific species of 128 at 1000 kPa, which decreased to 87 by 10 4 kPa and to 111 at 10 kPa.It is apparent then that a higher fraction of hydrogen sulphide in the biogas may make the performance of biogas cleaning via CaMOF less sensitive to pressure changes, especially below 1000 kPa.Altogether, both sets of multicomponent molecular simulations suggest that low-pressure or vacuum-pressure operation is desirable for the purposes of cleaning biogas using CaMOF on account of the high selectivity observed for non-calorific species, as well as the relatively low sensitivity to pressure changes and the magnitude of the uptake of the dominant species (i.e., methane and carbon dioxide).An additional biogas feedstock was also considered, as shown in Figure 6.This feed gas had a composition of 30 mol-% CO 2 , 14.9 mol-% N 2 , 50 mol-% CH 4 , 3 mol-% O 2 , 0.1 mol-% H 2 S, and 2 mol-% H 2 [5], therefore possessing a higher proportion of sour gas constituents.The selectivity trends were similar to the previous biogas feedstock, with a peak selectivity in favour of non-calorific species of 128 at 1000 kPa, which decreased to 87 by 10 4 kPa and to 111 at 10 kPa.It is apparent then that a higher fraction of hydrogen sulphide in the biogas may make the performance of biogas cleaning via CaMOF less sensitive to pressure changes, especially below 1000 kPa.Altogether, both sets of multicomponent molecular simulations suggest that low-pressure or vacuum-pressure operation is desirable for the purposes of cleaning biogas using CaMOF on account of the high selectivity observed for non-calorific species, as well as the relatively low sensitivity to pressure changes and the magnitude of the uptake of the dominant species (i.e., methane and carbon dioxide).The observed selectivity behaviour at peak adsorbent loading can be compared with several other adsorbents proposed for biogas purification, as outlined in Table 7.It is evident that CaMOF shows some potential for separating both H2S and CO2 from methane on account of it possessing substantial selectivity for hydrogen sulphide compared to other adsorbents that have been proposed and studied for this task.As stated previously, a significant limitation of the present study is its in silico nature, although the results suggest that proceeding with laboratory experiments may be a fruitful venture.However, the strength of the forces between H2S and the CaMOF, especially at low pressures, should be considered.In this regard, Figures 4 and 5 show that the adsorption sites of hydrogen sulphide at low pressure are strongly attractive, much more so than for CO2.This suggests that the regeneration of adsorption beds may be challenging when such systems are employed to remove hydrogen sulphide from biogas.The observed selectivity behaviour at peak adsorbent loading can be compared with several other adsorbents proposed for biogas purification, as outlined in Table 7.It is evident that CaMOF shows some potential for separating both H 2 S and CO 2 from methane on account of it possessing substantial selectivity for hydrogen sulphide compared to other adsorbents that have been proposed and studied for this task.As stated previously, a significant limitation of the present study is its in silico nature, although the results suggest that proceeding with laboratory experiments may be a fruitful venture.However, the strength of the forces between H 2 S and the CaMOF, especially at low pressures, should be considered.In this regard, Figures 4 and 5 show that the adsorption sites of hydrogen sulphide at low pressure are strongly attractive, much more so than for CO 2 .This suggests that the regeneration of adsorption beds may be challenging when such systems are employed to remove hydrogen sulphide from biogas.

Discussion
Firstly, the hydrogen sulphide sorption behaviour should be discussed.Figure 7 shows the density field of H 2 S (at 298 K and the highest pressure considered for hydrogen sulphide in this study, namely, 99.22 kPa) in the CaMOF lattice.It is apparent that there may be multiple types of hydrogen sulphide adsorption sites as there are sites located near Ca, N, and O atoms.This observation alone suggests that an adsorption isotherm model accounting for site heterogeneity such as the Freundlich model, or possibly a model incorporating its features such as the Redlich-Peterson model, would most accurately describe hydrogen sulphide adsorption in the CaMOF adsorbent.It may be noted, however, that it may be more likely that there are adsorption sites at the channel vertices where Ca and N atoms are located.The adsorption isotherm analysis presented in Section 3.1 would concur with this visual assessment, as it was observed that the best fit for H 2 S adsorption was obtained using the Redlich-Peterson mode, which has features of both the Langmuir and Freundlich isotherms, tending to Langmuir-type behaviour at low pressure with greater site heterogeneity at high pressures as it shifts towards Freundlich-type behaviour.An examination of the energy distribution for H 2 S adsorbed in CaMOF, as shown in Figure 7, may be instructive in this regard.
Figure 8 shows three adsorption site types for H 2 S defined in terms of the potential energy associated with them: −75 to −73 kJ/mol, −50 to −47 kJ/mol, and −32 to −30 kJ/mol.The latter two adsorption site types appear to be hard to distinguish at lower pressures, as evidenced by the significant degree of overlap between the energy distribution curves shown in Figure 4; this difficulty in distinguishing between the two higher energy sites is evident in the density field of adsorbed H 2 S molecules shown in Figure 4.The first of the identified adsorption site types (that are associated with a potential energy of −75 to −73 kJ/mol) may be related to the previously identified Ca/N regions of high affinity for H 2 S sorption.The change in the probability density of the different peaks shown in Figure 5 demonstrates the difficulties in determining the adsorption behaviour in terms of simplistic adsorption isotherm models such as the Langmuir and Freundlich models.The results suggest that, at low pressures (below 18 kPa), there may effectively be only one accessible type of adsorption site for H 2 S molecules, whereas at higher pressures, the shift is from one type of adsorption site to two site types at intermediate pressures (as shown by the clear distinction between the peaks at −75 to −73 kJ/mol and −50 to −47 kJ/mol at pressures from 18 kPa up to about 99 kPa), and then towards three site types from about 18 kPa onwards as the third energetic site type (associated with the peak at −32 to −30 kJ/mol) becomes noticeable.In terms of what this means in practical terms, if the modification or pre-treatment of the CaMOF is undertaken, the H 2 S sorption capacity of the adsorbent at predetermined pressure values can be tailored to some extent by targeting the different adsorption sites.Figure 8 shows three adsorption site types for H2S defined in terms of the potential energy associated with them: −75 to −73 kJ/mol, −50 to −47 kJ/mol, and −32 to −30 kJ/mol.The latter two adsorption site types appear to be hard to distinguish at lower pressures, as evidenced by the significant degree of overlap between the energy distribution curves shown in Figure 4; this difficulty in distinguishing between the two higher energy sites is evident in the density field of adsorbed H2S molecules shown in Figure 4.The first of the identified adsorption site types (that are associated with a potential energy of −75 to −73 kJ/mol) may be related to the previously identified Ca/N regions of high affinity for H2S sorption.The change in the probability density of the different peaks shown in Figure 5 demonstrates the difficulties in determining the adsorption behaviour in terms of simplistic adsorption isotherm models such as the Langmuir and Freundlich models.The results suggest that, at low pressures (below 18 kPa), there may effectively be only one accessible type of adsorption site for H2S molecules, whereas at higher pressures, the shift is from one type of adsorption site to two site types at intermediate pressures (as shown by the clear distinction between the peaks at −75 to −73 kJ/mol and −50 to −47 kJ/mol at pressures from 18 kPa up to about 99 kPa), and then towards three site types from about 18 kPa onwards as the third energetic site type (associated with the peak at −32 to −30 kJ/mol) becomes noticeable.In terms of what this means in practical terms, if the modification or pre-treatment of the CaMOF is undertaken, the H2S sorption capacity of the adsorbent at predetermined pressure values can be tailored to some extent by targeting the different For comparison, the energy distributions of the next two most strongly adsorbed species, CO 2 and CH 4 , are shown in Figures 9 and 10.For both species, it is apparent that increasing the gas pressure shifts the mean potential energy of the system, but in all cases, there remains a single energetic site type.Therefore, regardless of the pressure of the gas reservoir, the adsorption sites for a given pressure are homogeneous.This is in contrast to the case for hydrogen sulphide, where up to three energetic site types were observed.In terms of tailoring the CaMOF for targeted adsorption of carbon dioxide or methane, it is thus apparent that a single energetic site type needs to be considered, in contrast to the three site types observed for H 2 S. It is suggested that work on the selfdiffusion or binary diffusion of biogas constituents may be undertaken to probe gas-gas and gas-solid interactions from a dynamical perspective, which may complement the equilibrium-focussed work of the present study.For comparison, the energy distributions of the next two most strongly adsorbed species, CO2 and CH4, are shown in Figures 9 and 10.For both species, it is apparent that increasing the gas pressure shifts the mean potential energy of the system, but in all cases, there remains a single energetic site type.Therefore, regardless of the pressure of the gas reservoir, the adsorption sites for a given pressure are homogeneous.This is in contrast to the case for hydrogen sulphide, where up to three energetic site types were observed.In terms of tailoring the CaMOF for targeted adsorption of carbon dioxide or methane, it is thus apparent that a single energetic site type needs to be considered, in contrast to the three site types observed for H2S.It is suggested that work on the self-diffusion or binary diffusion of biogas constituents may be undertaken to probe gas-gas and gas-solid interactions from a dynamical perspective, which may complement the equilibrium-focussed work of the present study.In terms of potential design work, the results of the Monte Carlo simulations suggest that vacuum swing adsorption using CaMOF adsorbent may be pursued to yield a practicable gas separation system.Vacuum operation does, of course, entail its own challenges, such as higher capital costs for implementation, and a techno-economic assessment that balances capital and operating costs to determine the optimal vacuum pressure would be necessary [95].That being said, a favourable selectivity in favour of non-calorific versus calorific gases at 10 kPa suggests a promising pathway to developing an effective biogas upgrading system.For example, a quick estimate of the mass of adsorbent containing all of the hydrogen sulphide in the biogas mixture referred to in Figure 7 yields quantities of about 0.05 to 0.06 g of adsorbent per mole of feed gas.Hence, producing cleaner-burning biogas without sulphur may require low capital costs by virtue of smaller adsorbent bed sizes.With reference to CO 2 removal, it is also then clear that very small bed sizes may be sufficient to produce purified biogas with a higher heating value.
Apart from upgrading the calorific value of biogas for use as a fuel, CaMOF may be incorporated into biogas reforming processes for the removal of carbon dioxide from feedstock as part of biogas upgrading; for the purposes of this study, solar-driven biogas steam reforming [96] will be used as a reference case study due to its recent active development.As shown in Figure 11, a key step in this process is biogas upgrading wherein a feedstock comprising 60 mol-% methane with carbon dioxide is separated to produce pure carbon dioxide and a stream rich in methane with a carbon dioxide residue.Adsorption typically entails low energy costs compared to many separation schemes, and it can therefore be helpful to consider a CaMOF-based pressure swing adsorption system for the purpose of biogas upgrading.
Note that W is the mass of the adsorption bed.Due to the selective adsorption of carbon dioxide by CaMOF, the gas discharged by the pressure swing adsorption unit described by Equations ( 8) and ( 9) would be rich in methane.Solving the simultaneous material balance over a range of conditions yields the results presented in Figure 12.To replicate the process conditions in the literature [96] as closely as possible for the purposes of integration with the solar biogas reforming process, a bed mass of about 30-42 g/mole of feed gas with an operating pressure of approximately 70-78 kPa would be suitable for cleaning a gas mixture consisting of 60 mol-% methane with carbon dioxide.These conditions can yield a methane-rich stream comprising 93 mol-% methane and all of the methane that was in the feed, demonstrating that a CaMOF-based VSA system may be incorporated into a solardriven biogas reforming process.It can be noted that the secondary outlet stream consists solely of carbon dioxide, which fits with the requirements of the solar biogas reforming process, wherein the biogas upgrading step is intended to produce two gas streams: one containing only CO 2 , and another comprising about 93-94 mol-% CH 4 .Therefore, the proposed VSA employing CaMOF as the adsorbent may satisfy the requirements of this biogas reforming process.From an operational viewpoint, it is apparent that for a given bed mass, the vacuum pressure should be kept within a band of approximately 4 kPa from the setpoint.In terms of bed degradation such that the effective mass of accessibly adsorbent decreases, Figure 12 shows that the VSA process is less sensitive than it is to operating pressure fluctuations, as evidenced by the broader range of bed mass values for which the target concentration of 93-94 mol-% methane can be achieved.It can be noted that test calculations were performed for operations approaching 1000 kPa; however, the magnitude of methane uptake into CaMOF meant that achieving both outlet composition specifications-i.e., one gas stream with 100 mol-% CO 2 and another with 93-94 mol-% CH 4 -would be challenging in practice when employing a single-stage pressure swing adsorption system.In terms of potential design work, the results of the Monte Carlo simulations suggest that vacuum swing adsorption using CaMOF adsorbent may be pursued to yield a practicable gas separation system.Vacuum operation does, of course, entail its own challenges, such as higher capital costs for implementation, and a techno-economic assessment that balances capital and operating costs to determine the optimal vacuum pressure would be necessary [95].That being said, a favourable selectivity in favour of non-calorific versus calorific gases at 10 kPa suggests a promising pathway to developing an effective biogas upgrading system.For example, a quick estimate of the mass of adsorbent containing all of the hydrogen sulphide in the biogas mixture referred to in Figure 7 yields quantities of about 0.05 to 0.06 g of adsorbent per mole of feed gas.Hence, producing cleaner-burning biogas without sulphur may require low capital costs by virtue of smaller adsorbent bed sizes.With reference to CO2 removal, it is also then clear that very small bed sizes may be sufficient to produce purified biogas with a higher heating value.
Apart from upgrading the calorific value of biogas for use as a fuel, CaMOF may be incorporated into biogas reforming processes for the removal of carbon dioxide from feedstock as part of biogas upgrading; for the purposes of this study, solar-driven biogas steam reforming [96] will be used as a reference case study due to its recent active development.As shown in Figure 11, a key step in this process is biogas upgrading wherein a feedstock comprising 60 mol-% methane with carbon dioxide is separated to produce pure carbon dioxide and a stream rich in methane with a carbon dioxide residue.Adsorption typically entails low energy costs compared to many separation schemes, and it can therefore be helpful to consider a CaMOF-based pressure swing adsorption system for the purpose of biogas upgrading.Note that W is the mass of the adsorption bed.Due to the selective adsorption of carbon dioxide by CaMOF, the gas discharged by the pressure swing adsorption unit described by Equations ( 8) and ( 9) would be rich in methane.Solving the simultaneous material balance over a range of conditions yields the results presented in Figure 12.To replicate the process conditions in the literature [96] as closely as possible for the purposes of integration with the solar biogas reforming process, a bed mass of about 30-42 g/mole of feed gas with an operating pressure of approximately 70-78 kPa would be suitable for cleaning a gas mixture consisting of 60 mol-% methane with carbon dioxide.These If one considers the correction factors mentioned in Section 3.1 to adjust the quantity of CO 2 and CH 4 adsorbed, then the batch equilibrium modelling yields a qualitatively similar result, as shown in Figure 13.Here, it is apparent that a smaller bed size over a similar operating range yields performance such that the targeted methane content of 93-94 mol-% is achievable.The required bed size in this case is about 5-6 g/mole of feedstock.Therefore, when considering this correction factor to account for consistent deviations between the available experimental data and the predictions of the COMPASS forcefield, it is apparent that CaMOF can serve as a suitable adsorbent for the solar biogas reforming process.
achieved.It can be noted that test calculations were performed for operations approachi 1000 kPa; however, the magnitude of methane uptake into CaMOF meant that achievi both outlet composition specifications-i.e., one gas stream with 100 mol-% CO2 and a other with 93-94 mol-% CH4-would be challenging in practice when employing a sing stage pressure swing adsorption system.If one considers the correction factors mentioned in Section 3.1 to adjust the quant of CO2 and CH4 adsorbed, then the batch equilibrium modelling yields a qualitative similar result, as shown in Figure 13.Here, it is apparent that a smaller bed size ove similar operating range yields performance such that the targeted methane content of 9 94 mol-% is achievable.The required bed size in this case is about 5-6 g/mole of feedstoc Therefore, when considering this correction factor to account for consistent deviations tween the available experimental data and the predictions of the COMPASS forcefield is apparent that CaMOF can serve as a suitable adsorbent for the solar biogas reformi process.Regarding the outlook, this study integrated atomistic simulation with macro-sc batch equilibrium modelling with a view to integrating a gas separation unit with an isting chemical process design.This workflow can be employed in future work to scre multiple adsorbent materials to determine their suitability for integration into existi Regarding the outlook, this study integrated atomistic simulation with macro-scale batch equilibrium modelling with a view to integrating a gas separation unit with an existing chemical process design.This workflow can be employed in future work to screen multiple adsorbent materials to determine their suitability for integration into existing processes.In principle, new adsorbents can be designed in silico and then their potential for industrial application can be assessed immediately.The results obtained by such multiscale modelling are largely limited by the assumptions that go into interparticle potential energy modelling at the atomistic scale, although, in many cases, the fundamental behaviour of physisorption processes is well described by modern forcefields, and such computation screening can suggest pathways for future gas separation research.An advantage of this multiscale modelling approach is that highly variable feedstocks can be considered, which is a feature of biomass-derived products.Hence, this modelling approach can assist in rapidly adjusting existing processes or designing new ones to take advantage of the features of novel or inconsistent biogas or biomass-derived feedstocks.

Conclusions
Atomistic Monte Carlo simulations in the grand canonical ensemble were employed to study the sorption of biogas constituents on CaMOF to assess, in silico, its potential usefulness to purify biogas via sorption processes such as pressure or vacuum swing adsorption.The hydrogen sulphide displayed unusual behaviour in that it had two or three different adsorption site types within the CaMOF (based on an assessment of the potential energy of the simulation cell), complicating any description of hydrogen sulphide adsorption on CaMOF in terms of common adsorption isotherm models.With respect to the other species under consideration (i.e., nitrogen, oxygen, methane, hydrogen, and carbon dioxide), commonly used adsorption isotherm models adequately described the observed sorption behaviour.
Simulations of mixture adsorption were also undertaken, for which it was observed that the optimal removal of non-calorific components occurred at vacuum pressure (specifically at 10 kPa, the lowest pressure considered in this study).The strong selectivity of the CaMOF adsorbent in favour of hydrogen sulphide at low pressures suggests that future experimental work may yield a practicable vacuum swing adsorption system to produce cleaner-burning biogas free from sulphurous impurities and with reduced non-calorific content.In addition, low-sulphur biogas cleaned using vacuum swing adsorption in this way may be beneficial to steam reforming operations using renewable biogas feedstocks.
In addition to molecular simulations, batch equilibrium modelling was undertaken with a view to assessing the feasibility of the CaMOF as an adsorbent for the vacuum swing adsorption separation of a CH 4 /CO 2 mixture.This was undertaken within the context of solar biogas reforming, where the effects of operating pressure and bed size for a singlebed adsorber were examined.It was demonstrated that this unit may be incorporated into a solar biogas reforming process, as it can meet the requirements of the biogas feed preparation step in terms of separation power.A downside to consider, however, is that there may be challenges with adsorbent regenerability to be overcome, as indicated by an assessment of the isosteric heat of adsorption of H 2 S.

ChemEngineering 2024, 8 , 26 Figure 1 .
Figure 1.Unit cell of the CaMOF used in this study.Atoms are colour-coded according to the element: white = hydrogen, grey = carbon, red = oxygen, blue = nitrogen, and green = calcium.The MOF in question has a unit cell with parameters a = 3.18746 nm, b = 3.18746 nm, c = 0.89890 nm, α = β = 90.0°,and γ = 120° and the chemical formula C162H180N90O54Ca18.Coordinatively unsaturated metal sites, known as open-metal sites, are a common feature of MOF materials and are a key characteristic in terms of enhanced[31][32][33][34][35] and selective[36][37][38][39] gas adsorption at low pressures or fugacities.These sites may be unsaturated and made accessible via framework defects such as missing organic linkers[40][41][42] or, commonly, through chemical[43,44] or thermal[45][46][47][48][49] activation protocols.The CaMOF considered in this study has a honeycomb porous framework with regular 1D channels.The Ca 2+ cations have a pentagonal bipyramid coordinated geometry and are connected to one oxygen atom in DMA, three nitrogen atoms in the dtztp ligand, and three carboxylate oxygen atoms in the dtztp ligand.The aforementioned channels are lined with exposed Ca 2+ centres along with uncoordinated nitrogen atoms such that the effective window diameter is 1.5 nm and there is an accessible voidage of 68.9%[30].The aim of this study is to assess the feasibility of the newly developed CaMOF to separate undesirable sour components from biogas, specifically CO2 and H2S.Multicomponent molecular simulations of biogas blended with this adsorbent have not been published in the open literature.The batch equilibrium modelling of a pressure swing adsorption system using this adsorbent, another topic of analysis in this study, has also not been published in the open literature.This study also presents the first investigation into the incorporation of such a pressure swing adsorption system using the CaMOF into a biogas reforming process.

Figure 1 .
Figure 1.Unit cell of the CaMOF used in this study.Atoms are colour-coded according to the element: white = hydrogen, grey = carbon, red = oxygen, blue = nitrogen, and green = calcium.

Figure 2 .
Figure 2. Pure species adsorption isotherms for biogas constituents in CaMOF.Where they are not visible, the error bars (representing the standard deviation) are smaller than the symbols.The lines are guides for the eye.The open symbols represent experimental data [30] and the dashed line series represents an adjustment of the adsorption isotherms using a simple empirical factor, as outlined in the text.

Figure 3 .
Figure 3. Diffusivity plotted as a function of the molar mass of each gas species.

Figure 3 .
Figure 3. Diffusivity plotted as a function of the molar mass of each gas species.

Figure 4 .
Figure 4. Pure species isosteric heats of adsorption (h ad ) in the Ca-MOF adsorbent.

Figure 5 .
Figure 5. Mixture adsorption isotherms for biogas constituents in Ca-MOF for biogas composed of 30 mol-% CO 2 , 14.99 mol-% N 2 , 50 mol-% CH 4 , 3 mol-% O 2 , 0.01 mol-% H 2 S, and 2 mol-% H 2 .Where they are not visible, the error bars (representing the standard deviation) are smaller than the symbols.The lines are guides for the eye.

Figure 6 .
Figure 6.Mixture adsorption isotherms for biogas constituents in Ca-MOF for biogas composed of 30 mol-% CO 2 14.9 mol-% N 2 , 50 mol-% CH 4 , 3 mol-% O 2 , 0.1 mol-% H 2 S, and 2 mol-% H 2 .Where they are not visible, the error bars (representing the standard deviation) are smaller than the symbols.The lines are guides for the eye.

ChemEngineering 2024, 8 , 26 Figure 7 .
Figure 7. Density field of H2S adsorbed onto the unit cell of the Ca-MOF at 298 K and 99.22 kPa.Atoms are colour-coded according to the element: white = hydrogen, grey = carbon, red = oxygen, blue = nitrogen, and green = calcium.Hydrogen sulphide-containing regions are indicated in black.

Figure 7 .
Figure 7. Density field of H 2 S adsorbed onto the unit cell of the Ca-MOF at 298 K and 99.22 kPa.Atoms are colour-coded according to the element: white = hydrogen, grey = carbon, red = oxygen, blue = nitrogen, and green = calcium.Hydrogen sulphide-containing regions are indicated in black.

Figure 8 .
Figure 8. Energy distribution for H2S adsorbed in CaMOF.E is the potential energy derived from configurations sampled during the grand canonical Monte Carlo simulations, and p(E) is the probability of occurrence for each energy value.The lines are guides for the eye.

Figure 8 .
Figure 8. Energy distribution for H 2 S adsorbed in CaMOF.E is the potential energy derived from configurations sampled during the grand canonical Monte Carlo simulations, and p(E) is the probability of occurrence for each energy value.The lines are guides for the eye.

ChemEngineering 2024, 8 , 26 Figure 9 .
Figure 9. Energy distribution for CO2 adsorbed in CaMOF.E is the potential energy derived from configurations sampled during the grand canonical Monte Carlo simulations, and p(E) is the probability of occurrence for each energy value.The lines are guides for the eye.

Figure 9 .
Figure 9. Energy distribution for CO 2 adsorbed in CaMOF.E is the potential energy derived from configurations sampled during the grand canonical Monte Carlo simulations, and p(E) is the probability of occurrence for each energy value.The lines are guides for the eye.

Figure 9 .
Figure 9. Energy distribution for CO2 adsorbed in CaMOF.E is the potential energy derived from configurations sampled during the grand canonical Monte Carlo simulations, and p(E) is the probability of occurrence for each energy value.The lines are guides for the eye.

Figure 10 .
Figure 10.Energy distribution for CH4 adsorbed in CaMOF.E is the potential energy derived from configurations sampled during the grand canonical Monte Carlo simulations, and p(E) is the probability of occurrence for each energy value.The lines are guides for the eye.

Figure 10 .
Figure 10.Energy distribution for CH 4 adsorbed in CaMOF.E is the potential energy derived from configurations sampled during the grand canonical Monte Carlo simulations, and p(E) is the probability of occurrence for each energy value.The lines are guides for the eye.

Figure 12 .
Figure 12.Heat map of CO2 composition (in mol-%) in the outlet stream from the VSA system a function of bed mass and operating pressure ratio.

Figure 12 .
Figure 12.Heat map of CO 2 composition (in mol-%) in the outlet stream from the VSA system as a function of bed mass and operating pressure ratio.

Figure 13 .
Figure 13.Heat map of CO2 composition (in mol-%) in the outlet stream from the VSA system a function of bed mass and operating pressure ratio, when a correction factor is incorporated to count for deviations between experiments and COMPASS forcefield predictions.Details are p vided in the text in Section 3.1.

Figure 13 .
Figure 13.Heat map of CO 2 composition (in mol-%) in the outlet stream from the VSA system as a function of bed mass and operating pressure ratio, when a correction factor is incorporated to account for deviations between experiments and COMPASS forcefield predictions.Details are provided in the text in Section 3.1.

:
Effect of simulation cell size on adsorption of carbon dioxide in CaMOF at 298 K; Figure S2: Effect of van der Waals cut-off radius on adsorption of carbon dioxide in CaMOF at 298 K. Author Contributions: Conceptualisation, M.L. and K.T.; methodology, M.L. and V.T.A.; software, M.L.; validation, M.L. and K.T.; formal analysis, M.L., V.T.A. and K.T.; investigation, M.L. and V.T.A.; resources, M.L.; data curation, M.L.; writing-original draft preparation, M.L.; writing-review and editing, M.L., V.T.A. and K.T.; visualisation, M.L.; supervision, M.L. and K.T.; project administration, M.L.All authors have read and agreed to the published version of the manuscript.Funding: This research received no external funding.

Table 3 .
Fitted parameters and correlation coefficient (R 2 ) for the Freundlich model based on the results of grand canonical Monte Carlo simulations at T = 298 K.

Table 4 .
Fitted parameters and correlation coefficient (R 2 ) for the Redlich-Peterson model based on the results of grand canonical Monte Carlo simulations at T = 298 K.

Table 5 .
Fitted parameters and correlation coefficient (R 2 ) for the Temkin model based on the results of grand canonical Monte Carlo simulations at T = 298 K.

Table 6 .
Fitted parameters and correlation coefficient (R 2 ) for the Sips model based on the results of grand canonical Monte Carlo simulations at T = 298 K.

Table 7 .
Comparison of CaMOF selectivity with other adsorbent materials for the purpose of biogas purification.* These values are the selectivity factor qCO2/qCH4.

Table 7 .
Comparison of CaMOF selectivity with other adsorbent materials for the purpose of biogas purification.* These values are the selectivity factor q CO2 /q CH4 .