{"id":"1cec735f-1160-4aa8-aaee-0d1ef6d0e338","arxiv_id":"1908.07764","paper_version":1,"verdict":"ACCEPT","confidence":"HIGH","novelty_score":7.0,"correctness_risk":"low","formal_verification":"none","parameter_count":0,"one_line_summary":"For molecules interacting through an isotropic pair potential, the Doi volume reaction rate constant is given by a semi-analytic formula requiring only a Debye-type integral and a one-dimensional boundary-value problem.","lead":"This paper derives a mathematical formula that connects the microscopic reaction parameters used in particle-based reaction-diffusion simulations to the macroscopic reaction rate when molecules attract or repel each other. The result lets researchers calibrate simulation models against measured reaction rates, which matters for building quantitative models of chemistry inside cells.","discovery_kind":"extension","skeptic_critique":{"model":"deepseek-v4-flash","headline":"Singular-potential boundary condition at r=0 is not rigorously established, and the §IVD logarithmic check meant to support it is internally inconsistent.","rationale":"I agree with the reader that the weakest assumption is the r=0 boundary condition for divergent potentials, and I do not accuse the authors of any error beyond an internal inconsistency in the validation example. The central formula Eq. (32) is cleanly derived for regular potentials and is independently supported by the iPRD simulations, which is why I would not reject the paper. However, the manuscript states a logarithmic-potential validation that cannot be correct as written: for the quoted potential, g(r) is R^2/r^4 inside the reaction volume, not a step function, and the stated coefficient −4/r corresponds to that same potential while the step form would require the opposite sign. This means the singular BVP has not actually been validated by the analytic check shown in Fig. 2, and the extension to the Lennard-Jones case rests on a physical bound plus finite-difference convergence that is plausible but unproven. Since the LJ example is a central demonstration and the simulations lack error bars, I would make acceptance conditional on one independent check of the singular BVP rather than rejecting the work.","tokens_in":20637,"tokens_out":20878,"duration_ms":206379,"concrete_test":"For the Lennard-Jones potential (47) at the parameters of Table I, independently solve Eq. (29) with a method that treats r=0 as an irregular singular point, e.g., shooting from r=R while imposing the physical asymptotic ψ(r)=O(r^{m+2}e^{−βU(r)}) as r→0, or solving in the variable w=e^{βU}p. Compare ψ′(R) and the resulting k from Eq. (32) against the finite-difference result with h=R/10^4 for κR = 0.1, 1, and 10. If the relative difference exceeds 10^{−3}, the singular-potential extension of Eq. (32) is not verified; if it agrees, the concern is resolved.","verdict_should_be":"CONDITIONAL","load_bearing_attack":"The central result Eq. (32) depends on the inner boundary-value problem (29), specifically on the condition ψ(0)=0. For finite potentials this follows from Gauss' theorem via Eqs. (12)–(13). For potentials with U(r)→∞ at r=0, the paper excludes the origin and instead appeals to the bound p(r) ≤ c_A e^{−βU(r)} (Section III after Eq. (15); used in Section IVD). This bound is physically plausible and probably correct, but it is not proven, and it is load-bearing for the Lennard-Jones example, where the potential diverges at the origin. Moreover, the numerical validation intended to cover this regime is not self-consistent as printed: for the logarithmic potential in Eq. (42), U(r)=−2k_BT ln(r/R), a direct computation gives g(r)=e^{βU}/r^2 = R^2/r^4, not R^{−2}θ(R−r); obtaining the step form would require the opposite sign of U, which in turn gives βU′−2/r = 0 rather than −4/r. Thus the claimed analytic check does not clearly test the singular BVP used for the Lennard-Jones potential, and convergence for algebraic divergences (m=12) is asserted but not demonstrated.","agreement_with_reader":"partial"},"referee_report":{"model":"deepseek-v4-flash","summary":"The manuscript extends Doi's volume-reaction model to bimolecular reactions in which the two reactants interact through an arbitrary isotropic pair potential U(r). In the steady state, the authors derive a semi-analytical expression, Eq. (32), for the macroscopic association rate constant k in terms of the reactivity κ = √(λ/D), the reaction radius R, and the potential: k = 4πD [∫_R^∞ e^{βU(s)} s^{-2} ds + e^{βU(R)} R^{-2} ψ'(R)/κ^2]^{-1}, where ψ solves the one-dimensional boundary-value problem (29) with ψ(0)=0 and ψ(R)=1. The expression reproduces Debye's rate in the diffusion-limited limit and λV_eff in the well-mixed limit, and it decomposes uniquely into encounter and formation rates. The theory is checked against an analytically solvable logarithmic potential, against perturbation theory for slow reactions, and against 195 stochastic iPRD simulations for non-interacting, softly repulsive, and Lennard-Jones systems.","tokens_in":20849,"tokens_out":20091,"duration_ms":254111,"significance":"If the final result holds, this is a valuable and practical contribution. It bridges the well-mixed and diffusion-limited regimes for interacting particles and provides a direct route to calibrate the microscopic parameters λ and R in iPRD simulations. The derivation of Eq. (32) is clean and self-contained; the limiting cases correctly recover Smoluchowski, Debye, and Doi results; and the decomposition k^{-1} = k_e^{-1} + k_f^{-1} emerges naturally from the matching calculation. The perturbative result, Eq. (40), and the cheap finite-difference scheme with documented linear convergence are useful tools. The simulation campaign is extensive, and the agreement between theory and simulations in Figs. 4 and 8 is excellent. The paper also makes falsifiable predictions, such as the non-monotonic dependence of the formation rate on the reaction-boundary position in Fig. 5.","major_comments":[],"minor_comments":[{"comment":"The sentence \"With this, g(r) = R^{-2} θ(R−r) is a step function\" is incorrect. For the stated potential U(r) = −2k_BT ln(r/R) inside r<R, one has e^{βU} = R^2/r^2 and hence g(r) = R^2/r^4 inside, while g(r) = 1/r^2 outside. The analytic solution Eq. (43) is nevertheless the correct solution for this potential (for example, its slow-reaction limit matches the perturbative result Eq. (40) with V_eff = 4πR^3/5), so the numerical validation is not invalidated; however, the misleading statement about g should be corrected.","section":"IV.D, Eq. (42)"},{"comment":"The extension of the boundary condition j(0)=0 to potentials that diverge at r=0 rests on the bound p(r) ≤ c_A e^{−βU(r)}, which is physically motivated but not proven. Since this bound is load-bearing for the Lennard-Jones example via Eq. (29b) and for the boundedness argument in Section IV.D, the authors should either supply a short proof (e.g., by a maximum principle or a probabilistic killing argument) or explicitly state the bound as a regularity assumption on the potential.","section":"III, after Eq. (15); IV.D"},{"comment":"In the sentence \"The expression is proportional to [βU'(r) − 2/r] r^2 p(r)\", the proportionality constant is 4πκ^2D/K; making this explicit would help the reader verify the limiting argument that follows.","section":"IV.D, boundedness argument"},{"comment":"Reference 64 contains a typo in the journal name: \"Chem. Phys.ll\" should read \"Chem. Phys.\".","section":"References"}],"recommendation":"minor_revision","confidential_remarks":"The manuscript is technically sound and the simulations are convincing. The only substantive point I would raise to the authors is the status of the singular-potential boundary condition; as printed it is a physically argued assumption, and the authors should make that status explicit. I do not think the incorrect description of the logarithmic test (the g step statement) undermines the numerical validation, since Eq. (43) is consistent with the actual potential defined in Eq. (42)."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"The paper delivers a genuine extension of Doi's volume-reaction model to isotropic pair potentials, with Eq. (32) as a compact semi-analytic rate formula that interpolates cleanly between the well-mixed and diffusion-limited regimes. The derivation is transparent: outside the reaction radius it reduces to Debye, the inner ODE for ψ is exact, and the matching at R is continuous. The iPRD simulations are extensive and give good agreement for the overall rates and concentration profiles. I think the reader's ACCEPT verdict is basically right.\n\nBut the stress-test note is also right, and the reader missed it. The logarithmic potential check in Section IVD is internally inconsistent. The paper states that for U(r) = −2k_BT ln(r/R) (r<R), g(r)=R^{-2}θ(R−r). Direct computation gives g(r)=e^{βU}/r^2 = (R/r)^2 r^{-2} = R^2/r^4, not a step function. Getting the claimed step requires the opposite sign of U, which would change the ODE coefficient from −4/r to 0. So the analytic check used to validate the singular BVP for the Lennard-Jones potential does not actually test that BVP. The convergence plot in Fig. 2 is therefore not a valid test of the singular case, and the assertion of convergence for algebraic divergences (m=12) is not demonstrated.\n\nThat matters, because the boundary condition ψ(0)=0 for potentials diverging at the origin is an assumption, not a proven consequence of Gauss' theorem. The paper is transparent about this—it appeals to a physical bound p(r) ≤ c_A e^{−βU(r)}—but the one analytic check designed to support it is flawed. The LJ simulation agreement is encouraging evidence that the assumption is right, but it is weaker evidence than a correct singular analytic test would be, because the near-origin region is rarely sampled. The main formula Eq. (32) for finite potentials (like the harmonic repulsion) rests on rigorous footing, so the core result stands.\n\nOther soft spots are minor: no error bars in the simulation figures, and no code or data release. Self-citations are fine given the context.\n\nWho is this for? Practitioners of particle-based reaction–diffusion (iPRD) simulations who need to calibrate λ and R against macroscopic rate constants, and anyone working on diffusion-influenced reactions with pair interactions. It deserves a serious referee. Recommended course: engage the paper, but ask the authors to fix the logarithmic validation—either correct the sign in Eq. (42) and update the ODE and Eq. (43) accordingly, or replace it with a genuinely singular test case—and to provide error bars and data/code if possible.","headline":"Solid central result with a flawed validation detail: the logarithmic potential check in §IVD does not hold up, but Eq. (32) is still a worthwhile contribution.","tokens_in":21382,"tokens_out":5676,"would_cite":true,"duration_ms":96446,"reading_group":"yes","serious_thinker":"yes","would_accept_peer_review":true},"rs_alignment":null,"lean_confirmation":null,"pith_extraction":{"msc":[],"pacs":[],"model":"deepseek-v4-flash","headline":"A single formula now gives the reaction rate constant for diffusing molecules that attract or repel each other, valid from well-mixed to diffusion-limited kinetics.","keywords":["diffusion-influenced reactions","Doi volume reaction model","pair interaction potential","reaction rate constant","encounter rate","formation rate","reaction-diffusion equation","particle-based reaction-diffusion simulations"],"falsifier":"Measure the radial concentration profile $p(r)$ near the origin in a stochastic Brownian-dynamics simulation of a hard-core or strongly repulsive potential, e.g. $U(r) \\sim a r^{-12}$, with finite propensity $\\lambda$. If $p(r)$ exceeds $c_A e^{-\\beta U(r)}$ anywhere, the equilibrium upper bound used to extend $j(0)=0$ to divergent potentials is wrong and Eq. (32) would have to be revised. Alternatively, compare the theory's prediction for the diffusion-limited plateau, $k \\to 4\\pi D[\\int_R^\\infty e^{\\beta U(s)}s^{-2} ds]^{-1}$, against brute-force simulation at $\\kappa R = 10^3$ for a potential with a deep well inside the reaction sphere.","tokens_in":20440,"feed_emoji":"🧪","tokens_out":7313,"duration_ms":62697,"temperature":0.7,"pith_summary":"Chemical reactions in solution require two molecules to first meet by diffusion, and the forces between them—steric repulsion, electrostatic attraction—change how often they meet. This paper extends Doi's volume-reaction model, in which a pair reacts with propensity $\\lambda$ whenever their centres are closer than a reaction radius $R$, to molecules that interact through an arbitrary isotropic pair potential $U(r)$. It derives a closed, semi-analytical formula, Eq. (32), for the steady-state macroscopic rate constant $k$ in terms of $\\lambda$, $R$, the relative diffusion constant $D$, and $U$, valid across the whole range from well-mixed (reaction-limited) to diffusion-limited kinetics. The rate decomposes exactly into an encounter rate and a formation rate, reproducing Smoluchowski's and Debye's classical limits, and the predictions match extensive stochastic particle-based simulations. This gives a practical route to calibrate microscopic reaction parameters to experimental rate constants.","feed_headline":"One formula links reaction rates to intermolecular forces","feed_subtitle":"Extends the Doi reaction model to attracting and repelling molecules, from well-mixed to diffusion-limited kinetics.","key_machinery":"The load-bearing object is $\\psi(r)$, the dimensionless radial flux profile inside the reaction volume: it removes the unknown reaction frequency $K$ from the inner problem and turns the reaction–diffusion equation into a linear, one-dimensional boundary-value problem,\n$$\\psi'' + (\\$\\beta$ U' - 2/r)\\psi' - (\\$\\lambda$/D)\\psi = 0, \\quad \\psi(0)=0,\\ \\psi(R)=1.$$\nIts boundary derivative $\\psi'(R)$ fixes the concentration at the reaction surface, $p(R) = \\psi'(R)K/(4\\pi R^2\\lambda)$, and thereby the formation rate. Outside the sphere, flux conservation plus the Boltzmann-weighted gradient law give the Debye profile in terms of the single integral of $g(s)=e^{\\beta U(s)} s^{-2}$. Matching the two solutions at $r=R$ produces Eq. (32), so the whole argument is powered by this two-domain split and the linearity of the inner equation.","core_discovery":"The central result is that the steady-state reaction rate constant for an isotropic pair potential $U(r)$ is\n$$k = 4\\pi D\\left[\\int_R^\\infty $e^{{\\beta U(s)}}$ $s^{{-2}}$\\,ds + \\frac{$e^{{\\beta U(R)}}$ \\psi'(R)}{\\$\\lambda$/D}\\right]^{-1},$$\nwhere $\\psi(r)$, defined by $\\psi = -4\\pi r^2 j(r)/K$ inside the reactive sphere, solves $\\psi'' + (\\beta U' - 2/r)\\psi' - (\\lambda/D)\\psi = 0$ with $\\psi(0)=0$ and $\\psi(R)=1$; here $j(r)$ is the radial flux and $K$ is the reaction frequency. The first term in the bracket is the Debye encounter rate, the second is the formation rate, so the formula is an exact harmonic-mean decomposition $k^{-1}=k_e^{-1}+k_f^{-1}$. As the propensity $\\lambda$ grows, $k$ saturates at the Debye rate, and without any potential it reduces to Doi's and Smoluchowski's results. The same construction yields the full concentration profile $p(r)$, and a slow-reaction expansion gives $k \\approx \\lambda\\int_{|r|\\le R} e^{-\\beta U(r)} d^3r$, so the rate is the propensity times the Boltzmann-weighted accessible reaction volume. These predictions were compared with particle-based simulations for a repulsive harmonic potential, a Lennard-Jones potential, and no potential, with agreement across three decades of reactivity.","pith_inferences":["Because the inner equation is linear, the same matching logic should extend to time-dependent rates: a Laplace transform in time would turn Eq. (29) into a resolvent problem, giving the full autocorrelation of the rate, not just its steady-state value.","The non-monotonic formation rate seen when the reaction boundary moves through the potential well suggests that the reaction radius could be treated as a tunable design parameter: choosing $R$ near the potential minimum should maximize the total rate at fixed $\\lambda$.","In a crowded environment the relevant potential is a potential of mean force; if that quantity were measured or computed, Eq. (32) would predict how crowding shifts the apparent rate, an extension the paper explicitly leaves to future work.","Inverting Eq. (32) numerically yields $(\\lambda, R)$ pairs that reproduce a target experimental $k$; one could tabulate $\\lambda$ for a grid of $R$ and $U$, producing ready-to-use parameter tables for particle-based reaction–diffusion simulations."],"forward_implications":["For any isotropic potential, computing $k$ reduces to one quadrature and one one-dimensional boundary-value problem, so the rate can be evaluated essentially instantly for any choice of $\\lambda$ and $R$.","The exact decomposition $k^{-1}=k_e^{-1}+k_f^{-1}$ holds over the full reactivity range, with $k_e$ always Debye's rate and $k_f$ fixed by the boundary concentration; the diffusion-limited limit is the Debye/Smoluchowski rate.","A purely repulsive potential inside the reaction volume lowers both partial rates and is most damaging in the slow-reaction regime, while attraction outside the sphere raises the encounter rate and can raise the total rate by tens of percent.","In the well-mixed limit the rate reduces to $k = \\lambda V_{\\mathrm{eff}}$, with $V_{\\mathrm{eff}}$ the Boltzmann-weighted reaction volume, giving a direct way to measure accessible volume from reaction kinetics.","The matching of theory and particle-based simulation across three decades of $\\lambda$ means the model can be used to calibrate microscopic parameters to experimental rate constants without free fit parameters."],"supporting_citations":[{"why":"gives the potential-dependent encounter rate (Debye rate) that provides the outer-domain term in Eqs. (22) and (32).","marker":"[60]"},{"why":"introduces the propensity-based volume reaction model that this work extends to pair potentials.","marker":"[63]"},{"why":"defines the Doi volume-reaction model whose rate Eq. (26) is recovered in the non-interacting limit.","marker":"[64]"},{"why":"supplies the closed-form solution of the inner reacting-sphere problem used for the potential-free case.","marker":"[37]"},{"why":"provides the classical Smoluchowski diffusion-limited rate $4\\pi DR$ recovered as the zero-potential limit.","marker":"[67]"},{"why":"supplies the stochastic particle-based simulation code used for validation of the theory.","marker":"[54]"},{"why":"established the interacting-particle reaction-dynamics setting and the effective reaction volume used in the slow-reaction limit.","marker":"[48]"},{"why":"states the harmonic-mean composition of rates underlying the encounter/formation decomposition.","marker":"[58]"}],"fun_headline_variants":["One equation solves reaction rates for pair potentials","Pair forces dictate reaction rates from well-mixed to diffusion-limited","Exact steady-state rates for diffusive reactive pairs","Unified rate constant for interacting molecules"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The derivation requires the no-flux condition at the origin, $j(0)=0$, and for potentials that diverge at $r=0$ this condition is justified by an assumed physical bound on the concentration rather than a proof; if the near-origin behavior were different, Eq. (32) would change.","fun_headline_variants_meta":{"raw":{"variants":["One equation solves reaction rates for pair potentials","Pair forces dictate reaction rates from well-mixed to diffusion-limited","Exact steady-state rates for diffusive reactive pairs","Unified rate constant for interacting molecules"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.001296,"raw_usage":{"total_tokens":5361,"prompt_tokens":1087,"completion_tokens":4274,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":703,"completion_tokens_details":{"reasoning_tokens":4214}},"tokens_in":703,"tokens_out":4274,"duration_ms":33286,"temperature":1.0,"reasoning_tokens":4214,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-14T11:58:04.331057+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Measure the radial concentration profile $p(r)$ near the origin in a stochastic Brownian-dynamics simulation of a hard-core or strongly repulsive potential, e.g. $U(r) \\sim a r^{-12}$, with finite propensity $\\lambda$. If $p(r)$ exceeds $c_A e^{-\\beta U(r)}$ anywhere, the equilibrium upper bound used to extend $j(0)=0$ to divergent potentials is wrong and Eq. (32) would have to be revised. Alternatively, compare the theory's prediction for the diffusion-limited plateau, $k \\to 4\\pi D[\\int_R^\\infty e^{\\beta U(s)}s^{-2} ds]^{-1}$, against brute-force simulation at $\\kappa R = 10^3$ for a potential with a deep well inside the reaction sphere.","supporting_citations":[{"cited_title":"Debye ,\\ title title Reaction Rates in Ionic Solutions , \\ 10.1149/1.3071413 journal journal J","cited_arxiv_id":null,"evidence_quote":"gives the potential-dependent encounter rate (Debye rate) that provides the outer-domain term in Eqs. (22) and (32)."},{"cited_title":"Teramoto \\ and\\ author N","cited_arxiv_id":null,"evidence_quote":"introduces the propensity-based volume reaction model that this work extends to pair potentials."},{"cited_title":"Doi ,\\ title title Theory of diffusion-controlled reactions between non-simple molecules","cited_arxiv_id":null,"evidence_quote":"defines the Doi volume-reaction model whose rate Eq. (26) is recovered in the non-interacting limit."},{"cited_title":"von Smoluchowski ,\\ title title Versuch einer mathematischen Theorie der Koagulationskinetik kolloider Lösungen , \\ 10.1007/BF01427232 journal journal Z","cited_arxiv_id":null,"evidence_quote":"provides the classical Smoluchowski diffusion-limited rate $4\\pi DR$ recovered as the zero-potential limit."},{"cited_title":"Hoffmann , author C","cited_arxiv_id":null,"evidence_quote":"supplies the stochastic particle-based simulation code used for validation of the theory."},{"cited_title":"Fr \\\"o hner \\ and\\ author F","cited_arxiv_id":null,"evidence_quote":"established the interacting-particle reaction-dynamics setting and the effective reaction volume used in the slow-reaction limit."},{"cited_title":"Shoup \\ and\\ author A","cited_arxiv_id":null,"evidence_quote":"states the harmonic-mean composition of rates underlying the encounter/formation decomposition."}],"review_version":1}