Добавил:
Опубликованный материал нарушает ваши авторские права? Сообщите нам.
Вуз: Предмет: Файл:

Основы физической и квантовой оптики

..pdf
Скачиваний:
6
Добавлен:
05.02.2023
Размер:
3.95 Mб
Скачать

 

 

 

 

 

 

 

 

 

 

 

 

 

182

Pω1 2 ( t ) =

1

Pω3 1 2

exp[ i( ω + ω

t )] + c.c. =

 

i

 

 

2

 

i

 

 

 

 

1

2

 

= [d

ijk

E ω1

E ω2

exp[ i( ω + ω

2

)t ] + c.c]

 

 

 

j

k

 

 

 

1

 

 

 

 

или:

 

 

 

 

 

 

 

 

 

 

 

 

Pω3 1 + ω2

= 2d

ijk

E ω1 E

ω2

 

 

 

 

 

i

 

 

 

 

 

j

k

 

 

 

 

 

Отметим, что коэффициенты dijk ,

как правило, определяются из результатов

экспериментов. В большинстве случаев это эксперименты по генерации второй гармоники, когда ω1 = ω2 = ω . В этом случае:

P2ω = d

ijk

E ω E

ω .

 

 

 

 

 

 

 

 

 

 

 

 

 

i

 

j

k

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Рассмотрим теперь простейший одномерный случай, считая что

световое поле распространяется вдоль оси Z, а

x

=

y

= 0 . Полагаем, что

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

имеется три бегущих вдоль Z плоских световых волны с частотами ω1 , ω2 ,

ω3

с полями:

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

E( ω1 )( z,t ) =

1

[E ( z ) × exp[ i( w t - k

z )] + c.c.]

 

 

 

 

 

 

 

 

 

 

 

i

 

2

 

1i

 

1

1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Ek( ω2 )( z,t ) =

1

[E2k ( z ) ×exp[ i( w2t - k2 z )] + c.c.]

 

 

 

 

 

 

 

 

 

 

 

 

2

 

 

 

[E

 

 

 

 

z )] + c.c.]

 

 

 

 

 

E( ω3 )( z,t ) =

1

