{"id":"112cb0a8-e3e3-488b-a327-8227ac250abc","arxiv_id":"2507.00918","paper_version":1,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":6.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":1,"one_line_summary":"Static energetic disorder slows lower and upper polaritons, but for realistic organic microcavities the effect is too small to explain observed slowdown, pointing to phonon scattering as the dominant mechanism.","lead":"This paper calculates how static energetic disorder in molecular microcavities changes polariton group velocities. It finds the effect is usually modest, supporting the view that vibrational, phonon-driven scattering, not static disorder, explains recent slow-propagation experiments.","discovery_kind":"new_application","skeptic_critique":{"model":"deepseek-v4-flash","headline":"The large-slowdown predictions occur in the δq/q ≥ 1 regime where the paper itself says Eq. (4) is inapplicable, so the quantitative strength of the central claim is not established; the reader's conditional verdict is appropriate.","rationale":"I agree with the reader's weakest assumption: the load-bearing step is the identification of the physical group velocity with hbar^-1 d Re[EP(q)]/dq in the presence of static disorder. What must be true for the central claim is that Eq. (4) remains a valid measure of wavepacket velocity in the high-disorder, high-exciton regime where the largest slowdowns are reported. The paper itself states the opposite: modes with δq/q ≥ 1 are outside the applicability of Eq. (3), and the large-slowdown contours in Fig. 3 appear precisely in that excluded region. This makes the concern concrete and internal to the manuscript, not merely a disagreement with consensus. The weak-disorder perturbative analysis (Eqs. 7–9) supports a monotonic reduction of group velocity in the valid regime, and the perovskite-specific conclusion uses σ/ΩR < 0.235, which is more secure. However, the abstract's broad statement that static disorder appreciably slows polaritons when fluctuations approach the coupling strength is quantitatively unsupported where it matters most. The paper is honest about the boundary and hedges its language, so conditional acceptance remains appropriate. A direct spectral-function or exact-diagonalization check would settle whether the complex-branch derivative underestimates or overestimates the true wavepacket velocity, and would either rescue or remove the large-slowdown portion of the claim. The reader's conditional verdict should therefore stand unchanged.","tokens_in":7906,"tokens_out":6344,"duration_ms":85344,"concrete_test":"Build a finite 2D square lattice of N×N two-level molecules with Gaussian site energies (σ/ΩR = 0.5), periodic boundary conditions, and a cavity photon branch EC(q); diagonalize H for ~10^3 disorder realizations. From the disorder-averaged spectral function A(q,ω), compute (i) hbar^-1 d Re E_LP/dq from Eq. (3), and (ii) the spectrally weighted velocity v_sw(q) = ∫dω A(q,ω)(1/hbar)dEC(q)/dq / ∫dω A(q,ω), plus (iii) the wavepacket centroid velocity in real-time propagation initialized at q0. If (ii) or (iii) agrees with (i) within 20% for modes with δq/q ≥ 1, the concern is resolved; if they deviate, the large-slowdown claim is an artifact of using the complex-branch derivative outside its validity domain.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The paper's quantitative message is that static disorder can appreciably reduce polariton group velocities, especially for high-exciton-content modes when σ/ΩR is large (Abstract; Figs. 2–3). The operative definition, Eq. (4), identifies vg with hbar^-1 d Re[EP(q)]/dq, and the paper states (following Eqs. 4–5) that this is valid only if the polariton spectral function is narrow and δq/q < 1, with δq = Im[EP]/(hbar vg q). The largest predicted slowdowns are precisely in LP modes at large q, where the LP is predominantly excitonic. The text explicitly concedes that these modes have δq/q ≥ 1 and that Eq. (3) is no longer applicable there. Therefore the 'appreciable slowing' regime is outside the domain in which the group velocity is defined; the derivative of Re EP in that regime need not equal the velocity of a propagating wavepacket, which may instead be governed by spectral-weight transport, diffusive motion, or localization. This does not undermine the weak-disorder perturbative argument (Eqs. 7–9) or the conclusion that, for perovskite parameters (σ/ΩR < 0.235), static disorder is a minor effect; both live in the valid regime. But the abstract's general claim, and any quantitative use of the Fig. 3 red-contour region, require an independent check that the quasiparticle velocity is meaningful there.","agreement_with_reader":"agree"},"referee_report":{"model":"deepseek-v4-flash","summary":"The manuscript develops a nonperturbative scheme, based on the effective-medium dispersion relation of Litinskaya and Reineker, to compute how static Gaussian energetic disorder renormalizes the complex energies and group velocities of microcavity exciton-polaritons. The authors solve the resulting integral equation numerically and also derive a weak-disorder perturbative expansion. Using parameters from perovskite microcavity (Xu et al.) and BODIPY-BSW (Balasubrahmanyam et al.) experiments, they find that static disorder lowers the LP branch and raises the UP branch, reducing group velocities in both branches, with appreciable reductions only for strongly excitonic modes when the disorder width approaches half the Rabi splitting. The paper concludes that static disorder accounts for at most a minor fraction of the experimentally observed slowdown and that dynamical (phonon-assisted) disorder is the dominant mechanism.","tokens_in":8188,"tokens_out":6361,"duration_ms":69936,"significance":"The paper addresses a timely question in polariton transport: whether static energetic disorder can explain the slow propagation observed in molecular microcavities. Its strengths are the minimal parameter set, the absence of any fitted target quantities, and the internal consistency between the numerical solution and the perturbative expansion. The conclusion that static disorder plays a minor role for the specific perovskite parameters (σ/ΩR < 0.235) is well supported and lives in the regime where the group velocity definition is valid. However, the more dramatic quantitative claims of appreciable slowdown, emphasized in the abstract and Fig. 3, are made in a regime where the paper's own validity criterion δq/q < 1 is violated, so the strength of the central claim is not fully established by the present evidence.","major_comments":[{"comment":"The paper defines the renormalized group velocity as ħ⁻¹∂Re[EP(q)]/∂q and states that this is valid only when the disorder-induced broadening satisfies δq/q < 1. It then acknowledges that in LP modes at large q and UP modes at small q, δq/q ≥ 1 and the estimate is no longer applicable. Yet Fig. 3 presents quantitative velocity renormalizations exceeding 10% in exactly that regime, and the abstract's emphasis on σ/ΩR approaching 0.5 refers to conditions where many high-exciton modes fall outside the validity domain. The claim that static disorder can appreciably slow polaritons therefore rests on extrapolating Eq. (4) beyond its stated range. To support the quantitative message, the authors should either restrict the strong-renormalization claims to the δq/q < 1 region (where the effect is modest) or provide an independent check, such as wavepacket-propagation simulations or a spectral-function-based transport velocity, demonstrating that ∂Re[EP]/∂q remains the physical group velocity in the strongly broadened regime.","section":"Results and Discussion (Eqs. 4–5, Fig. 3)"},{"comment":"The perturbative argument leading to the conclusion that both vLP(q) < vLP⁽⁰⁾(q) and vUP(q) < vUP⁽⁰⁾(q) is presented as a derivation, but the reasoning for the UP branch is not pointwise: from δUP(q) ≥ 0 and δUP(q) → 0 as q → ∞, one can only infer that the derivative δ'UP(q) is negative on average, not at every q. Since the universal slowdown is a central claim, the authors should either quantify the derivative from Eq. (9) explicitly or clearly label this as a qualitative trend that is verified numerically in Fig. 2 rather than a proof.","section":"Results and Discussion (Eqs. 8–9)"}],"minor_comments":[{"comment":"The sentence 'these modes have δq/q ≥ 1 and therefore, Eq. 3 is no longer applicable' appears to contain a typo: Eq. (3) is the exact equation used to compute the complex energies and remains applicable for any δq. The statement likely refers to Eq. (4), the group-velocity definition, and should be corrected for clarity.","section":"Results and Discussion (after Eq. 5)"},{"comment":"The quantity PM(q), labeled 'exciton content' or 'exciton fraction', is used to define the axes of Fig. 3 and the crossover relation PM ∝ (σ/ΩR)⁻¹, but it is never defined explicitly in the text. Please state whether it is computed from the bare (σ = 0) Hopfield coefficients or from the renormalized dispersion, and specify the formula used.","section":"Fig. 3"},{"comment":"The introduction cites Ref. 44 and 45 as previous works showing that static disorder can slow polariton transport in one-dimensional wires; briefly stating the difference between those results and the present higher-dimensional calculation would help the reader understand the novelty.","section":"Introduction and Methods"},{"comment":"The statement 'data available from the corresponding author upon reasonable request' is acceptable, but depositing the numerical data underlying Figs. 1–3 in a repository would improve reproducibility.","section":"Data Availability"}],"recommendation":"major_revision","confidential_remarks":"The core message—that static disorder is not the primary cause of the observed polariton slowdown in the perovskite system—is likely correct and is supported by the valid-regime analysis. The main weakness is the mismatch between the paper's strong quantitative claims (large group-velocity renormalization at high exciton fraction and σ/ΩR ≈ 0.5) and its own validity criterion δq/q < 1. I recommend major revision rather than rejection because the issue is fixable by either qualifying the claims or adding a transport calculation in the strong-broadening regime. The manuscript would also benefit from correcting the apparent typo in the sentence about Eq. (3)/Eq. (4) and from defining PM(q)."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"Useful paper, with a real caveat. The genuinely new pieces are the application of the Litinskaya–Reineker effective-medium equation to 2D microcavity group velocities, the closed-form leading-order splitting renormalization (Eqs. 7–9), and the clean map of when static disorder matters as a function of σ/ΩR and exciton fraction. The calculation has no fitted parameters; all inputs come from experiment, and nothing is tuned to match the observed slowdown. That is worth underlining: the central claim that static disorder accounts for at most a minor fraction of the observed slowdown is a genuine prediction, not a fit.\n\nThe paper is also honest about its limits. It explicitly states that the group-velocity identification via d Re E/dq is only valid when δq/q < 1, and it marks that boundary in Fig. 3. A reader cannot accuse the authors of hiding it.\n\nThe soft spot is that the 'appreciable slowdown' region they emphasize—high exciton content, σ approaching half the Rabi splitting—is precisely where δq/q ≥ 1. In that regime the derivative of the real part of the pole does not, by their own criterion, correspond to a propagating wavepacket's velocity; spectral-weight transport or localization could set the actual speed. The abstract's 'particularly when the fluctuation approaches half the splitting' leans on that outside-domain region. This does not damage the weak-disorder perturbative result or the minor-role conclusion for realistic perovskite parameters, both of which live inside the valid domain, but it does mean the quantitative strength of their 'appreciable slowdown' message is not yet established.\n\nMinor: no code or data deposit, and no propagated error bars on measured parameters. For a theory paper of this type I would not weigh that heavily. The citation pattern is fair; the prior 1D work is theirs and is cited as such, and the effective-medium equation is properly credited to Litinskaya–Reineker.\n\nWho is this for: people working on polariton transport, especially experimentalists trying to interpret slow propagation. It deserves a serious referee. A good referee should ask for one consistency check in the high-δq region—e.g., a wavepacket propagation or a spectral-function computation—before the stronger claims are accepted. Even if that check softens the numbers, the qualitative conclusion that phonons dominate will stand.","headline":"Static disorder, applied to 2D microcavity polaritons, is shown to be a minor cause of observed slowdown; the paper's quantitative 'appreciable slowing' claims sit outside its own stated validity regime, but the weak-disorder analytics and the central conclusion hold up.","tokens_in":8713,"tokens_out":2989,"would_cite":true,"duration_ms":32829,"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":"Static disorder alone can account for only a minor fraction of the polariton slowdown reported in recent experiments, with phonon-assisted scattering the likely dominant mechanism.","keywords":["exciton-polaritons","group velocity","static disorder","inhomogeneous broadening","microcavity polaritons","Rabi splitting","phonon-assisted scattering","polariton transport"],"falsifier":"Simulate the real-space motion of a polariton wavepacket in a disordered microcavity by solving the full disordered microscopic model for the same $\\sigma/\\Omega_R$ values and compare its centroid velocity with $\\hbar^{-1}\\partial_q \\mathrm{Re}[E_P(q)]$; a mismatch in the $\\delta_q/q<1$ regime would invalidate the central method. On the experimental side, measuring low-temperature group velocities in two samples with identical Rabi splitting but different inhomogeneous broadening would directly test whether static disorder controls the slowdown.","tokens_in":7707,"feed_emoji":"🐌","tokens_out":15750,"duration_ms":164031,"temperature":0.7,"pith_summary":"Molecular exciton-polaritons are expected to travel across microcavities at speeds fixed by their dispersion, but recent experiments observe far slower propagation. This paper asks whether static energetic disorder—frozen sample inhomogeneities—can produce that slowdown. Using a nonperturbative effective-medium equation that needs only the exciton energy distribution, cavity dispersion, and Rabi splitting, the paper finds that static disorder does lower both lower- and upper-polariton group velocities, but appreciably only for strongly excitonic modes when the disorder width approaches half the Rabi splitting. With parameters from recent perovskite and BODIPY surface-wave experiments, the calculation accounts for at most a minor fraction of the measured velocity drop. The paper concludes that phonon-assisted dynamical disorder, not static disorder, dominates the reported polariton slowdown.","feed_headline":"Static disorder explains only a minor share of polariton slowdown","feed_subtitle":"A parameter-light calculation shows frozen disorder is not the dominant cause; phonons are.","key_machinery":"The load-bearing object is the nonperturbative effective-medium dispersion relation for complex polariton energies in a Gaussian disorder bath: $E_P(q)-E_C(q,m) = (\\Omega_R^2/4)\\int dE'\\,\\rho(E')/(E_P(q)-E')$, rewritten with the scaled complementary error function for numerical stability. Its real part supplies the renormalized dispersion whose derivative gives the group velocity, while its imaginary part gives the resonance-scattering lifetime and the momentum broadening $\\delta_q$ that marks the regime of validity. The weak-disorder limit of the same equation reduces to a quadratic dispersion with a disorder-enhanced Rabi splitting, yielding closed-form insight into why lower polaritons are redshifted and upper polaritons blueshifted with wave-vector-dependent magnitude. Every numerical result and the perturbative interpretation in the paper follows from this single equation.","core_discovery":"The central claim is that static energetic disorder renormalizes the polariton dispersion in a way that always lowers the group velocity on both branches, with the size of the effect controlled by the ratio of the inhomogeneous broadening $\\sigma$ to the Rabi splitting $\\Omega_R$ and by the exciton content of the mode. Solving the effective-medium dispersion relation with a Gaussian exciton density of states gives complex polariton energies $E_P(q)$; the paper identifies the renormalized group velocity as $v_g(q)=\\hbar^{-1}\\partial_q \\mathrm{Re}[E_P(q)]$ and the disorder-induced momentum broadening as $\\delta_q = \\mathrm{Im}[E_P(q)]/(\\hbar v_g(q))$, with the method valid while $\\delta_q/q<1$. A weak-disorder expansion shows the effective Rabi splitting grows to $\\tilde{\\Omega}_R \\approx \\Omega_R\\sqrt{1+\\sigma^2/(2[E_P^{(0)}(q)-E_M]^2)}$, pushing the lower polariton down and the upper polariton up with wave-vector-dependent strength, which is why both group velocities drop. Applied to perovskite microcavities and to BODIPY exciton-surface-wave polaritons under the parameters used in the recent experiments, the predicted slowdown is small except at high exciton content and large $\\sigma/\\Omega_R$, leading the paper to attribute most of the experimentally observed slowdown to dynamical, phonon-assisted disorder.","pith_inferences":["The paper's largest predicted slowdowns sit in the $\\delta_q/q \\ge 1$ regime where the group-velocity formula is declared invalid, so the true transport speed in those modes may differ from the derivative of the renormalized dispersion; a wavepacket simulation on a disordered lattice would settle the size and direction of that difference.","The 10%-renormalization boundary shown in the paper, approximately $P_M(q) \\propto (\\sigma/\\Omega_R)^{-1}$, can be read as a design rule: for a given disorder strength, keep the exciton content below that line to keep static-disorder velocity loss under 10%.","If static disorder is genuinely minor in typical samples, existing datasets should show polariton velocities correlating with temperature and vibronic coupling strength but not with sample-to-sample inhomogeneous linewidth; re-analyzing published ultrafast microscopy data along those lines is a direct extension."],"forward_implications":["Static disorder can be set aside as the dominant explanation for the polariton slowdown observed in perovskite and BODIPY surface-wave microcavities; phonon-assisted scattering is the more likely cause.","The renormalization is selective: lower-polariton velocities drop most at large in-plane wave vectors where exciton content is high, upper-polariton velocities drop most near zero wave vector, and both effects become significant only as $\\sigma$ approaches $\\Omega_R/2$.","The method offers a parameter-only route to predict disorder-renormalized group velocities from independently measurable quantities, so the same comparison can be extended to any new strongly coupled microcavity material.","At low temperature, where phonon populations are frozen out, static disorder should become the dominant slowdown mechanism in materials with weak vibronic coupling and strong inhomogeneous broadening."],"supporting_citations":[{"why":"It supplies the effective-medium dispersion relation (equation 1) that the entire nonperturbative calculation solves.","marker":"50"},{"why":"It provides the BODIPY surface-wave polariton parameters and the experimental slowdown data used as a comparison target.","marker":"34"},{"why":"It provides the perovskite microcavity parameters, the temperature-dependent group-velocity data, and the main experimental comparison.","marker":"35"},{"why":"It gives the experimental estimate of perovskite exciton inhomogeneous broadening ($\\sigma = 5$–$10$ meV) that sets the realistic relative disorder strengths.","marker":"57"},{"why":"It supports the spectral-function validity condition $\\delta_q/q<1$ under which the group-velocity formula applies.","marker":"56"},{"why":"It articulates the phonon-assisted disorder mechanism that the paper concludes dominates the observed slowdown.","marker":"39"}],"fun_headline_variants":["Static disorder only minor player in polariton slowdown","Phonons, not static disorder, drive polariton slowdown","Frozen disorder slows polaritons, but phonons matter more","Polariton slowdown: static disorder not the main cause"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The method assumes that a polariton still travels with a well-defined direction and speed, so that its velocity is just the slope of the energy curve; the paper's own validity check says this stops holding exactly in the strongly excitonic regime where the biggest slowdowns are predicted.","fun_headline_variants_meta":{"raw":{"variants":["Static disorder only minor player in polariton slowdown","Phonons, not static disorder, drive polariton slowdown","Frozen disorder slows polaritons, but phonons matter more","Polariton slowdown: static disorder not the main cause"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.00064,"raw_usage":{"total_tokens":2976,"prompt_tokens":1001,"completion_tokens":1975,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":617,"completion_tokens_details":{"reasoning_tokens":1907}},"tokens_in":617,"tokens_out":1975,"duration_ms":18427,"temperature":1.0,"reasoning_tokens":1907,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-06T21:04:07.723695+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Simulate the real-space motion of a polariton wavepacket in a disordered microcavity by solving the full disordered microscopic model for the same $\\sigma/\\Omega_R$ values and compare its centroid velocity with $\\hbar^{-1}\\partial_q \\mathrm{Re}[E_P(q)]$; a mismatch in the $\\delta_q/q<1$ regime would invalidate the central method. On the experimental side, measuring low-temperature group velocities in two samples with identical Rabi splitting but different inhomogeneous broadening would directly test whether static disorder controls the slowdown.","supporting_citations":[{"cited_title":"Litinskaya \\ and\\ author P","cited_arxiv_id":null,"evidence_quote":"It supplies the effective-medium dispersion relation (equation 1) that the entire nonperturbative calculation solves."},{"cited_title":"Balasubrahmaniyam , author A","cited_arxiv_id":null,"evidence_quote":"It provides the BODIPY surface-wave polariton parameters and the experimental slowdown data used as a comparison target."},{"cited_title":"Xu , author A","cited_arxiv_id":null,"evidence_quote":"It provides the perovskite microcavity parameters, the temperature-dependent group-velocity data, and the main experimental comparison."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"It gives the experimental estimate of perovskite exciton inhomogeneous broadening ($\\sigma = 5$–$10$ meV) that sets the realistic relative disorder strengths."},{"cited_title":"Suyabatmaz \\ and\\ author R","cited_arxiv_id":null,"evidence_quote":"It supports the spectral-function validity condition $\\delta_q/q<1$ under which the group-velocity formula applies."},{"cited_title":"Sokolovskii , author Y","cited_arxiv_id":null,"evidence_quote":"It articulates the phonon-assisted disorder mechanism that the paper concludes dominates the observed slowdown."}],"review_version":1}