{"id":"c8dcb926-5814-4143-87fd-0cd567149223","arxiv_id":"2505.21163","paper_version":1,"verdict":"REJECT","confidence":"MODERATE","novelty_score":6.0,"correctness_risk":"high","formal_verification":"none","parameter_count":1,"one_line_summary":"High laser power in a few-mode fiber inverts the modal population and drives the inferred optical temperature through infinity to negative values, an attractor the authors call spatial beam freezing.","lead":"In a few-mode optical fiber, raising the input laser power makes the output beam switch from the fundamental to the highest-order fiber modes, a regime called spatial beam freezing. The paper interprets this as a transition to negative absolute temperature, but that interpretation relies on a phenomenological model with a coefficient fitted to the data.","discovery_kind":"new_application","skeptic_critique":{"model":"deepseek-v4-flash","headline":"Fixed-UL calibration in Methods forces the temperature sign flip; since the measured modal inversion implies U_L(P) varies with power, the negative-temperature attractor is not established.","rationale":"The reader's weakest-assumption analysis identifies exactly the same load-bearing point: the fixed-U_L procedure is inconsistent with the observed power-dependent mode redistribution. My stress-test strengthens this by showing why the fixed-U_L step is not merely a quantitative inaccuracy but is likely what creates the sign flip in the model. With U_L/P constant (as expected for rescaled injections), the solution of Eq. (8) does not have the same U_L/P trajectory across the critical value; with U_L fixed, U_L/P decays as 1/P, forcing the T=±∞ divergence and the negative-temperature branch. The empirical observations of modal inversion and robustness are credible and independently supported by the cut-back and power-conservation checks, but they do not by themselves establish negative absolute temperature. The central theoretical claim depends on an equation of state that the authors admit is not derived from the conservation laws, and on a temperature inferred from a model whose key input is held constant in a way that contradicts the data. Therefore the reader's REJECT verdict is appropriate; I see no reason to change it.","tokens_in":10237,"tokens_out":13611,"duration_ms":157937,"concrete_test":"Re-analyze the Fig. 3a data: compute U_L(P) directly via Eq. (2) at every input power and test whether U_L(P)/P is constant. Then replace the fixed low-power U_L in Eq. (8) with the measured U_L(P) (or with U_L(P)=-PΣβ_i f_i from the input modal weights) and re-fit γ. If the recomputed T(P) no longer diverges and approaches 0− at high power, the negative-temperature attractor is an artifact of the fixed-U_L assumption; if the sign flip survives, the concern is refuted.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The most load-bearing condition for the central claim is in Methods, 'Determination of thermodynamic parameters': U_L is computed once from the mode occupancy at the lowest input power and 'kept as a constant for all power levels' when solving Eq. (8). This is not a harmless calibration. The paper's own observable is a strong power-induced redistribution from LP01 to LP21, so U_L(P)=-Σβ_i|c_i|^2 evaluated from the measured occupancies must change substantially with P. For a fixed input spatial profile rescaled in power, U_L(P)/P should remain near -β_1; holding U_L fixed instead makes U_L/P tend toward 0 as P grows, which drives U_L/P across the critical value -Σβ_i/M in Eq. (5) and produces the divergence of T and the subsequent T→0− branch in Eq. (8). The manuscript itself flags the proximate gap: Eq. (7) 'cannot be directly derived' from the RJ law and the conservation laws. The fixed-U_L step is where that unproved equation of state is converted into a sign flip. Because separate input-power runs are different initial conditions, U_L should scale with P, not remain constant. The empirical modal inversion may be real; the negative-absolute-temperature claim does not follow unless T is recomputed with a U_L(P) consistent with the measured distributions.","agreement_with_reader":"agree"},"referee_report":{"model":"deepseek-v4-flash","summary":"The paper reports experiments in a few-mode step-index fiber (SMF-28 operated below cutoff) in which increasing input power causes the output mode occupancy to redistribute from the fundamental mode to the highest-order modes, with the highest-order LP21 modes carrying about 90% of the power at the largest powers. The authors interpret this redistribution thermodynamically: they retain the Rayleigh-Jeans distribution, add a Kerr-type potential energy U_NL = γP²/(2M), modify the equation of state to U_L + U_NL − μP = MT, and infer from this model that the optical temperature rises with power, diverges at an equipartition point, and then approaches zero from below, defining a negative-absolute-temperature attractor. Supporting experiments include cut-back tests for thermalization, power-conservation checks, offset-injection robustness measurements, and observations at several wavelengths.","tokens_in":10466,"tokens_out":5568,"duration_ms":68629,"significance":"If the thermodynamic interpretation is correct, the paper would be a notable experimental demonstration of a negative optical temperature in a multimode waveguide and of a robust, all-optically controlled modal inversion, with potential interest for beam shaping and high-power fiber lasers. The experimental core is valuable: the modal decomposition is state-of-the-art, the cut-back and power-conservation checks support the claim of a thermalized, lossless propagation, and the offset-injection data provide a convincing demonstration of the robustness of the high-order-mode dominance. However, the central quantitative claim—that the system reaches a negative temperature approaching zero from below—is not directly measured but is inferred from a model that contains an unproved equation of state and a questionable calibration procedure. Those issues are load-bearing and require substantial revision before the negative-absolute-temperature attractor claim can be accepted.","major_comments":[{"comment":"The procedure of computing U_L once from the lowest-power mode occupancy and then 'keeping it as a constant for all power levels' when solving Eq. (8) is inconsistent with the definition U_L = −Σβ_i|c_i|² in Eq. (2). The measured modal redistribution from LP01 to LP21 directly implies that U_L/P changes substantially with power; for an input profile that is simply rescaled in power, U_L/P should remain near −β_1, whereas fixing U_L drives U_L/P toward zero as P grows. This forced crossing of the critical value −Σβ_i/M in Eq. (5) is what produces the temperature divergence and the subsequent T→0⁻ branch in the authors' Fig. 2d and Fig. 3b. Since the central claim of a negative-absolute-temperature attractor depends on this branch, the authors must either recompute T using a power-dependent U_L(P) obtained from the measured occupancies or provide a physical justification for why U_L, as used in Eq. (8), should be the low-power linear energy rather than the instantaneous kinetic energy.","section":"Methods, 'Determination of thermodynamic parameters'"},{"comment":"Equation (7), U_L + U_NL − μP = MT, is acknowledged in the manuscript to 'cannot be directly derived' from the Rayleigh-Jeans law and the conservation laws of P and U. This relation is not a minor technicality: it is the equation used to convert the measured mode occupancies into the optical temperature values plotted in Fig. 3b and to locate the T→0⁻ attractor. Because the equation of state is not derived, the quantitative temperature axis, the critical power at which T diverges, and the strength of the attractor are all model-dependent. The authors should either derive Eq. (7) from a controlled approximation of the nonlinear Hamiltonian or treat the model as a purely phenomenological fit and clearly separate the directly measurable quantities (mode occupancies, U_L) from the inferred thermodynamic parameters.","section":"Main text, 'Real photon gas'"},{"comment":"The nonlinear coefficient γ is obtained by fitting Eq. (8) to the highest-power data point, and the same point is the one that most strongly determines the negative-temperature branch. Consequently, the agreement between the theoretical curve and the experimental occupancies at the highest power in Fig. 3a is not an independent validation of the model. The authors should provide a sensitivity analysis showing how the inferred temperatures and the predicted critical powers vary with γ and with the choice of the reference power used to fix U_L, including confidence intervals on γ from the fit.","section":"Fig. 3 and Methods, 'Determination of thermodynamic parameters'"},{"comment":"The statement that the negative-temperature attractor occurs 'no matter the sign of the temperature at low power' is supported experimentally by the offset-injection measurements in Fig. 4, but those measurements are shown at only two power levels and are compared with the model only qualitatively. The theoretical statement of a universal attractor to T→0⁻ in Fig. 2c,d is computed from the unproved equation of state (7) and the fixed-U_L procedure. To make the attractor claim quantitative, the authors should report for the offset data the model-inferred temperatures as a function of power and verify that they indeed converge to the same negative near-zero value regardless of the low-power sign of T.","section":"Main text, 'Robustness of spatial beam freezing'"}],"minor_comments":[{"comment":"The paper uses 'demonstrate' for the negative absolute temperature regime, but T is always inferred from the model rather than directly measured. The wording should be softened, e.g., 'consistent with' or 'inferred from the thermodynamic model,' and the distinction between measured mode occupancies and inferred thermodynamic parameters should be made explicit in the abstract.","section":"Abstract and Main text"},{"comment":"The normalized nonlinear coefficient γ takes the value 0.05 in Fig. 2 but 2.04 in Fig. 3, with no statement of units or normalization. Please specify how γ is normalized, and report the fitted value with a confidence interval.","section":"Fig. 2 caption / Fig. 3"},{"comment":"The manuscript reports average input powers (e.g., 0.1 mW, 3.3 mW, 4 mW) for femtosecond pulses at 100 kHz repetition rate, but the thermodynamic description and the nonlinear potential energy depend on the instantaneous power or pulse energy. The authors should report peak powers or pulse energies and justify using average power in the thermodynamic fits.","section":"Methods, 'Experiments'"},{"comment":"The sign conventions in the denominator of Eq. (8) are easy to misread; a short derivation in the Methods section showing how Eq. (8) follows from Eqs. (3) and (7) would improve clarity.","section":"Eq. (8)"},{"comment":"The text states that U_NL remains about two orders of magnitude lower than U_L while also describing the gas as 'dense.' This apparent tension should be clarified: the non-ideal nature of the gas is invoked through the equation of state even when the potential energy is small compared with the linear energy.","section":"Main text, 'Real photon gas'"}],"recommendation":"major_revision","confidential_remarks":"The experimental work appears carefully executed, and the modal inversion itself is well supported. The main risk is that the negative-absolute-temperature attractor is an artifact of the fixed-U_L calibration combined with the unproved equation of state (7). I would urge the handling editor to require the authors to either recompute the temperatures directly from the measured occupancies via the Rayleigh-Jeans law alone or provide a derivation of Eq. (7) and a consistent power-dependent U_L. If those revisions are not feasible, the claims should be downgraded to a model-dependent interpretation rather than a demonstrated thermodynamic property."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"Quick take: the paper has a credible experimental discovery—power-induced modal inversion in a few-mode fiber, robust to offset injection—but the headline claim of a negative absolute temperature attractor is not supported. The thermodynamic analysis has a load-bearing calibration step that forces the sign flip.\n\nWhat's new: the observation that high power drives the output modal distribution from the fundamental to the highest-order modes, with the highest-order pair holding about 90% of power, is new and well documented. The cut-back, power conservation, and offset experiments are the right checks. The polychromatic freezing across spectral slices is a nice bonus. The authors are also honest: they explicitly state that Eq. (7), the modified equation of state, cannot be derived from the conservation laws.\n\nWhere it falls down: the determination of T in Methods. They compute U_L once from the lowest-power mode occupancy and keep it constant for all powers. But the experiment's own observable—modal redistribution—means U_L changes substantially. For a fixed input profile rescaled in power, U_L should scale with P, not stay fixed. Holding U_L constant forces U_L/P to drift toward zero as P grows, which pushes the system across the critical value in Eq. (5) and produces the divergence of T and the eventual T→0− branch. In other words, the sign flip is baked into the calibration. The equation of state is also a guess, not derived, and γ is fitted to the highest-power point, which is exactly where the attractor is claimed. Removing the fixed-U_L assumption would give different temperatures, and it is entirely possible the sign flip disappears or at least shifts. Since T is not measured, only inferred from this model, the negative-temperature claim is not established. The modal inversion could be real physics explained by something else (e.g., nonlinear mode coupling), but the paper does not make that case.\n\nBottom line: the experiments are worth engaging with, and a serious referee should see this. The paper needs major revision: either derive the equation of state, measure T directly, or at minimum recompute temperatures with a self-consistent U_L(P) from the measured mode occupancies. As it stands, I would not cite the negative-temperature claim, but I would cite the empirical freezing effect.\n\nRecommendation: send to peer review, with a note to the authors that the fixed-U_L step is the central issue.","headline":"Solid experiments on modal inversion, but the negative-temperature attractor is not supported: the sign flip is forced by a fixed-U_L calibration and an unproved equation of state.","tokens_in":11058,"tokens_out":3277,"would_cite":true,"duration_ms":40910,"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":"This paper claims that increasing input power in a few-mode step-index fiber inverts the mode population from the fundamental to the highest-order mode, which it interprets as the optical temperature passing through infinity to a negative…","keywords":["negative absolute temperature","photon gas","multimode fiber","few-mode fiber","mode population inversion","spatial beam freezing","Rayleigh-Jeans distribution","Kerr nonlinearity"],"falsifier":"Recompute the linear internal energy from the experimentally measured mode occupancies at each power level instead of carrying the lowest-power value through all powers, then solve the paper's Eq. (8) with that updated value; if the temperature no longer flips sign or the $T = 0^-$ branch disappears, the negative-temperature attractor would be an artifact of the fixed-energy assumption.","tokens_in":9964,"feed_emoji":"🌡️","tokens_out":7325,"duration_ms":73675,"temperature":0.7,"pith_summary":"The paper tries to show that a dense gas of photons inside a few-mode optical fiber behaves like a non-ideal, Van der Waals gas, and that this has a dramatic consequence: as the input laser power grows, the equilibrium distribution of light among the fiber modes inverts, with the highest-order mode taking over from the fundamental mode. The authors interpret this as the optical temperature of the photon gas increasing until it crosses infinity and flips sign, landing in a negative-absolute-temperature regime that approaches zero from below. They call this final state \"spatial beam freezing\" and argue that it is a stable attractor, so at high power the output beam shape is essentially independent of injection conditions and external perturbations. The significance is practical as well as fundamental, because it offers an all-optical way to control laser beam shape for high-power sources, nonlinear imaging, and optical communications.","feed_headline":"More power flips photon gas to negative absolute temperature","feed_subtitle":"In a few-mode fiber, high power makes the highest-order mode dominate, locking the output beam regardless of launch conditions.","key_machinery":"The carrying object is the modified equation of state $U_L + U_{NL} - \\mu P = M T$, paired with the Rayleigh-Jeans mode distribution $|c_i|^2 = -T/(\\mu + \\beta_i)$, where $\\beta_i$ are the propagation constants of the guided modes. The new term is the nonlinear potential energy $U_{NL} = \\frac{\\gamma}{2}\\frac{P^2}{M}$, which is quadratic in total power and inversely proportional to the number of modes, so it acts like a density-dependent interaction energy of a Van der Waals gas. Holding $U_L$ fixed at its low-power value while solving the resulting equation for $T/P$ makes the temperature a nonlinear function of power: it rises, diverges at the equipartition condition, and then approaches $0^-$ as the highest-order mode becomes the only populated state. The sign flip at $T = \\pm\\infty$ is not a singularity of the model, because the Lagrange multiplier conjugate to momentum conservation is $1/T$.","core_discovery":"The paper's central claim is that adding a Kerr-type potential energy to the equation of state of a photon gas in a few-mode fiber changes the power dependence of the optical temperature so drastically that the mode population spontaneously inverts. In the weakly nonlinear ideal-gas regime, the Rayleigh-Jeans distribution with conserved power and momentum gives an optical temperature whose sign is fixed by the linear internal energy; the paper's modified equation of state, $U_L + U_{NL} - \\mu P = M T$, with $U_{NL} = \\frac{\\gamma}{2}\\frac{P^2}{M}$, makes $T$ a nonlinear function of power. As power rises, $T$ increases, diverges at the equipartition point $U_L = U_{L,c}$, then returns from $T = 0^-$, which corresponds to the highest-order mode being the only populated state. In the experiment, the measured mode occupancy follows the Rayleigh-Jeans form at every power, with the highest-order guided modes carrying about 90% of the power at the highest input powers, and the same freezing is observed at other wavelengths where the number of guided modes differs.","pith_inferences":["The model's key quantitative step is to keep the linear internal energy $U_L$ fixed at its lowest-power value while solving for the temperature at every power; because the measured mode populations redistribute strongly with power, recomputing $U_L$ from the measured distribution at each power could shift or remove the sign flip.","A testable extension is to check whether the critical power for the inversion scales as predicted with the mode count $M$, since the added potential energy is proportional to $P^2/M$; varying the fiber core or wavelength to change $M$ would directly probe the density mechanism.","The attractor behavior may be a more general phenomenon than this specific fiber model: any nonlinearity that adds a positive, power-dependent term to the internal energy could push a Rayleigh-Jeans gas toward $T = 0^-$, so similar inversion should appear in other multimode nonlinear systems."],"forward_implications":["At low power the output is dominated by the fundamental mode; above a threshold the same fiber outputs a beam dominated by the highest-order guided modes, so beam quality degradation appears as a predictable thermodynamic transition rather than an uncontrolled defect.","The $T = 0^-$ state is stable: changing the launch offset, bending the fiber, or squeezing it does not change the final beam, which is what makes the effect usable as a form of spatial mode locking.","Because freezing occurs at each wavelength slice of a spectrally broadened pulse, the mechanism extends to polychromatic beams, with the frozen mode being the highest-order mode guided at that wavelength.","The inversion happens at power levels well below the self-focusing damage threshold, so it is accessible in ordinary step-index optical fibers.","Regardless of whether the low-power gas starts with positive or negative temperature, raising the density always drives it toward $T = 0^-$ rather than away from it."],"supporting_citations":[{"why":"Supplies the thermodynamic framework of a multimode waveguide as a box for a photon gas, with optical temperature and chemical potential defined from power and momentum conservation.","marker":"[1]"},{"why":"Provides the Rayleigh-Jeans law as the classical counterpart of the Bose-Einstein distribution used for the mode occupancies.","marker":"[2]"},{"why":"Establishes that optical temperatures can be positive or negative in multimode systems and defines $T = 0$ as the singular regime where the highest-order mode dominates.","marker":"[8]"},{"why":"Describes the beam self-cleaning effect in highly multimode fibers, the positive-temperature phenomenon that this paper's mode inversion extends and contrasts with.","marker":"[11]"},{"why":"Supports the claim that mode power redistribution in fibers results from thermalization toward a Rayleigh-Jeans equilibrium.","marker":"[12]"},{"why":"Introduces the Kerr-potential-energy correction in a waveguide-array photon gas via a Joule-Thomson analogue, the approach this paper adapts into its equation of state.","marker":"[18]"},{"why":"Proves uniqueness of the physical $(T,\\mu)$ solution for the ideal photon gas, which underlies the ideal-gas equation of state that the paper modifies.","marker":"[19]"},{"why":"Describes the holographic mode decomposition method used to measure the experimental mode occupancies.","marker":"[24]"}],"fun_headline_variants":["Photon gas flips to negative absolute temperature","Power-induced negative temperature in photon gas","Negative temperature photon gas attractor","Dense photon gas enters negative temperature regime","Photon gas's power-driven negative temperature"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The whole temperature-versus-power curve is computed while holding the linear part of the internal energy fixed at its lowest-power value, even though the measured mode populations redistribute strongly with power; if that energy value actually changes, the inferred temperatures, the sign flip, and the negative-temperature branch would not follow.","fun_headline_variants_meta":{"raw":{"variants":["Photon gas flips to negative absolute temperature","Power-induced negative temperature in photon gas","Negative temperature photon gas attractor","Dense photon gas enters negative temperature regime","Photon gas's power-driven negative temperature"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.00093,"raw_usage":{"total_tokens":4009,"prompt_tokens":998,"completion_tokens":3011,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":614,"completion_tokens_details":{"reasoning_tokens":2947}},"tokens_in":614,"tokens_out":3011,"duration_ms":24527,"temperature":1.0,"reasoning_tokens":2947,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-07T13:34:01.913509+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Recompute the linear internal energy from the experimentally measured mode occupancies at each power level instead of carrying the lowest-power value through all powers, then solve the paper's Eq. (8) with that updated value; if the temperature no longer flips sign or the $T = 0^-$ branch disappears, the negative-temperature attractor would be an artifact of the fixed-energy assumption.","supporting_citations":[{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Supplies the thermodynamic framework of a multimode waveguide as a box for a photon gas, with optical temperature and chemical potential defined from power and momentum conservation."},{"cited_title":"Zanaglia, J","cited_arxiv_id":null,"evidence_quote":"Provides the Rayleigh-Jeans law as the classical counterpart of the Bose-Einstein distribution used for the mode occupancies."},{"cited_title":"Baudin, J","cited_arxiv_id":null,"evidence_quote":"Establishes that optical temperatures can be positive or negative in multimode systems and defines $T = 0$ as the singular regime where the highest-order mode dominates."},{"cited_title":"Krupa, A","cited_arxiv_id":null,"evidence_quote":"Describes the beam self-cleaning effect in highly multimode fibers, the positive-temperature phenomenon that this paper's mode inversion extends and contrasts with."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Introduces the Kerr-potential-energy correction in a waveguide-array photon gas via a Joule-Thomson analogue, the approach this paper adapts into its equation of state."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Proves uniqueness of the physical $(T,\\mu)$ solution for the ideal photon gas, which underlies the ideal-gas equation of state that the paper modifies."},{"cited_title":"Gervaziev, I","cited_arxiv_id":null,"evidence_quote":"Describes the holographic mode decomposition method used to measure the experimental mode occupancies."}],"review_version":1}