{"id":"714f1064-d152-4272-88e1-8d33015f3436","arxiv_id":"2605.24227","paper_version":1,"verdict":"UNVERDICTED","confidence":"LOW","novelty_score":6.0,"correctness_risk":"unknown","formal_verification":"none","parameter_count":0,"one_line_summary":"In the Tavis-Cummings model, collective dynamics are linear at vanishing excitation density and nonlinear at finite density, with mean-field accurate for large N while cluster expansion captures finite-N effects.","lead":"This paper maps the regimes where mean-field and single-excitation approximations agree or diverge in collective light-matter systems by using two parameters: molecule number N and excitation density N_exc/N in the Tavis-Cummings model. A smart generalist might read it to understand when simplified models for molecular ensembles in optical cavities remain reliable versus when nonlinearity appears.","discovery_kind":"extension","skeptic_critique":{"model":"grok-4.3","headline":"Representativeness of Tavis-Cummings (even with vibronics) for real collective dynamics is untested","rationale":"The reader's weakest assumption is exactly the load-bearing step for extending the model-derived map to the broader claim about collective light-matter dynamics. Internal consistency of the MF/SE analysis within Tavis-Cummings is not challenged here; the concern is only the missing bridge to more complete descriptions or experiment.","tokens_in":1848,"tokens_out":323,"duration_ms":37711,"concrete_test":"Take the MF Duffing equation derived for finite N_exc/N, insert parameters from a Tavis-Cummings instance with added counter-rotating terms at N=200 and N_exc/N=0.3, and compare the resulting anharmonic frequency shift against exact diagonalization in the conserved-excitation subspace; if the shift differs by >15% the regime map does not transfer.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The central claim is a two-parameter regime map for when MF and SE descriptions are controlled. This map is derived inside the Tavis-Cummings Hamiltonian (and a variant with local vibronic coupling). The abstract states that the large-N limit at finite excitation density remains accurately described by MF and yields a Duffing equation. For the claim to control real light-matter systems, the model must capture the essential physics; no comparison to Hamiltonians with counter-rotating terms, disorder, or collective decay is supplied to establish this. The cluster-expansion and polaron-decoupling statements are internal to the chosen model.","agreement_with_reader":"agree"},"referee_report":{"model":"grok-4.3","summary":"The paper claims that in the Tavis-Cummings model, collective light-matter dynamics are governed by two parameters, N and the excitation density N_exc/N. When N is large and the excitation density approaches zero, mean-field (MF) and single-excitation (SE) approximations agree, yielding linear collective dynamics with harmonic Rabi oscillations. At finite excitation density, the large-N limit is described by MF dynamics that become nonlinear, leading to a Duffing equation for the cavity amplitude with anharmonic Rabi frequency. Cluster expansion is shown to restore finite-N correlations beyond MF. The same linear limit is reached when local vibronic interactions are included, with SE via polaron decoupling and MF via linearization. This provides a regime map for the validity of different theoretical descriptions.","tokens_in":1974,"tokens_out":385,"duration_ms":51853,"significance":"If the results hold, this work is significant for providing a systematic two-parameter characterization of approximation regimes in collective light-matter systems. It highlights the conditions under which linear vs nonlinear dynamics emerge and demonstrates how standard techniques like cluster expansion can be used to go beyond mean-field. The extension to vibronic interactions shows robustness of the linear limit. This could aid in selecting appropriate models for polariton chemistry and related fields. The analysis is internal to the Tavis-Cummings model and its variant, which is appropriate for the stated claims; the stress-test concern regarding representativeness for real systems does not land because the derivations and claims are explicitly scoped to the chosen Hamiltonian.","major_comments":[],"minor_comments":[{"comment":"Abstract: the central results about regime agreement and the Duffing equation are stated without referencing the specific sections or equations in the main text where the derivations, error analysis, or explicit calculations appear; adding such pointers would improve navigability.","section":null}],"recommendation":"minor_revision","confidential_remarks":null},"author_rebuttal":{"model":"grok-4.3","summary":"We thank the referee for their careful reading and positive assessment of our manuscript, including the accurate summary of our two-parameter regime characterization in the Tavis-Cummings model and its extension to vibronic interactions. The recommendation for minor revision is noted, though no specific major comments were provided in the report.","responses":[],"tokens_in":1361,"tokens_out":79,"duration_ms":17155,"standing_objections":[]},"desk_editor":{"model":"grok-4.3","letter":"The main thing to know is that this work draws a clean two-parameter map for when mean-field and single-excitation descriptions line up in the Tavis-Cummings model. At large N and vanishing excitation density both approximations recover linear Rabi oscillations. At finite density the large-N limit stays mean-field accurate but turns nonlinear, showing up as a Duffing equation for the cavity field. With added local vibronics the linear limit is reached by different routes in each approximation.\n\nThe paper does a straightforward job laying out these limits and showing how cluster expansion restores finite-N correlations beyond mean-field. The identification of the Duffing form at order-one density and the contrasting mechanisms for reaching linearity with vibronics are useful organizational points that were not spelled out this way before.\n\nThe soft spot is exactly the one flagged in the stress-test note. Everything is derived inside the Tavis-Cummings Hamiltonian. There is no check against versions that include counter-rotating terms, disorder, or collective decay, so it is not shown that the same regime boundaries control real light-matter systems. The claims about controlled descriptions therefore remain internal to the chosen model.\n\nThis is for theorists who routinely pick between mean-field and single-excitation treatments in cavity QED or polariton chemistry and want a practical guide for parameter choices. The reasoning is coherent and the scope is stated clearly, so the paper deserves a serious referee even though its reach is limited to one standard Hamiltonian.","headline":"The paper maps when mean-field and single-excitation approximations agree or diverge in the Tavis-Cummings model using N and excitation density as parameters.","tokens_in":2492,"tokens_out":368,"would_cite":false,"duration_ms":27672,"reading_group":"maybe","serious_thinker":"yes","would_accept_peer_review":true},"rs_alignment":null,"lean_confirmation":null,"pith_extraction":{"msc":[],"pacs":[],"model":"grok-4.3","headline":"In the Tavis-Cummings model, mean-field and single-excitation descriptions agree on linear collective dynamics only at near-zero excitation density; at finite density the large-N limit stays mean-field but turns nonlinear.","keywords":["Tavis-Cummings model","collective light-matter dynamics","mean-field approximation","single-excitation approximation","Rabi oscillations","Duffing equation","excitation density","polariton chemistry"],"falsifier":"A numerical simulation or cavity experiment at large but finite N and N_exc / N of order one that shows the cavity amplitude obeying the predicted Duffing equation rather than linear Rabi motion.","tokens_in":2721,"feed_emoji":"","tokens_out":749,"duration_ms":37444,"temperature":0.7,"pith_summary":"The paper maps the validity of common approximations for collective light-matter coupling by treating molecule number N and excitation number N_exc as independent controls. When N is large but excitations remain dilute, both mean-field and single-excitation treatments converge on the same linear response with ordinary harmonic Rabi oscillations. Once the excitation density reaches order one, mean-field continues to describe the system accurately yet the motion becomes nonlinear, appearing as a Duffing equation for the cavity amplitude. The same linear limit is recovered when local vibrations are added, though each approximation reaches it by a different route. This supplies a practical two-parameter chart for choosing controlled theoretical descriptions.","feed_headline":"Excitation density sets linear or nonlinear collective light-matter regimes","feed_subtitle":"Large molecule number plus dilute excitations yields harmonic Rabi motion; finite density produces Duffing nonlinearity while mean-field sta","key_machinery":"Two-parameter regime map in N and N_exc that separates the linear collective limit from the nonlinear finite-density regime inside the Tavis-Cummings model.","core_discovery":"In the Tavis-Cummings model, when N ≫ 1 and the excitation density N_exc / N → 0, MF and SE descriptions agree and yield linear collective dynamics, showing harmonic Rabi oscillations. At finite excitation density (N_exc / N ∼ O(1)), the large-N limit remains accurately described by MF dynamics but becomes nonlinear in N_exc / N, manifested by a Duffing equation for the cavity amplitude with anharmonic Rabi frequency. Cluster expansion systematically restores finite-N correlations beyond MF. When local vibronic interactions are included, the same linear collective limit is reached by both approximations, with SE reaching it through polaron decoupling and MF through linearization.","pith_inferences":["The regime map suggests that excitation-density tuning could be used experimentally to cross from linear to nonlinear collective response while keeping N large.","Similar density-dependent crossovers may appear in other collective models once both N and N_exc are treated as independent variables.","The Duffing description offers a concrete starting point for analytic or numeric studies of anharmonic effects in polariton systems at moderate excitation levels."],"forward_implications":["At low excitation density the linear collective response is recovered by both MF and SE.","Finite excitation density produces an anharmonic Rabi frequency through the Duffing dynamics of the cavity field.","Cluster expansion supplies controlled corrections that restore correlations missing from pure mean-field.","Local vibronic coupling does not destroy the linear collective limit but routes the two approximations to it differently."],"fun_headline_variants":["Excitation density splits collective dynamics into linear and nonlinear","Tavis-Cummings regimes depend on N and excitation density","Finite density leads to anharmonic Rabi in mean-field limit","Both MF and SE reach linear limit at vanishing excitation density"],"cache_read_input_tokens":2112,"weakest_assumption_plain":"The Tavis-Cummings model, with or without added local vibronic interactions, captures the essential physics of collective light-matter systems.","fun_headline_variants_meta":{"raw":{"variants":["Excitation density splits collective dynamics into linear and nonlinear","Tavis-Cummings regimes depend on N and excitation density","Finite density leads to anharmonic Rabi in mean-field limit","Both MF and SE reach linear limit at vanishing excitation density"]},"model":"grok-4.3","cost_usd":0.005892,"raw_usage":{"total_tokens":2835,"prompt_tokens":741,"num_sources_used":0,"completion_tokens":66,"cost_in_usd_ticks":58924500,"prompt_tokens_details":{"text_tokens":741,"audio_tokens":0,"image_tokens":0,"cached_tokens":256},"completion_tokens_details":{"audio_tokens":0,"reasoning_tokens":2028,"accepted_prediction_tokens":0,"rejected_prediction_tokens":0}},"tokens_in":741,"tokens_out":66,"duration_ms":23050,"temperature":1.0,"reasoning_tokens":2028,"cache_read_input_tokens":256,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-06-30T14:11:07.111087+00:00","model_set":{"reader":"grok-4.3"},"falsifier":"A numerical simulation or cavity experiment at large but finite N and N_exc / N of order one that shows the cavity amplitude obeying the predicted Duffing equation rather than linear Rabi motion.","supporting_citations":[],"review_version":1}