Two-color solitons in Kerr third harmonic generation model

Alejandro B Aceves; Alexey Sukhinin; Curtis Menyuk; Jean-Claude Diels; Natalia Litchinitser

REVIEW 4 major objections 3 minor 12 references

Two-color solitons in Kerr third harmonic generation model

T0 review · 4 major / 3 minor · reviewed 2026-08-01 · deepseek-v4-flash

Challenge this review Re-run · record.json Download PDF Read on arXiv ↗ Review JSON

Pith's one-line read This paper shows that in a resonantly coupled two-color Kerr medium, the fate of a beam—collapse or diffraction—is decided by the distribution of power between the harmonics, not by a universal critical power.

desk verdict Interesting numerical claim about power-ratio thresholds in resonant THG, but the evidence for the separatrix is still only a few runs; worth a referee with a request for convergence and basin-boundary checks. read the letter →

arxiv 2607.16435 v1 pith:HCGQ3DDM submitted 2026-07-17 physics.optics math-phmath.MPnlin.PS

Alexey Sukhinin , Curtis Menyuk , Natalia Litchinitser , Jean-Claude Diels , Alejandro B Aceves This is my paper · ORCID

classification physics.opticsmath-phmath.MPnlin.PS PACS 42.65.Jx42.65.Sf42.65.Tg52.38.Hb52.35.Sb

keywords self-focusingKerrnonlinearitysolitonscross-phasemodulationthirdharmonicgenerationcriticalpowercollapsetwo-colorbeams

verification ladder T0 review T1 audit T2 compute T3 formal

The pith

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

The reading

This paper studies two-color, two-dimensional localized light states in a Kerr medium where the second harmonic has exactly three times the frequency of the fundamental. It finds that, unlike single-component self-focusing or non-resonant two-color models, there is no global critical power that decides collapse versus diffraction. Instead, the outcome depends on which harmonic carries the larger share of the total power. Fundamental-dominated stationary states act as dynamical separatrices: a small power increase makes both beams collapse together, a small decrease makes both diffract. Third-harmonic-dominated states do not show such a sharp separatrix, and collapse is preceded by strong oscillations of the third harmonic, a mechanism not present in standard Kerr self-focusing.

What carries the argument

The central object is the coupled (2+1)D nonlinear Schrödinger system describing Kerr third-harmonic generation, with cubic self- and cross-phase modulation plus resonant coupling terms. The argument is carried by stationary two-color solutions obtained numerically via a Newton-conjugate-gradient method, and by direct pseudo-spectral integration of perturbed stationary states. The simulations map out the collapse–diffraction boundary and reveal the oscillatory third-harmonic precursor that characterizes resonant collapse.

What would settle it

Rerun the ±10% perturbation experiments of the fundamental-dominated state (Fig. 1) at higher resolution (e.g., 1024×1024 or adaptive mesh) and larger domain, and vary the perturbation shape (e.g., Gaussian vs. Townes-like). If the +10% case diffracts or the −10% case collapses, or if the oscillations disappear, the central claim fails.

Watch

Extended reading notes

Core claim

The central claim is that in the resonant third-harmonic generation model, the dynamics of two-color localized beams are governed by the power distribution between the harmonics rather than by a scalar critical power. Numerically, the authors identify two families of unstable stationary solutions of the coupled 2D nonlinear Schrödinger equations. For the fundamental-dominated family, the stationary state sits exactly on a boundary: a +10% amplitude perturbation of the fundamental triggers simultaneous self-focusing of both components, while a −10% perturbation causes joint diffraction. For the third-harmonic-dominated family, no such clean boundary exists—moderate perturbations on either sid

Load-bearing premise

The central claim rests on a small number of pseudo-spectral runs at 512×512 spatial resolution, with no reported timestep, domain size, or convergence checks; if those runs misclassify collapse versus diffraction, the fundamental-dominated state is not a true separatrix.

Editorial extensions

If this is right

Experimental efforts to create co-existing UV/IR filaments should select initial conditions based on the power ratio between harmonics, not just total power.
The Marburger formula for collapse distance, developed for single-component beams, cannot be directly extended to resonantly coupled multicolor beams.
Strong third-harmonic oscillations can serve as an observable early-warning signal of imminent resonant collapse.
Any generalized collapse criterion for resonantly coupled Kerr media must be a function of the power distribution, not a single scalar threshold.

Reading between the lines

Editorial extensions of the paper, not claims the author makes directly.

