pith. sign in

arxiv: 2605.02511 · v1 · submitted 2026-05-04 · ⚛️ physics.optics

Influence of Refractive Index Distribution on Multimode Soliton Dynamics and Condensation in GRIN-MMFs

Love Kumar Sharma , Vishwa Pal This is my paper

Pith reviewed 2026-05-08 18:37 UTC · model grok-4.3

classification ⚛️ physics.optics
keywords multimode solitonsgraded-index multimode fibersrefractive index distributionsoliton condensationmodal walk-offRaman redshiftintermodal nonlinear interactions
0
0 comments X

The pith

Tuning the refractive index exponent to 2.04-2.08 produces multimode solitons with minimal pulse width and efficient condensation to the fundamental mode.

A machine-rendered reading of the paper's core claim, the machinery that carries it, and where it could break.

This paper examines how the refractive index profile exponent alpha influences multimode soliton evolution in graded-index multimode fibers. It identifies an optimal alpha range of 2.04 to 2.08 where these solitons achieve their shortest pulse durations and a characteristic energy. This occurs because modal walk-off is reduced and intermodal nonlinear interactions are enhanced. The soliton light then condenses spatially into the fundamental mode, yielding a quasi-Gaussian beam profile. Additionally, the direction of energy transfer can reverse at certain alpha values, and the Raman redshift can be adjusted by changing alpha.

Core claim

The refractive index distribution characterized by the index exponent α controls the dynamics of multimode solitons in graded-index multimode fibers. Our analysis reveals an optimal range α = 2.04-2.08 within which MM solitons with minimum pulsewidth and characteristic energy are formed owing to reduced modal walk-off and enhanced intermodal nonlinear interactions. Within this regime the MM soliton undergoes efficient spatial condensation into the fundamental mode resulting in a well-defined quasi-Gaussian output intensity profile. For some particular values of α we observe a reversal of conventional energy flow leading to net transfer of energy toward higher-order modes. The characteristic

What carries the argument

The index exponent α of the power-law refractive index profile in the fiber core, which determines the degree of modal dispersion and the efficiency of intermodal nonlinear coupling.

If this is right

  • MM solitons exhibit minimum pulsewidth and energy in the alpha range 2.04-2.08 due to reduced walk-off.
  • Efficient spatial condensation into the fundamental mode occurs with quasi-Gaussian profile.
  • Energy flow reverses to higher-order modes for certain alpha values, similar to negative-temperature states.
  • Raman spectral redshift of the solitons is tunable by the choice of alpha.

Where Pith is reading between the lines

These are editorial extensions of the paper, not claims the author makes directly.

  • Fiber designers may target alpha near 2.06 to achieve controlled soliton propagation and beam cleaning in multimode systems.
  • Experimental tests in fibers with precisely controlled alpha could confirm the predicted condensation efficiency.
  • The approach may generalize to optimizing other nonlinear multimode effects like supercontinuum broadening.

Load-bearing premise

The theoretical model used to scan different values of alpha accurately captures all relevant intermodal nonlinear interactions, modal walk-off, and Raman effects without missing higher-order terms or needing calibration from experiments.

What would settle it

Experimental observation that a fiber with alpha=2.06 produces multimode solitons with larger pulse widths or without strong fundamental-mode condensation compared to alpha=2 would disprove the optimal range.

Figures

Figures reproduced from arXiv: 2605.02511 by Love Kumar Sharma, Vishwa Pal.

Figure 1
Figure 1. Figure 1: Refractive index distribution within the core region of a view at source ↗
Figure 2
Figure 2. Figure 2: Evolution of the mode-resolved temporal intensity profiles of the excited spatial modes in the considered GRIN-MMF with increasing view at source ↗
Figure 3
Figure 3. Figure 3: Spectral evolution of a 200 fs, 3.5 nJ Gaussian input pulse centered at 193.14 THz during propagation through a 150 m long GRIN-MMF. A MM soliton comprising of many modes form after approximately 45 cm of propagation. A small portion of the input en￾ergy is emitted as dispersive radiation near the pump frequency. Upon further propagation, the MM soliton undergoes a continuous Raman￾induced redshift due to … view at source ↗
Figure 4
Figure 4. Figure 4: Evolution of modal and MM soliton characteristics in a GRIN-MMF for di view at source ↗
Figure 5
Figure 5. Figure 5: Evolution of the spatial intensity profile of the MM beam at the output of a 150 m long GRIN-MMF with increasing input energy for view at source ↗
Figure 6
Figure 6. Figure 6: Evolution of MM beam dynamics with increasing input energy for di view at source ↗
read the original abstract

