REVIEW 3 major objections 5 minor 42 references
Two formalisms of stochastization of one-step models
T0 review · 3 major / 5 minor · reviewed 2026-08-14 · deepseek-v4-flash
Pith's one-line read For one-step processes, two stochastic formalisms give the same master equation.
desk verdict 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. 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 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.
What would settle it
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.
Extended reading notes
Core claim
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.
Load-bearing premise
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.
Editorial extensions
If this is right
- 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)$.
Reading between the lines
- 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).
Signed reviews
Editorial analysis
A structured set of objections, weighed in public.
Referee Report
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.
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 (3)
- [Section V, Eq. (15)] 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 IV, Eq. (10)] 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 V, introductory paragraph and Eq. (15)] 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.
minor comments (5)
- [Section II and Section III] 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.
- [Eq. (9)] 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^α.
- [Eq. (11)] 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 V, diagram rules] 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.
- [General structure] 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.
Circularity Check
No significant circularity: the operator/combinatorial equivalence is verified by direct calculation, not fitted or self-citation-dependent.
full rationale
The paper's central equivalence claim is that the combinatorial master equation and the operator Liouvillian produce the same one-step stochastic process. This is not a fitted prediction: the transition rates are specified once in Eq. (9) and then used identically in both formalisms. The operator Liouvillian (15) is constructed in the text from the reaction scheme (1) by explicit diagrammatic rules, so the earlier self-citations [5,6] are contextual rather than load-bearing. The Verhulst example is a genuine internal check: Eq. (18) is obtained by acting with the Liouvillian on number states, and the paper states that it 'completely coincides with the formula (17), obtained by the combinatorial method.' No parameter is fitted to a subset of data and then renamed a prediction, and no conclusion is defined in terms of the result it is supposed to derive. The mass-action form of Eq. (9) is an input assumption, not an output of the derivation; questioning its first-principles status is a modelling critique, not a circularity critique. Similarly, the multi-species issue raised about the printed sum in Eq. (15) is a mathematical correctness concern, not an instance of circular reasoning. Because the derivation is self-contained and the claimed equivalence is verified by explicit calculation, the circularity score is 0.
Assumptions & free parameters
assumptions (5)
- domain assumption Transition rates are falling factorials of reactant counts: s+/-_alpha = k+/-_alpha times the product over i of phi_i!/(phi_i - I/F_i_alpha)! (Eq. 9).
- domain assumption The Fokker-Planck equation is obtained by truncating the Kramers-Moyal expansion at second order.
- domain assumption A classical one-step process can be represented in Fock space with Bose commutation [a,pi]=1 and the exclusive inner product (11).
- domain assumption Probability conservation of the Liouville operator, <0|L = 0 (Eq. 13), is used to fix the subtraction term in the diagrams.
- domain assumption The operator treatment is one-dimensional and spatially homogeneous.
Cite this review
Pith. "Pith review of Two formalisms of stochastization of one-step models." pith.science (2026). https://pith.science/paper/HTRQA6WV
@misc{pith2026190804294,
author = {Pith},
title = {Pith review of: Two formalisms of stochastization of one-step models},
year = {2026},
howpublished = {\url{https://pith.science/paper/HTRQA6WV}},
note = {Machine review of arXiv:1908.04294}
}
read the original abstract
To construct realistic mathematical models from the first principles, the authors suggest using the stochastization method. In a number of works different approaches to stochastization of mathematical models were considered. In the end, the whole variety of approaches was reduced to two formalisms: combinatorial (state vectors) and operator (occupation numbers). In the article the authors briefly describe these formalisms with an emphasis on their practical application.
Figures
Figures from the paper (1 more)
Reference graph
Works this paper leans on
-
[1]
The abstract indices notation (see [11]) is used in this wo rk. Under this notation a tensor as a whole object is denoted just as an index (e.g., xi), components are denoted by underlined index (e.g., xi )
-
[2]
We will adhere to the following agreements. Latin indices from the middle of the alphabet ( i, j, k) will be applied to the space of the system state vectors. Latin indic es from the beginning of the alphabet ( a) will be related to the Wiener process space. Greek indices ( α) will set a number of different interactions in kinetic equations. III. INTERACTI...
-
[3]
RU DN University Program 5-100
· · ·(ϕ − (n − 1)) will be replaced by (ϕ)n: s+ fp α = k+ α n∏ i =1 (ϕi )I i α , s− fp α = k− α n∏ i =1 (ϕi )F i α . Then for the Fokker–Planck equation (6) we may obtain formul as for the coefficients: Ai :=Ai(ϕk) = riα [ s+ fp α − s− fp α ] , Bij :=Bij(ϕk) = riα rjα [ s+ fp α − s− fp α ] . (10) Using the relation (8), we may obtain the coefficients for the ...
-
[6]
W. Feller, Die Grundlagen der Volterraschen Theorie des Kampfes ums Dasein in wahrscheinlichkeitstheoretischer B ehand- lung, Acta Biotheoretica 5 (1) (1939) 11–40. doi:10.1007/BF01602932
-
[7]
D. Mollison, Dependence of epidemic and population velo cities on basic parameters, Mathematical Biosciences 107 ( 2) (1991) 255–287. doi:10.1016/0025-5564(91)90009-8
-
[8]
A. V. Korolkova, E. G. Eferina, E. B. Laneev, I. A. Gudkova , L. A. Sevastianov, D. S. Kulyabov, Stochastization Of One-Step Processes In The Occupations Number Representati on, Proceedings 30th European Conference on Modelling and Simulation (2016) 698–704 doi:10.7148/2016-0698
-
[10]
P. Grassberger, M. Scheunert, Fock-Space Methods for Id entical Classical Objects, Fortschritte der Physik 28 (10) (1980) 547–578. doi:10.1002/prop.19800281004
-
[11]
U. C. T¨ auber, Field-Theory Approaches to Nonequilibri um Dynamics, in: Ageing and the Glass Transition, Vol. 716, Springer Berlin Heidelberg, Berlin, Heidelberg, 2005, pp. 295–348. arXiv:0511743, doi:10.1007/3-540-69684-9_7
Show all 42 references
-
[12]
Janssen, U
H.-K. Janssen, U. C. T¨ auber, The field theory approach to percolation processes, Annals of Physics 315 (1) (2005) 147 –192. arXiv:0409670, doi:10.1016/j.aop.2004.09.011
2005 doi
-
[13]
Mobilia, I
M. Mobilia, I. T. Georgiev, U. C. T¨ auber, Fluctuations and correlations in lattice models for predator-prey inter action, Physical Review E 73 (4) (2006) 040903. arXiv:0508043, doi:10.1103/PhysRevE.73.040903
2006 doi
-
[14]
Penrose, W
R. Penrose, W. Rindler, Spinors and Space-Time: Volume 1, Two-Spinor Calculus and Relativistic Fields, Vol. 1, Cam bridge University Press, 1987. doi:10.1017/CBO9780511564048
1987 doi
-
[15]
Waage, C
P. Waage, C. M. Gulberg, Studies concerning affinity, J. C hem. Educ. 63 (12) (1986) 1044. doi:10.1021/ed063p1044
1986 doi
-
[16]
A. N. Gorban, G. S. Yablonsky, Three Waves of Chemical Dy namics, Math. Model. Nat. Phenom. Vol. 10 (5) (2015) 1–5. doi:10.1051/mmnp/201510501
2015
-
[17]
N. G. van Kampen, Stochastic Processes in Physics and Ch emistry, North-Holland Personal Library, Elsevier Scienc e, 2011
2011
-
[18]
C. W. Gardiner, Handbook of Stochastic Methods: for Phy sics, Chemistry and the Natural Sciences, Springer Series i n Synergetics, 1985. 9
1985
-
[19]
Doi, Second quantization representation for classi cal many-particle system, Journal of Physics A: Mathematic al and General 9 (9) (1976) 1465–1477
M. Doi, Second quantization representation for classi cal many-particle system, Journal of Physics A: Mathematic al and General 9 (9) (1976) 1465–1477. doi:10.1088/0305-4470/9/9/008
1976 doi
-
[20]
Doi, Stochastic theory of diffusion-controlled reac tion, Journal of Physics A: Mathematical and General 9 (9) (1 976) 1479–1495
M. Doi, Stochastic theory of diffusion-controlled reac tion, Journal of Physics A: Mathematical and General 9 (9) (1 976) 1479–1495. doi:10.1088/0305-4470/9/9/009
-
[21]
Peliti, Path integral approach to birth-death proce sses on a lattice, Journal de Physique 46 (9) (1985) 1469–148 3
L. Peliti, Path integral approach to birth-death proce sses on a lattice, Journal de Physique 46 (9) (1985) 1469–148 3. doi:10.1051/jphys:019850046090146900
1985 doi
-
[22]
P. F. Verhulst, Notice sur la loi que la population suit d ans son accroissement, Vol. 10, 1838
-
[23]
W. Feller, On the theory of stochastic processes, with p articular reference to applications, Proceedings of the [F irst] Berkeley Symposium on Mathematical Statistics and Probabi lity (1949) 403–432. arXiv:1908.04294v1 [cond-mat.stat-mech] 11 Aug 2019 Два формализма стохастиз...
1949 arXiv
-
[24]
В данной нотации тензор как целостный объект обозначается просто индексом (например, xi), компоненты обозначаются подчёркнутым индексом (напри- мер, xi )
В работе используется нотация абстрактных индексов [11] . В данной нотации тензор как целостный объект обозначается просто индексом (например, xi), компоненты обозначаются подчёркнутым индексом (напри- мер, xi )
-
[25]
Латинские и ндексы из середины алфавита ( i,j,k) будут относиться к пространству векторов состояний системы
Будем придерживаться следующих соглашений. Латинские и ндексы из середины алфавита ( i,j,k) будут относиться к пространству векторов состояний системы. Лат инские индексы из начала алфавита ( a) бу- дут относиться к пространству винеровского процесса. Греч еские индексы (α) бу...
-
[26]
V., Korolkova A
Demidova A. V., Korolkova A. V., Kulyabov D. S., Sevastya nov L. A. The Method of Constructing Models of Peer to Peer Pro tocols // 6th International Congress on Ultra Modern Telecommunicat ions and Control Systems and Workshops (ICUMT). — IEEE Computer Society, 2015. — P. 557...
2015 arXiv
-
[27]
R., Korolkova A
Velieva T. R., Korolkova A. V., Kulyabov D. S. Designing I nstallations for Verification of the Model of Active Queue Ma nagement Discipline RED in the GNS3 // 6th International Congress on Ultra Modern Telecommunicat ions and Control Systems and Workshops (ICUMT). — IEEE Compu...
2015 arXiv
-
[28]
Die Grundlagen der Volterraschen Theorie des K ampfes ums Dasein in wahrscheinlichkeitstheoretischer Be hand- lung // Acta Biotheoretica
Feller W. Die Grundlagen der Volterraschen Theorie des K ampfes ums Dasein in wahrscheinlichkeitstheoretischer Be hand- lung // Acta Biotheoretica. — 1939. — Bd. 5, H. 1. — S. 11–40
1939
-
[29]
Dependence of epidemic and population veloc ities on basic parameters // Mathematical Biosciences
Mollison D. Dependence of epidemic and population veloc ities on basic parameters // Mathematical Biosciences. — 19 91. — dec. — Vol. 107, no. 2. — P. 255–287
-
[30]
V., Eferina E
Korolkova A. V., Eferina E. G., Laneev E. B., Gudkova I. A. , Sevastianov L. A., Kulyabov D. S. Stochastization Of One-Step Processes In The Occupations Number Representation // Proceedings 30th European Conference on Modelling and Simu lation. — 2016. — jun. — P. 698–704
2016
-
[31]
G., Korolkova A
Hnatič M., Eferina E. G., Korolkova A. V., Kulyabov D. S., Sevastyanov L. A. Operator Approach to the Master Equation for the One-Step Process // EPJ Web of Conferences. — 2016. — V ol. 108. — P. 58–59. — 1603.02205
2016 arXiv
-
[32]
Fock-Space Methods for Ide ntical Classical Objects // Fortschritte der Physik
Grassberger P., Scheunert M. Fock-Space Methods for Ide ntical Classical Objects // Fortschritte der Physik. — 1980 . — Vol. 28, no. 10. — P. 547–578
1980
-
[33]
Täuber U. C. Field-Theory Approaches to Nonequilibrium Dyn amics // Ageing and the Glass Transition. — Berlin, Hei- delberg : Springer Berlin Heidelberg, 2005. — Vol. 716. — P. 2 95–348. — 0511743
2005
-
[34]
Janssen H.-K., Täuber U. C. The field theory approach to pe rcolation processes // Annals of Physics. — 2005. — jan. — Vol. 315, no. 1. — P. 147–192. — 0409670
2005
-
[35]
T., Täuber U
Mobilia M., Georgiev I. T., Täuber U. C. Fluctuations an d correlations in lattice models for predator-prey interac tion // Physical Review E. — 2006. — apr. — Vol. 73, no. 4. — P. 040903. — 0508043
2006
-
[36]
Спиноры и пространство-время
Пенроуз Р., Риндлер В. Спиноры и пространство-время. Два-с пинорное исчисление и релятивистские поля. /emdash.cyr Москва : Мир, 1987. /emdash.cyr Т. 1. /emdash.cyr С. 527
1987
-
[37]
Waage P., Gulberg C. M. Studies concerning affinity // J. Chem. Educ. — 1986. — Vol. 63, no. 12. — P. 1044
1986
-
[38]
N., Yablonsky G
Gorban A. N., Yablonsky G. S. Three Waves of Chemical Dyn amics // Math. Model. Nat. Phenom. Vol. — 2015. — Vol. 10, no. 5. — P. 1–5
2015
-
[39]
Ван-Кампен Н. Г. Стохастические процессы в физике и химии. /emdash.cyr М. : Высшая школа, 1990
1990
-
[40]
Гардинер К. В. Стохастические методы в естественных на уках. /emdash.cyr Мир, 1986
1986
-
[41]
Second quantization representation for classical ma ny-particle system // Journal of Physics A: Mathematical and General
Doi M. Second quantization representation for classical ma ny-particle system // Journal of Physics A: Mathematical and General. — 1976. — Vol . 9, no. 9. — P. 1465–1477
1976
-
[42]
Stochastic theory of diffusion-controlled react ion // Journal of Physics A: Mathematical and General
Doi M. Stochastic theory of diffusion-controlled react ion // Journal of Physics A: Mathematical and General. — 1976 . — Vol. 9, no. 9. — P. 1479–1495
1976
-
[43]
Path integral approach to birth-death proces ses on a lattice // Journal de Physique
Peliti L. Path integral approach to birth-death proces ses on a lattice // Journal de Physique. — 1985. — Vol. 46, no. 9 . — P. 1469–1483
1985
-
[44]
Verhulst P. F. Notice sur la loi que la population suit da ns son accroissement. — 1838. — Vol. 10. — P. 113–117
-
[45]
On the theory of stochastic processes, with pa rticular reference to applications // Proceedings of the [F irst] Berkeley Symposium on Mathematical Statistics and Probabi lity
Feller W. On the theory of stochastic processes, with pa rticular reference to applications // Proceedings of the [F irst] Berkeley Symposium on Mathematical Statistics and Probabi lity. — 1949. — P. 403–432
1949
Reviewed August 14, 2026 · model on record in the stance chip above.
Discussion (0). Continue with ORCID to comment.