Cubic quintic Ginzburg Landau equation as a model for resonant interaction of EM field with nonlinear media

In this work, we investigate the possibility of approximating saturable nonlinearity, which is commonly used in complex Ginzburg–Landau equation (CGLE) for modelling resonant interaction of an electromagnetic field with nonlinear media, with cubic-quintic (CQ) nonlinearity. To validate the suggested approximations, we use variational method to estimate 2D analytical solutions of the CGLE with both saturable and CQ nonlinearity. The paper compares three ways to determine parameters of the CQ approximation and discusses the obtained results in terms of accuracy.


Introduction
The class of CGLEs represents a good model for describing a wide variety of phenomena in dissipative systems, such as superconductivity (Ginzburg et al. 1950), nonequilibrium phenomena (Aranson and Kramer 2002), phase transitions (Nato Advanced Study Institute 1975), binary fluid convection (Kolodner 1992), laser generated spatiotemporal dissipative structures, soliton propagation in optical fiber systems (Kodama and Hasegawa 1992), biological systems (Morales et al. 2015), etc.
The CQ CGLE, in particular, is often utilized to describe competitive effects of linear and nonlinear gain and loss. One of the reasons of high popularity of the CQ CGLE is This article is part of the Topical Collection on Advanced Photonics Meets Machine Learning.
the fact that there are many well developed methods for solving it both analytically and numerically. In the domain of fiber lasers, the model commonly used is the "master equation" derived by Haus (Haus et al. 1991(Haus et al. , 1992, which, in fact, represents the stationary CGLE with cubic nonlinearity. The coefficients that appear in the Haus model were related to the physical parameters in a rather phenomenological way (Haus et al. 1994). The cubic CGLE has been investigated in sufficient details in the one-dimensional case. Analytical solutions of dissipative structures in cubic CGLE have been obtained, but all of them are unstable except the class of solutions with arbitrary amplitude in an "exotic" medium with zero linear loss/gain (Akhmediev et al. 1996).
To achieve the stabilization of the dissipative structures (solitons) in the CGLE model, higher-order nonlinear (quintic) terms have been introduced. The obtained equation is referred to as the quintic or CQ CGLE equation. The one-dimensional CQ CGLE, in the case of special relations between its parameters, has analytical solutions in the form of stable solitons with a fixed or arbitrary amplitude (Akhmediev and Afanasjev 1995;Akhmediev et al. 1997Akhmediev et al. , 1998. In the general case, the one dimensional and multidimensional CQ CGLE do not have known exact solutions and have been studied numerically, as well as using different analytical approximations, including the variational method (Skarka and Aleksić 2006;Skarka et al. 2014Skarka et al. , 2016Aleksic et al. 2007Aleksic et al. , 2012aAleksic et al. , b, 2014Aleksic et al. , 2015.
After the introduction of the cubic "master equation" (Haus et al. 1991(Haus et al. , 1992, there have been derived many refinements in the form of CQ CGLE, as a model of resonant interaction of electromagnetic field with nonlinear media. These refinements have been done assuming small intensities, which allows the expansion of the Bloch equations in the series of intensities.
In this work, we further discuss the possibility of using CQ CGLE for modelling resonant interaction of an electromagnetic field. We compare three CQ approximations, including the widely used series expansion, and the CGLE with saturable nonlinearity. More precisely, the comparison is performed in terms of the parameters of stationary two-dimensional fundamental solutions, obtained using variational approach.

Model
We consider the response of a medium where active (amplifier) and passive atoms (saturable absorber), with almost identical transition frequencies, are embedded in the background matrix (for example, in a waveguide). Due to quick relaxation of the upper energy levels (three or four), we can reduce the model to an effective two-level scheme. When the concentration of the resonance centers is small, the polarization can be considered as an additive quantity: consisting of the polarization of active (P ( ) ) , passive (P ( ) ) atoms, and non-resonant polarization of the background matrix (P (m) ).
From the Maxwell-Bloch equations, when the electromagnetic pulse duration is much greater than the relaxation time of the upper working level, we obtain the dimensionless CGLE for the electrical field: (1) P = P ( ) + P ( ) + P (m) where and ̄= , ̄= . The constants , are dimensionless detunings between the gain (amplification) and loss (absorption), with the respect to the spectral line center and the frequency of empty cavity mode. The coefficients and represent the linear amplification of active atoms and the absorption of passive atoms, (normalized to the absorption coefficient of the background medium), respectively (Fedorov et al. 2000). The linear part of the polarization determines the linear losses, temporal dispersion of the media, and spectral filtering ( 1 ). Under conditions of resonant interactions, when → 0 and → 0 , term N(I) in Eq. (2) can be neglected.
The function G(I) has a maximum for if the parameters satisfy the following relation: Outside this region the function G monotonically either decreases, when > b , or increases, when < ∕b . To establish the energy balance and the possibility of the existence of localized structures, the function G(I) should be positive in a finite interval of intensity, i.e. G(I m ) > 0 , and G(I = 0) < 0 , in order to have a stable background. These two conditions in the parameter space become: Under the previous conditions, function G(I) has two zeros: as shown by the middle curve in Fig. 1a. The lower curve ( ∕b < < min ) does not have positive regions, thus the function G provides only loss. The fact that in the upper curve ( max < < b) we have G(0) > 0 means the linear gain is bigger than the linear loss, which causes an unstable background. The aforementioned properties of the lower and upper curve prevent formation of stable localized structures. In order to simplify analytical estimation of Eq.
(2), we consider three approximations of G, by using a second-order polynomial: which satisfy the condition G a (0) = G(0) i.e.
The first approximation ( G 1 ) is a series expansion, the second ( G 2 ) reaches maximum value at I m , and the third ( G 3 ) has the value of zero at I 0+ , and G(I m ) at I m . Thus: 1. Approximation G 1 : the series expansion of Eq. (4) 2. Approximation G 2 : dG 2 (I m )∕dI m = 0 and G 2 (I m ) = G(I m ) 3. Approximation G 3 : G 3 (I 0+ ) = G(I 0+ ) = 0 and G 3 (I m ) = G(I m ) Figure 1b shows the saturable nonlinearity (G) and its three approximations ( G 1 , G 2 , G 3 ) as functions of intensity I, for parameters = 1 and = 1.8 . Clearly, the most widely used approximation G 1 deviates the most from G, while G 2 and G 3 better resemble G. In the next section we analyze each approximation in more details.
The coefficients , a and a account for linear loss, cubic gain, and quintic loss, respectively. If the loss saturates faster than the gain (b > 1) , then coefficients a and a are positive.
(10) G a (I) = − + a I − a I 2 , a = 1, 2, 3 Fig. 1 a G as a function of intensity I for different regions of parameter . b G(I) and its three approximations for = 1 and = 1.8

Variational approach
In order to validate the approximations, we use the variational method to solve the twodimensional CGLE: with saturable nonlinearity and its CQ approximations, in the form of localized fundamental solitons. We adopt the trial function in the form of a Gaussian: where the real variational parameters are the amplitude A = A(z) , width R = R(z) , wavefront curvature C(z), and phase = (z). Skipping technical details, along the lines of Skarka and Aleksić (2006), for the saturable nonlinearity of Eq. (4) we obtain the following system of the first-order evolution equations for the parameters of ansatz in Eq. (16): From the fixed points of Eqs.(17)- (19): we obtain the steady-state solutions, which are determined by the following system of equations: where e = ∕ √ 1 + 2 and Li 2 (x) is the polylogarithm function. In the case of CQ approximations Eq. (10), the stationary values of the trial function parameters are given by the following system of equations ( a = 1, 2, 3):  G 2 and G 3 . The series expansion G 1 , which is widely used in many research papers, represents a good approximation only for small values of the intensity. Approximation G 2 shows better accuracy for slightly larger values of the intensity, while G 3 gives values fairly close to the saturable nonlinearity over the whole intensity domain.
From Eqs. (22) , (25) and (26) we can see that in the case when the wave front curvature (C and C a ) approaches zero, the radius (R and R a ) approaches infinity. Therefore, the nonzero value of the wave front curvature is the essential condition in order to have a finite radius, i.e. the existence of localized solutions. The zero values of C and C a are marked with circles in Fig. 2.

Conclusion
A common model for resonant interaction of electromagnetic waves with nonlinear media is the CGLE with saturable nonlinearity. Due to limitations in the analytical analysis of that kind of a model, we approximate the saturable nonlinearity with CQ nonlinearity, whose analytical and numerical aspects have been extensively stuided.
Using the variational method, we analyze three different second-order polynomial approximations of the saturable nonlinearity. The obtained results show the differences of the steady-state solutions between different approximations and the saturable nonlinearity, and depending on the application, suggest the most accurate approach. The most adequate approximation is the second-order polynomial, which has the same value at the linear limit and the common distant zero with saturable nonlinearity.
In the parameter space of the resonant medium, the region of possible approximations is determined. The parameters of the polynomial approximation are related to the physical parameters of the resonant medium.