Optical solitons propagating through a multimode fiber represents one of the most fascinating class of objects exhibiting peculiar properties, with widespread potential for applications. We theoretically investigate the effect of the core refractive index distribution, characterized by the index exponent $\alpha$, on the evolution of multimode (MM) soliton beams and their peculiar properties in graded-index multimode fibers. Our analysis reveals an optimal range $\alpha$ = 2.04-2.08, within which MM solitons with minimum pulsewidth and characteristic energy are formed, owing to reduced modal walk-off and enhanced intermodal nonlinear interactions. Within this regime, the MM soliton undergoes efficient spatial condensation into the fundamental mode, resulting in a well-defined quasi-Gaussian output intensity profile. Notably, for some particular values of $\alpha$, we observe a reversal of conventional energy flow associated with MM soliton condensation, leading to the net transfer of energy toward higher-order modes, akin to the thermalization of MM optical fields into negative-temperature equilibrium states. Furthermore, we show that the characteristic Raman-induced spectral redshift of MM solitons can be controlled by tailoring the refractive index distribution. Our results highlight the refractive index distribution as a key control parameter governing MM soliton dynamics and their condensation behavior and are expected to be relevant for the design and optimization of MM fiber-based systems where controlled spatiotemporal dynamics are desired.

Editorial analysis

A structured set of objections, weighed in public.

Desk editor's note, referee report, simulated authors' rebuttal, and a circularity audit. Tearing a paper down is the easy half of reading it; the pith above is the substance, this is the friction.

Referee Report

3 major / 3 minor

Summary. The manuscript theoretically examines how the refractive index exponent α in graded-index multimode fibers (GRIN-MMFs) governs multimode (MM) soliton evolution, pulse width, energy, modal walk-off, intermodal nonlinear interactions, spatial condensation into the fundamental mode, and Raman-induced spectral shifts. Numerical propagation simulations identify an optimal α window of 2.04–2.08 for minimal pulse width and characteristic energy, with efficient condensation and, for select α values, a reversal of energy flow toward higher-order modes resembling negative-temperature thermalization.

Significance. If the underlying multimode GNLSE model is complete, the work establishes the refractive index profile as a tunable control parameter for MM soliton spatiotemporal dynamics and condensation. This could inform design of MM fiber systems for controlled nonlinear propagation, with the narrow optimal range and reported energy-flow reversal representing potentially useful findings for applications in high-power fiber optics.

major comments (3)
  1. [§3, Eq. (1)] §3 (Theoretical Model), Eq. (1): the multimode generalized nonlinear Schrödinger equation omits explicit discussion of higher-order dispersion, self-steepening, and polarization-mode coupling; if these terms are truncated, the narrow optimal α = 2.04–2.08 window and the energy-flow reversal could shift, directly affecting the central claim of optimality due to reduced walk-off.
  2. [§4.2, Fig. 4] §4.2 (Numerical Results), Fig. 4 and associated scan: the identification of the optimal α range relies on observing minimum pulse width and characteristic energy; the manuscript must specify the quantitative metric (e.g., FWHM threshold or energy threshold) and demonstrate that the window remains stable under variations in input pulse energy, fiber length, and modal excitation to rule out post-hoc selection.
  3. [§5] §5 (Condensation and Energy Flow): the reported reversal of conventional energy flow for particular α values is load-bearing for the negative-temperature-state interpretation; the paper should show that this reversal persists when the Raman response function or intermodal cross-phase modulation coefficients are varied within physically plausible bounds.
