{"id":"e69691c7-9d37-479b-90e0-b6274d49ae43","arxiv_id":"2608.05312","paper_version":1,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":7.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":7,"one_line_summary":"Off-diagonal exciton-phonon coupling creates a one-way, system-size-independent drain from dark states into the bright cavity mode, restoring near-unity transport.","lead":"This paper shows that vibrations of emitters inside an optical cavity can open a one-way door that releases energy trapped in 'dark states' and sends it back into the light channel. If the mechanism works in real materials, it could make cavity-based energy transport nearly lossless and give experimenters four clear signatures to test it.","discovery_kind":"new_application","skeptic_critique":{"model":"deepseek-v4-flash","headline":"Size-independence claim rests on flat non-Condon spectral density over polariton splitting scale; SM Eq. (S8) is marginal at N=6 and worsens with N, so the sum-rule escape rate need not be N-independent for realistic structured baths.","rationale":"Decomposing the central claim helps locate the fragile step. The assertion that transitions into dark states vanish identically is a structural property of the Redfield tensor (SM Eq. S6): every transition carries the target-state cavity weight w_k, and w_k=0 for dark states, so this part is robust for any spectral density at T=0. The second assertion, that the total escape rate from any dark state equals gamma_rec and is independent of N, is what makes the mechanism qualitatively different from dephasing. That assertion does not follow from the photonic-weight sum rule alone; it requires the replacement J(E_l - E_k) -> J(Delta) so that sum_k J(E_l-E_k) w_k factors as J(Delta) sum_k w_k = J(Delta). Once J varies on the scale sqrt(N) g, the escape rate is a weighted average of J over the bright-state manifold and depends on the dark state's energy and on N. The paper's SM Eq. (S8) states the needed condition and concedes it is marginal at N=6; at N=96 the polariton splitting is ~14.7 meV, comparable to or larger than typical phonon cutoffs in molecular systems, so the approximation cannot be assumed. The numerical demonstrations use the site-basis Lindblad model, which already contains the flat-J reduction; hence they confirm the reduced model's behavior, not the validity of the reduction. I therefore agree with the reader that this is the load-bearing approximation. I would not move the verdict: the mechanism is physical and the unidirectionality is exact in the clean T=0 limit, but the headline size-independence should be stated as conditional on a spectrally flat non-Condon bath, or demonstrated with the full Redfield tensor for realistic J. The disorder caveat (w_k small but nonzero) is secondary and, given the small delta t used, does not change this conclusion. Verdict remains CONDITIONAL.","tokens_in":17501,"tokens_out":9223,"duration_ms":75485,"concrete_test":"Recompute the dark-manifold dynamics using the full eigenbasis Redfield tensor of SM Eq. (S6) with a structured non-Condon spectral density, e.g., super-Ohmic J(omega)=lambda omega^3 exp(-omega/omega_c) with omega_c = 5 meV and the main-text parameters (g=1.5 meV, delta t=0.5 meV, gamma_lead=0.5 meV, T=0), for N=6, 32, 96. Compare the total dark decay rate (slope of ln p_D(t)) and peak efficiency eta* to the flat-J site-basis results in Fig. 1(c). If the decay rate or eta* shifts by more than ~10% at any N, the headline size independence fails.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The central size-independence result, W_{k<-l}=gamma_rec w_k(1-w_l) (Eq. 4), is exact only after the full Redfield tensor (SM Eq. S6) is reduced to a site-basis Lindblad form by assuming J_i(omega_lk) ≈ J_i(Delta) for all transitions. The microscopic rate from a dark state l to a bright state k is proportional to J(E_l - E_k) w_k |<i|psi_l>|^2; summing over k and i yields sum_k J(E_l - E_k) w_k (1 - w_l), which depends on the bright-state energy distribution and on N via the polariton splitting ~ sqrt(N) g. The paper's own validity condition, J'(Delta) sqrt(N) g / J(Delta) << 1 (SM Eq. S8), is marginal for N=6 (cutoff ~5 meV) and is increasingly violated at N=96 (sqrt(96)*1.5 ≈ 14.7 meV) for any J with a peak or cutoff on that scale. The SM even concedes O(30%) modifications to gamma_rec in the worst case at N=6. Thus the quantitative headline—total dark escape rate exactly gamma_rec, independent of N—is not established for structured phonon baths; the unidirectionality itself (transitions into w=0 dark states vanish) does survive for any J, but the size-independent rate is an artifact of the flat approximation. The numerics in Figs. 1(c) and 2 are generated with the site-basis Lindblad model, i.e., they assume the reduction under test, so they cannot validate it. A secondary issue is that static disorder (delta t=0.5 meV) gives dark states small but nonzero w_k, making the strict one-way property approximate in the simulated system; the paper acknowledges this, but it means the numerics test near-unidirectionality, not exact unidirectionality.","agreement_with_reader":"agree"},"referee_report":{"model":"deepseek-v4-flash","summary":"The paper considers N two-level emitters coupled to a single cavity mode and to a phonon bath, and identifies a non-Condon (off-diagonal) exciton-phonon coupling that modulates the light-matter coupling. In the zero-temperature, single-excitation limit this yields a jump operator L_rec,i = sqrt(gamma_rec)|cav><i|, whose eigenbasis transition rates are W_{k<-l} = gamma_rec w_k (1 - w_l). Because dark states have zero photonic weight in the clean uniform-coupling limit, transitions into dark states vanish and the total escape rate out of any dark state equals gamma_rec, independent of N. The authors support this claim with a Born-Markov/Redfield derivation in the Supplemental Material, closed-form rate equations, direct Lindblad propagation in cavity-drain and site-drain geometries, a QCLE benchmark, and four proposed experimental signatures: single-exponential dark decay, a size-scaling efficiency gap, a detailed-balance temperature crossover, and resonant enhancement of the rescue rate at vibrational modes.","tokens_in":17936,"tokens_out":6188,"duration_ms":56686,"significance":"If the central reduction is valid, the paper identifies a genuinely new mechanism: a unidirectional dark-to-bright valve protected by a photonic-weight sum rule, structurally distinct from the bidirectional dark-bright exchange of standard ENAQT. The derivation of Eq. (4) from the site-basis Lindblad model is clean, the photonic-weight sum rule is exact in the clean limit, and the paper offers concrete falsifiable predictions with clear experimental discriminators. The authors also provide reproducible numerical infrastructure, disorder-averaged Lindblad simulations, closed-form two-exponential fits, and an independent QCLE consistency check. The main caveat is that the headline size-independence rests on a flat-spectral-density reduction whose validity is marginal even at the parameters used in the main text; because the numerics assume that reduction, the size-scaling plots do not by themselves establish the claim for structured phonon baths.","major_comments":[{"comment":"The headline result that the total dark escape rate is exactly gamma_rec and independent of N is obtained only after replacing the Redfield rates J_i(omega_lk) by J_i(Delta) in the eigenbasis Redfield tensor. Without that replacement, the rate from a dark state l to a bright state k carries J(E_l-E_k), so the summed escape rate depends on the bright-state energy distribution and, through the polariton splitting sqrt(N) g, on N. The paper's own validity condition, SM Eq. (S8), is described as marginal at N=6 (cutoff ~5 meV) and becomes harder to satisfy as N grows to 96, where sqrt(96) g is about 14.7 meV. The numerics in Figs. 1(c) and 2 are generated with the site-basis Lindblad model, i.e., they assume the very reduction whose validity is at issue; they therefore cannot validate the size independence. The claim should either be restricted to the flat-spectral-density regime with a quantitative error estimate, or tested directly for a structured J(omega) using the full Redfield tensor.","section":"SM Eqs. (S6)-(S8); main-text Eq. (4), Figs. 1(c) and 2"},{"comment":"The statement that transitions into dark states vanish identically relies on w_k=0 for all dark states. The numerical simulations, however, include static hopping disorder delta_t=0.5 meV, which the paper itself notes lifts dark states to small but finite w_k. The 'strictly unidirectional' property and the exact photonic-weight sum rule are therefore clean-limit statements, not properties of the disordered system actually simulated. The paper acknowledges this in passing, but the abstract and main text do not qualify the exactness, and the size-scaling curves in Fig. 1(c) are computed in the disordered system. Please state explicitly whether the reported efficiencies use the clean-limit prediction or the disordered eigenstates with nonzero w_k, and quantify how the finite-w_k corrections scale with N and delta_t.","section":"Main-text Eq. (4); SM Eq. (S6); Figs. 1(c), 2 and SM Table I"},{"comment":"There is an inconsistency in the finite-temperature detailed-balance ratio. Main-text Eq. (3) defines gamma_rec as the T=0 emission rate 2J(Delta), and Fig. 3 states gamma_abs = e^{-Delta/k_B T} gamma_rec. However, SM Eq. (S4) gives gamma_abs = 2J(Delta) n(Delta,T) and gamma_em(T) = 2J(Delta)[1+n(Delta,T)], so gamma_abs/gamma_em(T) = n/(1+n) = e^{-Delta/k_B T}. If gamma_rec in Fig. 3 is intended to mean the finite-temperature emission rate gamma_em(T), the notation conflicts with Eq. (3); if it means the T=0 rate, the equality gamma_abs = e^{-Delta/k_B T} gamma_rec is only approximate for k_B T << Delta. The qualitative crossover survives, but the formula should be stated consistently.","section":"Main-text Fig. 3; SM Eq. (S4) and Eq. (S17)"}],"minor_comments":[{"comment":"The caption contains typos: 'relxation' should be 'relaxation' and 'avergae' should be 'average'.","section":"Fig. 1 caption"},{"comment":"The sentence 'In the cavity-drain configuration (Gamma_D=0, F_D=0)' introduces the symbol F_D without defining it; the second condition appears to be a typo.","section":"SM Eq. (S16)"},{"comment":"The statement that for pure dephasing the dark population relaxes to p_D -> p_B is a degenerate-manifold statement; it may benefit from a brief qualification that this holds in the symmetric limit, since the SM later discusses density-of-states-weighted equilibration.","section":"Main-text, rate-equation paragraph after Eq. (5)"},{"comment":"The claim that |lambda^{(o)}/lambda^{(d)}|^2 is generically of order unity in molecular systems is plausible but cited only to spectroscopic observations; a short derivation or explicit molecular estimate would strengthen the parameter-scaling argument for gamma_rec ~ gamma_deph.","section":"SM, 'Microscopic origin' section"}],"recommendation":"major_revision","confidential_remarks":"The central mechanism is interesting and the clean-limit algebra is sound, but the paper's most striking claim, exact size independence of the dark-to-bright escape rate, is currently established only under a spectral-flatness assumption that the numerics do not independently test. I would like to see either a structured-J calculation with the full Redfield tensor or a clear reframing of the result as valid in the flat-bath regime, before publication."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"Short version: the non-Condon rescue channel is real and the unidirectionality is robust, but the size-independence claim is conditional on a flat spectral density that gets harder to justify as N grows. That is the main thing to know.\n\nWhat the paper does well: it takes a known ingredient (non-Condon, off-diagonal exciton-phonon coupling) and applies it to the collective dark manifold, where it opens a strictly one-way dark-to-bright route because transitions into dark states vanish identically (w_k = 0). That unidirectionality holds for any bath at T=0, as the Redfield tensor in the SM shows. The photonic-weight sum rule then gives a total escape rate that is independent of the dark state's identity and of N, provided the flat-bath reduction is valid. The four signatures (single-exponential dark decay, size-scaling efficiency gap, temperature behavior, resonance at bath modes) are well-chosen and experimentally distinguishable from dephasing-based transport. The paper is well-written and honest about the validity condition.\n\nThe soft spots: the size-independence is exact only when J(omega) is flat over the polariton splitting scale sqrt(N) g. The paper's own condition, SM Eq. (S8), is marginal at N=6 (cutoff ~5 meV) and clearly violated at N=96 (scale ~15 meV). For any structured spectral density, the escape rate will depend on the dark state energy and on N, so the headline quantitative claim is not established for realistic phonon baths. The numerics in Figs. 1(c) and 2 use the site-basis Lindblad model, i.e., they assume the reduction under test and therefore cannot validate it. Static disorder (delta_t = 0.5 meV) gives dark states small nonzero photonic weight, making the strict one-way property approximate in the simulated system; the paper acknowledges this, but it is a real caveat. Also, no code or repository link appears in the arXiv version, which hampers independent reproduction.\n\nThese are addressable rather than fatal. The unidirectionality is exact and new; the size-independence is a conditional bonus that needs a stated regime of validity and a quantitative error estimate. I would send this to a serious referee, with the request that the referee press for a clear statement of when the flat-bath reduction holds and a code release. The paper deserves review, not desk rejection.","headline":"Genuine new rescue channel with a clean sum-rule argument, but the headline size-independence rests on a flat-bath approximation that will not survive structured spectral densities at large N.","tokens_in":18492,"tokens_out":2705,"would_cite":true,"duration_ms":24859,"reading_group":"yes","serious_thinker":"yes","would_accept_peer_review":true},"rs_alignment":null,"lean_confirmation":null,"pith_extraction":{"msc":[],"pacs":["42.50.Pq","71.36.+c"],"model":"deepseek-v4-flash","headline":"Off-diagonal exciton–phonon coupling opens a one-way dark-to-bright rescue channel whose rate is independent of the emitter count.","keywords":["polariton transport","dark states","non-Condon coupling","cavity quantum electrodynamics","environment-assisted quantum transport","Lindblad master equation","photonic weight","Herzberg-Teller polaritons"],"falsifier":"A time-resolved photoluminescence or transient-absorption experiment on arrays with $N$ from about 4 to 64 could test the claim: if the dark-state decay time grows with $N$, or the decay is not single-exponential, or the extracted escape rate $\\gamma_{\\mathrm{rec}}$ changes when the detuning is moved across a phonon mode in a way not captured by $2J(\\Delta)$, the unidirectional-valve picture would be contradicted. In particular, measuring the escape rate for two different dark states with the same detuning but different energies within a structured spectral density would reveal whether the sum-rule collapse holds.","tokens_in":17282,"feed_emoji":"⚛️","tokens_out":6841,"duration_ms":50461,"temperature":0.7,"pith_summary":"This paper argues that the off-diagonal (non-Condon) part of the exciton–phonon coupling creates a strictly one-way relaxation channel from the dark manifold into the bright cavity state in a Tavis–Cummings array. Because dark eigenstates carry zero photonic weight, transitions into them vanish, while the photonic-weight sum rule fixes the total dark-to-bright escape rate at $\\gamma_{\\mathrm{rec}} = 2J(\\Delta)$, independent of the number of emitters $N$. If correct, cavity-coupled transport no longer degrades as the dark manifold grows: the dark population decays single-exponentially and the drained efficiency approaches unity at every system size. The paper also identifies four observable signatures that distinguish this valve from dephasing-based environment-assisted transport.","feed_headline":"Dark states escape at a rate that never shrinks with size","feed_subtitle":"Non-Condon phonon coupling drains the dark manifold at near-unity transport efficiency regardless of N.","key_machinery":"The load-bearing object is the non-Condon spectral density $J(\\omega) = (\\pi/2)\\sum_k |\\lambda_k^{(o)}|^2 \\delta(\\omega-\\omega_k)$, which enters the rates $\\gamma^{\\mathrm{em}} = 2J(\\Delta)[1+n_{\\Delta,T}]$ and $\\gamma^{\\mathrm{abs}} = 2J(\\Delta) n_{\\Delta,T}$. At $T=0$ the absorption channel switches off and the dissipator reduces to the single jump operator $\\hat{L}_{\\mathrm{rec},i}=\\sqrt{\\gamma_{\\mathrm{rec}}}\\,|\\mathrm{cav}\\rangle\\langle i|$ with $\\gamma_{\\mathrm{rec}}=2J(\\Delta)$. The argument then runs on the photonic-weight sum rule: transitions into a dark state vanish because $w_k=0$, and the escape rate from any dark state is pinned to $\\gamma_{\\mathrm{rec}}$ because the bright states carry total weight one.","core_discovery":"The central claim is that the non-Condon (transition-dipole-modulating) exciton–phonon coupling, treated in a Born–Markov rotating-wave approximation, produces the jump operator $\\hat{L}_{\\mathrm{rec},i} = \\sqrt{\\gamma_{\\mathrm{rec}}} |\\mathrm{cav}\\rangle\\langle i|$ at zero temperature. In the eigenbasis the secular rate is $W_{k\\leftarrow l} = \\gamma_{\\mathrm{rec}} w_k (1-w_l)$, where $w_k$ is the cavity weight of eigenstate $k$; since $w_k=0$ for all $N-1$ dark states, the reverse transitions are identically zero and the total rate out of any dark state is $\\sum_{k\\in B} \\gamma_{\\mathrm{rec}} w_k = \\gamma_{\\mathrm{rec}}$ by the sum rule. The paper shows this unidirectionality survives in the full Redfield tensor, not just the secular site-basis reduction, and that in the cavity-drain configuration it yields transport efficiency $\\eta \\to 1$ on a timescale set by $\\max(\\gamma_{\\mathrm{rec}}^{-1}, \\gamma_{\\mathrm{lead}}^{-1})$ with no residual $N$ dependence.","pith_inferences":["If the mechanism is generic, it turns the dark manifold from a liability into a resource: engineered spectral densities $J(\\omega)$ with peaks at chosen detunings would allow targeted, size-independent extraction, which could be exploited in light-harvesting and dissipation-based quantum control.","The same unidirectional structure may apply beyond single-excitation transport to finite pumping and multi-excitation manifolds, where the photonic-weight prefactor would still suppress transitions into dark configurations.","A natural experimental extension is to cold-atom cavity-QED arrays where $\\gamma_{\\mathrm{rec}}$, $\\gamma_{\\mathrm{deph}}$, and $N$ are independently tunable, allowing a direct test of the $N$-independent decay before moving to molecular systems.","In structured molecular baths where the spectral density is not flat, the paper's own validity condition suggests the size independence will fail progressively; measuring $J(\\omega)$ independently would let one predict the $N$ at which the rescue rate starts to drift."],"forward_implications":["In cavity-drain transport, the dark manifold is evacuated at rate $\\gamma_{\\mathrm{rec}}$ regardless of $N$, so the transport efficiency reaches near unity for any array size, whereas dephasing-only transport leaves a dark plateau that grows with $N$.","The same photonic-weight argument predicts four signatures: single-exponential dark-state decay, a size-scaling efficiency gap, temperature dependence set by detailed balance with crossover at $k_B T/\\Delta \\sim 1$, and resonant enhancement of $\\gamma_{\\mathrm{rec}}$ when the detuning matches a vibrational bath mode.","At finite temperature the absorption channel reappears with ratio $\\gamma^{\\mathrm{abs}}/\\gamma^{\\mathrm{rec}} = e^{-\\Delta/k_B T}$, refilling the dark manifold and degrading the valve; the crossover temperature for molecular polariton detunings of tens of meV is around 100 K.","The unidirectional property is robust to the secular approximation: the full Redfield tensor has $w_k$ prefactors on every transition, so transitions into dark states vanish identically at $T=0$ for any bath spectral density.","The rescue advantage is specific to the cavity-drain geometry; with a site-$N$ drain, dephasing slightly outperforms rescue at the peak, as the paper's parameter sweep shows."],"supporting_citations":[{"why":"Supplies the hierarchical-equations treatment of Herzberg–Teller/non-Condon effects used to motivate the non-Condon channel.","marker":"[40]"},{"why":"Reports light emission from vibronic polaritons in coupled metalloporphyrin–multimode cavity systems, the experimental platform where non-Condon coupling is active.","marker":"[42]"},{"why":"Provides evidence distinguishing phonon-modulated emission and absorption in porphyrin cavity polariton samples, supporting the resonance signature.","marker":"[43]"},{"why":"Baseline Condon-coupling theory showing dark–bright relaxation in polariton systems, which the paper contrasts with its unidirectional rescue.","marker":"[27]"},{"why":"Bosonic-mapping analysis of radiative pumping versus vibrational relaxation that the paper compares against its $1/N$-independent rescue rate.","marker":"[28]"},{"why":"Defines the ENAQT dephasing-assisted transport framework whose bidirectional dark–bright exchange is the benchmark the rescue channel is designed to beat.","marker":"[22]"},{"why":"Supplemental derivation of the Born–Markov reduction, the Redfield-tensor unidirectionality, and the validity condition Eq. (S8).","marker":"[39]"},{"why":"Standard open-quantum-systems reference used for the operator identity that collapses the Redfield sum to the site-basis Lindbladian.","marker":"[68]"}],"fun_headline_variants":["One-way dark-state escape, rate fixed by photonic weight","Dark states drain at a size-independent rate in cavities","Non-Condon phonon coupling unlocks dark states for transport","Size no longer matters: dark states escape at constant rate","Cavity transport: dark states drain at fixed rate"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The size-independent rescue rate assumes the non-Condon spectral density is nearly flat on the scale of the polariton splittings, so the site-basis Lindblad collapse is valid, and that dark states have essentially zero photonic weight.","fun_headline_variants_meta":{"raw":{"variants":["One-way dark-state escape, rate fixed by photonic weight","Dark states drain at a size-independent rate in cavities","Non-Condon phonon coupling unlocks dark states for transport","Size no longer matters: dark states escape at constant rate","Cavity transport: dark states drain at fixed rate"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.00106,"raw_usage":{"total_tokens":4462,"prompt_tokens":975,"completion_tokens":3487,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":591,"completion_tokens_details":{"reasoning_tokens":3406}},"tokens_in":591,"tokens_out":3487,"duration_ms":22086,"temperature":1.0,"reasoning_tokens":3406,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-08T15:32:38.819010+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"A time-resolved photoluminescence or transient-absorption experiment on arrays with $N$ from about 4 to 64 could test the claim: if the dark-state decay time grows with $N$, or the decay is not single-exponential, or the extracted escape rate $\\gamma_{\\mathrm{rec}}$ changes when the detuning is moved across a phonon mode in a way not captured by $2J(\\Delta)$, the unidirectional-valve picture would be contradicted. In particular, measuring the escape rate for two different dark states with the same detuning but different energies within a structured spectral density would reveal whether the sum-rule collapse holds.","supporting_citations":[{"cited_title":"Seibt and T","cited_arxiv_id":null,"evidence_quote":"Supplies the hierarchical-equations treatment of Herzberg–Teller/non-Condon effects used to motivate the non-Condon channel."},{"cited_title":"Takahashi, K","cited_arxiv_id":null,"evidence_quote":"Reports light emission from vibronic polaritons in coupled metalloporphyrin–multimode cavity systems, the experimental platform where non-Condon coupling is active."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Provides evidence distinguishing phonon-modulated emission and absorption in porphyrin cavity polariton samples, supporting the resonance signature."},{"cited_title":"Xiang, R","cited_arxiv_id":null,"evidence_quote":"Baseline Condon-coupling theory showing dark–bright relaxation in polariton systems, which the paper contrasts with its unidirectional rescue."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Bosonic-mapping analysis of radiative pumping versus vibrational relaxation that the paper compares against its $1/N$-independent rescue rate."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Supplemental derivation of the Born–Markov reduction, the Redfield-tensor unidirectionality, and the validity condition Eq. (S8)."}],"review_version":1}