OSA's Digital Library

Optics Express

Optics Express

  • Editor: J. H. Eberly
  • Vol. 7, Iss. 9 — Oct. 23, 2000
  • pp: 299–304
« Show journal navigation

Analytical solutions for three and four diffraction orders interaction in Kerr media

N. Korneev  »View Author Affiliations


Optics Express, Vol. 7, Issue 9, pp. 299-304 (2000)
http://dx.doi.org/10.1364/OE.7.000299


View Full Text Article

Acrobat PDF (154 KB)





Browse Journals / Lookup Meetings

Browse by Journal and Year


   


Lookup Conference Papers

Close Browse Journals / Lookup Meetings

Article Tools

Share
Citations

Abstract

The analytical solution for the interaction of three diffraction orders in the Kerr medium is obtained by reducing the problem to the completely integrable Hamiltonian task. Intensities of all waves are periodic with propagation length and linearly related, the amplitudes are quasi-periodic and expressed in elliptic functions. Symmetrical four-order interaction also admits an analytical solution.

© Optical Society of America

1 Introduction.

The interaction of two plane waves in a Kerr medium (two wave mixing) is important in dynamic holography and nonlinear beam propagation studies. If the angle between beams is big, they do not produce higher diffraction orders, and their evolution is described using the Kogelnik diffraction theory [1

1. H. Kogelnik, “Coupled wave theory for thick hologram gratings,” Bell Syst.Tech.J 48,2909–2947 (1969).

]. If the angle is small enough, higher diffraction orders appear [2

2. R. Magnusson and T.K. Gaylord “Analysis of multiwave diffraction of thick gratings,” J.Opt.Soc.Am. 67,1165–1170(1977). [CrossRef]

]. Three orders must be taken into account in the first approximation for the interaction of a strong beam with a weak one, when the additional order is formed symmetrically with respect to the strong beam. At least four orders are needed when the two initial beams have equal intensities. Higher diffraction orders are easily observable, in particular, in photorefractive crystals [3

3. L. Solymar, D.J. Webb, and A. Grunnet-Jepsen, The Physics and applications of photorefractive crystals (Calderon, Oxford, 1996).

], and recently it was shown that the non-Bragg diffraction can produce in these materials big gains [4

4. A. Apolinar-Iribe, N. Korneev, and J.J. Sanchéz -Mondragon “Beam amplification resulting from non-Bragg wave mixing in photorefractive strontium barium niobate,” Opt.Lett. 23,1877–79 (1998). [CrossRef]

].

2 The Hamiltonian for interacting diffraction orders.

In a normalized form the nonlinear Schrödinger equation (NLSE) describing the propagation of wave in a Kerr medium is [7

7. V.E. Zakharov and A.B. Shabat, “Exact theory of two-dimensional self-focusing and one-dimensional self-modulation of waves in nonlinear media,” Sov.Physics JETP ,61,62–69(1972).

]:

izψ+12xxψ+κψ2ψ=0.
(1)

Here ψ is the amplitude of wave depending on one transversal coordinate x and the propagation coordinate z; κ is the nonlinear coefficient, which is positive for self-focusing media and negative for self-defocusing ones. We are interested in the initial value problem when ψ(x, z) is known for z=z 0, and we need to determine it for z>z 0. The N-order interaction corresponds to the amplitude written as:

ψ(x,z)=kSk(z)exp(ikx),
(2)

where Sk are the amplitudes of orders and the sum is taken over appropriate wavevectors k. For the odd number of orders we take for k values …-2K,-K, 0,K, 2K…. If the number of orders is even, they are … - 3/2K,-1/2K, 1/2K, 3/2K…. Here K=2π/Λ, where Λ is the interference fringes period. By substituting Eq.(2) into NLSE and combining the coefficients for Fourier harmonics, the infinite system of coupled nonlinear differential equations follows. We can, nevertheless, truncate this infinite system by assuming that higher diffraction orders are small and forcing the corresponding amplitudes to be zero. The coupled differential equations for the general case are written down in [5

5. N. Korneev, A. Apolinar-Iribe, and J.J. Sanchéz -Mondragon, “Theory of multiple beam interaction in photorefractive media,” J. Opt. Soc. Am. B 16,580–586(1999). [CrossRef]

]. For three orders we have :

