American Journal of Computational Mathematics
Vol.05 No.03(2015), Article ID:59545,8 pages

Computer Simulation of Transition Regimes of Solitons in Stimulated Raman Scattering with Excitation of Polar Optical Phonons

Galyna Feshchenko1, Vladimir Feshchenko2

1Vanier College, Montreal, Canada

2Dawson College, Montreal, Canada


Copyright © 2015 by authors and Scientific Research Publishing Inc.

This work is licensed under the Creative Commons Attribution International License (CC BY).

Received 7 August 2015; accepted 7 September 2015; published 11 September 2015


The system of nonlinear equations modeling the process of nonstationary stimulated Raman scattering (SRS) in noncentrosymmetric crystals for the waves on laser, Stokes, polariton, and phonon frequencies is investigated by using the numerical methods. The general case for amplitudes of waves that resulted in doubling of the number of equations is considered. It is shown that the application of the methods of finite differences to the computer simulation of transition regimes is completely consistent with the analytical results found for the asymptotical solutions in form of solitons. The obtained results also indicate that the laser pulses of Gaussian shape appearing at the boundary of nonlinear medium tend to become solitons of Lorentzian shape. It was also found that the formation of solitons occurs when the vibrations of optical phonons and that of electromagnetic wave were either in or out of phase. It is shown that all electromagnetic waves entering the medium with different speeds become solitons having the same speed. In the second part of the paper we considered the computer simulation of soliton stability with respect to small (weak) perturbations of all interacting waves. In the present paper we considered the case of evolution of those disturbances in the vicinity of peaks of solitons. The numerical analysis showed that in wide range of parameters the solitons were stable.


Finite Differences Method, Abbreviated Maxwell’s Equations, Stimulated Raman Scattering, Stability

1. Introduction

A fast progress in producing ultrashort pulses stimulates the constantly growing interest in obtaining the additional information on both the optical characteristics of matter (the longitudinal and transverse relaxation times, the constants of the transition, etc.) and the possibility of creation of super short pulses propagating in a medium without energy losses [1] [2] . One of the effects that results in solitons is SRS under the conditions of the interaction of ultrashort pulses of exciting and Stokes radiation in a medium [3] - [5] . The effect of SRS is modeled by the system of abbreviated Maxwell’s equations for the amplitudes of the electromagnetic fields and the equations of motion for the averaged polarization [6] . In many cases due to the complexity only the asymptotic behavior of those systems can be found [7] . However, to make clear, whether the SRS pulses go over to the soliton regime, it is necessary to solve the complete set of nonstationary equations which describes the temporal evolution of all interacting waves from the boundary of the medium. It is also important to consider the stability of the SRS solitons [8] . Using different numerical methods in computer simulation of space-time evolution of the interacting waves is the effective approach in analysis of the equations analytical solutions of which could not be found [9] - [11] . In this paper we use the assumption that solitons can exist if, at least, three conditions are satisfied: 1) the asymptotical solution in form of solitons is found; 2) the results of transition regime are consistent with the ones found in the previous part; 3) the obtained results are stable with respect to small (weak) perturbations. In our previous paper we established the possibility of existence of the solitons in nonstationary SRS with excitation of polar optical phonons [12] . In the present work it has been shown that the results of the computer simulation of the system of nonlinear equations (see [12] ) are consistent with the asymptotic solutions. The computer simulation was provided in Sections 3 and 4. The stability of soliton propagation was studied in Section 5.

2. Basic Equations

Let us assume that the pump laser wave with frequency and Stokes wave excite both the polariton radiation and the polar optical phonon, in the vicinity of which falls ( the vibration is supposed to be dipole-active and solitary). The two optical impulses with propagate in noncentrosymmetric crystal at the angles with respect to z-axis (perpendicular to the crystal plane). The vibration, due to its combinational activity, is excited by the cubic nonlinear polarization and electric field on. The electric strengths are:



where: are the wave vectors in linear theory; are the unit vectors ( we assume that only waves of certain fixed polarization are effectively involved into interaction); are the amplitudes of interacting waves. The phonon vibration is:


where is collinear with. We neglect the pumping of energy into the anti-Stokes and higher Stokes components, since we consider the case in which the intense fields of the exciting (laser) and first Stokes waves are present at the input of the medium. The basic system of nonlinear equations was found in [7] (the first three equations are abbreviated Maxwell’s whereas the fourth one is the equation that simulates the behavior of the nonlinear medium):





where: and are the corresponding tensor contractions of nonresonance quadratic and cubic nonlinear

polarizabilities with unit vectors of polarization of interacting waves;; is the oscillator

strength of the oscillation; is the ratio of effective charge to effective mass for given vibration;

is the nonresonance part of dielectric permeability on frequency;;;

is the tensor of combinational scattering per one elementary cell ;;; is constant of phonon attenuation; -z-components of velocities of waves on.

The asymptotic solution of the system (4)-(7) in form of solitons was found in [12] . Soliton amplitudes and phases were related with as follows

, (8)

, (9)

where the real amplitudes and phases are considered to be functions of only v is the soliton velocity; are the indexes of refraction in linear theory. The soliton solution was found to be of Lorentzian shape:






3. Finite-Difference Equivalent of Basic System of Equations

