{"id":"a795a54a-c2dd-4ad5-b84f-58be651b1854","arxiv_id":"2509.03067","paper_version":2,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":6.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":1,"one_line_summary":"Vibrational coupling in organic emitters can preserve, and for red-detuned cavities even enhance, dynamical superradiance, with a characteristic asymmetry in the photon-rise time vs detuning.","lead":"This paper develops theory for how molecular vibrations affect collective light emission, superradiance, in organic materials inside a cavity. It predicts that for certain cavity frequencies, vibrations can actually speed up the collective emission, and it offers a clear experimental signature.","discovery_kind":"new_application","skeptic_critique":{"model":"deepseek-v4-flash","headline":"Mean-field closure for local vibrational modes is uncontrolled in HTC; predicted risetime asymmetry may be an artifact.","rationale":"The reader's verdict correctly flags the extrapolation from validated TC benchmarks to HTC results. I agree that the θ=π/4 validation does not directly cover θ=10^-3π, but the more load-bearing issue is the mean-field factorization of the local vibrational mode. In the TC model, mean-field becomes exact in the N→∞ limit because the cavity field is a collective bosonic mode whose fluctuations vanish; benchmarking at θ=π/4 and seeing cumulants approach mean-field with N supports this. In the HTC model, each molecule carries its own vibrational mode. The fluctuations of that mode around its mean are independent of N: they are O(1) for every molecule. The Holstein term √S(b+b†)σ_z is linear in the bosonic displacement and commutes with the electronic population; the resulting ground state is a polaron with electronic-vibrational entanglement. Mean-field discards these correlations by factorizing ⟨σ_z b⟩. For S=0.2 and ων=0.15 eV, the polaron shift (Sων ≈ 0.03 eV) is comparable to the detunings where the enhancement is claimed (|Δ| = 0.15–0.3 eV), and the Franck-Condon sidebands invoked in the physical explanation are entirely absent from the mean-field equations. N→∞ cannot help because the vibrational modes are local and not collective. Therefore the central asymmetry prediction may be an artifact of the closure rather than a real vibronic effect. The concrete second-order cumulant test directly addresses this: if including ⟨σ_z b⟩ correlations changes the τ(Δ) shift, the mean-field result is unreliable. I do not recommend rejecting the paper; the exact TC method is a solid contribution, and the HTC prediction is falsifiable. The verdict should remain conditional on this additional validation.","tokens_in":17900,"tokens_out":11031,"duration_ms":129573,"concrete_test":"Run the second-order cumulant equations for the HTC model (same QuantumCumulants.jl framework as Sec. II) at N=100 and N=1000, with θ=10^-3π and parameters of Fig. 6, and extract the risetime τ(Δ) for S=0.2. If the cumulant result does not converge to the mean-field curve—in particular, if the minimum of τ does not shift to negative Δ or the enhancement disappears—the mean-field prediction underlying the central claim is not reliable.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The central claim that moderate vibrational coupling enhances dynamical superradiance for negative detunings is obtained solely from the mean-field solution of the HTC model in Sec. III. The approximation factorizes each molecule's electronic and vibrational degrees of freedom, i.e. ⟨σ_z(b+b†)⟩ ≈ ⟨σ_z⟩⟨b+b†⟩. This is not justified by N→∞ because the vibrational mode is local to each molecule; its quantum fluctuations are O(1) and do not vanish in the thermodynamic limit. For S=0.1–0.4 and ων=0.15 eV, such fluctuations produce polaron dressing and Franck-Condon sidebands, and the paper's own explanation of the enhancement (Sec. III.B: 'cavity frequency matches transitions from the excited state manifold with zero vibrational excitations to the ground-state manifold with vibrational excitations') appeals to these discrete vibronic states. A mean-field calculation, which treats vibrations only through their mean displacement, cannot represent this multi-level vibronic structure; it yields at most a classical Stark shift of the electronic transition. The benchmark in Sec. II validates mean-field only for the TC model (no vibrations) and only for θ=π/4, so it does not establish the validity of the factorization for HTC at θ=10^-3π. If the risetime asymmetry disappears or changes sign once second-order vibrational correlations are included, the headline prediction is an artifact of the mean-field closure.","agreement_with_reader":"partial"},"referee_report":{"model":"deepseek-v4-flash","summary":"The paper studies dynamical superradiance in a single-mode cavity coupled to N emitters, comparing two models for vibrational effects: Markovian pure dephasing and explicit Holstein–Tavis–Cummings (HTC) vibrational dressing. It introduces an exact numerical method for the dissipative Tavis–Cummings (TC) master equation up to N ≈ 140, exploiting weak permutation and weak U(1) symmetries, and benchmarks mean-field and second-order cumulant approximations against it. The manuscript then uses mean-field theory to simulate the HTC model at N = 10^8, predicting that moderate vibrational coupling can enhance superradiance for negative cavity detunings and produce an asymmetric photon-number risetime as a function of detuning, in contrast to pure dephasing. The proposed asymmetry is presented as an experimentally accessible signature of vibrationally assisted superradiance.","tokens_in":18163,"tokens_out":7765,"duration_ms":94200,"significance":"If the HTC mean-field prediction is robust, the predicted risetime asymmetry is a concrete, experimentally testable signature for organic microcavity systems, and the paper's exact numerical solver for dissipative TC models is a valuable technical contribution with code availability. However, the headline HTC result rests on a mean-field closure that is not benchmarked for the HTC model, and the paper's own mechanism explanation invokes vibronic structure that the closure cannot represent. The significance of the paper therefore depends on additional validation of the HTC mean-field treatment.","major_comments":[{"comment":"The mean-field approximation factorizes each molecule's electronic and vibrational operators, e.g., ⟨σ_z(b+b†)⟩ ≈ ⟨σ_z⟩⟨b+b†⟩. For local vibrational modes this closure is uncontrolled: N→∞ does not suppress O(1) local vibrational quantum fluctuations, and the benchmark in Sec. II validates mean-field only for the TC model (S=0), not for HTC. Moreover, the proposed enhancement mechanism in Sec. III.B ('cavity frequency matches transitions from the excited state manifold with zero vibrational excitations to the ground-state manifold with vibrational excitations') is a vibronic Franck–Condon mechanism that a factorized mean-field treatment cannot represent. To support the central claim, the authors should benchmark HTC mean-field against a second-order cumulant or exact small-N calculation, or derive the result in a controlled limit (e.g., polaron frame, or classical cavity field with exact","section":"Sec. III.A–III.B, Eq. (9)–(11)"},{"comment":"Mean-field is validated against exact results only for SR initial angle θ=π/4 (Fig. 2d), while the HTC simulations use θ=10^{-3}π. This is a significant extrapolation: for θ=0 (SF), mean-field gives identically zero photon emission, so the small-θ regime is singular. The benchmark at θ=π/4 does not establish that mean-field remains reliable at θ=10^{-3}π. The authors should either provide a TC benchmark at small θ, or compare HTC mean-field to cumulant/exact results at small θ, before using this parameter regime for the central prediction.","section":"Sec. II.C vs. Sec. III.B"},{"comment":"The definition of the risetime τ is not sufficiently precise. The text says τ is extracted from 'an exponential fit to the linear regions' of ⟨n⟩/N on a log scale, but the fitting interval, the criterion for what counts as a linear region, and any statistical uncertainty are not given. Since the asymmetry in Fig. 6(b) is the paper's main experimental signature, the extraction procedure should be specified quantitatively and error bars or confidence intervals should be provided, at least for representative parameter sets.","section":"Sec. III.B and Fig. 6"}],"minor_comments":[{"comment":"The caption says '(c), (d), (f) Same as (a), (b), (c)' but should presumably read '(d), (e), (f) Same as (a), (b), (c)'.","section":"Fig. 4 caption"},{"comment":"The y-axis label in the inset appears garbled ('10□1'); please correct the typographical error.","section":"Fig. 5 inset"},{"comment":"The Lamb-shift term H_LS = i(γν/4)(b_i†² − b_i²) is unusual; a brief physical justification beyond the citation to Ref. [82] would improve readability.","section":"Eq. (10)"},{"comment":"It would be helpful to state explicitly that only the ν=ν′ diagonal block is implemented in the released PIBS code, and whether this limits which observables can be computed with the public version.","section":"Sec. II.B"}],"recommendation":"major_revision","confidential_remarks":"The exact TC solver and the cumulant benchmarks are solid contributions. The main risk is the uncontrolled mean-field closure for the HTC model: the central asymmetry prediction is not supported by the current benchmarks, and the reported mechanism is inconsistent with a factorized treatment. I would be willing to accept the paper after the authors supply either a controlled HTC benchmark (e.g., second-order cumulant comparison or exact small-N results) or a derivation that justifies the closure in the relevant parameter regime, and after clarifying the risetime extraction."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"The reader's conditional verdict is sensible, and the stress-test concern lands. The exact numerical method (PIBS), combining permutation and weak U(1) symmetry to push dissipative TC dynamics up to N=140, is a real contribution. Code and data are shared, and the benchmarks against mean-field and second-order cumulants for the TC model are clean and useful. The qualitative contrast between pure dephasing and explicit vibrations, such as the Bloch vector length dynamics, would be interesting even if the enhancement claim were wrong.\n\nThe soft spot is exactly where the reader's weakest assumption sits. Section III uses mean-field with N=10^8 and theta=10^-3 pi. The validation in Section II covers only the TC model (no vibrations) and only theta=pi/4. The cumulant-to-mean-field convergence in Fig. 3 is also TC-only. More importantly, each local vibrational mode has O(1) quantum fluctuations regardless of N; there is no thermodynamic-limit argument that justifies factorizing <sigma_z(b+b^dagger)>. The paper's own explanation in Section III.B invokes transitions from the excited state with zero vibrational quanta to the ground state with vibrational quanta. That is a discrete vibronic mechanism, not something a mean-field that replaces the vibration by its average displacement can produce. So the risetime asymmetry in Fig. 6 is at risk of being an artifact of the closure.\n\nThis is not a case where the paper is incoherent or dishonest; the authors explicitly state they use mean-field for the HTC model. But the conclusion overclaims when it says the TC benchmarks 'justified the application of mean-field theory' to the vibrational case. The leap is unstated there.\n\nWho should read this: people working on superradiance in organic cavities will want the exact method and the dephasing-versus-vibration comparison. People planning to use the enhancement prediction as a guide to experiments should wait until the HTC mean-field is checked, for example by a second-order cumulant or a small-N exact solution of the vibrational model.\n\nRecommendation: send to peer review. A good referee can ask for that benchmark, and the exact method deserves to be published regardless.","headline":"Useful exact method for the dissipative Tavis-Cummings model, but the headline claim of vibrationally assisted superradiance rests on an unvalidated mean-field closure.","tokens_in":18691,"tokens_out":4433,"would_cite":false,"duration_ms":50484,"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":"Vibrational coupling in organic emitters can enhance, not just suppress, dynamical superradiance when the cavity is tuned below the molecular transition.","keywords":["dynamical superradiance","superfluorescence","Holstein-Tavis-Cummings model","vibrational coupling","cavity quantum electrodynamics","permutational invariance","mean-field approximation","photon-number risetime"],"falsifier":"Measure the early exponential photon-number risetime in an organic microcavity as the cavity frequency is swept across the molecular transition. If the risetime remains symmetric about zero detuning at all Huang–Rhys couplings, or if no parameter regime shows a negative-detuning risetime shorter than the S=0 curve, the vibrationally assisted enhancement is falsified. A complementary calculation: an exact small-N solution of the Holstein–Tavis–Cummings model at θ=10^-3π that disagrees with mean-field would undercut the extrapolation.","tokens_in":17786,"feed_emoji":"⚛️","tokens_out":6542,"duration_ms":68495,"temperature":0.7,"pith_summary":"This paper develops the theory of dynamical superradiance—collective emission from many excited molecules into a single cavity mode—for organic emitters whose electronic states are dressed by molecular vibrations. The authors show that treating vibrations as mere dephasing is insufficient: in the Holstein–Tavis–Cummings model, moderate vibrational coupling can actually speed up the photon-number rise when the cavity frequency sits below the molecular transition. They support this with an exact numerical method for the dissipative Tavis–Cummings model up to 140 emitters, use it to validate mean-field and second-order cumulant approximations, then apply mean-field theory to macroscopic emitter numbers. If correct, the work provides an experimentally accessible signature—an asymmetry in the photon-number risetime versus cavity detuning—that distinguishes coherent vibronic coupling from Markovian dephasing.","feed_headline":"Vibrational coupling can speed up superradiance in organic cavities","feed_subtitle":"When the cavity sits below the molecular resonance, moderate vibronic coupling shortens the photon-number risetime—a signature experiments c","key_machinery":"The load-bearing objects are two complementary constructions. First, the Holstein–Tavis–Cummings Hamiltonian (two-level electronic states coupled to a common cavity mode and, per molecule, to a local harmonic vibrational mode with Huang–Rhys coupling S) supplies the physical mechanism: negative detuning lets the cavity pick out vibrationally dressed transitions. Second, an exact solution method for the dissipative Tavis–Cummings model combines weak permutation symmetry (emitters are interchangeable) with weak U(1) symmetry (the master equation is invariant under rotating photon and dipole phases together); this block structure lets the authors integrate the excitation-conserving diagonal sec","core_discovery":"The central claim is that vibrational coupling in organic emitters changes dynamical superradiance in a way that goes beyond simple dephasing. In the Holstein–Tavis–Cummings model, the photon-number rise time grows with coupling for small negative detunings, but for large negative detunings moderate coupling (Huang–Rhys parameter S between about 0.1 and 0.4) shortens the risetime below its vibration-free value. The explanation is that a negatively detuned cavity can resonantly connect the excited-state manifold with zero vibrational quanta to the ground-state manifold with vibrational quanta; increasing S makes one-vibration transitions more probable, so the cavity emission is assisted rathe","pith_inferences":["A natural experimental next step is to scan cavity detuning in an organic microcavity while fitting the early exponential rise; the predicted asymmetric risetime would be evidence for coherent vibronic dressing, while a symmetric curve would point to pure dephasing.","Because the exact solver is formulated for two-level emitters and a single cavity mode, extending the same block structure to multi-level emitters or multi-mode cavities could test whether vibrational assistance survives those complications.","The paper validates mean-field on the vibration-free model at θ=π/4 but applies it to the vibration-coupled model at θ=10^-3π; a finite-N exact or cumulant check at the smaller angle would tighten the quantitative risetime predictions.","If the enhancement mechanism is correct, varying the vibrational frequency relative to the detuning should produce a resonant peak in the risetime speed-up, which could be used to calibrate the Huang–Rhys parameter in situ."],"forward_implications":["Dynamical superradiance survives in organic molecules at realistic Huang–Rhys parameters (up to roughly S=0.1), so organic microcavities are viable platforms for observing it.","For negative cavity detunings, moderate vibrational coupling can shorten the exponential photon-number risetime compared with the vibration-free case, an enhancement that pure-dephasing models cannot produce.","The risetime-versus-detuning curve is symmetric for Markovian dephasing and asymmetric for coherent vibronic coupling; measuring that asymmetry distinguishes the two microscopic pictures.","S beyond about 0.3 suppresses the exponential rise and inhibits electronic-coherence buildup, so strong vibrational dressing destroys dynamical superradiance.","The exact block solver extends benchmark-quality solutions to 140 emitters, roughly five times larger than the standard permutation-symmetry limit, providing a tool to test other approximations."],"supporting_citations":[{"why":"Supplies the mean-field solution of the Tavis–Cummings model, the baseline dynamical-superradiance dynamics this work extends.","marker":"[20]"},{"why":"Supplies the weak-permutation-symmetry compression that underpins the exact solver's emitter-state representation.","marker":"[57–63]"},{"why":"Supplies the weak-U(1) block structure that lets the solver integrate sequentially and reach N≈140 emitters.","marker":"[64]"},{"why":"Supplies the Holstein–Tavis–Cummings Hamiltonian used to model vibrational dressing of the emitters.","marker":"[47, 56]"},{"why":"Justifies mean-field accuracy in the thermodynamic limit, motivating the large-N vibration-coupled calculations.","marker":"[65–67]"},{"why":"Supplies the software used to generate and integrate the mean-field and second-order cumulant equations.","marker":"[80]"},{"why":"Provides the finite-N quantum-correction estimate used to argue that early-time mean-field agreement extends to macroscopic emitter numbers.","marker":"[34]"}],"fun_headline_variants":["Vibrational coupling boosts superradiance in organic cavities","Organic superradiance aided by vibrations at negative detuning","Vibronic coupling shortens superradiance risetime in organic systems","Superradiance enhanced by vibrations in organic materials","Vibrational assistance speeds superradiance in organic cavities"],"cache_read_input_tokens":2688,"weakest_assumption_plain":"The argument depends on mean-field theory, validated only for the vibration-free Tavis–Cummings model at a large initial coherence angle, also being accurate for the vibration-coupled model at a much smaller initial coherence angle and macroscopic emitter number.","fun_headline_variants_meta":{"raw":{"variants":["Vibrational coupling boosts superradiance in organic cavities","Organic superradiance aided by vibrations at negative detuning","Vibronic coupling shortens superradiance risetime in organic systems","Superradiance enhanced by vibrations in organic materials","Vibrational assistance speeds superradiance in organic cavities"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000425,"raw_usage":{"total_tokens":2036,"prompt_tokens":782,"completion_tokens":1254,"prompt_tokens_details":{"cached_tokens":256},"prompt_cache_hit_tokens":256,"prompt_cache_miss_tokens":526,"completion_tokens_details":{"reasoning_tokens":1172}},"tokens_in":526,"tokens_out":1254,"duration_ms":9356,"temperature":1.0,"reasoning_tokens":1172,"cache_read_input_tokens":256,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-05T11:05:42.729595+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Measure the early exponential photon-number risetime in an organic microcavity as the cavity frequency is swept across the molecular transition. If the risetime remains symmetric about zero detuning at all Huang–Rhys couplings, or if no parameter regime shows a negative-detuning risetime shorter than the S=0 curve, the vibrationally assisted enhancement is falsified. A complementary calculation: an exact small-N solution of the Holstein–Tavis–Cummings model at θ=10^-3π that disagrees with mean-field would undercut the extrapolation.","supporting_citations":[{"cited_title":"Bonifacio and G","cited_arxiv_id":null,"evidence_quote":"Supplies the mean-field solution of the Tavis–Cummings model, the baseline dynamical-superradiance dynamics this work extends."},{"cited_title":"Plankensteiner, C","cited_arxiv_id":null,"evidence_quote":"Supplies the software used to generate and integrate the mean-field and second-order cumulant equations."},{"cited_title":"Keeling, Quantum corrections to the semiclassical col- lective dynamics in the Tavis-Cummings model, Physical Review A 79, 053825 (2009)","cited_arxiv_id":null,"evidence_quote":"Provides the finite-N quantum-correction estimate used to argue that early-time mean-field agreement extends to macroscopic emitter numbers."}],"review_version":1}