izS1=12K2S1κ(S12S1+2S02S1+2S12S1+S02S1*),
(3)
izS0=κ(S02S0+2S12S0+2S12S0+2S0*S1S1),
(4)
izS1=12K2S1κ(S12S1+2S02S1+2S12S1+S02S1*).
(5)

It is well known that NLSE is formally the equation of motion for the Hamiltonian system with infinite number of degrees of freedom. The corresponding Hamiltonian function is [7

7. V.E. Zakharov and A.B. Shabat, “Exact theory of two-dimensional self-focusing and one-dimensional self-modulation of waves in nonlinear media,” Sov.Physics JETP ,61,62–69(1972).

] :

H(ψ,ψ*)=12(xψ2κψ4)dx.
(6)

The NLSE is given by the equation with Fréchet derivative i∂zψ=δH/δψ*.

We substitute Eq.(2) in Eq.(6), and obtain a Hamiltonian which depends on Sk and S*k. For soliton problems, the integration in Eq.(6) is done over the whole x axis. In our case we integrate over one period of interference fringes and normalize the result.

H(S,S*)=12kk2SkSk*12κa1+a2=a3+a4Sa1Sa2Sa3*Sa4*,
(7)

where every a can take any possible value of k. The Hamiltonian is real and symmetric with respect to S and its complexconj ugate.

The Hamiltonian equations for amplitudes are:

izSk=HSk*.
(8)

To obtain them in more familiar canonical form of mechanics, it is necessary to introduce real coordinates qk and momenta pk with Sk =1/√2(qk+ipk ). For the optical task, the complexno tation is more natural.

The truncated Hamiltonian Eq.(7) gives Eqs.(35), this is verified by direct substitution. It is also not difficult to show that for N orders, Eqs.(8,7) give the same result as Eqs.(2,1).

As a completely integrable task, NLSE has an infinite number of integrals of motion which are called Cn in the original work of Zakharov and Shabat [7

7. V.E. Zakharov and A.B. Shabat, “Exact theory of two-dimensional self-focusing and one-dimensional self-modulation of waves in nonlinear media,” Sov.Physics JETP ,61,62–69(1972).

]. Of them, apart from the Hamiltonian itself, we will utilize the intensity : I=∑k SkS* k , and the momentum: M=∑ kkSkS*k. Their conservation for our case does not follow automatically from the argument of [7

7. V.E. Zakharov and A.B. Shabat, “Exact theory of two-dimensional self-focusing and one-dimensional self-modulation of waves in nonlinear media,” Sov.Physics JETP ,61,62–69(1972).

]. It is necessary to prove that truncated sums are integrals of motion for the truncated Hamiltonian. For any function U (S, S*) the condition for its conservation is the equality to zero of the Poisson bracket :

izU={U,H}=kUSkHSk*USk*HSk
(9)

The conservation of Hamiltonian itself is evident from Eq.(9). The conservation of momentum and intensity can be proved for any number of orders by manipulations with Eq.(7), Eq.(9) and the quantity in question. It is also seen that I and M are in involution, i.e. {I,M}=0.

3 Solution for three orders.

If there are three integrals of motion for the Hamiltonian system with three degrees of freedom and they are in involution, it follows from the general theorem of mechanics that the trajectories of this system in the phase space lie on invariant 3-dimensional tori. The resulting movement is quasi-periodic, and the system is completely integrable [6

6. V.I. ArnoldMathematical Methods of Classical Mechanics (2-nd edition, Springer-Verlag1989).

],[8

8. D. ter HaarElements of Hamiltonian mechanics.(Pergamon Press, 1971).

].

