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

Голографические фотонные структуры в наноструктурированных материалах

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

 

M j (r,t)e

iK

j

r

+ M120 (r,t)e

iK

20

r

+ M121

(r,t)eiK

21r +

 

 

 

N

 

 

 

 

 

 

 

 

 

 

 

M = 0.5 j =0

 

 

 

 

 

 

 

 

 

 

 

 

,

 

 

 

+ M102 (r,t)eiK 02 r + M112 (r,t)eiK12 r + c.c.

 

 

 

 

 

 

 

.

(3.11)

 

 

N

 

 

 

 

 

 

 

 

 

 

 

 

 

 

iK

r

+ n120 (r,t)eiK 20 r + n121 (r,t)eiK

 

 

 

 

n j (r,t)e

 

 

j

 

21r +

 

 

n = nst + 0.5 j =0

 

 

 

 

 

 

 

 

 

 

 

.

 

 

 

 

+ n102 (r,t)eiK 02 r + n112 (r,t)eiK12 r + c.c.

 

 

 

 

 

 

 

 

где Mj, nj – амплитуды

гармоник

 

концентрации

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

преломления

основной решетки, M20, M21, M02, M12,

n20, n21, n02, n12

амплитуды первых гармоник

концентрации мономера и показателя преломления ДДР.

 

 

 

 

Принимая во внимание малость амплитуд ДДР, можно пренебречь их влиянием на гармоники основной решетки и разделить рассматриваемую задачу на две части. Это позволяет в первой части найти решение для пространственных гармоник основной решетки и, используя его, решить задачу формирования ДДР.

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

El = iG (E n (τ)ei K 0l y + E n (τ)ei 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 0l y

+ E1nl −1 (τ)e

i K

0l 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)}ei K

 

 

 

 

 

 

 

 

 

H(τ,p)=L{H(τ,y)}=

 

 

 

 

 

 

где введены

обозначения

 

= L M 0 (t) {nl (t) +

 

0 l ydy¢ ,

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

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

()eF1τ′

dt¢ -

 

2

H (, p)eF1τ′ d.

 

(3.15)

 

 

dn p

 

 

 

p

dnp

 

 

 

 

 

 

 

 

 

 

 

0

 

 

 

0

 

 

 

 

 

Подставляя полученное решение (3.15) во второе уравнение из (3.14) и используя обратное интегральное преобразование Лапласа по пространственной координате у, решение для n1l0(t,y) запишем в виде:

 

τ

d y

 

 

nl 0

 

Q(, y¢) × H 0

 

(3.16)

(t, y) = F2 Q(, y) +

(, t, y - y¢) d,

 

 

0

 

 

 

0

 

 

 

 

 

 

2k

τ

H (, y)

eF1τ′ 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, y) = 1+

 

R(t¢¢)dt¢¢×J1 2

F2

(d - y) ∫

R(t¢¢)dt¢¢

,

 

 

 

 

d - y τ′

 

 

 

τ′

 

 

 

R(t) = C

n

2

 

2k

 

F1τ

τ

M

0

()

 

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 (τ))ei K0l ydy, J1(x) – функция Бесселя.

0

Таким образом, выражение (3.16) определяет временную динамику пространственного распределения амплитуды ДДР с вектором решетки Kl0, образованной записывающими волнами и дифрагированной на l– той пространственной гармонике ОГДР. Из решения видно, что ДДР в отличии от гармоник основной решетки имеет пространственно неоднородные распределения амплитуды и фазы вдоль глубины решетки даже в отсутствии поглощения.

Проведя аналогичную процедуру решения для ДДР с вектором решетки Kl1 можно получить следующие решения:

71

l1

 

τ

(t, y) = F2

 

n1

Q(, y)

 

 

 

 

 

0

 

F

m

τ

где H 0 (, t, y) = 1+

2

0

R(t¢¢)dt¢¢×

d - y

 

τ′

d y

 

 

+ Q(, y¢) × H 0

 

(3.17)

(, t, y - y¢) d,

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(, y) + Q(, y¢) × H0

(, t, y - y¢) d,

 

 

 

0

 

 

 

 

0

 

 

 

F

m

τ

 

 

 

τ

 

 

 

 

где H 0 (, 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 (τ)dy(nl −1 (τ) +

 

nl (τ))ei K0l ydy′ .

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 × (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 ×(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]:

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

 

 

 

 

i

( y, Dq

, Dl)

 

 

 

 

,Dl) y

 

 

 

1

 

1

 

 

 

= G1i n1 ( y) × Å 0i ( y, Dq0

, Dl) ×e-i× K1 (Dq1

 

 

 

 

dy

 

 

 

,

(3.23)

 

 

 

 

 

 

 

 

 

 

............................................................

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

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

(yn

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 ∫ ∑nj ( y′)Tji,2k −1 ( y′, θ j , λ) ei K j ( θ j , λ) ydy′,

 

 

 

 

0 j =1

 

 

(3.26)

 

 

 

 

 

d y

 

 

 

 

 

 

 

 

 

 

Tji,2k +1 ( y, θ j , λ) = Gij

nj ( y′)T0,i 2k ( y′, θ0 , λ) ei K j ( θ j , λ) ydy′,

 

 

 

 

 

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

= kj + Pj θ j kj в ряд Тейлора вблизи направлений

kj в (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мкм, θ01=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мкм,

θ01=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

Соседние файлы в предмете [НЕСОРТИРОВАННОЕ]