In general, the pulses (including solitons) are randomly polarized. However, it is desirable to have pulses with a certain (better yet, predominantly defined polarization). It stands to a reason that the polarization could also be a carrier of the information about processes that took place during the interaction. Especially it is important in the case of optical fibers . One of the features of the propagation of solitons in the randomly birefringent optical fibers is that they are randomly polarized whereas there is a need for the solitons with well-defined polarization. In , it was both theoretically and experimentally shown that in special cases when circular birefringence was induced when the fibers were twisted, the random linear birefringence was mitigated. In , it was shown that an induced process of modulational instability may be exploited for the generation of THz train of vector dark solitons. The technique of frequency-resolved optical gating was used to completely characterize the intensity and phase of the dark solitons. In the above article, the authors analyzed the experiments that showed that with dual-frequency orthogonal polarization, one may achieve the simultaneous suppression of modulational instability. The second nonlinear phenomenon that is significant for the generation of ultrashort optical pulses is stimulated Raman scattering (SRS). In , by using a nanosecond laser as the pump source, polarization modulation instability (PMI) and SRS were experimentally investigated in an optical fiber. It was shown that such investigation on the PMI gain in the As2S5 optical fiber may provide a reference for the development of fiber amplifiers and fiber lasers. The detailed analysis of the polarization of the solitons generated by PMI was provided in  in which was found that the circular polarization of the pump pulse in the twisted fiber produces circularly polarized solitons with a high grade of polarization whereas in the fiber without twisting, the soliton polarization was random. The calculation of a polarization-dependent effective SRS gain factor considering the random birefringence character of the fiber and the relative mismatch between the continuous pump and the signal was carried out in . There were also achievements in the theory of solitons and SRS . The numerical simulation of the transient stimulated Raman scattering based on using the inverse scattering transform technique is reported in . A new type of solitons—dissipative Raman solitons—was found based on the nonlinear Ginzburg-Landau equation in . In past years, there was also significant progress related to the development of the micro-resonators which convert laser light into a series of ultrashort pulses traveling around the resonator’s circumference . These pulses can propagate while maintaining their shape, therefore, becoming solitons. One of the promising applications of solitons is, of course, spectroscopy. The comparative investigation of characteristics of spontaneous and stimulated Raman scattering (SRS) in different alkali-earth tungstate and molybdate crystals at both high and low-frequency anionic group vibrations was presented, for example, in . It has been found that, among these crystals, the SrMoO4 and SrWO4 crystals are the most perspectives for SRS generation on both stretching and bending modes of internal anionic group vibrations with the strongest SRS pulse shortening under synchronous laser pumping. The significant progress in applying the methods of SRS and CARS (Coherent Anti-Stokes Raman Spectroscopy) in recent years was achieved in biology and medicine  - . For example, in  it was reported a new procedure for alignment-free frequency modulation SRS utilizing polarization encoding since the SRS microscopy is a powerful method for imaging molecular distributions based on their intrinsic vibrational contrast. The successful application of polarization-sensitive SRS imaging in treating skin, nail, and fungal infections was reported in .
However, in our opinion, there is a need for studying the processes of the generation of solitons in nonlinear processes with predetermined polarization that would allow generating the stable ultrashort pulses traveling through the medium not only at different frequencies but different polarization (for example, the soliton at one frequency but with two perpendicular polarizations). Such creation of ultrashort stable pulses could significantly increase the resolution of SRS microscopy of imaging of the molecular vibrations. In this paper, we theoretically considered the formation of solitons of such type in transient SRS by polaritons  - . Also, in our article principal attention was paid to a calculation and analysis of the gain of both stationary SRS and spontaneous Raman scattering. We considered the case of nonstationary SRS in cubic crystals and showed that the theory developed is consistent with experimental results.
2. Basic Principles and Equations
In this paper, we carry out our analysis assuming that the pump field is a linearly polarized plane electromagnetic wave. It is also assumed that the nonlinear medium takes the form of a layer bounded by the planes z = 0 and z = L. The pump wave
propagates along the z-axis. The subscripts l, s, and p henceforth denote the pump (laser), Stokes, and polariton wave fields; are the frequencies, n and are the refractive indices and the wave vectors in the unpumped medium, and are the real unit vectors of electromagnetic fields. The medium is nonmagnetic and transparent at frequencies . We use the Stokes and polariton fields in the form
where: , , , , , , , , .
The longitudinal component of the Stokes wave can be neglected, but this cannot be done for the polariton wave in the phonon region. It has been shown in  that with a further advance into this region all three amplitudes first become comparable, after which becomes dominant, provided, of course, the excitation of the longitudinal waves is allowed by the selection rules. The phase shift of the polariton wave is determined by the vector and not by ( , , is the dielectric constant at the frequency ).
The fields are interrelated via the nonlinear part of the polarization . The latter quantity has at the frequencies the following forms
where , , , .
The shortened equations for the amplitudes are obtained from Maxwell’s equations by the standard procedure  and take the form
Note, that in (4) .
Given the strong absorption we have 
and we can, therefore, neglect in (6) and (7) the terms with the derivatives after which these equations yield
Substituting the obtained expressions in (6) and (7), we arrive at a system of two differential equations for
, , ,
is the oscillator strength of the o-f transition.
3. Asymptotic Solutions: Polarization Solitons
To do that we bring the system (9, 10) to unitless form.
where , , , , , , , , , , , , , is the characteristic time related to the laser field (pump).
We are looking for stationary solutions as
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 (13)-(15):
where , , .
We are looking for the solitary (asymptotic) solution for solitons at as following:
where , and .
The introduction of Q allows to reduce the system of nonlinear equations (17)-(22) to
where , ,
, , . The space-time evolution of the normalized intensities is shown in Figure 1.
Figure 1. The space-time evolution of the normalized intensities of the pump and Stokes pulses for: (a) and , (b) and , and (c) and , respectively.
4. Gain Factor
Now we show that the system of Equations (14) and (15) is consistent with the experimental results presented, for example, in . To do that we first bring the system of Equations (14) and (15) to unitless form and change the variables z, t to variables , (we assume that and consider the preset-pump approximation):
where , , , , , , , , , is the characteristic time related to the laser field (pump).
The theoretical consideration of the gain factor for SRS by polaritons is based on the modeling of the quasi-stationary solutions of the coupled wave equations for the different polarizations of the Stokes (for the complete analysis see  ). Therefore, we seek the solutions of (25) in the form , assuming and to be independent of z. We then obtain the system of algebraic equations for . Choosing in a plane perpendicular to a two-dimensional coordinate system with axes along the unit vectors, we represent the equations for in the form of a tensor relation
Equating the determinant of the system (26) to zero, we obtain the solutions for
We will need the explicit expressions for the tensors and . They can be found within the framework of the microscopic theory in the dipole approximation based on the perturbation theory states . The resultant expressions are
The summation in (28) and (29) is over all dipole-active phonons, the frequencies of which are equal , where are the attenuation constants. For example in a cubic crystal, the dipole-active phonons are triply degenerate  so that the number of the mutually degenerated oscillations we introduce the index ( is a triad of real unit vectors denoting the vibrations along the edges of the unit cube. Furthermore, is the dipole moment of the transition 0-fv for the unit cell with its volume ; is the tensor of the phonon spontaneous scattering per cell ; is the number of cells in the crystal. The tensor represents the contribution to the remote electronic states. The tensor determines the contribution due to the electronic states as well. It is convenient to represent the tensors (28) and (29) in the simplified form as follows
where , , is the Raman differential cross-section per unit cell (cm−1/sr).
We introduce the principal axes of the tensor as a whole. If we denote its principal values as we obtain from (26) . Finally, we introduce the gain which can be expressed as
where , , is the pump intensity, are the principal values of the tensor , is the scattering angle (the angle between and ( , )).
Formula (31) denotes two gain coefficients for Stokes waves polarized along .
To verify (31), we were using the parameters of crystals widely used in optical display and storage, optical communication network, optical detection, etc. such as ZnO  -  and ZnS    . Figure 2 shows the intensity as a function of the polariton frequency in zinc blende ZnS in the range 200 - 400 cm−1. The red dots represent the experimental points .
Figure 2. The graph of intensity (in arbitrary units) versus the polariton frequency (in 1/cm) in ZnS in the range 200-400 1/cm.
In this paper, we theoretically showed that in the case of nonstationary SRS by polaritons, there is a possibility of occurrence of simultaneously propagating ultrafast stable pulses (simultons) not only at different frequencies but with different polarizations as well. This can be used in optoelectronics when creating polarization filters.
 Rosas, A.F., Hernández, O.D., Arceo, R., Santos, G.J.E., Vázquez, S.M., álvarez, E.R., Flores, C.I.E. and Kuzin, E. (2018) Polarization Properties, Nonlinear Optics—Novel Results in Theory and Application.
 Millot, G., TchofoDinda, P., Seve, E. and Wabnitz, S. (2001) Modulational Instability and Stimulated Raman Scattering in Normally Dispersive Highly Birefringent Fibers. Optical Fiber Technology, 7, 170-205.
 Cheng, T.L., Wang, Q.M., Yan, X., Wang, F., Zhang, X.N., Li, S.G., Suzuki, T. and Ohishi, Y. (2020) Experimental Investigation of the Polarization Modulation Instability and Stimulated Raman Scattering in a Chalcogenide Optical Fiber. Journal of Applied Physics, 128, Article ID: 193103.
 Flores-Rosas, A., Kuzin Evgeny, A., Arceo, R., Díaz-Hernán, O., Ruiz-Pérez, V.I. and Escalera-Santos, G.J. (2017) Polarization Properties of Vector Solitons Generated by Modulation Instability in Fiber with Circular Birefringence. Optical Engineering, 56, Article ID: 036115.
 Ortiz-Mora, A., Rodríguez, P., Díaz-Soriano, A., Martínez-Muñoz, D. and Dengra, A. (2020) Method of Moments Optimization of Distributed Raman Amplification in Fibers with Randomly Varying Birefringence. Photonics, 7, 86.
 Frank, M., Smetanin, S.N., Jelinek, M., Vyhlidal, D., Shukshin, V.E., Ivleva, L.I., Dunaeva, E.E., Voronina, I.S., Zverev, P.G. and Kubecek, V. (2019) Stimulated Raman Scattering in Alkali-Earth Tungstate and Molybdate Crystals at Both Stretching and Bending Raman Modes under Synchronous Picosecond Pumping with Multiple Pulse Shortening Down to 1 ps. Crystals, 9, 167.
 Hill, A.H., Munger, E., Francis, A.T., Manifold, B. and Fu, D. (2019) Frequency Modulation Stimulated Raman Scattering Microscopy through Polarization Encoding. The Journal of Physical Chemistry B, 123, 8397-8404.
 Xiong, H.Q., Qian, N.X., Miao, Y.P., Zhao, Z.L., Chen, C. and Min, W. (2021) Super-Resolution Vibrational Microscopy by Stimulated Raman Excited Fluorescence. Light: Science & Applications, 10, Article No. 87.
 Li, Y.P., Shen, B.L., Zou, G.J., Wang, S.Q., Qu, J.L., Hu, R. and Liu, L.W. (2021) Fast Denoising and Lossless Spectrum Extraction in Stimulated Raman Scattering Microscopy. Journal of Biophotonics, 14, e202100080.
 Jiang, J., Grass, D., Zhou, Y., Warren, W.S. and Fischer, M.C. (2021) Beyond Intensity Modulation: New Approaches to Pump-Probe Microscopy. Optics Letters, 46, 1474.
 Lee, M., Herrington, C.S., Ravindra, M., Sepp, K., Davies, A., Hulme, A.N. and Brunton, V.G. (2021) Recent Advances in the Use of Stimulated Raman Scattering in Histopathology. The Analyst, 146, 789-802.
 Xiong, H.Q., Qian, N.X., Zhao, Z.L., Shi, L.Y., Miao, Y.P. and Min, W. (2020) Background-Free Imaging of Chemical Bonds by a Simple and Robust Frequency-Modulated Stimulated Raman Scattering Microscopy. Optics Express, 28, 15663-15677.
 Dong, P.-T., Zong, C., Dagher, Z., Hui, J., Li, J.J., Zhan, Y.W., Zhang, M., Mansour, M.K. and Cheng, J.-X. (2021) Polarization-Sensitive Stimulated Raman Scattering Imaging Resolves Amphotericin B Orientation in Candida Membrane. Science Advances, 7, eabd5230.
 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.
 Fonoberov, V.A. and Balandin, A.A. (2004) Interface and Confined Polar Optical Phonons in Spherical ZnO Quantum Dots with a Wurtzite Crystal Structure. Physica Status Solidi (c), 1, 2650-2653.
 Alim, K.A., Fonoberov, V.A. and Balandin, A.A. (2005) Origin of the Optical Phonon Frequency Shifts in ZnO Quantum Dots. Applied Physics Letters, 86, Article ID: 053103.
 Serrano, J., Manjon, E.J., Romero, A.H., Ivanov, A. Cardona, M., Lauck, R., Bosak, A. and Krisch, M. (2010) Phonon Dispersion Relations of Zinc Oxide: Inelastic Neutron Scattering and ab Initio Calculations. Physical Review B, 81, Article ID: 174304.
 Phan, T.L., Yu, S.C., Nghia, N.X. and Lam, V.D. (2010) Resonant Raman Scattering in ZnO Nanostructures Annealed at Different Temperatures. Journal of the Korean Physical Society, 57, 1569-1573.
 Mallet, E., Reveret, F., Disseix, P., Shubina, T.V. and Leymarie, J. (2014) Influence of Excitonic Oscillator Strengths on the Optical Properties of GaN and ZnO. Physical Review B, 90, Article ID: 045204.
 Milekhin, A.G., Yerykov, N.A., Sveshnikova, L.L., Duda, T.A., Zenkevitch, E.I., Kosolobov, S.S., Latyshev, A.V., Himenski, C., Surovtsev, N.V., Adichchev, S.V., Feng, Z.C., Wu, C.C., Wuu, D.S. and Zahn, D.R.T. (2011) Surface-Enhanced Raman Scattering of Light by ZnO Nanostructures. Journal of Experimental and Theoretical Physics, 113, 983-991.
 Zeng, X., Yan, S., Cui, J., Liu, H., Dong, J., Xia, W., Zhou, M. and Chen, H. (2015) Size- and Morphology-Dependent Optical Properties of ZnS: Al One-Dimensional Structures. Journal of Nanoparticle Research, 17, 188.
 Serrano, J., Cantarero, A., Cardona, M., Garro, N., Lauck, R., Tallman, R.F., Ritter, T.M. and Weinstein, B.A. (2004) Raman Scattering in β-ZnS. Physical Review B, 69, Article ID: 014301.