“…Let ψ x,y = (L d /L s I 0 ) 1/2 e ±iβz/4 E x,y , where L d = k x a 2 ≈ k y a 2 is the diffraction length for close wavenumbers of two components k x , k y , β = 2|k x −k y |L d is the normalized phase mismatch, L s = 2c/(ω 0 n 2 I 0 ) is the self-action length for nonlinear coefficient n 2 , a and I 0 are the characteristic transverse scale and intensity, respectively. Then, ψ x,y are governed by the system [25]:…”