If the separatrix claim holds, a continuous family of power-ratio thresholds likely exists, forming a critical manifold that interpolates between the fundamental- and third-harmonic-dominated regimes; the paper only exhibits two isolated separatrix states.
The same power-ratio logic may apply to higher-order resonant processes (e.g., fourth or fifth harmonic generation), suggesting that resonances add internal degrees of freedom that break the universality of collapse thresholds.
The observed third-harmonic oscillations may be a signature of periodic energy exchange between the modes that, when nonlinearity wins, destabilizes into collapse; a reduced two-mode oscillator model could be derived to predict the oscillation frequency and collapse onset.
A natural experimental test is to launch two-color beams with controlled power ratios and monitor whether the collapse threshold shifts along the predicted separatrix in the (P1, P2) plane.

Share X Bluesky LinkedIn Reddit HN

Editorial analysis

A structured set of objections, weighed in public.

Desk editor's note, referee report, and a circularity audit.

Desk Editor's Note

read the letter

Colleague,

Short version: this paper finds something genuinely new in a clean model — in resonant two-color Kerr collapse, the old single critical power appears replaced by dependence on how power is split between the harmonics, and fundamental-dominated stationary states look like separatrices. The model derivation is solid, the stationary states are computed with a proper Newton-CG method, and the taxonomy of fundamental- vs third-harmonic-dominated families is a useful organizing device. Credit where due: this is not a rehash of [9,10].

The soft underbelly is the dynamical evidence. The separatrix claim for Fig. 1 rests on exactly two amplitude perturbations: +10% collapses, -10% diffracts. For the third-harmonic-dominated family, three runs. No timestep, domain-size, or resolution convergence is reported; the statements refer to "up to 512x512" with no numerical error bars. The stress-test note is on point: these states are linearly unstable saddles, and a saddle automatically has trajectories on both sides of its stable manifold. Two adjacent trajectories do not prove the stationary state itself lies on the basin boundary. That requires either a bisection along a perturbation family, or a projection onto the stable manifold, or at minimum more perturbation shapes and longer runs.

I also think the abstract overstates the body. The body says "we do not claim a rigorous proof of singularity formation" and allows that the critical power notion "may not apply," but the abstract says unconditionally there is no universal critical power and the fundamental-dominated states act as dynamical separatrices. That gap should be closed in revision, either by softening the abstract or strengthening the evidence.

There is genuine value here for nonlinear optics and filamentation. The model is standard, the question is well-posed, and the potential implication — that UV/IR filament formation depends on the initial power ratio, not just total power — is concrete and testable. The missing pieces are all numerical hygiene: convergence tests, domain checks, a small parameter sweep around the stationary states, and ideally a basin-boundary estimate. I don't see a fatal flaw in the model or the classification.

Who is this for? Anyone working on multi-frequency self-focusing, THG, or collapse. It deserves a serious referee — send it to review, but with the expectation of a major revision report demanding more dynamics. I would cite it with a caveat as numerical evidence, not as settled law, and it's a good reading-group discussion of how thin the line is between suggestive numerics and a separatrix claim.

Best

Referee Report

4 major / 3 minor

Summary. The paper studies two-color localized states in a resonant Kerr third-harmonic generation model. Starting from the quasi-monochromatic approximation with ω2=3ω1, the authors derive a coupled (2+1)D NLS system (Eqs. (1)–(2)), rescale it, and solve for stationary two-color solitons using a Newton-conjugate-gradient method. They classify solutions into a 'fundamental-dominated' family (P1>P2, Figs. 1–4) and a 'third-harmonic-dominated' family (P2>P1, Figs. 5–6), and show the linear stability spectra, all of which contain eigenvalues with positive real parts. Direct simulations with perturbed stationary states as initial conditions are used to argue that fundamental-dominated states act as 'dynamical separatrices' between simultaneous collapse and joint diffraction, while third-harmonic-dominated states do not. The paper further identifies strong third-harmonic oscillations as a precursor to collapse in the fundamental-dominated family, and concludes that the concept of a universal critical power does not apply to this resonant system, with dynamics instead governed by the power distribution between harmonics.

Significance. If the central claims are correct, the paper identifies a genuinely new mechanism in multi-frequency Kerr self-focusing: a rational power-ratio threshold manifold rather than a scalar critical power, and an oscillatory resonant-collapse channel absent from the single-component NLS and from non-resonant two-color models. The model derivation is transparent, the stationary-state computations are standard and reproducible in principle, and the paper explicitly avoids fitted parameters, which is a strength. The results are potentially important for filamentation and for the theory of multi-frequency collapse. However, the dynamical claims currently rest on a very small number of numerical runs with no reported convergence study, and the evidence is not yet strong enough to support the qualitative conclusions as stated.

