Основы физической и квантовой оптики
..pdf
|
|
|
|
|
|
|
|
|
|
|
|
|
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), |
||
w×( 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 |
|
||
откуда получим: |
|
|
n−2 ( ω ) − n−2 ( 2ω ) |
|
||||||||
sin2 θs |
= |
|
||||||||||
|
|
0 |
0 |
|
(8.34). |
|||||||
n |
− 2 |
( 2ω ) − n− 2 ( 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=e– el=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. Пространственные оптические солитоны в среде с