There is a general procedure for obtaining angle-action variables from known integrals, but we just use conserved quantities to obtain the equations for order intensities. Let us denote the order intensities as I -1, I 0, I 1. Then I -1+I 0+I 1=I, and I 1-I -1=M/K are conserved. It follows: I1=12(I+MKI0), and I1=12(IMKI0), thus all intensities can be linearly expressed in terms of the zeroth order one. Using the equations for order intensities, the integral given by the Hamiltonian is expressed as:

H=(12K2κI0)(II0)12κI214κ((II0)2M2K2)2κRe(S0*2S1S1).
(10)

We rewrite Eq.(4) for the zeroth order intensity:

zI0=4κIm(S0*2S1S1).
(11)

Comparing Eq.(11) and Eq.(10), and taking into account that

(Im(S0*2S1S1))2+(Re(S0*2S1S1))2=I02I1I1=14I02((II0)2M2K2),
(12)

it is seen that zI0 is expressed as a square root of the fourth degree polynomial of I 0, and it follows that the intensity is an elliptic function. So, all order intensities are periodic with a propagation length. It is not true for amplitudes, because phases generally do not return to their initial value after a period of intensity. The solutions for amplitudes can be obtained by writing down the equations for i∂z ln(Sk )=i∂zSk/Sk . For example, for S 0 we have

izln(S0)=κ(S02+2S12+2S12+2S0*2S1S1S02).
(13)
Fig. 1. The trajectory in the complex plane for S 1. It is seen that the amplitude gains a constant phase shift after one period of intensity. The numerical solution of equations (3–5) for K=1.8, κ=0.4, and initial conditions S -1=0.3, S 0=1, S 1=0.1+0.24i.
Fig. 2. The maximal intensity transferred to the initially weak side order depending on the nonlinearity parameter q. The central beam intensity I 0=1.

In the right-hand side we have order intensities and the combination S0*2 S 1 S -1, and all of them can be expressed as functions of I 0 using Eqs.(1012). It follows that the phases are sums of periodic and linear functions of propagation length. The phase gains over the intensity period for 3 orders are not independent, because the combination S0*2 S 1 S -1 is periodic. In Figure 1 we present the computer calculation for a typical trajectory of the diffraction order. The formal solution in elliptic functions requiers the knowledge of the roots of the general fourth order polynomial, and thus it is quite cumbersome. Practically, the simplest way is the direct computer calculation. The most important property of three orders interaction is the quasi-periodical character of amplitudes and periodicity of intensities. So, it is possible to calculate amplitudes over one period of intensities, and to expand the solution for arbitrary propagation length.

4 Symmetrical four orders interaction.

To solve for four-beam interaction in the general case it is necessary to find the fourth integral of motion. The procedure of deriving it from the fourth NLSE integral given in the paper of Zakharov and Shabat [7

7. V.E. Zakharov and A.B. Shabat, “Exact theory of two-dimensional self-focusing and one-dimensional self-modulation of waves in nonlinear media,” Sov.Physics JETP ,61,62–69(1972).

] fails: because of the trunctaion we do not obtain the conserved quantity here. We can, nevertheless, solve a task if initial conditions are such that S -2=S 2 and S -1=S 1. These relations are maintained with propagation. Resulting evolution equations include only two orders :

izS1=18K2S1κ((I+S12+2Re(S1S2*))S1+(S12+4Re(S1S2*))S2),
(14)
izS2=98K2S2κ((I+S22)S2+(S12+4Re(S1S2*))S1).
(15)

The Hamiltonian as an integral of motion is given by

H=94K2S22+14K2S1212κ(I2+2S14+2S24+8S12Re(S1S2*)+16(Re(S1S2*))2)
(16)

Since the intensity I=2|S 1|2+2|S 2|2 is also the conserved quantity, there are two integrals of motion. This permits us to state that the task is completely integrable. The equation for the second order intensity is:

zI2=2κ(S12+4Re(S1S2*))Im(S1S2*),
(17)