major comments (4)

[Collapse Events, Figs. 7–8] The central 'dynamical separatrix' claim for the fundamental-dominated family rests on exactly two perturbed runs: E1→1.1E1 collapses and E1→0.9E1 diffracts (Fig. 7). Since Fig. 1(d) shows the stationary state is linearly unstable, any saddle point has a codimension-one stable manifold, so one collapsing and one diffracting neighboring trajectory is generic and does not establish that the stationary state lies on the basin boundary. A proper separatrix identification requires a basin-boundary computation, e.g., bisection in perturbation amplitude/power for several perturbation shapes, including perturbations of E2 and of both components, plus confirmation that the threshold converges under resolution/domain/timestep refinement. This is load-bearing for the abstract's main claim.
[Collapse Events, Fig. 8] The conclusion that the third-harmonic-dominated family does not act as a separatrix is drawn from three propagations: +10% E1 collapses, -10% E1 collapses, and both -20% diffract. The absence of a detected threshold in three runs is not evidence for the absence of a threshold; a cliff between -10% and -20% in E1 amplitude is still a threshold. Moreover, no perturbations of E2 or mixed perturbations are explored. Systematic threshold and perturbation-shape scans are needed before the different dynamical role of this family can be accepted.
[Collapse Events, first paragraph] No convergence or numerical-validation data are reported. The text states 'spatial resolutions up to 512×512' but gives no domain size, timestep, temporal integrator, dealiasing/filtering procedure, or resolution study. In collapse simulations, peak amplitudes grow until limited by the grid, so without comparing e.g. 256², 512² and 1024² runs on a fixed larger domain, the distinction between genuine collapse-like dynamics and grid-scale artifacts cannot be evaluated. This numerical evidence underpins every qualitative claim about collapse and separatrix behavior.
[Abstract and Conclusions] The phrase 'no universal critical power' is stronger than the evidence presented. The perturbed runs in Figs. 7–8 compare states with different total powers and different profiles; the observation that a lower-power state (Fig. 8(a), total ≈ 11.5) collapses while a higher-power state (Fig. 7(c), total ≈ 15.3) diffracts is suggestive of a non-monotone total-power threshold, but only if those runs are reliable. The paper does not directly test power-ratio dependence at fixed total power or at fixed profile shape. A targeted numerical experiment isolating the power-distribution variable would substantiate (or refute) the title claim.

minor comments (3)

[Conclusions vs. Abstract] The Conclusions appropriately hedge ('we do not claim a rigorous proof of singularity formation'), but the Abstract and Introduction state the separatrix and no-universal-critical-power conclusions unconditionally. The wording should be aligned with the actual strength of the evidence.
[Figures 1–6] The figure captions do not give axis labels, spatial scales, or grid parameters. The radial profiles would be much more useful with explicit r-units and the amplitude normalization specified. This is needed for reproducibility.
[Search for Resonant Critical Power] This section presents the Townes-soliton background but does not connect it to the two-color problem. It could be shortened or explicitly linked to the two-component equations, e.g., by stating why the single-component critical power does not apply.

Circularity Check

0 steps flagged · score 1.0 of 10

No significant circularity: the numerical results are self-contained and the self-citations are not load-bearing.

full rationale

The paper's derivation chain is not circular. The coupled model (1)-(2) is obtained from the stated quasi-monochromatic ansatz and the cubic Kerr polarization under the resonance condition ω2 = 3ω1; no fitted parameter or target observable is inserted into the model. The stationary states are solved numerically from (3)-(4) using the Newton-conjugate-gradient method, and the claimed collapse/diffraction behaviors are forward simulations of (1)-(2) with ±10% perturbations of those states. These outcomes are not imposed by construction: the +/− perturbation runs are independent dynamical evolutions, not fits to a desired conclusion. The Townes constant α≈1.86225 is quoted as background from the single-component NLSE literature and is not used as an input to the two-color simulations. The self-citations [9,10] are used only to contrast the non-resonant case and do not carry the load of the paper's new claim. The manuscript does contain an explicit limitation—'we do not claim a rigorous proof of singularity formation'—and the reported evidence is sparse (a few pseudo-spectral runs with no convergence checks), but this is a correctness/evidence concern rather than a circularity. Under the stated rubric, no step reduces to its own input by definition or by fitted prediction.

