{"id":"45e87c37-a066-4899-a89a-cc4ac6593bcd","arxiv_id":"1908.01225","paper_version":2,"verdict":"REJECT","confidence":"MODERATE","novelty_score":6.0,"correctness_risk":"high","formal_verification":"none","parameter_count":0,"one_line_summary":"The authors derive a general product formula for the product of m multiple integrals with respect to the compensated Poisson random measure of a pure jump Levy process.","lead":"This paper proves a formula for the product of several multiple stochastic integrals driven by a pure jump process, expressing it as a sum of other such integrals. The result gives a compact tool for working with nonlinear functionals of jump processes, relevant to chaos expansions and Malliavin calculus.","discovery_kind":"extension","skeptic_critique":{"model":"deepseek-v4-flash","headline":"Theorem 2.2's stated L2 domain fails: for m=2, n=m=1, the formula requires I1(f^2) for f in L2, but f^2 need not be L2, so RHS terms are undefined.","rationale":"The reader's weakest_assumption (contraction maps not L2-preserving) matches the main issue. I agree. The concern is load-bearing because it strikes the theorem's stated domain: not an edge case of the proof but a counterexample at m=2. The formal derivation with exponential vectors only works for kernels like p^{⊗q} where p is sufficiently integrable (or where all products occurring are L2), and polarization plus 'routine limiting' cannot produce a finite sum of multiple integrals if the product is not square integrable. This is not a disagreement with consensus; the known m=2 formula for Poisson random measures in the literature imposes integrability conditions on the contraction kernels. The result might be recoverable by restricting the domain, e.g., assuming f_k are bounded or that all relevant contractions lie in L2, which would make the paper a conditional contribution. As written, REJECT is appropriate.","tokens_in":8683,"tokens_out":10787,"duration_ms":114836,"concrete_test":"Run the following analytic check. Let η be a Poisson process with jump size 1, T=1, and set f(t)=(1-t)^(-1/4). Verify f∈L2([0,1]) and f^2∉L2([0,1]). Use the m=2 formula (2.13) for I1(f)^2: it contains the term I1(f^2), which is not defined because its kernel is not square-integrable. Alternatively compute E[I1(f)^4]; the fourth cumulant is ∫_0^1 f(t)^4 dt = ∞, so the LHS is not an L2 random variable, but every multiple integral on the RHS is L2. Either route shows Theorem 2.2/2.3 fails on its stated domain.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The central claim (Theorem 2.2) asserts the product formula for all f_k in symmetric L2 spaces. The proof reaches it by comparing coefficients of a formal power series in u's and then says the identity extends by a routine limiting argument. The load-bearing defect is that the extension target is not merely unproved; it is false as stated. The contraction operators in (2.12) do not preserve L2. For instance, take m=2, q1=q2=1, f=g. The formula (Theorem 2.3, (2.13)) gives I1(f)^2 = I2(f⊗f) + I1(f^2) + ∫ f^2 dν dt. The term I1(f^2) is a well-defined multiple integral only if f^2 ∈ L2(T). But f ∈ L2(T) does not imply f^2 ∈ L2(T). For a unit Poisson process on [0,1] with ν=δ_1 and f(t)=(1-t)^(-1/4), one has f∈L2 while f^2∉L2, so the RHS contains an undefined term; moreover E[I1(f)^4] = ∞ because the fourth cumulant ∫ f^4 dt diverges, so the LHS is not in L2. Thus the identity cannot hold in the stated domain. The same issue occurs for the general V_j^ν terms with |j|≥3 in Theorem 2.2: products of several L2 functions with a shared variable need not be L2. Consequently the 'routine limiting argument' is impossible; any correct statement needs extra integrability assumptions on the kernels (or on the Lévy measure) guaranteeing all contracted kernels lie in L2.","agreement_with_reader":"agree"},"referee_report":{"model":"deepseek-v4-flash","summary":"The paper proposes a general product formula for the product of m multiple stochastic integrals with respect to the compensated Poisson random measure of a pure-jump Lévy process. The formula expresses the product as a finite sum of multiple integrals whose kernels are obtained by a family of contraction and multiplication operations. The proof is based on exponential vectors: the authors expand an exponential functional in two ways, compare coefficients, and then extend by polarization and a claimed routine limiting argument. The m=2 case is stated as Theorem 2.3 and is shown to reduce to the classical Brownian product formula. The central claim is Theorem 2.2, equation (2.12).","tokens_in":9049,"tokens_out":5534,"duration_ms":52821,"significance":"The product formula for multiple integrals is a standard tool in stochastic analysis; a compact general formula for m factors in the Lévy setting would be useful. The exponential-vector method is elegant, the coefficient comparison is explicit, and the paper's formal derivation is a strength. If the formula were valid on the stated L2 domain, it would unify and extend known results (e.g., Shigekawa and Lee–Shih). However, the stated domain is not correct: as detailed below, the right-hand side can contain terms that are not well-defined multiple integrals for admissible L2 kernels, and the proof's final density argument is invalid. Thus the main result as stated cannot be accepted; the useful part of the computation would require a corrected statement with additional integrability assumptions.","major_comments":[{"comment":"The formula is not well-defined for all f in L2. In the case m=2, q1=q2=1, Eq. (2.13) reads I_1(f)^2 = I_2(f⊗f) + I_1(f^2) + ∫_T f^2 dλdν. For f∈L2(T,dλ×dν) there is no guarantee that f^2∈L2(T,dλ×dν), so the term I_1(f^2) is undefined. For example, take T=[0,1], ν=δ_1, and f(t)=(1−t)^{−1/4}. Then f∈L2, but ∫ f^4 dλ = ∞, so f^2∉L2; consequently the right-hand side is not defined, and in fact E[I_1(f)^4]=∞, so the left-hand side does not even belong to L2. Hence Theorem 2.2 cannot hold for the stated domain.","section":"§2, Theorem 2.3 (Eq. (2.13))"},{"comment":"The proof only establishes the identity for kernels that are finite linear combinations of symmetric tensor products p_1⊗...⊗p_q with p_i∈L2. The extension to all f_k∈L2 requires the contraction maps in (2.9)–(2.10) to be continuous on L2 and to preserve square-integrability. They do not: for a multi-index j with |j|≥3, V_j^ν(f_1,...,f_m) contains a product of several L2 functions evaluated at the same variables (t_1,z_1),...,(t_ν,z_ν), and such a product need not be integrable; even in the two-factor case the map f↦f^2 shows the problem. Therefore the 'routine limiting argument' is not a technical gap but an impossibility. A correct statement would need to impose extra integrability conditions on the kernels or on the Lévy measure and then prove the identity for that restricted domain; as written, the argument cannot be repaired by a denseness argument.","section":"§3, last paragraph ('routine limiting argument')"}],"minor_comments":[{"comment":"There is a typographical error in the second exponential factor: an unmatched closing parenthesis appears after 'm−1'. The intended expression is exp{∫_T ( e^{Σρ_k} − Σ e^{ρ_k} + m − 1 ) ν(dz)ds }.","section":"§3, Eq. (3.8)"},{"comment":"The notation In(fn)Im(gm) is confusing; in Theorem 2.3 the symbols f and g denote kernels of order n and m, not functions named f_n and g_m. Also, the kernel definition in (2.14) loses some arguments in the displayed formula; it should show the dependence on (s_1,z_1),...,(s_{n+m−k−2l}, z_{n+m−k−2l}) explicitly.","section":"§2, Theorem 2.3 (2.13)–(2.14)"},{"comment":"The symbol T is used both for the time horizon and for the domain [0,T]×R0. This overload is confusing; a different letter for the product domain would improve readability.","section":"§2, after Eq. (2.2)"},{"comment":"There is a typo in the exponent of the last term: 'u^{lliκm}_{iκm}' should read 'u^{liκm}_{iκm}'.","section":"§3, Eq. (3.12)"}],"recommendation":"reject","confidential_remarks":"The counterexample to Theorem 2.2 is elementary and robust: it uses a unit Poisson process and a square-integrable integrand whose square is not square-integrable. The formal exponential-vector computation may be correct for sufficiently integrable kernels, and a revised statement with such assumptions might be worth considering, but the present theorem as stated is false. I therefore recommend rejection."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"Dear colleague,\n\nYou should know that this paper has a nice proof idea and a genuinely new-looking formula for products of m≥3 multiple integrals with respect to the compensated Poisson random measure. But Theorem 2.2 as stated is false, not merely unproved: the RHS can contain terms that are not defined for arbitrary square-integrable kernels. The stress-test example is right. For m=2, n=m=1, the formula gives I(f)^2 = I_2(f⊗f) + I_1(f^2) + ∫ f^2 ν dt. If ν=δ_1, T=1, f(t)=(1-t)^{-1/4}, then f∈L^2 but f^2∉L^2, so I_1(f^2) is undefined (and the fourth moment diverges). So the 'routine limiting argument' at the end of Section 3 cannot work.\n\nWhat is genuinely good: the exponential-vector/chaos-expansion derivation is short and elegant. It reduces correctly to Shigekawa's Brownian formula and to the m=2 result of Lee-Shih when the kernels are nice. For m≥3, the general formula appears new, and the bookkeeping with multi-indices is clear enough to follow. The proof works verbatim for kernels that are finite linear combinations of tensor products of bounded functions, say, or whenever all contracted kernels remain in L^2.\n\nThe soft spots are exactly where the reader's report points. The contraction operators V^ν_j and ⊗̂^μ_i do not preserve L^2 when |j|≥2 or |i|≥2, and no attempt is made to check integrability of the output. Also, the exponential functional (3.1) requires something like e^ρ-1 ∈ L^2; for ρ=log(1+up) with p∈L^2, B contains ∫ p_i p_j ν dt which may be infinite unless p_i p_j ∈ L^1. That is the same root issue.\n\nThis is a load-bearing flaw, but a fixable one. The paper should state the theorem for kernels satisfying explicit sufficient conditions (e.g., all needed contractions in L^2, or work with bounded functions and compact support and then close under a norm that controls the contractions). The m=2 case with the extra condition fg∈L^2 is known and fine; the general m version, with such conditions, would be a useful contribution.\n\nWho should read it: researchers in stochastic analysis who want product formulas for Lévy multiple integrals. It deserves a serious referee rather than a desk reject, because the idea is sound and the gap is a standard domain issue, not an error of principle. But my own verdict: as written, the theorem is not true on its stated domain.","headline":"Elegant formal computation of a general product formula for Lévy multiple integrals, but the stated L2 domain is indefensible; the theorem needs extra integrability conditions before it can be used.","tokens_in":9540,"tokens_out":5412,"would_cite":false,"duration_ms":53464,"reading_group":"maybe","serious_thinker":"yes","would_accept_peer_review":true},"rs_alignment":null,"lean_confirmation":null,"pith_extraction":{"msc":["60H05","60G51","60H30"],"pacs":[],"model":"deepseek-v4-flash","headline":"The product of any finite number of multiple stochastic integrals of a pure-jump Lévy process is a finite sum of single multiple integrals of symmetrically contracted kernels.","keywords":["Lévy process","nonlinear functional of Lévy process","multiple integrals","chaos expansion","product formula","exponential vector","polarization technique"],"falsifier":"Take a Lévy process with Lévy measure $\\nu(dz)=\\mathbf{1}_{(0,1)}(z)\\,dz$, set $m=3$, and let $f_1=f_2=f_3=h$ with $h(t,z)=z^{-1/4}$ on $(0,1)^2$. Then $h\\in L^2(dt\\,d\\nu)$ but $h^3\\notin L^2(dt\\,d\\nu)$. The right-hand side of (2.12) contains the term $I_1(h^3)$ (the triple contraction with one shared variable), which is not defined, while the left-hand side $I_1(h)^3$ is a well-defined random variable. Hence the theorem is false as stated for arbitrary $L^2$ kernels; checking just this term would settle the claimed domain.","tokens_in":8465,"feed_emoji":"🧮","tokens_out":21916,"duration_ms":198443,"temperature":0.7,"pith_summary":"The paper gives a single closed-form formula for the product of any finite number of multiple integrals with respect to the compensated Poisson random measure of a purely discontinuous Lévy process. The formula expands the product into a finite sum of multiple integrals whose kernels are obtained by symmetrized contractions of the original kernels, with explicit factorial coefficients. This generalizes the classical Brownian product formula and unifies earlier two-integral results. The proof is short: it compares the chaos expansions of a product of exponential vectors in two ways and reads off coefficients.","feed_headline":"All Lévy multiple integrals multiply via one contraction sum","feed_subtitle":"A compact contraction identity covers all pure-jump Lévy processes with a short exponential-vector proof.","key_machinery":"The exponential vector $\\mathcal{E}(\\rho)=\\exp\\{\\int_{\\mathbb{T}}\\rho\\,d\\tilde N - \\int_{\\mathbb{T}}(e^\\rho-1-\\rho)\\,d\\nu dt\\}$, whose chaos expansion has kernels $(e^\\rho-1)^{\\hat\\otimes n}$. The proof expands the product of $m$ such exponential vectors in two ways—directly as a product of independent chaos expansions, and as a single exponential vector with an extra deterministic correction factor—and then polarizes by taking $\\rho_k=\\log(1+u_k p_k)$. Comparing coefficients of $u_1^{q_1}\\cdots u_m^{q_m}$ yields the contraction formula. The combinatorial core is the pair of operators $\\hat\\otimes_i^\\mu$ and $V_j^\\nu$ that glue or integrate shared time–jump variables.","core_discovery":"Theorem 2.2 asserts that for $m\\ge2$ and symmetric kernels $f_k \\in (L^2([0,T]\\times \\mathbb{R}_0, dt\\otimes \\nu))^{\\hat\\otimes q_k}$, the product $\\prod_{k=1}^m I_{q_k}(f_k)$ equals a finite sum, indexed by multi-indices $\\vec l,\\vec n$, of single multiple integrals $I_{|q|+|\\vec n|-|\\chi(\\vec l,\\vec n)|}(\\hat\\otimes_{\\vec i}^{\\vec l}\\hat\\otimes V_{\\vec j}^{\\vec n}(f_1,\\dots,f_m))$, with weights $\\prod_k q_k! /(\\prod_\\alpha l_{i_\\alpha}! \\prod_\\beta \\mu_{j_\\beta}! \\prod_k (q_k-\\chi(k,\\vec l,\\vec n))!)$. In words, the product of several multiple jump integrals is a linear combination of one multiple jump integral of a symmetrized kernel in which groups of the original variables are glued together or integrated out. For $m=2$ the formula reduces to a two-integral expansion; in the Brownian limiting case it reduces to the classical contraction formula.","pith_inferences":["The stated domain over all $L^2$ kernels is probably too broad: the contraction operator $V_j^\\nu$ for $|j|\\ge3$ does not preserve $L^2$, so the theorem likely needs an added integrability condition on the kernels or on the Lévy measure.","The exponential-vector derivation suggests an analogous product formula for multiple integrals with respect to other random measures that admit a Wick-type exponential and a chaos decomposition, as long as the deterministic correction factor is computed correctly.","A concrete check of the formula for $m=3$ and a Poisson process with unit jumps, comparing $I_1(1)^3$ against the direct polynomial $(N(T)-T)^3$, would pin down the contraction coefficients and expose any missing terms."],"forward_implications":["For any $m\\ge2$, the product of $m$ multiple integrals with respect to the jump measure is a finite sum of multiple integrals; no infinite series appears beyond the finite indexing by subsets.","The $m=2$ case recovers the known formula for products of two Lévy integrals and, in the Brownian limit, the classical contraction rule.","Because square-integrable functionals of the Lévy process have a unique chaos expansion, the formula allows products of arbitrary chaos expansions to be re-expanded into a single chaos series.","The compact form of the coefficients makes the formula suitable for moment computations and limit theorems for statistics built from jump processes."],"supporting_citations":[{"why":"Supplies the classical Brownian product formula that the new formula generalizes and must reduce to for $m=2$.","marker":"[9]"},{"why":"Gives the earlier two-integral product formula for Lévy integrals that this paper extends to $m\\ge2$ with respect to the jump measure.","marker":"[3]"},{"why":"Provides the polarization technique and the chaos-expansion framework used throughout the proof.","marker":"[2]"},{"why":"Used to justify the stochastic differential rule for the exponential vector, which yields its chaos expansion.","marker":"[7]"}],"fun_headline_variants":["Lévy integrals multiply: one explicit formula","Lévy integrals: one sum contracts all products","Product of Lévy integrals: a single contraction sum","All Lévy multiple integrals reduce to one finite expansion"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The proof shows the identity for kernels that are finite sums of tensor products of one-variable functions, then asserts it extends to all square-integrable kernels by a routine limiting argument; that extension requires the contracted kernels on the right-hand side to remain square-integrable, which need not happen when a term glues together three or more kernels.","fun_headline_variants_meta":{"raw":{"variants":["Lévy integrals multiply: one explicit formula","Lévy integrals: one sum contracts all products","Product of Lévy integrals: a single contraction sum","All Lévy multiple integrals reduce to one finite expansion"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.001003,"raw_usage":{"total_tokens":4166,"prompt_tokens":793,"completion_tokens":3373,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":409,"completion_tokens_details":{"reasoning_tokens":3311}},"tokens_in":409,"tokens_out":3373,"duration_ms":22046,"temperature":1.0,"reasoning_tokens":3311,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-14T15:20:57.453059+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Take a Lévy process with Lévy measure $\\nu(dz)=\\mathbf{1}_{(0,1)}(z)\\,dz$, set $m=3$, and let $f_1=f_2=f_3=h$ with $h(t,z)=z^{-1/4}$ on $(0,1)^2$. Then $h\\in L^2(dt\\,d\\nu)$ but $h^3\\notin L^2(dt\\,d\\nu)$. The right-hand side of (2.12) contains the term $I_1(h^3)$ (the triple contraction with one shared variable), which is not defined, while the left-hand side $I_1(h)^3$ is a well-defined random variable. Hence the theorem is false as stated for arbitrary $L^2$ kernels; checking just this term would settle the claimed domain.","supporting_citations":[{"cited_title":"L\\'evy processes and infinitely divisible distributions","cited_arxiv_id":null,"evidence_quote":"Supplies the classical Brownian product formula that the new formula generalizes and must reduce to for $m=2$."},{"cited_title":"and ksendal, B","cited_arxiv_id":null,"evidence_quote":"Provides the polarization technique and the chaos-expansion framework used throughout the proof."},{"cited_title":"The Malliavin calculus and related topics","cited_arxiv_id":null,"evidence_quote":"Used to justify the stochastic differential rule for the exponential vector, which yields its chaos expansion."}],"review_version":1}