{"id":"d83a86db-845f-4d4b-88e2-7bcf793ff004","arxiv_id":"1908.04294","paper_version":1,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":2.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":0,"one_line_summary":"Two standard formalisms for stochastizing one-step models, combinatorial state vectors and operator occupation numbers, are shown to give the same master equation, with a diagram technique for building the Liouville operator.","lead":"This paper lays out two equivalent ways to make one-step random processes into equations: working directly with state probabilities, or using creation and annihilation operators in Fock space. It is useful to modelers of populations, chemical reactions, and network traffic who need master or Langevin equations.","discovery_kind":"review","skeptic_critique":{"model":"deepseek-v4-flash","headline":"For multi-species schemes, Eq. (15) sums over species instead of taking products, so the claimed general Liouvillian does not reproduce the multivariate master equation.","rationale":"I focused on Eq. (15) rather than the reader's weakest assumption because the mass-action rate (9) is a modeling hypothesis common to both formalisms and does not threaten the claimed equivalence; it only limits the domain of applicability. The product/sum error is internal to the operator formalism: if left uncorrected, the operator method does not reproduce the combinatorial master equation for the interaction schemes defined in Section III. The reader did flag that Eq. (15) is asserted rather than proven; my concern makes that assertion concretely false in the multi-species case as printed. The Eq. (10) diffusion-coefficient sign error noted by the reader is real but localized to the Fokker-Planck reduction and does not affect the equivalence of the two master equations. Since the defect is fixable by replacing the sum with a product, or by restricting the paper to scalar one-step processes, the CONDITIONAL verdict remains appropriate; I do not move it.","tokens_in":15169,"tokens_out":12468,"duration_ms":132211,"concrete_test":"Apply Eq. (15) to the two-species annihilation scheme A+B->0 with I=(1,1), F=(0,0) and rate k n_A n_B, and compare the resulting operator action on |n_A,n_B> term by term with the master equation from Eq. (5) using s(n_A,n_B)=k n_A n_B. If the printed formula is used, the two disagree: Eq. (15) predicts independent death rates n_A and n_B, whereas the correct joint process has pair-annihilation rate n_A n_B. Replacing the sum over i in Eq. (15) by a product over i should restore exact agreement, distinguishing a typographical slip from a substantive defect.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The central claim requires Eq. (15) to be the Liouvillian of interaction scheme (1). As printed, Eq. (15) is L = sum over alpha,i of [ k+_alpha ((pi_i)^{F_i alpha} - (pi_i)^{I_i alpha}) (a_i)^{I_i alpha} + k-_alpha ((pi_i)^{I_i alpha} - (pi_i)^{F_i alpha}) (a_i)^{F_i alpha} ], i.e. a sum over species. This is correct only for a one-dimensional state space. In the multi-species setting of Section III, one reaction alpha changes several components phi_i simultaneously through the vector r_i alpha; its jump operator must be a product (tensor product) over the participating species, roughly k+_alpha (prod_i pi_i^{F_i alpha} - prod_i pi_i^{I_i alpha}) prod_i a_i^{I_i alpha}. The printed sum replaces one joint reaction by a sum of independent one-species reactions. Concretely, take A+B -> 0 with I=(1,1), F=(0,0) and rate k n_A n_B. The master equation from Eqs. (5) and (9) has pair-annihilation gain k(n_A+1)(n_B+1)p(n_A+1,n_B+1) and loss k n_A n_B p(n_A,n_B). The correct Liouvillian is k(1 - pi_A pi_B) a_A a_B. Eq. (15) as written gives k[(1 - pi_A)a_A + (1 - pi_B)a_B], whose master equation is two independent death processes, not pair annihilation. The Verhulst example is one-dimensional and cannot reveal this. Thus the general equivalence stated around Eq. (15) is not established as printed, independent of the mass-action assumption in Eq. (9).","agreement_with_reader":"partial"},"referee_report":{"model":"deepseek-v4-flash","summary":"The paper aims to present two formalisms for stochastizing one-step Markov models: the combinatorial approach based on the master equation for state vectors, and the operator approach based on creation and annihilation operators in an occupation-number (Fock) space. The authors derive transition rates from combinatorial counting, give Fokker–Planck and Langevin coefficients, propose a diagram technique for constructing the Liouville operator, and demonstrate the purported equivalence of the two formalisms on the Verhulst model. The central claim is that an arbitrary one-step interaction scheme (1) corresponds to the general Liouvillian operator (15), and that the combinatorial master equation and the operator Liouvillian produce identical dynamics.","tokens_in":15454,"tokens_out":3098,"duration_ms":35992,"significance":"If the claimed general equivalence were correct, the paper would provide a useful practical recipe for translating interaction schemes into stochastic master equations. The Verhulst example is computed carefully, and the algebraic check between Eq. (17) and Eq. (18) is correct, illustrating the operator method in one dimension. However, the general statement is not established as printed: Eq. (15) is incorrect for multi-species schemes, and the Fokker–Planck diffusion coefficient in Eq. (10) has a sign error. These defects are load-bearing because the paper's stated purpose is to provide a general stochastization method. The paper is essentially expository, summarizing ideas from the authors' earlier works, and would need correction before it can serve as a reliable reference.","major_comments":[{"comment":"The Liouvillian in Eq. (15) is written as a sum over species i, but a single reaction that simultaneously changes several species must be represented by a tensor product over those species. For the reaction A+B -> 0 with stoichiometry I_A=I_B=1 and F_A=F_B=0, Eq. (15) gives k[(1-pi_A)a_A + (1-pi_B)a_B], whose master equation describes two independent one-species death processes. The correct Liouvillian is k(1-pi_A pi_B)a_A a_B, corresponding to pair-annihilation transitions with rate k n_A n_B. This discrepancy is not visible in the one-dimensional Verhulst example, so the general equivalence claimed around Eq. (15) is not established.","section":"Section V, Eq. (15)"},{"comment":"The diffusion coefficient in Eq. (10) is stated as B_ij = r_iα r_jα [s+_fpα - s-_fpα]. The Kramers–Moyal expansion of the one-step master equation (5) yields B_ij = Σ_α r_iα r_jα [s+_fpα + s-_fpα], with a plus sign. With the printed minus sign, a pure death process (s- = 0) would have zero diffusion, which is incorrect. This error propagates through the relation (8) and the resulting Langevin equations, so the stochastic differential equations obtained from this framework are not reliable.","section":"Section IV, Eq. (10)"},{"comment":"The text explicitly says 'For simplicity, we consider the one-dimensional version' at the start of Section V, yet Eq. (15) is presented as the Liouvillian for the general scheme (1), which is n-dimensional. The paper should either prove a correctly product-ordered multi-species version of Eq. (15) or explicitly restrict its claimed equivalence to one-dimensional models. As written, the reader cannot determine the intended scope of the central claim.","section":"Section V, introductory paragraph and Eq. (15)"}],"minor_comments":[{"comment":"The notation I_i^j is introduced as an operator, but then used as a vector via I_j^α = I_i^j δ_i. This transition is confusing; define the diagonal reduction explicitly and consistently.","section":"Section II and Section III"},{"comment":"The falling-factorial expression ϕ_i!/(ϕ_i - I_i^α)! is only meaningful when ϕ_i ≥ I_i^α; the text should state that the rate is zero for ϕ_i < I_i^α.","section":"Eq. (9)"},{"comment":"In the second inner product, the summation index is denoted k while factorial moments are also denoted n_k; this is confusing and appears to be a typo. Clarify the index convention.","section":"Eq. (11)"},{"comment":"The diagrams in Figures 5 and 6 illustrate the subtraction of the denominator term, but the text does not define a precise rule for reading the multiplicity of lines in the normal-ordered product. A worked example with two identical incoming lines would help.","section":"Section V, diagram rules"},{"comment":"The manuscript contains a complete Russian translation after the English text. For an English-language journal submission, this should be removed or moved to supplementary material.","section":"General structure"}],"recommendation":"major_revision","confidential_remarks":"The paper is largely a summary of the authors' previous work, and the claimed novelty lies in the diagram technique. The central equivalence is not correct as printed for multi-species schemes, and the Fokker–Planck coefficient error further weakens the technical content. I believe the issues are fixable within the scope of the manuscript, so I recommend major revision rather than rejection. The paper also contains a full Russian translation, which is unusual for a journal submission and should be addressed editorially."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"Dear X,\n\nThe headline: the Verhulst example is correct, but the general Liouvillian (Eq. 15) is wrong for multi-species reactions and the Fokker-Planck diffusion term has a sign error. Those need fixing before the paper can be trusted as a general recipe.\n\nThe one-species Verhulst check is the best part: Eqs. (17) and (18) do coincide, and the calculation is straightforward and correct. The diagram notation is a reasonable mnemonic for single-species birth-death processes, and the authors correctly cite the original Doi-Peliti and Grassberger-Scheunert work. It is not a novel formalism; it is a compact restatement of known material.\n\nThe soft spots are significant. The big one is Eq. (15). As printed, it sums the Liouvillian over species. For a reaction that changes several components at once, the operator must be a product over the participating species, not a sum. Take A+B -> 0 at rate k n_A n_B. The correct Liouvillian is k(1 - pi_A pi_B) a_A a_B. Eq. (15) gives k(1 - pi_A) a_A + k(1 - pi_B) a_B, which describes two independent death processes. So the claimed general correspondence between interaction scheme (1) and Eq. (15) is not established; it is false for multi-species reactions. The Verhulst example is one-dimensional and cannot expose this.\n\nSecond, Eq. (10) gives the Fokker-Planck diffusion coefficient as r_i r_j (s+ - s-). The Kramers-Moyal expansion of Eq. (5) gives r_i r_j (s+ + s-). The sign error makes the Langevin/Fokker-Planck part unreliable as stated.\n\nThird, Eq. (9) is a mass-action assumption, not a first-principles derivation. That is acceptable, but it should be presented as an assumption because the whole construction inherits it.\n\nThe paper is a tutorial, not a research advance. If you work on single-species birth-death processes, the Verhulst example is a useful worked exercise. But the general claims need corrections before the paper can be trusted.\n\nI would send this to peer review, not desk-reject, because the errors are concrete and fixable, and a corrected version would be a readable consolidation. It needs major revision before publication.\n\nBest,","headline":"The Verhulst example works, but the general Liouvillian (Eq. 15) is wrong for multi-species reactions and the Fokker-Planck diffusion term has a sign error.","tokens_in":16095,"tokens_out":6357,"would_cite":false,"duration_ms":53521,"reading_group":"no","serious_thinker":"yes","would_accept_peer_review":true},"rs_alignment":null,"lean_confirmation":null,"pith_extraction":{"msc":["60J27","82C31"],"pacs":["05.40.-a","02.50.Ey"],"model":"deepseek-v4-flash","headline":"For one-step processes, two stochastic formalisms give the same master equation.","keywords":["stochastization","one-step processes","master equation","occupation numbers representation","Fock space","Dirac notation","Liouville operator","diagram technique"],"falsifier":"Apply both formalisms to a two-species one-step scheme such as $A+B \\rightleftharpoons C$ with rates $k n_A n_B$ and $k' n_C$ as prescribed by Eq. (9), and compare the master equation obtained from Eq. (5) with the one obtained from the Liouvillian (15). The paper's claim implies they agree term by term; a single mismatch in the two-species cross terms or in the ordering of $\\pi_A a_A$ and $\\pi_B a_B$ would refute the general recipe.","tokens_in":14888,"feed_emoji":"🎲","tokens_out":12433,"duration_ms":119009,"temperature":0.7,"pith_summary":"The paper tries to establish that two routes for making a one-step model stochastic—the combinatorial state-vector formalism and the operator occupation-number formalism—lead to the same master equation. It provides a general rule that converts any interaction scheme into a Liouville operator, plus a diagram technique that produces the same operator directly from the reaction drawings. The Verhulst (logistic growth) model is worked out both ways, and the two master equations coincide term by term. If the claim is right, a modeler can choose either formalism for convenience and still describe the same stochastic process.","feed_headline":"One-step processes: two formalisms, same master equation","feed_subtitle":"Combinatorial state vectors and occupation-number operators agree, giving a diagram recipe for stochastic models.","key_machinery":"The load-bearing object is the Liouville operator of Eq. (15), assembled from creation operators $\\pi_i$ and annihilation operators $a_i$ with powers given by the initial and final stoichiometric coefficients of each reaction. It is built to satisfy the probability-conservation condition $\\langle 0|L = 0$, which fixes the subtraction of the interaction terms. Its companion is the rate ansatz of Eq. (9): a reaction runs at the rate constant times the number of ordered arrangements of the reactants, i.e., a falling factorial. Together these turn an interaction scheme into either a state-space master equation or a Fock-space Liouvillian, and the paper's diagram technique reads the Liouvillian off the diagrams without returning to the state space.","core_discovery":"Under the falling-factorial transition rates of Eq. (9), an arbitrary one-step interaction scheme $I_\\alpha \\rightleftharpoons F_\\alpha$ is shown to yield one and the same master equation whether it is derived combinatorially from Eq. (5) or from the Liouville operator $L = \\sum_{\\alpha,i} \\bigl[ k^+_\\alpha ((\\pi_i)^{F_{i\\alpha}} - (\\pi_i)^{I_{i\\alpha}})(a_i)^{I_{i\\alpha}} + k^-_\\alpha ((\\pi_i)^{I_{i\\alpha}} - (\\pi_i)^{F_{i\\alpha}})(a_i)^{F_{i\\alpha}} \\bigr]$, where $\\pi_i$ creates and $a_i$ annihilates a particle of species $i$. The paper's demonstration is the Verhulst model: the combinatorial master equation (17), with birth rate $\\lambda n$, death rate $\\beta n$, and competition rate $\\gamma n(n-1)$, is exactly the master equation (18) recovered from the Liouvillian. The diagrammatic rules of Section V reconstruct the same Liouvillian, so the operator approach is a transcription of the same stochastic process, not an approximation of it.","pith_inferences":["The paper leaves implicit that the agreement between formalisms is built into the shared falling-factorial rate law, so a mismatch with real data would implicate the rate ansatz rather than the choice of formalism.","A natural extension not implemented here is to apply the diagram rules to a multi-species one-step reaction such as $A+B \\rightleftharpoons C$ and check the cross terms term by term.","A further testable consequence is that rates not expressible as finite falling factorials, such as saturating or Hill-type rates, fall outside both formalisms as stated and would require extending Eq. (15)."],"forward_implications":["Any one-step model written as an interaction scheme can be converted into a master equation, a Fokker–Planck equation, or a Langevin equation without ad-hoc guessing.","The two formalisms agree for every scheme that obeys the falling-factorial rate law, so a discrepancy between them would indicate a mistake in applying the rules rather than a difference between the formalisms.","The diagram technique gives a direct route to the Liouvillian: incoming lines become annihilation operators, outgoing lines become creation operators, and the interaction line carries the rate constant.","For the Verhulst model, the resulting stochastic process has birth rate $\\lambda n$, death rate $\\beta n$, and competition rate $\\gamma n(n-1)$."],"supporting_citations":[{"why":"It supplies the Fock-space formalism for identical classical objects that the operator approach adopts.","marker":"[7]"},{"why":"It introduces second quantization for classical many-particle systems, the basis of the occupation-number Liouvillian.","marker":"[16]"},{"why":"It applies the same occupation-number method to diffusion-controlled reactions, supporting the operator route.","marker":"[17]"},{"why":"It extends the operator formalism to path integrals for birth-death processes on a lattice.","marker":"[18]"},{"why":"It provides the master-equation framework that the combinatorial approach starts from.","marker":"[14]"},{"why":"It supplies the Kramers–Moyal expansion used to pass from the master equation to Fokker–Planck and Langevin forms.","marker":"[15]"},{"why":"It is the original Verhulst logistic model used as the worked comparison example.","marker":"[19]"}],"fun_headline_variants":["Two formalisms, one master equation for one-step models","Combinatorial and operator approaches agree exactly","Same stochastic equation from state vectors or operators","Stochastization dual paths converge in one-step case","Diagram rules prove equivalence of two formalisms"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The load-bearing premise is that each reaction's probability per unit time equals a rate constant times the number of ordered ways to pick the reacting individuals from the current population; if real rates depend on the state in any other way, both formalisms describe the wrong stochastic process even though they still agree with each other.","fun_headline_variants_meta":{"raw":{"variants":["Two formalisms, one master equation for one-step models","Combinatorial and operator approaches agree exactly","Same stochastic equation from state vectors or operators","Stochastization dual paths converge in one-step case","Diagram rules prove equivalence of two formalisms"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000299,"raw_usage":{"total_tokens":1682,"prompt_tokens":849,"completion_tokens":833,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":465,"completion_tokens_details":{"reasoning_tokens":761}},"tokens_in":465,"tokens_out":833,"duration_ms":9212,"temperature":1.0,"reasoning_tokens":761,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-14T13:58:06.784921+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Apply both formalisms to a two-species one-step scheme such as $A+B \\rightleftharpoons C$ with rates $k n_A n_B$ and $k' n_C$ as prescribed by Eq. (9), and compare the master equation obtained from Eq. (5) with the one obtained from the Liouvillian (15). The paper's claim implies they agree term by term; a single mismatch in the two-species cross terms or in the ordering of $\\pi_A a_A$ and $\\pi_B a_B$ would refute the general recipe.","supporting_citations":[{"cited_title":"Mollison, Dependence of epidemic and population velo cities on basic parameters, Mathematical Biosciences 107 ( 2) (1991) 255–287","cited_arxiv_id":null,"evidence_quote":"It supplies the Fock-space formalism for identical classical objects that the operator approach adopts."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"It introduces second quantization for classical many-particle systems, the basis of the occupation-number Liouvillian."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"It applies the same occupation-number method to diffusion-controlled reactions, supporting the operator route."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"It extends the operator formalism to path integrals for birth-death processes on a lattice."},{"cited_title":"Waage, C","cited_arxiv_id":null,"evidence_quote":"It supplies the Kramers–Moyal expansion used to pass from the master equation to Fokker–Planck and Langevin forms."}],"review_version":1}