Например, Бобцов

Особенности дифракции импульсного излучения фемтосекундной длительности на объемной пропускающей голографической решетке

ГОЛОГРАФИЯ

УДК 535.317.1; 778.38
ОСОБЕННОСТИ ДИФРАКЦИИ ИМПУЛЬСНОГО ИЗЛУЧЕНИЯ ФЕМТОСЕКУНДНОЙ ДЛИТЕЛЬНОСТИ НА ОБЪЕМНОЙ ПРОПУСКАЮЩЕЙ ГОЛОГРАФИЧЕСКОЙ РЕШЕТКЕ

© 2009 г. Н. В. Ионина
Санкт-Петербургский Государственный университет информационных технологий, механики и оптики, Санкт-Петербург
Е-mail: shurch@yandex.ru

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

Ключевые слова:

Коды OCIS: 090.0090, 090.2900

Поступила в редакцию 10.02.2009

В работе численно и экспериментально исследованы некоторые особенности временной и спектральной структуры импульсов немонохроматического излучения при дифракции на объемной полимерной фазовой пропускающей голографической решетке.
При широком использовании дифракционных решеток разного рода в оптике сверхкоротких импульсов, лазерной технике, волноводной оптике [1, 2] представляется, что описания особенностей работы пропускающих решеток в случае сверхкоротких световых импульсов далеко не исчерпывают разнообразия возможных условий и результатов применения. В данной работе рассматривалась объемная пропускающая фазовая голограмма-решетка, формирующаяся в результате модуляции показателя преломления n светочувствительной среды. Такие голограммы обладают высокой угловой селективностью и высокой дифракционной эффективностью, которая для монохроматического излучения может достигать в пределе 100%. При этом по сравнению с отражательными объемными голограммамирешетками объемные пропускающие голографические ре-шетки менее чувствительны к изменению сред-него показателя преломления регистрирующей среды и, соответственно, имеют более высокую стабильность голографических параметров.
Для такой голограммы периодическое изменение показателя преломления в пространстве можно описать зависимостью вида

n(r) = n– + n1cos(K·r),

(1)

где n– – средний показатель преломления регистрирующей среды, n1 – амплитуда модуляции показателя преломления, K – вектор решетки; |K| = 2π/d; d – период решетки, r – радиус-вектор точки среды. Предполагается также, что параметры голограммы-решетки не изменяются при прохождении светового пучка.
При анализе процесса дифракции на голограмме сверхкоротких импульсов излучения использовалось спектральное представление светового поля, что продуктивно для сверхкоротких импульсов, ширина спектра которых не позволяет считать их монохроматическими. Соответственно комплексная амплитуда падающей на голограмму-решетку волны представлялась в виде
¥
òEin(r,t) = Ein(ω)exp(-jr¢(ω)⋅ r)exp(jωt)dω, (2) -¥
где Ein(ω) – спектр Фурье комплексной амплитуды падающей волны Ein(r, t), ω – круговая частота излучения, r′(ω) – волновой вектор падающего на голограмму излучения.
Рассматривались собственно дифракционные явления (влияние преломления и френелевского отражения от поверхностей исключалось), для чего средний показатель преломления голограм-

34 “Оптический журнал”, 76, 6, 2009

мы полагался равным показателю преломления окружающей среды.
Для описания воздействия фазовой объемной решетки вида (1) на проходящий световой импульс использовалось представление частотных передаточных функций объемных голограмм, полученных в теории связанных волн [3]. При этом поле внутри голограммы представлялось в виде суммы волн нулевого порядка R(r, t) и дифрагированной S(r, t), каждая из которых, в свою очередь, описывалась суперпозицией плоских монохроматических волн, распространяющихся в одном направлении [4],

E(r,t) = R(r,t) + S(r,t) =
¥
= ò R(ω,z)exp(-jr(ω)⋅ r)exp(jωt)dω + -¥ ¥
+ ò S(ω,z)exp(-js(ω)⋅ r)exp(jωt)dω, -¥

(3)

где R(ω, z), S(ω, z) – комплексные амплитуды плоских компонентов волн нулевого порядка и

дифрагированной; z – координата по оси оz, нормальной к поверхности голограммы; r(ω) – волновой вектор волны нулевого порядка, |r(ω)| = = β(ω) = ωn/c, с – скорость света в вакууме, s(ω) – волновой вектор дифрагированной волны. Вектор Е нормален плоскости XOZ. Внутри голограммы [3] s(ω) = r(ω) – K (см. рис. 1).
Учитывая линейность волновых уравнений
теории связанных волн для объемных голограмм,
комплексные амплитуды плоских компонентов
можно записать в виде [5]

