{"id":"4acd3e80-8ca8-4c30-9a26-7679250d4323","arxiv_id":"2505.01213","paper_version":1,"verdict":"REJECT","confidence":"MODERATE","novelty_score":4.0,"correctness_risk":"high","formal_verification":"none","parameter_count":8,"one_line_summary":"Catalysis is re-expressed as a spatiotemporal integral of site density, modulation, and local rate, but the examples calculate rather than predict and the claimed unification is a re-labeling of existing kinetic expressions.","lead":"This preprint proposes a 'unified theory of catalysis' by writing turnover frequency as an integral over space and time of site density, modulation, and local rate. The core formula is a definitional restatement of standard kinetics, and the worked examples use the same values as inputs and outputs.","discovery_kind":"extension","skeptic_critique":{"model":"deepseek-v4-flash","headline":"The enzyme demonstration contradicts the central integral: Eq. 3 is additive, but Eq. 17 replaces the sum with a minimum, so the reported BOD TOF does not come from Eq. 3.","rationale":"The reader's verdict is REJECT, and my independent stress-test identifies a concrete internal inconsistency that reinforces rejection rather than changing it. The load-bearing concern is not merely that the factorization in Eq. 2 assumes independence; it is that the paper's own enzyme example violates the additive structure of the central integral. Eq. 3 defines TOF as an integral of local contributions, implying additivity over volume elements. Yet the BOD analysis first writes Eq. 16 as a sum of two volumetric integrals, then asserts Eq. 17, a minimum over the two site contributions. These two expressions are numerically and structurally incompatible. This matters because the paper's strongest claim is that the ACS framework is a unified, quantitative theory from which all catalytic behavior, including enzymatic series kinetics, emerges. If the central integral cannot produce the min constraint without an extra nonlocal postulate, then Eq. 3 is not the universal predictive formula claimed. The rest of the applications are mostly bookkeeping: constants chosen to reproduce literature rates, with no independent validation. The TST 'limiting case' is obtained by inserting the Eyring expression for k_local and grouping factors into P‡, so it is an assumption rather than a derivation. None of these observations require questioning the author's intent; they concern the argument's support. The proposed concrete test settles the internal inconsistency directly: recompute the BOD TOF from Eq. 16 and compare to Eq. 17. The discrepancy will either confirm the problem or reveal a hidden assumption that must be stated. Since the central claim already fails on this one worked example, the reader's REJECT verdict stands.","tokens_in":11655,"tokens_out":3249,"duration_ms":33583,"concrete_test":"Take the BOD parameters exactly as stated: ρ=1, f_T1=1.0, k_T1=674.8 s⁻¹, f_TNC=0.85, k_TNC=96.18 s⁻¹. Evaluate Eq. 16 as written: 674.8 + 81.76 = 756.6 s⁻¹. Compare with Eq. 17's min = 81.76 s⁻¹. If the two disagree, then the enzyme example does not use Eq. 3, and the universal TOF integral is not the predictive formula for that case. Optionally, attempt to derive the min operation from Eq. 3 by adding a series constraint term; if no local/volume-integral expression reproduces min(S1,S2), the integral cannot encode the stated rate-determining-step coupling.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The central claim is that TOF(t)=∫ρ f k_local dV is a universal descriptor from which TST and all catalytic regimes emerge. The BOD application is internally inconsistent with that claim. After writing Eq. 16 as the sum of T1 and TNC volumetric integrals, the paper states that the sites act in series and sets TOF_enzyme = min(ρ_T1 f_T1 k_T1, ρ_TNC f_TNC k_TNC) (Eq. 17). Plugging their own numbers gives a sum of 674.8 + 81.76 = 756.6 s⁻¹, not the reported 81.76 s⁻¹. Thus the enzyme TOF is not the value of Eq. 3; it is the output of a non-additive series constraint. A theory whose central equation is an additive volume integral cannot claim to encompass a system whose own analysis requires a global minimum over sites, unless that constraint is derived from the integral or added as a separate postulate. It is not derived. This is not a parameter error; it is a structural mismatch between the proposed universal equation and one of its three headline demonstrations. The same pattern—choosing field functions to reproduce known rates—appears in the Pt and Ru examples, but the enzyme case is the cleanest demonstration that Eq. 3 is not doing the predictive work.","agreement_with_reader":"partial"},"referee_report":{"model":"deepseek-v4-flash","summary":"The paper proposes a 'theory of the active catalytic space' (ACS) in which the turnover frequency is written as a volume integral over three fields: catalytic site density ρ(r), a dynamic modulation function f(r,t), and a local per-site rate k_local(r,t) (Eq. 3). The authors claim that this integral unifies heterogeneous, homogeneous, and enzymatic catalysis, that transition state theory emerges as a limiting case, and that the framework 'predicts turnover performance with high fidelity.' Applications are presented for a single platinum atom and for Pt{111}/{100} facets, for the enzyme bilirubin oxidase (BOD), and for a homogeneous ruthenium water-oxidation catalyst. The central mathematical object is a definition, and the paper's demonstrative examples largely consist of choosing functional forms and parameters to reproduce known rate ranges.","tokens_in":12017,"tokens_out":3991,"duration_ms":42192,"significance":"If the framework were actually predictive and able to subsume transition state theory as a limiting case, a unified spatiotemporal description of catalysis would be a valuable conceptual contribution. However, as written, Eq. (3) is a bookkeeping identity that has not been shown to have predictive content. The applications either reduce to normalization (single-atom Pt), to multiplication of chosen constants (Pt facets), to an ad hoc min() operation that contradicts the additive integral (BOD), or to fitted functions whose agreement with experiment is asserted rather than tested (Ru water oxidation). The paper does compile a useful table of experimental techniques for measuring spatial and temporal catalytic fields, and the dimensional analysis in the Methods is internally consistent, but those positive elements do not compensate for the absence of a derivation of the central claims. The manuscript does not provide machine-checked proofs, reproducible code, parameter-free derivations, or falsifiable predictions.","major_comments":[{"comment":"Eq. (16) writes the BOD turnover frequency as the sum of two volumetric integrals over the T1 and TNC sites, which is the direct application of Eq. (3). Eq. (17) then replaces that sum with the minimum of the two contributions, TOF_enzyme = min(ρT1 fT1 kT1, ρTNC fTNC kTNC). With the paper's own numbers, the sum is 674.8 + 81.76 = 756.6 s⁻¹, not the reported 81.76 s⁻¹. The min() operation is not derived from Eq. (3); it is an additional series constraint imposed by hand. Thus the BOD demonstration does not evaluate the central equation, and it cannot be used as evidence that Eq. (3) describes enzymatic catalysis.","section":"Application for Enzyme Catalysis, Eqs. (16)-(18)"},{"comment":"The homogeneous example assigns a Gaussian ρ(r), a sinusoidal f(t) = 1 + 0.6 sin(2πt), and an exponential k_local(r) = 5 exp(-r/4), with no statement of how these forms or their parameters were determined from the system. The integration yields TOF_mean ≈ 723 s⁻¹, which is then said to be 'consistent with reported experimental ranges' for ruthenium water-oxidation catalysts. This is curve fitting, not prediction: the functional forms and coefficients are free parameters, and the agreement with experiment is asserted without any statistical or uncertainty quantification. The example therefore provides no falsifiable test of the ACS framework.","section":"Application for Homogeneous Catalysis, Eqs. (19)-(22)"},{"comment":"The single-atom platinum result TOF(t) = k0 follows entirely from the imposed normalization ∫ρ dV = 1. Substituting ρ = 1/V into Eq. (3) makes the integral equal to k0 by construction. This is tautological and does not test Eq. (3). In addition, modeling a single atom as a uniform density inside a sphere of arbitrary radius R is physically unmotivated; the reader is given no criterion for choosing R, so the 'dimensionally consistent integral framework' claim is not supported by this example.","section":"Application for Heterogenous Catalysis, Eqs. (5)-(9)"},{"comment":"The paper states that 'a full derivation is provided in the SI' and then asserts that after grouping the density, modulation, and infinitesimal volume into P‡, the ACS expression becomes TOF = P‡ (kBT/h) exp(-ΔG‡/RT), 'formally equivalent' to TST. No derivation appears in the manuscript, and the SI is not available to the reader. The grouping of three independent fields into a single 'effective population' is not shown to follow from Eq. (3), and the treatment omits standard TST ingredients such as transmission coefficients and the distinction between Gibbs energy of activation and enthalpy/entropy contributions. As written, the claim that TST emerges as a limiting case is an assertion, not a derived result.","section":"Transition State Theory (TST) as a Limiting Case of ACS Framework"},{"comment":"Eq. (17) sets ρT1 = ρTNC = 1 and omits the volume element dV, so the quantities fT1 kT1 and fTNC kTNC have units of s⁻¹ and are per-site rates, not volume integrals over the catalytic space. This is dimensionally incompatible with Eq. (3) unless an implicit unit volume is introduced. Similarly, Eq. (10) writes three-dimensional volume integrals over V111 and V100 but evaluates them as area integrals with A111 and A100 (sites/m² × s⁻¹ × m²), mixing 3D and 2D domains. These inconsistencies obscure what quantity the theory actually predicts.","section":"Equation (17) and dimensional analysis"}],"minor_comments":[{"comment":"There are two different items labeled 'Figure 3': the platinum TOF map and the BOD site-resolved TOF schematic. The figures must be renumbered and cross-references updated.","section":"General / Figures"},{"comment":"Reference 9 is cited as the basis for the 'active catalytic space framework,' but that reference concerns heterogeneities of individual catalyst particles and does not define the ACS concept. The framework appears to be introduced in this manuscript, so the citation does not support the stated provenance.","section":"References"},{"comment":"The Methods state that 'The final TOF(t) values were benchmarked against known literature values for each case to ensure model consistency,' but no quantitative comparison, error bars, or benchmark data are reported anywhere in the main text. The reader cannot assess the claimed fidelity.","section":"Methods / Validation"},{"comment":"Eq. (14) uses Marcus theory parameters λ and ΔG without defining them, and Eq. (15) uses [O2] without specifying the concentration value used in the BOD calculation. These definitions and values are needed to reproduce the numbers.","section":"Equations (14)-(15)"},{"comment":"The conclusion states that the framework was 'demonstrated across five distinct case studies,' but the main text presents only three applications (Pt, BOD, Ru). The additional cases (carbonic anhydrase, lactate dehydrogenase, perovskite photoelectrode) are said to be in the SI, which is not part of the submitted manuscript; the count should be corrected or the additional cases included.","section":"Conclusion / Scope"},{"comment":"There are numerous typographical errors, e.g., 'integration overtime' (should be 'over time'), 'Fig, 2' (should be 'Fig. 2'), and inconsistent spacing around commas and mathematical symbols. A careful copyedit is needed.","section":"Typos"}],"recommendation":"reject","confidential_remarks":"The manuscript is best characterized as a perspective or position paper rather than a research article with testable claims. The central equation is a definition, and the three main applications either reduce to chosen inputs or contradict the definition. The internal inconsistency between Eq. (16) and Eq. (17) is the clearest indicator that the unifying claim is not supported by the analysis. If the authors wish to pursue this direction, they would need to (i) derive the series-constraint min() from the integral, (ii) provide out-of-sample predictions with measured or independently computed fields, and (iii) supply the promised SI derivation of TST as a limiting case. As it stands, the paper is not suitable for publication in a research journal."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"Short version: Eq. 3 is a definition, not a theory. Writing catalytic rate as the integral of site density times accessibility times local rate is fine, but it's a restatement of how we already think about distributed catalysts. The paper then dresses it up as the 'active catalytic space' and claims it predicts turnover with high fidelity. It doesn't predict; it fits. The single-atom example reduces to TOF = k0 by normalization. The two-facet platinum example computes an extensive rate (9.3×10^13 reactions/s) and calls it a TOF, which is a per-site quantity. The homogeneous ruthenium example chooses a Gaussian, a sine, and an exponential, and reports that the output lands in experimental ranges—that's curve fitting, not validation.\n\nThe enzyme example is the real problem, and the stress-test note is correct. Eq. 3 is an additive volume integral. BOD's T1 and TNC sites act in series, so the paper quietly replaces the sum with a minimum (Eq. 17). The reported 81.76 s⁻¹ comes from that minimum, not from the integral. The paper even acknowledges the switch. That means the central equation is not generating the headline result for one of the three flagship cases. It's an internal contradiction, not a parameter error.\n\nThe paper does earn some credit. It compiles a useful table mapping experimental techniques (cryo-EM, SECM, transient absorption, impedance) onto the three fields. The pitch that catalysis is spatiotemporally distributed and should be studied with operando tools is reasonable and well-supported by the cited literature. As a perspective, that could be a nice piece. But the mathematical formalism is trivial and the examples don't demonstrate anything beyond arithmetic.\n\nOther soft spots: the SI with the 'full derivation' is referenced but not available, and no code or data are provided. The TST limiting case is true by construction—they group ρ·f·dV into P‡ and call it the effective population reaching the transition state. That's relabeling, not emergence.\n\nMy recommendation: if this is submitted as a research article claiming a unified predictive theory, reject. It could be acceptable as an invited perspective or viewpoint if the author drops the predictive claims and the misleading 'unified theory' title. A desk rejection is defensible because the internal inconsistency in the enzyme section is enough to sink the central claim.","headline":"Definitional framework with unvalidated 'predictions'; the enzyme example contradicts its own central integral.","tokens_in":12495,"tokens_out":3149,"would_cite":false,"duration_ms":35593,"reading_group":"no","serious_thinker":"no","would_accept_peer_review":false},"rs_alignment":null,"lean_confirmation":null,"pith_extraction":{"msc":[],"pacs":[],"model":"deepseek-v4-flash","headline":"Catalysis can be described by one spatiotemporal volume integral, with transition-state theory as a limiting case.","keywords":["active catalytic space","turnover frequency","turnover number","transition state theory","enzyme catalysis","homogeneous catalysis","heterogeneous catalysis","spatiotemporal fields"],"falsifier":"Measure turnover frequency for a catalytic surface whose site density is varied continuously, for example by diluting active sites with an inert spacer or changing surface coverage while holding the chemistry fixed. The ACS integral predicts TOF scales linearly with the integrated site density; if the per-site rate changes measurably with coverage or crowding, the prediction from Equation 3 will deviate from the measured rate. The size of that deviation directly quantifies the missing coupling terms that the exact multiplicative factorization omits.","tokens_in":11413,"feed_emoji":"🧪","tokens_out":10798,"duration_ms":91339,"temperature":0.7,"pith_summary":"This paper tries to establish that catalysis is one phenomenon, not three: enzymatic, homogeneous, and heterogeneous catalysis are all described by a single field equation for the turnover frequency (TOF), the number of reaction cycles per unit time. The equation is an integral over the catalytic volume of three coexisting fields: the density of active sites $\\rho(\\mathbf r)$, a dimensionless dynamic modulation $f(\\mathbf r,t)$, and the local per-site rate constant $k_{\\mathrm{local}}(\\mathbf r,t)$. An enzyme with gated conformational changes and a platinum surface with differently active crystal facets become two different profiles of the same integrand rather than two separate theories. The paper claims that, under the usual assumptions, the familiar transition-state-theory rate expression emerges from this integral as a special case, and that computed TOFs match known experimental ranges for a platinum single atom, platinum facets, an oxidase enzyme, and a ruthenium water-oxidation complex. A reader should care because, if correct, this gives a common quantitative language for comparing any catalysts and for plugging spatially resolved imaging data directly into rate predictions.","feed_headline":"Three catalysis fields collapse into one volume integral","feed_subtitle":"Enzymes, metal surfaces, and molecular catalysts share one equation; transition state theory is its limiting case.","key_machinery":"The load-bearing object is the active catalytic space: a volume $V \\subset \\mathbb R^3$ equipped with the three fields $\\rho(\\mathbf r)$, $f(\\mathbf r,t)$, and $k_{\\mathrm{local}}(\\mathbf r,t)$. The workhorse is the product identity $d\\mathrm{TOF} = \\rho\\, f\\, k_{\\mathrm{local}}\\, dV$, integrated over the volume to give $\\mathrm{TOF}(t)$ and once more over time to give the turnover number $\\mathrm{TON}(T)$. This integral does the unifying work: it lets an enzyme's allostery be encoded in $f$, a metal's crystal facets in a piecewise-constant $\\rho$ and $k_{\\mathrm{local}}$, and a solution's microheterogeneity in a Gaussian $\\rho$ with a decaying $k_{\\mathrm{local}}$. The derivation of transition state theory as a limit is made by localizing the integrand at one configuration $\\mathbf r^*$ and grouping $\\rho f\\,dV$ into the effective activated-state population $P^\\ddagger$, leaving $\\mathrm{TOF} = P^\\ddagger (k_{\\mathrm B}T/h)\\, e^{-\\Delta G^\\ddagger/RT}$.","core_discovery":"The central claim, stated on the paper's own terms, is that the instantaneous global turnover frequency of any catalyst is exactly $\\mathrm{TOF}(t) = \\int_V \\rho(\\mathbf r)\\, f(\\mathbf r,t)\\, k_{\\mathrm{local}}(\\mathbf r,t)\\, dV$, with $\\rho$ the local density of catalytically competent sites, $f$ a space- and time-dependent factor describing accessibility, conformation, or environmental gating, and $k_{\\mathrm{local}}$ the intrinsic rate of a single site. Each classical kinetic law is absorbed into one of these fields: saturating enzyme kinetics or electron-transfer rate expressions enter through $k_{\\mathrm{local}}$, structural data enter through $\\rho$, and fluctuating environments enter through $f$. The paper further asserts that transition state theory is the limiting case in which the integrand collapses to a single dominant configuration $\\mathbf r^*$, and that the energetic-span TOF$^\\circ$ formula used for benchmarking is derivable from this limit. In five worked systems, a single platinum atom, a two-facet platinum catalyst, bilirubin oxidase, a solvated ruthenium water-oxidation catalyst, and additional systems in the supplement, the framework yields TOF values that the paper reports as consistent with experimentally known ranges. The point of the framework is not to discard existing kinetics but to embed every existing rate law in a common spatiotemporal integral.","pith_inferences":["Editorial inference: The multiplicative form treats each infinitesimal volume element as an independent catalytic population, so the framework is likely to fail for strongly coupled systems such as substrate channeling or surface spillover; measuring TOF while continuously varying site density would reveal the missing coupling terms.","Editorial inference: The same structure could be applied to selectivity by promoting $k_{\\mathrm{local}}$ to a vector of branch rates, making product distributions volume integrals over the same $\\rho$ and $f$ fields, an extension the paper does not develop.","Editorial inference: Because the paper parameterizes $f$ with sinusoidal or gated forms, the framework can be tested dynamically by applying a modulated perturbation and comparing the predicted TOF waveform's phase and amplitude with time-resolved experimental turnover.","Editorial inference: With $\\rho$, $f$, and $k_{\\mathrm{local}}$ each independently measurable, the integrated equation acts as a consistency check on experimental data; large discrepancies between integrated and measured TOF would indicate missing cross-terms or a mis-assigned field."],"forward_implications":["Any catalyst can be scored by the same three fields, so an enzyme and an electrode can be compared directly in the same units rather than through field-specific metaphors.","Spatially resolved measurements, such as imaging, scanning probes, and operando spectra, can be inserted directly into the integral, turning qualitative maps of reactive hot spots into quantitative TOF predictions.","If the integral is exact, the standard transition-state and energetic-span formulas are diagnostic special cases, and deviations between them and the full integral identify when single-barrier thinking fails.","The formula decomposes a catalyst's performance into three independent levers, raising local site density, improving temporal accessibility, or increasing intrinsic per-site rate, each with a calculable effect on total TOF.","Because turnover number is defined as the time integral of the same expression, the framework connects instantaneous activity directly to long-term stability and deactivation behavior."],"supporting_citations":[{"why":"Introduces the active catalytic space notion that the paper formalizes into the field integral.","marker":"9"},{"why":"Supplies the standard turnover-frequency definitions whose energetic-span formula the paper derives as the transition-state limiting case.","marker":"10"},{"why":"Provides the multicopper oxidase system and its kinetic context used in the enzyme case study.","marker":"11"},{"why":"Provide the water-oxidation ruthenium complex and experimental turnover ranges used to benchmark the homogeneous case.","marker":"12,13"},{"why":"Provide platinum facet site densities, local rate constants, and edge-site reactivity used in the heterogeneous case.","marker":"14-17"}],"fun_headline_variants":["One volume integral unifies all three catalysis fields","From enzymes to metal surfaces: one TOF equation","Spatiotemporal integral folds catalysis into one framework","Unified catalysis: single integral predicts turnover frequency"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The load-bearing premise is that the three factors in $d\\mathrm{TOF}=\\rho\\, f\\, k_{\\mathrm{local}}\\, dV$ are independent and purely multiplicative, meaning the per-site rate $k_{\\mathrm{local}}$ is assumed not to change with local site density, crowding, or the modulation field; if real catalysts have site-to-site coupling or structure-dependent local rates, the volume integral does not equal the true turnover frequency.","fun_headline_variants_meta":{"raw":{"variants":["One volume integral unifies all three catalysis fields","From enzymes to metal surfaces: one TOF equation","Spatiotemporal integral folds catalysis into one framework","Unified catalysis: single integral predicts turnover frequency"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000218,"raw_usage":{"total_tokens":1427,"prompt_tokens":920,"completion_tokens":507,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":536,"completion_tokens_details":{"reasoning_tokens":446}},"tokens_in":536,"tokens_out":507,"duration_ms":5042,"temperature":1.0,"reasoning_tokens":446,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-16T04:24:14.214298+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Measure turnover frequency for a catalytic surface whose site density is varied continuously, for example by diluting active sites with an inert spacer or changing surface coverage while holding the chemistry fixed. The ACS integral predicts TOF scales linearly with the integrated site density; if the per-site rate changes measurably with coverage or crowding, the prediction from Equation 3 will deviate from the measured rate. The size of that deviation directly quantifies the missing coupling terms that the exact multiplicative factorization omits.","supporting_citations":[{"cited_title":"& Weckhuysen, B","cited_arxiv_id":null,"evidence_quote":"Introduces the active catalytic space notion that the paper formalizes into the field integral."},{"cited_title":"Turning over","cited_arxiv_id":null,"evidence_quote":"Supplies the standard turnover-frequency definitions whose energetic-span formula the paper derives as the transition-state limiting case."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Provides the multicopper oxidase system and its kinetic context used in the enzyme case study."}],"review_version":1}