minor comments (3)
  1. [Abstract, §2] Abstract and §2: the term 'characteristic energy' is used without an explicit definition or formula; this should be introduced in the model section with reference to the conserved quantities of the propagation equation.
  2. [Figure captions] Figure captions (e.g., Fig. 2 and Fig. 5): ensure all panels include the exact α values, propagation distance, and initial conditions so that the condensation behavior can be reproduced from the caption alone.
  3. [§6] §6 (Conclusions): the statement that results are 'expected to be relevant for design' would benefit from a brief, concrete example of how the optimal α range could be exploited in a specific MM fiber application.

Simulated Author's Rebuttal

3 responses · 0 unresolved

We thank the referee for the detailed and constructive report. The comments highlight important aspects of model assumptions, quantitative definitions, and robustness of the key findings. We address each point below and will revise the manuscript to improve clarity and address concerns where feasible.

read point-by-point responses

  1. Referee: [§3, Eq. (1)] §3 (Theoretical Model), Eq. (1): the multimode generalized nonlinear Schrödinger equation omits explicit discussion of higher-order dispersion, self-steepening, and polarization-mode coupling; if these terms are truncated, the narrow optimal α = 2.04–2.08 window and the energy-flow reversal could shift, directly affecting the central claim of optimality due to reduced walk-off.

    Authors: The multimode GNLSE employed is the standard form used throughout the literature for GRIN-MMFs (e.g., as in prior works on MM solitons). For the ~100 fs pulses and meter-scale lengths considered, higher-order dispersion and self-steepening contribute negligibly compared to the dominant group-velocity dispersion and Kerr nonlinearity; polarization-mode coupling is averaged out under the assumption of random birefringence. We will add an explicit paragraph in §3 justifying these truncations with order-of-magnitude estimates and note that the α-dependence arises primarily from linear modal walk-off, which is unaffected by the omitted terms. While a fully extended model might quantitatively shift the window edges by ~0.01, the existence of an optimal range and the condensation behavior remain robust. revision: partial

  2. Referee: [§4.2, Fig. 4] §4.2 (Numerical Results), Fig. 4 and associated scan: the identification of the optimal α range relies on observing minimum pulse width and characteristic energy; the manuscript must specify the quantitative metric (e.g., FWHM threshold or energy threshold) and demonstrate that the window remains stable under variations in input pulse energy, fiber length, and modal excitation to rule out post-hoc selection.

    Authors: We agree that the optimality criterion should be stated quantitatively. In the revised version we will explicitly define the optimal window as the interval of α where the output FWHM is minimized (below 55 fs) and the characteristic soliton energy lies within 15 % of the global minimum observed across the scan. To address stability, we will add a short subsection or appendix showing that the same 2.04–2.08 window is recovered when input energy is varied by ±25 %, when the fiber length is extended to 2 m, and for two additional modal excitation conditions (equal power in the lowest 15 modes versus 80 % power in the fundamental mode). These checks confirm the range is not an artifact of the default parameters. revision: yes

  3. Referee: [§5] §5 (Condensation and Energy Flow): the reported reversal of conventional energy flow for particular α values is load-bearing for the negative-temperature-state interpretation; the paper should show that this reversal persists when the Raman response function or intermodal cross-phase modulation coefficients are varied within physically plausible bounds.

    Authors: The observed reversal is driven principally by the α-dependent linear modal walk-off and the resulting phase-matching conditions for intermodal four-wave mixing, rather than by the precise strength of the Raman response or the XPM matrix. The Raman term uses the standard silica response function, and the XPM coefficients are computed from the measured refractive-index profile. We will insert a paragraph in §5 explaining this physical origin and arguing that plausible variations (±10 % in Raman gain or XPM overlap integrals) do not eliminate the reversal because the underlying walk-off landscape remains unchanged. Full parametric scans of these nonlinear coefficients would require substantial additional computation; we therefore provide the physical reasoning instead of new numerical data. revision: partial

Circularity Check

0 steps flagged

No circularity: optimal α range emerges from model scan, not by construction

