{"id":"36a3a931-3c2c-451b-a339-bcf97a853782","arxiv_id":"1908.06594","paper_version":2,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":5.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":3,"one_line_summary":"Two cavity-QED dissipative protocols deterministically prepare the maximally discordant two-qubit mixed state as a unique steady state.","lead":"This paper proposes two dissipative cavity-QED schemes that prepare a two-qubit separable state with maximal quantum discord. Both schemes are designed so this non-entangled but non-classically correlated state is the unique steady state of the engineered dynamics, with simulated fidelities above 99%.","discovery_kind":"new_application","skeptic_critique":{"model":"deepseek-v4-flash","headline":"First-scheme switching interval is below 1/κ for the advertised 87Rb parameters, so the 99.41% fidelity may not follow from the Trotter-derived master equation.","rationale":"The central claim rests on the effective master equations (13) and (22) being faithful descriptions of the physical sequences. The second scheme is supported by full-vs-effective agreement (Fig. 7) and does not rely on rapid phase switching, so I do not see a comparably sharp problem there. The first scheme's Trotter step is the weakest link: the authors themselves state the interval must be much longer than 1/κ, and the advertised 87Rb parameters put the default horizon below that threshold. This is not a dispute with the general dissipative-preparation idea or with the algebra of the steady state; it is a concrete parameter-regime inconsistency that determines whether the quoted fidelity demonstrates the claim. The paper's numerical comparisons (Figs. 2,3) are credible within their stated parameters, but those parameters (κ=0.1g, etc.) are different from the experimental κ≈0.046g used in Sec. V. The reader's weakest assumption identified the same Trotter timescale issue, and this check makes it quantitative. Because the second scheme and the general formalism remain sound, the appropriate verdict stays CONDITIONAL; the concern does not overturn the paper but adds a specific condition on the first-scheme feasibility claim.","tokens_in":15886,"tokens_out":15567,"duration_ms":160963,"concrete_test":"Take the 87Rb parameters from Sec. V for the first scheme (g=2π×14.4 MHz, κ=2π×0.66 MHz, Ω=0.3g, Δ=76g, N=200) and rerun the full master equation built from Hamiltonian (6) with total time T=40000/g so that the switching interval T/(2N)=100/g exceeds 1/κ≈21.8/g by a factor ~4.6 and remains below 1/γ≈735/g. Compare the final super-fidelity G(ρ(t),σ) with the 99.41% value quoted in Sec. V and with the prediction of the effective master equation (13). If the full-model fidelity is still >99% and matches Eq. (13), the concern is resolved; if it drops or disagrees, the first-scheme experimental claim is not supported by the derived effective dynamics.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The load-bearing step is the Trotterization in Sec. III: Eq. (12) represents each phase interval by e^{L_x T/2N} and e^{L_y T/2N}, which requires the interval I = T/(2N) to satisfy 1/κ ≪ I ≪ 1/γ (stated below Eq. (12)). For the 87Rb numbers quoted in Sec. V, g = 2π×14.4 MHz, κ = 2π×0.66 MHz, Ω = 0.3g, Δ = 76g, and N = 200. Then G = gΩ/Δ ≈ 0.00395g, γ = 4G^2/κ ≈ 1.36×10^-3 g, so 1/γ ≈ 735/g and 1/κ ≈ 21.8/g. The rest of the paper uses gt = 8000 for its simulations; with that horizon I = 20/g, so κI ≈ 0.92, not ≫1. The cavity field has not relaxed within each interval, so the physical process is not the product of the two effective Lindblad semigroups used in Eq. (12). Consequently the 99.41% fidelity quoted for the first scheme may be obtained in a parameter window where Eq. (13) is not justified, rather than demonstrating the claimed dissipative preparation. The coupled-cavity scheme does not have this Trotter step and appears unaffected; the concern is specific to the first scheme's experimental-feasibility claim.","agreement_with_reader":"agree"},"referee_report":{"model":"deepseek-v4-flash","summary":"The paper proposes two cavity-QED schemes for the deterministic dissipative preparation of the two-qubit maximally discordant mixed state ρ = (|Ψ+⟩⟨Ψ+| + |00⟩⟨00| + |11⟩⟨11|)/3. In the first scheme, two four-level atoms in a single lossy cavity are driven by classical fields whose phases are switched alternately, and the authors use a Trotter product formula to derive the effective master equation (13), which is a sum of two collective Lindblad dissipators Sx and Sy with equal weights. In the second scheme, two atoms are placed in a lossy coupled-cavity array, and delocalized modes are used to realize the same two Lindblad terms simultaneously, yielding Eq. (22). The target state is shown to be a steady state of these Lindblad equations, and numerical simulations of the full master equations are compared with the effective dynamics for various parameters, atomic spontaneous emission rates, initial states, and phase mismatches. The paper also provides an experimental feasibility discussion based on 87Rb parameters and reports fidelities above 99% for both schemes.","tokens_in":16220,"tokens_out":9338,"duration_ms":94068,"significance":"The work is significant if the derivations hold, because it offers concrete cavity-QED implementations of a known but experimentally elusive target state, using dissipation as a resource rather than a nuisance. The paper has clear strengths: the full-model simulations reproduce the effective dynamics in Figs. 2 and 7, the microscopic Hamiltonians are given explicitly, no parameter is fitted to the target state, and the phase-mismatch analysis adds a useful practical perspective. The coupled-cavity scheme in Sec. IV is the more robust part of the paper, since it avoids the Trotterization step that is problematic in the first scheme. However, the first scheme's experimental-feasibility claim rests on a parameter regime where the Trotter timescale separation is violated, and the uniqueness of the steady state is asserted rather than proved. These issues are load-bearing for the paper's central claims and require a major revision.","major_comments":[{"comment":"The Trotterization step that produces Eq. (13) is not justified in the experimental parameter window quoted for the first scheme. The text below Eq. (12) requires each switching interval I = T/(2N) to satisfy 1/κ ≪ I ≪ 1/γ. For the 87Rb numbers in Sec. V (g = 2π×14.4 MHz, κ = 2π×0.66 MHz, Ω = 0.3g, Δ = 76g, N = 200, and gt = 8000), one obtains G ≈ 0.00395g, γ = 4G²/κ ≈ 1.36×10⁻³g, so 1/κ ≈ 21.8/g, while I = 20/g. Thus κI ≈ 0.92, which is not ≫1; the cavity field has not relaxed within each phase interval, so the physical map is not a concatenation of e^{L_x I} and e^{L_y I} as in Eq. (12). The 99.41% fidelity quoted for the first scheme therefore does not demonstrate the dissipative-preparation mechanism of Eq. (13) in that parameter window. The authors should either choose parameters (for example, a longer total evolution time or a smaller N) that satisfy κI ≫ 1 and recompute the fidelity, or explicitly present the first scheme as a stroboscopic simulation of the full Hamiltonian and derive the corresponding validity conditions.","section":"Sec. III and Sec. V"},{"comment":"The uniqueness claim that underlies the deterministic-preparation statement is asserted but not proved. The appendix states that Eq. (A.3) 'can be examined both numerically and analytically' to have a unique steady state except for δϕ1 = −δϕ2 = ±0.5π, but no analytic argument is given and the numerical check is not shown. Since independence of the initial state (except for the singlet |Ψ−⟩) is a central advertised advantage of both schemes, the authors should provide a complete proof of uniqueness for Eq. (23), and likewise for Eq. (13) where uniqueness is also stated without proof. Citing a known theorem for the uniqueness of steady states of Lindblad equations with a given set of jump operators would be a sufficient resolution.","section":"Appendix"},{"comment":"The parameter scans include regimes where the adiabatic-elimination condition κ ≫ G is violated. With Δ = 100g and Ω = 0.5g, one has G = 0.005g, and the curve for κ = 0.01g in Fig. 3(a) therefore corresponds to κ/G = 2, so the effective master equation (10) is not valid in that regime. If the full master equation still reaches high fidelity there, the mechanism cannot be attributed to the dissipative dynamics described by Eq. (10), and the interpretation of that curve as a test of the effective model is misleading. The parameter discussion should be restricted to κ ≫ G, or the non-adiabatic regime should be analyzed separately.","section":"Sec. III, Fig. 3(a)"}],"minor_comments":[{"comment":"Equation (A.2) contains a stray ρ inside the Lindblad operator: the term '−ie^{iδϕ2}|2⟩⟨3|ρ' should read '−ie^{iδϕ2}|2⟩⟨3|', since the state ρ should only appear outside the Lindblad generator.","section":"Appendix, Eq. (A.2)"},{"comment":"In the Introduction, 'Glave et al.' should be 'Galve et al.', as the reference is to F. Galve, G. L. Giorgi, and R. Zambrini.","section":"Sec. I"},{"comment":"The symbol γ is used for both the collective decay rate (e.g., γy = 4G²/κ in Eq. (10)) and the atomic spontaneous emission rate (γ = 2π×3 MHz in Sec. V). This overloading is confusing in Figs. 3(b), 5, and 8(d), and the authors should distinguish the two quantities, for example by using Γ for the collective rate.","section":"Secs. III and V"},{"comment":"The phrase 'V on neumann entrophy' in Sec. II is a typo for 'Von Neumann entropy'.","section":"Sec. II"}],"recommendation":"major_revision","confidential_remarks":"The paper contains a sound and useful coupled-cavity scheme, but the first scheme's main numerical claim is made in a parameter regime where the Trotter-derived master equation is not justified. The uniqueness proof is also missing. These issues can be fixed within the scope of a revision, so I recommend major revision rather than rejection."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"Dear colleague,\n\nThe short version: this is a legitimate, workmanlike contribution to dissipative preparation of non-entangled quantum correlations. It provides two concrete cavity-QED implementations for the known maximally discordant mixed state (MDMS), with full-model numerics supporting the effective master equations. The coupled-cavity scheme is clean and convincing. The first scheme has a questionable parameter choice for its headline experimental numbers.\n\nWhat's actually new: the specific atom-cavity realizations of the Sx/Sy Lindblad dynamics, including the alternating-phase Trotterized sequence and the coupled-cavity version. The target state and Eq. (4) come from earlier work (Galve et al., López et al.), so the novelty lies in the physical implementations, which is legitimate. The derivations are standard adiabatic elimination, and the full-model vs effective-model comparisons in Figs. 2 and 7 are the right evidence. The sensitivity analysis for atomic spontaneous emission and phase mismatch is useful and honest.\n\nSoft spots, in order of importance:\n\n1. The stress-test on the first scheme's 87Rb parameters is correct. The paper states that the switching interval I must satisfy 1/κ ≪ I ≪ 1/γ. With g=2π×14.4 MHz, κ=2π×0.66 MHz, Ω=0.3g, Δ=76g, N=200, and T=8000/g, I=20/g, while 1/κ≈21.8/g. So I is actually smaller than 1/κ, not much larger. The authors choose N=200 and report 99.41% fidelity for these parameters, but that window does not meet their own criterion for Trotterized adiabatic elimination. The high fidelity is therefore not adequately explained by the effective master equation (13). It might still be a genuine physical result, but the paper doesn't analyze why it works outside the intended regime. The coupled-cavity scheme does not have this problem.\n\n2. The appendix claims uniqueness of the steady state for Eq. (23) but only offers a 'numerically and analytically' assertion with a messy matrix equation, not a proof. That's a gap, though not fatal.\n\n3. No comparison with Ref. [42] (Altintas et al.), which appears to have considered MDMS steady states. The reader flagged this, and it's worth addressing to establish what's genuinely new.\n\nThese are all addressable in revision. The math and numerics otherwise hold up; the effective-model derivations are reproducible, and the full-model simulations are the right check.\n\nWho this is for: people working on dissipative state preparation in cavity QED, especially those interested in non-entangled correlations. It deserves a serious referee, not a desk rejection. I'd recommend sending it to review, with the expectation that the first-scheme timescale issue gets fixed or acknowledged, and the comparison to Ref. [42] is added.\n\nBest,\n[Your name]","headline":"Solid dissipative-state-preparation paper with one real flaw: the first scheme's 87Rb parameters violate the paper's own interval condition, so the 99.41% fidelity is not a clean demonstration of the Trotterized effective dynamics.","tokens_in":16727,"tokens_out":7928,"would_cite":true,"duration_ms":74965,"reading_group":"maybe","serious_thinker":"yes","would_accept_peer_review":true},"rs_alignment":null,"lean_confirmation":null,"pith_extraction":{"msc":[],"pacs":["03.67.-a","42.50.Pq"],"model":"deepseek-v4-flash","headline":"Two cavity-QED protocols use engineered dissipation to drive two qubits into the maximally discordant mixed state, the unique steady state in the subspace orthogonal to the singlet, with simulated fidelity above 99 percent.","keywords":["quantum discord","maximally discordant mixed state","dissipative state preparation","cavity quantum electrodynamics","collective decay operators","Lindblad master equation","super-fidelity","separable quantum correlations"],"falsifier":"Solve the full coupled-cavity master equation without the adiabatic-elimination approximation for the parameters of Fig. 7 ($\\kappa=0.1g$, $\\Omega=0.2g$, $\\Delta=100g$) from the initial state $|00\\rangle|00\\rangle_c$: the central claim predicts a late-time super-fidelity above 99 percent to $\\rho=(|\\Psi^+\\rangle\\langle\\Psi^+|+|00\\rangle\\langle00|+|11\\rangle\\langle11|)/3$; a substantial deviation, or a steady-state discord below $1/3$, would falsify it.","tokens_in":15692,"feed_emoji":"⚛️","tokens_out":18732,"duration_ms":169325,"temperature":0.7,"pith_summary":"The paper aims to show that the maximally discordant mixed state of two qubits, $\\rho=(|\\Psi^+\\rangle\\langle\\Psi^+|+|00\\rangle\\langle00|+|11\\rangle\\langle11|)/3$, can be generated deterministically by treating cavity loss as an engineered resource rather than as an error. Because discord survives in separable mixed states and is more resistant to decoherence than entanglement, a fixed dissipative route to the maximally discordant state would supply a practical non-classical resource. The paper derives two effective master equations, one from a phase-switched single-cavity setup and one from a lossy coupled-cavity setup, and shows both have this $\\rho$ as their unique steady state in the subspace orthogonal to the singlet $|\\Psi^-\\rangle$. Numerical simulation of the full models with currently available cavity-QED parameters gives super-fidelities above 99 percent, with steady-state concurrence zero and quantum discord $1/3$.","feed_headline":"Cavity-QED schemes prepare a maximally discordant state by dissipation","feed_subtitle":"Both designs steer any initial state except the singlet to a discord-1/3 steady state with fidelities above 99 percent.","key_machinery":"The load-bearing object is the pair of collective decay operators $S_x=\\sigma_x^1+\\sigma_x^2$ and $S_y=\\sigma_y^1+\\sigma_y^2$: the target state is the common steady state of the Lindblad generators $L[S_x]$ and $L[S_y]$. The paper realizes these operators in two physical arrangements. In the single-cavity scheme, two four-level atoms driven by phase-controlled classical fields and coupled to a lossy cavity are reduced by large detunings to an effective coupling $G(J_-+J_+)a^\\dagger$; strong cavity decay converts this into the collective decay $L[S_y]$ (or $L[S_x]$ with appropriate drive phases), and alternating the phases at a fast rate superposes the two generators through the Trotter product formula, giving Eq. (13). In the coupled-cavity scheme, the transformation to delocalized modes $m_1=(a_1-a_2)/\\sqrt{2}$ and $m_2=(a_1+a_2)/\\sqrt{2}$ makes one mode couple to $S_y$ and the other to $S_x$; adiabatic elimination of these lossy modes gives both Lindblad terms at once, with $\\gamma_x=\\gamma_y=4G^2/\\kappa$, yielding Eq. (22).","core_discovery":"The central discovery is that collective dissipative dynamics generated by the two Lindblad operators $S_x=\\sigma_x^1+\\sigma_x^2$ and $S_y=\\sigma_y^1+\\sigma_y^2$ select the rank-3 maximally discordant mixed state $\\rho=(|\\Psi^+\\rangle\\langle\\Psi^+|+|00\\rangle\\langle00|+|11\\rangle\\langle11|)/3$ as their unique steady state in the subspace orthogonal to the singlet $|\\Psi^-\\rangle$, and that these Lindblad terms can be realized physically. The single-cavity scheme alternates the phases of the two classical drives; the Trotter limit of fast switching yields $\\dot\\rho=\\frac{1}{2} L_{\\gamma_x}[S_x]\\rho+\\frac{1}{2} L_{\\gamma_y}[S_y]\\rho$ (Eq. 13). The coupled-cavity scheme uses delocalized cavity modes to couple to $S_x$ and $S_y$ simultaneously, giving $\\dot\\rho=L_{\\gamma_x}[S_x]\\rho+L_{\\gamma_y}[S_y]\\rho$ after adiabatic elimination (Eq. 22). In both cases the cavity decay rate $\\kappa$ enters through $\\gamma=4G^2/\\kappa$ and is the engine of the protocol. Full master-equation simulations show the steady state has concurrence 0, quantum discord 1/3, and super-fidelity above 99 percent under realistic parameters, and is reached from any initial state that does not populate $|\\Psi^-\\rangle$.","pith_inferences":["A natural testable extension is to apply the same Trotterized phase-switching construction to other pairs of collective decay operators, whose common fixed points would form new families of steady-state quantum-correlated states.","The delocalized-mode trick suggests a general recipe: a cavity-hopping term can be converted into an engineered collective dissipator by working in symmetric and antisymmetric modes, a strategy that could carry over to multi-qubit and higher-dimensional systems.","One could also add a third dissipative channel that damps the singlet, which would extend preparation to genuinely all initial states and remove the paper's stated restriction."],"forward_implications":["Any initial state with zero overlap on the singlet $|\\Psi^-\\rangle$ converges to the same maximally discordant mixed state, so neither precise state preparation nor operation-time control is needed; this removes the main limitation of unitary MDMS generation.","The prepared steady state has zero concurrence but quantum discord $1/3$, providing a deterministic source of non-classical correlations that do not rely on entanglement.","Because the protocol is driven by cavity decay through $\\gamma=4G^2/\\kappa$, faster cavity loss accelerates convergence, while atomic spontaneous emission only weakly degrades the final fidelity (still above 99 percent at $\\gamma\\sim g$).","Both schemes work with currently available cavity-QED parameters, with reported super-fidelities from 99.10 percent to 99.67 percent across three experimental platforms.","The phase-mismatch analysis shows the target state remains the unique steady state for a broader family of collective Lindblad operators, except for the singular combination $\\delta\\phi_1=-\\delta\\phi_2=\\pm0.5\\pi$, so the drive phases need not be perfectly tuned."],"supporting_citations":[{"why":"defines the rank-2 and rank-3 maximally discordant mixed states and identifies $\\epsilon=1/3$, $x=1/2$ as the target state used in the paper.","marker":"[12]"},{"why":"rewrites the target state in the Bell basis and supplies the master equation whose steady state the paper aims to engineer.","marker":"[13]"},{"why":"provides the four-level atom-cavity configuration and the large-detuning route to collective dissipative terms.","marker":"[31]"},{"why":"supplies the Trotter product formula (Corollary 5.8) that justifies the fast phase-switching limit and yields Eq. (13).","marker":"[33]"},{"why":"introduces the delocalized cavity modes $m_1,m_2$ used to obtain the coupled-cavity effective Hamiltonian (20).","marker":"[47]"},{"why":"defines super-fidelity, the metric used to compare the simulated mixed state with the target.","marker":"[40]"},{"why":"supplies projected high-cooperativity cavity parameters used in the experimental-feasibility fidelity estimates.","marker":"[52]"},{"why":"gives microscopic optical-resonator parameters used for the second experimental-feasibility estimate.","marker":"[53]"}],"fun_headline_variants":["Dissipation steers qubits to a maximally discordant steady state","Deterministic discord without entanglement, via cavity decay","Cavity loss creates maximally discordant mixed states with 99% fidelity","Dissipative protocols yield unique discord-1/3 steady state","Max discord from dissipation: cavity QED achieves 99% fidelity"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"Both derivations collapse unless the cavity loses photons much faster than the atoms exchange excitations with it, and the phase-switched scheme additionally needs each dwell time to be longer than the cavity decay time yet shorter than the collective decay time, so the two decay mechanisms cleanly alternate.","fun_headline_variants_meta":{"raw":{"variants":["Dissipation steers qubits to a maximally discordant steady state","Deterministic discord without entanglement, via cavity decay","Cavity loss creates maximally discordant mixed states with 99% fidelity","Dissipative protocols yield unique discord-1/3 steady state","Max discord from dissipation: cavity QED achieves 99% fidelity"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000606,"raw_usage":{"total_tokens":2905,"prompt_tokens":1102,"completion_tokens":1803,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":718,"completion_tokens_details":{"reasoning_tokens":1721}},"tokens_in":718,"tokens_out":1803,"duration_ms":11409,"temperature":1.0,"reasoning_tokens":1721,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-14T12:40:57.157382+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Solve the full coupled-cavity master equation without the adiabatic-elimination approximation for the parameters of Fig. 7 ($\\kappa=0.1g$, $\\Omega=0.2g$, $\\Delta=100g$) from the initial state $|00\\rangle|00\\rangle_c$: the central claim predicts a late-time super-fidelity above 99 percent to $\\rho=(|\\Psi^+\\rangle\\langle\\Psi^+|+|00\\rangle\\langle00|+|11\\rangle\\langle11|)/3$; a substantial deviation, or a steady-state discord below $1/3$, would falsify it.","supporting_citations":[{"cited_title":"Experimental one-way quantum computing,","cited_arxiv_id":null,"evidence_quote":"defines the rank-2 and rank-3 maximally discordant mixed states and identifies $\\epsilon=1/3$, $x=1/2$ as the target state used in the paper."},{"cited_title":"The quantum dissipation characterized by a Lindblad gen- erator in Markovian quantum master equations originates from the weak coupling between quantum systems and envi- ronment","cited_arxiv_id":null,"evidence_quote":"rewrites the target state in the Bell basis and supplies the master equation whose steady state the paper aims to engineer."},{"cited_title":"Unconventional rydberg pumping and applications in quantum information processing,","cited_arxiv_id":null,"evidence_quote":"provides the four-level atom-cavity configuration and the large-detuning route to collective dissipative terms."},{"cited_title":"Scheme for entanglement generation in an atom-cavi ty system via dissipation,","cited_arxiv_id":null,"evidence_quote":"supplies the Trotter product formula (Corollary 5.8) that justifies the fast phase-switching limit and yields Eq. (13)."},{"cited_title":"Sub- and super-ﬁdelity as bounds for quantum ﬁdelity,","cited_arxiv_id":null,"evidence_quote":"introduces the delocalized cavity modes $m_1,m_2$ used to obtain the coupled-cavity effective Hamiltonian (20)."},{"cited_title":"Engel and R","cited_arxiv_id":null,"evidence_quote":"defines super-fidelity, the metric used to compare the simulated mixed state with the target."},{"cited_title":"Multimode entanglement in cou- pled cavity arrays,","cited_arxiv_id":null,"evidence_quote":"supplies projected high-cooperativity cavity parameters used in the experimental-feasibility fidelity estimates."},{"cited_title":"Quantum many-body phenomena in coupled cavity arrays,","cited_arxiv_id":null,"evidence_quote":"gives microscopic optical-resonator parameters used for the second experimental-feasibility estimate."}],"review_version":1}