R(ω,z) = ÏR (ω,z)Ein (ω), S(ω,z) = ÏS (ω,z)Ein (ω).

(4)

Здесь ПR(ω, z), ПS(ω, z) – частотные передаточные функции объемной пропускающей фазовой

голограммы [6]:

ÏR (ω,z) = exp(-jξ(ω)) êéëêcos ëêéê ξ(ω)2 + ν(ω)2 úúùû + + j éêëêξ(ω)sin ξ(ω)2 + ν(ω)2 ξ(ω)2 + ν(ω)2 úûùúùúûú,

(5)

ΠS (ω,z) = -j ρz (ω)/σz (ω) exp(-jξ(ω))ν(ω)sin ξ(ω)2 + ν(ω)2 ξ(ω)2 + ν(ω)2 ,

(6)

где брэгговская расстройка и сила решетки определяются соответственно как

ξ(ω) = ëéêβ(ω)2- σ(ω) 2 úûùT/(4 σz (ω) ),
( )ν(ω) = ωβ(ω)n1T 2c σz (ω) ρz (ω) .

(7)

Здесь подстрочный индекс z означает проекцию вектора на ось ОZ, Т – толщина голограммы.
Численные эксперименты проводились для условий, соответствующих реальной голографической регистрирующей среде (полиметилметакрилат с фенантренхиноном) [7, 8] и параметрам излучения лазера на сапфире с титаном (длина волны излучения λ = 810–815 нм, длительность импульсов излучения 20–30 фс). Спектр Фурье комплексной амплитуды падающей волны задавался в виде
( )Ein = exp êêëé- ω -2,32´1015 2/2(0,075´1015 )2 úûúù (8)
и соответствовал импульсу с центральной длиной волны λ0  812 нм и длительностью τ  27 фс по уровню 1/е максимальной интенсивности (см. рис. 1б).

“Оптический журнал”, 76, 6, 2009

Интенсивность, отн. ед.

(a)
N 1 Ti:S

X K
О

Z 2 –K

3 ASP-100M

2000 1000

1

(б) 2

4 ПК

0 800

850

Длина волны, нм

Рис. 1. а – cхема эксперимента. 1 – фемтосекундный лазер, 2 – объемная пропускающая голограмма с векторной диаграммой реконструкции, 3 – автоматизированный спектрометр ASP-100M, 4 – компьютер, N – нормаль к поверхности голограммы; б – спектр падающего на голограмму излучения. 1 – экспериментальный, 2 – расчетный.

35

Интенсивность, отн. ед. Интенсивность, отн. ед.

2000 1000

21

(a)

0 770,092
2000

823,261

876,443 (б)

1000
0 749,892

1 2

803,055

856,232

Длина волны, нм

Рис. 2. Спектр излучения в нулевом по-
рядке дифракции на выходе голограммы.
n1 = 0,0001, d = 0,00295 мм, Т = 2,52 мм, λ0 = 812 нм. а – угол падения θ = 0,094 рад, б – θ = 0,092 рад, 1 – расчет, 2 – экспе-
римент.

Средний показатель преломления регистрирующей среды принимался равным n~  1,49 [8].
При расчетах учитывалась материальная диспер-
сия среды как дисперсия полиметилметакрилата в области частот [2×1015, 4,1×1015] с–1 [9]

n(ω) =1,48474-0,00441ω /1015 +
( )+ 0,00191 ω /1015 2 .

(10)

При указанных условиях проведено численное моделирование эволюции формы и спектра сверхкороткого (немонохроматического) импульса при дифракции излучения на объемной пропускающей фазовой решетке.
Полученные для различных значений параметров n1 и Т зависимости говорят о сложном характере связи между параметрами голограммы и формой и спектром импульсов как дифрагированного, так и прямопроходящего излучения (см. рис. 2, 3). Как видно из графиков, спектр импульсов может претерпевать существенные изменения (и в нулевом, и в первом порядках дифракции), которые определяются характером частотной передаточной функции голограммырешетки. Изменения спектров сопровождаются сложным изменением формы огибающей импульсов излучения.

Ir(t), отн. ед. 0,4
0,2
0 1,2

1,23

Ir(ω), отн. ед. (а) 0,8
0,4

1,26