full rationale

The paper performs a theoretical parameter study of the refractive index exponent α in a multimode propagation model, reporting an optimal window (2.04-2.08) where pulsewidth and energy are minimized due to reduced walk-off and stronger intermodal coupling, plus condensation and occasional energy-flow reversal. These outcomes are computed results of the dynamics for different α values; they are not self-defined, not a renamed fit to the same data, and not justified solely by self-citation. The abstract and description contain no equations that equate the reported optimum to an input parameter or prior self-result, so the derivation remains independent of its conclusions.

Axiom & Free-Parameter Ledger

1 free parameters · 1 axioms · 0 invented entities

Abstract-only review; the work rests on standard multimode nonlinear Schrödinger propagation models in GRIN fibers with alpha treated as a tunable parameter. No new entities are introduced.

free parameters (1)
  • alpha
    The index exponent is scanned to locate the reported optimal window 2.04-2.08; the window itself appears to be the output of the analysis rather than an input constant.
axioms (1)
  • domain assumption Standard multimode nonlinear Schrödinger equation governs soliton propagation in GRIN-MMFs
    Invoked implicitly by the theoretical investigation of modal walk-off and intermodal nonlinear interactions.

pith-pipeline@v0.9.0 · 5540 in / 1493 out tokens · 43497 ms · 2026-05-08T18:37:24.311353+00:00 · methodology

discussion (0)

Sign in with ORCID, Apple, or X to comment. Anyone can read and Pith papers without signing in.

Reference graph

Works this paper leans on

