{"id":"44428e95-8f1b-474d-a4cc-64ef0e3341b2","arxiv_id":"2607.18392","paper_version":2,"verdict":"ACCEPT","confidence":"MODERATE","novelty_score":7.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":0,"one_line_summary":"A completely-positive numerical method computes the optimal time-averaged (Robinet) quantum trajectory from a system–transmission-line dilation to arbitrarily high order, plus a proof that finitely many extra binned records cannot purify beyond order Δt³.","lead":"Quantum states tracked by continuous detectors must be reconstructed from digitized, time-averaged signals. This paper gives a numerical method that computes the optimal such reconstruction to very high order, keeps the state update exactly physically valid, and proves a limit on how much extra signal information can help.","discovery_kind":"new_method","skeptic_critique":{"model":"deepseek-v4-flash","headline":"Appendix C's infinite-mode rank argument does not by itself establish the Δt^4 purity obstruction: non-proportionality of the unmeasured tail Kraus operator is also needed.","rationale":"The reader's weakest_assumption focuses on the quantum-stochastic-calculus manipulations at the dt scale; that concern is largely mitigated by the independent cascaded-systems benchmarks and the exact CP construction. However, the reader also states that the Appendix C proofs 'are rigorous and directly establish the paper's most interesting theoretical claim.' That assessment is too generous. The infinite-rank one-body correlation kernel shows that the line state cannot be expressed with finitely many modes, but purity of the system after measuring finitely many modes is controlled by the span of the resulting Kraus operators. A rank-one tail operator A with an infinite-mode line state contributes, after tracing the unmeasured rest, a single effective Kraus operator proportional to A; whether that causes impurity depends on whether A is proportional to the measured-mode operator B. Appendix C provides no such proportionality argument. The central numerical claim — an exactly CP, arbitrarily high-order, external-solver-free scheme — remains well supported by the construction and by the order-10 benchmarks, so rejection is not warranted. But the highlighted purity/no-extension theorem should be accepted only conditionally on a corrected proof or an explicit generic-proportionality argument.","tokens_in":70040,"tokens_out":31700,"duration_ms":277578,"concrete_test":"Compute the order-4 CP instrument after adding the WCW first-mode quadrature, using the explicit Kraus operators (98)–(101), for a driven qubit with L=σ_- and H=(ω/2)σ_x, where [L,[L,G]]=iωL is rank-one but B(x,y) contains the identity term ⟨x|0⟩⟨y|0⟩. Check whether the conditional map has two linearly independent Kraus operators for generic (x,y). Then repeat with g'_6 replaced by a single-mode two-photon state carrying the same effective tail operator A; if the same mixedness appears, the infinite-mode rank is not the operative obstruction and Appendix C's inference is invalid.","verdict_should_be":"CONDITIONAL","load_bearing_attack":"The paper's headline theoretical byproduct — that no finite number of integrated-signal statistics can extend the WCW purity guarantee to order Δt^4 — is supported by Appendix C through an inference that does not follow. Appendix C proves that the line states g'_6 and |ψ⊥> have one-body correlation kernels with no zero eigenvalue, hence are not representable with finitely many modes. But the purity of the reconstructed system state after measuring a finite set of modes is decided by the number of linearly independent Kraus operators acting on the system after tracing the unmeasured modes — not by the internal mode entanglement of the line. Concretely, at order 2N=4 the problematic contributions appear in Eqs. (100)–(101) as P2(x) and P3(x), both proportional to A=[L,[L,G]]. If A and the measured-mode operator B(x,y)=P0(x)⟨y|0⟩+P1(x)⟨y|1⟩ were proportional for all outcomes, the order-4 conditional state would remain pure despite the infinite tail; measuring one or several tail modes would merely re-partition the common operator A. The paper never proves generic non-proportionality of B and A. Thus the infinite-rank result establishes only that the tail cannot be fully measured with finitely many modes, not that the unmeasured tail causes impurity. The conclusion 'there is no way to preserve purity' does not follow from Appendix C alone; an additional argument about the span of the Kraus operators is required. This is a fixable proof gap, not a demonstrated counterexample, but it undermines the claim that Appendix C rigorously proves the most interesting theoretical result.","agreement_with_reader":"disagree"},"referee_report":{"model":"deepseek-v4-flash","summary":"The paper develops a completely positive (CP) discretization scheme for reconstructing time-averaged diffusive quantum trajectories (Robinet states). The construction dilates the stochastic master equation into a unitary system-plus-transmission-line dynamics, isolates a “zero mode” of the line whose homodyne measurement yields the Robinet state, expands the Dyson series, Gram–Schmidt orthogonalizes the resulting line states, and obtains explicit Kraus operators P_μ(x) such that the averaged map is CPTP to any order up to Δt^{10}. The scheme is benchmarked on a random 7-level system, with convergence verified against a cascaded-quantum-systems reference solver. The paper also claims a no-go result: the purity-preserving extension of Wonglakhon–Chantasri–Wiseman cannot reach order Δt^4 with finitely many integrated-signal statistics, because the required line state g'_6 is kinked and has an infinite one-body mode rank; an analogous claim is made for two non-commuting jump operators.","tokens_in":70323,"tokens_out":11566,"duration_ms":109064,"significance":"If the central construction is correct, the paper provides a practical, high-order, CPTP numerical integrator for the Robinet state that does not require an external solver, together with a transparent zero-mode picture of why purity is lost under time binning. The paper is explicit and reproducible: the symbolic derivations are automated, the code is archived, and the numerical verification is careful, using a random 7-level system and geometric-mean errors over 100 trajectories. The physical insight that the WCW purity result is tied to the zero-mode/first-mode truncation is appealing and likely correct in broad terms. However, one of the advertised theoretical byproducts—the no-go against Δt^4 purity with finitely many integrated statistics—is supported by an incomplete argument, as detailed below. This does not undermine the numerical scheme, but it does affect the strength of a headline claim and therefore needs correction.","major_comments":[{"comment":"The claim that no finite number of integrated-signal statistics can preserve purity at order Δt^4 is not established by the proof given. Appendix C 1 shows that the line state g'_6 (Eq. (61)) has a one-body correlation kernel with no zero eigenvalue (Eqs. (C4)–(C8)), hence that the line cannot be represented with finitely many modes. But conditional purity of the system after a finite-mode measurement is not controlled by the internal mode rank of the line state; it is controlled by the linear dependence of the residual Kraus operators on the system. At order 4 the problematic contributions are P2(x) and P3(x) in Eqs. (100)–(101), both proportional to A=[L,[L,G]]. If, for a given finite set of measured modes, the operator produced by the measured part were proportional to A for every outcome, the conditional state would remain pure even though the tail has infinitely many modes. Appendix","section":"I.C and Appendix C 1"},{"comment":"The same logical gap appears in the two-jump-operator case. The paper claims that when L1 and L2 do not commute, purity is lost at order 2 and no additional statistics can restore it. The supporting Appendix C 2 shows that the antisymmetric state |ψ⊥> (Eq. (147)) has a one-body correlation kernel with no zero eigenvalue (Eqs. (C11)–(C13)), so it populates infinitely many modes. As in the single-jump case, this only establishes that the line state is genuinely multimode. It does not exclude the possibility that all unmeasured tail Kraus operators are proportional to the operators obtained by measuring the zero modes plus any finite additional set. In that case the conditional state would remain pure despite the infinite tail. The distinction matters because the tail operator is M⊥ (Eq. (160)), and the paper does not prove that M⊥ is linearly independent of the operators generated by any f","section":"VII.D and Appendix C 2"}],"minor_comments":[{"comment":"The BCH manipulations at the dt scale with [a,a†]=1/dt are formal. A brief remark connecting these steps to the standard Ito/quantum-stochastic-calculus justification would improve rigor and accessibility. The numerical benchmarks mitigate the practical risk, but the formal presentation is faster than is ideal.","section":"Sec. II A, Eqs. (28)–(31)"},{"comment":"The geometric-mean error used in Fig. 2b is defined in the text but not in the captions. Adding a one-sentence definition in the caption would be helpful, especially since the exact error convention affects the apparent scatter.","section":"Figs. 1 and 2"},{"comment":"The printed expansion is extremely long. It may be preferable to move full symbolic expressions to the archived code or a supplementary file and retain only representative terms or recurrence relations in the appendix, to improve readability.","section":"Appendix F"},{"comment":"The Cholesky-based S^{-1/2} normalization is described briefly. A sentence clarifying the domain of validity (S strictly positive) and the non-uniqueness of the polar decomposition in degenerate cases would be useful, although these issues do not arise in the reported numerical experiments.","section":"Sec. VI"}],"recommendation":"major_revision","confidential_remarks":"The central CP reconstruction scheme appears sound and well-verified, and the numerical evidence is convincing. My main reservation is the proof of the purity no-go claim: the infinite-mode-rank argument in Appendix C does not, by itself, imply impurity under finite-mode measurements. This is a fixable proof gap rather than a demonstrated counterexample, but it affects a headline theoretical byproduct of the paper. I recommend major revision, with the expectation that the authors can supply the missing non-proportionality argument or appropriately weaken the claim."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"Worth a careful read. The main result holds up: a genuinely new dilation-based construction that gives an exactly CP, exactly-CPTP-average, arbitrary-order scheme for the Robinet state, with no external solver, and it is tested to order 10 against an independent cascaded-system solver. The zero-mode idea is the real innovation: splitting the transmission line so that measuring one quadrature yields K_I in Kraus form, then Gram-Schmidt orthogonalizing the rest. The authors are honest about the physics-level rigor in the QSC steps, and the numerics cover that risk. The paper also ships code.\n\nThe soft spot is the purity obstruction in Sec. I C and Appendix C. The appendix proves that the two-photon state g'_6 has infinite one-body mode rank, hence cannot be fully measured with finitely many modes. But that alone does not imply loss of purity: if the unmeasured tail modes all carry system operators proportional to the measured ones, the conditional state can still be pure. The paper's own escape clause—the obstruction vanishes when [L,[L,G]]=0—shows they were alert to operator structure, but they never prove generic non-proportionality. The stress-test scenario (P2 and P3 both proportional to A) is a concrete instance where infinite mode rank does not hurt purity. So the claim that no finite set of integrated records can extend WCW to order Δt^4 is not rigorously established by the argument given. I expect the claim is true for generic L,H, but it needs a linear-independence argument, not just a spectral argument on the line state.\n\nMinor issues: the floating-point saturation in Fig. 1b is asserted without evidence, and the extensions (heterodyne, finite η, time-dependent) are derived but not benchmarked.\n\nThis is a serious paper. The CP scheme itself is a structural upgrade over previous expansions and should be published. I would send it to peer review with a request to fix the Appendix C gap or weaken the claim.","headline":"Zero-mode dilation yields a solid, well-tested CP discretization of the Robinet state; the companion purity theorem has a real but fixable proof gap.","tokens_in":70945,"tokens_out":5352,"would_cite":true,"duration_ms":51168,"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":"This paper shows that the Robinet state—the best estimate of a quantum system given only time-binned measurement records—can be reconstructed using an explicitly completely positive Kraus expansion that is accurate to arbitrarily high order","keywords":["Robinet state","continuous quantum measurement","Kraus operators","completely positive maps","stochastic master equation","transmission line dilation","high-order discretization","time-averaged trajectories"],"falsifier":"Independently run the explicit order-10 Kraus sum on a high-dimensional random system, checking (i) the Choi matrix of K_I is positive semidefinite for many random x and (ii) the geometric-mean error against a small-time-step SME reference scales as Δt^{10}; a single negative Choi eigenvalue or a fitted exponent below 10 would disprove the scheme's structural claims. Separately, to test the purity barrier, include integrals of the signal against all polynomial functionals up to degree, say, 6 and check whether the purity of an initially pure state decays as Δt^4 for a generic system; if it doe","tokens_in":69784,"feed_emoji":"⚛️","tokens_out":6208,"duration_ms":56790,"temperature":0.7,"pith_summary":"This paper aims to make the Robinet state—the optimal quantum-state estimate conditioned only on time-binned, digitized measurement records—numerically computable to arbitrarily high accuracy in the bin size Δt while preserving the mathematically essential property of complete positivity. The authors achieve this by dilating the stochastic master equation into a unitary system-plus-transmission-line model, expanding the joint wavefunction in a Dyson series, Gram–Schmidt-orthogonalizing the resulting line states, and isolating the 'zero mode' whose quadrature is the binned signal. Measuring that mode yields the Robinet map as an explicit finite Kraus sum, accurate to order 2N+1; the average map is made exactly trace-preserving by an S^{-1/2} normalization. They verify the predicted scaling numerically up to order 10. The same machinery gives a sharp structural result: past the order at which a kinked, mode-entangled line state appears, no finite number of integrals of the signal against additional functions can keep the reconstructed state pure.","feed_headline":"Time-binned quantum trajectories rebuilt to order 10","feed_subtitle":"Kraus-operator scheme keeps the reconstruction completely positive and drops external solvers at any order in Δt.","key_machinery":"The central mechanism is the zero-mode decomposition H_line = H_0 ⊗ H_rest of the transmission line, where H_0 is spanned by the constant temporal mode g^(1)_1(t) = 1/√Δt. Measuring a quadrature of this mode is exactly the time-averaged homodyne signal that defines the Robinet state. The paper's construction expands the unitary evolution in a Dyson series, uses BCH normal-ordering and vacuum contractions to control the singular commutation relations [â(t),â†(u)] = δ(t−u), orthogonalizes the multi-photon line states by Gram–Schmidt, splits them into zero mode plus 'rest', and finally forms Kraus operators P_μ(x) by projecting onto the quadrature states ⟨x|k⟩. The workhorse identity is the til","core_discovery":"At the core of the paper is the construction, for any half-integer N, of the map K_I(ρ)=Σ_μ P_μ(x)ρP_μ(x)†+O(Δt^{2N+1}), with P_μ(x)=Σ_k ⟨x|k⟩P_{k,μ} obtained by projecting a Dyson expansion of the system–line state onto a Gram–Schmidt-orthogonalized basis of the line. Because the sum is a genuine Kraus form, complete positivity holds by construction even after truncation, and the normalization S^{-1/2} restores an exactly trace-preserving average. This makes the Robinet state directly computable without an external solver, and the paper validates convergence on a random 7-level system, observing the expected orders through 2N=10. In the same framework, the purity-preservation property previ","pith_inferences":["An adaptive integrator could be built directly from the paper's data: the joint system–zero-mode density at successive orders is computable before sampling the record, so the difference across orders gives a bias-free error estimate and suggests automatic step-size control.","The same Gram–Schmidt/zero-mode machinery likely transfers to discretizing continuous matrix product states: choosing non-constant bin shapes (wavelets, finite elements) changes which line states are 'zero modes' and might bypass the infinite-mode obstruction at the price of a different measurement interpretation.","The kink-state argument suggests a general diagnostic: for any truncation of the line to finitely many modes, the accuracy floor is set by the lowest-order Kraus operator whose line state has an infinite-rank correlation kernel; systems where that operator vanishes are permissive, while generic systems are not."],"forward_implications":["A single drop-in replacement for first-order SME discretization: reconstructing experimental time-binned records at low order is already principled, and stepping up to order 10 matches the exact quadrature/cascaded reference without an external solver.","The probability distribution of binned signals can be sampled directly at any order (as a Gaussian times a polynomial), enabling physical, law-accurate sampling of Robinet trajectories and pure-state unravelings with a discrete index plus continuous record.","Purity is lost at order 4 for a single monitored channel whenever [L,[L,G]] does not vanish, and no finite set of extra signal integrals can beat order 3; the obstruction is an infinite-rank kink state, not a matter of cleverly chosen functionals.","For two non-commuting jump operators, purity is lost already at order 2, and no additional measurement statistics can restore it; the same formalism identifies exactly when commuting pairs such as heterodyne detection evade this loss."],"fun_headline_variants":["Robinet state computed to 10th order without solvers","Completely positive Robinet reconstruction from binned data","Kraus method rebuilds quantum trajectories to order 10","High-order positive quantum reconstruction no external solver"],"cache_read_input_tokens":2304,"weakest_assumption_plain":"The whole coefficient hierarchy rests on the standard physics-level treatment of singular field commutators [â(t),â†(u)]=δ(t−u), the BCH normal-ordering at the dt scale with [â,â†]=1/dt, and the vacuum contraction that produces the tilted Liouvillian; if those manipulations hide a systematic error at some order, all Kraus coefficients inherit it while the order-10 numerics might still look convergent.","fun_headline_variants_meta":{"raw":{"variants":["Robinet state computed to 10th order without solvers","Completely positive Robinet reconstruction from binned data","Kraus method rebuilds quantum trajectories to order 10","High-order positive quantum reconstruction no external solver"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000944,"raw_usage":{"total_tokens":3889,"prompt_tokens":784,"completion_tokens":3105,"prompt_tokens_details":{"cached_tokens":256},"prompt_cache_hit_tokens":256,"prompt_cache_miss_tokens":528,"completion_tokens_details":{"reasoning_tokens":3040}},"tokens_in":528,"tokens_out":3105,"duration_ms":20579,"temperature":1.0,"reasoning_tokens":3040,"cache_read_input_tokens":256,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-01T15:33:54.245097+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Independently run the explicit order-10 Kraus sum on a high-dimensional random system, checking (i) the Choi matrix of K_I is positive semidefinite for many random x and (ii) the geometric-mean error against a small-time-step SME reference scales as Δt^{10}; a single negative Choi eigenvalue or a fitted exponent below 10 would disprove the scheme's structural claims. Separately, to test the purity barrier, include integrals of the signal against all polynomial functionals up to degree, say, 6 and check whether the purity of an initially pure state decays as Δt^4 for a generic system; if it doe","supporting_citations":[],"review_version":1}