Solitons or self-reinforcing solitary waves can emerge spontaneously in a physical system in which some energy is fed in, for instance as thermal energy or by an excitation with an electromagnetic wave, even if the excitation does not match exactly the soliton solution. Therefore, if a system possesses the necessary properties to allow the existence of solitons, it is highly likely that any large excitation will indeed lead to their formation [1] - [3] . The field of solitons and related nonlinear phenomena has been substantially advanced and enriched by research and discoveries in nonlinear optics [4] - [7] .

In our previous research [8] , we established the possibility of the existence of simultons (simultaneously propagating solitons at different frequencies) in the case of nonstationary Raman scattering with excitation of polar optical phonons under the conditions of the interaction of ultrashort pulses of exciting and Stokes radiation in nonlinear crystals. The relevance of this study is connected both with the fact that one can extract additional information on the optical characteristics of matter, and with the possibility of obtaining of ultrashort pulses.

The second topical problem in modern nonlinear optics is the production of coherent and frequency-tunable radiation in the far ultraviolet (UV) and infrared (IR). In these spectral areas, solid materials have broad absorption bands and this narrows down the application of nonlinear crystals for the generation of electromagnetic radiation. Possible ways of overcoming those difficulties are related with the utilization of nonlinear phenomena in gases and metal vapors. The resonant four-photon interaction (RFPI) in the case of two-photon resonance is one of them. Among the advantages of gases are the presence of narrow resonances and possibility of continuous variation of density, width of spectral line, length of the medium, etc. [9] - [11] . Ultrashort pulse propagation in the case of two-photon resonance was first examined in [12] where the two-photon self-induced transparency effect was predicted. This prediction was subsequently confirmed experimentally [13] and by numerical studies [14] . Third-harmonic generation (THG) in media exhibiting resonance behaviour has also attracted considerable attention [15] - [21] . However, RFPIs that are not frequency degenerate are of no less interest; they can be used to transfer the tuning of radiation from one range to another [22] [23] .

The present paper is devoted to the computer simulation of transition regimes of RFPI solitons in the case of two-photon resonance. The basic equations describing this process are given in Section 2 [12] [24] . The results of computer simulation are shown in Section 3. The stability of solitons is considered in Section 4.

Let us assume that two optical pulses with frequencies propagate in the nonlinear medium at the angles with respect to the z-axis. The value of is close to the frequency of resonant transition between levels 2 and 1 in the medium. The nonlinear interaction between and the medium results in parametric generation of and. The values of are considered to be in the transparent range of frequencies. We also assume that all electromagnetic waves have the same polarizations.

To find the system of equations that governs the processes of propagation of optical pulses with frequencies in the medium we take the standard system of equations for the amplitudes of probability of finding the system in state with energy [25]

where (2)

is the dipole moment of the transition; are the frequencies, wave vectors, real “slowly-varying amplitudes”, and phases of the interacting waves, respectively.

We next use (1) and the theory of perturbations [26] to find (the perturbation coefficient is of

order)

To obtain the system of equations for we introduce the expression (3) into the Equation (1), which becomes

where

.

The expression for the polarization induced by the superposition of nonlinear waves is defined by

We introduce (3) into (7) and find that the expression for the induced polarization becomes

where , ,

The system of Equations (4) and (5) can now be rewritten in terms of and n as follows

where , , is the maximum pulse amplitude; is the pulse width.

To make the system (8) - (10) complete we add Maxwell’s equations for the all real “slowly-varying amplitudes” and their phases. We obtain

where

, N is the number of molecules in cm³,.

We carried out the computer simulation of the system (8) - (10) and (11) - (18) for the following parameters of electromagnetic radiation and medium (gases): and [27] - [29] . The optical pulses on the pump and trigger frequencies were chosen to be of Gaussian shape. The accuracy of numerical results was based upon monitoring the conservation of energy of the system at every cross-sectional area in the medium. We considered the following trends: the phase locking and confinement; the length of medium needed for soliton formation; the distribution of energy during the transition regime; the energy limitations on solitons formation; the effect of slowing down of solitons in resonant systems; the relationship between the soliton speed and characteristics of incoming electromagnetic waves and nonlinear medium; the conditions leading to formation solitary wave at one frequency (instead of generation of sequence of them); the connection between the amplitude of soliton and parameters upon consideration. The space-time evolution of the normalized intensities (, , , and) is shown in Figures 1-4.

To investigate the stability of solitons we perform the summation of the Equations (12) (14) (16) and (18) for phases and transform them to the equations for. We assume that both processes occur at the conditions of synchronism, so that and. This condition is usually satisfied in gases [29] . Moreover, we also suggest that the phase differences are locked to or due to the nonlinear effects [30] . In this case (). Let be or, where are some small phase fluctuations. Finally, the modified system of (12) (14) (16) and (18) can be written in terms of as follows

where ,

is the simulation velocity,

The behaviour of the latter system is analyzed in terms of phase planes. As an example, Figure 5 shows the phase plane of the following system

at 11 different initial conditions for.

The space-time evolution of the optical pulses by using the computer simulation of transition regimes of four- photon resonant parametric processes in case of two-photon resonance is investigated. The computer simulation was based on application of the finite difference methods to the system of nonlinear equations modeling the foregoing interactions. It is shown that at certain boundary conditions (those result from the “area theorem” (see, e.g. [30] )) the incoming laser pulses at frequencies first generate new waves at, and then all become simultons of Lorentzian shape. It has also been shown that upon the conditions of phase locking (

or) and synchronism (and) in wide range of typical values of polarizabilities, the simultons are stable. These results could be useful for the applications related with designing the lossless communication systems using the tunable frequencies ranging from IR () to UV ().

M. Arif,J. Hussain,I. Hussain,S. Kumar,G. Bhati,Vladimir Feshchenko,Galina Feshchenko, (2016) Computer Simulation of Transition Regimes of Solitons in Four-Photon Resonant Parametric Processes in Case of Two-Photon Resonance. American Journal of Computational Mathematics,06,267-274. doi: 10.4236/ajcm.2016.63028