( z ) × exp[ i( w t - k

 

 

(8.16),

 

 

 

 

3

 

 

 

 

j

 

 

 

 

2

 

3 j

 

3

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

где индексы i, j, k означают декартовы координаты, и каждый из индексов может принимать значения x или y.

Подставим соотношения (8.13) в нелинейное волновое уравнение (8.9). Очевидно, что нужно провести дифференцирование полей по координате z,

поскольку они однородны вдоль

поперечных

координат (

x

=

y

= 0 ).

 

 

 

 

 

 

 

 

 

 

 

 

 

Тогда, взяв i

ю компоненту, получим:

 

 

 

 

 

 

 

Ñ2 E( ω1 )( z,t ) =

2

E( ω1 )( z,t ) =

1

2

{ E ( z ) × exp[ i( w t - k

z )] + c.c.} .

 

 

 

 

 

 

 

 

 

i

 

z2

i

2

 

z2

1i

1

1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Проведя дифференцирование, учтем, что амплитуды взаимодействующих волн слабо изменяются на расстояниях порядка длины волны. Соответственно, можно пренебречь второй производной поля по z:

2 E

E

 

 

 

 

 

 

1i << k

1i . В итоге получим:

 

 

 

 

z2

1 z

 

 

 

E1i ( z ) ] × exp[ i( w t - k z )] + c.c. (8.17).

Ñ2 E( ω1 )( z,t ) = -

1

[ k 2 E

( z ) + 2ik

 

 

 

i

2

1 1i

 

1

z

1

1

 

 

 

 

Аналогичные соотношения получим для других составляющих светового поля. Учтем, что нелинейное волновое уравнение должно удовлетворяться

183

для каждой частотной компоненты, запишем соответствующее соотношение для поля E1( ω1 )( z,t ) :

 

k 2

 

 

E

 

 

 

 

 

 

 

[

1

E + ik

 

1i ] × exp[ i( w t - k

z )] + c.c. =

 

 

 

 

 

 

 

2

 

1i

 

1

z

1

1

 

 

 

 

(8.18).

 

 

1

 

 

 

 

 

 

 

 

2

 

= [

w2m

eE

× exp[ i( w t - k

z )] + c.c - m

 

[ P

( z,t )]

 

 

t 2

2

1 0

 

1i

1

1

 

 

0

NL

i

Подставим в это уравнение выражение для нелинейной поляризации, учтем,

что w2m

e = k 2

, и, проведя

 

несложные

преобразования,

получим три

1

0

1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

уравнения:

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

E1i

 

 

 

 

 

 

 

d`

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

= -iw

 

m0

 

E

 

 

E*

 

× exp[ -i( k

 

- k

 

- k

 

)z ]

 

(8.19),

 

 

 

z

 

 

1

 

 

e1

ijk

 

 

3 j

 

2k

 

 

 

3

 

 

 

2

 

 

 

1

 

 

 

 

 

 

E*

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

= iw

 

 

m

0

 

d`

E E*

×exp[ -i( k

 

- k

 

+ k

 

 

 

 

 

 

 

 

 

2k

 

 

 

 

 

 

 

 

 

)z ]

 

(8.20)

 

 

 

z

 

 

e2

 

 

 

 

 

 

 

 

 

2

 

 

 

kij

 

 

1i

 

3 j

 

 

 

1

 

 

 

3

 

 

 

2

 

 

 

 

 

 

 

 

E3 j

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

= -iw

 

 

 

 

m

 

 

`

 

 

 

 

E ×exp[-i( k

 

+ k

 

- k )z ]

 

 

 

 

 

 

 

 

 

 

 

 

0 d

 

 

E

 

 

 

(8.21).

 

 

 

z

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

3

 

 

 

 

e3

 

 

jik

 

1i

 

2k

 

 

 

 

 

1

 

 

 

2

3

 

Полученные уравнения могут использоваться для анализа различных взаимодействий в среде с квадратичной оптической нелинейностью. Остановимся лишь на случае генерации второй гармоники. При этом ω1 = ω2 ,

а ω3 = 2ω1 . Соответственно, нам необходимы лишь два из полученных

уравнений – первое или второе и последнее. При анализе считаем взаимодействие слабым, чему соответствует условие dE1i dz = 0 , т.е.

пренебрегаем изменением амплитуды волны основной частоты в процессе генерации второй гармоники. Поэтому мы можем рассматривать лишь последнее из этих трех уравнений - уравнение (8.21):

E3 j

 

 

 

 

 

 

 

= -iw

m

0 d`

E

E ×exp[ i × Dk × z ]

(8.22),

z

 

 

e

jik

1i

1k

 

где w = w1 = ω3 2 ; ε = ε3 ; Dk = k3j - k1i - k1k . Здесь k1i - волновое число для волны с частотой ω1 при поляризации света вдоль направления i.

Положим, что i, j, k могут принимать значения x и y вдоль направления поляризации собственных мод в кристаллической среде. Световая волна с частотой ω1 представляет собой, в общем случае, суперпозицию собственных

мод, поэтому может иметь и x, и y составляющие. Величина волнового числа k зависит от состояния поляризации, что отмечено верхними индексами.

Легко получить решение уравнения (8.22) для случая отсутствия второй гармоники на входе нелинейной среды, при E3 j ( 0 ) = 0 . Для кристалла

длиной L это решение имеет вид:

184

E

( L ) = -iw

m0

d`

E

E

×

exp[ i × Dk × L ] - 1

 

(8.23).

e

iDk

3 j

 

jik

1i

1k

 

 

Интенсивность поля второй гармоники на выходе среды найдем в виде:

 

 

 

4m0

 

 

 

 

 

sin2 [

1

Dk × L ]

 

 

 

 

 

 

 

 

 

 

E*

3 j ( L ) × E

( L ) = -w2

( d`

)2 E 2

E 2

L2

×

 

2

 

(8.24).

 

 

 

 

 

3 j

 

e

jik

1i

1k

 

[

1

Dk × L ] 2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Выражение для выходной мощности найдем, используя соотношение:

P( 2ω )

 

 

 

 

 

 

 

1

 

 

e

*

 

 

=

 

 

 

 

 

E3 j E 3 j

(8.25),

 

2

 

 

m0

S

 

 

 

 

 

где S – площадь поперечного сечения пучка. Подставляя в это выражение

(8.25), получим:

P( 2ω )

 

 

 

 

 

 

 

 

 

sin2 [

1

Dk × L ]

 

= 2

 

m0

 

× w2 ( d`

)2 E 2

E 2

L2

×

2

(8.26).

 

 

 

 

 

S

 

 

e

 

jik

1i

1k

 

[

1

Dk × L ] 2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

2

 

 

 

 

Запишем это выражение через мощность первой гармоники:

P( 2ω )

m

 

3 2

 

w2

( d`jik )2

L2

 

P( ω )P( ω )

 

sin2 [

1

Dk × L ]

 

0

 

 

 

 

 

 

= 8

 

 

×

 

 

 

×

i k

×

2

 

(8.27).

 

 

 

 

 

 

 

 

 

S

 

e

 

 

 

n3

 

 

S 2

 

[

1

Dk × L ] 2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

2

 

 

 

 

Или, при P( ω ) = P( ω ) = P( ω ) / 2 оно примет вид:

 

 

 

 

 

i

 

k

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

( 2ω )

 

 

 

 

3

 

2 `

2 2

 

 

( ω )

 

sin

2

[

1

Dk × L ]

 

 

P

 

m

0

 

2

 

w ( d jik

) L

 

P

 

 

 

 

 

 

2

 

 

 

 

 

= 2

 

 

 

×

 

 

×

 

 

×

 

 

 

 

 

(8.28).

 

P( ω )

e

 

n3

 

 

 

 

1

 

 

 

 

 

 

 

 

 

 

 

S

[

Dk × L ] 2

 

 

 

 

 

 

 

 

2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

8.4. Условия фазового синхронизма при генерации второй гармоники

Из (8.28) вытекает, что эффективная генерация второй гармоники

требует

выполнения условия

фазового

синхронизма Dk=0.

Поскольку

w3 = 2 w,

а

w1 = w2 = w, то

требование

Dk=0 эквивалентно

условию

k( 2ω ) = 2k ( ω ) .

Действительно,

причиной

возникновения излучения с

удвоенной частотой является нелинейная поляризация, наводимая в среде световой волной с частотой ω . Волна наведенной поляризации распространяется в среде со скоростью распространения индуцирующего

излучения. Если k ( 2ω ) ¹ 2k ( ω ) , то между излучением с удвоенной частотой и наводимой поляризацией возникает фазовый сдвиг, нарастающий с

185

расстоянием, так что в некоторой точке он может достигать величины p. Расстояние, на котором фазовый сдвиг изменяется от нуля до p, называется "когерентной длиной":

lc

=

2π

=

 

2π

(8.29).

 

 

ω ) - 2k ( ω )

 

 

Dk k( 2

 

Ее величина служит оценкой максимальной длины кристалла, которую можно эффективно использовать для генерации второй гармоники. В области нормальной дисперсии показатель преломления растет с увеличением частоты, поэтому:

Dk = k ( 2ω ) - 2k( ω ) =

2ω

( n( 2ω ) - n( ω ) )

(8.30).

 

 

c

 

Здесь использовано соотношение k=wn/c. Тогда когерентная длина может быть выражена в форме:

 

=

 

πc

 

=

 

λ

lc

 

 

 

 

 

(8.31),

( n( 2

ω )

 

 

 

 

 

- n( ω ) ) 2( n( 2

ω ) - n( ω ) )

где l – длина волны падающего света.

Для обеспечения условий фазового синхронизма используется оптическая анизотропия кристаллов. Из (8.30) видим, что для этого

необходимо выполнить условие n( 2ω ) = n( ω ) , т.е. условие равенства показателей преломления среды на основной частоте и на частоте второй гармоники. В области нормальной дисперсии показатели преломления как обыкновенной, так и необыкновенной волн, распространяющихся в заданном направлении, увеличиваются с ростом частоты. Таким образом, выполнить условие равенства показателей преломления невозможно, если волны с частотами w и 2w относятся к одному типу (обыкновенные или необыкновенные). Однако условие фазового синхронизма может быть выполнено при использовании волн разных типов.

В качестве примера рассмотрим зависимости величин показателей преломления необыкновенной и обыкновенной волн в одноосном,

отрицательном кристалле (пример такого

кристалла –

кристалл ниобата

 

 

 

лития),

от

угла

q между

направлением

 

 

 

распространения и оптической осью (осью Z). Эти

 

 

 

зависимости схематично изображены на рис. 8.3.

 

 

 

Можно видеть, что для обыкновенной волны

 

 

 

основной частоты и необыкновенной волны с

 

 

 

 

 

 

удвоенной

частотой

имеются

направления

 

 

 

распространения

света,

удовлетворяющие

 

 

 

 

Рис. 8.3.

условию фазового синхронизма (одно из них

 

отмечено стрелкой).

 

 

 

 

 

 

 

 

 

186

Величина необыкновенного показателя преломления одноосного кристалла в зависимости от угла в плоскости, содержащей оптическую ось, определяется соотношением:

 

1

 

=

cos2 θ

+

sin2 θ

(8.32),

n2

( θ )

n2

n2

 

 

 

e

 

 

0

 

e0

 

где θ – угол между направлением распространения света и оптической осью;

n0 и ne0 – главные значения

показателей преломления кристалла.

Таким

образом, направления синхронизма могут быть найдены из условия:

 

1

 

=

cos2 θs

+

sin2 θs

 

(8.33),

 

 

 

 

 

 

n2 ( ω ) n2

( 2ω ) n2 ( 2ω )

 

 

o

 

 

 

 

0

 

 

e0

 

откуда получим:

 

 

n2 ( ω ) n2 ( 2ω )

 

sin2 θs

=

 

 

 

0

0

 

(8.34).

n

2

( 2ω ) n2 ( 2ω )

 

 

 

 

 

 

 

 

 

 

e

0

 

 

8.5. Распространение светового пучка в нелинейно оптической среде

При анализе нелинейной поляризации, наводимой световым полем в среде с кубичной нелинейностью, мы выяснили, что возникает составляющая поляризации на частоте падающей световой волны, величина которой пропорциональна квадрату ампитуды светового поля. Это значит, что световое поле с неоднородным распределением интенсивности наводит в такой среде оптические неоднородности, которые, в свою очередь, могут изменять пространственную структуру светового поля. В этом случае говорят о пространственном самовоздействии светового поля или пучка. Рассмотрим основные особенности такого самовоздействия и некоторые его результаты.

Пусть световой пучок падает на границу нелинейно –

оптической

среды. Полагаем его поле в виде:

 

E( x, z,t ) = A( x, z ) × exp( ikz - iwt )

(8.35),

где x и z – поперечная и продольная координаты. В данном случае предполагается, что размер пучка в направлении x значительно превышает длину волны света, а в направлении y его поле однородно. Соответственно, можно использовать параксиальное приближение, т.е. считать, что амплитуда A(x,z) изменяется в направлении z намного медленнее, чем в направлении x. Уравнение, описывающее поведение пучка в нелинейной диэлектрической среде, получим, исходя из стандартной системы уравнений Максвелла:

rotH = D

t

rotE = − B

t

187

Для этого используем обычную процедуру, применяя к обеим частям второго из этих уравнений операцию rot:

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

D

 

 

 

 

 

 

 

rot( rotE ) = -

 

 

rotB = -m

rotH = -m

,

 

 

 

 

 

 

t

 

t 2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

t

 

 

 

 

 

 

 

 

 

 

= -m 2

 

 

 

 

 

 

 

 

 

 

 

) - Ñ2

 

 

 

D

. Поскольку для диэлектрической среды

или grad( divE

E

 

 

 

 

 

 

 

 

 

 

t 2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

) = 0 , то получаем:

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

grad( divE

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

= m 2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Ñ2

 

 

D

= -mew2

 

 

 

 

 

 

 

 

(8.36).

 

 

 

 

E

E

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

t 2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

С учетом (8.35) выражения для пространственных

производных поля

E

 

имеют вид:

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

2 E

=

2 A

× exp( ikz - iwt )

 

x2

x2

 

 

 

 

2 E =

2 A exp( ikz - iwt ) + 2ik

A exp( ikz - iwt ) - k 2 A × exp( ikz

z 2

z2

 

 

 

z

В параксиальном приближении можно положить

второй производной 2 A , тогда (8.36) примет вид:

z2

2 A + 2ik A - k 2 A = -mew2 A x2 z

2 A << A z2 z

(8.37a),

-iwt ) (8.37б).

ипренебречь

(8.38).

Здесь мы ограничились скалярным приближением, считая световой пучок линейно поляризованным. Поскольку рассматриваемая среда является оптически нелинейной, представим выражение для ее диэлектрической проницаемости e в форме:

e = el

+ enl (

 

E

 

2 )

(8.39),

 

 

где el и enl – линейная и

нелинейная

составляющие диэлектрической

проницаемости. Но mel w2 = k 2 , поэтому:

2 A + 2ik

A

 

 

 

 

 

 

 

 

 

 

 

 

 

 

x2

 

 

z

A i

2 A

 

imenl w2

 

 

 

 

 

z -

 

×

x2

 

=

 

 

 

 

 

 

A

 

 

2k

 

 

 

2k

 

 

 

Правую часть этого уравнения запишем в форме:

 

 

ime

nl

w2e

l

=

ik 2e

nl

 

=

ik

 

×

e

nl

 

 

2kel

 

 

 

2kel

 

2

 

el

 

 

 

 

 

 

 

 

 

 

 

= -menl w2 A или:

(8.40).

(8.41).

Но enl можно представить как enl=eel=De. Т.к. e=n2, то: De = ( n2 )¢ = 2n × Dn и

enl

=

2n × Dnnl

=

2Dnnl

. Таким образом, (8.40) примет вид:

el

n2

 

 

 

n

 

 

 

 

 

188

 

 

 

A

-

i

×

2 A

=

ikDn

nl

A

(8.42).

z

2k

x2

n

 

 

 

 

 

 

 

Это уравнение для амплитуды светового пучка, распространяющегося в нелинейно – оптической среде. В правой части данного уравнения величина нелинейной добавки к показателю преломления среды nnl является

функцией пространственных координат и зависит от локальной интенсивности светового поля. Рассмотрим различные механизмы оптической нелинейности, приводящие к изменениям показателя преломления среды на частоте распространяющегося в ней светового поля.

1. Керровская оптическая нелинейность, обусловленная нелинейной связью наведенной поляризации в среде и напряженности электрического поля в световой волне. Это изменение показателя преломления на частоте падающей волны в среде с кубичной нелинейностью (см. 8.7). Величина изменения показателя преломления определяется соотношением:

Dn( к ) = n

 

× I( x, z )

(8.43),

nl

( 2 )

 

 

где n(2) – нелинейный показатель

 

преломления

среды. Для обычных

материалов (газы, жидкости, стекло, кристаллы) величина n(2) очень мала, так что заметные изменения показателя преломления, приводящие к эффектам пространственного самовоздействия световых пучков, наблюдаются при интенсивностях света в сотни МВт/см2 и более. Данный механизм характеризуется очень высокой скоростью отклика (она ограничивается лишь инерционностью электронов).

Отметим, что знак нелинейного показателя преломления n(2) определяет знак нелинейности среды. В случае n(2)>0 среда называется самофокусирующей, при n(2)<0 ее называют самодефокусирующей.

2. Термооптическая нелинейность – изменение показателя преломления при изменении температуры среды, обусловленном поглощением света. Величина изменения показателя преломления в этом случае определяется соотношением:

 

 

Dn( TO ) =

n

× dT [ I( x, z )]

(8.44),

 

 

 

 

nl

T

 

 

 

 

 

 

где

n

- температурный коэффициент показателя преломления, а dT[I(x,z)] –

 

 

T

 

 

 

локальное изменение температуры среды вследствие поглощения света. Термооптическая нелинейность может быть значительно более сильной, чем керровская, однако она и значительно более медленная.

3. Оптическая нелинейность, обусловленная фоторефрактивным эффектом. Как отмечалось в п. 7.3, такая нелинейность обусловлена изменением показателя преломления вследствие линейного электрооптического эффекта при условии индуцирования светом в фоторефрактивном кристалле электрического поля пространственного

189

заряда. Она может быть поистине гигантской, приводя к значительным возмущениям показателя преломления даже при микроваттных интенсивностях света. Однако она также значительно медленнее в сравнении с керровской и во многих случаях медленнее термооптической. Изменения показателя преломления для разных механизмов фоторефракции определяются выражениями:

Для фотовольтаического механизма:

 

Dn( phv ) = A × I 1

×

 

I( x, z ) Id

 

(8.45),

 

 

 

nl

d

1

+ I( x, z ) Id

 

 

 

 

где А – константа,

зависящая

от физических

характеристик и

электрооптических коэффициентов материала; Id – так называемая темновая интенсивность, определяемая условием σd = σ ph I , т.е. интенсивность света,

при которой величина фотопроводимости материала равна его темновой проводимости.

Для дрейфового механизма:

Dn( drift ) = B ×

 

Ibg

× E

 

(8.46),

 

 

 

ext

nl

Ibg

+ I( x, z )

 

 

 

 

 

 

где Ibg – интенсивность фоновой подсветки; Eext

напряженность внешнего

электрического поля, приложенного к кристаллу; B - постоянная, зависящая от физических характеристик и электрооптических коэффициентов материала.

Для диффузионного механизма:

 

 

 

 

 

 

 

 

Dndif

= -

1

n3r × E

 

( z ) = -

1

n3r ×

kBT

 

1

 

×

dn( x )

 

(8.47).

 

 

 

e n( x )

dx

nl

 

2 l

D

 

2 l

 

 

Отметим, что в случае диффузионного механизма транспорта электрических зарядов изменение показателя преломления носит нелокальный характер, т.е. оно определяется величиной градиента концентрации фотовозбужденных носителей заряда. Соответственно, эффекты самовоздействия световых пучков проявляются в данном случае не в самофокусировке или самодефокусировке пучков, а в их самоискривлении.

Уравнение (8.42) позволяет проследить, как изменяется профиль светового пучка в нелинейной среде при различных механизмах и величине оптической нелинейности, а также различных параметрах светового пучка. Его решение может быть получено в аналитической форме лишь в ряде частных случаев. В большинстве же ситуаций возможно только его численное решение.

Особый интерес представляют решения уравнения, отвечающие режимам сохранения поперечного профиля светового пучка, называемым режимами пространственных оптических солитонов. В таком случае решение (8.42) ищем в виде:

190

A( x,z ) = U( x )×exp( igz )

Тогда:

 

A( x,z ) = ig × A( x,z ) = ig ×U ( x )×exp( igz ),

 

z

 

 

 

 

 

 

 

 

 

 

 

 

 

2 A( x,z )

=

2U × exp( igz ) .

 

 

 

 

 

 

x2

 

x2

 

 

 

Подставляя эти соотношения в (8.42), получим:

 

 

i

 

 

 

 

2

 

 

 

 

ikDnнл

 

 

igU × exp( igz ) -

 

×

 

U × exp( igz )

=

×U ( x ) × exp( igz ),

 

 

 

 

 

2k

 

 

x2

 

 

 

n

 

 

 

 

 

 

 

 

что после упрощения дает:

 

 

 

 

 

 

 

 

 

 

 

 

U ¢¢ - 2kgU +

2k 2 Dnнл

×U = 0 или

 

 

 

 

 

 

 

 

 

 

 

 

n

 

 

 

 

 

2k 2 Dnнл

 

 

 

 

 

 

U ¢¢ +

 

 

 

 

 

-

2kg ×U = 0

 

 

 

 

 

 

n

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

где U ¢¢ =

2U .

 

 

 

 

 

 

 

 

 

 

 

 

x2

 

 

 

 

 

 

 

 

 

 

 

(8.48).

(8.49),

Это уравнение для пространственных солитонов, которое называют также нелинейным уравнением Шредингера. Действительно, в данной форме записи аналогия с соответствующим уравнением очевидна. В то же время это уравнение совпадает по форме с волновым уравнением для градиентного диэлектрического волновода. Соответственно, его решения представляют собой собственные моды оптического волновода, индуцированного в нелинейной среде самим световым пучком. Следует отметить, что нелинейная добавка к показателю преломления среды Dnnl может быть как положительной, так и отрицательной. Соответственно, нелинейная добавка g к постоянной распространения положительна либо отрицательна. В первом случае (Dnnl>0) среда является самофокусирующей, а решение уравнения (8.49) представляет собой светлый солитон, т.е. световой пучок, распространяющийся без дифракционного расплывания в нелинейно – оптической среде. Противоположный случай соответствует самодефокусирующей среде. Собственное решение (8.49) при этом называют темным солитоном, под которым понимается как бы бездифракционное поведение неосвещенной области в световом поле.

Рассмотрим отдельно решения нелинейного уравнения Шредингера (или солитонного уравнения) для различных механизмов оптической нелинейности среды.

8.5.1. Пространственные оптические солитоны в среде с

191

керровской нелинейностью

В данном случае нелинейное возмущение показателя преломления среды определяется нелинейной восприимчивостью третьего порядка и, как уже отмечалось, может быть записано в форме:

Dnnl = n( 2 ) × I( x, z ) = n( 2 ) ×U 2 .

Солитонное уравнение (8.49) примет вид:

 

2k

2 n

 

 

U ¢¢ +

 

( 2 )

×U 2

- 2kg

×U = 0

 

 

 

 

n

 

 

 

 

 

 

Далее, обозначая kn2 = a получим окончательно: n

U -

1

×U ¢¢ =

α

×U 3

2kg

g

 

 

 

Коэффициент a может быть как положительным, первом случае, для самофокусирующей среды, солитона имеет вид:

(8.50).

(8.51).

так и отрицательным. В решение для светлого

 

2g 1 / 2

 

U ( x ) =

 

 

sec h( x )

(8.52),

 

 

a

 

 

где x = ( 2kg )12 x – нормированная поперечная координата.

В случае самодефокусирующей среды решение дает профиль темного солитона в виде:

g 1 / 2

 

 

 

 

 

U ( x ) =

 

 

tanh( x / 2 )

(8.53),

 

a

 

 

 

 

где x = ( -2kg )12 x .

Таким образом, в случае среды с керровской нелинейностью профили световых пучков, соответствующих светлым и темным пространственным солитонам, выражаются в аналитической форме и определяются соотношениями (8.52) и (8.53). Рисунок 5.4 иллюстрирует поперечные профили интенсивности таких солитонов.