{"id":"cb855d3b-8c66-44f9-bd63-342804f9ffab","arxiv_id":"1908.01708","paper_version":1,"verdict":"CONDITIONAL","confidence":"HIGH","novelty_score":6.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":0,"one_line_summary":"The optical temperature and chemical potential of a thermalized nonlinear multimode system are uniquely fixed by the conserved power and energy, with the sign of the temperature set by whether the normalized energy lies above or below the average spectrum.","lead":"This paper proves that for any given multimode optical system, the optical temperature and chemical potential at thermal equilibrium are uniquely determined by the input power and energy. It then gives simple conditions on the average spectrum that tell whether the equilibrium temperature is positive or negative.","discovery_kind":"extension","skeptic_critique":{"model":"deepseek-v4-flash","headline":"Uniqueness is proven only within the two-constraint Rayleigh-Jeans ensemble; the paper's universal claim that any weakly nonlinear multimode system reaches this equilibrium is unsupported because systems with extra conserved quantities can have different equilibria.","rationale":"The paper's mathematical derivation is a clean proof of uniqueness of the temperature for a fixed linear spectrum, power, and internal energy, under the explicit assumption that the equilibrium occupancies follow the Rayleigh-Jeans distribution with only two conserved quantities. The reader's weakest_assumption identifies exactly this point, and my stress-test concurs: the proof is conditional, while the abstract and conclusion claim universal validity for 'any' weakly nonlinear multimode system. The strongest reason to doubt the universal claim is the existence of additional conserved quantities in symmetric systems. A circular fiber conserves angular momentum, and a periodic lattice conserves quasimomentum; entropy maximization with those extra constraints yields a distribution with an additional chemical potential, so the two-parameter uniqueness theorem does not apply. The paper provides no argument that such symmetries are either absent or ineffective, and its numerical evidence covers only a single non-symmetric lattice geometry. This is not an internal inconsistency in the proof, but a scope-overreach in the central claim as advertised. The verdict should remain CONDITIONAL: the theorem is correct within its stated assumptions, but the paper must temper the universality claim or extend the proof to systems with additional conserved quantities. A secondary, smaller issue is the S=0 transition point where the polynomial degree drops and the physical root is at infinite temperature; Table 1 mentions the limit, but the proof's degree-N-1 statement is not valid there. This does not alter the overall verdict.","tokens_in":7071,"tokens_out":16978,"duration_ms":175475,"concrete_test":"Run the DNLS simulation of Eq. (5) on a 20x20 lattice with periodic boundary conditions so that total quasimomentum is conserved, using an input concentrated on modes with a single nonzero Bloch wavevector; after thermalization, fit the ensemble-averaged occupancies to the two-parameter RJ form |c_i|^2 = -T/(E_i+mu) and compare with the Eq. (2) prediction from the conserved P and U. If a third chemical potential (for quasimomentum) is needed for a statistically acceptable fit, the universal uniqueness claim is refuted.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The central theorem (Eq. 2 and the root-counting argument after it) is internally sound: for a fixed spectrum, power P, and linear energy U, the function f(x) has exactly one root outside (rho_1, rho_N) when the equilibrium occupancies are exactly the Rayleigh-Jeans form |c_i|^2 = -T/(E_i+mu). The load-bearing assumption is the step 'the thermalized average power levels conveyed by each mode are found to obey a Rayleigh-Jeans distribution,' which follows from maximizing entropy with only P and U as constraints. The abstract and conclusion elevate this to a universal statement: 'any complex weakly nonlinear highly multimode system can reach thermodynamic equilibrium that is characterized by a unique temperature and chemical potential.' That universality is not established. A weakly nonlinear system with a continuous symmetry has an additional conserved quantity: a circular multimode fiber conserves angular momentum, and a periodic lattice conserves quasimomentum. For such systems the equilibrium distribution obtained by entropy maximization carries an additional Lagrange multiplier, e.g. |c_i|^2 = -T/(E_i + mu + nu L_i), so the temperature is no longer uniquely determined by (P,U) alone, and Eq. (2) does not apply. The paper gives no argument that these invariants are absent or irrelevant, and the numerical support (Fig. 3) is confined to one finite, anisotropic, boundary-broken 20x20 lattice with four initial conditions. Consequently, the proof establishes uniqueness conditional on reaching the two-constraint RJ equilibrium, but not the abstract's claim that every weakly nonlinear multimode system reaches that equilibrium. A secondary technical gap: the statement that f(x) is a degree N-1 polynomial fails when S = sum g_i rho_i = 0 (the -U/P = E_bar transition); there the finite root moves to T = ±infinity, which the proof does not cover except as a limit in Table 1.","agreement_with_reader":"agree"},"referee_report":{"model":"deepseek-v4-flash","summary":"The paper considers weakly nonlinear multimode optical systems with M modes, a linear spectrum E_i, and two conserved quantities: total power P and linear internal energy U. Assuming the equilibrium modal occupancies are exactly the Rayleigh-Jeans distribution |c_i|^2 = -T/(E_i + mu), the authors rewrite the power constraint as a rational equation f(x)=0 in the variable x = M T (Eq. (2)). They prove that f(x) is equivalent to a degree N-1 polynomial, that it has N-2 real roots in the intervals between distinct eigenvalue levels, and that exactly one additional real root lies outside the interval (rho_1, rho_N); this exterior root is the only one giving positive modal occupancies and hence uniquely determines the temperature T and chemical potential mu for fixed (P, U, spectrum). Table 1 gives explicit inequalities for positive versus negative temperature, and the paper tests the predicted equilibrium distributions against numerical simulations in a 20x20 nonlinear waveguide array (Fig. 3). The abstract and conclusion claim that the results are universal for any weakly nonlinear highly multimode optical system.","tokens_in":7356,"tokens_out":10676,"duration_ms":110240,"significance":"The root-counting argument is mathematically sound and gives a concise, parameter-free proof of uniqueness of the temperature inside the two-constraint Rayleigh-Jeans ensemble. The explicit conditions in Table 1, and the associated physical interpretation in terms of positive and negative temperatures, are genuinely useful for interpreting beam self-cleaning and related multimode thermalization experiments. However, the central theorem is conditional on an equilibrium distribution that maximizes entropy subject only to P and U; the manuscript does not establish the universal reachability claim in its abstract and conclusion. If the claims are properly restricted, the paper is a solid contribution to the thermodynamics of multimode photonic systems.","major_comments":[{"comment":"The uniqueness proof assumes that the equilibrium occupancies are exactly the two-constraint Rayleigh-Jeans form |c_i|^2 = -T/(E_i + mu), obtained by maximizing entropy with only P and U conserved. The paper does not prove that every weakly nonlinear multimode system thermalizes to this ensemble; the conclusion that the results 'universally apply to any optical nonlinear multimode system' is therefore not supported. Systems with additional conserved quantities (for example, angular momentum in a circular multimode fiber or quasimomentum in a periodic lattice) would have equilibrium distributions with an extra Lagrange multiplier, such as |c_i|^2 = -T/(E_i + mu + nu L_i), for which Eq. (2) and the uniqueness conclusion do not follow. The numerical evidence in Fig. 3 is restricted to a single finite anisotropic 20x20 lattice whose boundary conditions break translation invariance, so it does not test the universal claim. I recommend restricting the claims to systems whose thermodynamic equilibrium is described by the two-constraint Rayleigh-Jeans distribution, or supplying a substantive argument that no other conserved quantity constrains the equilibrium.","section":"Abstract, Eq. (2), and Conclusion"},{"comment":"The exterior-root argument assumes sum g_i rho_i is nonzero. When -U/P = Ebar exactly, sum g_i rho_i = 0, and the asymptotics in Eq. (4) give f(x) -> 0 at both infinities, so there is no finite exterior root; the only physical solution is the infinite-temperature state. The text mentions T -> +/-infinity only as a limit when approaching the transition point. The uniqueness statement should explicitly address this measure-zero case, for example by stating that T = +/-infinity is the unique physical solution, rather than leaving it as a limiting statement.","section":"Table 1 and paragraph after Eq. (4)"}],"minor_comments":[{"comment":"The number of random-phase realizations used for the ensemble averages in Fig. 3 is not stated; the authors should provide this number and an estimate of the statistical spread of the reported average occupancies.","section":"Fig. 3 and numerical methods"},{"comment":"The sentence 'If in other hand the quantity -U/P approaches from above or below the transition point Ebar' contains a typo ('If in other hand' should be 'If, on the other hand'); the equality case -U/P = Ebar should also be listed in Table 1.","section":"Discussion after Table 1"},{"comment":"The definition U = -<Psi|H_L|Psi> and the ordering E_1 <= ... <= E_M are specific to the waveguide convention; the authors note that the ordering is reversed in cavity arrangements but do not give the corresponding form of Eq. (2) or state whether Table 1 applies unchanged under that convention.","section":"Definitions of U and eigenvalue ordering"},{"comment":"Reference [15] is cited as 'in press'; the authors should update the citation with the published details or provide a preprint identifier.","section":"References"},{"comment":"The statement that 'the system is expected to reach thermal equilibrium by maximizing its entropy' omits the standard caveat that this holds in the weak-nonlinearity, random-phase approximation; the authors should explicitly cite the conditions from refs. [16,17] under which the Rayleigh-Jeans distribution is the correct equilibrium.","section":"Entropy maximization statement"}],"recommendation":"major_revision","confidential_remarks":"The mathematical core of the paper is sound and likely publishable after modification. The main concern is the gap between the conditional theorem (unique T given P, U, and a two-constraint RJ equilibrium) and the universal claims in the abstract and conclusion. The authors should be asked to either narrow those claims or provide evidence that additional conserved quantities do not alter the equilibrium in the systems they target. A secondary concern is that the paper is written as a letter for a photonics venue, but the primary contribution is a mathematical uniqueness proof; the authors should make the scope and conditions of that proof very explicit for the intended readership."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"Read this as a conditional uniqueness result, not as the universal thermalization theorem the abstract advertises. For a fixed spectrum, fixed P and U, and assuming the two-constraint Rayleigh-Jeans equilibrium, the algebra is right: f(x) reduces to a degree N-1 polynomial when S=sum g_i rho_i is nonzero, Bolzano gives N-2 interior roots, and the asymptotics plus positivity just outside (rho_1,rho_N) force exactly one exterior physical root. Table 1 is a useful, explicit set of conditions for positive/negative T, and the four numerical examples in the 20x20 lattice match the predicted distributions well. That is a solid extension of the earlier thermodynamic multimode program, and it earns a serious referee.\n\nThe soft spot is the scope of the claim. The derivation simply assumes the equilibrium occupancies are |c_i|^2 = -T/(E_i+mu), which follows from maximizing entropy with only P and U conserved. Any additional conserved quantity—angular momentum in a circular fiber, quasimomentum in a periodic lattice—adds a Lagrange multiplier and moves the occupancy away from that simple form. Eq. (2) then has no reason to describe the equilibrium. The paper offers no argument that such invariants are absent or negligible, and the numerics, one anisotropic open-boundary lattice, cannot test the universal language because that lattice has no continuous symmetry. The conclusion's 'any optical nonlinear multimode system' is the part I would ask the authors to retire unless they add a proper treatment of extra invariants.\n\nTwo minor notes. At the S=0 transition the physical root goes to infinity; the finite-root proof silently assumes S != 0, though Table 1 covers the transition as a limit. And the equation of state is taken from the authors' own earlier paper; that is a dependence, not circularity, since they are extending that formalism rather than fitting it.\n\nVerdict: this is a competent, useful paper for people working on multimode thermodynamics, beam self-cleaning, and multimode lasers. It should go to review, with a clear request to narrow the claims.","headline":"A clean conditional uniqueness theorem with an overbroad universality claim; referee it, but push the authors to either prove or drop the 'any system' language.","tokens_in":7997,"tokens_out":4457,"would_cite":true,"duration_ms":46965,"reading_group":"yes","serious_thinker":"yes","would_accept_peer_review":true},"rs_alignment":null,"lean_confirmation":null,"pith_extraction":{"msc":[],"pacs":[],"model":"deepseek-v4-flash","headline":"This paper proves that the optical temperature and chemical potential of a thermalized multimode system are uniquely set by the spectrum, power, and energy.","keywords":["optical thermodynamics","Rayleigh-Jeans distribution","negative temperature","multimode nonlinear optics","beam self-cleaning","chemical potential","mode thermalization"],"falsifier":"Take a fixed multimode array with known spectrum and launch two different input distributions that have identical total power $P$ and internal energy $U$; if the equilibrium modal occupancies measured after long propagation differ measurably from each other or from the predicted Rayleigh-Jeans form, the uniqueness claim fails. More sharply, in a translationally invariant array, manufacturing initial conditions with the same $P$ and $U$ but different conserved momentum would reveal whether an additional invariant controls the final state.","tokens_in":6866,"feed_emoji":"🌡️","tokens_out":4806,"duration_ms":46088,"temperature":0.7,"pith_summary":"This paper proves that a weakly nonlinear multimode optical system that reaches thermodynamic equilibrium has one and only one physically meaningful temperature and one chemical potential for a given linear spectrum, total power, and internal energy. The physically meaningful requirement is that every mode carries positive power, which forces the temperature root to lie outside the interval set by the spectrum and energy. The paper gives explicit conditions: temperature is positive when the ratio $-U/P$ exceeds the average spectral value, negative when it falls below, and diverges through infinity exactly at the average. It demonstrates the result in a $20 \\times 20$ nonlinear waveguide array, where numerically thermalized modal distributions match the predicted Rayleigh-Jeans form and temperature to within a few percent.","feed_headline":"Every heated multimode beam has one unique temperature","feed_subtitle":"Input power, internal energy, and the linear spectrum alone decide the optical temperature and chemical potential.","key_machinery":"The load-bearing object is the rational function $f(x) = \\sum_i g_i x/(x-\\rho_i) - M$, whose real roots are the candidate temperatures multiplied by $M$. Each pole $\\rho_i = U + P\\mathcal{E}_i$ combines the internal energy and the $i$-th eigenvalue, and $g_i$ is the degeneracy. The distribution formula $|c_i|^2 = -PT/(\\rho_i - x)$ shows that positive modal powers require the root $x$ to stay outside the whole interval $(\\rho_1, \\rho_N)$. Counting signs across the $N-2$ intervals between poles, then using the asymptotic behavior of $f$ at $\\pm\\infty$, proves that exactly one physically admissible root exists. A final one-line equation of state, $\\mu = (U - MT)/P$, then fixes the chemical potential.","core_discovery":"At thermal equilibrium, the modal powers obey a Rayleigh-Jeans distribution, $|c_i|^2 = -T/(\\mathcal{E}_i + \\mu)$. Substituting the conserved power $P$ and internal energy $U$ turns the normalization condition into an algebraic equation for $x = MT$ whose roots are all real: $f(x) = \\sum_i g_i x/(x-\\rho_i) - M = 0$, with $\\rho_i = U + P\\mathcal{E}_i$. Sign changes between consecutive $\\rho_i$ give $N-2$ roots by Bolzano's theorem, and the asymptotic sign of $f$ at $\\pm\\infty$ supplies one more root in one of the two outer intervals. Only one root lies outside $(\\rho_1, \\rho_N)$, and only that root keeps every modal occupancy positive; hence the temperature is unique. The corresponding chemical potential follows from the equation of state $U - \\mu P = MT$. The paper states the threshold explicitly: $T > 0$ when $-U/P > \\bar{\\mathcal{E}}$ and $T < 0$ when $-U/P < \\bar{\\mathcal{E}}$, with the domains $T > (U + P\\mathcal{E}_N)/M$ and $T < (U + P\\mathcal{E}_1)/M$ respectively.","pith_inferences":["An implication left implicit is that the uniqueness proof would break down in systems with additional conserved quantities (for example transverse momentum in a translationally invariant array): the equilibrium would then be set by a generalized Gibbs distribution, and the argument would need more constraints than just $P$ and $U$.","A testable extension would be to launch multiple input realizations that share the same $P$ and $U$ but differ in random phases; the paper's numerical ensembles already do this, so a natural experiment is to check whether the fitted $T$ is identical up to fluctuations for every such ensemble.","If the nonlinearity is increased until the linear Hamiltonian no longer dominates, the conserved quantity $U$ is no longer fixed and the Rayleigh-Jeans form is expected to be replaced by a different equilibrium; the uniqueness statement should be read as restricted to the weakly nonlinear regime."],"forward_implications":["For any multimode structure whose spectrum is known, an experimentalist can predict the equilibrium temperature and chemical potential directly from the input power and the energy set by the launch condition.","The beam self-cleaning effect, where power settles into the low-order modes, corresponds to positive optical temperatures; negative temperatures instead concentrate power in high-order modes.","At $-U/P = \\bar{\\mathcal{E}}$ the temperature diverges and all modes become equally populated; approaching the spectral edges makes $T$ tend to $0^+$ or $0^-$, the condensate limits.","Since negative temperatures lie above infinite temperature on the thermodynamic scale, energy flows from a negative-temperature system to a positive-temperature one, which determines the direction of mode redistribution.","The results apply to continuous and discrete multimode systems alike, including multimode fibers, cavities, waveguide arrays, and coupled-resonator structures."],"supporting_citations":[{"why":"Provides the Rayleigh-Jeans distribution $|c_i|^2 = -T/(\\mathcal{E}_i+\\mu)$ that the thermal equilibrium ansatz assumes.","marker":"[15-17, 20]"},{"why":"Supplies the equation of state $U - \\mu P = MT$ connecting the temperature to the conserved quantities.","marker":"[15]"},{"why":"Establishes the continuous and discrete multimode structures (waveguides, arrays, coupled resonators) to which the argument applies.","marker":"[18, 19]"},{"why":"Bolzano's theorem, used to count the $N-2$ real roots of $f(x)$ between consecutive poles.","marker":"[21]"},{"why":"Supports the interpretation of negative temperatures as hotter than positive temperatures, used in the paper's thermodynamic discussion.","marker":"[22]"}],"fun_headline_variants":["One temperature for every heated multimode beam","Input power and spectrum dictate optical temperature","Unique thermal equilibrium for nonlinear multimode systems","Positive or negative temperature? Power and spectrum say","Multimode light reaches a single thermal state"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The derivation assumes that the modal occupancies at equilibrium are exactly the Rayleigh-Jeans distribution obtained by maximizing entropy with only total power and total energy conserved; if a real system has extra conserved quantities or nonlinear corrections change the distribution, the uniqueness proof no longer applies.","fun_headline_variants_meta":{"raw":{"variants":["One temperature for every heated multimode beam","Input power and spectrum dictate optical temperature","Unique thermal equilibrium for nonlinear multimode systems","Positive or negative temperature? Power and spectrum say","Multimode light reaches a single thermal state"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000553,"raw_usage":{"total_tokens":2593,"prompt_tokens":863,"completion_tokens":1730,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":479,"completion_tokens_details":{"reasoning_tokens":1663}},"tokens_in":479,"tokens_out":1730,"duration_ms":12946,"temperature":1.0,"reasoning_tokens":1663,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-14T15:06:06.448811+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Take a fixed multimode array with known spectrum and launch two different input distributions that have identical total power $P$ and internal energy $U$; if the equilibrium modal occupancies measured after long propagation differ measurably from each other or from the predicted Rayleigh-Jeans form, the uniqueness claim fails. More sharply, in a translationally invariant array, manufacturing initial conditions with the same $P$ and $U$ but different conserved momentum would reveal whether an additional invariant controls the final state.","supporting_citations":[{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Supplies the equation of state $U - \\mu P = MT$ connecting the temperature to the conserved quantities."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Bolzano's theorem, used to count the $N-2$ real roots of $f(x)$ between consecutive poles."},{"cited_title":"Braun, J","cited_arxiv_id":null,"evidence_quote":"Supports the interpretation of negative temperatures as hotter than positive temperatures, used in the paper's thermodynamic discussion."}],"review_version":1}