Soliton Energy Ratio is a normalized measure quantifying the energy carried by soliton modes relative to the total energy, with definitions varying by governing equations.
It spans diverse contexts—from spectral fractions in ocean waves to peak density ratios in soliton collisions—displaying bounded, linear, or super-unit scaling based on the setup.
This metric aids in diagnosing soliton-dominated regimes and evaluating energy transfer and amplification in both conservative and dissipative systems.
In the cited literature, the expression soliton energy ratio denotes a family of normalized energy measures attached to nonlinear coherent structures. The numerator may be the energy carried by soliton modes in a spectral decomposition, the energy contained in a soliton-dominated frequency band, the peak local energy density generated in a multi-soliton collision, or the energy of a moving, dissipative, or quantum-corrected soliton state relative to a reference configuration. Accordingly, the quantity may be bounded in [0,1], may scale linearly with a control parameter, or may exceed unity; its precise meaning is determined by the governing equation and by the decomposition used to identify the soliton sector (Costa et al., 2014, Lee et al., 6 Oct 2025, Saadatmand et al., 2015, Kalashnikov et al., 2023).
1. Principal definitions and scope
Across the cited works, the same label is used for several distinct ratios. In ocean-wave applications, it usually means the fraction of total wave energy carried by soliton modes. In integrable field theory, it may instead compare the maximum local energy density attained during a collision to the total energy of the incoming solitons. In optical-fibre, dissipative, and relativistic settings, it frequently denotes a transport or parameter-continuation quotient rather than a spectral fraction.
Setting
Definition
Characteristic behavior
Shallow-water KdV waves
Rsoliton=Esoliton/Etotal with Esoliton=∫0ωcSsoliton(ω)dω (Costa et al., 2014)
Controlled by linear-in-Rsoliton=Esoliton/Etotal9 slopes and offsets
A common source of ambiguity is that only some of these ratios are literal “fractions of total energy.” In the NLS-NFT sea-state formulation the ratio is explicitly constrained to Esoliton=∫0ωcSsoliton(ω)dω0, whereas in the sine-Gordon collision problem the numerator is a peak density and the ratio grows as Esoliton=∫0ωcSsoliton(ω)dω1, so values above unity are expected rather than anomalous.
2. Shallow-water KdV formulation and finite-gap extraction
In the shallow-water analysis of Currituck Sound, the total mean wave energy per unit span is taken, in linear theory, to be proportional to the integral of the measured power spectrum,
Esoliton=∫0ωcSsoliton(ω)dω2
The soliton contribution is defined as the low-frequency component attributed to the KdV soliton band,
Esoliton=∫0ωcSsoliton(ω)dω3
The measured spectrum exhibits Esoliton=∫0ωcSsoliton(ω)dω4 for Esoliton=∫0ωcSsoliton(ω)dω5 and a higher-frequency wind-wave cascadeEsoliton=∫0ωcSsoliton(ω)dω6; in practice Esoliton=∫0ωcSsoliton(ω)dω7 is chosen at the observed spectral minimum separating the KdV band from the higher-frequency wind-wave band (Costa et al., 2014).
The dynamical basis is the Korteweg-deVries equation
Esoliton=∫0ωcSsoliton(ω)dω8
A single-soliton solution in a moving frame is written as
Esoliton=∫0ωcSsoliton(ω)dω9
and its energy is defined by the Rsoliton≃0.6±0.10 norm
Rsoliton≃0.6±0.11
Under periodic or quasiperiodic boundary conditions, the exact KdV solution is characterized by Rsoliton≃0.6±0.12 real branch points Rsoliton≃0.6±0.13. In the soliton limit each spectral band-gap collapses to repeated eigenvalues, yielding Rsoliton≃0.6±0.14 discrete real eigenvalues Rsoliton≃0.6±0.15, each mapped one-to-one onto a soliton amplitude Rsoliton≃0.6±0.16, speed Rsoliton≃0.6±0.17, and energy Rsoliton≃0.6±0.18. Numerically, finite-gap theory is applied to the low-frequency portion of the time series to recover the set Rsoliton≃0.6±0.19.
The empirical significance of the ratio in this setting is twofold. First, finite-gap analysis shows that the low-frequency Rsoliton=Esol/Etotal=(Etotal−Erad)/Etotal∈[0,1]0 region is soliton dominated. Second, the solitons have random FGT phases, supporting a soliton random phase approximation. The probability density of the solitons demonstrates that they are dense in time and highly non-Gaussian. Although the study did not tabulate Rsoliton=Esol/Etotal=(Etotal−Erad)/Etotal∈[0,1]1 explicitly, the supplied reconstruction states that all of the low-frequency power is carried by approximately Rsoliton=Esol/Etotal=(Etotal−Erad)/Etotal∈[0,1]2 solitons in each Rsoliton=Esol/Etotal=(Etotal−Erad)/Etotal∈[0,1]3 min record near storm peak, and that integrating the observed Rsoliton=Esol/Etotal=(Etotal−Erad)/Etotal∈[0,1]4 law up to Rsoliton=Esol/Etotal=(Etotal−Erad)/Etotal∈[0,1]5 against the full measured spectrum gives typically
so that over half of the surface-wave energy is contained in a dense gas of KdV solitons at peak storm conditions in shallow water with Rsoliton=Esol/Etotal=(Etotal−Erad)/Etotal∈[0,1]7 (Costa et al., 2014).
3. Deep-ocean soliton gases and the nonlinear Fourier transform
In deep water, the same phrase is given a more explicitly spectral meaning through the nonlinear Schrödinger nonlinear Fourier transform. Let Rsoliton=Esol/Etotal=(Etotal−Erad)/Etotal∈[0,1]8 be the complex envelope of a normalized surface-elevation record, extended by zero outside the measurement interval. The total energy is the standard Rsoliton=Esol/Etotal=(Etotal−Erad)/Etotal∈[0,1]9 norm,
Rsoliton≥0.900
and the nonlinear Parseval relation decomposes it into discrete and continuous parts,
Rsoliton≥0.901
Here Rsoliton≥0.902 is the energy carried by discrete eigenvalues, each corresponding to an NLS soliton, and Rsoliton≥0.903 is the energy carried by the continuous spectrum. The soliton energy ratio is then
Rsoliton≥0.904
For the normalization used in the study,
Rsoliton≥0.905
where Rsoliton≥0.906 are discrete Zakharov-Shabat eigenvalues and Rsoliton≥0.907 is the reflection coefficient (Lee et al., 6 Oct 2025).
The extraction pipeline is applied to each Rsoliton≥0.908-minute record sampled at Rsoliton≥0.909. The raw elevation is detrended and windowed, the carrier frequency Rsoliton>0.500 is estimated and removed by multiplication with Rsoliton>0.501, and the envelope is normalized to unit-coefficient focusing NLS scaling. The forward scattering problem
Rsoliton>0.502
is then solved to locate all discrete eigenvalues with Rsoliton>0.503 and to compute the continuous reflection coefficient.
In this formulation, the soliton energy ratio is used as a criterion for the closeness of a sea state to a soliton gas. A record is flagged as “very high soliton energy” when Rsoliton>0.504. These cases are associated with wave steepness Rsoliton>0.505 exceeding Rsoliton>0.506 and Benjamin-Feir index Rsoliton>0.507 exceeding Rsoliton>0.508. Because directional interference can artificially increase the ratio, the study applies a probabilistic directional-filtering method: a retention angle Rsoliton>0.509 about the mean direction is selected, off-axis bands are removed from R(N)=emax(N)/Etot(N)0, random uniform phases are assigned, R(N)=emax(N)/Etot(N)1 synthetic realizations are reconstructed, and a distribution of R(N)=emax(N)/Etot(N)2 is formed. After this correction, R(N)=emax(N)/Etot(N)3 is used to declare a soliton gas in the principal propagation direction.
The field results are unusually specific. The dataset comprises R(N)=emax(N)/Etot(N)4 deep-water, unimodal-spectrum R(N)=emax(N)/Etot(N)5-min records from three buoys in Taiwan waters, with R(N)=emax(N)/Etot(N)6–R(N)=emax(N)/Etot(N)7, R(N)=emax(N)/Etot(N)8, R(N)=emax(N)/Etot(N)9, R(N)=N/40, and directional spread R(N)=N/41. Eleven events had raw R(N)=N/42, corresponding to R(N)=N/43 of the dataset. These records had typically R(N)=N/44, R(N)=N/45–R(N)=N/46, R(N)=N/47, R(N)=N/48, nearly symmetric statistics with skewness approximately zero and kurtosis R(N)=N/49, and no large rogue-wave indices. One Eluanbi example at 2019-05-16 14:00 UTC had N0, N1, N2, N3, N4, and N5 discrete eigenvalues reaching N6 in amplitude. After directional filtering, three Eluanbi cases retained mean N7 for both N8 and N9, and are therefore reported as confirmed deep-ocean soliton gases (Lee et al., 6 Oct 2025).
4. Collision, amplification, and transport ratios in optical systems
In dispersion-mapped optical fibres, the relevant quantity is an energy-transfer ratio between colliding solitons. Two well-separated fundamental solitons are launched into a fibre governed by a generalized nonlinear Schrödinger equation with piecewise-constant group-velocity dispersion, anomalous outside a short segment of length Rsoliton=Esoliton/Etotal00 and normal within it. Each input pulse has energy
Rsoliton=Esoliton/Etotal01
with Rsoliton=Esoliton/Etotal02. After a collision in the normal-dispersion section, the output energies are Rsoliton=Esoliton/Etotal03 and Rsoliton=Esoliton/Etotal04. The transfer into soliton 1 is quantified by
Rsoliton=Esoliton/Etotal05
and the corresponding gain ratio is
Rsoliton=Esoliton/Etotal06
For fixed physical parameters, the transfer is fitted semi-analytically by
Rsoliton=Esoliton/Etotal07
with phase offset Rsoliton=Esoliton/Etotal08, and amplitudes scaling as power laws in Rsoliton=Esoliton/Etotal09, Rsoliton=Esoliton/Etotal10, Rsoliton=Esoliton/Etotal11, Rsoliton=Esoliton/Etotal12, Rsoliton=Esoliton/Etotal13, and Rsoliton=Esoliton/Etotal14. The quoted best-fit parameters are
Rsoliton=Esoliton/Etotal15
Rsoliton=Esoliton/Etotal16
At Rsoliton=Esoliton/Etotal17, Rsoliton=Esoliton/Etotal18, Rsoliton=Esoliton/Etotal19, Rsoliton=Esoliton/Etotal20, Rsoliton=Esoliton/Etotal21, and Rsoliton=Esoliton/Etotal22, the maximum transfer is approximately Rsoliton=Esoliton/Etotal23 at Rsoliton=Esoliton/Etotal24, while the minimum is approximately Rsoliton=Esoliton/Etotal25 at Rsoliton=Esoliton/Etotal26. A representative single collision yields Rsoliton=Esoliton/Etotal27, and a second engineered collision raises this to Rsoliton=Esoliton/Etotal28, corresponding to a cumulative gain of approximately Rsoliton=Esoliton/Etotal29. The mechanism relies on the dispersion-sign change and does not require third-order dispersion or Raman terms (Savojardo et al., 2017).
A different optical meaning appears in lossy fibres, where the ratio tracks propagation-induced depletion or compensation. The fundamental-mode energy is
Rsoliton=Esoliton/Etotal30
From the perturbed CPDE model,
Rsoliton=Esoliton/Etotal31
Assuming a Rsoliton=Esoliton/Etotal32-shaped fundamental soliton,
Rsoliton=Esoliton/Etotal33
one obtains Rsoliton=Esoliton/Etotal34 and therefore
Rsoliton=Esoliton/Etotal35
with Rsoliton=Esoliton/Etotal36 and Rsoliton=Esoliton/Etotal37. The corresponding ratio satisfies
Rsoliton=Esoliton/Etotal38
with logistic-form solution
Rsoliton=Esoliton/Etotal39
where Rsoliton=Esoliton/Etotal40 and Rsoliton=Esoliton/Etotal41. The subcases are explicit: Rsoliton=Esoliton/Etotal42 for pure linear loss and Rsoliton=Esoliton/Etotal43 for pure cubic gain or loss. For simulations with Rsoliton=Esoliton/Etotal44, the first Rsoliton=Esoliton/Etotal45 where Rsoliton=Esoliton/Etotal46 occurs at approximately Rsoliton=Esoliton/Etotal47, Rsoliton=Esoliton/Etotal48, Rsoliton=Esoliton/Etotal49, Rsoliton=Esoliton/Etotal50, and Rsoliton=Esoliton/Etotal51 for Rsoliton=Esoliton/Etotal52, Rsoliton=Esoliton/Etotal53, Rsoliton=Esoliton/Etotal54, Rsoliton=Esoliton/Etotal55, and Rsoliton=Esoliton/Etotal56, respectively. In this sense, a gain of Rsoliton=Esoliton/Etotal57 triples the propagation distance when the dissipation rate is Rsoliton=Esoliton/Etotal58 (DalľAgnol et al., 2021).
5. Relativistic and dissipative formulations
In the reduced-QED soliton-like construction, the energy ratio compares a moving state to its rest state rather than separating solitonic and radiative sectors. After neglecting transverse photons and adopting a self-consistent mean-field ansatz in the Coulomb gauge, the rest-state energy is obtained from coupled nonlinear Dirac equations and a Poisson-type self-consistent potential. The lowest-energy solution is reported as
Rsoliton=Esoliton/Etotal59
A canonical separation of total momentum then yields the relativistic dispersion relation
Rsoliton=Esoliton/Etotal60
The corresponding ratio is
Rsoliton=Esoliton/Etotal61
with the physical positive branch written as Rsoliton=Esoliton/Etotal62. The small-momentum expansion is
Rsoliton=Esoliton/Etotal63
while in the ultra-relativistic limit Rsoliton=Esoliton/Etotal64 one has Rsoliton=Esoliton/Etotal65 (Skoromnik et al., 2016).
In the cubic-quintic Ginzburg-Landau equation, the energy ratio is instead a parameter-continuation quotient on a master diagram. The dissipative-soliton energy is
Rsoliton=Esoliton/Etotal66
and in the adiabatic strongly chirped limit with Rsoliton=Esoliton/Etotal67 the spectrum has Rayleigh-Jeans form
Rsoliton=Esoliton/Etotal68
which integrates to
Rsoliton=Esoliton/Etotal69
The family is controlled by the universal parameter
Rsoliton=Esoliton/Etotal70
At the vacuum-stability threshold Rsoliton=Esoliton/Etotal71,
Rsoliton=Esoliton/Etotal72
The soliton-energy ratio between two parameter sets is then
Rsoliton=Esoliton/Etotal73
Dissipative soliton resonance is identified by the condition Rsoliton=Esoliton/Etotal74 while the peak power remains finite at Rsoliton=Esoliton/Etotal75, which occurs at
Rsoliton=Esoliton/Etotal76
At this threshold, Rsoliton=Esoliton/Etotal77, Rsoliton=Esoliton/Etotal78, and Rsoliton=Esoliton/Etotal79. The same formulation admits a thermodynamic interpretation with inverse spectral temperature Rsoliton=Esoliton/Etotal80, chemical potential Rsoliton=Esoliton/Etotal81, entropy Rsoliton=Esoliton/Etotal82, and free energy Rsoliton=Esoliton/Etotal83 (Kalashnikov et al., 2023).
6. Topological charge, quantum corrections, and comparative interpretation
For two-component solitons in one spatial dimension, the relevant ratios are organized by the topological charge Rsoliton=Esoliton/Etotal84. Both the classical energy and the vacuum polarization energy are fitted to very high numerical accuracy by straight lines,
Rsoliton=Esoliton/Etotal85
For the parameter choice Rsoliton=Esoliton/Etotal86, Rsoliton=Esoliton/Etotal87, the reported values are
Rsoliton=Esoliton/Etotal88
With loop-counting parameter Rsoliton=Esoliton/Etotal89, the total energy becomes
Rsoliton=Esoliton/Etotal90
Two ratios are then natural:
Rsoliton=Esoliton/Etotal91
and
Rsoliton=Esoliton/Etotal92
Because both energies are linear in Rsoliton=Esoliton/Etotal93, the binding properties are controlled by the offsets rather than curvature. The critical condition is set by solving
Rsoliton=Esoliton/Etotal94
For small Rsoliton=Esoliton/Etotal95, all higher-charge solitons are bound; once Rsoliton=Esoliton/Etotal96 exceeds the critical value, all higher-charge states become unbound. The reported numerical behavior is that Rsoliton=Esoliton/Etotal97 is typically of order Rsoliton=Esoliton/Etotal98 for the parameter sets studied (Graham et al., 6 Mar 2025).
A common misconception is to treat soliton energy ratio as a single invariant fraction of soliton content. The cited literature does not do so. It uses the same term for spectral soliton fractions in KdV and NLS sea states, for localization efficiency in sine-Gordon collisions, for propagation and amplification quotients in optical fibres, for relativistic dispersion in a reduced-QED soliton-like state, and for parameter- or charge-dependent quotients in dissipative and quantum-corrected soliton theories. This suggests that the term identifies a problem-dependent diagnostic rather than a universal observable.