REVIEW 3 major objections 5 minor 40 references
On the use of equilibrium models to describe dynamic adsorption regimes
T0 review · 3 major / 5 minor · reviewed 2026-08-01 · deepseek-v4-flash
Pith's one-line read PFO kinetics coupled to the Sips isotherm is structurally inconsistent: it yields abrupt breakthrough for α<1 and a rate constant that depends on inlet concentration.
desk verdict Solid analytical work on the PFO-Sips travelling wave model; the empirical indictment of PFO is suggestive but depends on an unverified small-dispersion reduction. read the letter →
The pith
A machine-rendered reading of the paper's core claim, the machinery that carries it, and where it could break.
The reading
What carries the argument
The travelling-wave reduction of the advection–dispersion–adsorption system. In scaled coordinates the Sips-PFO kinetic equation becomes dC/dη = γ g_α(C) with g_α(C) = [v C^α − (C^α + v − 1)C] / [v(C^α + v − 1)], where α = m/n is the Sips exponent and v = 1 + κ^{1/n} is the wave speed. This function decides the qualitative behaviour: for α ≤ 1 it has the correct sign on the interval (0,1), so a heteroclinic travelling wave exists; for α < 1 the integral ∫_0^{1/2} ds/g_α(s) converges, producing the finite cut-off η*_α and the abrupt breakthrough; for α > 1 the steady state at C = 0 becomes unstable and an extra root can block the wave entirely. The same function is then fitted to experimental
What would settle it
Run a column experiment in the α<1 regime, such as the Cu(II) system, with high-time-resolution sampling near breakthrough: the Sips-PFO model predicts outlet concentration identically zero until t_b and then a jump, while the Sips model predicts a smooth continuous rise; a measured smooth non-zero tail before the predicted t_b would falsify the structural claim. Alternatively, fix all operating conditions and vary inlet concentration by a factor of ten; if the fitted kLDF does not change, the central inconsistency is absent.
Extended reading notes
Core claim
The paper's central discovery, on its own terms, is that the pseudo-first-order model combined with the Sips isotherm is not a mild approximation of consistent Sips kinetics but a different dynamical class. The travelling-wave analysis shows that for α=m/n<1 the Sips-PFO equation has a finite cut-off time η*_α: before that time the outlet concentration is exactly zero, so breakthrough is abrupt rather than smooth. Fitting to three datasets then shows an independent inconsistency: the fitted Sips rate ka varies by at most factors of 1.3, 1.1, 1.6 across inlet concentrations, while the fitted PFO rate kLDF increases by factors of 2.5, 2.0, 13.5. The paper concludes that PFO is structurally unf
Load-bearing premise
The analytical derivation assumes the Damköhler and inverse Péclet numbers are so small that accumulation and dispersion terms can be dropped from the mass-balance equation, and the paper does not verify this smallness for the three experimental datasets before applying the analytical breakthrough curves.
Editorial extensions
If this is right
- If the Sips-PFO model is used in the α<1 regime, the predicted breakthrough curve is discontinuous: it stays at zero until t_b and then jumps, so any dataset with a smooth onset cannot be represented by the model's functional form.
- The PFO rate constant kLDF extracted at one inlet concentration will not carry over to another concentration; predictions for a new experiment require refitting, so the model is not predictive.
- The consistent Sips model gives a rate constant ka that remains stable across inlet concentrations, meaning a single ka can describe multiple breakthrough curves, which the PFO model cannot do.
- For α>1, Sips-PFO travelling waves may fail to exist unless the wave speed is high enough, so the model can fail even to provide a bounded wave solution in some parameter regimes.
- In the Henry limit (K_S c_in << 1) both models reduce to the same linear form, hiding the structural flaw precisely in the regime where the fitted curves look similar; agreement in that limit should not be read as validation of PFO.
Reading between the lines
- Because the Sips isotherm is mathematically equivalent to the Hill equation, the same structural objection would transfer to any Hill-type kinetic system approximated by a linear driving force in q with a borrowed equilibrium curve; this extends beyond the adsorption examples tested here.
- The observed concentration dependence of kLDF could be used as a cheap consistency diagnostic for any proposed kinetic model: fit at two inlet concentrations and check whether the rate constant moves; no single-curve R2 value can substitute for that test.
- Published studies that report PFO fits for column adsorption could be re-examined to see whether their fitted kLDF values drift with inlet concentration or bed length; if they do, the curves are interpolation rather than physical prediction.
Editorial analysis
A structured set of objections, weighed in public.
Referee Report
Summary. The paper analyzes a Pseudo-First-Order (PFO) kinetic model coupled with the Sips isotherm (S-PFO) for fixed-bed column adsorption, as an alternative to the thermodynamically consistent Sips kinetic model. Using a traveling-wave ansatz and assuming Da, Pe^{-1} \ll 1, the authors derive implicit analytical breakthrough curves for the S-PFO model for \alpha = 1, 1/2, 1/3, 2/3, and show that for \alpha < 1 the reduced first-order PDE yields a finite breakthrough time \hat{t}_b with an abrupt jump in outlet concentration. They compare these solutions with the Sips traveling-wave results and fit both models to three experimental datasets (toluene, Cu(II), Hg(II)). They report that the S-PFO model gives consistently worse fits and that the fitted k_{LDF} increases monotonically with inlet concentration, while the Sips rate constant k_a remains approximately constant. From this they conclude that the PFO model is 'inherently flawed and structurally unfit' for column adsorption modeling.
Significance. If correct, the paper would be a significant contribution to adsorption modeling: it provides closed-form traveling-wave solutions for a widely used but thermodynamically inconsistent kinetic formulation, and it challenges nearly a century of PFO-based column modeling with a clear, falsifiable observation (the concentration dependence of fitted k_{LDF}). The mathematical derivations in Section 4 and the appendices are careful and reduce correctly to known limits, and the experimental datasets cover three different systems. However, the central empirical and qualitative claims rest on the small-Da and small-Pe^{-1} reduction that is not verified for the experimental conditions. Because the paper's broad conclusion depends on that reduction, the current evidence is not yet sufficient; the analysis is sound in principle but needs either validation of the reduction or a more cautious interpretation of the results.
major comments (3)
- [Section 4.2-4.3, Eq. (39)-(41), Figures 4-5] The experimental validation and the central comparison between Sips and S-PFO rely on the reduction from Eq. (12) to the first-order PDE (21), which assumes Da, Pe^{-1} \ll 1. The paper states in Section 3.3 that these parameters 'are typically considered small' but never estimates them for the datasets in Table 1. For dataset 1, using the reported operating conditions and fitted k_a from Table 2, one estimates \ell \approx u\epsilon/(\rho_b q_{max} k_a) \approx 1.1\times 10^{-4} m; with a typical axial dispersion coefficient in the range D \sim 10^{-5}--10^{-4} m^2/s, Pe^{-1}=D/(\ell u) is of order 1--10, not small. If Pe^{-1} is not negligible, Eq. (21), and hence the analytical breakthrough curves (59) and (61), are not applicable to those experiments. The fitted k_{LDF} and k_a values in Table 2 would then be effective parameters that absorb neglected dispersion and accumulation, and
- [Section 5.5, Appendix D, Table 2] The predicted finite cut-off time arises from the reduced first-order PDE and would be regularized by dispersion; the structural distinction claim needs qualification.
- [Section 5.5, Appendix D, Table 2] The (m,n) selection is partly data-driven; a robustness check is recommended.
minor comments (5)
- [Eq. (46)] In the denominator of the expression for g_{1/3}, the term s^{1/2} appears; it should be s^{1/3} to be consistent with the definition of g_\alpha and with the subsequent change of variable u=s^{1/3}.
- [Table C.3 vs Table 2] The R^2 values for the S-PFO linear fit in Table C.3 differ substantially from the R^2 values reported in Table 2 for the same fits (e.g., dataset 3 at the highest c_{in}: 0.7176 vs 0.9853). While the text explains that the two tables use different definitions, the large discrepancy may confuse readers. Please state explicitly in the captions which quantity is plotted or tabulated and why the linear-fit R^2 differs from the breakthrough-curve R^2.
- [Section 5.1] Typo: 'effet' should be 'effect' in the sentence beginning 'The effet on the travelling wave profiles...'.
- [Section 3.3] The statement that Da and Pe^{-1} are 'typically considered small' would benefit from a reference or a concrete order-of-magnitude estimate for the systems of interest. As discussed in the major comments, this is not just a presentation issue for dataset 1.
- [References] Reference [12] is listed as 'under review' and reference [38] as 'submitted'; if possible, update these with publication or preprint identifiers. Also, reference [29] contains a typo in the DOI ('0.1016' instead of '10.1016').
Circularity Check
No significant circularity: the PFO travelling-wave results and the model comparison are derived from the stated equations, not from fitted outputs or self-citations.
full rationale
The central qualitative claims are derived in a self-contained way. The travelling-wave existence condition (alpha <= 1) and the finite cut-off time for alpha < 1 follow from the sign analysis of g_alpha in Eq. (27) and from the convergence of Phi_alpha(0+) in Appendix B; these are mathematical consequences of the stated PFO kinetic equation coupled to the Sips isotherm, not of fitted parameters. The Sips benchmark is not imported as a black box: the paper derives the Sips isotherm from F=0 in Eq. (6), and the comparison in Figs. 5-6 uses the explicit ODEs (26) and (27). The experimental section fits only the rate constants k_a and k_LDF after fixing isotherm parameters from equilibrium data; the reported variation of k_LDF with c_in is a fitted-parameter comparison, not a prediction forced by construction, and the same (m,n) configuration is used for both models, so the model-order selection does not bias the relative comparison. The main weakness noted in Section 3.4 is that Da, Pe^{-1} << 1 is assumed but not verified for the experimental datasets; however, this is an unverified approximation condition and a correctness risk, not a circularity, because Eqs. (21)/(59) are not defined in terms of the experimental outcomes they are used to judge. The self-citations ([11], [21], [25]) are to prior published derivations and datasets used as benchmarks or supporting context; none substitutes for the derivation performed in this paper. No exhibited equation reduces to another by construction, and no fitted quantity is renamed as a prediction.
Assumptions & free parameters
free parameters (5)
- qmax =
4.059 (dataset 1), 0.1208 (dataset 2), 0.3789 (dataset 3) mol/kg
- KS =
406.09, 4.0016, 0.2688 (m^3/mol)^α
- ka =
Varies per dataset and cin; e.g., 1.12e-1 to 1.07e-1 m^3/(mol s) for dataset 1
- kLDF =
Varies per dataset and cin; e.g., 6.05e-4 to 1.53e-3 s^-1 for dataset 1
- (m,n) =
(1,1), (1,2), (1,3) for datasets 1-3
assumptions (5)
- domain assumption Trace adsorbate assumption: fluid velocity u is constant along the column (Section 2.1).
- domain assumption Da, Pe^-1 ≪ 1: dispersion and accumulation terms negligible relative to advection/adsorption (Section 3.4).
- ad hoc to paper Traveling wave ansatz: solutions depend on x and t through η = x - L - v(t - th) (Section 4.1).
- domain assumption Sips kinetic model (3) is the physically consistent model; the isotherm is its steady state (Section 2.3.1).
- domain assumption Isotherm consistency requirement: a kinetic model should reduce to an isotherm at steady state (Section 1).
Cite this review
Pith. "Pith review of On the use of equilibrium models to describe dynamic adsorption regimes." pith.science (2026). https://pith.science/paper/I6U5NUIQ
@misc{pith2026260717941,
author = {Pith},
title = {Pith review of: On the use of equilibrium models to describe dynamic adsorption regimes},
year = {2026},
howpublished = {\url{https://pith.science/paper/I6U5NUIQ}},
note = {Machine review of arXiv:2607.17941}
}
read the original abstract
We present a column adsorption model that couples a Pseudo-First-Order (PFO) kinetic formulation with the Sips isotherm framework. Using a traveling wave approximation, we derive analytical solutions for specific operating conditions. Qualitatively, these solutions deviate significantly from their pure Sips counterparts: instead of a smooth, continuous increase in concentration at the column outlet, the PFO-Sips model predicts an abrupt, sudden breakthrough. We validate these analytical solutions against diverse experimental datasets from the literature. The results reveal that the PFO-based model consistently underperforms compared to the original Sips formulation. Furthermore, this validation exposes fundamental inconsistencies within the PFO framework. We demonstrate that despite its widespread use in the literature for almost a century, the PFO model is inherently flawed and structurally unfit for describing column adsorption dynamics.
Figures
Figures from the paper (10 more)
Reference graph
Works this paper leans on
-
[1]
D. Tefera, Z. Hashisho, J. Philips, J. Anderson, M. Nichols, Modeling competitive adsorp- tion of mixtures of volatile organic compounds in a fixed-bed of beaded activated carbon, Environmental Science & Technology 48 (2014) 5108–5117.doi:10.1021/es404667f. 32
-
[2]
V. Tzanakopoulou, M. Pollitt, D. Castro-Rodriguez, D. Gerogiorgis, Dynamic model val- idation and simulation of acetone–toluene and benzene–toluene systems for industrial volatile organic compound (VOC) abatement, Industrial & Engineering Chemistry Re- search 63 (2024) 7281–7299.doi:10.1021/acs.iecr.4c00251
-
[3]
N. Casas, J. Schell, R. Pini, M. Mazzotti, Fixed bed adsorption of CO2/H2 mixtures on activated carbon: experiments and modeling, Adsorption 18 (2012) 143–161.doi: 10.1007/s10450-012-9389-z
- [4]
-
[5]
S. Liu, Y. Chen, B. Yue, C. Wang, B. Qin, Y. Chai, G. Wu, J. Li, X. Han, I. da Silva, P.Manuel, S.J.Day, S.P.Thompson, N.Guan, S.Yang, L.Li, Regulatingextra-framework cations in faujasite zeolites for capture of trace carbon dioxide, Chemistry – A European Journal 28 (50) (2022) e202201659.doi:10.1002/chem.202201659
- [6]
-
[7]
A. H. Sulaymon, B. A. Abid, J. A. Al-Najar, Removal of lead copper chromium and cobalt ions onto granular activated carbon in batch and fixed-bed adsorbers, Chemical Engineering Journal 155 (3) (2009) 647–653.doi:10.1016/j.cej.2009.08.021
-
[8]
D. Ruthven, Principles of Adsorption and Adsorption Processes, Wiley-Interscience pub- lication, Wiley, 1984.doi:10.1002/aic.690310335
Show all 40 references
-
[9]
Sips, On the structure of a catalyst surface, Journal of Chemical Physics 16 (5) (1948) 490–495.doi:10.1063/1.1746922
R. Sips, On the structure of a catalyst surface, Journal of Chemical Physics 16 (5) (1948) 490–495.doi:10.1063/1.1746922
1948 doi
-
[10]
I.Langmuir, Theadsorptionofgasesonplanesurfacesofglass, micaandplatinum., Journal of the American Chemical Society 40 (9) (1918) 1361–1403.doi:10.1021/ja02242a004
1918 doi
-
[11]
Aguareles, E
M. Aguareles, E. Barrabés, T. Myers, A. Valverde, Mathematical analysis of a sips-based model for column adsorption, Physica D: Nonlinear Phenomena 448 (2023) 133690.doi: 10.1016/j.physd.2023.133690
2023
-
[12]
J. A. Lloveras, M. Aguareles, E. Barrabés, Analysis of travelling wave equations in sorption processes, under review. Preprint available athttps://doi.org/10.48550/arXiv.2507. 16404
-
[13]
A. V. Hill, The possible effects of the aggregation of the molecules of haemoglobin on its dissociation curves, The Journal of Physiology 40 (1910) 4–7
1910
-
[14]
Goutelle, M
S. Goutelle, M. Maurin, F. Rougier, X. Barbaut, L. Bourguignon, M. Ducher, P. Maire, Thehillequation: areviewofitscapabilitiesinpharmacologicalmodelling, Fundamental& Clinical Pharmacology 22 (6) (2008) 633–648.doi:10.1111/j.1472-8206.2008.00633.x. 33
2008
-
[15]
K. H. Chu, J. Debord, M. Harel, J.-C. Bollinger, Mirror, mirror on the wall, which is the fairest of them all? Comparing the Hill, Sips, Koble–Corrigan, and Liu adsorption isotherms, Industrial & Engineering Chemistry Research 61 (19) (2022) 6781–6790.doi: 10.1021/acs.iecr.2c00507
2022 doi
-
[16]
Lagergren, Zur Theorie der Sogenannten Adsorption Gelöster Stoffe, Bihang till K
S. Lagergren, Zur Theorie der Sogenannten Adsorption Gelöster Stoffe, Bihang till K. Svenska Vet-Akad Handlingar 24 (II(4)) (1898) 1–39. URLhttps://weburn.kb.se/metadata/824/EOD_2758824.htm
-
[17]
Fenti, S
A. Fenti, S. Salvestrini, Analytical solution of the langmuir-based linear driving force model for adsorption kinetics, Journal of Environmental Chemical Engineering 6 (2) (2018) 2893– 2900.doi:10.1016/j.jece.2018.04.002
2018 doi
-
[18]
M. S. Shafeeyan, W. M. A. Wan Daud, A. Shamiri, A review of mathematical modeling of fixed-bed columns for carbon dioxide adsorption, Chemical Engineering Research and Design 92 (5) (2014) 961–988.doi:10.1016/j.cherd.2013.08.018
2014 doi
-
[19]
Y. Ho, G. McKay, Pseudo-second order model for sorption processes, Process Biochemistry 34 (5) (1999) 451–465.doi:10.1016/S0032-9592(98)00112-5
1999 doi
-
[20]
H. C. Thomas, Heterogeneous ion exchange in a flowing system, Journal of the American Chemical Society 66 (10) (1944) 1664–1666.doi:10.1021/ja01238a017
1944 doi
-
[21]
T. G. Myers, A. Cabrera-Codony, A. Valverde, On the development of a consistent math- ematical model for adsorption in a packed column (and why standard models fail), In- ternational Journal of Heat and Mass Transfer 202 (2023) 123660.doi:10.1016/j. ijheatmasstransfer.2022.123660
2023
-
[22]
T. G. Myers, Is it time to move on from the Bohart-Adams model for column adsorption?, International Communications in Heat and Mass Transfer 159 (2024) 108062.doi:10. 1016/j.icheatmasstransfer.2024.108062
2024
-
[23]
Chuang, P
C. Chuang, P. Chiang, E. Chang, Modeling vocs adsorption onto activated carbon, Chemo- sphere 53 (1) (2003) 17–27.doi:10.1016/S0045-6535(03)00357-6
2003 doi
-
[24]
Myers, A
T. Myers, A. Valverde, M. Aguareles, M. Calvo-Schwarzwalder, F. Font, Modelling mass transfer from a packed bed by fluid extraction, International Journal of Heat and Mass Transfer 188 (2022) 122562.doi:10.1016/j.ijheatmasstransfer.2022.122562
2022
-
[25]
Valverde, A
A. Valverde, A. Cabrera-Codony, M. Calvo-Schwarzwalder, T. G. Myers, Investigating the impact of adsorbent particle size on column adsorption kinetics through a mathematical model, International Journal of Heat and Mass Transfer 218 (2024) 124724.doi:10.1016/ j.ijheatmasstrans...
2024
-
[26]
W. L. McCabe, J. C. Smith, P. Harriott, Unit operations of chemical engineering, 5th Edition, McGraw Hill, Inc., Singapore, 1993
1993
-
[27]
A. H. Sulaymon, S. A. Yousif, M. M. Al-Faize, Competitive biosorption of lead mercury chromium and arsenic ions onto activated sludge in fixed bed adsorber, Journal of the Taiwan Institute of Chemical Engineers 45 (2) (2014) 325–337.doi:10.1016/j.jtice. 2013.06.034. 34
2014 doi
-
[28]
T. T. Tran, R. Severac, A. Becker, M. Douard, V. Papaefthymiou, G. Schlatter, A. Hébraud, S. Le Calvé, Toluene adsorption studies on carbon-based materials under con- trolled temperature, humidity, flow rate, and toluene concentration, using a high-precision dynamic adsorption...
2026
-
[29]
X. Yang, H. Yi, X. Tang, S. Zhao, Z. Yang, Y. Ma, T. Feng, X. Cui, Behaviors and kinetics of toluene adsorption-desorption on activated carbons with varying pore structure, Journal of Environmental Sciences 67 (2018) 104–114.doi:0.1016/j.jes.2017.06.032
2018
-
[30]
Ben Amar, M
M. Ben Amar, M. Mallek, A. Valverde, H. Monclús, T. G. Myers, V. Salvadó, A. Cabrera- Codony, Competitive heavy metal adsorption on pinecone shells: Mathematical modelling of fixed-bed column and surface interaction insights, Science of The Total Environment 917 (2024) 170398....
2024
-
[31]
R. I. Soria, S. A. Rolfe, M. P. Betancourth, S. F. Thornton, The relationship between properties of plant-based biochars and sorption of Cd(II), Pb(II) and Zn(II) in soil model systems, Heliyon 6 (11) (2020).doi:10.1016/j.heliyon.2020.e05388
2020 doi
-
[32]
X. Xu, A. Schierz, N. Xu, X. Cao, Comparison of the characteristics and mechanisms of hg(ii) sorption by biochars and activated carbon, Journal of Colloid and Interface Science 463 (2016) 55–60.doi:10.1016/j.jcis.2015.10.003
2016 doi
-
[33]
Yousif, A
S. Yousif, A. Sulaymon, M. Alfaize, Experimental and theoretical investigations of lead mercury chromium and arsenic biosorption onto dry activated sludge from wastewater, International Review of Chemical Engineering 5 (1) (2013) 30–40
2013
-
[34]
van Veenhuyzen, S
B. van Veenhuyzen, S. Tichapondwa, C. Hörstmann, E. Chirwa, H. Brink, High capac- ity pb(ii) adsorption characteristics onto raw- and chemically activated waste activated sludge, Journal of Hazardous Materials 416 (2021) 125943.doi:10.1016/j.jhazmat. 2021.125943
2021
-
[35]
Giraldo, I
S. Giraldo, I. Robles, A. Ramirez, E. Flórez, N. Acelas, Mercury removal from wastewater using agroindustrial waste adsorbents, SN Applied Sciences 2 (6) (2020) 1029.doi:10. 1007/s42452-020-2736-x
2020
-
[36]
B.Deng, C.-Y.Ren, Z.-H.Ren, Z.-J.Zhang, Q.Li, Q.-X.Yuan, H.-P.Zhao, Bioremediation of high concentrations of mercury and nickel in activated sludge by black soldier fly lar- vae: Accumulation, transformation, and sustainability assessment, Journal of Hazardous Materials 508 (2...
2026
-
[37]
Z. He, J. Shen, Y. Zhao, Y. Ru, D. Zhang, X. Pan, Microbial antagonistic mechanisms of hg(ii) and se(iv) in efficient wastewater treatment using granular sludge, Water Research 253 (2024) 121311.doi:10.1016/j.watres.2024.121311
2024
-
[38]
Myers, U
T. Myers, U. Thorat, M. Calvo-Schwarzwalder, A. Valverde, A. Cabrera-Codony, F. Font, S. Sahu, B. Wetton, Modelling the simultaneous adsorption of two competing species, submitted to International Communications in Heat and Mass Transfer, May 2026. 35
2026
-
[39]
L. C. Auton, M. Aguareles, A. Valverde, T. G. Myers, M. Calvo-Schwarzwalder, An analytical investigation into solute transport and sorption via intra-particle diffusion in the dual-porosity limit, Applied Mathematical Modelling 130 (2024) 827–851.doi: 10.1016/j.apm.2024.03.023
2024 doi
-
[40]
L. C. Auton, S. S. Ravuru, S. De, T. G. Myers, A. Valverde, Development and experi- mental validation of a mathematical model for fluoride-removal filters comprising chemi- cally treated mineral rich carbon, Journal of Water Process Engineering 79 (2025) 108914. doi:10.1016/j....
2025
Reviewed August 1, 2026 · model on record in the stance chip above.
Discussion (0). Continue with ORCID to comment.