{"id":"36ff0c8c-eb34-4415-9f35-927e9af72d86","arxiv_id":"2506.23386","paper_version":1,"verdict":"REJECT","confidence":"MODERATE","novelty_score":3.0,"correctness_risk":"high","formal_verification":"none","parameter_count":0,"one_line_summary":"A hybrid qubit-field Wigner function for the Jaynes-Cummings model is constructed, reproducing standard Rabi and inversion results, but the one-mode Wigner function and the reduced-field purity formula are incorrect as printed.","lead":"The paper constructs a full phase space Wigner function for the Jaynes-Cummings atom-field model and re-derives the model's Rabi oscillation, collapse-revival, and field-purity dynamics in that language. The main closed-form Wigner function appears internally consistent, but the one-mode Wigner function (Eq. 28) and the purity formula (Eq. 40) contain errors that fail elementary consistency checks.","discovery_kind":"extension","skeptic_critique":{"model":"deepseek-v4-flash","headline":"The purity formula (40) fails the Fock check: with C_n = δ_{n,0} it returns ξ = 1, whereas the exact reduced-field purity is cos^4(gt)+sin^4(gt), equal to 1/2 at t=π/(4g). Entanglement conclusions are unsupported.","rationale":"The reader's weakest_assumption identifies Eq. (40) and Eq. (28) as fragile premises; my focused check confirms Eq. (40) is the more load-bearing one because the paper's entanglement conclusions rest on it. The Fock-state test is decisive and does not depend on any approximation or external convention: Eq. (40) contradicts the exact two-level evolution for an initial state that lies inside the model's own Hilbert space. I also agree with the reader that Eq. (23) itself passes the normalization check and reproduces the standard Rabi/inversion results, which weakens any global claim that the whole construction is unsound. Still, since the central advertised outcome—using the Wigner-function purity to measure entanglement—is quantitatively wrong as printed, the REJECT verdict is justified. The paper may be salvageable if Eq. (40) is corrected and Eq. (28) is properly derived, but neither fix is cosmetic.","tokens_in":12479,"tokens_out":3092,"duration_ms":34020,"concrete_test":"Compute ξ(t) for C_n = δ_{n,0} in two ways: (i) insert the delta coefficients into Eq. (40); (ii) form the exact reduced field from Eq. (22) and compute tr(ρ_f^2). Evaluate both at t = π/(4g). Any correct purity expression must give 1/2; Eq. (40) as printed gives 1. If the discrepancy appears, re-derive Eq. (40) from the reduced Wigner function (38) using the orthogonality relation (26) in a truncated Fock subspace (n=0,1) to isolate where the cross terms are lost.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The central advertised application is the entanglement dynamics obtained from the purity ξ(t) of the reduced field. Equation (40) is the load-bearing premise, and it is directly falsified by a two-dimensional Fock-state check. For an initially excited atom with field in |0⟩, C_n = δ_{n,0}. The exact JC evolution (20) gives |ψ(t)⟩ = cos(gt)|e,0⟩ − i sin(gt)|g,1⟩, so the reduced field is cos^2(gt)|0⟩⟨0| + sin^2(gt)|1⟩⟨1|, whose purity is cos^4(gt)+sin^4(gt), i.e., 1/2 at t=π/(4g). Inserting C_n = δ_{n,0} into Eq. (40) gives, for all t, ξ(t)=1: the double sum contains only n=m=0 and equals 1/2 + 1/2 cos(2gt(√1−√1)) = 1, while the second sum vanishes because C_{−1}=0. Thus Eq. (40) cannot be the purity of the reduced field. This also invalidates Eq. (41), since for this initial state the exact purity is periodic in gt and revives at T_rev; it does not tend to a t→∞ limit. A further internal inconsistency is Eq. (28), which omits the cross-field Laguerre terms L_{r,r+1}, L_{r+1,r} and the ground-state factor W_{r+1} that follow from specializing (23) to C_n = δ_{n,r}; hence the one-mode illustration is also not a valid reduction of the central formula. These are not stylistic issues: the abstract and Section 3 explicitly claim that ξ(t) tracks the entanglement dynamics of the JC model, and Eq. (40) is the quantitative basis of that claim.","agreement_with_reader":"agree"},"referee_report":{"model":"deepseek-v4-flash","summary":"The paper constructs a Stratonovich-Weyl Wigner function for the resonant Jaynes-Cummings (JC) model with an initially excited atom and a coherent field, presenting the full hybrid qubit-bosonic quasi-probability distribution in Eq. (23). The authors claim that this Wigner function simplifies to the one-mode Fock-state expression (28), that it reproduces Rabi oscillations and atomic inversion (Eqs. (30)-(35)), and that it yields the purity of the reduced field (Eqs. (39)-(41)) as a measure of entanglement dynamics. The construction of Eq. (23) appears internally consistent and the Rabi/inversion results match the standard JC solution. However, the claimed one-mode reduction and the purity formula are not correct, and these errors directly affect the paper's advertised applications.","tokens_in":12611,"tokens_out":15436,"duration_ms":144232,"significance":"If correct, the explicit closed-form hybrid Wigner function would be a valuable addition to phase-space methods for light-matter interaction, and the purity-based entanglement analysis would provide a practical tool for cavity QED and related settings. The paper deserves credit for a systematic Stratonovich-Weyl construction, for verifying normalization of the full Wigner function, and for reproducing standard JC Rabi oscillations. However, the novel claims in Section 3 are not supported: Eq. (28) is not the specialization of Eq. (23) that the text claims, and Eq. (40) fails the elementary vacuum Fock-state check. These are load-bearing errors because the abstract and conclusions explicitly advertise the purity calculation as the main tool for investigating entanglement dynamics.","major_comments":[{"comment":"The claimed purity of the reduced field is incorrect. For the initial state |e,0> (C_n = δ_{n,0}), the exact JC evolution gives |ψ(t)> = cos(gt)|e,0> − i sin(gt)|g,1>, so the reduced field is cos²(gt)|0><0| + sin²(gt)|1><1| and its purity is cos⁴(gt) + sin⁴(gt), which equals 1/2 at t = π/(4g). Substituting C_n = δ_{n,0} into Eq. (40) gives ξ(t) = 1 for all t, because the double sum contains only the n = m = 0 term and the second sum vanishes. Hence Eq. (40) does not equal tr(ρ_f²) as claimed, and the entanglement analysis built on it is invalid.","section":"Section 3, Eq. (40)"},{"comment":"The asserted reduction of the full Wigner function (23) to the one-mode Fock-state case C_n = δ_{n,r} is not correct. Inserting C_n = δ_{n,r} into (23) yields, in addition to the L_{r,r} term retained in (28), off-diagonal field terms proportional to L_{r,r+1}(2β,2β*) and L_{r+1,r}(2β,2β*) from the atom-field coherence terms, together with a term containing the (r+1)-Fock Wigner function. These extra terms do not vanish; therefore (28) is not the Wigner function of the evolved state |ψ(t)> = cos(tg√(r+1))|e,r> − i sin(tg√(r+1))|g,r+1>. Although the final Rabi probabilities (30)-(31) are correct, their derivation via (28) is not.","section":"Section 3, Eq. (28)"},{"comment":"The long-time limit of the purity is not a valid consequence of the JC dynamics. For the vacuum initial state, the exact purity is cos⁴(gt) + sin⁴(gt), which is periodic and does not tend to a limit; more generally, coherent-state revival dynamics prevent the purity from converging as t → ∞. The claimed asymptotic value ½ + ½ e^{−2|α|²} I₀(2|α|²) therefore cannot be correct, independent of the algebraic error in Eq. (40).","section":"Section 3, Eq. (41)"}],"minor_comments":[{"comment":"The phase factor e^{−itω(E_n−E_m)} has incorrect dimensions because E_n already contains ω; it should be e^{−it(E_n−E_m)} as in Eq. (23).","section":"Eq. (22)"},{"comment":"In the second line, one term contains e^{itE_{n−1}} while the other contains e^{−itE_{n−1}}; the sign appears inconsistent with the rest of the derivation and should be checked.","section":"Eq. (20)"},{"comment":"The revival time expression is typeset ambiguously; the standard formula is T_rev ≈ 2π√⟨N⟩/g.","section":"Eq. (36)"},{"comment":"The identity for the two-dimensional Laguerre polynomial should state explicitly that L_{r,r}(2β,2β*) = (−1)^r r! L_r(4|β|²).","section":"After Eq. (28)"},{"comment":"There are several typographical issues: 'hatDelta' in Eq. (11), 'dyanmics' in the Section 3 heading, 'avors' in Section 3, and 'quentum' in reference [32].","section":"Throughout"}],"recommendation":"reject","confidential_remarks":"I recommend rejection. The errors in Eqs. (28) and (40) are central to the paper's advertised contributions, and correcting them would change the paper's conclusions rather than merely polish the presentation. The authors could, in a future manuscript, rederive the reduced-field purity correctly and reconsider the one-mode Wigner function, but that would require substantial new calculation. The Rabi/inversion results are standard textbook material and do not compensate for the invalid entanglement analysis."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"Quick take: the paper builds a hybrid qubit-field Wigner function for the Jaynes-Cummings model, and the main object (23) is a legitimate closed-form construction. But the two showcase results—the one-mode reduction (28) and the reduced-field purity (40)—are both wrong as printed, and the purity error sinks the entanglement-dynamics claim that the abstract sells.\n\nWhat's actually good: the Stratonovich-Weyl machinery is applied cleanly, the normalization (27) checks out, and the Rabi oscillation and inversion formulas (30)-(35) reproduce the textbook JC results. Equation (23) itself appears to be internally consistent with the density operator (22), modulo a spurious omega in the phase factor of (22) that does not match (23). So there is a real mathematical object here, and the phase-space route to the JC model is worth having.\n\nThe soft spots are not cosmetic. Equation (28) is not the C_n = delta_{n,r} specialization of (23). Inserting the delta's actually produces cross-field Laguerre terms L_{r,r+1}, L_{r+1,r} and a W_{r+1} piece; none appear in (28). So the 'one-mode oscillation' illustration is a different object. More seriously, Eq. (40) fails the Fock check: with C_n = delta_{n,0}, the exact reduced-field purity is cos^4(gt)+sin^4(gt), but (40) returns 1 for all t. That is not a minor slip—it is the quantitative basis for the paper's claim to capture entanglement dynamics. The t→∞ statement (41) is also misleading, because the purity for a coherent state keeps reviving; there is no strict limiting value. These errors do not invalidate (23), but they do invalidate the paper's advertised application.\n\nNovelty is also modest. Reference [24] already built an informationally complete Wigner function for the Tavis-Cummings model, which contains the JC model as the single-atom case. The explicit closed form (23) may still be useful, but the paper overstates what is new.\n\nBottom line: this is a salvageable paper—the core Wigner function is real—but the current version is not ready. A serious referee should see it, but it should be rejected in its present form. If the authors fix the purity computation and the reduction, and recalibrate the claims, it could become a decent contribution to phase-space quantum optics.","headline":"A real closed-form Wigner function for the JC model, but the two advertised applications—the one-mode reduction and the purity formula—are wrong as printed, and the entanglement claims do not hold.","tokens_in":13443,"tokens_out":3346,"would_cite":false,"duration_ms":31551,"reading_group":"no","serious_thinker":"yes","would_accept_peer_review":true},"rs_alignment":null,"lean_confirmation":null,"pith_extraction":{"msc":["81S30","81V80","81R30"],"pacs":[],"model":"deepseek-v4-flash","headline":"The paper derives an explicit Wigner function for the resonant Jaynes-Cummings model on the product phase space of the atomic Bloch sphere and the field's complex plane, and uses it to reproduce Rabi oscillations, collapse and revival…","keywords":["Jaynes-Cummings model","Wigner function","phase space quantum mechanics","Stratonovich-Weyl correspondence","hybrid quantum systems","entanglement dynamics","quasiprobability distribution","coadjoint orbit method"],"falsifier":"A decisive check is to set the initial field to the Fock vacuum, $C_0=1$ and all other $C_n=0$, and evaluate the purity (40) at $t=\\pi/(4g)$. The formula yields $\\xi=1$, while direct computation from the time-evolved state (20) gives a reduced field $\\tfrac{1}{2}(|0\\rangle\\langle0|+|1\\rangle\\langle1|)$ with purity $1/2$; the mismatch would settle the validity of the purity derivation.","tokens_in":11995,"feed_emoji":"⚛️","tokens_out":11694,"duration_ms":96299,"temperature":0.7,"pith_summary":"The paper aims to put the Jaynes-Cummings model, the standard description of a two-level atom coupled to a single quantized field mode, into the phase-space (Wigner function) formalism. It constructs a Wigner function on the product of the atomic Bloch sphere and the field's complex phase plane and claims this quasiprobability distribution is informationally complete: it fully encodes the evolving state of the hybrid system. The authors then recover from it the model's well-known features—Rabi oscillations, the collapse and revival of atomic inversion, and the purity of the reduced field state as a measure of atom-field entanglement. A sympathetic reader should care because, if the construction is right, one closed-form function replaces the full density matrix for this canonical model, potentially simplifying how entanglement is tracked in cavity QED and similar experiments.","feed_headline":"Closed-form Wigner function now captures full Jaynes-Cummings dynamics","feed_subtitle":"A single quasiprobability distribution encodes Rabi oscillations, collapse-revival, and atom–field entanglement.","key_machinery":"The central object is the Stratonovich–Weyl kernel for hybrid qubit–bosonic systems, $\\hat\\Delta(\\Omega)=\\hat\\Delta_q(\\theta,\\phi)\\otimes \\hat\\Delta_f(\\beta,\\beta^*)$, where $\\hat\\Delta_q$ is the spin-$1/2$ kernel built from the SU(2) rotation and parity operator on the Bloch sphere, and $\\hat\\Delta_f$ is the displaced-parity kernel $\\tfrac{2}{\\pi}\\hat D \\hat\\Pi_f \\hat D^\\dagger$ on the field phase plane. The Wigner function is $W_{\\hat\\rho}=\\mathrm{tr}[\\hat\\rho \\hat\\Delta]$. The computations are carried by the Laguerre 2D polynomials $L_{n,m}(\\beta,\\beta^*)$ and their orthogonality relation (26), which turn the trace integrals into the normalization check and the purity formula.","core_discovery":"The authors construct the Wigner function for the resonant Jaynes-Cummings model, with the atom initially excited and the field in a coherent state, by tracing the time-evolved density operator $\\hat\\rho(t)$ against the tensor-product Stratonovich–Weyl kernel $\\hat\\Delta(\\Omega)=\\hat\\Delta_q(\\theta,\\phi)\\otimes \\hat\\Delta_f(\\beta,\\beta^*)$. The resulting quasiprobability distribution, Eq. (23), is a double sum over Fock indices with Laguerre 2D polynomials $L_{n,m}(2\\beta,2\\beta^*)$ and trigonometric Jaynes–Cummings factors; the paper shows it is real and normalized via the identities (25)–(27). From this distribution, the Rabi probability and atomic inversion (33)–(35) match the standard solution, including collapse and revival. Integrating out the spin variables gives the reduced field Wigner function $W_{\\hat\\rho_f}(\\beta,\\beta^*)$ (38), and the purity $\\xi(t)=\\pi\\int|W_{\\hat\\rho_f}|^2 d^2\\beta$ (39) evaluates to (40), whose long-time limit is $\\tfrac{1}{2}+\\tfrac{1}{2} e^{-2|\\alpha|^2}I_0(2|\\alpha|^2)$. The central claim is that this phase-space object is informationally complete: it carries the full dynamics of the light–matter interaction, including entanglement dynamics, without reconstructing the density matrix.","pith_inferences":["Since the Stratonovich–Weyl kernel is fixed by the symmetry groups (SU(2) for the qubit and the Heisenberg–Weyl group for the field), the same construction should extend to other hybrid systems—for example, multi-mode fields or the Tavis–Cummings model—by replacing only the time-evolved density operator.","The purity expression depends on the Laguerre orthogonality relation; a direct finite-sum check with two-Fock initial states or a very small coherent amplitude would make the derivation's validity transparent and could be done with standard symbolic computation.","If the Wigner function is truly informationally complete, it should support a star product for the hybrid algebra, turning the Jaynes-Cummings model into a fully deformation-quantized system; the paper anticipates this direction but does not construct it.","The phase-space purity measure could be compared experimentally with the von Neumann entropy of the field mode in cavity QED, since the field Wigner function is directly measurable, providing a practical entanglement witness."],"forward_implications":["The Wigner function (23) is real and normalized by construction, so expectation values of atom-field observables can be computed by phase-space integrals instead of operator traces.","The Rabi probabilities and atomic inversion extracted from the Wigner function, Eqs. (33)–(35), reproduce the standard Jaynes-Cummings result, including the collapse–revival times $T_{\\mathrm{rev}}\\approx 2\\pi k/(g\\sqrt{\\langle \\hat N\\rangle})$.","The reduced field Wigner function (38) and purity (40) allow entanglement dynamics to be read off from the field's quasiprobability distribution, with the long-time purity approaching $\\tfrac{1}{2}+\\tfrac{1}{2} e^{-2|\\alpha|^2}I_0(2|\\alpha|^2)$.","In the single-mode (Fock-state) limit, the Wigner function reduces to the known number-state Wigner function multiplied by a factor describing Rabi oscillations at frequency $2g\\sqrt{r+1}$, isolating the elementary excitation-exchange step."],"supporting_citations":[{"why":"Defines the Jaynes–Cummings model and its standard solution, which the Wigner function must reproduce.","marker":"[19, 20]"},{"why":"Supplies the Stratonovich–Weyl postulates that justify calling the kernel-based map a valid phase-space representation.","marker":"[30]"},{"why":"Provides the spin-1/2 Stratonovich–Weyl kernel on the Bloch sphere used for the qubit part.","marker":"[31, 32]"},{"why":"Establishes the hybrid qubit–bosonic Wigner function construction and the negativity-volume entanglement witness that motivates the purity analysis.","marker":"[24, 25]"},{"why":"Gives the Laguerre 2D-function identity and orthogonality relation used for normalization and purity evaluation.","marker":"[41]"},{"why":"Supplies the resonant Jaynes–Cummings evolution operator and the standard Rabi/collapse-revival results used as the consistency benchmark.","marker":"[23]"},{"why":"Defines purity in terms of the Wigner function integral, which is the basis for the entanglement measure.","marker":"[46]"},{"why":"Defines the Laguerre 2D polynomials appearing in the Wigner function's closed form.","marker":"[40]"}],"fun_headline_variants":["Exact Wigner function for Jaynes-Cummings dynamics","Full phase-space Wigner function captures atom-field entanglement","Quasiprobability encodes Rabi oscillations and purity","Wigner function solves light-matter interaction in phase space","Phase-space Wigner function reveals full Jaynes-Cummings dynamics"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The load-bearing assumption is that inserting Fock-state amplitudes into the full Wigner function yields the single-mode expression (28) and that the stated polynomial orthogonality turns the purity integral into the closed form (40); if either step is wrong, the claimed results collapse.","fun_headline_variants_meta":{"raw":{"variants":["Exact Wigner function for Jaynes-Cummings dynamics","Full phase-space Wigner function captures atom-field entanglement","Quasiprobability encodes Rabi oscillations and purity","Wigner function solves light-matter interaction in phase space","Phase-space Wigner function reveals full Jaynes-Cummings dynamics"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000158,"raw_usage":{"total_tokens":1240,"prompt_tokens":974,"completion_tokens":266,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":590,"completion_tokens_details":{"reasoning_tokens":183}},"tokens_in":590,"tokens_out":266,"duration_ms":3209,"temperature":1.0,"reasoning_tokens":183,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-06T21:50:54.652588+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"A decisive check is to set the initial field to the Fock vacuum, $C_0=1$ and all other $C_n=0$, and evaluate the purity (40) at $t=\\pi/(4g)$. The formula yields $\\xi=1$, while direct computation from the time-evolved state (20) gives a reduced field $\\tfrac{1}{2}(|0\\rangle\\langle0|+|1\\rangle\\langle1|)$ with purity $1/2$; the mismatch would settle the validity of the purity derivation.","supporting_citations":[{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Supplies the Stratonovich–Weyl postulates that justify calling the kernel-based map a valid phase-space representation."},{"cited_title":"W¨ unsche, Laguerre 2D-functions and their application in quantum optics , J","cited_arxiv_id":null,"evidence_quote":"Gives the Laguerre 2D-function identity and orthogonality relation used for normalization and purity evaluation."},{"cited_title":"Meystre and M","cited_arxiv_id":null,"evidence_quote":"Supplies the resonant Jaynes–Cummings evolution operator and the standard Rabi/collapse-revival results used as the consistency benchmark."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Defines purity in terms of the Wigner function integral, which is the basis for the entanglement measure."},{"cited_title":"W¨ unsche, Generating Functions for Products of Special Laguerre 2D and Hermite 2D Polynomials , Appl","cited_arxiv_id":null,"evidence_quote":"Defines the Laguerre 2D polynomials appearing in the Wigner function's closed form."}],"review_version":1}