Assumptions & free parameters 2 free parameters · 5 assumptions · 0 invented entities

No data-fitting parameters are used; the central numbers are outputs of simulations. The two families are separated by the hand-chosen criterion P1>P2 versus P2>P1, and the dynamics depend on hand-chosen perturbation sizes. No new physical entities are postulated; 'dynamical separatrix' is a classification label.

free parameters (2)

Stationary-state scale β (third-harmonic propagation constant) = not reported
Each computed localized state corresponds to some β via E_i=√βE_i(s), s=√βr, and the power values used for classification depend on it; β values are not listed, so the family branches cannot be directly reconstructed.
Perturbation amplitudes in collapse runs = 1.1, 0.9, 0.8 (multiplicative); also 1.1 only on E1
The separatrix classification is inferred from these hand-picked perturbations of one or two states in Figs. 7-8; no scan or justification is provided.

assumptions (5)

domain assumption Quasi-monochromatic slowly-varying-envelope approximation; only transverse diffraction included (paraxial (2+1)D model).
Used at the start to reduce Maxwell's equations to the coupled NLSEs (1)-(2); nonparaxial and dispersive corrections are neglected.
domain assumption Exact resonance/phase matching ω2=3ω1 and k2≈3k1; coherent THG terms E2(E1*)² and E1³ are retained.
Defines the resonant coupling that the whole paper contrasts with non-resonant XPM; enters through the nonlinear polarization expansion.
domain assumption Ideal instantaneous cubic Kerr nonlinearity with no saturation, dispersion, losses, or higher-order nonlinear terms.
Model equations (1)-(2) contain only cubic terms; if saturation/dispersion matters before collapse, the resonant collapse scenario may not occur.
domain assumption Newton-CG solutions of (3)-(4) are converged localized stationary states on the 512×512 grid.
No residuals, continuation parameters, or convergence criteria are reported; the two families are defined by these computed solutions.
domain assumption Pseudo-spectral Fourier propagation resolves collapse/diffraction accurately; no grid or boundary artifacts influence the outcomes.
No timestep, domain size, or resolution study is given, yet the separatrix conclusion depends on the collapse/diffraction classification of individual runs.

how reviews work

0 comments

Cite this review

Pith. "Pith review of Two-color solitons in Kerr third harmonic generation model." pith.science (2026). https://pith.science/paper/HCGQ3DDM