0 2,1

2,3

(в) 2,5

Is (t), отн. ед. 0,8

Is(ω), отн. ед. (б) 1

(г)

0,4 0,5

00

1,2

1,23 t 1011, c 1,26

2,2

2,35 

10–15, c–1

2,5

Рис. 3. Форма и спектр импульса на выходе голограммы. а, в – в нулевом порядке дифракции; б, г – в первом порядке дифракции, n1 = 0,0005, d = 0,00295 мм, Т = 2,45 мм.
36 “Оптический журнал”, 76, 6, 2009

Однако при некоторых значениях n1и Т в структуре огибающей импульса в первом поряд-
ке дифракции Is(t) появляется определенная регулярность. А именно, формируются два фемто-
секундных максимума (пичка) с длительностью (τout), несколько превышающей длительность импульса падающего излучения, и имеющих
один относительно другого временную задержку Δt (см. рис. 4а).Спектр импульса Is(t) в этом случае приобретает характерный вид с глубокой
модуляцией по интенсивности (рис. 4в). Как
показывают расчеты, задержка (временной интервал) Δt изменяется при изменении угла падения излучения на голограмму (θ). График зависимости Δt(θ) приведен на рис. 5.
Для оценки корректности расчетной модели
был проведен контрольный эксперимент по диф-
ракции фемтосекундного импульса на объемной
пропускающей фазовой голограмме-решетке с известными параметрами [8, 10]: Т  2,52 мм, n1  10–4, n~  1,49, d  2,95×10–3 мм. Условия дифракции были выбраны так, чтобы получить
характерные и ожидаемые изменения параме-
тров импульса, фиксируемые прямыми измере-
ниями:
а) минимум в спектре нулевого порядка,
определяемый контуром спектральной селектив-
ности объемной голограммы и соответствующий

условиям Брэгга для заданной длины волны и заданной решетки [11];
б) “смещение” этого минимума по спектру при изменении угла падения излучения на голограмму и соответствие этого “смещения” изменению условий Брэгга.
Схема экспериментальной установки приведена на рис. 1. Источником излучения в эксперименте служил фемтосекундный лазер 1 на сапфире с титаном “Фемос-2” (Новосибирского государственного университета). Средняя мощность его излучения 80–100 мВт, энергия в импульсе до 100 нДж, частота повторения импульсов длительностью 15–50 фс порядка 80 МГц. Пучок излучения диаметром примерно 2 мм в поперечном сечении имеет гауссово распределение интенсивности. Расходимость пучка – порядка десятых долей миллирадиана. Излучение лазера с центральной длиной волны λ0 ≈ 812 нм и длительностью импульса τ ≈ 27 фс (по уровню 1/е максимальной интенсивности) падает на объемную пропускающую фазовую голограммурешетку 2, установленную на поворотном столике, обеспечивающем изменение углового положения с точностью до 0,05 мрад. Излучение нулевого порядка дифракции с выхода голограммы подается на вход автоматизированного спектрометра ASP-100M (ООО “Авеста-Проект”),

Is(t), отн. ед. 0,4

Is(ω), отн. ед. 03
(а)

(в)

0,2

0 1,25
Ir (t), отн. ед. 0,4

1,27

1,29 (б)

0 2,2 Ir(ω), отн. ед.
0,8

2,36

2,5 (г)

0,2 0,4

0 1,25

1,27 t 10–11, c 1,29

0 2,1

2,3 

1015, c–1

2,5

Рис. 4. Формирование в первом порядке дифракции пары фемтосекундных пичков. а), б) – форма и

спектр импульса на выходе голограммы в первом порядке дифракции; в), г) – форма и спектр им-

пульса на выходе голограммы в нулевом n~  1,49, d  0,00295 мм, θ  0,093 рад.

порядке

дифракции;

λ0



812

нм,

Т



2,52

мм,

n1



10



4,

“Оптический журнал”, 76, 6, 2009

37

t, фс 202

194

186 0,087

0,090

0,093

0,096 , рад

Рис. 5. Зависимость временного интервала Δt от угла падения излучения на голограмму.

