^{1}

^{*}

^{2}

One of the significant problems of molecular spectroscopy is the determination and detailed analysis of how molecular vibrations are dephased. The dephasing of infrared-active (IR-active) vibrations of molecules was investigated by IR absorption spectroscopy. Pulse methods were used to investigate IR-vibrations as well. These methods revealed such coherent nonstationary effects as optical nutation, damping of the free polarization, photon echo, etc. New means of studying dephasing processes were uncovered by the method of nonstationary (time-domain) coherent anti-Stokes Raman scattering (CARS) spectroscopy. However, there are some aspects of CARS that still are not fully covered. One of them is related to Raman scattering by polaritons in dipole-active crystals whereas the second one is the increase of efficiency of CARS (minimization of the wave mismatch, the relationship between pulse width and the relaxation time, etc.). The purpose of the present research to study the case of “extreme” coherency between all interacting pulses (the duration of each pulse is smaller than characteristic times and those pulses are traveling with the same speed) in dipole-active crystals. In this research, we analyzed the process of simultaneous propagation of three waves (anti-Stokes, Stokes, and the pump) under CARS by polaritons. We have found some solutions modeling such simultaneous propagation. We also found the expression for the gain factor for such scattering. The gain factor was evaluated under the assumption of a given stationary pump field. It was shown that the typical values of the relative intensities were consistent with the experimental results.