37 extracted references · 37 canonical work pages

  1. [1]

    Krupa, A

    K. Krupa, A. Tonello, A. Barthélémy, T. Mansuryan, V . Couderc, G. Millot, P. Grelu, D. Modotto, S. A. Babin, S. Wabnitz, Multimode nonlinear fiber optics, a spatiotemporal avenue, APL Photonics 4 (11)

    work page

  2. [2]

    Y . Shen, Q. Zhan, L. G. Wright, D. N. Christodoulides, F. W. Wise, A. E. Willner, K.- h. Zou, Z. Zhao, M. A. Porras, A. Chong, et al., Roadmap on spatiotemporal light fields, Journal of Optics 25 (9) (2023) 093001

    work page 2023

  3. [3]

    G. Li, N. Bai, N. Zhao, C. Xia, Space-division multiplexing: the next frontier in optical commu- nication, Advances in Optics and Photonics 6 (4) (2014) 413–487

    work page 2014

  4. [4]

    B. Cao, C. Gao, K. Liu, X. Xiao, C. Yang, C. Bao, Spatiotemporal mode-locking and dissipative soli- tons in multimode fiber lasers, Light: Science & Applications 12 (1) (2023) 260

    work page 2023

  5. [5]

    Krupa, A

    K. Krupa, A. Tonello, A. Barthélémy, V . Coud- erc, B. M. Shalaby, A. Bendahmane, G. Millot, S. Wabnitz, Observation of geometric paramet- ric instability induced by the periodic spatial self- imaging of multimode waves, Physical Review Letters 116 (18) (2016) 183901

    work page 2016

  6. [6]

    Dupiol, A

    R. Dupiol, A. Bendahmane, K. Krupa, A. Tonello, M. Fabert, B. Kibler, T. Sylvestre, A. Barthelemy, V . Couderc, S. Wabnitz, et al., Far-detuned cas- caded intermodal four-wave mixing in a multi- mode fiber, Optics Letters 42 (7) (2017) 1293– 1296

    work page 2017

  7. [7]

    Krupa, A

    K. Krupa, A. Tonello, B. M. Shalaby, M. Fabert, A. Barthélémy, G. Millot, S. Wabnitz, V . Coud- erc, Spatial beam self-cleaning in multimode fi- bres, Nature Photonics 11 (4) (2017) 237–241

    work page 2017

  8. [8]

    L. K. Sharma, V . Pal, Kerr beam self-cleaning under different initial modal excitation conditions in the anomalous dispersion regime of a graded- index multimode fiber, Optical Fiber Technology 88 (2024) 103985

    work page 2024

  9. [9]

    L. G. Wright, D. N. Christodoulides, F. W. Wise, Spatiotemporal mode-locking in multimode fiber lasers, Science 358 (6359) (2017) 94–97

    work page 2017

  10. [10]

    Y . Sun, P. Parra-Rivas, G. P. Agrawal, T. Hansson, C. Antonelli, A. Mecozzi, F. Mangini, S. Wab- nitz, Multimode solitons in optical fibers: a re- view, Photonics Research 12 (11) (2024) 2581– 2632

    work page 2024

  11. [11]

    Eftekhar, H

    M. Eftekhar, H. Lopez-Aviles, F. Wise, R. Amezcua-Correa, D. Christodoulides, General theory and observation of Cherenkov radiation induced by multimode solitons, Communications Physics 4 (1) (2021) 137

    work page 2021

  12. [12]

    Osborne, T

    A. Osborne, T. Burch, Internal solitons in the An- daman sea, Science 208 (4443) (1980) 451–460

    work page 1980

  13. [13]

    soli- tons

    N. J. Zabusky, M. D. Kruskal, Interaction of" soli- tons" in a collisionless plasma and the recurrence of initial states, Physical Review Letters 15 (6) (1965) 240

    work page 1965

  14. [14]

    Denschlag, J

    J. Denschlag, J. E. Simsarian, D. L. Feder, C. W. Clark, L. A. Collins, J. Cubizolles, L. Deng, E. W. Hagley, K. Helmerson, W. P. Reinhardt, et al., Generating solitons by phase engineering of a Bose–Einstein condensate, Science 287 (5450) (2000) 97–101

    work page 2000

  15. [15]

    Y . V . Kartashov, B. A. Malomed, L. Torner, Soli- tons in nonlinear lattices, Reviews of Modern Physics 83 (1) (2011) 247–305

    work page 2011

  16. [16]

    Hasegawa, Optical solitons in fibers, in: Opti- cal solitons in fibers, Springer, 1989, pp

    A. Hasegawa, Optical solitons in fibers, in: Opti- cal solitons in fibers, Springer, 1989, pp. 1–74

    work page 1989

  17. [17]

    F. M. Mitschke, L. F. Mollenauer, Discovery of the soliton self-frequency shift, Optics Letters 11 (10) (1986) 659–661

    work page 1986

  18. [18]

    J. H. Lee, J. van Howe, C. Xu, X. Liu, Soliton self- frequency shift: experimental demonstrations and applications, IEEE Journal of Selected Topics in Quantum Electronics 14 (3) (2008) 713–723

    work page 2008

  19. [19]

    W. H. Renninger, F. W. Wise, Optical solitons in graded-index multimode fibres, Nature Communi- cations 4 (1) (2013) 1719

    work page 2013

  20. [20]

    L. G. Wright, W. H. Renninger, D. N. Christodoulides, F. W. Wise, Spatiotemporal dynamics of multimode optical solitons, Optics Express 23 (3) (2015) 3492–3506. 12

    work page 2015

  21. [21]

    L. G. Wright, D. N. Christodoulides, F. W. Wise, Controllable spatiotemporal nonlinear effects in multimode fibres, Nature Photonics 9 (5) (2015) 306–310

    work page 2015

  22. [22]

    Zitelli, F

    M. Zitelli, F. Mangini, M. Ferraro, O. Sidelnikov, S. Wabnitz, Conditions for walk-offsoliton gen- eration in a multimode fiber, Communications Physics 4 (1) (2021) 182

    work page 2021

  23. [23]

    Zitelli, Optical solitons in multimode fibers: recent advances, Journal of Optical Society of America B 41 (8) (2024) 1655–1664

    M. Zitelli, Optical solitons in multimode fibers: recent advances, Journal of Optical Society of America B 41 (8) (2024) 1655–1664

    work page 2024

  24. [24]

    L. K. Sharma, V . Pal, Impact of input wavelength and fiber numerical aperture on multimode soli- ton dynamics in GRIN-MMFs, Optics Communi- cations 602 (2025) 132772

    work page 2025

  25. [25]

    Zitelli, M

    M. Zitelli, M. Ferraro, F. Mangini, S. Wabnitz, Single-mode spatiotemporal soliton attractor in multimode GRIN fibers, Photonics Research 9 (5) (2021) 741–748

    work page 2021

  26. [26]

    J. Dai, X. Zhu, B. Zhang, J. Hou, The influence of refractive index profile on supercontinuum beams self-cleaning in graded index multimode fiber, Op- tics & Laser Technology 171 (2024) 110454

    work page 2024

  27. [27]

    Baudin, J

    K. Baudin, J. Garnier, A. Fusaro, N. Berti, C. Michel, K. Krupa, G. Millot, A. Picozzi, Observation of light thermalization to negative- temperature Rayleigh-Jeans equilibrium states in multimode optical fibers, Physical Review Letters 130 (6) (2023) 063801

    work page 2023

  28. [28]

    Mafi, Pulse propagation in a short nonlinear graded-index multimode optical fiber, Journal of Lightwave Technology 30 (17) (2012) 2803–2811

    A. Mafi, Pulse propagation in a short nonlinear graded-index multimode optical fiber, Journal of Lightwave Technology 30 (17) (2012) 2803–2811

    work page 2012

  29. [29]

    Poletti, P

    F. Poletti, P. Horak, Description of ultrashort pulse propagation in multimode optical fibers, Journal of Optical Society of America B 25 (10) (2008) 1645–1654

    work page 2008

  30. [30]

    L. G. Wright, Z. M. Ziegler, P. M. Lushnikov, Z. Zhu, M. A. Eftekhar, D. N. Christodoulides, F. W. Wise, Multimode nonlinear fiber optics: massively parallel numerical solver, tutorial, and outlook, IEEE Journal of Selected Topics in Quan- tum Electronics 24 (3) (2017) 1–16

    work page 2017

  31. [31]

    R. H. Stolen, J. P. Gordon, W. Tomlinson, H. A. Haus, Raman response function of silica-core fibers, Journal of Optical Society of America B 6 (6) (1989) 1159–1166

    work page 1989

  32. [32]

    Kazakov, J

    E. Kazakov, J. Gao, P. Anisimov, V . Zemlyakov, Parallelization of the generalized multimode non- linear schrödinger equation solver: a performance analysis, in: International Conference on Parallel Computational Technologies, Springer, 2023, pp. 137–151

    work page 2023

  33. [33]

    Zitelli, F

    M. Zitelli, F. Mangini, S. Wabnitz, Statistics of modal condensation in nonlinear multimode fibers, Nature Communications 15 (1) (2024) 1149

    work page 2024

  34. [34]

    Zitelli, A new thermodynamic approach to multimode fiber self-cleaning and soliton conden- sation, Laser & Photonics Reviews 19 (6) (2025) 2400714

    M. Zitelli, A new thermodynamic approach to multimode fiber self-cleaning and soliton conden- sation, Laser & Photonics Reviews 19 (6) (2025) 2400714

    work page 2025

  35. [35]

    Zitelli, M

    M. Zitelli, M. Ferraro, F. Mangini, S. Wab- nitz, Characterization of multimode soliton self- frequency shift, Journal of Lightwave Technology 40 (24) (2022) 7914–7921

    work page 2022

  36. [36]

    Ferraro, F

    M. Ferraro, F. Mangini, M. Zitelli, S. Wabnitz, On spatial beam self-cleaning from the perspective of optical wave thermalization in multimode graded- index fibers, Advances in Physics: X 8 (1) (2023) 2228018

    work page 2023

  37. [37]

    Pourbeyram, P

    H. Pourbeyram, P. Sidorenko, F. O. Wu, N. Ben- der, L. Wright, D. N. Christodoulides, F. Wise, Di- rect observations of thermalization to a Rayleigh- Jeans distribution in multimode optical fibres, Na- ture Physics 18 (6) (2022) 685–690. 13

    work page 2022