связанного с персональным компьютером (ПК). Результаты измерения распределения интенсивности в исследуемом спектре излучения выводятся в виде графика на монитор ПК, а также формируют соответствующий файл с массивом численных данных, позволяющий проводить обработку результатов измерений. Результаты измерений и расчетов распределения интенсивности в спектрах нулевого порядка дифракции фемтосекундного импульса излучения представлены на рис. 2.
Приведенные графики показывают, что в экспериментальных спектрах нулевого порядка дифракции наблюдаются, как и предполагалось, минимум интенсивности и смещение положения этого минимума в сторону более высоких или низких частот в зависимости от угла падения излучения на голограмму. При этом экспериментальным кривым соответствуют расчетные зависимости.
По данным экспериментальных графиков и расчетных кривых (рис. 2.) были определены положения спектральных минимумов (ωmin) в зависимости от угла падения излучения на голограмму (см. таблицу). Как видно из приведенных в таблице данных расчетные значения соответствуют экспериментальным.

Таким образом, при численном анализе воздействия объемной фазовой пропускающей полимерной голограммы-решетки на фемтосекундный (немонохроматический) световой импульс продемонстрировано, что временная и спектральная структуры импульса на выходе голограммы, как в 0-м, так и в 1-м порядках дифракции в целом сложным образом зависят от параметров решетки. При этом соответствие результатов расчетов и данных контрольного эксперимента позволяют говорить о корректности расчетной модели.
В то же время при моделировании установлено, что для определенных значений параметров решетки возможно сформировать в 1-м порядке дифракции импульс с определенной формой огибающей в виде пары фемтосекундных пичков, временной интервал между которыми меняется в зависимости от угла падения излучения на голограмму. Такая особенность может быть, например, использована в системах обработки, передачи и кодирования информации, экспериментах по исследованию свойств различных оптических сред.
Автор выражает благодарность С.Э. Путилину за помощь в проведении экспериментов и Ю.Л. Корзинину за полезные обсуждения.
ЛИТЕРАТУРА
1. Ахманов С.А., Выслоух В.А., Чиркин А.С. Оптика фемтосекундных лазерных импульсов. М.: Наука, 1988. 312 с.
2. Оптика наноструктур. СПб.: Недра, 2005. 361с.
3. Kogelnik H. Coupled wave theory for thick hologram grating // Bell Syst. Techn. 1969. V. 48. № 9. P. 2909–2947.
4. Борн М., Вольф Э. Основы оптики. М.: Наука, 1973. 720 с.
5. Корзинин Ю.Л., Киба И.А. Канализация энергии в объемных пропускающих фазовых решетках // Сб. трудов и официальных материалов V Междунар. конфер. “Голография ЭКСПО-2008”. СПб.: 2008. С. 123–127.

Зависимость положения спектральных минимумов ωmin от угла падения излучения θ

θ,рад

0,0905 0,0919 0,0926 0,0932 0,0939 0,0949 0,0955

Экспериментальные значения ωmin ×10–15, с–1
Расчетные значения ωmin ×10–15, с–1

2,22 2,22

2,23 2,24

2,24 2,25

2,27 2,27

2,29 2,29

2,31 2,31

2,33 2,33

0,0963 2,35
2,35

0,097 2,39
2,39

38 “Оптический журнал”, 76, 6, 2009

6. Solymar L., Cook D.J. Volume Holography and Volume Gratings. London: Academic Press, 1988.
7. Вениаминов А.В., Седунов Ю.Н., Попов А.П., Бандюк О.В. Постэкспозиционное поведение голограмм под влиянием диффузии макромолекул // Опт. и спектр.1996. Т. 81. В. 4. С. 676–680.
8. Андреева О.В., Бандюк О.В., Парамонов А.А., Черкасов А.С., Гаврилюк Е.Р., Андреев П.В. Объемные пропускающие голограммы в полимерной среде с фенантренхиноном // Оптический журнал. 2000. Т. 67. № 12. С. 27–33.
9. Физические величины. Справ. / Под ред. Григорьева И.С., Мейлихова Е.З. М.: Энергоатомиздат, 1991.

10. Ионина Н.В., Андреева О.В., Бандюк О.В., Парамонов А.А. Предэкспозиционный контроль образцов полимерной регистрирующей среды // Сб. трудов V Международной конференции “Прикладная оптика-2002”, СПб., 2002. Т. 2. С. 144–146.
11. Андреева О.В., Беспалов В.Г., Васильев В.Н., Городецкий А.А., Кушнаренко А.П., Лукомский Г.В., Парамонов А.А. Исследование спектральной селективности объемных голограмм с помощью импульсного излучения фемтосекундной длительности // Опт. и спектр. 2004. Т. 96. № 2. С. 190–196.

“Оптический журнал”, 76, 6, 2009

39