In this paper, we apply the method of finite differences to show that space-time evolution of impulses is consistent with asymptotic behaviour in form of solitons. For the time derivative the forward difference was applied [11] :


For the space derivative the backward difference was applied:


The system (4)-(7) can be rewritten as follows










The chosen mesh was

4. Transition Regime of the Stimulated Raman Scattering Solitons

The pulses on laser and Stokes frequencies were chosen to be of Gaussian shape (their amplitudes were proportional to (the coefficient was 0.159155; this was done in order to minimize the transition re-

gime and to get the area of each pulse close to). The two-dimensional arrays A for the amplitudes of waves were A [10] [1000] which resulted in 10,000 points for each wave. The conservation of energy (in our case the area under the pulse shape) was monitored by calculating the time integrals for each cross-sectional area in a medium. Since we suggested that the normalized areas of incoming waves (area divided by the amplitude) were equal to, all those integrals showed the evolution of pulses of Gaussian shape to those of Lorentzian shape with normalized area of ( [13] ). The results are shown in Figure 1.

Figure 1. The space-time evolution of the normalized intensities of laser, Stokes, and polariton pulses.

5. Stability of the Stimulated Raman Scattering Solitons

Let the new amplitudes (amplitudes with disturbances) in unitless form be



where: are the unitless amplitudes of waves in presence of disturbance ();

are the small disturbances; -is the peak amplitude of the incoming laser pulse;

;;;; c is the speed of light in vacuum; is the characteristic time (=

1 fs). When we substitute (22) and (23) into the system (4)-(7) we get the system of equations that can be significantly simplified if we consider the evolution of disturbances in the vicinity of the maximum value of soliton amplitude that is in the case when. In that case the simplified system of equations (the disturbances are considered to be small compared to the amplitudes of solitons) becomes that of four linear equations of second order:









The system of equations (24)-(27) was solved numerically. The results are presented in Figure 2.

Figure 2. The time evolution of normalized disturbances (red solid line; green dashed line, blue dotted line) at

6. Conclusion

The system of eight nonlinear equations describing the SRS with excitation of polar optical phonons was numerically solved. The results of the computer simulation are completely consistent with the asymptotical solutions in form of solitons. It is numerically shown that with SRS on dipole-active phonons the process of propagation of pulses of Gaussian shape on laser and Stokes frequencies results in formation of solitons of Lorentzian shape at all frequencies of interacting waves. Those solitons are stable in wide range of parameters that characterize both the electromagnetic waves and medium.

Cite this paper

GalynaFeshchenko,VladimirFeshchenko, (2015) Computer Simulation of Transition Regimes of Solitons in Stimulated Raman Scattering with Excitation of Polar Optical Phonons. American Journal of Computational Mathematics,05,336-344. doi: 10.4236/ajcm.2015.53031


  1. 1. Rolle, J., Berge, L., Duchateau, G. and Skupin, S. (2014) Filamentation of Ultrashort Laser Pulses in Silica Glass and KDP Crystal: A Comparative Study. Physical Review A, 90, 1-8.

  2. 2. Chen, Z., Segev, M. and Christodoulides, D.N. (2012) Optical Spatial Solitons: Historical Overview and Recent Advances. Reports on Progress in Physics, 75, 1-20.

  3. 3. Prokopovich, I.P. and Khrushchinskii, A.A. (1997) Highly Efficient Generation of Attosecond Pulses in Coherent Stimulated Raman Self-Scattering of Intense Femtosecond Laser Pulses. Laser Physics, 7, 305-308.

  4. 4. Prokopovich, I.P. and Peatross, J. (1999) Direct Shaping and Amplifying of High-Intense Single Attosecond Pulses from High-Intense Femtosecond Optical Pulses in Inert Gases. Laser Physics, 9, 588-591.

  5. 5. Kaplan, A.E. and Shkolnikov, P.L. (1996) Subfemtosecond Pulses in the Multicascade Stimulated Raman Scattering. Journal of the Optical Society of America B, 13, 347-354.

  6. 6. Makhviladze, T.M. and Sarychev, M.E. (1976) Soliton Regimes of Stimulated Raman Scattering. Soviet Physics— JETP, 44, 471-477.

  7. 7. Makhviladze, T.M., Sarychev, M.E. and Shelepin L.A. (1975) Raman Scattering of Light by an Excited Medium. Soviet Physics—JETP, 42, 255-261.

  8. 8. Kivshar, Y. and Agrawal, G. (2003) Optical Solitons: From Fibers to Photonic Crystals. Academic Press, Waltham.

  9. 9. Ames, W. (1992) Numerical Methods for Partial Differential Equations. 3rd Edition, Academic Press, Waltham.

  10. 10. Farlow, S. (1993) Partial Differential Equations for Scientists and Engineers. Dover Books, Mineola.

  11. 11. Kreyszig, E. (1993) Advanced Engineering Mathematics. 7th Edition, John Wiley & Sons Ltd., Hoboken.

  12. 12. Marchevskii, F.N., Strizhevskii, V.L. and Feshchenko, V.P. (1984) Soliton Generation by Stimulated Raman Scattering with Excitation of Polar Optical Phonons. Quantum Electronics, 11, 277-282.

  13. 13. Boyd, R. (2008) Nonlinear Optics. 3rd Edition, Academic Press, Waltham.