{"id":"4496557b-622a-40d1-bd4b-db68f90705c6","arxiv_id":"1909.01002","paper_version":2,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":7.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":0,"one_line_summary":"The paper derives the exact distribution of the Wigner-Smith time-delay matrix for a chaotic cavity with arbitrary absorption and N channels, reducing large-N statistics to two coupled Coulomb gases.","lead":"This mathematical physics paper derives the full probability distribution of the Wigner-Smith time-delay matrix for waves scattering in a chaotic cavity with absorption. It gives exact formulas for any absorption strength and shows that large-N statistics follow from two interacting Coulomb gases.","discovery_kind":"extension","skeptic_critique":{"model":"deepseek-v4-flash","headline":"Abstract overstates beta=1 scope: Eq. (56) is derived only for even N in the orthogonal class, so the claimed joint eigenvalue distribution for all three symmetry classes is not established for odd N.","rationale":"I read the paper as a technical advance: for beta=2, the full-matrix distribution (47) is derived from the fictitious-channel model with several internal consistency checks (weak- and strong-absorption limits, N=1 reduction), and the beta=1,4 eigenvalue formula (56) reduces correctly at gamma=0. The mathematical steps I checked — the unitary-group integral reduction, the N_phi→∞ limit, and the change of variables to Gamma — are coherent. The single most load-bearing weakness is not an algebraic error but a scope mismatch: the central claim is advertised for all three symmetry classes, while the beta=1 derivation is inherently even-N because of the half-integer integral dimension. This does not invalidate the even-N result, but it means the paper does not fully deliver the claimed generality. The reader's conditional verdict is appropriate; the authors should either supply the odd-N beta=1 derivation or explicitly restrict the abstract and main summary. I do not see a more serious internal inconsistency, and I do not regard the fictitious-channel model itself as a load-bearing weakness, since it is a standard and internally consistent modeling step.","tokens_in":21283,"tokens_out":12129,"duration_ms":103552,"concrete_test":"Compute the N=1, beta=1 distribution of Gamma directly from the known exact single-channel absorbing-cavity result (Savin and Sommers, Phys. Rev. E 68, 036211 (2003)) and attempt to cast it in a form analogous to Eq. (56). If no representation with N_t=1/2 can be defined, the odd-N beta=1 gap is genuine and the abstract must be amended. Alternatively, check whether the input Eq. (53) has an odd-N extension in Refs. [39,67]; if it does not, the restriction must be stated in the abstract and the summary.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The paper's advertised central result — the joint eigenvalue distribution of the Wigner-Smith matrix in absorbing cavities for beta=1, 2 and 4 — is not proven for beta=1 with odd N. The derivation of Eq. (56) in Section 2.2 starts from the finite-N reflection eigenvalue distribution Eq. (53) of Refs. [39,67]. That input has an integral of dimension N_t = beta N/2; for beta=1 and odd N, N_t is half-integer, so the representation is not defined. The paper acknowledges this in Section 2.2 ('the dimension N_t = N/2 of the integral in (53) restricts the number of channels to even numbers') and in the Conclusion, but the abstract and the summary of main results in Section 1.1 state without qualification that the joint distribution is derived in the other symmetry classes. The even-N computation may be correct, but the claim as advertised is narrower than proven. The gap is load-bearing because beta=1 odd N includes the physically relevant single-channel orthogonal case N=1, which is covered by none of the paper's new formulas.","agreement_with_reader":"partial"},"referee_report":{"model":"deepseek-v4-flash","summary":"The paper addresses the joint distribution of the Wigner-Smith time-delay matrix Q for a chaotic cavity with uniform absorption rate γ, perfectly coupled to N channels, within random matrix theory. Working with Γ = (NQ)^{-1}, the author derives for β = 2 an exact matrix-integral representation (Eq. (47)) and for β = 1, 2, 4 a joint eigenvalue distribution (Eq. (56)) by taking a fictitious-channel limit of the known reflection-eigenvalue distribution (Eq. (53)). The large-N analysis is recast as two interacting Coulomb gases, and saddle-point methods give the cumulants of the Wigner time delay in the weak- and strong-absorption limits (Eqs. (5), (6), (85)). The paper checks the γ → 0 Wishart-Laguerre limit, the strong-absorption limit, and the known N = 1 result.","tokens_in":21414,"tokens_out":13415,"duration_ms":124347,"significance":"If the derivations are correct, this is a significant contribution: it provides the first exact distribution of the full Wigner-Smith matrix for an absorbing chaotic cavity at arbitrary absorption, generalizing the zero-absorption result of Brouwer, Frahm and Beenakker. The β = 2 formula (47) is compact, and the Coulomb-gas treatment yields concrete predictions for cumulants of the Wigner time delay. The derivation is transparent and the consistency checks are convincing. The main caveat is the scope for β = 1: Eq. (56) is established only for even N because the input Eq. (53) has half-integer dimension for odd N; the abstract and Section 1.1 currently overstate the coverage.","major_comments":[{"comment":"The abstract and Section 1.1 claim, without qualification, that the joint eigenvalue distribution is derived for β = 1, 2, and 4. However, Eq. (56) is obtained from Eq. (53), whose integral dimension N_t = βN/2 is half-integer for β = 1 and odd N; the author acknowledges in Section 2.2 that the dimension N_t = N/2 restricts the number of channels to even numbers, and the Conclusion states that the odd-N β = 1 case remains open. Since β = 1 odd N includes the physically relevant single-channel case N = 1, the advertised claim is narrower than what is proven. Please qualify the abstract and Section 1.1 (for example, by stating that for β = 1 the result holds for even N) and explicitly note that odd-N β = 1 remains open.","section":"Abstract; §1.1; §2.2"}],"minor_comments":[{"comment":"The product over i < j in Eq. (50) appears to be a typographical error: it should read ∏_{i<j}(t_i - t_j)^2 ∏_{i,j=1}^N (Γ_i - t_j), rather than ∏_{i<j}((t_i - t_j)^2(Γ_i - t_j)), which omits most of the Γ-t cross terms and is inconsistent with the Andréief form in Eq. (51).","section":"Eq. (50)"},{"comment":"Eq. (4) should carry the same even-N restriction for β = 1 as Eq. (56), or a footnote should point to Section 2.2 explaining that the integral dimension N_t = N/2 restricts the formula to even N in the orthogonal class.","section":"§1.1, Eq. (4)"},{"comment":"There is a typo in the phrase introducing Eq. (3): 'Kroenecker' should be 'Kronecker'.","section":"§1.1, below Eq. (3)"},{"comment":"The index convention is inconsistent in Eqs. (49)-(50): the eigenvalues of Γ are denoted Γ_n in the text but Γ_i appears in the integral over t; please make the notation uniform.","section":"Eqs. (49)-(50)"},{"comment":"The sentence in the Conclusion that the derivation covers 'most situations' is vague; please specify there that the β = 1 odd-channel case is excluded, matching the caveat in Section 2.2.","section":"Conclusion"}],"recommendation":"minor_revision","confidential_remarks":"The paper is technically sound in the parts I checked, and the body already contains the correct even-N caveat for β = 1. The main issue is the unqualified claim in the abstract and Section 1.1; a revised abstract and summary should resolve it. I would not require the author to solve the open odd-N β = 1 problem for acceptance."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"Short version: this paper fills a real gap. The distribution of the Wigner-Smith matrix for a chaotic cavity with absorption was previously known only in limits or for small channel numbers; here it is given exactly for any absorption rate. For beta=2, Eq. (47) is a compact matrix integral for the full matrix; for beta=1 and 4 (and beta=2 as a check) Eq. (56) gives the joint eigenvalue distribution. The derivation from the Poisson kernel is a clean sequence of unitary matrix integrals, and it reproduces the zero-absorption Wishart-Laguerre limit, the strong-absorption limit, and the known N=1 case. The two-gas Coulomb gas representation, with different couplings inside and between the gases, is new and gives a sensible handle on the large-N problem.\n\nI agree with the stress-test note: the advertised beta=1 coverage is overstated. Eq. (56) is derived from the reflection-eigenvalue distribution of Refs. [39,67], which is only defined for even N in the orthogonal class. The paper discloses this in Section 2.2 and in the conclusion, but the abstract and the Section 1.1 summary state without qualification that the joint distribution is derived in all three symmetry classes. Since beta=1 with odd N includes the single-channel orthogonal case, the mismatch is not cosmetic. This is a fixable presentation problem, not a fatal flaw: the even-N computation is what the math supports, and the author should say so up front.\n\nThe other soft spot is that the large-N cumulant expansions are not numerically benchmarked. They are consistent with the known mean in both limits, and the saddle-point setup is standard, but a single numerical check would raise confidence. Correctness risk is moderate, not high.\n\nThe citation pattern is honest: the paper leans on [39,67] for the input and on the author's earlier thermodynamic identity for the Coulomb gas step. No circularity. The self-citations are appropriate.\n\nWho should read this: random matrix theorists working on scattering, and experimentalists extracting Wigner-Smith time delays from microwave or coherent-transport data. It is a solid technical advance within an established framework, not a change of worldview. Send it to a serious referee. Ask the referee to push for the abstract to match the even-N beta=1 scope, and suggest a numerical check if feasible.","headline":"Solid exact RMT result for the Wigner-Smith matrix with absorption, but the beta=1 coverage is proven only for even N and the abstract should say so.","tokens_in":22027,"tokens_out":3057,"would_cite":true,"duration_ms":29584,"reading_group":"yes","serious_thinker":"yes","would_accept_peer_review":true},"rs_alignment":null,"lean_confirmation":null,"pith_extraction":{"msc":["15B52","60B20","81Q50"],"pacs":[],"model":"deepseek-v4-flash","headline":"This paper derives the exact distribution of the Wigner-Smith time-delay matrix for chaotic cavities with uniform absorption, covering all symmetry classes and any absorption rate.","keywords":["Wigner-Smith matrix","time delay","chaotic cavity","random matrix theory","absorption","Coulomb gas","Wishart-Laguerre distribution"],"falsifier":"Simulate the scattering matrix of an absorbing chaotic cavity directly, for instance by drawing the unitary matrix from the Poisson kernel while including fictitious channels (or by an imaginary energy shift), and compare the empirical distribution of $\\Gamma = (NQ)^{-1}$ for small $N$ (say $N=2$, $\\beta=2$) with the integrated form of Eq. (56); a mismatch in the eigenvalue density would falsify the claimed exact distribution.","tokens_in":20982,"feed_emoji":"⏱️","tokens_out":8515,"duration_ms":72615,"temperature":0.7,"pith_summary":"This paper derives the full statistical distribution of the Wigner-Smith time-delay matrix, a Hermitian matrix that encodes how long waves spend inside a scattering cavity, when the cavity absorbs energy uniformly. The result covers the three symmetry classes (orthogonal, unitary, symplectic) and any absorption rate, and reduces to the known Wishart-Laguerre distribution in the zero-absorption limit. For broken time-reversal symmetry the distribution takes the compact form of a matrix integral; in the other classes the eigenvalues are governed by a joint distribution with an auxiliary integral. The paper also shows that, for many channels, the statistics are described by two coupled Coulomb gases, and uses this to compute the first cumulants of the Wigner time delay in the weak- and strong-absorption regimes.","feed_headline":"Exact distribution found for time delays in absorbing chaotic cavities","feed_subtitle":"A matrix integral captures all absorption rates, generalizing the zero-loss Wishart-Laguerre law.","key_machinery":"The argument is carried by two ingredients. First, absorption is modelled by $N_\\phi$ fictitious channels with tunnel probability $T$, in the limit $N_\\phi \\to \\infty$ with $N_\\phi T = \\gamma N$; this converts the absorbing cavity into a larger unitary scattering problem governed by the Poisson kernel. Second, for $\\beta = 2$ the distribution of the reflection matrix is obtained by explicit integration over the unitary group, while for general $\\beta$ the eigenvalue distribution follows from the known finite-$N$ reflection eigenvalue distribution of Jarosz, Vidal and Kanzieper. A change of variables $\\Gamma = \\gamma (1 - r^\\dagger r)^{-1}$ then yields the time-delay distribution. For large $N$, the joint distribution is reinterpreted as two logarithmically repelling Coulomb gases on $[0,\\gamma]$ and $[\\gamma,\\infty)$, with the inter-gas repulsion half as strong as the intra-gas repulsion, whose saddle-point equations are solved with Tricomi's theorem in the weak- and strong-absorption limits.","core_discovery":"The central discovery is that, for a chaotic cavity with uniform absorption rate $\\gamma$ and $N$ perfectly coupled channels, the inverse time-delay matrix $\\Gamma = (NQ)^{-1}$ has the distribution, for $\\beta = 2$, $P(\\Gamma) \\propto e^{-N \\operatorname{tr}\\Gamma} \\int_0^{\\gamma 1_N} dT \\det(1_N \\otimes \\Gamma - T \\otimes 1_N) e^{-N \\operatorname{tr}T}$, with $\\Gamma > \\gamma 1_N$. For the other symmetry classes, the joint distribution of the eigenvalues is given by an integral over $\\beta N/2$ auxiliary variables $t_n \\in [0,\\gamma]$. These expressions extend the zero-absorption Wishart-Laguerre law of Brouwer, Frahm and Beenakker to arbitrary absorption, and imply that the eigenvalues of $\\Gamma$ are bounded below by $\\gamma$, so the proper time delays are bounded above by $1/(N\\gamma)$.","pith_inferences":["The edge constraint $\\Gamma > \\gamma$ suggests that the density of eigenvalues of $\\Gamma$ develops a universal edge behaviour at $\\gamma$; checking this in a microwave-cavity experiment would be a direct test of the distribution.","For the orthogonal class with odd $N$, the paper leaves the distribution open; an inference is that a Pfaffian or a modified auxiliary ensemble may be needed, and the even-$N$ duality hints at a closed form waiting to be found.","The two-gas representation with unequal inter-gas repulsion may transfer to other random-matrix settings where two coupled spectral sets interact, such as eigenvalues above a threshold in Wishart matrices, with modified exponents; this could be tested numerically."],"forward_implications":["The full distribution of the time-delay matrix (or its eigenvalues) is now known for any absorption rate, so all linear statistics of the proper time delays become computable at finite $N$.","The hard bound $\\Gamma > \\gamma 1_N$ means proper time delays cannot exceed $1/(N\\gamma)$, quantifying how absorption eliminates long-lived resonances.","The Coulomb-gas description yields a large-deviation function for the Wigner time delay, producing the mean and variance in the weak- and strong-absorption limits and matching the known mean $\\langle \\tau_W \\rangle = 1/(N(\\gamma+1))$.","In the weak-absorption regime, all cumulants of the Wigner time delay can be expressed in terms of the cumulants at zero absorption via Eq. (85), extending earlier zero-absorption results."],"supporting_citations":[{"why":"Gives the zero-absorption Wishart-Laguerre distribution of the inverse time-delay matrix that this paper generalizes to γ > 0.","marker":"[10,11]"},{"why":"Supplies the finite-N reflection-eigenvalue distribution (Eq. (53)) used as the starting point for the β = 1 and 4 results.","marker":"[39,67]"},{"why":"Provides the fictitious-channel model of absorption and the change of variables on the unitary group used in the β = 2 derivation.","marker":"[8]"},{"why":"Establishes the relation r†r = 1 - γ N Q and gives the N = 1 distribution used to check the new result.","marker":"[59]"},{"why":"Gives the large-N mean density of reflection eigenvalues that the Coulomb-gas typical density must reproduce.","marker":"[60]"},{"why":"Provides the unitary-group integral formula needed to evaluate the β = 2 reflection-matrix distribution.","marker":"[28]"},{"why":"Derives the Wigner time-delay distribution at zero absorption whose cumulant structure is extended to weak absorption.","marker":"[64]"}],"fun_headline_variants":["Exact time-delay law found for absorbing chaotic cavities","Time delays in absorbing cavities: exact distribution","Matrix integral solves time-delay distribution with absorption","Coulomb gas approach yields exact time-delay law","Beyond Wishart-Laguerre: time delays with absorption"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The derivation assumes that a cavity with $N_\\phi$ fictitious absorbing channels of tunnel probability $T$ reproduces a uniform absorption rate $\\gamma$ in the limit $N_\\phi \\to \\infty$ with $N_\\phi T = \\gamma N$; if that equivalence fails, the exact distributions do not follow.","fun_headline_variants_meta":{"raw":{"variants":["Exact time-delay law found for absorbing chaotic cavities","Time delays in absorbing cavities: exact distribution","Matrix integral solves time-delay distribution with absorption","Coulomb gas approach yields exact time-delay law","Beyond Wishart-Laguerre: time delays with absorption"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000437,"raw_usage":{"total_tokens":2191,"prompt_tokens":884,"completion_tokens":1307,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":500,"completion_tokens_details":{"reasoning_tokens":1231}},"tokens_in":500,"tokens_out":1307,"duration_ms":11175,"temperature":1.0,"reasoning_tokens":1231,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-14T05:30:09.690200+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Simulate the scattering matrix of an absorbing chaotic cavity directly, for instance by drawing the unitary matrix from the Poisson kernel while including fictitious channels (or by an imaginary energy shift), and compare the empirical distribution of $\\Gamma = (NQ)^{-1}$ for small $N$ (say $N=2$, $\\beta=2$) with the integrated form of Eq. (56); a mismatch in the eigenvalue density would falsify the claimed exact distribution.","supporting_citations":[{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Provides the fictitious-channel model of absorption and the change of variables on the unitary group used in the β = 2 derivation."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Establishes the relation r†r = 1 - γ N Q and gives the N = 1 distribution used to check the new result."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Gives the large-N mean density of reflection eigenvalues that the Coulomb-gas typical density must reproduce."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Provides the unitary-group integral formula needed to evaluate the β = 2 reflection-matrix distribution."},{"cited_title":"Texier and S","cited_arxiv_id":null,"evidence_quote":"Derives the Wigner time-delay distribution at zero absorption whose cumulant structure is extended to weak absorption."}],"review_version":1}