{"id":"bb241015-4e60-490f-a835-75d5248de186","arxiv_id":"2607.17941","paper_version":1,"verdict":"CONDITIONAL","confidence":"HIGH","novelty_score":7.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":5,"one_line_summary":"The PFO model coupled with the Sips isotherm predicts finite-time abrupt breakthrough and yields a concentration-dependent effective rate constant, indicating it is unfit for column adsorption dynamics.","lead":"An analysis of the widely used pseudo-first-order (PFO) adsorption model shows it predicts an abrupt, nonphysical breakthrough and that its fitted rate constant changes with inlet concentration, unlike the consistent Sips model. The authors conclude that the century-old PFO model is structurally unsuited to describing column adsorption.","discovery_kind":"extension","skeptic_critique":{"model":"deepseek-v4-flash","headline":"The empirical evidence for PFO's structural flaw rests on an unverified small-Da/Pe^{-1} reduction; if Pe^{-1} is not small for the datasets, the fitted k_LDF variation may be an artifact.","rationale":"The paper's central claim—that PFO is structurally unfit—rests on two pillars: the mathematical prediction of a discontinuous breakthrough for α<1, and the empirical demonstration that the fitted k_LDF varies with inlet concentration while the Sips k_a does not. The reader correctly identifies the model reduction in Section 3.4 as the weakest assumption: both derivations assume Da, Pe^{-1} ≪ 1. The paper asserts these are 'typically small' but never verifies them for the three datasets. My order-of-magnitude estimates indicate Pe^{-1} may indeed be order 1 or larger for these small-scale columns, because the adsorption length ℓ can be very small (~10^{-4} m) while axial dispersion is appreciable. If Pe^{-1} is not small, the analytical breakthrough formulas (59) and (61) do not describe the experimental data, and the fitted parameters lose physical meaning. The observed trend in k_LDF could then be an artifact of systematically increasing approximation error with concentration, not evidence of an inherent model flaw. Similarly, the finite cut-off time is a property of the reduced ODE; with dispersion the front is smoothed, so the claimed 'structural distinction' is at least partially an artefact of the approximation. This concern is concrete and testable: estimating Pe^{-1} from the published conditions and, if necessary, re-fitting with a full numerical solution would settle it. The paper's analytical work is careful and the mathematical results for the reduced system appear sound, but their applicability to the experimental validation is not established. The reader's conditional verdict is appropriate; I do not find a more severe internal inconsistency. The same concern was the reader's weakest assumption, so I agree with that identification.","tokens_in":21116,"tokens_out":8247,"duration_ms":70491,"concrete_test":"Estimate Pe^{-1} for each dataset: compute ℓ = u ε c_in τ/(ρ_b q_max) with τ = 1/(k_a c_in^m q_max^{n-1}) from Table 2 values, and D from the Chung–Wen or Gunn correlation for axial dispersion (or D ≈ 0.73 D_m + 0.5 u d_p). If any Pe^{-1} is not ≪1, fit the full PDE (1) with boundary conditions (10) and PFO kinetics (9) numerically (e.g., method of lines) using the same isotherm parameters, and compare to the analytical curve (59). If the numerical breakthrough deviates from (59) (e.g., by >5% in c_b or t_b), re-fit k_LDF and k_a to the numerical solution; if the k_LDF vs c_in trend weakens or vanishes, the paper's central conclusion is unsupported.","verdict_should_be":"UNCHANGED","load_bearing_attack":"Section 3.4 obtains the first-order PDE (21) and the traveling-wave breakthrough formulas (59)–(61) by neglecting Da and Pe^{-1}. For the experimental datasets (Table 1), the paper does not estimate these numbers. Da = ε c_in/(ρ_b q_max) is indeed tiny, but Pe^{-1} = D/(ℓ u) with ℓ = u ε c_in τ/(ρ_b q_max) can be order 1 or larger: using the fitted k_a values (Table 2), ℓ for dataset 1 is ~10^{-4} m, so with typical axial dispersion coefficients Pe^{-1} ~ 1–10. If Pe^{-1} is not ≪1, Eq. (21) is not a valid approximation; the analytical breakthrough curves (59) and (61) are then not applicable, and the fitted k_LDF and k_a are effective parameters that absorb neglected dispersion/accumulation. The reported monotonic increase of k_LDF with c_in (Figure 12) could therefore reflect increasing approximation error rather than an intrinsic property of the PFO model. The predicted discontinuous breakthrough at t_b for α<1 is also a property of the reduced ODE; finite dispersion regularizes the singularity, weakening the claimed structural distinction.","agreement_with_reader":"agree"},"referee_report":{"model":"deepseek-v4-flash","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.","tokens_in":21450,"tokens_out":8470,"duration_ms":77145,"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":[{"comment":"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":"Section 4.2-4.3, Eq. (39)-(41), Figures 4-5"},{"comment":"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":"Section 5.5, Appendix D, Table 2"},{"comment":"The (m,n) selection is partly data-driven; a robustness check is recommended.","section":"Section 5.5, Appendix D, Table 2"}],"minor_comments":[{"comment":"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}.","section":"Eq. (46)"},{"comment":"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":"Table C.3 vs Table 2"},{"comment":"Typo: 'effet' should be 'effect' in the sentence beginning 'The effet on the travelling wave profiles...'.","section":"Section 5.1"},{"comment":"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.","section":"Section 3.3"},{"comment":"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').","section":"References"}],"recommendation":"major_revision","confidential_remarks":"The core mathematical derivations are solid, and if the authors can convincingly establish the small-Pe^{-1} regime for the experimental datasets, the paper would be a valuable contribution. The stress-test concern about Pe^{-1} is real: for dataset 1 the estimates suggest Pe^{-1} of order 1 or larger, which would invalidate the empirical comparison. I recommend major revision, not rejection, because the issue is fixable by either recomputing with the full PDE or by restricting the claims to parameter regimes where the reduction is justified. The abstract's sweeping statement that the PFO model is 'inherently flawed and structurally unfit' should be softened or explicitly conditioned on the validity of the asymptotic reduction."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"Marc — this one is worth a look if you care about adsorption kinetics. The authors derive travelling wave solutions for the pseudo-first-order model coupled with the Sips isotherm, and they show that for Sips exponent α<1 the PFO model gives a discontinuous breakthrough with a finite cut-off time, unlike the smooth Sips waves. That existence condition (α≤1) and the analytical implicit solutions for α=1,1/2,1/3,2/3 are new as far as I know, and the derivations are careful. The paper also fits both models to three literature datasets and finds that the PFO rate constant k_LDF varies with inlet concentration while the Sips rate constant k_a stays roughly constant. Nice demonstration.\n\nThe main soft spot is that the travelling wave reduction relies on Da, Pe^{-1} << 1. The paper asserts these are 'typically small' but never estimates them for the datasets used in Section 5.5. The Damköhler number is indeed tiny (ε c_in/(ρ_b q_max) ~10^{-4}), but the inverse Péclet number can be order 1 for these columns if axial dispersion is typical. If Pe^{-1} isn't negligible, equation (21) isn't valid, the analytical breakthrough curves (59) are not applicable, and the fitted k_LDF values are effective parameters absorbing neglected dispersion. The abrupt breakthrough would also be regularized by dispersion, so the claimed structural distinction may be an artifact of the reduction. This doesn't necessarily invalidate the mathematical results, but it weakens the empirical claim. The authors need to estimate Pe for their datasets and show the reduction is justified, or present a dispersion-aware comparison.\n\nA second, minor issue: the (m,n) orders are selected partly by comparing the breakthrough fits themselves, and the isotherm data alone don't pin down α. They justify the choices chemically, but the model-selection step has an inherent circularity. Also, no error bars on the fitted parameters.\n\nThat said, the paper is honest about its assumptions and the derivations check out. If the authors address the dispersion question, this could be a strong piece. It deserves a serious referee. I'd send it out.","headline":"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.","tokens_in":21910,"tokens_out":2717,"would_cite":true,"duration_ms":24714,"reading_group":"maybe","serious_thinker":"yes","would_accept_peer_review":true},"rs_alignment":null,"lean_confirmation":null,"pith_extraction":{"msc":["35C07","76S05","80A30"],"pacs":[],"model":"deepseek-v4-flash","headline":"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.","keywords":["adsorption","pseudo-first order model","Sips isotherm","travelling waves","breakthrough curve","packed column","kinetic modelling","linear driving force"],"falsifier":"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.","tokens_in":1451,"feed_emoji":"🧪","tokens_out":6526,"duration_ms":89664,"temperature":0.7,"pith_summary":"The paper tries to settle whether the pseudo-first-order (PFO) kinetic model, a workhorse of adsorption-column modelling for nearly a century, can be trusted as a physical description of breakthrough curves. It shows that when PFO is coupled to the Sips isotherm, the mathematics becomes structurally different from the underlying Sips kinetics: for Sips exponent α<1 the predicted outlet concentration stays exactly zero until a finite breakthrough time and then jumps, whereas the consistent Sips model gives a smooth rise. Fitting both models to experimental data for three different adsorbent–contaminant systems, the paper finds the PFO rate constant kLDF climbs with inlet concentration while the Sips adsorption rate ka stays essentially constant. The conclusion is that apparent good fits of PFO are curve-fitting artefacts, not evidence of a valid mechanism, so parameters extracted with PFO are not transferable across operating conditions.","feed_headline":"PFO adsorption model is structurally unfit for columns","feed_subtitle":"Its fitted rate constant shifts with inlet concentration; the consistent Sips rate holds steady across three datasets.","key_machinery":"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","core_discovery":"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","pith_inferences":["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."],"forward_implications":["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."],"fun_headline_variants":["PFO breakthrough is abrupt, unlike smooth Sips predictions","PFO rate constant shifts 13-fold; Sips stays steady","Abrupt breakthrough exposes PFO's fundamental flaw","PFO's structural flaw: rate constant swings with feed","PFO breakthrough is a step function, not Sips smooth rise"],"cache_read_input_tokens":23296,"weakest_assumption_plain":"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.","fun_headline_variants_meta":{"raw":{"variants":["PFO breakthrough is abrupt, unlike smooth Sips predictions","PFO rate constant shifts 13-fold; Sips stays steady","Abrupt breakthrough exposes PFO's fundamental flaw","PFO's structural flaw: rate constant swings with feed","PFO breakthrough is a step function, not Sips smooth rise"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000711,"raw_usage":{"total_tokens":3005,"prompt_tokens":680,"completion_tokens":2325,"prompt_tokens_details":{"cached_tokens":256},"prompt_cache_hit_tokens":256,"prompt_cache_miss_tokens":424,"completion_tokens_details":{"reasoning_tokens":2242}},"tokens_in":424,"tokens_out":2325,"duration_ms":14376,"temperature":1.0,"reasoning_tokens":2242,"cache_read_input_tokens":256,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-01T16:32:39.106487+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"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.","supporting_citations":[],"review_version":1}