{"id":"b84899c3-9f4e-4801-b632-7d1e8e046927","arxiv_id":"1908.07745","paper_version":1,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":4.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":3,"one_line_summary":"In a coupled-mode simulation, random coupling among degenerate modes is shown to reproduce both the speckled output and the Kerr beam self-cleaning observed in graded-index multimode fiber experiments.","lead":"This paper tests several models of random mode coupling in a simulation of light traveling through a multimode fiber. The model that couples only modes of the same family reproduces the speckle and beam self-cleaning seen in experiments.","discovery_kind":"extension","skeptic_critique":{"model":"deepseek-v4-flash","headline":"Eq. (6) writes the random-coupling term as -(1/2)C instead of -(i/2)C, so the reported evolution is not unitary and does not follow from Eq. (4); the claimed agreement with experiment rests on an unverified or misprinted model equation.","rationale":"I agree with the reader that the random-coupling model is under-specified, particularly the undisclosed variance and correlation length of C, and that no quantitative comparison to experiment is shown. However, I see a more fundamental and concrete issue: Eq. (6) does not follow from Eq. (4) unless an i is inserted in the coupling term. This is not merely a missing parameter; it changes the unitarity of the linear propagator and could invalidate the reported numerical results. The paper provides no code, no data, and no experimental overlay, so the reader cannot determine whether the printed equation was actually solved or whether the missing i is a typographical artifact. Because the central claim is a claim of agreement with experiment, the absence of quantitative validation is already load-bearing; the Eq. (6) inconsistency makes the situation worse. The concern is potentially fixable: a corrected equation, a statement of the coupling parameters, and a direct comparison with Ref. [14] would allow the claim to be evaluated. Thus I keep the reader's CONDITIONAL verdict, but with the added condition that the propagation equation be corrected or justified and that the missing model parameters be reported. This is why I mark agreement as partial: the reader identified the under-specification, which is real, but not the more precise equation-level inconsistency that I find most decisive.","tokens_in":5550,"tokens_out":7040,"duration_ms":74513,"concrete_test":"Independently re-derive the matrix form: divide Eq. (4) by 2i and check whether the coupling term is -(i/2)C A. Then, with access to the code or a rerun, simulate the equal-n coupling model using both the corrected -iC/2 term and the printed -C/2 term, monitoring total power vs propagation distance. The concern is settled if the corrected simulation conserves total power and reproduces the linear speckle and nonlinear self-cleaning at the same parameter values, while the printed form does not. Additionally, request the variance and correlation length of C and a quantitative comparison (e.g., radial intensity profiles or mode-power spectra) with the experimental data in Ref. [14].","verdict_should_be":"UNCHANGED","load_bearing_attack":"The central claim is that the equal-mode-number random-coupling model reproduces the linear speckle and nonlinear self-cleaning data. The decisive weak point is the matrix equation actually used. Eq. (4), after division by 2i, gives dA/dζ = ... - (i/2) C A ..., but Eq. (6) states M = ... - (1/2) C A ... . Since C is declared Hermitian to preserve total power, the printed -C/2 term is Hermitian and makes the linear evolution non-unitary (exponential growth or decay of the total power), contradicting the stated power-preservation requirement. If the leapfrog scheme in Eq. (9) solved Eq. (6) literally, the speckle and self-cleaning could be artifacts of gain/loss rather than the proposed coupling physics; if it solved Eq. (4), then Eq. (6) is a typo and the manuscript still withholds the variance and correlation length of C for the equal-n case and shows no quantitative side-by-side comparison with the experimental data of Ref. [14]. In either branch, the 'complete agreement' claim is not currently verifiable from the text.","agreement_with_reader":"partial"},"referee_report":{"model":"deepseek-v4-flash","summary":"The manuscript reports numerical simulations of Kerr beam self-cleaning in a graded-index multimode fiber using a coupled-mode model. Random linear coupling between modes is added to the base model of Ref. [14], and three coupling scenarios are compared: coupling between all modes, coupling between neighboring modes satisfying n1=n±1 and m1=m±1, and coupling only among degenerate modes with equal mode number n. The authors find that only the degenerate-mode coupling reproduces a speckled output in the linear regime and a stable self-cleaned beam in the nonlinear regime, and they claim that the results are in complete agreement with their own experimental data from Ref. [14].","tokens_in":5825,"tokens_out":6454,"duration_ms":69851,"significance":"If the model were fully specified and the agreement with experiment quantified, this would be a useful contribution: it offers a computationally efficient coupled-mode description and identifies degenerate-mode random coupling as the minimal ingredient that preserves the self-cleaning mechanism while generating realistic speckle. The paper's strengths are the systematic comparison of three coupling models, the direct calculation of modal overlap integrals, and the clear separation of linear and nonlinear regimes. The significance is currently limited, however, because the central claim of complete agreement is not backed by a single quantitative comparison, the random-coupling parameters are not reported, and the printed evolution equation appears to contain an error that breaks power conservation.","major_comments":[{"comment":"The matrix form in Eq. (6) reads dA/dζ = M A with M = -iD L/2 - (1/2) C - i p NL/2, but substituting Eq. (4) gives a coupling term of -(i/2) C A. Because C is declared Hermitian in order to preserve total power, the printed -C/2 term is Hermitian with the wrong sign of the imaginary unit, making the linear evolution non-unitary and capable of introducing spurious exponential growth or decay. If the leapfrog scheme in Eq. (9) solves Eq. (6) literally, the reported speckle and self-cleaning could be numerical artifacts of this gain or loss; if it solves Eq. (4), then Eq. (6) is a typo that must be corrected. In either case the text needs a clear statement of which equation is actually integrated.","section":"Section 2, Eqs. (4) and (6)"},{"comment":"The parameters of the random coupling matrix C are never fully specified. The text states that each element is normally distributed with zero mean, but it does not give the standard deviation or variance of that distribution, and it gives a correlation length of 10 cm only for the n±1 model, not for the equal-mode-number model that is claimed to reproduce the experiments. The all-modes model is said to have coefficients that change at every integration step, which would make the correlation length equal to the step size rather than a physical value. Since the speckle contrast and the self-cleaning dynamics depend on the coupling strength, the model is under-specified and the reported agreement cannot be reproduced or checked.","section":"Section 3, random-coupling models"},{"comment":"The statement that the numerical results are in complete agreement with the experimental data of Ref. [14] is not supported by any quantitative comparison in the manuscript. There is no side-by-side plot of simulated and measured output intensity profiles, no metric such as speckle contrast or output beam radius versus power, and no error bars. The conclusion repeats the claim without evidence. Please add at least one quantitative comparison, for example a plot of fundamental-mode fraction or output beam radius versus input power for the same fiber parameters and launch conditions as in Ref. [14].","section":"Section 3 and Conclusion"}],"minor_comments":[{"comment":"The sentence 'the beam will not change during propagation in the linear regime' is not correct for Eq. (4) with C=0: the D(n+1)^2 term produces power-conserving phase evolution among modes and hence self-imaging of the beam. Please either account for this phase evolution in the discussion or state explicitly that D is neglected in that argument.","section":"Section 3, first paragraph of Numerical results"},{"comment":"The symbol p is used both for normalized power (p = P/Psf) and as a mode index in A_{p,m} and in the subscripts of the overlap integrals, which makes equations such as (4) and (8) unnecessarily hard to read. Please use separate notation, for example P0 for normalized power.","section":"Throughout"},{"comment":"There are several typographical and grammatical errors that should be corrected in a revision: 'Intoduction' in the section heading, 'profle' instead of 'profile', 'we start out investigation', and 'completely agreement' in the abstract-like phrasing.","section":"Throughout"},{"comment":"The numerical scheme is presented as an explicit leapfrog finite-difference method, but the text does not discuss the chosen integration step size, numerical stability, or convergence checks. A sentence reporting these details would be useful, especially because the random coefficients vary along the propagation distance.","section":"Section 3, Eq. (9)"}],"recommendation":"major_revision","confidential_remarks":"The central claim of the paper is, at present, a self-referential check: both the base coupled-mode model and the experimental benchmark come from Ref. [14] by the same group, and the random-coupling parameters are not disclosed. The apparent sign error in Eq. (6) is the most serious technical issue because it directly affects what the numerical code may actually be solving. I believe the work is within the scope of the journal and the main idea is plausible, but the manuscript needs a corrected equation, full reporting of the coupling parameters, and a genuine quantitative comparison with the experimental data before it can be accepted."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"The paper you asked about is a short numerical study of Kerr beam self-cleaning in graded-index multimode fiber. The thing worth knowing: the authors test three random-coupling models and find that only coupling among degenerate modes (equal mode number n) reproduces both the linear speckle and the nonlinear self-cleaning observed in their earlier experiment. That negative result is a useful contribution to model selection.\n\nWhat the paper does well is the systematic comparison of coupling rules. The all-modes coupling gives rapid energy oscillations and no stable speckle; the n±1 coupling distorts but doesn't speckle; the equal-n coupling, with a slightly offset input, gives a realistic speckle and then self-cleaning at high power. The narrative is clear and the logic is honest about why the centered input fails (symmetry selection).\n\nThe soft spots are real, and one is load-bearing. Equation (6), the matrix form actually solved, has -1/2 C where Eq. (4) gives -i/2 C. With C Hermitian, -1/2 C is Hermitian, so the linear evolution is not unitary: total power would grow or decay, contradicting the paper's own requirement that C be Hermitian to preserve power. If the code solved Eq. (6) literally, the speckle could be a gain/loss artifact. If it solved Eq. (4), then Eq. (6) is a typo. Either way the manuscript can't be verified as written. This is not a nitpick; it's the central model.\n\nThe other issues are more standard. The \"complete agreement\" with experiment is supported by no quantitative comparison: no error bars, no side-by-side figure, no metric. The variance of the random coupling is never reported for the equal-n case, only a correlation length (10 cm) that appears in the nearest-neighbor model. No code or data is released. So the reader cannot reproduce the figures or judge whether the coupling strength was tuned.\n\nThe citation pattern is fine; the benchmark is the authors' own prior experiment, which is legitimate but means independent confirmation is still open.\n\nThis paper is for people actively modeling multimode fiber nonlinearities. It doesn't open a new direction, but it gives a concrete, computationally cheap model that explains a known effect. I'd send it to peer review, but a serious referee should demand a corrected Eq. (6), full parameter reporting, and a quantitative comparison before it can be considered verified.","headline":"A numerically useful model showing that degenerate-mode random coupling explains Kerr self-cleaning, but a sign error in the central equation and missing quantitative support make the current claims unverifiable.","tokens_in":6323,"tokens_out":5541,"would_cite":false,"duration_ms":48399,"reading_group":"maybe","serious_thinker":"unclear","would_accept_peer_review":true},"rs_alignment":null,"lean_confirmation":null,"pith_extraction":{"msc":[],"pacs":[],"model":"deepseek-v4-flash","headline":"Random coupling among modes with equal mode number, added to a coupled-mode model, reproduces both the linear speckle and the nonlinear beam self-cleaning observed in a graded-index multimode fiber.","keywords":["beam self-cleaning","multimode optical fiber","graded-index fiber","Kerr effect","random mode coupling","coupled-mode model","degenerate modes","speckle pattern"],"falsifier":"Launch a single higher-order mode into a real 62.5-micrometer graded-index fiber and measure how its power leaks into modes with different mode numbers over sub-meter scales; if significant non-degenerate coupling appears at the 10 cm scale, or if the coupling strength required to fit the speckle statistics is far outside measured fiber perturbations, the degenerate-mode coupling model would be ruled out.","tokens_in":5381,"feed_emoji":"🔆","tokens_out":5439,"duration_ms":91371,"temperature":0.7,"pith_summary":"The paper argues that Kerr beam self-cleaning in a graded-index multimode optical fiber can be understood within a coupled-mode model once random linear coupling between degenerate spatial modes is added. Without such coupling, the model produces no linear speckle and the fundamental mode power oscillates instead of stabilizing. With coupling restricted to modes sharing the same mode number n, and with a slightly off-center input beam, the model reproduces both the speckled output in the linear regime and the collapse to a near-fundamental-mode beam at high power. The authors report complete agreement with their experimental data, making degenerate-mode coupling the operative ingredient rather than coupling among all modes or among neighboring modes only.","feed_headline":"Equal-mode random coupling explains fiber beam self-cleaning","feed_subtitle":"A coupled-mode model with degenerate-mode coupling reproduces both the speckled linear output and the cleaned nonlinear beam seen in…","key_machinery":"The central object is a random Hermitian coupling matrix C inserted into the coupled-mode equations for a graded-index fiber. Each off-diagonal element is normally distributed with zero mean and changes randomly along propagation with a 10 cm correlation length; the particular variant that works restricts coupling to modes with equal mode number n, i.e., degenerate modes. This matrix models fiber imperfections such as bending, tilting, and fabrication irregularities. It randomizes the phases and energy distribution among higher-order modes to create linear speckle, while leaving the fundamental mode energetically isolated, which the Kerr term then exploits to condense the beam into a stable bell-shaped output.","core_discovery":"The central claim is that random linear coupling between degenerate modes, meaning modes with equal mode number n = 2p + |m|, is the key mechanism that lets Kerr nonlinearity clean a noisy multimode beam. In the linear regime this coupling scrambles energy among modes of the same n and, for an input beam slightly displaced from the fiber axis, produces a realistic speckle pattern; in the nonlinear regime the fundamental mode, being the only n = 0 mode, cannot lose energy through this coupling, so its power stabilizes and dominates, giving an output beam close to the fundamental mode. Other coupling prescriptions, such as coupling all modes or only neighboring modes, fail to reproduce the observed speckle and self-cleaning. The paper concludes that the degenerate-mode coupling model is in complete agreement with the available experimental data.","pith_inferences":["If degenerate-mode coupling is indeed the operative perturbation, fiber designs that minimize non-degenerate perturbations, such as precisely controlled micro-bending, should show a sharper self-cleaning threshold, which could be tested experimentally.","The paper does not report the variance or full statistics of the coupling matrix; a natural next step is to extract these statistics from single-mode-launch experiments and test whether the model's quantitative predictions survive without tuning.","The equal-n coupling rule suggests a connection to mode condensation: the Kerr nonlinearity acts on a beam whose higher-order modes have already been randomized by disorder, which may link this work to hydrodynamic descriptions of multimode fibers."],"forward_implications":["The coupled-mode model with degenerate-mode random coupling reproduces the linear speckled output observed in experiments, provided the input beam is slightly off-center.","In the nonlinear regime, the same model yields robust beam self-cleaning to a near-fundamental-mode output at 10 kW input power.","Coupling all modes, or only neighboring modes with n1 = n ± 1 and m1 = m ± 1, does not reproduce the experimental speckle or self-cleaning; only equal-n coupling works.","Because the fundamental mode is the sole n = 0 mode, its power is insensitive to degenerate coupling and can stabilize under Kerr nonlinearity.","Using the coupled-mode model with a large integration step reduces computation time compared with full three-dimensional nonlinear Schrödinger propagation."],"supporting_citations":[{"why":"Supplies the experimental data against which the degenerate-mode coupling model is compared; the paper reports complete agreement with it.","marker":"[14]"},{"why":"Reported the spatial beam self-cleaning effect in multimode fibers that this model aims to reproduce.","marker":"[11]"},{"why":"Demonstrated Kerr self-cleaning of femtosecond-pulsed beams in graded-index multimode fiber, providing another experimental baseline.","marker":"[12]"},{"why":"Observed self-organized instability in graded-index multimode fibers, relevant to the nonlinear regime modeled here.","marker":"[13]"}],"fun_headline_variants":["Degenerate-mode coupling drives fiber beam self-cleaning","Random coupling among equal modes cleans fiber beams","Key to Kerr beam cleaning: degenerate-mode coupling","Equal-mode coupling scrambles and cleans fiber beams"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The model's agreement depends on real fiber imperfections behaving like random coupling that mixes only modes with the same mode number, reshuffled along a 10 cm length, with a strength that the paper does not report.","fun_headline_variants_meta":{"raw":{"variants":["Degenerate-mode coupling drives fiber beam self-cleaning","Random coupling among equal modes cleans fiber beams","Key to Kerr beam cleaning: degenerate-mode coupling","Equal-mode coupling scrambles and cleans fiber beams"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000407,"raw_usage":{"total_tokens":2029,"prompt_tokens":771,"completion_tokens":1258,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":387,"completion_tokens_details":{"reasoning_tokens":1199}},"tokens_in":387,"tokens_out":1258,"duration_ms":114657,"temperature":1.0,"reasoning_tokens":1199,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-14T11:56:39.978043+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Launch a single higher-order mode into a real 62.5-micrometer graded-index fiber and measure how its power leaks into modes with different mode numbers over sub-meter scales; if significant non-degenerate coupling appears at the 10 cm scale, or if the coupling strength required to fit the speckle statistics is far outside measured fiber perturbations, the degenerate-mode coupling model would be ruled out.","supporting_citations":[{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Supplies the experimental data against which the degenerate-mode coupling model is compared; the paper reports complete agreement with it."},{"cited_title":"Krupa, A","cited_arxiv_id":null,"evidence_quote":"Reported the spatial beam self-cleaning effect in multimode fibers that this model aims to reproduce."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Demonstrated Kerr self-cleaning of femtosecond-pulsed beams in graded-index multimode fiber, providing another experimental baseline."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Observed self-organized instability in graded-index multimode fibers, relevant to the nonlinear regime modeled here."}],"review_version":1}