Parameter Estimation of Breakthrough Curve Models in the Adsorption Process of H 2 S and CO 2 Using the Markov Chain Monte Carlo Method

: The increase in emissions of toxic gasses such as hydrogen sulfide (H 2 S) and carbon dioxide (CO 2 ), resulting from growing urbanization and industrialization, has caused environmental and public health problems, making the implementation of air purification techniques through adsorption important. Thus, modeling the gas adsorption process is fundamental for good agreement with experimental data, employing mathematical models that enable the prediction of adsorption capacity. In this way, the present work aimed to compare different analytical breakthrough curve models (Thomas, Yoon–Nelson, Adams–Bohart, and Yan) for the adsorption of H 2 S and CO 2 in fixed-bed columns, using experimental data from the literature, estimating the curve parameters through the Markov Chain Monte Carlo (MCMC) method with the Metropolis–Hastings algorithm, and ranking using the determination coefficients (R 2 and R 2Adjusted ) and the Bayesian Information Criterion (BIC). The models showed better agreement using the estimation of maximum adsorption capacity (q s , N 0 ) and the constants of each model (k th , k yn , and k ba ). In the adsorption of H 2 S, the Yan model stood out for its precision in estimating q s . For the adsorption of CO 2 , the Adams–Bohart model achieved better results with the estimation of N 0 , along with the Yoon–Nelson model. Furthermore, the use of this method allows for a reduction in computational effort compared to models based on complex differential equations.


Introduction
Increasing urbanization and industrialization worldwide have heightened concerns about air quality due to rising emissions of atmospheric pollutants.With technological advancements and the expansion of industrial activities, the release of substances such as hydrogen sulfide or hydrogen sulfide gas (H 2 S), nitric oxide (NO), nitrogen dioxide (NO 2 ), carbon dioxide (CO 2 ), and sulfur dioxide (SO 2 ) has become more frequent.These pollutants are known for their adverse impact on both the environment and human health [1].
Hydrogen sulfide is released primarily during industrial processes such as petroleum refineries and water and sewage treatment plants, and it occurs naturally in biological processes such as organic matter decomposition in waste deposits and wetlands.Furthermore, it is a substance known for its toxicity and can pose significant risks to human health when present in high concentrations [2].
On the other hand, carbon dioxide (CO 2 ), a colorless, odorless, and non-flammable gas, is produced on a large scale by the industry, mainly as a byproduct in the production of ammonia and hydrogen, as well as in the burning of fossil fuels and petroleum refining.Although it can be used in some chemical reactions, CO 2 is known as the primary greenhouse gas, contributing significantly to climate change and global warming, becoming a globally recognized environmental concern [3][4][5].
In this context, gas adsorption has therefore become an area of great interest due to the need to develop more effective air purification systems.This process involves the adherence of atoms, molecules, or ions of substances (adsorbates) to a solid surface (adsorbent), resulting in the formation of a molecular or atomic layer on this surface.Adsorption has various applications, including water purification, separation of chemical compounds, and especially the capture of polluting gasses [6][7][8][9].
The use of fixed-bed columns for adsorption is preferred in the industry due to their ability to treat large quantities of pollutants in a continuous flow system, enabling their implementation on a large scale.These columns provide breakthrough curves that describe the variation in pollutant concentration over time.Therefore, mathematical modeling of this system is important for understanding, analyzing, and optimizing separation and purification processes [10][11][12].
However, there are challenges in modeling continuous adsorption systems due to their non-linearity.Parameter estimation of the models often relies on linear regression techniques, including the method of least squares; however, such methods can distort parameter interpretation and compromise model assumptions [13].
Bayesian inference emerges as a promising alternative to least squares fitting, surpassing its limitations.This technique provides a method to combine available experimental data with mathematical modeling and the prior probability distribution of parameters.In general, it involves a statistical analysis of the posterior probability density, which represents the conditional probability of the parameters based on observations.The posterior probability density integrates the prior density of parameters-the information available before measurements-with the likelihood function, which expresses the probability of observations given the parameters in question [14,15].
The posterior probability can be approximated through the Markov Chain Monte Carlo (MCMC) method, a computational tool for approximating integrals and generating samples from a posterior probability.The Metropolis-Hastings algorithm, widely used in MCMC, is based on the acceptance-rejection method and is effective in handling complex probability distributions, enabling random sampling from any probability distribution [16,17].
Therefore, the main objective of the study was to compare different breakthrough curve models for the adsorption of H 2 S and CO 2 in a fixed-bed column system using data obtained from the literature, estimating the parameters of the breakthrough curves through the Markov Chain Monte Carlo method (MCMC) with the application of the Metropolis-Hastings algorithm.

