{"id":"67d6daa6-ef43-4504-89b5-443f95587c46","arxiv_id":"2608.07034","paper_version":1,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":5.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":5,"one_line_summary":"Adding more phonon modes between cavity-magnon polaritons monotonically boosts their stationary entanglement, while phase synchronization falls and amplitude synchronization turns into collective squeezing.","lead":"This paper proposes that coupling two cavity-magnon polariton modes to multiple vibrational (phonon) modes increases their steady-state entanglement as the number of phonon channels grows. It also reports that the same phonon scattering produces amplitude squeezing while degrading phase synchronization, which may help design quantum-correlation control in hybrid magnomechanical devices.","discovery_kind":"extension","skeptic_critique":{"model":"deepseek-v4-flash","headline":"The scaling with k assumes all phonon modes have equal bare coupling and near-resonant frequencies; this is stated, not justified, and no k=1 baseline is given.","rationale":"The paper's main new claim is the monotonic enhancement of polariton entanglement with the number of phonon modes k, together with the quadrature-selective synchronization tradeoff. Both results are computed under the explicit assumption that all k modes have comparable bare magnetostrictive couplings G0j and are placed near resonance by a small frequency spacing Δω/2π = 2 MHz. The reader identified exactly this assumption as the weakest point, and I agree: it is the single load-bearing condition for the claimed scaling. The drift-matrix construction and Lyapunov solution are internally consistent; I checked the signs and block structure of the coupling matrices and found no algebraic error. The synchronization measures are used correctly: Sc ≤ 1 follows from the Heisenberg uncertainty relation for the difference quadratures, and Sx > 1 correctly signals collective amplitude squeezing. What is missing is any demonstration that the equal-coupling, near-degenerate phonon ensemble is physically realizable or that the result survives realistic disorder. Since the paper is careful to phrase the claim as 'in the parameter regime explored here,' the appropriate verdict remains CONDITIONAL: the numerics are sound, but the headline scaling rests on a stated but unverified assumption. A single numerical experiment—varying G0j across modes—would settle whether the assumption is load-bearing. The reader's CONDITIONAL verdict already reflects this, so no verdict change is needed.","tokens_in":20188,"tokens_out":9394,"duration_ms":774418,"concrete_test":"Recompute Fig. 4(b) (Emax_N vs k) with the same parameters but replace G0j = G0 by a realistic disordered set, e.g., draw G0j from a log-normal distribution with 50% spread or set G0j = G0/(1+(j-1)/10), while keeping ωj = ω1+(j-1)Δω and all other parameters fixed. If Emax_N(k) is no longer strictly increasing (e.g., k=50 value falls below k=25), the identical-mode assumption is load-bearing and the central claim needs to be qualified to the fine-tuned equal-coupling regime. As a complementary check, add the k=1 point to Fig. 4(b) under identical parameters to establish the single-mode baseline.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The central scaling result—Figs. 2(d), 4(b), and 8—rests on treating all k phonon modes as identical parallel channels: same bare coupling G0j = G0, same damping κj, and frequencies ωj within a few MHz of the polariton resonance (Δω/2π = 2 MHz). The paper states this assumption explicitly in Sec. III (discussion of Fig. 5): 'we assume comparable bare magnetostrictive couplings G0j ≃ G0 for the selected phonon modes.' This is not a harmless normalization: with Gj ∝ G0j⟨c⟩, the collective phonon-mediated rate entering the Lyapunov equation is a sum over modes weighted by G_j^2 (and detunings), so a realistic spread in G0j—which YIG eigenmodes inevitably have—changes the effective mode count. If G0j decrease with mode index or fluctuate, the marginal contribution of each added mode is not equal, and the monotonic EN(k) enhancement shown in Fig. 4(b) may saturate or reverse. The paper provides no evidence that k=50 modes with comparable G0j and near-degenerate frequencies can be selected in a 250-μm YIG sphere; the cited experiments (Refs. [10,29]) demonstrate strong coupling to only a few phonon modes. A k=1 baseline is also absent, so the claimed 'beyond single-mode schemes' is inferred from k=2 rather than demonstrated.","agreement_with_reader":"agree"},"referee_report":{"model":"deepseek-v4-flash","summary":"The paper studies a cavity magnomechanical system in which two cavity-magnon polaritons are coupled to k phonon modes, and analyzes their stationary quantum correlations using linearized quantum Langevin equations, the Lyapunov equation for the covariance matrix, logarithmic negativity for bipartite entanglement, residual contangle for tripartite entanglement, and the Mari et al. synchronization measures. The central claim is that, in the parameter regime explored, increasing k monotonically enhances the stationary polariton-polariton entanglement (from about 0.14 at k=2 to 0.33 at k=50), while the amplitude synchronization Sx exceeds unity and grows with k, and the phase synchronization Sp decreases with k. The paper also reports mode-resolved tripartite entanglement and a robustness analysis with respect to temperature and damping.","tokens_in":20544,"tokens_out":5556,"duration_ms":54587,"significance":"If the central scaling result is correct, the paper would establish multimode phonon mediation as a practical route to stronger polariton entanglement and quadrature-selective synchronization, going beyond single-phonon protocols. The manuscript uses standard and internally consistent methods: the linearization, Lyapunov solution, stability checks, and Gaussian entanglement measures are appropriate, and the parameter table is explicit and experimentally anchored. The paper also contains useful numerical results on mode-resolved tripartite entanglement and the trade-off between entanglement and phase synchronization. However, the decisive claim of monotonic k-enhancement currently rests on an idealized identical-channel assumption and on the absence of a k=1 baseline, so the significance is conditional on a robustness analysis that is not yet provided.","major_comments":[{"comment":"The monotonic enhancement of EN with k is computed under the explicit assumption that all selected phonon modes have comparable bare magnetostrictive couplings (G0j ≃ G0), the same damping κj, and near-resonant frequencies spaced by Δω/2π = 2 MHz. Since the coupling blocks in the drift matrix, Eqs. (9)-(12), enter as sums over modes weighted by Gj, a realistic spread in G0j or in the phonon frequencies changes the effective number of participating channels. The manuscript provides no evidence that a 250-μm YIG sphere offers many modes with comparable magnetostrictive strength within a few MHz of the polariton resonance, and the cited experiments demonstrate strong coupling only to a few modes. I ask the authors to either justify the mode-selection assumption experimentally or add a systematic robustness study with disordered G0j, κj, and ωj, showing that the monotonic EN(k) behavior survives realistic spreads. The closing sentence of the Conclusion, which lists \"disorder in phonon spectra\" as future work, indicates that this missing test is recognized by the authors.","section":"§III, discussion of Fig. 5(b), underpinning Figs. 2, 4, 6, 8"},{"comment":"No k=1 baseline is provided. The abstract and introduction claim that the multimode scheme goes \"beyond\" conventional single-mode protocols, but the comparison is made only against k=2. Since the claimed monotonic scaling cannot be assessed without the single-channel reference, the authors should compute EN for k=1 in the same setup, or provide the analytic single-mode limit, and show that the k=2 value is not already saturating the accessible entanglement. If the single-mode entanglement is comparable to the k=2 value, the phrase \"monotonic enhancement\" would require reinterpretation.","section":"§III, Fig. 4(b)"},{"comment":"The scaling result is obtained after optimizing θ (and selecting Δω/2π = 2 MHz) for each k, with all modes treated as identical. Because θopt shifts with k and because every added mode has the same effective coupling, it is not clear whether the monotonicity is a physical property of the multimode mechanism or an artifact of re-optimizing a collective parameter under an identical-channel assumption. The authors should either (i) provide an analytic argument for the scaling of the effective collective coupling with k, or (ii) show monotonicity at a fixed, experimentally chosen θ and at detunings away from the optimal resonance condition. Without this, the central claim that \"each additional phonon mode opens a new scattering pathway\" with a monotonic effect is not yet established beyond the specific selected parameters.","section":"§III, Fig. 4(a)-(b) and Fig. 8"}],"minor_comments":[{"comment":"The line \"show that , in the parameter regime\" contains an extra space before the comma; it should read \"show that, in the parameter regime\".","section":"Abstract"},{"comment":"The block matrix notation uses \"0k\" in the off-diagonal phonon blocks; this should be written as 02×2 or replaced by an explicit zero matrix, since each block is two-dimensional.","section":"Eq. (9)"},{"comment":"The table caption reads \"Tab. ( I)\" instead of \"Table I\", and the mechanical damping is denoted κb in Fig. 7 while the main text uses κj; please unify the notation.","section":"Table I caption and Fig. 7"},{"comment":"The difference quadratures are introduced as X̃− and P̃−, but the arbitrary quadrature Xϕ is written without a tilde; define all quadrature variables consistently.","section":"Eq. (28) and surrounding text"},{"comment":"There are several typographical ligature issues such as \"efficient\" for \"efficient\"; a final proofreading pass is needed.","section":"Throughout"}],"recommendation":"major_revision","confidential_remarks":"The manuscript's central claim is interesting and the numerical framework is sound, but the current evidence for the monotonic k-scaling is not yet convincing because it relies on identical near-resonant phonon channels and lacks a k=1 baseline. In my view, this is a load-bearing issue rather than a presentation issue: the abstract and conclusion state the monotonic enhancement as a general result, while the underlying assumption is explicitly local to the example. If the authors can add the requested robustness analysis and single-mode baseline, the revised manuscript would be suitable for reconsideration. I would also ask the editor to ensure that the novelty relative to Ref. [28] is clearly delineated, since the two-mode case is a direct extension of that work."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"Two things to know. First, the core result—stationary polariton–polariton entanglement grows monotonically with phonon mode count k, from EN≈0.14 (k=2) to ≈0.33 (k=50)—is real in the model, but it is not a robust prediction about any specific experiment. It rests on assuming all k phonon modes are identical, equally coupled, near-resonant channels. The paper states this plainly, but it is load-bearing for the headline. Real YIG spheres have discrete mechanical modes with different frequencies, damping, and magnetostrictive strengths; no evidence is given that 50 modes with comparable G0j can be selected. Second, the genuinely new content is the trade-off between entanglement and synchronization: increasing k strengthens amplitude squeezing (Sx>1) while degrading phase synchronization. That quadrature-selective behavior is a nice, non-obvious observation.\n\nWhat the paper does well: the QLE linearization, Lyapunov solution, and CV entanglement/synchronization measures are standard and executed carefully, with stability checks for every parameter set. The robustness scans over temperature, damping, and detuning are genuinely useful. The mode-resolved tripartite entanglement—each phonon mode entangles most efficiently near its own resonance—is a multimode feature absent in single-phonon schemes (Ref. [28]).\n\nSoft spots, in order. (1) No k=1 baseline. The claim of going beyond single-mode schemes is inferred from k=2, never computed. This is trivial to fix and should be required. (2) The identical-mode assumption is stated, not justified. The authors should offer a physical selection argument or run a sensitivity scan with a spread in G0j. If monotonic EN(k) survives realistic disorder, the claim becomes much stronger. (3) The k-scaling is purely numerical; a simple analytic estimate of how the collective coupling enters EN would add real trust. (4) Minor: a stray comma in the abstract (\"we show that ,\") and occasional loose wording, but nothing damaging.\n\nThe stress-test concern lands. Still, the central argument holds up as a theoretical proposal under its stated assumptions. This paper deserves a serious referee, not a desk reject. Ask for the k=1 baseline and a disorder scan; the core mechanism and the synchronization trade-off are defensible.","headline":"Standard-formalism paper whose k-scaling result is genuine but rests on an idealized identical-mode assumption; worth refereeing with clear requests for a k=1 baseline and a robustness scan over mode couplings.","tokens_in":21074,"tokens_out":4324,"would_cite":true,"duration_ms":38687,"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":"More phonon modes monotonically strengthen stationary entanglement between cavity-magnon polaritons, while phase synchronization weakens and amplitude squeezing grows.","keywords":["cavity magnomechanics","polariton entanglement","multimode phonon mediation","quantum synchronization","Stokes and anti-Stokes scattering","logarithmic negativity","continuous-variable Gaussian states","YIG sphere"],"falsifier":"Take a real YIG sphere with its measured discrete mechanical eigenfrequencies and bare magnetostrictive couplings, and compute the steady-state logarithmic negativity as more phonon modes are coupled while keeping the drive fixed; if $E_N$ stops growing monotonically with $k$ once the comparable-coupling assumption is replaced by realistic uneven couplings, the central claim is falsified.","tokens_in":19991,"feed_emoji":"🧲","tokens_out":6357,"duration_ms":54746,"temperature":0.7,"pith_summary":"This paper claims that in a cavity magnomechanical system, coupling the two cavity-magnon polaritons to many vibrational (phonon) modes turns each mode into a parallel scattering channel, and the stationary entanglement between the polaritons grows monotonically with the number of channels $k$: the logarithmic negativity rises from about $0.14$ at $k=2$ to about $0.33$ at $k=50$. The same Stokes and anti-Stokes scattering processes that build entanglement also produce quadrature-selective synchronization: amplitude synchronization rises above the vacuum level ($S_x>1$, growing with $k$), while phase synchronization decreases with $k$. If correct, this gives a quantitative route to stronger bipartite entanglement and collective squeezing in hybrid magnomechanical platforms without external squeezing drives. The paper also shows that genuine tripartite entanglement between the two polaritons and each phonon mode peaks at that mode's own frequency, so the multimode configuration acts as a frequency-selective correlation engine.","feed_headline":"Adding phonon modes doubles polariton entanglement","feed_subtitle":"In a magnomechanical cavity, every extra vibration channel strengthens quantum correlations but weakens phase locking.","key_machinery":"The load-bearing object is the multimode drift matrix $\\Lambda$ of the linearized quantum Langevin equations, obtained by expanding about large steady-state polariton amplitudes produced by a strong magnon drive. Each phonon mode $b_j$ contributes a pair of $2\\times 2$ coupling blocks $A_{j\\pm}$ and $B_{j\\pm}$ that connect the polariton quadratures to the mechanical quadratures with effective couplings $G_j = 2i G_{0j}\\langle c\\rangle$, so when the bare couplings are comparable ($G_{0j}\\simeq G_0$) all channels share a common strength set by the drive. Solving the Lyapunov equation $\\Lambda V + V\\Lambda^T = -\\Gamma$ for the $(4+2k)\\times(4+2k)$ covariance matrix yields the entanglement and synchronization quantifiers; the resonance condition $\\Delta_+ \\simeq \\omega_1$, $\\Delta_- \\simeq -\\omega_1$ places the lower polariton on the Stokes sideband and the upper polariton on the anti-Stokes sideband, which is the physical route by which correlations are generated and transferred.","core_discovery":"The central claim is that multimode phonon mediation is not merely a small correction to the single-phonon protocol: in the parameter regime explored, each additional phonon mode adds a near-resonant Stokes and anti-Stokes scattering pathway, and the correlations transferred along these parallel pathways accumulate. As a result, the steady-state logarithmic negativity $E_N$ between the upper and lower polariton modes increases monotonically with the number of phonon modes $k$, from $E_N\\simeq 0.14$ at $k=2$ to $E_N\\simeq 0.33$ at $k=50$. The same mechanism drives quantum synchronization between the polaritons but with opposite scaling in the two quadratures: the phase-synchronization measure $S_p$ falls as $k$ grows, while the amplitude measure $S_x$ exceeds unity and rises from about $1.15$ to $1.30$, meaning the polariton difference mode is collectively squeezed in amplitude without any injected squeezing. The paper further reports that genuine tripartite entanglement among the polaritons and a given phonon mode is activated only above a $k$-dependent frequency-spacing threshold and peaks when the polariton splitting matches that mode's frequency.","pith_inferences":["If the monotonic scaling persists under realistic mode spectra, YIG spheres or phononic crystals with many near-degenerate mechanical modes could be used to increase polariton entanglement without raising the drive power.","The reported tradeoff suggests a control knob: selecting $k$ tunes the quadrature in which synchronization is strongest, which could be used to route quantum correlations in a frequency-selective way.","The mode-resolved tripartite entanglement indicates that the same platform could act as a multimode quantum interface where each mechanical mode is individually addressable by detuning.","A direct experimental test would be to sweep the number of actively coupled phonon modes and check whether both the entanglement curve and the reduction in $S_p$ match the predicted scaling."],"forward_implications":["Stationary polariton-polariton entanglement increases with phonon number $k$, reaching $E_N\\simeq 0.33$ at $k=50$, more than twice the $k=2$ value.","The optimal hybridization angle shifts toward the instability boundary as $k$ grows, so the multimode ensemble relaxes the asymmetry between Stokes and anti-Stokes channels required for stable entanglement.","Amplitude synchronization $S_x$ exceeds unity and grows with $k$, demonstrating collective squeezing of the polariton difference mode without an external squeezing drive.","Phase synchronization $S_p$ decreases as $k$ grows, so entanglement enhancement and phase synchronization compete in the multimode regime.","Each phonon mode mediates genuine tripartite entanglement peaked at its own frequency, enabling mode-resolved, frequency-selective correlation engineering."],"supporting_citations":[{"why":"This single-phonon protocol is the baseline that the paper extends to $k$ parallel scattering channels.","marker":"[28]"},{"why":"It supplies the cavity magnomechanics model and the experimentally feasible parameter values used throughout the paper.","marker":"[10]"},{"why":"It establishes the drive-enhanced effective magnomechanical coupling and the entanglement-generation approach in cavity magnomechanics.","marker":"[12]"},{"why":"It defines the continuous-variable synchronization measures $S_c$, $S_p$, and $S_x$ used to quantify quadrature-selective synchronization.","marker":"[19]"},{"why":"It provides the Lyapunov-equation method for computing steady-state covariance matrices of linearized Gaussian systems.","marker":"[37]"},{"why":"It demonstrates strong coupling of cavity-magnon polaritons to phonons, grounding the polaromechanical platform in experiment.","marker":"[29]"},{"why":"It documents the experimentally accessible range of cavity-magnon coupling strengths needed to re-optimize the mixing angle as $k$ grows.","marker":"[47]"},{"why":"It establishes strong cavity-magnon coupling and the polariton hybridization used as the base platform.","marker":"[5]"}],"fun_headline_variants":["More phonon modes boost polariton entanglement","Extra phonon channels multiply polariton entanglement","More phonons: stronger entanglement, weaker phase sync","Multimode phonon mediation entangles more, syncs less","Phonon-rich cavities entangle polaritons, trade sync"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The monotonic scaling assumes that every added phonon mode has comparable bare magnetostrictive coupling and sits close enough to the polariton resonance that the microwave drive uniformly raises all effective couplings to one common value.","fun_headline_variants_meta":{"raw":{"variants":["More phonon modes boost polariton entanglement","Extra phonon channels multiply polariton entanglement","More phonons: stronger entanglement, weaker phase sync","Multimode phonon mediation entangles more, syncs less","Phonon-rich cavities entangle polaritons, trade sync"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000971,"raw_usage":{"total_tokens":4166,"prompt_tokens":1017,"completion_tokens":3149,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":633,"completion_tokens_details":{"reasoning_tokens":3071}},"tokens_in":633,"tokens_out":3149,"duration_ms":20935,"temperature":1.0,"reasoning_tokens":3071,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-10T16:00:43.288411+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Take a real YIG sphere with its measured discrete mechanical eigenfrequencies and bare magnetostrictive couplings, and compute the steady-state logarithmic negativity as more phonon modes are coupled while keeping the drive fixed; if $E_N$ stops growing monotonically with $k$ once the comparable-coupling assumption is replaced by realistic uneven couplings, the central claim is falsified.","supporting_citations":[{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"This single-phonon protocol is the baseline that the paper extends to $k$ parallel scattering channels."},{"cited_title":"Zuo, Z.-Y","cited_arxiv_id":null,"evidence_quote":"It supplies the cavity magnomechanics model and the experimentally feasible parameter values used throughout the paper."},{"cited_title":"Li, S.-Y","cited_arxiv_id":null,"evidence_quote":"It establishes the drive-enhanced effective magnomechanical coupling and the entanglement-generation approach in cavity magnomechanics."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"It defines the continuous-variable synchronization measures $S_c$, $S_p$, and $S_x$ used to quantify quadrature-selective synchronization."},{"cited_title":"Y., Qian, H., Shen, R","cited_arxiv_id":null,"evidence_quote":"It demonstrates strong coupling of cavity-magnon polaritons to phonons, grounding the polaromechanical platform in experiment."},{"cited_title":"D., Jusserand, B., Perrin, B, Strong Optical-Mechanical Coupling in a Vertical GaAs/AlAs Microcavity for Subterahertz Phonons and Near-Infrared Light","cited_arxiv_id":null,"evidence_quote":"It documents the experimentally accessible range of cavity-magnon coupling strengths needed to re-optimize the mixing angle as $k$ grows."},{"cited_title":"Zhang, C.-L","cited_arxiv_id":null,"evidence_quote":"It establishes strong cavity-magnon coupling and the polariton hybridization used as the base platform."}],"review_version":1}