Голографические фотонные структуры в наноструктурированных материалах
..pdf
|
∑ M j (r,t)e |
−iK |
j |
r |
+ M120 (r,t)e− |
iK |
20 |
r |
+ M121 |
(r,t)e−iK |
21r + |
|
|
|||
|
N |
|
|
|
|
|
|
|
|
|
|
|
||||
M = 0.5 j =0 |
|
|
|
|
|
|
|
|
|
|
|
|
, |
|
|
|
|
+ M102 (r,t)e−iK 02 r + M112 (r,t)e−iK12 r + c.c. |
|
|
|
|
|||||||||||
|
|
|
. |
(3.11) |
||||||||||||
|
|
N |
|
|
|
|
|
|
|
|
|
|
|
|
||
|
|
−iK |
r |
+ n120 (r,t)e−iK 20 r + n121 (r,t)e−iK |
|
|
||||||||||
|
|
∑ n j (r,t)e |
|
|
j |
|
21r + |
|
|
|||||||
n = nst + 0.5 j =0 |
|
|
|
|
|
|
|
|
|
|
|
. |
|
|
||
|
|
+ n102 (r,t)e−iK 02 r + n112 (r,t)e−iK12 r + c.c. |
|
|
|
|||||||||||
|
|
|
|
|
||||||||||||
где Mj, nj – амплитуды |
гармоник |
|
концентрации |
мономера и показателя |
преломления |
|||||||||||
основной решетки, M20, M21, M02, M12, |
n20, n21, n02, n12 – |
амплитуды первых гармоник |
||||||||||||||
концентрации мономера и показателя преломления ДДР. |
|
|
|
|
||||||||||||
Принимая во внимание малость амплитуд ДДР, можно пренебречь их влиянием на гармоники основной решетки и разделить рассматриваемую задачу на две части. Это позволяет в первой части найти решение для пространственных гармоник основной решетки и, используя его, решить задачу формирования ДДР.
ДДР формируются в результате интерференции падающих записывающих волн с волнами, дифрагированными на высших пространственных гармониках основной решетки. Принимая во внимание малость дифрагированных волн, можно не учитывать изменения падающих записывающих волн и ограничить рассмотрение дифрагированными волнами на рассматриваемой высшей пространственной гармонике El, E′l
∂El = iG (E n (τ)e−i K 0l y + E n (τ)e−i K0 l y + E nl 0 (τ, y) + E nl1 (τ, y)) |
|
|
||||||||||||
|
∂y |
l |
0 l |
|
|
1 l −1 |
|
|
|
0 1 |
1 1 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||
|
∂E′ |
|
(E0 nl (τ)e |
|
|
|
|
|
|
|
|
(τ, y)) |
, |
(3.12) |
|
l |
= iGl |
−i K 0′ l y |
+ E1nl −1 (τ)e |
−i K |
0′ l y |
|
0l |
1l |
|
|
|||
|
∂y |
|
|
|
+ E0 n1 |
(τ, y) + E1n1 |
|
|
||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
где Gl=π/(λcosθl), λ – длина волны света в ФПМ, nl(τ) – |
l– тая гармоника n, определенная в |
|||||||||||||
решении для основной решетки с учетом αd≈0 (l>1), n1j(τ,у) – |
первая гармоника j– |
той ДДР |
||||||||||||
показателя преломления n, индекс j=l0, |
l1, 0l, 1l – |
соответствует рассматриваемой ДДР с |
||||||||||||
волновым вектором Kj, формируемой одной из записывающих волн и волной, дифрагированной на l– той пространственной гармонике ОГДР, K0l=| K0l| – модуль вектора фазовой расстройки (см. рис.3.8), τ=t/Tm – относительное время, Tm=1/(K12Dm) – время диффузии, Dm – коэффициент диффузии, K1=|K1|.
Учитывая малость амплитуд ДДР, можно пренебречь их взаимовлиянием. Используя методику из подраздела 3.1.3 [64,65], где кинетические уравнения записи дополняются дифракционным уравнением в приближении заданного поля и без учета поглощения,
запишем интегро– |
дифференциальные кинетические уравнения записи для ДДР с вектором |
||||||||||||||||||||||||||||||
Kl1: |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
∂ |
|
l 0 |
|
2 |
l 0 |
|
|
|
2k |
|
|
|
E E* (τ, y) |
|
|
l 0 |
|
|
|
|
|
|||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
0 |
l |
|
|
|
|
|
|
|
|
|
|||||||
|
|
|
M1 |
(τ, y) = −rl 0 M1 |
(τ, y) − |
|
|
|
2kM 0 (τ) |
|
|
|
|
|
+ M1 |
(τ, y) |
|
|
|
||||||||||||
|
∂τ |
b |
|
E |
2 |
+ E 2 |
|
, |
(3.13) |
||||||||||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
0 |
|
1 |
|
|
|
|
|
|
|||||
|
|
∂ |
|
|
|
|
|
2k |
M 0 (τ) |
|
E0 El* (τ, y) |
|
|
|
M1l 0 (τ, y) |
|
|
|
M1l 0 (τ, y) |
|
|
||||||||||
|
|
nl 0 (τ, y) = δn |
|
2k |
|
|
+ |
− δn r 2 |
|
|
|||||||||||||||||||||
|
|
∂τ |
p b |
|
|
|
|
|
|
||||||||||||||||||||||
|
1 |
|
|
|
|
M |
n |
|
|
E2 |
+ E2 |
|
|
|
|
M |
n |
|
|
|
i l 0 |
M |
n |
|
|
||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
0 |
1 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|||
где El(τ,y) – |
определяется решением системы (3.12) в приближении заданного поля и учетом |
||||||||||||||||||||||||||||||
|n1l1|<|nl– 1 |; |
М0(τ) |
– |
нулевая гармоника M, |
определенная в решении для основной решетки |
|||||||||||||||||||||||||||
выражениями |
с |
учетом αd≈0; |
M1l0(τ,y), |
– первая |
гармоника |
дополнительной |
решетки |
||||||||||||||||||||||||
концентрации мономера с вектором Kl0, rl0=Kl0/K1, Kl0=|Kl0|.
Для решения (3.13) используем интегральное преобразование Лапласа по пространственной координате у
70
|
¶ |
|
|
|
|
|
|
|
|
|
|
|
|
|
2 |
k |
|
|
|
|
|
|
|
|
|
|
|
|
||
|
|
M |
1l 0 (t, p) = -rl20 M1l 0 |
(t, p) - |
|
M1l 0 (t, p) + |
|
|
|
|||||||||||||||||||||
|
|
|
|
|
|
|
|
|||||||||||||||||||||||
|
¶t |
|
|
|
|
|
|
|
|
|
|
|
|
b |
|
|
|
|
|
|
|
|
|
|
|
|
||||
|
|
|
|
|
|
|
2k |
|
|
|
|
|
k |
|
|
|
|
n1l 0 (t, p) |
|
|
|
|||||||||
|
|
|
|
|
+ 2i |
|
|
|
|
M 0 (t) |
|
|
|
|
|
|
|
Gl |
|
|
|
|
+ H (t, p) |
|
||||||
|
|
|
|
|
|
|
1+ m0 |
|
p |
|
||||||||||||||||||||
|
|
|
|
|
|
|
b |
|
|
|
|
|
|
|
|
|
|
|||||||||||||
|
¶ l 0 |
|
|
|
|
|
|
l 0 |
(t, p) |
|
|
|
|
2 |
k |
|
l 0 |
(r, t) |
, |
(3.14) |
||||||||||
|
|
2 M1 |
|
|
|
|
|
|
|
M1 |
|
|
||||||||||||||||||
|
|
n1 |
|
(t, p) = -dni rl 0 |
|
|
|
|
|
|
+ dnp |
|
|
|
|
|
- |
|
||||||||||||
¶t |
|
|
|
M n |
|
|
|
b |
|
M n |
|
|
||||||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
2k |
M 0 (t) |
|
|
|
|
k |
|
|
|
|
n1l 0 (t, p) |
|
|
||||||||||
|
|
|
|
|
- 2i |
|
|
|
|
|
|
|
|
|
|
|
Gl |
|
|
|
|
+ H (t, p) |
|
|||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||||||||||
|
|
|
|
|
|
|
b M n |
1+ m0 |
|
|
|
p |
|
|
|
|||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
d − y |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
m0 nl −1 (t)}e−i K |
||||
|
|
|
|
|
|
|
|
|
H(τ,p)=L{H(τ,y)}= |
|
|
|
|
|
|
||||||||
где введены |
обозначения |
|
= L M 0 (t) ∫{nl (t) + |
|
0 l y′dy¢ , |
||||||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
0 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
2k |
|
2k |
2iGl kdnp |
|
2 |
2 |
|
|
|
|
|
|
|
|
|
|
||||||
F1 = rl20 + |
|
, F2 |
= - |
|
|
|
|
|
|
|
, m0 = E1 |
E0 |
. |
|
|
|
|
|
|
|
|
|
|
|
b |
|
1+ m0 |
|
|
|
|
|
|
|
|
|
|||||||||||
|
b |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||||
Учитывая, что функция n1l0(t,у) является медленно меняющейся по сравнению с |
|||||||||||||||||||||||
M0×exp(– |
F1×τ) [64,65], и используя начальное условие М1l0(τ=0)=0 и теорему о среднем для |
||||||||||||||||||||||
|
|
|
|
|
|
|
b |
|
|
|
|
|
|
b |
|
|
|
|
|
|
|
|
|
определенного |
интеграла ∫ j(x) f (x)dx = j(x)∫ f (x)dx , |
где |
a ≤ ξ ≤ b , |
получим |
следующее |
||||||||||||||||||
|
|
|
|
|
|
|
a |
|
|
|
|
|
|
a |
|
|
|
|
|
|
|
|
|
решение для первого уравнения из (3.14): |
|
|
|
|
|
|
|
|
|
|
|
||||||||||||
|
|
|
|
|
|
|
F eF1τ |
|
nl 0 |
(t, p) |
τ |
|
|
|
F eF1τ τ |
|
|
|
|
|
|||
|
M1l 0 (t, p) = - |
2 |
1 |
|
|
∫ M 0 |
(t¢)e− F1τ′ |
dt¢ - |
|
2 |
∫ H (t¢, p)e− F1τ′ dt¢. |
|
(3.15) |
||||||||||
|
|
dn p |
|
|
|
p |
dnp |
|
|||||||||||||||
|
|
|
|
|
|
|
|
|
|
0 |
|
|
|
0 |
|
|
|
|
|
||||
Подставляя полученное решение (3.15) во второе уравнение из (3.14) и используя обратное интегральное преобразование Лапласа по пространственной координате у, решение для n1l0(t,y) запишем в виде:
|
τ |
d − y |
|
|
nl 0 |
|
∫Q(t¢, y¢) × H 0 |
|
(3.16) |
(t, y) = F2 ∫ Q(t¢, y) + |
(t¢, t, y - y¢) dt¢ , |
|||
|
|
0 |
|
|
|
0 |
|
|
|
|
|
|
2k |
τ |
H (t¢, y) |
e− F1τ′ d |
|
|
H (t, y) |
|
|
|
|
|
|
где Q( , y) |
C r 2 |
|
|
eF1τ |
|
|
¢ |
|
|
, C |
= |
n / |
n |
, |
||
- |
|
∫ |
- |
|
||||||||||||
t |
= |
n l 0 |
b |
|
|
t |
M n |
n |
|
d i |
d p |
|
||||
|
|
|
|
|
0 |
M n |
|
|
|
|
|
|
|
|
||
|
|
|
F2 |
τ |
|
|
|
|
|
|
|
|
|
|
|
|
τ |
|
|
||||
H 0 |
(t¢, t, y) = 1+ |
|
∫ |
R(t¢¢)dt¢¢×J1 2 |
F2 |
(d - y) ∫ |
R(t¢¢)dt¢¢ |
, |
|||
|
|||||||||||
|
|
|
d - y τ′ |
|
|
|
τ′ |
|
|
|
|
R(t) = C
n
2 |
|
2k |
|
F1τ |
τ |
M |
0 |
(t¢) |
|
− F1τ′ |
|
- |
|
|
×e |
∫ |
|
|
|
||||
rl 0 |
b |
|
|
M n |
e |
|
|||||
|
|
|
|
|
0 |
|
|
||||
dt¢ - M 0 (t) ,
M n
H (τ, y) = M 0 (τ)d − y (nl (τ) + 
m0 nl−1 (τ))e−i K0l y′dy′, J1(x) – функция Бесселя.
∫
0
Таким образом, выражение (3.16) определяет временную динамику пространственного распределения амплитуды ДДР с вектором решетки Kl0, образованной записывающими волнами и дифрагированной на l– той пространственной гармонике ОГДР. Из решения видно, что ДДР в отличии от гармоник основной решетки имеет пространственно неоднородные распределения амплитуды и фазы вдоль глубины решетки даже в отсутствии поглощения.
Проведя аналогичную процедуру решения для ДДР с вектором решетки Kl1 можно получить следующие решения:
71
l1 |
|
τ |
(t, y) = F2 |
|
|
n1 |
∫ Q(t¢, y) |
|
|
|
|
|
|
0 |
|
F |
m |
τ |
|
где H 0 (t¢, t, y) = 1+ |
2 |
0 |
∫ R(t¢¢)dt¢¢× |
|
d - y |
||||
|
τ′ |
|||
d − y |
|
|
+ ∫Q(t¢, y¢) × H 0 |
|
(3.17) |
(t¢, t, y - y¢) dt¢ , |
||
0 |
|
|
|
|
|
|
|
|
|
|
|
|
|
τ |
|
|||
m0 |
||||||
J1 2 F2 |
(d - y) ∫ |
R(t¢¢)dt¢¢ , H(τ,y) полностью |
||||
|
|
τ′ |
|
|
||
совпадает с приведенным в (3.16), а Q(τ,y) и R(τ,y) совпадают с приведенными в (3.16) с точностью до замены rl0 на rl1=Kl1/K1.
Проведя аналогичную процедуру решения для ДДР с вектором решетки Kl1 можно получить следующие решения:
l1 |
|
τ |
y |
|
|
(t, y) = F2 |
|
|
|
(3.18) |
|
n1 |
∫ Q(t¢, y) + ∫Q(t¢, y¢) × H0 |
(t¢, t, y - y¢) dt¢ , |
|||
|
|
|
0 |
|
|
|
|
0 |
|
|
|
|
F |
m |
τ |
|
|
|
τ |
|
|
|
|
||||||||
где H 0 (t¢, t, y) = 1+ |
2 |
0 |
∫ |
R(t¢¢)dt¢¢×J1 2 F2 |
m0 |
y ∫ |
R(t¢¢)dt¢¢ , H(τ,y) полностью совпадает |
||
y |
|
||||||||
|
|
τ′ |
|
|
τ′ |
|
|
||
с приведенным в (3.16), а Q(τ,y) и R(τ,y) совпадают с приведенными в (3.16) с точностью до замены rl0 на rl1=Kl1/K1.
Решения для ДДР с векторами K0l совпадают с решениями (3.16) для ДДР с векторами
K0l с точностью до замены rl0 на r0l=K0l/K1 и H (τ, y) = M 0 (τ)d∫− y(nl −1 (τ) + |
|
nl (τ))e−i K0′ l y′dy′ . |
m0 |
||
0 |
|
|
Аналогично, решения для ДДР с векторами K1l совпадают с решениями (3.17) для |
||
ДДР с векторами K1l с точностью до замены rl1 |
на r1l=K1l/K1 и |
|
d − y( ) −i K ′ y′ ′ .
H (τ, y) = M 0 (τ) ∫ nl −1 (τ) +
m0 nl (τ) e 0 l dy
0
На рис.3.9 представлены результаты расчета модуля и аргумента пространственного профиля решетки с волновым вектором K20 при Cn=2, dnp=0.01, q0=800, q1=850 (в воздухе), q2=740, d=85мкм, DK02=71.4.
Из рис.3.9(а,б) видно, что решетка n120(t,y) имеет неравномерный квазипериодический амплитудно– фазовый пространственный профиль. Наличие фазовой составляющей профиля ДДР приводит к образованию дополнительного максимума угловой селективности в угловом положении соответствующем порождающей гармонике ОГДР (рис. 3.9б).
ед. |
0,020 |
|
|
b=0.25 |
|
|
|
|
|
|
|||
|
|
|
b=5 |
|
||
. |
0,015 |
|
|
|
|
|
, отн |
|
|
|
|
|
|
0,010 |
|
|
|
|
|
|
p |
|
|
|
|
|
|
n |
|
|
|
|
|
|
(τ,y)/δ |
0,005 |
|
|
|
|
|
20 |
|
|
|
|
|
|
1 |
|
|
|
|
|
|
n |
|
|
|
|
|
|
|
0,000 |
|
|
|
|
|
|
0,0 |
0,2 |
0,4 |
0,6 |
0,8 |
1,0 |
|
|
y/d, отн. ед. |
|
|
||
а)
. |
1,5 |
|
|
|
b=0.25 |
|
), рад |
|
|
|
|
|
|
|
|
|
|
b=5 |
|
|
|
|
|
|
|
|
|
p |
|
|
|
|
|
|
δn |
|
|
|
|
|
|
,y)/ |
0,0 |
|
|
|
|
|
(τ |
|
|
|
|
|
|
20 |
|
|
|
|
|
|
1 |
|
|
|
|
|
|
arg( n |
-1,5 |
|
|
|
|
|
|
|
|
|
|
|
|
|
0,0 |
0,2 |
0,4 |
0,6 |
0,8 |
1,0 |
|
|
y/d, отн. ед. |
|
|||
б)
Рисунок 3.9
Амплитуда ДДР может быть сравнима с амплитудой порождающей гармоники при
72
малых величинах угла записи (θ0, θ1) и толщины материала (d). С увеличением d и θ0, θ1 период осцилляции профиля ДДР вдоль у уменьшается, как и амплитуда ДДР. Необходимо отметить, что неоднородность фазового фронта приводит к пространственной неоднородности направления вектора решетки.
Расчеты также показали, что для решетки с вектором K12 амплитудно– фазовый профиль n112(τ,y) совпадает с приведенными на рис. 3.9 (а,б).
Для иллюстрации на рис.3.10 приведем двумерные контурные графики пространственных профилей, как гармоник основной решетки, так и первых гармоник ДДР образованных на второй пространственной гармонике.
|n20(x,y)| |
|n21(x,y)| |
x
x
|
K20 |
|
y |
|
K21 |
|
|
|
|n2(x,y)| |
K2 |
y |
K1 |
|n1(x,y)| |
x |
x |
y
y
Рисунок 3.10
Необходимо отметить основные особенности пространственных профилей ДДР. Во– первых, профили ДДР имеют двумерное квази периодическое амплитудное распределение, а, во– вторых, неоднородность фазовой составляющей профиля ДДР приводит к повороту вектора решетки и искажению фазового фронта.
3.2 Дифракционные свойства пространственно неоднородных ОГДР
Рассмотрим процесс считывания ДР произвольно поляризованным монохроматическим световым пучком в пренебрежении остаточным поглощением ФПМ
|
|
|
|
|
Er (t, r) = ∑ ∫ e0i Eri (k0 ) × exp[i × (w×t - k0 × r)]dkτ0 + ê.ñ. , |
(3.19) |
|||
|
|
|
|
|
i = s , p |
|
|
|
|
где Ei |
(k |
0 |
) |
– |
угловой спектр (УС), k |
= k∆θ |
0 |
– тангенсальная компонента волнового вектора |
|
r |
|
|
|
τ0 |
|
|
|
||
|
|
|
|
|
|
|
|
73 |
|
k0, связанная с углом отклонения ∆q0 плосковолновой компоненты углового спектра от оси
пучка k0¢.
Пространственная и векторная геометрия процесса считывания представлена на рис.3.11 для двух дифракционных порядков.
x
|
|
|
ep |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
es |
E0 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
N |
|
|
|
|
|
|
|
|
|
|
|
θ0 γ |
θ2 |
|
y |
|||
|
|
|
|
|
|
|||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
k2′ |
θ1 |
|
|
||
|
|
|
|
|
k0′ |
|
|
|
|
|
|
|
|
E2(d) |
|
k1′ |
|
|
|
||
|
|
|
|
|
|
|
|
|
||
|
|
|
E1(d) |
|
d |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
а) |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
x |
|
|
|
|
|
|
|
|
|
|
|
|
γ |
|
|
|
|
|
|
|
|
|
|
|
θ1 |
|
|
|
|
|
|
|
|
|
|
θ2 |
|
|
θ0 |
|
y |
|
|
|
|
θ2 |
|
|
|
|||
|
|
|
|
|
|
|
|
|
||
|
|
|
k2 |
|
θ1 |
θ0 |
|
|
||
K02 |
|
K2 |
k1 |
|
k0′ |
|
|
|||
|
P2 |
|
|
|
||||||
|
|
|
k2′ |
k0 |
|
|
||||
|
|
|
P1 |
|
k1′ |
|
|
|||
|
K2 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
K1 |
|
|
|
|
|
|
|
|
N2 |
|
|
K1 |
|
|
|
|
|
|
|
|
N1 |
|
|
P0 |
|
|
|||
|
|
|
|
|
|
|
||||
|
|
|
|
|
|
|
|
|
N0 |
б) |
|
|
|
|
|
|
|
|
|
|
|
Рисунок 3.11
Световое поле E в области ФПМ в силу дифракции считывающего пучка Er (3.19) на пространственных гармониках решетки можно записать в виде суммы N+1 световых пучков
Ej
N |
|
1 |
|
N |
|
|
|
E (t, r) = ∑ E j |
= |
|
∑∑ ∫ eij E ij (r, k j ) × exp[i ×(w×t - k j × r )]dkτj |
+ ê.ñ. , |
(3.20) |
||
|
|||||||
j=0 |
2 |
i=s , p j=0 |
|
|
|||
каждый из которых представлен двумя составляющими вектора напряженности Ej в соответствующем ортогональном поляризационном базисе, заданном двумя ортами ejp и ejs, лежащими в плоскости, перпендикулярной оси пучка Ej. Здесь Eji(r, kj) – медленно меняющиеся амплитуды плосковолновых составляющих угловых спектров компонент Ej×eji, j=0 соответствует проходящему пучку, j=1..N – дифрагированному пучку на решетке с Kj=
j·K1, причем ejp лежит в плоскости дифракции XOY, а ej |
s – перпендикулярен ей (см. |
74 |
|
рис.3.11).
Напряженность электрического поля E(t,r) в области взаимодействия описывается векторным волновым уравнением, следующим из уравнений Максвелла
|
|
|
|
¶2 |
|
|
, |
|
(3.21) |
|
rot rot E (t, r ) = -m0 e |
0 ¶t 2 [e(t, r) × E (t, r)] |
|
||||||
|
|
|
|
||||||
где возмущение диэлектрической проницаемости ε представляется в виде |
|
|
|||||||
|
|
|
|
N |
|
|
|
||
e(t, r) = e0 |
+ De(t, r) = e0 |
+ 0.5nst |
n0 (t, r) + ∑n j (t, r)eiK j ×r + ê.ñ. , |
(3.22) |
|||||
|
|
|
|
j =1 |
|
|
|
||
где ε 0 = nst2 , nj(y) – |
определены |
решениями, |
полученными ранее |
при |
рассмотрении |
||||
процесса записи ОГДР.
Подставляя выражения (3.20) и (3.22) в (3.21) и следуя методу медленно– меняющихся амплитуд (ММА), получим две независимых системы уравнений связанных волн для амплитуд плосковолновых составляющих УС перпендикулярных Еs и тангенсальных Еp компонент поля Е, определяемых с точностью до Δε/ε0 [79]:
dÅ |
0 ( y, Dq0 , Dl) = G0i ∑ n j ( y) × Å ij ( y, Dq j , Dl) ×ei× K j (Dq j ,Dl) y |
|
|||||||||||
|
i |
|
|
|
|
|
|
N |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||
|
|
|
dy |
|
|
|
|
j =1 |
|
|
|
|
|
dÅ |
i |
( y, Dq |
, Dl) |
|
|
|
|
,Dl) y |
|
||||
|
|
1 |
|
1 |
|
|
|
= G1i n1 ( y) × Å 0i ( y, Dq0 |
, Dl) ×e-i× K1 (Dq1 |
|
|||
|
|
|
dy |
|
|
|
, |
(3.23) |
|||||
|
|
|
|
|
|
|
|
|
|
||||
............................................................ |
|
|
|
|
|||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
dÅ |
i |
|
( y, Dq |
N |
, Dl) |
= GNi nN ( y) × Å 0i ( y, Dq0 |
, Dl) ×e-i× K N |
(DqN ,Dl) y |
|
||||
|
|
N |
|
|
|
|
|
||||||
|
|
|
dy |
|
|
|
|
|
|||||
|
|
|
|
|
|
|
|
|
|
|
|
||
где индекс i=s |
соответствует взаимодействию s– |
компонент, i=p – |
|
р– |
|
компонент; |
||||||||||||||||
Gsj |
= i k /(2cos(θ j )), |
G jp |
= −i k /(2 cos2 (θ j )), |
Δθj=Δθj(Δθ0 , |
λ ), отклонение |
длины волны |
||||||||||||||||
считывания от условий Брэгга. |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||||||
|
Принимая |
во |
внимание |
граничные |
|
условия |
E0s, p ( y = 0, Dq0 , Dl) = Ers, p (Dq0 , Dl) , |
|||||||||||||||
Es, p (y = d, θ , |
|
λ) = 0 |
и |
учитывая |
малость |
амплитуд гармоник n |
(y)δn |
p |
<< |
n |
st |
, решение |
||||||||||
j |
|
0 |
|
|
|
|
|
|
|
|
|
|
|
|
|
j |
|
|
|
|
||
системы (3.23) запишем с помощью метода возмущений [80] в матричном виде: |
|
|
|
|||||||||||||||||||
|
|
|
|
|
|
|
|
E j ( y, Dq j , Dl) = Tj ( y, Dq j , Dl)Er (Dq0 , Dl) , |
|
|
|
|
|
(3.24) |
||||||||
где |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
E jp ( y, θ j , |
λ) |
|
|
Erp (Dq0 , Dl) |
|
|
|
|
|
|
|||||
|
E j ( y, θ j , λ) = |
|
s |
|
|
|
|
|
|
s |
|
|
|
|
|
|
|
|||||
|
|
θ j , |
, Er |
(Dq0 ) = |
, |
|
|
|
|
|
|
|||||||||||
|
|
|
|
|
|
|
E j ( y, |
λ) |
|
|
|
Er (Dq0 , Dl) |
|
|
|
|
|
|
||||
|
T |
( y, θ |
|
, |
|
T p ( y, θ |
|
, λ) |
0 |
|
|
, j=0..N; компоненты Tj |
p и Tj |
s являются |
||||||||
|
|
λ) = |
j |
0 |
j |
|
T s ( y, |
θ |
|
|
||||||||||||
|
j |
|
j |
|
|
|
|
|
|
j |
, λ) |
|
|
|
|
|
|
|
|
|||
|
|
|
|
|
|
|
|
|
|
|
j |
|
|
|
|
|
|
|
|
|
|
|
когерентными передаточными функциями для p и s – |
компонент и задаются выражениями |
|||||||||||||||||||||
¥ |
)2k , |
|
T0i ( y, θ0 , λ) = ∑T0i,2k ( y, θ0 , λ) (δnp / nst |
||
k =0 |
|
(3.25) |
¥ |
|
|
|
)2k +1 , |
|
Tji ( y, θ j , λ) = ∑Tji,2k +1 ( y, θ j , λ) (δnp / nst |
||
k =0
75
а коэффициенты T i |
, T i |
определяются рекуррентными соотношениями |
0,2k |
j,2k +1 |
|
|
|
|
|
d − y N |
|
|
|
|
T0,i 2k ( y, θ0 , λ) = G0i ∫ ∑n′j ( y′)Tji,2k −1 ( y′, θ j , λ) ei K j ( θ j , λ) y′ dy′, |
||||||||
|
|
|
|
0 j =1 |
|
|
(3.26) |
|
|
|
|
|
|
d − y |
|
|
|
|
|
|
|
|
|
|
|
|
Tji,2k +1 ( y, θ j , λ) = Gij |
∫ n′j ( y′)T0,i 2k ( y′, θ0 , λ) e−i K j ( θ j , λ) y′ dy′, |
|||||||
|
|
|
|
|
0 |
|
|
|
T i |
( y, θ |
0 |
, λ) = 1, |
n′ |
( y) = n |
( y) / δn |
p |
. |
0,0 |
|
|
j |
j |
|
|
||
Таким образом, соотношения (3.24)– (3.26) являются решением самосогласованной дифракционной задачи и позволяют исследовать дифракционную эффективность и селективные свойства ПГДР при произвольных амплитудах и количестве гармоник пространственно– неоднородной решетки показателя преломления, а также при произвольных поляризации и амплитудно– фазовом распределении считывающего пучка.
3.2.1 Дифракционная эффективность и селективные свойства
Дифракционную эффективность в j– том порядке определим как отношение потока энергии j– го дифрагированного пучка к потоку энергии считывающего пучка в направлении нормали к границе раздела сред (вдоль оси у). Учитывая, что каждой монохроматической волне с комплексной векторной амплитудой Ej можно сопоставить вектор Пойтинга S j = c /(2p) N j ∫ (E j × E j * )d (Dq j ) , запишем выражение для дифракционной эффективности в
виде
hd j = (S j × y0 ) /(S0 × y0 ) .
где у0 – единичный вектор вдоль оси y, Nj – нормаль вдоль оси j– го пучка (см. рис.3.11). Для исследования зависимости дифракционной эффективности в j– том (j=1..N)
дифракционном порядке от угла падения считывающего пучка и селективных свойств дифракционной решетки необходимо получить соотношения для модулей векторов фазовой расстройки ∆Kj = ∆K0j +∆Kj′ в (3.23) от θ0 и ∆θ0 в явном виде. Для этого запишем векторные соотношения фазового синхронизма для центральной угловой компоненты УС (∆θ0=0) j– го дифрагированного пучка (см. рис. 3.11, для j=1,2)
|
|
|
′ |
′ |
|
K0 j , |
(3.27) |
|
|
k j = k0 + jK1 + |
|||||
и текущей компоненты УС (∆θ0≠0) |
|
|
|
|
|||
|
|
k j = k0 + jK1 + DK0 j + DK ′j . |
(3.28) |
||||
Далее запишем (3.27) для проекций векторов на оси координат. Из уравнения для |
|||||||
проекций на ось х, определим угловые положения центральных компонент УС |
|
||||||
|
p |
k0 x¢ |
+ K jx |
|
|
|
|
q j = |
2 |
- arccos |
|
|
|
. |
(3.29) |
|
k |
||||||
|
|
|
|||||
|
|
|
|
|
|
|
|
q j (q0 , g) = p / 2 - arccos[( j -1) sin(q0 )- j sin(q0 - 2g)] |
|
||||||
Подставляя (3.29) в уравнение для проекций на ось у, найдем модули векторов начальной фазовой расстройки:
76
K |
0 j |
= k′ |
− jK |
1y |
− k′ |
, |
|
|
|
jy |
|
0 y |
|
|
|
||
DK |
0 j (q0 , g) = k{cos(p - q j ) - j cos(p - q0 + 2g) + ( j -1) cos(q0 )}. |
(3.30) |
||||||
Используя разложение векторов k j |
= k′j + Pj θ j k′j в ряд Тейлора вблизи направлений |
|||||||
k′j в (3.28) и учитывая (3.27) |
и (3.29), |
определим фазовые расстройки |
K ′j и угловые |
|||||
положения ∆qj текущих компонент УС дифракционных порядков, взаимодействующих с угловой компонентой считывающего пучка, отклоненной от его оси на угол ∆q0:
DK ¢j |
(Dq0 ) = -Dq0 k0¢ |
(P0 × N j ) |
, |
DK ¢j |
(Dq0 ) = -Dq0 k0¢ |
cos(p / 2 - q0 - q j ) |
, |
(3.31) |
|
|
|||||||
|
|
(q × N j ) |
|
|
cos(g + q j ) |
|
||
Dq |
|
(Dq |
|
) = Dq |
|
P0 x |
, |
Dq |
|
(Dq |
|
) = Dq |
|
cos(q0 ) |
. |
(3.32) |
||
|
|
|
|
|
|
|
||||||||||||
|
j |
|
0 |
|
0 P |
|
j |
|
0 |
|
0 cos(q |
j |
) |
|
|
|||
|
|
|
|
|
|
jx |
|
|
|
|
|
|
|
|
|
|
||
где Рj – единичный вектор, указывающий направление смещения волнового вектора j– ой
волны от начального направления k ′j |
(см. рис.3.11б), q=K1/|K1|. |
|
||||
При отклонении длины волны считывания от условий Брэгга на Dl используем, как и |
||||||
ранее, разложение в ряд Тейлора вблизи начальных направлений k ′j |
|
|||||
k0 = k0¢ + 2pnN 0 |
λ |
, |
k j = k¢j + 2pnN j |
λ |
+ Pj Dq j k ¢j , |
(3.33) |
|
|
|||||
|
l2 |
|
l2 |
|
||
где j=1..N.
Тогда, подставляя (3.33) в (3.28) и учитывая (3.27) и (3.29), можно получить следующие соотношения
|
DK ¢ |
(Dl) = k |
|
¢ |
Dl |
N |
|
{1- ( N |
|
× N |
|
)}, |
|
|
|
|
|||||||
|
0 |
|
|
0 |
|
|
|
|
|
||||||||||||||
|
|
j |
|
|
|
|
l |
|
|
|
j |
|
|
|
|
j |
|
|
|
|
(3.34) |
||
|
|
|
|
|
|
Dl |
|
|
|
|
|
cos(p - q0 - q j ) |
|
|
|
||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||||||||
|
DK ¢ |
(Dl) = - |
¢ |
|
1 |
- |
|
|
|
|
|||||||||||||
|
|
|
k |
|
|
|
|
|
|
|
|
|
|||||||||||
|
|
|
|
|
|
|
|
|
|
, |
|
|
|
||||||||||
|
|
j |
|
|
|
l |
|
0 |
|
|
|
cos(p - q j ) |
|
|
|
|
|||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||||||
|
Dl N 0 x - N jx |
|
|
|
|
|
|
|
|
|
|
Dl cos(p + q0 ) + cos(p + q j ) |
|
||||||||||
Dq j (Dl) = |
|
|
, |
|
Dq j |
(Dq0 ) = |
l |
|
|
|
|
(3.35) |
|||||||||||
l |
|
P |
|
|
|
|
cos(q ) |
. |
|||||||||||||||
|
|
|
jx |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
j |
||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||
3.2.2 Дифракционные свойства ОГДР в поглощающем ФПМ
На практике информация о кинетике профиля решеток определяется из кинетики дифракционной эффективности hd(t), которая, являясь интегральной характеристикой, учитывает неоднородность профиля n j ( y) .
На рис. 3.12 приведены графики динамики дифракционных характеристик: эффективности дифракции ηd(τ) (а,б), ширины полосы пропускания 2D0,5 (в,г), уровня первого бокового лепестка hd_SL [dB] (д,е). Все характеристики рассчитаны для Cn=2, s=1,
δnp=0.004, d=20мкм, θ0=θ1=100, b=0.25 (а,в,д), b=5 (б,г,е). На рис. 3.12 использованы следующие обозначения: кривая 1 – ad »0Неп и m0=1, кривая 2 – ad »2Неп и m0=1, кривая 3
– ad=4Неп и m0=1.
Увеличение поглощения в области b<1 при m0=1, как видно из рис. 3.12 а (кривые 1– 3), приводит к увеличению времен достижения максимального и стационарного значений hd записываемой ДС и величины hd на стационарном уровне, в то время как величина
77
максимального значения ηd уменьшается. Уменьшение максимального значения обусловлено |
||||||||||
снижением локального контраста my, что приводит к снижению амплитуды ДС, а увеличение |
||||||||||
стационарного уровня можно объяснить увеличением bу, приводящее к затягиванию |
||||||||||
кинетики и увеличению вклада диффузии мономера в процесс полимеризации. |
|
|
||||||||
При b>1 диффузия мономера происходит достаточно быстро, поэтому увеличение |
||||||||||
поглощения приводит к дальнейшему увеличению локального by>b, что сказывается на |
||||||||||
замедлении полимеризации. |
|
|
|
|
|
|
|
|||
|
0,04 |
|
2 |
|
0,45 |
|
|
|
|
1 |
|
|
|
|
|
|
|
|
|
|
|
ед. |
0,03 |
|
3 |
ед. |
0,30 |
|
|
|
|
2 |
|
|
|
|
|
|
|
||||
отн. |
0,02 |
|
1 |
отн. |
|
|
|
|
|
3 |
|
|
|
|
|
|
|
||||
τ), |
0,01 |
|
|
τ), |
0,15 |
|
|
|
|
|
( |
|
|
( |
|
|
|
|
|
|
|
d |
|
|
d |
|
|
|
|
|
|
|
η |
|
|
|
η |
|
|
|
|
|
|
|
0,00 |
1 |
2 |
3 |
0,00 |
10 |
20 |
30 |
40 |
5 |
|
0 |
0 |
||||||||
|
|
τ, отн. ед. |
|
|
|
τ, отн. ед. |
|
|
||
а) |
|
|
|
б ) |
|
|
|
|
|
|
|
8,0 |
|
|
|
8,0 |
|
|
|
3 |
|
|
|
|
3 |
|
|
|
|
|
|
|
. |
7,5 |
|
. |
7,5 |
|
|
|
|
|
|
|
|
|
|
|
|
|
||||
ед |
|
|
ед |
|
|
|
|
|
||
|
|
|
|
|
|
|
2 |
|
||
отн. |
7,0 |
|
|
отн. |
7,0 |
|
|
|
|
|
|
|
|
|
|
|
|
||||
|
|
|
|
|
|
|
|
|
||
), |
6,5 |
|
2 |
), |
6,5 |
|
|
|
1 |
|
(τ |
|
(τ |
|
|
|
|
||||
|
|
|
|
|
|
|
|
|
||
0.5 |
6,0 |
|
|
0.5 |
6,0 |
|
|
|
|
|
2 |
|
1 |
2 |
|
|
|
|
|
||
5,5 |
|
5,5 |
|
|
|
|
|
|||
|
1 |
2 |
3 |
10 |
20 |
30 |
40 |
5 |
||
|
0 |
0 |
||||||||
|
|
τ, отн. ед. |
|
|
|
τ, отн. ед. |
|
|
||
в) |
|
|
|
г) |
|
|
|
|
|
|
|
-10 |
1 |
|
|
-10 |
|
|
1 |
|
|
|
|
2 |
|
|
|
|
|
2 |
|
|
дБ |
-20 |
|
|
дБ |
-20 |
|
|
|
|
|
|
|
|
|
|
|
|
||||
|
3 |
|
|
|
|
|
|
|
||
(τ), |
|
|
(τ), |
|
|
|
|
|
|
|
-30 |
|
|
-30 |
|
|
|
|
|
||
d SL |
|
|
d SL |
|
|
3 |
|
|
||
|
|
|
|
|
|
|
|
|||
η |
|
|
|
η |
|
|
|
|
|
|
|
-40 |
1 |
2 |
3 |
-40 |
10 |
20 |
30 |
40 |
50 |
|
0 |
0 |
||||||||
|
|
τ, отн. ед. |
|
|
|
τ, отн. ед. |
|
|
||
д) |
|
|
|
е) |
|
|
|
|
|
|
|
|
|
Рисунок 3.12 – |
Дифракционные характеристики ОГДС |
|
|
|
|||
Для m0>0.5 для отражающих решеток, поглощение приводит к уменьшению дифракционной эффективности, особенно для b>1, за счет увеличению влияния поглощения на локальный контраст записываемой интерференционной картины.
Для b<1 (рис3.12 в,д) все изменения дифракционных характеристик обусловлены в основном только неоднородностью профиля, т.к. влияние амплитуды ДС и следовательно
78
ηd(τ,0) на данные характеристики, как видно из зависимости, показанной кривыми 1 на рис..3.12 (в,д) (равномерный профиль, т.к. αd=0 Неп), пренебрежимо мало в сравнении с влиянием неравномерности профиля. Из чего можно заключить, что неоднородность профиля приводит к существенному увеличению полосы пропускания и снижению уровня первого бокового лепестка. Первое обусловлено уменьшением эффективной длины ДС, за счет уменьшения амплитуды ДС от y/d=0.5 к y/d=0 и от y/d=0.5 к у/d=1, что также можно характеризовать как уменьшение скорости нарастания амплитуды на периферийных участках ДС, приводящее в свою очередь и к уменьшению уровня боковых лепестков.
Присутствие ЖК при b<1 приводит к уменьшению спада ДЭ на стационарном уровне,
атакже отсутствию локальных максимумов (минимумов) в 2 0,5(τ), Δηd_SL(τ).
Вотсутствии поглощения для b>1 (рис.3.12 (г,е), кривые 1) можно проследить влияние амплитуды ДС на ее селективные свойства, т.к. пространственный профиль в данном случае является равномерным. Сравнение кривых 2, 3 с 1 на рис.3.12 (г,е) показывает, что поглощение приводит к дальнейшему увеличению полосы пропускания и снижению уровня первого бокового лепестка, что можно связать с локализацией профиля ДС в области 0.3<y/d<0.7 (см. рис. 3.12 б), приводящей к уменьшению эффективной длины взаимодействия света с решеткой. Данное влияние поглощения на селективные свойства аналогично случаю с b<1, т.к. длина область эффективного взаимодействия остается практически без изменения для всех b.
Разрыв в зависимости бокового лепестка при αd=4Неп (рис. 3.12 (е), кривая 3), обусловлен тем, что первый лепесток уменьшается по величине и ширине под влиянием расширяющегося основного лепестка и неподвижного второго (см. сплошная линия на рис. 3.12 (в). В некоторый момент, второй боковой лепесток становится первым, как видно из сравнения сплошной и штрих линий на рис. 3.12 (в), а затем проявляется и заплывание минимумов характеристики (сплошная линия на рис. 3.12 в).
Перейдем к рассмотрению влияния фотоиндуцированного изменения поглощения на дифракционные характеристики ОГДС . Как и для ПГДС промоделируем случай, когда за время записи затухание изменяется от 4 Неп до нуля. На рис.3.13 представлены дифракционные характеристики ηd(τ), 2 0,5(τ), Δηd_SL(τ) для Cn=2, s=1, δnp=0.004, d=20мкм,
θ0=θ1=100, α2d≈0 и α1d=0 (кривые 1), α2d=4 Неп и α1d=0 (кривые 2), αd=4 Неп и α1d=1.4 Неп
(кривые 3).
η (τ), отн.ед. d
a)
0,20 |
|
|
|
8,0 |
|
2 |
|
-10 |
|
1 |
|
|
7,5 |
|
|
|
|
0,15 |
|
),отн.ед. |
|
|
|
|
||
3 |
|
|
|
3 |
τ),(дБ |
-20 |
||
|
|
7,0 |
|
|||||
|
|
|
|
|
||||
0,10 |
2 |
|
|
|
|
|
|
|
|
|
6,5 |
|
|
|
|
||
|
|
|
(τ |
|
|
SLd |
-30 |
|
0,05 |
|
|
|
|
|
|||
|
|
0.5 |
|
|
1 |
η |
||
|
|
6,0 |
|
|
||||
|
|
|
2 |
|
|
|
|
|
0,00 |
|
|
5,5 |
|
|
|
-40 |
|
5 |
10 |
15 |
5 |
10 |
|
|||
0 |
0 |
|
0 |
|||||
|
τ, отн. ед. |
|
|
τ, отн. ед. |
|
|
||
|
|
|
б) |
|
|
|
в) |
|
|
|
|
|
Рисунок 3.13 |
|
|
|
|
1
3
2
5 10
τ, отн. ед.
Как видно из рис.3.13, фотоиндуцированное изменение поглощения приводит все характеристики к некоторому среднему между случаем с постоянным затуханием и случаем без затухания. Данное поведение присуще для всей области параметра b.
3.2.3 Дифракционные свойства с учетом самодифракции при малом контрасте
На рис.3.4 приведены модуль (а) и аргумент (б) нормированного пространственного профиля ДС, рассчитанные по полученным выражениям На рис. 3.4 (в) приведены кривые
79