Analytical Models of Breakthrough Curves
When conducting dynamic adsorption experiments in fixed-bed columns, the breakthrough curve is generated as a response.This curve demonstrates the interaction between kinetic and thermodynamic effects of the solute with the adsorbent along the column as the gas or liquid flows.Throughout the process, factors such as adsorption capacity, selectivity, adsorption rate, initial concentration, and fluid velocity can significantly influence the position and shape of the breakthrough curve [11].
Figure 1 shows a schematic of the fixed-bed adsorption process, where the adsorbate (gas), with initial concentration C 0 , transfers from the solution to the adsorbent surface.Initially, the concentration at the outlet is zero, but it gradually increases as the sites become saturated.The breakthrough curve, which describes this process, is obtained by relating the concentration of the component at the inlet and outlet of the column (C/C 0 ) over time.
The region of highest component removal, known as the mass transfer zone (MTZ), can vary with temperature, component concentration, flow rate, and mass transfer rates [18,19].
become saturated.The breakthrough curve, which describes this process, is obtained by relating the concentration of the component at the inlet and outlet of the column (C/C0) over time.The region of highest component removal, known as the mass transfer zone (MTZ), can vary with temperature, component concentration, flow rate, and mass transfer rates [18,19].The mass transfer region is typically limited by values between 0.95 and 0.05 of the relative concentration (C/C0).Breakthrough (tr) occurs when the relative concentration reaches 0.05.If adsorption continues beyond this point, the outlet concentration increases until it reaches the inlet concentration (C/C0 = 1.0), indicating bed saturation (ts) [18].
Modeling these curves is essential for predicting variables and determining optimal conditions in fixed-bed columns, whether through a theoretical model using mass transfer coefficients and diffusivity or an analytical approach involving experiments with columns of varying lengths and parameters [20,21].Thus, for fitting to experimental data, the analytical breakthrough curve models of Thomas, Yoon-Nelson, Adams-Bohart, and Yan were used.

Thomas Model
The Thomas model, initially proposed in 1944 and later refined in 1948, is an approach used to describe adsorption in fixed-beds.In the initial formulation, the capture rate is considered an ion exchange reaction between ionic species "A" in the fluid phase and ionic species "B" in the adsorbent phase, governed by reaction rate constants for ion exchange and regeneration [22,23].
The model is based on Langmuir kinetics for adsorption-desorption and does not consider axial dispersion within the column.The adsorption rate is governed by reversible second-order reaction kinetics, and the model assumes a constant separation factor, applicable to both favorable and unfavorable isotherms [23].The mass transfer region is typically limited by values between 0.95 and 0.05 of the relative concentration (C/C 0 ).Breakthrough (t r ) occurs when the relative concentration reaches 0.05.If adsorption continues beyond this point, the outlet concentration increases until it reaches the inlet concentration (C/C 0 = 1.0), indicating bed saturation (t s ) [18].
Modeling these curves is essential for predicting variables and determining optimal conditions in fixed-bed columns, whether through a theoretical model using mass transfer coefficients and diffusivity or an analytical approach involving experiments with columns of varying lengths and parameters [20,21].Thus, for fitting to experimental data, the analytical breakthrough curve models of Thomas, Yoon-Nelson, Adams-Bohart, and Yan were used.

Thomas Model
The Thomas model, initially proposed in 1944 and later refined in 1948, is an approach used to describe adsorption in fixed-beds.In the initial formulation, the capture rate is considered an ion exchange reaction between ionic species "A" in the fluid phase and ionic species "B" in the adsorbent phase, governed by reaction rate constants for ion exchange and regeneration [22,23].
The model is based on Langmuir kinetics for adsorption-desorption and does not consider axial dispersion within the column.The adsorption rate is governed by reversible second-order reaction kinetics, and the model assumes a constant separation factor, applicable to both favorable and unfavorable isotherms [23].

