{"id":"0a01b28a-2185-4ec0-9c3b-98f237efd1bd","arxiv_id":"2504.15559","paper_version":1,"verdict":"CONDITIONAL","confidence":"HIGH","novelty_score":4.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":6,"one_line_summary":"Numerical simulations show magnon blockade, g^(2)(0) ≈ 0.04, in a strongly dispersive qubit-magnon system, caused by qubit-induced anharmonicity.","lead":"This paper shows that a superconducting qubit can block magnons, the quantized magnetic excitations in a magnetic sphere, from appearing in pairs when the qubit is far detuned from the magnon mode. The result suggests a path to single-magnon control for quantum sensing, but it depends on model parameters that are more optimistic than current experiments.","discovery_kind":"new_application","skeptic_critique":{"model":"deepseek-v4-flash","headline":"Effective dispersive Hamiltonian may be used outside its validity regime; missing g_qm and Δ_qm leaves the central blockade claim unverified.","rationale":"The reader's weakest assumption identifies the same load-bearing issue: the paper operates in a parameter range where χ_qm/γ reaches 45 but never reports the underlying g_qm and Δ_qm, even though Eq. (4) is derived under Δ_qm ≫ g_qm. My independent reading of the model confirms this is the most consequential gap. The rest of the analysis is internally consistent: the rotating-frame Hamiltonian in Eq. (9) is standard, the resonance conditions in Fig. 3 are compatible with the dressed-level structure of the dispersive model, and the numerical method (Lindblad master equation via QuTiP) is appropriate for the stated problem. The thermal-noise discussion is also internally coherent, with m_th ≈ 0.0035 at T ≈ 72 mK for ω_m/2π = 8.5 GHz. I do not see an internal mathematical contradiction that would force rejection; the concern is that the effective model may not describe the physical qubit–magnon system at the quoted strong-dispersive parameters. This is addressable by reporting g_qm and Δ_qm and by comparing against the full Jaynes–Cummings model at the same χ_qm. If the comparison confirms the effective-model results, the remaining issue is experimental feasibility rather than correctness; if it does not, the central claim is substantially weakened. The conditional verdict is therefore appropriate, and my stress-test pass does not change it.","tokens_in":11860,"tokens_out":12039,"duration_ms":120130,"concrete_test":"Recompute the steady-state g^(2)(0) for the points in Figs. 2 and 3 using the full Jaynes–Cummings interaction H_int = g_qm(m†σ_− + mσ_+) together with the same free terms, drives, and dissipators, with g_qm and Δ_qm chosen to reproduce each quoted χ_qm. For example, set Δ_qm/γ = 10 g_qm/γ and g_qm/γ = sqrt(χ_qm Δ_qm/γ²); for χ_qm/γ = 45 this gives g_qm/γ = 450, Δ_qm/γ = 4500. Repeat with a second, less dispersive ratio Δ_qm = 3g_qm to bracket the edge of the expansion. If the optimal blockade depth or the detuning positions shift by more than one linewidth κ_m/γ = 1.4 relative to the effective-model curves, Eq. (4) is being used outside its validity regime and the headline claim needs to be qualified.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The central numerical claim, g^(2)(0) → 0.04 for χ_qm/γ = 20–45, rests on Eq. (4), the second-order dispersive Hamiltonian derived under the explicit condition Δ_qm ≡ |ω_q − ω_m| ≫ g_qm (Section II, just before Eq. 4). The paper never reports g_qm or Δ_qm, only χ_qm = g_qm²/Δ_qm. Because χ_qm/γ is scanned up to 45 while the dissipative scales are κ_m/γ = 1.4 and κ_q/γ = 1.2, the strong-dispersive requirement χ_qm ≫ κ forces the underlying ratio Δ_qm/g_qm to be modest unless g_qm is very large. Concretely, χ_qm/γ = 45 with Δ_qm = 10g_qm requires g_qm/γ = 450; with the experimentally cited g_qm ~ 10γ, χ_qm/γ = 45 is mathematically impossible because χ_qm = g_qm²/Δ_qm ≤ g_qm when Δ_qm ≫ g_qm. Higher-order corrections to Eq. (4) are of relative order (g_qm/Δ_qm)²; at Δ_qm/g_qm = 3 this is ~11% of χ_qm, i.e., about 5γ for the largest plotted point, comparable to or larger than the linewidths that set the blockade line shapes. The blockade depth and resonance positions in Figs. 2 and 3 are precisely governed by level shifts of order χ_qm, so a 10–25% correction to the effective Hamiltonian is not an innocuous renormalization. Without reporting g_qm and Δ_qm, the paper does not establish that the predicted blockade is a property of the full qubit–magnon system rather than an artifact of using Eq. (4) outside its derivation regime.","agreement_with_reader":"agree"},"referee_report":{"model":"deepseek-v4-flash","summary":"The manuscript studies magnon blockade in a hybrid superconducting-qubit/yttrium-iron-garnet system operating in the dispersive regime. Starting from a cavity-mediated qubit-magnon interaction, the authors adopt the dispersive Hamiltonian H_disp = χ_qm(m†m + 1/2)σ_z (Eq. 4) with χ_qm = g_qm²/Δ_qm, add separate drives on the magnon and qubit, and solve the Lindblad master equation for the steady state. The main results are: (i) the second-order correlation function g^(2)(0) reaches about 0.04 for χ_qm/γ ≈ 20–40 at appropriate magnon drive detunings; (ii) the blockade is attributed to the qubit-induced anharmonicity that suppresses two-magnon transitions; (iii) the resonance conditions obtained by diagonalizing the driven qubit part match the numerically observed dips; and (iv) thermal magnon noise of mth ≈ 0.0035 (corresponding to T ≈ 72 mK for ω_m/2π = 8.5 GHz) destroys the blockade, while qubit thermal noise is less influential.","tokens_in":12293,"tokens_out":15085,"duration_ms":134468,"significance":"If the predictions are correct, the paper extends magnon blockade into the dispersive regime, which is the regime used for quantum non-demolition magnon readout in several recent experiments. The work therefore suggests that a single architecture can both read out and manipulate single magnons. The manuscript provides standard, internally consistent master-equation simulations and gives explicit analytic resonance formulas that are compared with the numerics; the thermal-robustness analysis yields a concrete, testable temperature threshold. The central reservation is that the effective dispersive Hamiltonian is used for parameter values whose consistency with the derivation condition |ω_q−ω_m| ≫ g_qm is not established, so the quantitative prediction of g^(2)(0) ≈ 0.04 may be an artifact of the model rather than a property of the physical system.","major_comments":[{"comment":"The paper does not report the Fock-space truncation dimension used in the QuTiP master-equation simulations, nor does it provide any convergence test. The quantity g^(2)(0) is computed from the two-magnon population and is sensitive to the truncation of high-number magnon states, particularly in the finite-temperature runs of Fig. 4 where mth > 0 populates higher Fock states. Please state the maximum magnon number retained for each figure and confirm that g^(2)(0) and P_2 are converged with respect to increasing the truncation cutoff.","section":"Section III, numerical simulations"}],"minor_comments":[{"comment":"The text lists the single-magnon resonance positions as ≈ ±20 and ≈ ±56 for the parameters of Fig. 3 (χ_qm/γ = 45, Δ_q/γ = −20, Ω_s/γ = 15). Evaluating the displayed formulas with these parameters gives ≈ ±23.3 and ≈ ±48.3, respectively; please reconcile this discrepancy or correct the quoted numbers/formulas.","section":"Section III, resonance positions"},{"comment":"The sentence preceding Eq. (3) states that the qubit and Kittel mode are nearly resonant when |ω_q−ω_m| ≪ g_cq, g_cm. This condition compares a frequency detuning to coupling strengths rather than to a spectral width; please clarify the intended dimensionless statement or provide a more explicit derivation of the cavity-mediated qubit–magnon coupling.","section":"Section II, Eq. (3)"},{"comment":"The symbol γ is used both for the gyromagnetic ratio (γ/2π = 28 GHz/T in Section II) and as the frequency scale (γ = 2π × 1 MHz in the figure captions and parameter values). This dual use is confusing; please use a different symbol, such as κ_0 or Γ, for the scaling unit in the figures.","section":"Notation"}],"recommendation":"major_revision","confidential_remarks":"The parameter ranges used (χ/γ up to 45, i.e., χ/2π up to ≈45 MHz) appear far larger than the dispersive shifts of order MHz reported in the qubit–magnon experiments cited by the authors. The manuscript would benefit from a concrete statement of which experimental parameters could realize the predicted regime, or an explicit discussion of whether such large χ values are achievable."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"The authors numerically study magnon blockade in a qubit-magnon system where the qubit is detuned from the magnon, so the interaction is the standard dispersive shift χ m†m σ_z. They solve the Lindblad master equation with QuTiP and show g^(2)(0) can drop to ~0.04 for χ/γ ~ 20–45, with a matching resonance structure. Credit where due: this is a clean, honest numerical study. The resonance conditions in Eq. 11 do match the dips in Fig. 3, the thermal noise analysis is sensible, and the parameters are chosen from real experiments. This is a legitimate extension of magnon blockade from the resonant regime (Ref 42) to the dispersive regime, and the mapping of the anharmonic level structure is useful.\n\nThe soft spots are real and one is load-bearing. The paper uses H_disp = χ(m†m + 1/2)σ_z, derived from the Jaynes-Cummings interaction under Δ = |ω_q − ω_m| ≫ g. But the text never reports g or Δ, only χ = g²/Δ. From that relation, the plotted χ/γ = 45 with Δ/g = 10 requires g/γ = 450, far above anything in the cited experiments (g/γ ~ 10). If Δ/g is only 3, the neglected higher-order terms are (g/Δ)² ~ 11% of χ, which is larger than the linewidths. Since the blockade minima sit exactly at the dispersive level shifts, a few-percent error in the Hamiltonian can move or wash them out. The paper needs to state g and Δ, justify the validity regime for every plotted χ, and preferably verify the blockade in the full qubit-magnon model without adiabatic elimination.\n\nTwo smaller issues: there are no Fock-space truncation or convergence checks, and the claim of 'experimentally viable' based on 50 mK base temperature glosses over the fact that χ/γ = 20–45 has not been demonstrated; the experiments cited operate at smaller dispersive shifts. These are fixable.\n\nWho is this for? Theorists working on quantum magnonics and circuit QED. It is not a transformative result, but it is a solid extension with a clear gap in parameter validation. I would send it to a serious referee, but the referee should ask for the missing validity check and convergence data before acceptance.","headline":"Good numerical extension of magnon blockade to the dispersive regime, but the paper never checks whether the effective Hamiltonian is valid at the plotted couplings; the blockade may be an artifact.","tokens_in":12826,"tokens_out":2528,"would_cite":false,"duration_ms":23325,"reading_group":"maybe","serious_thinker":"yes","would_accept_peer_review":true},"rs_alignment":null,"lean_confirmation":null,"pith_extraction":{"msc":[],"pacs":[],"model":"deepseek-v4-flash","headline":"Dispersive qubit coupling drives magnon blockade to 0.04.","keywords":["magnon blockade","dispersive regime","superconducting qubit","yttrium-iron-garnet","second-order correlation function","single-magnon source","qubit-induced anharmonicity","quantum magnonics"],"falsifier":"Take the parameters the paper uses, recover $g_{qm}$ and $\\Delta_{qm}$ from $\\chi_{qm}=g_{qm}^2/\\Delta_{qm}$, and solve the same master equation with the full linear interaction $g_{qm}(m^{\\dagger}\\sigma_- + m\\sigma_+)$ instead of the dispersive Hamiltonian; if $g^{(2)}(0)$ does not dip near 0.04 under those drive conditions, the blockade is an artifact of the effective model.","tokens_in":11633,"feed_emoji":"🧲","tokens_out":8648,"duration_ms":72516,"temperature":0.7,"pith_summary":"The paper claims that a magnon blockade—suppression of simultaneous two-magnon events—can be produced by a dispersive, not resonant, interaction between a superconducting qubit and the Kittel mode of a yttrium-iron-garnet sphere. Solving the Lindblad master equation for the driven system, the authors find that the second-order correlation drops to $g^{(2)}(0)\\approx0.04$ when the dispersive coupling is strong and the drive sits on a single-magnon resonance. The mechanism is qubit-induced anharmonicity: each magnon shifts the qubit transition by $2\\chi_{qm}$, pushing two-magnon states off resonance. This matters because it extends single-magnon control to the dispersive regime, where qubit readout is straightforward, and the authors show the blockade survives at dilution-refrigerator temperatures.","feed_headline":"Dispersive qubit coupling drives magnon blockade to 0.04","feed_subtitle":"In a YIG sphere coupled to a superconducting qubit, two-magnon events nearly vanish, enabling single-magnon sources.","key_machinery":"The machine is the dispersive interaction $H_{\\rm disp}=\\frac{1}{2}[2\\chi_{qm}(m^{\\dagger}m+\\tfrac{1}{2})]\\sigma_z$, a qubit-state-dependent shift of the magnon ladder that grows with magnon number. This shift makes the $|g,2\\rangle$ and $|e,2\\rangle$ levels leave the harmonic ladder, forbidding two-magnon transitions from the driven low states; the paper names this qubit-induced anharmonicity. The observable that carries the argument is the equal-time second-order correlation $g^{(2)}(0)=\\langle m^{\\dagger}m^{\\dagger}mm\\rangle/\\langle m^{\\dagger}m\\rangle^2$, with $g^{(2)}(0)<1$ marking antibunching, and the companion magnon-number probabilities $P_1$ and $P_2$, which show $P_2\\ll P_1$ at every blockade point.","core_discovery":"The central claim is that strong dispersive coupling alone turns a magnon mode into a single-magnon emitter. In the rotating frame, the system is governed by the effective Hamiltonian $H' = \\Delta_m m^{\\dagger}m + \\frac{1}{2}\\Delta_q\\sigma_z + \\frac{1}{2}(2\\chi_{qm}m^{\\dagger}m)\\sigma_z + \\Omega_s(\\sigma_+ + \\sigma_-) + \\Omega_d(m^{\\dagger}+m)$, together with Lindblad dissipation for the qubit and magnon baths. Numerically solving the resulting master equation, the paper reports that at $\\chi_{qm}/\\gamma=20$ and $40$—and up to $45$—there is a window of driving detuning in which $g^{(2)}(0)\\to0.04$, with the two-magnon probability $P_2$ far below the single-magnon probability $P_1$; outside these windows the same parameters give bunching with $g^{(2)}(0)\\to100$. The blockade minima occur at single-magnon resonances, while the bunching peaks occur at two-magnon resonances, matching the level-shift picture quantitatively.","pith_inferences":["One direct test would be to fix the bare qubit-magnon coupling and sweep the detuning, checking whether the blockade minimum tracks the predicted single-magnon resonance; this would confirm the shift mechanism rather than only the fit.","Because the mechanism is generic level anharmonicity, the same dispersive blockade should appear for other bosonic modes—phonons or microwave photons—when a qubit is dispersively coupled to them, although the required coupling range will differ.","The reported optimum $g^{(2)}(0)\\approx0.04$ is not shown to be a fundamental lower bound; the paper leaves open whether still larger $\\chi_{qm}/\\gamma$ pushes the correlation closer to zero or whether dissipation sets a floor.","Existing single-shot dispersive magnon detection suggests that the magnon-number distribution itself, not just the correlation function, could be measured directly in the regime studied here."],"forward_implications":["A dispersively coupled YIG-qubit device can serve as a tunable single-magnon source, with the blockade switched on and off by choosing the drive detuning.","Because the effect appears in the already-demonstrated strong dispersive regime, it can be combined with dispersive qubit readout and single-shot magnon detection.","At dilution-refrigerator base temperatures near 46 to 48 mK the blockade remains observable; the dominant threat is thermal population of the magnon mode, and the paper quantifies the destruction threshold near $m_{\\rm th}\\sim0.0035$.","The same level-shift mechanism gives a controlled transition from antibunching ($g^{(2)}(0)\\approx0.04$) to strong bunching ($g^{(2)}(0)\\approx100$) as the drive is swept across single- and two-magnon resonances."],"supporting_citations":[{"why":"Supplies the dispersive-regime magnon-qubit platform and the dissipation-based single-magnon sensing parameters used in the numerical simulations.","marker":"[27]"},{"why":"Provides the experimental demonstration of single-magnon detection in the strong dispersive regime and the cryogenic operating conditions quoted for viability.","marker":"[28]"},{"why":"The prior resonant-regime magnon blockade result that this paper extends to the dispersive regime and uses as the comparison for blockade physics.","marker":"[42]"},{"why":"Derives the Jaynes-Cummings interaction for superconducting circuits, from which the dispersive Hamiltonian of Eq. (4) is obtained.","marker":"[57]"},{"why":"Derives qubit-photon dispersive shifts and number splitting, the physics used here to justify the qubit-induced anharmonicity.","marker":"[58]"},{"why":"Supplies the Lindblad master equation and the definition and interpretation of the second-order coherence function used to quantify blockade.","marker":"[62]"},{"why":"Provides the numerical master-equation solver used to obtain the steady state and hence the reported $g^{(2)}(0)$ values.","marker":"[63]"}],"fun_headline_variants":["Dispersive coupling squeezes magnon blockade to 0.04","Strong qubit coupling enables single-magnon emission","Magnon blockade via dispersive qubit: g2 = 0.04","Superconducting qubit suppresses two-magnon events","Dispersive magnon blockade: single magnons at 0.04"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The prediction rests on the dispersive Hamiltonian $H_{\\rm disp}$ remaining valid at the strong couplings used ($\\chi_{qm}/\\gamma$ up to 45), but the paper does not give the underlying coupling $g_{qm}$ and detuning $\\Delta_{qm}$ needed to check the derivation condition $\\Delta_{qm}\\gg g_{qm}$.","fun_headline_variants_meta":{"raw":{"variants":["Dispersive coupling squeezes magnon blockade to 0.04","Strong qubit coupling enables single-magnon emission","Magnon blockade via dispersive qubit: g2 = 0.04","Superconducting qubit suppresses two-magnon events","Dispersive magnon blockade: single magnons at 0.04"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000876,"raw_usage":{"total_tokens":3789,"prompt_tokens":948,"completion_tokens":2841,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":564,"completion_tokens_details":{"reasoning_tokens":2749}},"tokens_in":564,"tokens_out":2841,"duration_ms":19317,"temperature":1.0,"reasoning_tokens":2749,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-16T11:23:43.674616+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Take the parameters the paper uses, recover $g_{qm}$ and $\\Delta_{qm}$ from $\\chi_{qm}=g_{qm}^2/\\Delta_{qm}$, and solve the same master equation with the full linear interaction $g_{qm}(m^{\\dagger}\\sigma_- + m\\sigma_+)$ instead of the dispersive Hamiltonian; if $g^{(2)}(0)$ does not dip near 0.04 under those drive conditions, the blockade is an artifact of the effective model.","supporting_citations":[{"cited_title":"Blais, R.-S, Huang, A","cited_arxiv_id":null,"evidence_quote":"Derives the Jaynes-Cummings interaction for superconducting circuits, from which the dispersive Hamiltonian of Eq. (4) is obtained."},{"cited_title":"Gambetta, A","cited_arxiv_id":null,"evidence_quote":"Derives qubit-photon dispersive shifts and number splitting, the physics used here to justify the qubit-induced anharmonicity."},{"cited_title":"Walls and G","cited_arxiv_id":null,"evidence_quote":"Supplies the Lindblad master equation and the definition and interpretation of the second-order coherence function used to quantify blockade."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Provides the numerical master-equation solver used to obtain the steady state and hence the reported $g^{(2)}(0)$ values."}],"review_version":1}