{"id":"5acf7ed5-80e5-445a-a255-63589fdba5c5","arxiv_id":"2608.09903","paper_version":1,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":6.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":2,"one_line_summary":"The Fokker-Planck density of a nonlinear SDE is the diagonal of a Lindblad-evolved quantum state, so classical expectations and two-time correlations become exact quantum channel observables.","lead":"This paper shows that the probability law of a broad class of nonlinear stochastic differential equations can be represented exactly as the position diagonal of a quantum density matrix evolving under a Lindblad master equation, with classical diffusion acting as decoherence. The result gives a mathematically controlled bridge from classical SDE statistics to open-quantum-system simulation, and the numerics validate it on double-well Langevin and noisy Lorenz-63 models.","discovery_kind":"unification","skeptic_critique":{"model":"deepseek-v4-flash","headline":"The theorem's standing flow/domain assumptions are not verified for the paper's own non-globally-Lipschitz examples; this gap is load-bearing but testable, so the conditional verdict stands.","rationale":"The reader's weakest_assumption correctly identifies the gap between the standing assumptions and the non-globally-Lipschitz examples as the key load-bearing point. I agree that the continuum derivation is mathematically standard given those assumptions, and the numerical evidence is consistent with the claimed identity. The concern is not that the result is false, but that the theorem's hypotheses are not verified for the paper's own test cases, so the exactness claim is broader than what is proven. A concrete analytical check on the double-well and Lorenz examples would settle whether the gap is merely technical or whether additional conditions (e.g., Lyapunov-type growth or refined integrability) are needed. Since this is exactly the reader's stated condition for the CONDITIONAL verdict, my stress test does not move the verdict; it reinforces it. The other issues (Galerkin error theorem's unverified H^s regularity, code not yet shipped) are secondary to the central claim; the flow/domain assumption gap is the single most load-bearing concern.","tokens_in":35083,"tokens_out":8482,"duration_ms":81419,"concrete_test":"For the double-well Langevin SDE (64) with Gaussian initial condition, verify the standing assumptions analytically: (a) prove global existence of the C^1 stochastic flow via the dissipative Lyapunov function E(q,v)=v^2/2+(q^2-1)^2/4, which controls no-blow-up for all initial conditions; and (b) prove the moment bound E∫_0^T (||H_0 ψ_t||^2 + ||H_j ψ_t||^2) dt < ∞ using the known exponential/super-exponential tail decay of the Fokker–Planck solution for the quartic potential. If (b) fails, the trace-class Ito formula in Theorem 2 is not justified for this example; if both hold, repeat the same check for the Lorenz–63 SDE (82).","verdict_should_be":"UNCHANGED","load_bearing_attack":"Theorem 2 is conditional on (i) a global C^1 stochastic flow of diffeomorphisms for each Brownian realization and (ii) unbounded-operator domain conditions justifying the trace-class Ito formula used in Appendix B. The only stated sufficient condition (Section II) is global Lipschitz vector fields with linear growth and bounded derivatives. The paper's own experiments use a cubic double-well force in Eq. (64b) and a quadratic Lorenz drift in Eq. (82), neither globally Lipschitz. No proof is supplied that these examples satisfy the standing assumptions: neither the global two-sided flow (no finite-time explosion for every initial condition and almost every Brownian path) nor the integrability condition E∫_0^T (||H_0 ψ_t||^2 + Σ_j ||H_j ψ_t||^2) dt < ∞ required for the Hilbert-space Ito formula on the common test domain. The numerics strongly suggest the results hold, but the central claim is stated as exact; as written, the theorem does not formally cover the headline examples. Closing this gap is necessary before the KLM mapping can be claimed for the class of nonlinear SDEs advertised in the abstract.","agreement_with_reader":"agree"},"referee_report":{"model":"deepseek-v4-flash","summary":"The paper introduces a Kolmogorov–Lindblad mapping (KLM) that lifts the pathwise probability density of an Itô SDE to a half-density in L², shows that each Brownian realization generates a unitary half-density flow, and proves that averaging the resulting pure states yields a Lindblad equation whose position diagonal equals the Fokker–Planck density. It further establishes forward and backward intertwining between the KLM semigroup and the Kolmogorov semigroups, an exact quantum-regression identity for two-time correlations, and a structure-preserving Galerkin projection that retains Lindblad form at finite dimension. Numerical tests on double-well Langevin dynamics and noisy Lorenz–63 are compared with independent Fokker–Planck references.","tokens_in":35156,"tokens_out":6206,"duration_ms":63090,"significance":"If the stated assumptions are met, this is a clean and potentially useful exact law-level embedding: it avoids the L¹/L² normalization mismatch of direct amplitude encoding, produces a CPTP channel before any quantum-backend choice, and gives a concrete path from classical SDE statistics to Lindblad simulation. Strengths of the manuscript include a self-contained derivation with no parameters fitted to the reference solution, an explicit separation of approximation cost from coherent-access cost, and numerical comparisons against independent Fokker–Planck solvers with reported extrapolation uncertainties.","major_comments":[{"comment":"Theorems 1–3 and Corollary 1 are proved only under the standing assumptions of Section II, whose only stated sufficient condition is global Lipschitz vector fields with linear growth. The paper's headline examples violate this: Eq. (64b) contains the cubic term -q_t^3+q_t and Eq. (82) contains quadratic products x_t y_t and x_t z_t, neither of which is globally Lipschitz. The manuscript does not verify the global C¹ two-sided stochastic flow for these coefficients, nor the domain/integrability condition E∫_0^T (‖H_0ψ_t‖² + Σ_j‖H_jψ_t‖²) dt < ∞ used for the Hilbert-space Itô formula in Appendix B. This gap is load-bearing because the abstract advertises the mapping for general nonlinear stochastic dynamics and the numerical sections apply the theorems to exactly these examples. I ask for explicit verification for the tested models, for instance via Lyapunov functions and a priori weighted-Sobolev estimates, or a restatement of the hypotheses so that the theorems are not stronger than their verified domain.","section":"Section II (standing assumptions) and Eqs. (64b), (82)"},{"comment":"Appendix D ends with the sentence 'The preceding calculation is summarized by the following theorem,' but no theorem is stated. The Galerkin residual estimate in this appendix is the stated basis for the weak-bias condition (59) used in Theorem 4, so the missing statement is a missing piece of support rather than a typographical issue. Please supply the theorem with its precise hypotheses and proof, or delete the sentence and clearly label the appendix as an informal derivation sketch.","section":"Appendix D"}],"minor_comments":[{"comment":"The sentence after Eq. (24) refers to 'the generator on the right-hand side of (26)' before Eq. (26) has been displayed; the equation order or cross-reference should be adjusted.","section":"Eqs. (24)–(26)"},{"comment":"The symbol bΦ_Q is reused from the mass-matrix transformation, and its role in the quadrature-level Galerkin matrix H_{j,h}^{(Q)} should be spelled out so that the reader can distinguish the raw basis, the isometry, and the sampled coefficient matrices.","section":"Eq. (52)"},{"comment":"The notation P_t^* is used in the proof of Theorem 3 before its definition has been fully separated from the later display; a one-sentence definition immediately before the proof would improve readability.","section":"Appendix C"},{"comment":"The convergence-rate fits in Eqs. (92)–(93) use only five cutoff values each; the authors correctly call them preasymptotic, but the discussion should state explicitly that these fits are not used as a substitute for an asymptotic error bound in any of the theorem statements.","section":"Section IV.C"}],"recommendation":"major_revision","confidential_remarks":"To the editor: the central construction appears mathematically sound within its stated hypotheses, and I found no internal inconsistency in the derivation itself. The main risk is that the abstract and introduction promise a mapping for general nonlinear SDEs while the theorem hypotheses are not checked for the paper's own numerical examples. This is fixable either by additional analysis for the double-well and Lorenz models or by careful narrowing of the claims. The omitted Appendix D theorem should also be supplied."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"Short version: the core result is real. For SDEs with a global two-sided C^1 stochastic flow, averaging the half-density transport gives a Lindblad equation whose position diagonal is exactly the Fokker-Planck density. The derivation is standard Stratonovich/Ito calculus done carefully, and Theorem 3's intertwining and Corollary 1's regression identity are clean, exact statements. The numerical section is honest: independent Fokker-Planck reference, no fitting, convergence in K. The structure-preserving Galerkin projection is a sensible practical interface. All of that deserves credit.\n\nWhat is genuinely new is not the half-density lift itself (Applebaum had it) but the law-level reading: the averaged pure states form a CPTP channel, the classical semigroup is recovered as the diagonal, and two-time correlations come out as an exact regression formula. That combination, with the discretization, is a useful contribution.\n\nSoft spots, in order of weight. First, the standing assumptions are not verified for the paper's own examples. The stated sufficient condition is global Lipschitz with linear growth; the double-well force is cubic and the Lorenz drift is quadratic. The theorems are asserted for the flow-regular class, but no proof is given that these specific SDEs have global two-sided flows and satisfy the integrability condition for the trace-class Ito formula. The numerics strongly suggest it is true, and for these particular systems it is probably provable, but as written the headline examples fall outside the proven regime. This has to be fixed before the claimed scope is justified. Second, Appendix D's Galerkin error bound relies on an a priori H^s regularity assumption and a two-power loss; the bound is conditional, not a theorem about these examples. That is stated openly, so it is a limitation, not a flaw. Third, no code is shipped; only a promise.\n\nThe citation pattern looks fine: Applebaum and Kunita are properly credited, and the single author-overlapping citation is contextual. No fitted parameters, no circularity. This is a serious paper with a clear mathematical core. It deserves a proper referee. I would send it out, with the request that the authors either prove the flow/domain conditions for their examples or explicitly narrow the theorems to the verified class. If that gap is closed, I would be happy to cite it.","headline":"A correct and useful exact mapping from nonlinear SDE laws to Lindblad dynamics, with one load-bearing regularity gap that is closable but must be closed.","tokens_in":35837,"tokens_out":1845,"would_cite":true,"duration_ms":17592,"reading_group":"yes","serious_thinker":"yes","would_accept_peer_review":true},"rs_alignment":null,"lean_confirmation":null,"pith_extraction":{"msc":["60H10","60H15","35Q84","81P68"],"pacs":[],"model":"deepseek-v4-flash","headline":"The paper proves that, under flow regularity and domain assumptions, the statistics of a nonlinear SDE are exactly the position diagonal of a state evolving by a Lindblad equation.","keywords":["Kolmogorov–Lindblad mapping","stochastic differential equations","Fokker–Planck equation","Lindblad equation","half-density","stochastic Schrödinger equation","quantum channel","Galerkin discretization"],"falsifier":"Take a globally Lipschitz nonlinear diffusion, for instance $dX_t=\\sin(X_t)\\,dt+dW_t$, and prepare two trace-one encodings with the same initial position density $p_0$ but different phases; evolve both under the KLM equation for a time where an independent Fokker–Planck reference is accurate. Any difference in $\\mathrm{Tr}(M_\\varphi\\Gamma(t))$ for a bounded $\\varphi$, or any deviation of the reconstructed diagonal from the reference density, would falsify the claimed exact intertwining.","tokens_in":34709,"feed_emoji":"⚛️","tokens_out":12428,"duration_ms":108566,"temperature":0.7,"pith_summary":"The paper aims to establish an exact bridge between classical stochastic dynamics and open quantum systems: for a broad class of nonlinear stochastic differential equations, the probability law can be encoded as the position diagonal of a density operator that evolves by a Lindblad equation. The construction works by taking, for each Brownian realization, the square root of the transported probability density; this half-density satisfies a linear stochastic Schrödinger equation, and averaging the resulting pure states produces a completely positive, trace-preserving channel. On the diagonal, the channel reproduces the Fokker–Planck density exactly, so expectations of bounded functions become quantum traces without the unknown normalization required by direct amplitude encoding. Two-time correlations and, by iteration, ordered multi-time correlations obey exact quantum regression identities, and a structure-preserving Galerkin projection keeps the finite-dimensional equation in Lindblad form. If the theorems hold, the main obstruction to quantum simulation of nonlinear SDEs shifts from representing the nonlinearity to approximation and coherent operator access.","feed_headline":"Nonlinear SDE laws become exact quantum channels","feed_subtitle":"A half-density lift maps Fokker–Planck laws to Lindblad channels, removing the normalization mismatch.","key_machinery":"The central object is the pathwise half-density $\\psi_t(x)=\\psi_0((\\Phi^\\omega_t)^{-1}(x))\\,|\\det D(\\Phi^\\omega_t)^{-1}(x)|^{1/2}$, whose modulus square is the noise-conditional density. Its generator is $K_V\\psi=V\\cdot\\nabla\\psi+\\frac{1}{2}(\\nabla\\cdot V)\\psi$; the added divergence term is the square-root Jacobian correction that makes each Brownian realization act unitarily on $L^2$. With $H_j=-iK_{V_j}$, the Itô correction of the stochastic Schrödinger equation turns the Brownian quadratic variation into the double-commutator dissipator, and the diagonal map $D$ plus the multiplication-operator identity $E^\\dagger_t(M_f)=M_{P_t f}$ carry the exact intertwining with the classical semigroup.","core_discovery":"On the paper's own terms, the central discovery is Theorem 2. For an Itô SDE $dX_t=b(X_t)\\,dt+\\sigma(X_t)\\,dW_t$ rewritten in Stratonovich form with vector fields $V_0,\\dots,V_m$, define the Hermitian half-density operators $H_j=-i(V_j\\cdot\\nabla+\\frac{1}{2}\\nabla\\cdot V_j)$ on $L^2(\\mathbb{R}^d)$. If $\\psi_t=U^\\omega_t\\psi_0$ is the pathwise half-density transported by the stochastic flow $\\Phi^\\omega_t$, then $\\Gamma(t)=\\mathbb{E}_\\omega[|\\psi_t\\rangle\\langle\\psi_t|]$ obeys $$\\frac{d\\Gamma}{dt}=-i[H_0,\\Gamma]-\\frac{1}{2}\\sum_{j=1}^m[H_j,[H_j,\\Gamma]],$$ and its position diagonal is exactly the Fokker–Planck density, $D\\Gamma(t)=p(t,\\cdot)$, so $E[\\varphi(X_t)]=\\mathrm{Tr}(M_\\varphi\\Gamma(t))$. The paper further claims exact forward and backward semigroup intertwining, an exact quantum regression formula for two-time correlations, a mixed-unitary structure for the channel, and a Galerkin projection that keeps the finite-dimensional equation in Lindblad form. Classical diffusion enters as decoherence generated by the Hermitian jump operators $H_j$, and the zero-noise limit is deterministic unitary transport of the half-density.","pith_inferences":["If the continuum identities are exact, then any mismatch between the projected channel's diagonal and an accurate Fokker–Planck reference is attributable to Galerkin, quadrature, or time-stepping error; the experiments can therefore be read as a clean benchmark for approximation spaces and operators.","For SDEs whose Stratonovich diffusion fields commute, the dissipator becomes a Gaussian twirl of unitary translations and the channel may be simulable more cheaply than generic Lindblad simulation; the paper mentions Gaussian twirls but does not exploit this into an algorithm.","The backward intertwining invites a goal-oriented weak-error analysis: regularity of the solution to the backward Kolmogorov equation, rather than strong pathwise regularity of the half-density, should control observable errors. The paper sketches this direction but leaves the estimates open.","A direct hardware test of the gauge freedom would prepare two different phase or purification encodings of the same initial density and verify identical statistics for multiplication observables; the paper does not perform such a test."],"forward_implications":["Bounded classical observables of the SDE become ordinary quantum expectations $\\mathrm{Tr}(M_\\varphi\\Gamma(t))$ with no unknown time-dependent normalization, so the normalization mismatch of direct amplitude encoding disappears.","Every trace-one initial encoding with the same position density produces the same classical law and the same multiplication-observable statistics; off-diagonal coherences are a gauge freedom for the law.","Two-time correlations obey the exact regression identity $\\mathbb{E}[\\varphi_1(X_t)\\varphi_2(X_s)]=\\mathrm{Tr}[M_{\\varphi_1}E_{t-s}(M_{\\varphi_2}\\Gamma(s))]$, and recursive application yields ordered multi-time correlations.","The structure-preserving Galerkin projection yields a finite-dimensional mixed-unitary, completely positive, trace-preserving channel before any quantum simulator is applied.","In the deterministic limit $m=0$ the Lindblad equation reduces to unitary Hamiltonian evolution for the half-density, recovering the classical transport lift."],"supporting_citations":[{"why":"supplies the stochastic-flow theory behind the standing regularity and completeness assumptions.","marker":"[17]"},{"why":"constructs the half-density operator action and unitary Brownian evolution used for the pathwise lift.","marker":"[19]"},{"why":"gives the forward Kolmogorov equation that the diagonal must reproduce.","marker":"[36]"},{"why":"provides the stochastic calculus in infinite dimensions used to derive the master equation from the stochastic Schrödinger equation.","marker":"[37]"},{"why":"supplies the characterization of Markovian master-equation generators used to identify the dissipative form.","marker":"[40]"},{"why":"identifies the Lindblad form of completely positive semigroup generators used for the double-commutator dissipator.","marker":"[41]"},{"why":"supplies the ancilla-dilation construction that turns a finite Lindblad step into a CPTP map without postselection.","marker":"[51]"}],"fun_headline_variants":["SDEs to Lindblad: exact mapping via half-density lift","Half-density lift turns SDE laws into Lindblad channels","Exact quantum channels for nonlinear SDEs","From Fokker–Planck to Lindblad without approximation","Kolmogorov–Lindblad: SDEs as open quantum systems"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The load-bearing premise is that the SDE's vector fields generate a global, two-sided stochastic flow of smooth diffeomorphisms, with all initial points moved by the same Brownian realization and with enough far-field decay to make the half-density operators self-adjoint; note that the paper's own numerical examples use a cubic force and a quadratic drift that do not satisfy the stated global-Lipschitz sufficient condition.","fun_headline_variants_meta":{"raw":{"variants":["SDEs to Lindblad: exact mapping via half-density lift","Half-density lift turns SDE laws into Lindblad channels","Exact quantum channels for nonlinear SDEs","From Fokker–Planck to Lindblad without approximation","Kolmogorov–Lindblad: SDEs as open quantum systems"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000283,"raw_usage":{"total_tokens":1762,"prompt_tokens":1128,"completion_tokens":634,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":744,"completion_tokens_details":{"reasoning_tokens":562}},"tokens_in":744,"tokens_out":634,"duration_ms":5106,"temperature":1.0,"reasoning_tokens":562,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-11T04:43:36.806478+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Take a globally Lipschitz nonlinear diffusion, for instance $dX_t=\\sin(X_t)\\,dt+dW_t$, and prepare two trace-one encodings with the same initial position density $p_0$ but different phases; evolve both under the KLM equation for a time where an independent Fokker–Planck reference is accurate. Any difference in $\\mathrm{Tr}(M_\\varphi\\Gamma(t))$ for a bounded $\\varphi$, or any deviation of the reconstructed diagonal from the reference density, would falsify the claimed exact intertwining.","supporting_citations":[{"cited_title":"Tennie and L","cited_arxiv_id":null,"evidence_quote":"supplies the stochastic-flow theory behind the standing regularity and completeness assumptions."},{"cited_title":"Gnanasekaran, A","cited_arxiv_id":null,"evidence_quote":"constructs the half-density operator action and unitary Brownian evolution used for the pathwise lift."},{"cited_title":"Da Prato and J","cited_arxiv_id":null,"evidence_quote":"supplies the characterization of Markovian master-equation generators used to identify the dissipative form."},{"cited_title":"Tronci and I","cited_arxiv_id":null,"evidence_quote":"identifies the Lindblad form of completely positive semigroup generators used for the double-commutator dissipator."}],"review_version":1}