Yoon-Nelson Model
The study by Yoon and Nelson published in 1984 was developed with the aim of applying gas adsorption kinetics principles to develop a new theoretical model related to the adsorption and breakthrough of vapors or contaminant gasses on solid adsorbents, specifically activated carbon [24].
The model proposed by Yoon and Nelson is based on the idea that the rate at which the probability of adsorption for each adsorbate molecule decreases over time is proportional to both the initial probability of adsorption and the probability of the adsorbate advancing into the interior of the adsorbent.A distinctive feature of this model is its simplification compared to more complex formulations involving mass balance along the adsorption bed, eliminating the need for detailed information on contaminant characteristics, adsorbent type, or physical bed properties [24,25].

Adams-Bohart Model
In 1920, Adams and Bohart conducted the first investigation of fixed-bed dynamics by analyzing the adsorption of chlorine gas in columns filled with activated carbon.Within this context, the proposed model considers that the adsorption rate is directly linked to the remaining capacity of the adsorbent and the concentrations of the adsorbed species.Additionally, it is recognized that adsorption equilibrium is not instantaneous and there is no axial dispersion in the adsorption bed [26].
The model also posits that the adsorbate is adsorbed irreversibly, with a local removal rate proportional to both the residual capacity of the adsorbent and the adsorbate concentration in the gas phase.This implies that the adsorption rate varies with the availability of adsorption sites on the adsorbent and the concentration of the adsorbate in the gas phase [26,27].

Yan Model
The model developed by Yan, known as the "dose-response" model, aimed to enhance the accuracy of mathematical fitting compared to the Thomas model, with a special focus on scenarios characterized by very long or very short operating times.The model considers the absence of axial dispersion, highly reduced internal and external diffusion resistances, second-order reversible reaction kinetics (pseudo-second order), and Langmuir isotherm [28].

Methodology
The estimation of parameters for the analytical models (Thomas, Yoon-Nelson, Adams-Bohart, and Yan) of fixed-bed adsorption for hydrogen sulfide (H 2 S) and carbon dioxide (CO 2 ) gasses was conducted using experimental data from [29], who analyzed the efficiency of babassu biochar as an adsorbent for H 2 S removal from air under low-pressure conditions in a fixed-bed column, and [30], who conducted an experimental and simulation study using a fixed-bed column for CO 2 adsorption using activated carbon.The authors employed the Linear Driving Force (LDF) model for parameter estimation, which is based on a system of partial differential equations that requires greater complexity for equation solving and subsequent estimation.The development of the model used by the authors is detailed in the Supplementary Material [31].
In this perspective, analytical breakthrough curve models of Thomas, Yoon-Nelson, Adams-Bohart, and Yan (dose-response) were used, aimed at assessing the models' ability to achieve good agreement with experimental data, reducing computational costs, enabling faster analysis, and avoiding the complexity of partial differential equations.This approach is advantageous when working with large datasets or when multiple simulations are required [32].
Each model involves unknown parameters that cannot be directly determined from the observed or available data.Thus, these parameters were inferred through Bayesian inference using the Markov Chain Monte Carlo method (MCMC) with the Metropolis-Hastings algorithm.
In Table 1 are the operating conditions used in the adsorption process; Table 2 presents a summary of the models used with their respective equations, parameters, and assumptions; and Table 3 presents the models with the parameters to be estimated.The model assumes a constant separation factor and is applicable to both favorable and unfavorable isotherms. [23]

Yoon-Nelson
τ-Time required to reach 50% breakthrough (min); t-Flow time (min); k yn -Yoon-Nelson kinetic constant (1/min). - The probability of adsorption of each adsorbate molecule decreasing over time is directly proportional to both the initial probability of adsorption and the probability of adsorbate's progress within the adsorbent. [24] Adams-Bohart  [28] Table 3. Parameters to be estimated.

