Lecture 12
Lecture 12
Pi = ε 0 χij E j + χ ijk
( 2)
E j Ek + !
A 3×6 matrix χil is introduced, χ11 χ12 χ13 χ14 χ15 χ16
where l = 1…6 is a contracted index, χil = χ 21 χ 22 χ 23 χ 24 χ 25 χ 26
and i =1…3. χ
31 χ32 χ33 χ34 χ35 χ36
Three-wave mixing
Coupling between two optical waves ω1 and ω2:
{
E ω1 (t ) = Re E ω1 eiω1t = } 12 (E ω1
eiω1t + ( E ω1 )∗ e −iω1t = ) 12 (E ω1
eiω1t + c.c. )
E ω2 (t ) = Re{E ω2
}= 12 (E
e iω 2 t ω2
) = 12 (E
e iω 2 t + ( E ω 2 ) ∗ e − iω 2 t ω2
eiω2t + c.c. )
The total field:
E = E ω1 (t ) + E ω2 (t ) = (
1 ω1 iω1t
2
E e + E ω2 eiω2t + c.c. )
Linear part of the polarization PL:
(
PL = ε 0 χ(ω1 ) E ω1 (t ) + χ(ω2 ) E ω2 (t ) )
Nonlinear part of the polarization PNL:
PNL = χ ( 2) E E =
1 ( 2) ω1 ω1 2iω1t
4
χ E E e ( + E ω2 E ω2 e 2iω2t + 2 E ω1 E ω2 ei (ω1 +ω2 )t
The polarization PNL, when introduced into the Maxwell equations, becomes
the source of the radiation at frequencies 2ω1, 2ω2, ω1 + ω2 et ω1 − ω2
It causes an energy transfer between the fundamental and the mixed spectral
components
Three wave mixing: two initial components (ω1 and ω2) give raise to a third
one (ω3)
A phase matching condition has to be fulfilled : at most one efficient energy
transfer channel is in general possible
Second harmonic generation (SHG)
Coupling between ω and 2ω — other spectral components are omitted:
E = E ω (t ) + E 2ω (t ) = (
1 ω iω t
2
)
E e + E 2ω e i 2ωt + c.c.
Pi ω, NL =
2
(
1 ( 2 ) 2 ω ω ∗ iω t
χ ijk E j ( Ek ) e + c.c. )
ω
Pi ,2NL = χ ijk
4
(
1 ( 2 ) ω ω 2 iω t
E j Ek e + c.c. )
The wave equation in the time domain then reads:
∂2E ∂ 2 PNL
∇ E = µ 0ε 2 + µ 0
2
∂t ∂t 2
Absorption can be taken into account in ε; however, we neglect it here
The waves are supposed to propagate along z; their amplitudes do not
depend on x and y.
SHG: continued
E ωj (z , t ) =
1 ω
2
( )
E j ( z ) ei (ωt − k1z ) + c.c. E 2j ω (z , t ) = (
1 2ω
2
)
E j ( z ) e i (2ωt − k2 z ) + c.c.
The energy transfer between the two waves is assumed to be very small in
the scale of the wavelength:
dE ωj d 2 E ωj dE 2j ω d 2 E 2j ω
k1 >> 2
k 2 >>
dz dz dz dz 2
We finally obtain:
dE ωj iωη0 ( 2) 2ω ω ∗ −i (k2 −2 k1 )z
=− χ jkl Ek ( El ) e
dz 2nω
dE 2j ω iωη 0 ( 2) ω ω −i (2 k1 −k2 )z
=− χ Ek El e
dz 2n2ω jkl
Constant field approximation
The fundamental wave is supposed not to be depleted:
dE 2j ω iωη 0 ( 2) ω ω −i (2 k1 −k2 )z
=− χ jkl Ek El e
dz 2n2ω
Solution:
E 2j ω ( z ) = B − A ei∆kz
with
∆k = k 2 − 2k1
2) ω ω
ω η0 χ (jkl Ek El
A=
2 n2ω ∆k
I 2ω
n
= 2ω E 2j ω ( L)
2 1
= η30
( ) 2
ω2 χ (eff2) L2
I ω2
sin 2 ( 12 ∆kL)
2η0 2 n2ω nω2 (12 ∆kL )2
Character of the solution depends critically on the value of ∆k
∆k ≠ 0
Both waves do not propagate with the same phase velocity: they are not
constantly in phase, but become periodically out-of-phase. This leads to
a modulation of I2ω with the period (called coherence length):
2π 2π λ
lc = = =
∆k k 2 − 2k1 2(n2ω − nω )
I 2ω =
1
η30
( )
( 2) 2
ω χ eff
2
I ω2 L2
2 n2ω nω2
ne, 2ω (θ ) =
1
2
(no,ω + ne,ω (θ ))
no,ω
ne,ω
ne,2ω(θ)
The choice of the polarizations ne,2ω ne,ω(θ)
depends on the available coeffi-
cients of χijk (e.g. χ111 couples only
parallel polarizations and thus can
never allow the phase matching)
Three-wave mixing: summary
General equations of three-wave mixing
ω1 ± ω2 ± ω3 = 0 (frequency transformation)
k1 ± k 2 ± k3 = 0 (phase matching condition)
Sum and difference frequency generation (SFD, DFD):
•Input: two strong beams ω1 and ω2 ω1 ± ω2 = ω3
•Output: strong beam ω3 k1 ± k 2 = k3
Parametric generation (amplification of weak beams):
•Input: strong ω3 + weak ω1 ω3 − ω1 = ω2
•Output: medium ω2 + medium ω1 k3 − k1 = k 2
Up-conversion
•Input: strong ω1 + weak ω2 ω1 + ω2 = ω3
•Output: weak ω3 k1 + k 2 = k3
Four-wave mixing
Third order effect:
ω4 ( 3) ω1 ω2 ω3
PNL = χ ijkl E j Ek El
Wave equation:
(
2 2 2
ω n c − k A − 2ik
2
)
dA i (ωt −kz )
e = −3µ 0ω2 χ (3) A2 A∗ei (ωt −kz )
dz