Maker and Terhune were first who demonstrated the CARS technique [

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

E → l ( r → , t ) = e ^ l A l ( z , t ) exp [ i ( k l z z − ω l t ) ] + c . c . (1)

propagates along the z-axis. The subscripts a, l, s, and p denote the anti-Stokes, pump (laser), Stokes, and polariton wave fields, ω a , l , s , p are the frequencies, n a , l , s , p and k → a , l , s , p are the refractive indices and the wave vectors in the unpumped medium, and e ^ a , l , s , p are the real unit vectors of electromagnetic fields. The nonlinear medium is assumed to be nonmagnetic and transparent at the frequencies ω a , l , s . We use the anti-Stokes, Stokes, and polariton fields in the form

E → a ( r → , t ) = e ^ a A a ( z , t ) exp [ i ( k a z z − ω a t ) ] + c . c . , (2)

E → s ( r → , t ) = e ^ s A s ( z , t ) exp [ i ( k s z z − ω s t ) ] + c . c . , (3)

E → p ( r → , t ) = e ^ p A p ( z , t ) exp [ i ( W z z − ω p t ) ] + c . c . , (4)

where k a , s = q a , s n a , s ; q a , s = ω a , s / c ; W z = k l z − k s z ; ω p = ω l − ω s .

In the process of CARS, the nonlinear interaction of two electromagnetic waves ω l , s results in the generation of anti-Stokes and polariton waves. The system of shortened equations for the amplitudes A a , l , s , p is obtained from Maxwell’s equations by using the standard approximation of slowly-varying amplitudes [

∂ A a ∂ z + 1 v a z ∂ A a ∂ t = i 2 π ω a c n a cos ( θ a z ) { χ a A l A p e i Δ k z z + γ a ( | A l | 2 + | A s | 2 ) A a } , (5)

∂ A l ∂ z + 1 v l z ∂ A l ∂ t = i 2 π ω l c n l cos ( θ l z ) { χ l 1 A s A p + χ l 2 A a A p * e − i Δ k z z } , (6)

∂ A s ∂ z + 1 v s z ∂ A s ∂ t = i 2 π ω s c n s cos ( θ s z ) { χ s A l A p * + γ s ( | A l | 2 + | A s | 2 ) A s } , (7)

∂ A p * ∂ z + 1 v p z ∂ A p * ∂ t = i q p 2 ε p ∞ 2 W z ( W 2 q p 2 ε p ∞ − 1 ) A p * − i 2 π q p 2 W z { χ p 1 A l * A s + χ p 2 A a * A l e i Δ k z z } , (8)

where χ a , χ l 1 , l 2 , χ s , χ p 1 , p 2 , γ a , s are the corresponding tensor contractions of non-resonance quadratic and cubic nonlinear polarizabilities with unit vectors of polarization of interacting waves; e p ∞ is the non-resonance part of dielectric permeability at frequency ω p ; v a , l , s , p z are z-components of velocities of waves on ω a , l , s , p ; Δ k z ≡ k l z + W z − k a z is the wave mismatch between the pump, polariton, and anti-Stokes waves.

Given the strong polariton absorption we have [

| ∂ A p * ∂ z | ≈ | 1 v p z ∂ A p * ∂ t | ≪ q p 2 ε p ∞ 2 W z ( W 2 q p 2 ε p ∞ − 1 ) A p * , (9)

so that we can neglect in (8) the terms with the derivatives after which this equation yields

A p * = 4 π q p 2 ( W 2 − q p 2 ε p ∞ ) ( χ p 1 A l * A s + χ p 2 A a * A l e i Δ k z z ) . (10)

If we insert the obtained expression for the amplitude of polariton wave in (5)-(6), we get a system of 3 differential equations for A a , l , s as follows:

∂ A a ∂ z + 1 v a z ∂ A a ∂ t = i 2 π ω a c n a cos ( θ a z ) { 4 π q p 2 χ a χ p 1 * ( W 2 − q p 2 ε p ∞ ) A l 2 A s * e i Δ k z z + γ a 1 | A l | 2 A a + γ a | A s | 2 A a } , (11)

∂ A l ∂ z + 1 v l z ∂ A l ∂ t = i 2 π ω l c n l cos ( θ l z ) { 4 π q p 2 ( χ l 1 χ p 2 * + χ l 2 χ p 1 ) ( W 2 − q p 2 ε p ∞ ) A a A l * A s e − i Δ k z z + 4 π q p 2 χ l 2 χ p 2 ( W 2 − q p 2 ε p ∞ ) | A a | 2 A l + γ l 1 | A s | 2 A l + γ l | A l | 2 A l } , (12)

∂ A s ∂ z + 1 v s z ∂ A s ∂ t = i 2 π ω s c n s cos ( θ s z ) { 4 π q p 2 χ s χ p 2 ( W 2 − q p 2 ε p ∞ ) A l 2 A a * e i Δ k z z + γ s 1 | A l | 2 A s + γ s | A s | 2 A s } , (13)

where q p ≡ ω p / c , γ a 1 ≡ γ a + 4 π q p 2 χ a χ p 2 * ( W 2 − q p 2 ε p ∞ ) , γ s 1 ≡ γ s + 4 π q p 2 χ s χ p 1 ( W 2 − q p 2 ε p ∞ ) , and γ l 1 ≡ γ l + 4 π q p 2 χ l 1 χ p 1 * ( W 2 − q p 2 ε p ∞ ) .

The system (11)-(13) can be simplified if we use new variables

A ′ a ≡ A a e − i Δ k z z 2 (14)

and A ′ s ≡ A s e − i Δ k z z 2 . (15)

The system (11)-(13) in terms of A ′ a , s can be written as follows:

∂ A ′ a ∂ z + 1 v a z ∂ A ′ a ∂ t + i Δ k z 2 A ′ a = i 2 π ω a c n a cos ( θ a z ) { 4 π q p 2 χ a χ p 1 * ( W 2 − q p 2 ε p ∞ ) A l 2 A ′ s * + γ a 1 | A l | 2 A ′ a + γ a | A s | 2 A ′ a } , (16)

∂ A l ∂ z + 1 v l z ∂ A l ∂ t = i 2 π ω l c n l cos ( θ l z ) { 4 π q p 2 ( χ l 1 χ p 2 * + χ l 2 χ p 1 ) ( W 2 − q p 2 ε p ∞ ) A ′ a A ′ s A l * + γ l 1 | A ′ s | 2 A l + γ l | A l | 2 A l + 4 π q p 2 χ l 2 χ p 2 ( W 2 − q p 2 ε p ∞ ) | A ′ a | 2 A l } , (17)

∂ A ′ s ∂ z + 1 v s z ∂ A ′ s ∂ t + i Δ k z 2 A ′ s = i 2 π ω s c n s cos ( θ s z ) { 4 π q p 2 χ s χ p 2 ( W 2 − q p 2 ε p ∞ ) A l 2 A ′ a * + γ s 1 | A l | 2 A ′ s + γ s | A ′ s | 2 A ′ s } . (18)

And, finally, if we assume a “week” wave mismatch at Stokes and anti-Stokes frequencies, that is

| ∂ A ′ a , s ∂ z + 1 v a , s z ∂ A ′ a , s ∂ t | ≫ Δ k z 2 A ′ a , s , (19)

then the final system of equations simulating CARS can be expressed as

∂ A ′ a ∂ z + 1 v a z ∂ A ′ a ∂ t = i 2 π ω a c n a cos ( θ a z ) { 4 π q p 2 χ a χ p 1 * ( W 2 − q p 2 ε p ∞ ) A l 2 A ′ s * + γ a 1 | A l | 2 A ′ a + γ a | A ′ s | 2 A ′ a } , (20)

∂ A l ∂ z + 1 v l z ∂ A l ∂ t = i 2 π ω l c n l cos ( θ l z ) { 4 π q p 2 ( χ l 1 χ p 2 * + χ l 2 χ p 1 ) ( W 2 − q p 2 ε p ∞ ) A ′ a A ′ s A l * + 4 π q p 2 χ l 2 χ p 2 ( W 2 − q p 2 ε p ∞ ) | A ′ a | 2 A l + γ l 1 | A ′ s | 2 A l + γ l | A l | 2 A l } , (21)

∂ A ′ s ∂ z + 1 v s z ∂ A ′ s ∂ t = i 2 π ω s c n s cos ( θ s z ) { 4 π q p 2 χ s χ p 2 ( W 2 − q p 2 ε p ∞ ) A l 2 A ′ a * + γ s 1 | A l | 2 A ′ s + γ s | A ′ s | 2 A ′ s } . (22)

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 z 0 / A 0 ( A 0 and τ 0 are the peak amplitude and characteristic duration of the pump, z 0 = c τ 0 ). After that, the system (20)-(22) can be reduced to

∂ A ˜ ′ a ∂ z ˜ + 1 v ˜ a z ∂ A ˜ ′ a ∂ t ˜ = i { C a 1 A ˜ l 2 A ˜ ′ s * + C a 2 | A ˜ l | 2 A ˜ ′ a + C a 3 | A ˜ ′ s | 2 A ˜ ′ a } , (23)

∂ A ˜ l ∂ z ˜ + 1 v ˜ l z ∂ A ˜ l ∂ t ˜ = i { C l 1 A ˜ ′ a A ˜ ′ s A ˜ l * + C l 2 | A ˜ ′ s | 2 A ˜ l + C l 3 | A ˜ l | 2 A ˜ l + C l 4 | A ˜ ′ a | 2 A ˜ l } , (24)

∂ A ˜ ′ s ∂ z ˜ + 1 v ˜ s z ∂ A ˜ ′ s ∂ t ˜ = i { C s 1 A ˜ l 2 A ˜ ′ a * + C s 2 | A ˜ l | 2 A ˜ ′ s + C s 3 | A ˜ ′ s | 2 A ˜ ′ s } , (25)

where A ˜ ′ a , s ≡ A ′ a , s A 0 , A ˜ l ≡ A l A 0 , t ˜ ≡ t τ 0 , C a 1 ≡ 2 π ω a z 0 c n a cos ( θ a z ) 4 π q p 2 χ a χ p 1 * A 0 2 ( W 2 − q p 2 ε p ∞ ) ,

C a 2 ≡ 2 π ω a z 0 c n a cos ( θ a z ) γ a 1 A 0 2 , C a 3 ≡ 2 π ω 0 z 0 c n a cos ( θ a z ) γ a A 0 2 , C l 2 ≡ 2 π ω l z 0 c n l cos ( θ l z ) γ l 1 A 0 2 ,

C l 3 ≡ 2 π ω l z 0 c n l cos ( θ l z ) γ l A 0 2 , C l 1 ≡ 2 π ω l z 0 c n l cos ( θ l z ) 4 π q p 2 ( χ l 1 χ p 2 * + χ l 2 χ p 1 ) A 0 2 ( W 2 − q p 2 ε p ∞ ) ,

C l 4 ≡ 2 π ω l z 0 c n l cos ( θ l z ) 4 π q p 2 χ l 2 χ p 2 A 0 2 ( W 2 − q p 2 ε p ∞ ) , C s 1 ≡ 2 π ω s z 0 c n s cos ( θ s z ) 4 π q p 2 χ s χ p 2 A 0 2 ( W 2 − q p 2 ε p ∞ ) ,

C s 2 ≡ 2 π ω s z 0 c n s cos ( θ s z ) γ s 1 A 0 2 , C s 3 ≡ 2 π ω s z 0 c n s cos ( θ s z ) γ s A 0 2 .

We are looking for stationary solutions as

A ˜ ′ a , s ( z ˜ , t ˜ ) ≡ B a , s ( ξ ˜ ) e i Φ a , s ( ξ ˜ ) and A ˜ l ( z ˜ , t ˜ ) ≡ B l ( ξ ˜ ) e i Φ l ( ξ ˜ ) , (26)

where ξ ˜ ≡ t ˜ − z ˜ / ν ˜ z ; ν ˜ z is the velocity of simultaneously propagating waves at the frequencies ω a , l , s ; B a , l , s and Φ a , l , s 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):

d B a d ξ ˜ = − κ a C a 1 B l 2 B s sin ( Φ ) , (27)

d Φ a d ξ ˜ = κ a { C a 1 B l 2 B s B a cos ( Φ ) + C a 2 B l 2 + C a 3 B s 2 } , (28)

d B l d ξ ˜ = κ l C l 1 B a B s B l sin ( Φ ) , (29)

d Φ l d ξ ˜ = κ l { C l 1 B a B s cos ( Φ ) + C l 2 B s 2 + C l 3 B l 2 + C l 4 B a 2 } , (30)

d B s d ξ ˜ = − κ s C s 1 B l 2 B a sin ( Φ ) , (31)

d Φ s d ξ ˜ = κ s { C s 1 B l 2 B a B s cos ( Φ ) + C s 2 B l 2 + C s 3 B s 2 } , (32)

where κ a . s . l ≡ v a , s , l z v z / ( v z − v a , s , l z ) , Φ ≡ 2 Φ l − Φ s − Φ a .

If we introduce the amplitude of simultaneously propagated waves as

Q ≡ − B a 2 κ a C a 1 = − B s 2 κ s C s 1 = B l 2 κ l C l 1 (33)

we could reduce the system above to

d Q d ξ ˜ = α Q 2 sin ( Φ ) , (34)

d Φ d ξ ˜ = 2 α Q cos ( Φ ) + β Q , (35)

where

λ a 2 = − κ a C a 1 , λ l 2 = κ l C l 1 , λ s 2 = − κ s C s 1 , α ≡ 2 λ a λ s λ l 2 , β ≡ 2 κ l ( C l 2 λ s 2 + C l 3 λ l 2 + C l 4 λ a 2 ) − κ s ( C s 2 λ l 2 + C s 3 λ s 2 ) − κ a ( C a 2 λ l 2 + C a 3 λ s 2 ) . (36)

The system (34)-(35) can be further simplified as follows

d Q d x = Q 2 sin ( Φ ) , (37)

d Φ d x = Q ( 2 cos ( Φ ) + β ˜ ) , (38)

where x ≡ ξ ˜ α , β ˜ = β α .

We can reduce the number of equations by using the integral of motion

Q = 1 2 cos ( Φ ) + β ˜ , (39)

where Q > 0 , β ˜ > 2 .

If we express the phase Φ as the function of Q in (37) we get

∫ d Q ( 4 − β ˜ 2 ) Q 4 + 2 β ˜ Q 2 − 1 = 1 2 x (40)

The integral on the left can be found as follows:

∫ d Q ( 4 − β ˜ 2 ) Q 4 + 2 β ˜ Q 2 − 1 = − ( i ( 1 − ( β ˜ − 2 ) Q 2 ) ( 1 − ( β ˜ + 2 ) Q 2 ) × F ( i sinh − 1 ( ( 2 − β ˜ ) Q | β ˜ − 2 β ˜ + 2 ) ) ) / ( 2 − β ˜ ) ( ( 4 − β ˜ 2 ) Q 4 + 2 β ˜ Q 2 − 1 ) , (41)

where F ( x ˜ | m ) is the elliptic integral of the first kind with the parameter m = k 2 , sinh − 1 ( x ˜ ) is the inverse hyperbolic sine function.

In ^{−}^{3} cm/MW [^{2} - 10^{3} MW and z = 1 cm , then g ≈ 1 .

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

∂ A ˜ s ∂ z ˜ = i C s 1 A ˜ a * e i Δ k z z (42)

∂ A ˜ a * ∂ z ˜ = − i C a 1 A ˜ s e − i Δ k z z (43)

After introducing new variables

F s ≡ A ˜ s e − i Δ k z z / 2 and F a * ≡ A ˜ a * e i Δ k z z / 2 (44)

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 F s ( z ) ) as

∂ 2 F s ∂ z 2 + ( Δ k z 2 ) 2 F s = C a 1 C s 1 F s . (45)

We solve this equation by adopting a trial solution for F s in the form

F s ( z ) = F s ( 0 ) e g z (46)

where g represents a gain factor. Then we substitute (46) into (45) to obtain the approximate value for g as

g ≈ ( C a 1 C s 1 ) 1 2 ≈ C (47)

(here we assumed that the pump was strong enough to provide C a 1 C s 1 ≫ ( Δ k z 2 ) 2 ).

Finally, the expression for g can be reduced to

g ≈ C ≈ 8 π 2 ω z 0 χ 2 A 0 2 / ( c n ) (48)

( g ≈ ( C a 1 C s 1 ) 1 2 = ( 2 π ω a z 0 c n a cos ( θ a z ) 4 π q p 2 χ a χ p 1 * A 0 2 ( W 2 − q p 2 ε p ∞ ) 2 π ω s z 0 c n s cos ( θ s z ) 4 π q p 2 χ s χ p 2 A 0 2 ( W 2 − q p 2 ε p ∞ ) ) 1 / 2 ≈ ( 8 π 2 ω z 0 / ( ε ′ n c ) χ 2 A 0 2 ) ≈ 8 π 2 ω z 0 χ 2 A 0 2 / ( c n ) ).

As the experimental data for this gain, we used the following [^{+} laser ; 30 ps, the peak output power ≈ 2.5 kW, the wavelength was 514.5 nm, the cross-section ≈ 10^{−18} cm^{−2}, γ f ≈ 10 cm − 1 , and χ ≈ 10 − 8 esu . In [^{−1}. Both the experimental results for the gain factor in [

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.

The authors declare no conflicts of interest regarding the publication of this paper.

Feshchenko, V. and Feshchenko, G. (2020) Nonstationary CARS by Polaritons. Journal of Applied Mathematics and Physics, 8, 1949-1958. https://doi.org/10.4236/jamp.2020.89146