Models Parameters
Markov Chain Monte Carlo (MCMC) In parameter estimation using the MCMC method with the Metropolis-Hastings algorithm, a uniform prior probability distribution was employed with a lower limit of zero (since all parameters are physically non-negative) and an upper limit chosen to be high (10 times the literature value), as described in Equation ( 1).This choice is justified by the lack of prior knowledge of such parameters; therefore, a uniform prior probability distribution is chosen because it is equiprobable.
Regarding likelihood, it is assumed that the experimental data follow a Gaussian probability distribution and that the experimental information is uncorrelated (Equation ( 2)): where σ 2 T is the variance of the measurement uncertainty; F is the solution of the direct model calculated with the parameters to be estimated; Y meas is the measured value; and P is the value of the parameters of interest.
The sample acceptance/rejection algorithm used was Metropolis-Hastings.In this way, the number of states (n) is defined, and then the interaction counter i = 0 is initiated, and an initial value P (0) is defined.A candidate value P * is generated from the distribution P(P * |P (i) ), according to Equation (3), with w being the search step [15,16,33,34].
Subsequently, the acceptance probability α(P (i) |P * ) of the candidate value is calculated, according to Equation ( 4).An auxiliary random sample is generated from a uniform distribution u~U (0,1), so that if u ≤ α(P * |P (i) ), the new value is accepted and set as P (i+1) = P * ; otherwise, P (i+1) = P (i) .Finally, the counter i is incremented to i + 1 and returns to the third step, where a candidate value P * is generated from the distribution P(P * |P (i) ).
For the selection of models with the best agreement to the experimental data, classical statistical metrics were applied, such as the coefficient of determination (R 2 ) and the adjusted coefficient of determination (R 2 Adjusted ), presented in Equations ( 5) and ( 6), as well as the Bayesian Information Criterion (BIC), represented by Equation (7).Accordingly, a good fit is indicated by R 2 values close to 1, reflecting a high capacity of the model to represent the experimental data, while lower BIC values suggest a better model fit.
In the equations for the coefficients of determination (5-6), SSE represents the sum of squared errors, SST the total sum of squares, Y the observed data from i = 1 to the number of observations in the sample (n), Ŷ the estimated data, Y the mean of the observed data, and NP the number of parameters to be estimated by the model.In the BIC equation, p( Y meas |P) represents the likelihood, NP the number of parameters, and J the number of observations in the sample.