and it is again possible with a help of H and I to express the right-hand side of this equation as a function of I 2 only, reducing the solution to integration. The integral, nevertheless, is more complicated than the elliptic one. Again, the intensities are periodic with propagation distance and amplitudes are quasi-periodic.

5 Representation of solutions on the Poincaré sphere and their properties.

The Poincaré sphere representation is useful to figure out how beam intensities change with propagation. Instead of the two amplitudes S 1, S 2 we introduce Stokes parameters given by: A=I 1-I 2=S 1 S*1-S 2 S*2, B=-i(S 1 S*2-S*1 S 2), C=S 1 S*2+S*1 S 2. The intensity conservation leads to the identity:

A2+B2+C2=I24,
(18)

thus the trajectories A(z),B(z), C(z) lie on a sphere. To determine their form one has to rewrite the Hamiltonian in terms of Stokes parameters. This gives the condition:

2K2A+κ(A2+IC+2AC+4C2)=const,
(19)

which means that projections of trajectories on A-C plane are concentric ellipses, and trajectories are intersections of elliptic cylinders parallel to B axis, and the sphere. In Fig.3, we show these ellipces built for illustration purposes with the direct numerical solution of motion equations. The ellipticity and orientation of axes of the ellipses do not depend on q parameter or the ellipse size. Their centers are given by A0=I6(14q) and C0=I6(q1), thus with changes in q, the center position moves along the straight line. The A parameter gives the intensity of beams (I 1,2=I/4±A/2), thus the Figure 3 contains a wealth of information about behavior of order intensities in particular situations. For example, trajectories with one maximum per period or with two of them separated by one minimum are possible, and the maximal values can be obtained by solving algebraic ecuations. The self-similar solutions for which which order intensities do not change with propagation are also well seen. For the situation of the Fig.3, there are sixo f them, and two are unstable.

Fig. 3. Phase portraits of trajectories with different initial conditions on the Poincaré sphere for a four-beam interaction. The nonlinearity parameter is q=0.5.

6 Discussion and conclusions.

The Hamiltonian method is a powerful tool for solving nonlinear propagation problems with a finite number of degrees of freedom. Once the problem in question is expressed in Hamiltonian form, the non-trivial integral of motion is obtained. If we are interested in the interaction of two nonlinear modes, and the intensity is conserved, we immediately have two integrals needed for a complete integrability with far-reaching consequences about the general behavior of solutions. An example is the symmetrical four-orders system discussed above. The three-order problem is solved by noting that the momentum is conserved as well.

For the interaction of bigger number of modes, additional integrals are needed. The direct substitution does not work for the fourth integral of Zakharov-Shabat [7

7. V.E. Zakharov and A.B. Shabat, “Exact theory of two-dimensional self-focusing and one-dimensional self-modulation of waves in nonlinear media,” Sov.Physics JETP ,61,62–69(1972).

]: a conserved quantity is not obtained. Even if the N-orders task is not completely integrable, truncated higher integrals of motion of Zakharov-Shabat can be used when orders beyond N are weak. It is exactly when truncated equations of motion are valid. Physically, the condition can be small nonlineartity (big q parameter), or limited propagation length when initially the energy was concentrated in a small number of central orders. It seems probable that the approximation to the analytical solution (action-angle variables) can still be obtained by using truncated integrals of motion, though some of these integrals are approximations.

In conclusion, we have shown that the propagation of three and symmetrical four diffraction orders in a Kerr medium can be reduced to a completely integrable Hamiltonian task. Integrals of motion can be obtained from known integrals of nonlinear Schödinger equation. In both cases that we have considered the order intesities are periodic functions of propagation length and the amplitudes are quasi-periodic.

Acknowledgments

The work was performed within the CONACyT project 211290-5-28495E.

References and links

1.

H. Kogelnik, “Coupled wave theory for thick hologram gratings,” Bell Syst.Tech.J 48,2909–2947 (1969).

2.

R. Magnusson and T.K. Gaylord “Analysis of multiwave diffraction of thick gratings,” J.Opt.Soc.Am. 67,1165–1170(1977). [CrossRef]

