Maker and Terhune were first who demonstrated the CARS technique . In  was also shown that the efficiency of the CARS generation is the function of the third-order susceptibilities which were extensively studied in the experiments of the nonlinear properties of solids and liquids   . Since then CARS spectroscopy has become a powerful technique in many fields of knowledge such as physics, biology, chemistry, healthcare, etc.     . The CARS imaging proved its efficiency in cancer diagnosis as well  . Begley et al. were among the researches who summarized the important advantages of vibrational spectroscopy based on nonlinear anti-Stokes generation   . This technique went to the next level when the ultrashort laser pulses resulted in the possibility of the coherent excitation of multiple Raman modes   . For example, the method of nonstationary (time-domain) CARS spectroscopy permitted direct observation of vibration dephasing in an ensemble of atoms or molecules or even in the simplest system-molecular hydrogen . The further increase in efficiency would result from simultaneous propagation in the medium of all interacting waves. In    we considered the cases of Raman scattering by polaritons in dipole-active crystals. In this paper, we considered the theoretical modeling of the processes of nonstationary CARS by polaritons in dipole-active crystals.
2. Basic Principles and Equations
In this paper, we consider the nonlinear interaction of four electromagnetic waves: anti-Stokes, Stokes, pump (laser), and polariton. Those waves are assumed to be linearly polarized plane waves. It is also assumed that the nonlinear medium takes a form of a layer bounded by the planes z = 0 and z = L. The pump wave
propagates along the z-axis. The subscripts a, l, s, and p denote the anti-Stokes, pump (laser), Stokes, and polariton wave fields, are the frequencies, and are the refractive indices and the wave vectors in the unpumped medium, and are the real unit vectors of electromagnetic fields. The nonlinear medium is assumed to be nonmagnetic and transparent at the frequencies . We use the anti-Stokes, Stokes, and polariton fields in the form
where ; ; ; .
In the process of CARS, the nonlinear interaction of two electromagnetic waves results in the generation of anti-Stokes and polariton waves. The system of shortened equations for the amplitudes is obtained from Maxwell’s equations by using the standard approximation of slowly-varying amplitudes  and takes the form
where are the corresponding tensor contractions of non-resonance quadratic and cubic nonlinear polarizabilities with unit vectors of polarization of interacting waves; is the non-resonance part of dielectric permeability at frequency ; are z-components of velocities of waves on ; is the wave mismatch between the pump, polariton, and anti-Stokes waves.
Given the strong polariton absorption we have 
so that we can neglect in (8) the terms with the derivatives after which this equation yields
If we insert the obtained expression for the amplitude of polariton wave in (5)-(6), we get a system of 3 differential equations for as follows:
where , , , and .
The system (11)-(13) can be simplified if we use new variables
and . (15)
The system (11)-(13) in terms of can be written as follows:
And, finally, if we assume a “week” wave mismatch at Stokes and anti-Stokes frequencies, that is
then the final system of equations simulating CARS can be expressed as
3. Asymptotic Solutions in a Form of Simultaneously Propagating Waves at Frequencies
Since we will conduct the numerical analysis of the system (20)-(22) we bring it to unitless form first. To do that we multiply both the left and right part of each equation by the factor ( and are the peak amplitude and characteristic duration of the pump, ). After that, the system (20)-(22) can be reduced to
where , , , ,
, , ,
We are looking for stationary solutions as
and , (26)
where ; is the velocity of simultaneously propagating waves at the frequencies ; and are the real amplitudes and phases of the waves, respectively. Such a standard procedure of presenting the complex amplitudes of waves in terms of real and imaginary parts results in duplication of the system of (23)-(25):
where , .
If we introduce the amplitude of simultaneously propagated waves as
we could reduce the system above to
The system (34)-(35) can be further simplified as follows
where , .
We can reduce the number of equations by using the integral of motion
where , .
If we express the phase as the function of Q in (37) we get
The integral on the left can be found as follows:
where is the elliptic integral of the first kind with the parameter , is the inverse hyperbolic sine function.
In Figure 1 it is shown that the solution of (37) and (38) exists in the form of pulses. The duration of those pulses can be easily evaluated as follows: first, we assume, that (in the next topic it is shown that where g is the gain factor of Raman scattering) so that the coefficient . The typical values of the gain factor in crystals are of order 10−3 cm/MW . Hence, if we consider the pump of the intensity of 102 - 103 MW and , then .
4. Gain Factor g
To show that the system of Equations (5)-(8) is consistent with experimental results for CARS by polaritons we consider the stationary solutions of the coupled wave equations in the constant pump approximation. The system for Stokes and anti-Stokes (20), (22) under the above suggestions can be expressed as
Figure 1. Q versus x.
After introducing new variables
the system of differential equations of the first order can be readily transformed to the single differential equation of the second order (for example, for ) as
We solve this equation by adopting a trial solution for in the form
where g represents a gain factor. Then we substitute (46) into (45) to obtain the approximate value for g as
(here we assumed that the pump was strong enough to provide ).
Finally, the expression for g can be reduced to
As the experimental data for this gain, we used the following : pulse width of the pulsed Ar+ laser ; 30 ps, the peak output power ≈ 2.5 kW, the wavelength was 514.5 nm, the cross-section ≈ 10−18 cm−2, , and . In  the nonlinear medium was zinc blende ZnS, in which the polariton frequencies were in the rage 200 - 400 cm−1. Both the experimental results for the gain factor in  and calculations based on (48) have resulted in g ≈ 1.
In this paper, we have found the system of differential equations that model the process of coherent anti-Stokes Raman scattering by polaritons in crystals. We have also found the asymptotic solutions of that system that correspond to the simultaneous propagation of all waves participating in the process of Raman scattering. And, lastly, we showed that the value of such an important feature as the gain factor resulted from that system is consistent with the experimental results.
 Maker, P.D. and Terhune R.W. (1965) Study of Optical Effects Due to an Induced Polarization of Third Order in the Electric Field Strength. Physical Review, 137, A801-A818.
 Tolles, W.M., Nibler, J.W., McDonald, J.R. and Harvey, A.B. (1977) A Review of the Theory and Application of Coherent Anti-Stokes Raman Spectroscopy. Applied Spectroscopy, 31, 253.
 Zheltikov, A.M. (2000) Coherent Anti-Stokes Raman Scattering: From Proof-of-the Principle Experiments to Femtosecond CARS and Higher-Order Wave-Mixing Generalizations. Journal of Raman Spectroscopy 31, 653-667.
 Shen, Y., Voronine, D.V. Sokolov, A.V., and. Scully, M.O. (2015) A Versatile Setup Using Femtosecond Adaptive Spectroscopic Techniques for Coherent Anti-Stokes Raman Scattering. Review of Scientific Instruments, 86, Article ID: 083107.
 Potcoava, M.C., Futia, G.L., Aughenbaugh, J., Schlaepfer, I.R. and Gibson, E.A. (2014) Raman and Coherent Anti-Stokes Raman Scattering Microscopy Studies of Changes in Lipid Content and Composition in Hormone-Treated Breast and Prostate Cancer Cells. Journal of Biomedical Optics, 19, Article ID: 111605.
 Harris, I.E. and McIlwain, M.E. (1981) Coherent Anti-Stokes Raman Spectroscopy in Propellant Flames. In: Capellos, C. and Walker, R.F., Eds., Fast Reactions in Energetic Systems, Springer, Dordrecht, 473-484.
 Mukamel, S. (1990) Femtosecond Optical Spectroscopy: A Direct Look at Elementary Chemical Events. Annual Review of Physical Chemistry, 41, 647-681.
 Nelson, K.A. and Ippen, E.P. (1989) Femtosecond Coherent Spectroscopy. In: Prigogine, I. and Rice, S.A., Eds., Advances in Chemical Physics, John Wiley & Sons, Inc., Hoboken, 1-35.
 Dyakov, Y.S., Krikunov, S.A., Magnitskii, C.A., Nikitin, S.Y. and Tunkin, V.G. Time-Domain Coherent Anti-Stokes Raman Spectroscopy for Hydrogen Gas in the Dicke Narrowing Region. Zh. Exp. Teor. Fiz, 84, 2013-2025.
 Feshchenko, G. and Feshchenko, V. (2015) Computer Simulation of Transition Regimes of Solitons in Stimulated Raman Scattering with Excitation of Polar Optical Phonons. American Journal of Computational Mathematics, 5, 336-344.
 Feshchenko, V. and Feshchenko, G. (2018) Nonstationary Stimulated Raman Scattering by Polaritons in Continuum of Dipole-Active Phonons. Journal of Applied Mathematics and Physics, 6, 405-417.
 Feshchenko, V. and Feshchenko, G. (2019) Nonstationary Stimulated Raman Scattering by Polaritons in Cubic Crystals. Journal of Applied Mathematics and Physics, 7, 2122-2129.
 Cheng, Y.C., Jin, C.Q., Gao, F., Wu, X.L., Zhong, W., Li, S.H., and Chu, P.K. (2009) Raman Scattering Study of Zinc Blende and Wurtzite ZnS. Journal of Applied Physics, 106, Article ID: 123505.