@misc{pith2026260716435,
  author       = {Pith},
  title        = {Pith review of: Two-color solitons in Kerr third harmonic generation model},
  year         = {2026},
  howpublished = {\url{https://pith.science/paper/HCGQ3DDM}},
  note         = {Machine review of arXiv:2607.16435}
}

read the original abstract

We investigate two-color, two dimensional spatially localized light modes in a resonant Kerr third-harmonic generation model. Using computational tools, we identify two distinct families of localized states. Unlike the single-component 2D nonlinear Schr\"odinger equation and previously studied non-resonant two-color systems, the dynamics are not dictated by a universal critical power, but depend on the distribution of power between the harmonics. The "fundamental-dominated" family acts as a "dynamical separatrix" between simultaneous collapse and joint diffraction, whereas the "third-harmonic-dominated" family does not. We further identify resonant collapse events accompanied by strong oscillations of the third harmonic, revealing a collapse mechanism absent from standard Kerr self-focusing.

Figures

Figures reproduced from arXiv: 2607.16435 by the authors.

**Figure 3.** FIG. 3. 3D profile of field [PITH_FULL_IMAGE:figures/full_fig_p003_3.png] view at source ↗

**Figure 2.** FIG. 2. 3D profile of field [PITH_FULL_IMAGE:figures/full_fig_p003_2.png] view at source ↗

Figures from the paper (2 more)

**Figure 7.** Figure 7: FIG. 7. (a) Propagation of peak amplitudes of perturbed [PITH_FULL_IMAGE:figures/full_fig_p004_7.png]

**Figure 6.** Figure 6: FIG. 6. 3D profile of field [PITH_FULL_IMAGE:figures/full_fig_p004_6.png]

Discussion (0). Continue with ORCID to comment.

Reference graph

Works this paper leans on

12 extracted references

[1]

R. Y. Chiao, E. Garmire, and C. H. Townes, Self-trapping of optical beams, Physical review letters13, 479 (1964)

1964
[2]

A. L. Gaeta, Catastrophic collapse of ultrashort pulses, Physical Review Letters84, 3582 (2000)

2000
[3]

K. Moll, A. L. Gaeta, and G. Fibich, Self-similar optical wave collapse: observation of the townes profile, Physical review letters90, 203902 (2003)

2003
[4]

Chen and C.-L

C.-A. Chen and C.-L. Hung, Observation of universal quench dynamics and townes soliton formation from modulational instability in two-dimensional bose gases, Physical Review Letters125, 250401 (2020)

2020
[5]

Bakkali-Hassani, C

B. Bakkali-Hassani, C. Maury, Y.-Q. Zou, ´E. Le Cerf, R. Saint-Jalm, P. C. M. Castilho, S. Nascimbene, J. Dal- ibard, and J. Beugnon, Realization of a townes soliton in a two-component planar bose gas, Physical Review Let- ters127, 023603 (2021)

2021
[6]

Ak¨ ozbek, A

N. Ak¨ ozbek, A. Iwasaki, A. Becker, M. Scalora, S. Chin, and C. Bowden, Third-harmonic generation and self- channeling in air using high-power femtosecond laser pulses, Physical review letters89, 143901 (2002)

2002
[7]

Berg´ e, S

L. Berg´ e, S. Skupin, C. K¨ ohler, I. Babushkin, and J. Her- rmann, 3d numerical simulations of thz generation by two-color laser filaments, Physical Review Letters110, 073901 (2013)

2013
[8]

V. Y. Fedorov and S. Tzortzakis, Extreme thz fields from two-color filamentation of midinfrared laser pulses, Phys- ical Review A97, 063842 (2018)

2018

Show all 12 references

[9]

Sukhinin, A

A. Sukhinin, A. B. Aceves, J.-C. Diels, and L. Arissian, Collapse events of two-color optical beams, Physical Re- view A95, 031801 (2017)

2017
[10]

Sukhinin, A

A. Sukhinin, A. Aceves, and J.-C. Diels, Two-color pulse collapse events in dispersive media, JOSA B36, G62 (2019)

2019
[11]

J. H. Marburger, Self-focusing: theory, Progress in quan- tum electronics4, 35 (1975)

1975
[12]

Yang, Newton-conjugate-gradient methods for solitary wave computations, Journal of Computational Physics 228, 7007 (2009)

J. Yang, Newton-conjugate-gradient methods for solitary wave computations, Journal of Computational Physics 228, 7007 (2009)

2009

Pith tools

Reviewed August 1, 2026 · model on record in the stance chip above.

[1] [1]

R. Y. Chiao, E. Garmire, and C. H. Townes, Self-trapping of optical beams, Physical review letters13, 479 (1964)

1964

[2] [2]

A. L. Gaeta, Catastrophic collapse of ultrashort pulses, Physical Review Letters84, 3582 (2000)

2000

[3] [3]

K. Moll, A. L. Gaeta, and G. Fibich, Self-similar optical wave collapse: observation of the townes profile, Physical review letters90, 203902 (2003)

2003

[4] [4]

Chen and C.-L

C.-A. Chen and C.-L. Hung, Observation of universal quench dynamics and townes soliton formation from modulational instability in two-dimensional bose gases, Physical Review Letters125, 250401 (2020)

2020

[5] [5]

Bakkali-Hassani, C

B. Bakkali-Hassani, C. Maury, Y.-Q. Zou, ´E. Le Cerf, R. Saint-Jalm, P. C. M. Castilho, S. Nascimbene, J. Dal- ibard, and J. Beugnon, Realization of a townes soliton in a two-component planar bose gas, Physical Review Let- ters127, 023603 (2021)

2021

[6] [6]

Ak¨ ozbek, A

N. Ak¨ ozbek, A. Iwasaki, A. Becker, M. Scalora, S. Chin, and C. Bowden, Third-harmonic generation and self- channeling in air using high-power femtosecond laser pulses, Physical review letters89, 143901 (2002)

2002

[7] [7]

Berg´ e, S

L. Berg´ e, S. Skupin, C. K¨ ohler, I. Babushkin, and J. Her- rmann, 3d numerical simulations of thz generation by two-color laser filaments, Physical Review Letters110, 073901 (2013)

2013

[8] [8]

V. Y. Fedorov and S. Tzortzakis, Extreme thz fields from two-color filamentation of midinfrared laser pulses, Phys- ical Review A97, 063842 (2018)

2018