3.

L. Solymar, D.J. Webb, and A. Grunnet-Jepsen, The Physics and applications of photorefractive crystals (Calderon, Oxford, 1996).

4.

A. Apolinar-Iribe, N. Korneev, and J.J. Sanchéz -Mondragon “Beam amplification resulting from non-Bragg wave mixing in photorefractive strontium barium niobate,” Opt.Lett. 23,1877–79 (1998). [CrossRef]

5.

N. Korneev, A. Apolinar-Iribe, and J.J. Sanchéz -Mondragon, “Theory of multiple beam interaction in photorefractive media,” J. Opt. Soc. Am. B 16,580–586(1999). [CrossRef]

6.

V.I. ArnoldMathematical Methods of Classical Mechanics (2-nd edition, Springer-Verlag1989).

7.

V.E. Zakharov and A.B. Shabat, “Exact theory of two-dimensional self-focusing and one-dimensional self-modulation of waves in nonlinear media,” Sov.Physics JETP ,61,62–69(1972).

8.

D. ter HaarElements of Hamiltonian mechanics.(Pergamon Press, 1971).

9.

V.I. Bespalov and V.I. Talanov “Filamentary structure of light beams in nonlinear liquids,” JETP Lett. 3, 307–310(1966).

OCIS Codes
(090.7330) Holography : Volume gratings
(190.5330) Nonlinear optics : Photorefractive optics
(190.7070) Nonlinear optics : Two-wave mixing

ToC Category:
Research Papers

History
Original Manuscript: August 25, 2000
Published: October 23, 2000

Citation
N. Korneev, "Analytical solutions for three and four diffraction orders interaction in Kerr media," Opt. Express 7, 299-304 (2000)
http://www.opticsinfobase.org/oe/abstract.cfm?URI=oe-7-9-299


Sort:  Journal  |  Reset  

References

  1. H. Kogelnik, "Coupled wave theory for thick hologram gratings," Bell Syst. Tech. J 48, 2909-2947 (1969).
  2. R. Magnusson and T. K. Gaylord "Analysis of multiwave diffraction of thick gratings," J. Opt. Soc. Am. 67,1165 -1170 (1977). [CrossRef]
  3. L. Solymar, D. J. Webb and A. Grunnet-Jepsen, The Physics and applications of photorefractive crystals (Calderon, Oxford, 1996).
  4. A. Apolinar-Iribe, N. Korneev and J. J. Sanchez-Mondragon "Beam amplification resulting from non-Bragg wave mixing in photorefractive strontium barium niobate," Opt. Lett. 23, 1877-79 (1998). [CrossRef]
  5. N. Korneev, A. Apolinar-Iribe,and J. J. Sanchez -Mondragon, "Theory of multiple beam interaction in photorefractive media," J. Opt. Soc. Am. B 16, 80- 86(1999). [CrossRef]
  6. V. I. Arnold Mathematical Methods of Classical Mechanics ( 2-nd edition, Springer-Verlag 1989).
  7. V. E. Zakharov and A. B. Shabat, "Exact theory of two-dimensional self-focusing and one-dimensional self-modulation of waves in nonlinear media," Sov. Physics JETP, 61,62-69 (1972).
  8. D. ter Haar, Elements of Hamiltonian mechanics, (Pergamon Press, 1971).
  9. V. I. Bespalov and V. I. Talanov "Filamentary structure of light beams in nonlinear liquids," JETP Lett. 3, 307-310 (1966).

Cited By

Alert me when this paper is cited

OSA is able to provide readers links to articles that cite this paper by participating in CrossRef's Cited-By Linking service. CrossRef includes content from more than 3000 publishers and societies. In addition to listing OSA journal articles that cite this paper, citing articles from other participating publishers will also be listed.

Figures

Fig. 1. Fig. 2. Fig. 3.
 

« Previous Article  |  Next Article »

OSA is a member of CrossRef.

CrossCheck Deposited