Results and Discussion
In the development of the parameter estimation method, experimental data from [29] for H 2 S adsorption and from [30] for CO 2 adsorption were utilized.The Markov Chain Monte Carlo method with the Metropolis-Hastings algorithm was applied to obtain the parameters.The results obtained for the estimation of maximum adsorption capacity (q s , N 0 ) were compared to those obtained deterministically in the literature.Subsequently, based on selection metrics, it was determined which models showed better agreement with the experimental data.
Thus, the number of states of the Markov chain (n) was set in a range from 2000 to 10,000, and the search step (w) used ranged from 0.003 to 0.03, ensuring the stabilization of the chains.For uncertainties, a variance of 1% of the maximum experimental measurement was applied, which is related to the ratio between concentrations (C/C 0 ) in the breakthrough curves.
The elements used for parameter estimation such as the number of states (n), search step (w), and heating states (aq) for each model, with different initial concentrations, are included in the Supplementary Material.Table 4 presents the parameter estimates of the analytical models for H 2 S gas adsorption with and without estimation of the maximum adsorption capacity.In the fitting analysis of the Thomas, Yoon-Nelson, Adams-Bohart, and Yan models for H 2 S gas adsorption, it was observed that all models, except Yoon-Nelson, showed good agreement in estimating the maximum adsorption capacity (q s , N 0 ) for the lowest initial concentration (C 0 ) of H 2 S (1.327 mg/L).Although the Yoon-Nelson model does not include q s in its equation for this estimation, it also demonstrated a good fit for the respective concentration.
Among the analyzed models, the Yan model showed the best fit using the estimation of q s , with an R 2 of 0.9982, R 2  Adjusted of 0.9981, and BIC of 11.2732.On the other hand, the Adams-Bohart model was the worst fit without the estimation of N 0 , with an R 2 e R 2 Adjusted of 0.9331 and 0.9284, respectively, and a BIC of 3.6145×10 3 .The Thomas model, commonly used to describe adsorption in fixed-bed systems, showed a higher estimated q s for the lowest concentration (1.327 mg/L), suggesting that the model is more effective in describing the H 2 S adsorption process at lower concentrations.The kinetic constant of the model (k th ) decreased with increasing concentration, indicating that the adsorption rate decreases as the initial concentration increases.With the q s , estimation, the model exhibited high R 2 and R 2  Adjusted values, suggesting a good fit to the experimental data.The BIC values were relatively low, and the coefficients of determination remained high even without q s estimation, but they were lower compared to when q s was estimated.
The Yoon-Nelson model also showed a good fit with high R 2 values for both concentrations, although the model does not include the parameter q s .It demonstrated a good representation of the adsorption data.The kinetic constant k yn was higher for the higher concentration (2.577 mg/L), indicating a higher probability of adsorption per unit time, while τ was lower at higher concentrations, suggesting faster saturation of the adsorbent bed.
For the Adams-Bohart model, it was observed that the constant k ba decreased with increasing C 0 (whether estimating N 0 or not).This behavior suggests that the adsorption rate tends to be higher at lower concentrations, which can be explained by the saturation of adsorption sites.At lower initial concentrations, available adsorption sites fill up more quickly, thereby increasing the adsorption rate represented by k ba [21,22,35].Furthermore, when estimating N 0 , the k ba constant adjusted to reflect the change in C 0 , reinforcing the inverse relationship between the kinetic constant and the initial concentration of the adsorbate.
In the Yan model, it was observed that a y decreased with increasing initial adsorbate concentration.This behavior is associated with the possibility that, for H 2 S adsorption, the maximum adsorption capacity decreases as the initial concentration increases.Additionally, a decrease in (qs) with increasing (C 0 ) was observed, which may occur due to competition among adsorbate molecules for adsorption sites [36].
Figure 2a,b presents breakthrough curves simulated using the estimation of maximum adsorption capacity for the Thomas (q s ), Adams-Bohart (N 0 ) and Yan (q s ). Figure 3a,b displays the simulated breakthrough curves with deterministic q s and N 0 obtained from the study by [29], both including the corresponding experimental data curve for the Thomas, Yoon-Nelson, Adams-Bohart, and Yan models.In general, the models showed better fit when estimating the maximum adsorption capacity, reflecting a more accurate description of the adsorption behavior.Among the models analyzed, the Yan model stood out for providing the best fit.This model was developed to reduce the uncertainty associated with the Thomas model, particularly effective in situations involving extremely long or short operating periods.The Yan model considers the absence of axial dispersion, extremely low internal and external diffusion resistances, reversible second-order reaction kinetics (pseudo-second-order), and Langmuir isotherm [28].
Additionally, Table 5 presents parameter estimates of analytical models for carbon dioxide (CO2) adsorption with and without estimation of maximum adsorption capacity (qs, N0).In general, the models showed better fit when estimating the maximum adsorption capacity, reflecting a more accurate description of the adsorption behavior.Among the models analyzed, the Yan model stood out for providing the best fit.This model was developed to reduce the uncertainty associated with the Thomas model, particularly effective in situations involving extremely long or short operating periods.The Yan model considers the absence of axial dispersion, extremely low internal and external diffusion resistances, reversible second-order reaction kinetics (pseudo-second-order), and Langmuir isotherm [28].
Additionally, Table 5 presents parameter estimates of analytical models for carbon dioxide (CO 2 ) adsorption with and without estimation of maximum adsorption capacity (q s , N 0 ).For CO 2 adsorption, the Thomas, Adams-Bohart, and Yan models showed good fits in estimating the maximum adsorption capacity for the highest concentration (35 mg/L), as indicated by the determination coefficients (R 2 e R 2 Adjusted ) and the Bayesian Information Criterion (BIC).The Yoon-Nelson model, although it does not estimate q s , also demonstrated a satisfactory fit for this concentration.
By analyzing the selection metrics (R 2 e R 2 Adjusted e BIC), the Adams-Bohart models (estimating N 0 ) and Yoon-Nelson (without estimating q s ) showed the best fits for the 35 mg/L, with R 2 and R 2 Adjusted of 0.9951 and BIC varying between 153.3194 and 153.3209.In contrast, the Adams-Bohart model, without estimating N 0 , exhibited the worst fit for the lowest concentration (18 mg/L), with an R 2 of 0.2020 and R 2 Adjusted of 0.1522.Generally, the maximum adsorption capacity for CO 2 concentrations increased with the rise of (C 0 ).In the paper [37], the authors state that the concentration gradient between the adsorbate and the solute is higher at elevated concentrations, increasing mass transfer and consequently the adsorption capacity.
The Thomas model provided q s estimates close to the deterministic q s .Moreover, it was observed that, in general, the kinetic constant (k th ) decreased with increasing C 0 , suggesting that the adsorption rate tends to decrease at higher concentrations due to adsorbent saturation, a dynamic also reported by [38].
On the other hand, in the Yoon-Nelson model, the constant k yn increased with the rise of C 0 .Thus, it was possible to observe a decrease in the parameter τ, indicating a shorter contact time required for adsorption to reach half of its capacity, relating to a more dynamic behavior of the adsorbent at higher concentrations.
In the Adams-Bohart model, when estimating N 0, it was noted that the constant k ba exhibited non-linear behavior in relation to C 0 .This behavior may occur due to the adsorption dynamics at different concentrations and the possible interference of interactions between CO 2 molecules at high concentrations.Furthermore, the Yan model provided similar maximum adsorption capacity estimates to those of the Thomas model, also indicating higher q s at higher concentrations.It was observed that the constant a y decreased with increasing concentration, both when estimating q s and using deterministic q s .Figure 4a-c show the breakthrough curves simulated using the estimated maximum adsorption capacity for the Thomas (q s ), Adams-Bohart (N 0 ) and Yan (q s ) models.Figure 5a-c present the breakthrough curves estimated with the deterministic values of q s and N 0 obtained from the study by [30], covering the Thomas, Yoon-Nelson, Adams-Bohart, and Yan models, including the curve corresponding to the experimental data for proper comparison.
As in the analysis for H 2 S gas, the breakthrough curves for CO 2 gas adsorption also showed superior fitting when estimating the maximum adsorption capacity.Among the models examined, the Adams-Bohart model stood out for providing the best fit with the estimate of maximum adsorption capacity at the highest concentration (35 mg/L).However, this model performed the worst in describing the adsorption process without this estimate.
The Adams-Bohart model assumes that the adsorption rate is directly related to the residual capacity of the adsorbent and the concentrations of the adsorbed species, recognizing that adsorption equilibrium does not occur instantly and that there is no axial dispersion in the adsorption bed.Additionally, the model suggests that the adsorbate is irreversibly adsorbed, with a local removal rate proportional to the residual capacity of the adsorbent and the concentration of the adsorbate in the gas phase.Thus, the adsorption rate varies according to the availability of adsorption sites on the adsorbent and the concentration of the adsorbate in the gas phase [26,27].
Although the Adams-Bohart model stood out, the Thomas model also showed good fit for the concentration of 35 mg/g.This model assumes that adsorption kinetics follow a first-order reaction with respect to the adsorbate concentration and that the adsorption isotherm follows the Langmuir model, implying a homogeneous adsorption surface with energetically equivalent adsorption sites.Furthermore, the Thomas model neglects internal and external mass transfer resistances, assuming that the adsorption rate is limited by the adsorption reaction at the active site [21,23].
When comparing the adsorption process of the two gasses, H 2 S and CO 2 , models that considered the estimation of maximum adsorption capacity provided a good description of adsorption processes in fixed-bed columns, better reflecting column dynamics.Conversely, treating this parameter as deterministic tends to result in larger differences between estimated and experimental values.This discrepancy arises because estimating it as a random variable allows for a more flexible adjustment, providing an optimized mathematical fit that better represents the reality of the adsorption system [19].also indicating higher qs at higher concentrations.It was observed that the constant ay decreased with increasing concentration, both when estimating qs and using deterministic qs.
Figure 4a-c show the breakthrough curves simulated using the estimated maximum adsorption capacity for the Thomas (qs), Adams-Bohart (N0) and Yan (qs) models.Figure 5a-c present the breakthrough curves estimated with the deterministic values of qs and N0 obtained from the study by [30], covering the Thomas, Yoon-Nelson, Adams-Bohart, and Yan models, including the curve corresponding to the experimental data for proper comparison.Moreover, it was observed that the accuracy of the models can be influenced by the initial concentrations of the adsorbates.During the adsorption of H 2 S at low concentrations, the models presented a more precise fit, whereas in the adsorption of CO 2 at low concentrations, the models showed a poorer fit to the experimental data.This may occur because analytical models rely on simplified assumptions, such as bed uniformity, the absence of interactions between adsorbed molecules, and the absence of axial dispersion, which may not fully capture the complexity of real systems.
Table 6 presents estimates of the parameters of analytical models of breakthrough curve adsorption of H 2 S and CO 2 gasses in fixed-bed systems found in other reports in the literature.These estimates, although derived from different operational conditions, are roughly consistent with the order of magnitude obtained from the experimental data of [29,30].As in the analysis for H2S gas, the breakthrough curves for CO2 gas adsorption also showed superior fitting when estimating the maximum adsorption capacity.Among the models examined, the Adams-Bohart model stood out for providing the best fit with the estimate of maximum adsorption capacity at the highest concentration (35 mg/L).However, this model performed the worst in describing the adsorption process without this estimate.
The Adams-Bohart model assumes that the adsorption rate is directly related to the residual capacity of the adsorbent and the concentrations of the adsorbed species, recognizing that adsorption equilibrium does not occur instantly and that there is no axial dispersion in the adsorption bed.Additionally, the model suggests that the adsorbate is irreversibly adsorbed, with a local removal rate proportional to the residual capacity of the adsorbent and the concentration of the adsorbate in the gas phase.Thus, the adsorption rate varies according to the availability of adsorption sites on the adsorbent and the concentration of the adsorbate in the gas phase [26,27].

Conclusions
This study showed that the use of Bayesian inference for model fitting in continuous adsorption systems is promising, enabling better understanding and optimization of gas separation and purification processes.
Overall, models that included the estimation of maximum adsorption capacity provided more accurate descriptions of adsorption system behavior.On the other hand, it was noted that when maximum adsorption capacity is treated as a deterministic parameter, the difference between estimated and experimental values tends to be larger.
The analysis of determination coefficients (R 2 and R 2 Adjusted ) and Bayesian Information Criterion (BIC) reinforced the efficiency of the fitted models.For H 2 S, the Yan model best described the gas adsorption, obtaining R 2 values above 0.998 and a lower BIC value of 11.2732, whereas for CO 2 , the standout models were Adams-Bohart and Yoon-Nelson, with R 2 values close to 0.995 and a BIC of 153.3.These results confirm the effectiveness of the analytical models in representing fixed-bed adsorption processes, offering a viable alternative with lower computational effort compared to models based on complex differential equations.
Thus, the study highlights the importance of modeling in fixed-bed column adsorption, contributing to the development of more efficient systems for capturing atmospheric pollutants, such as H 2 S e CO 2 , with the potential for large-scale implementation.Additionally, the application of Bayesian inference provides a promising alternative to overcome the limitations of traditional parameter estimation methods, offering improved fits/descriptions and enhanced reliability in analyses.

Figure 1 .
Figure 1.Schematic representation of the adsorption process in a fixed bed.

Figure 1 .
Figure 1.Schematic representation of the adsorption process in a fixed bed.

Figure 2 .Figure 2 .
Figure 2. (a) Breakthrough curves for concentration of 1.327 mg/L and (b) concentration of 2.577 mg/L with estimated qs and N0.

Figure 2 .Figure 3 .
Figure 2. (a) Breakthrough curves for concentration of 1.327 mg/L and (b) concentration of 2.577 mg/L with estimated qs and N0.

Figure 3 .
Figure 3. (a) Breakthrough curves for concentration of 1.327 mg/L and (b) concentration of 2.577 mg/L with deterministic q s and N 0 .

Figure 4 .Figure 4 .
Figure 4. (a) Concentration of 18 mg/L, (b) concentration of 29 mg/L, and (c) concentration of 35 mg/L with the estimation of qs and N0.

Figure 5 .
Figure 5. (a) Concentration of 18 mg/L, (b) concentration of 29 mg/L, and (c) concentration of 35 mg/L with deterministic q s and N 0 .

Table 2 .
Summary of the models employed in the study.

Table 4 .
Parameter estimation of analytical models for H 2 S gas.

Table 5 .
Parameter estimation of analytical models for CO 2 gas.

Table 6 .
Parameter estimates of analytical models obtained from the literature.