CC BY
43
Литература
1. Dahmen J., Hektor J., Perrey R., Ney H. Automatic Classification of Red Blood Cells Using Gaussian Mixture Densities // Proc. Bildverarbeitung für die Medizin. - 2000. - P. 331-335.
2. Costrarido L. Medical Image Analysis Methods: Evaluation Strategies for Medical-image Analysis. - Taylor & Francis, United Stated of America, 2005. - P. 433-471.
3. Kumar B.R., Joseph D.K., Teager T.V.S. Energy Based Blood Cell Segmentation // 14th International Conference on Digital Signal Processing. - DSP, 2002. - 1-3 July. - Santorini, Greece. - V. 2. - P. 619-622.
4. Bamford P. Empirical Comparison of Cell Segmentation Algorithms Using an Annotated Dataset // Proc. IEEE International Conference on Image Processing. - 2003. - V. 2. - P. 1073-1077.
5. Mukheijee D.P., Ray N., Acton S.T. Level Set Analysis for Leukocyte Detection and Tracking // IP. - 2004.
- V. 13. - № 4. - P. 562-572.
6. Park J., Keller J.M. Snakes on the Watershed // IEEE Transactions on Pattern Analysis and Machine Intelligence. - PAMI, 2001. - V. 23. - № 10. - P. 1201-1205.
7. Sinha N., Ramakrishnan A.G. Blood Cell Segmentation Using EM Algorithm // Proc. Third Indian Conference on Computer Vision. Graphics Image Processing (ICVGIP), 2002. - Ahmadabad, India, 2002, December 16-18. - P. 376-382.
8. Kumar R.S., Verma A., Singh J. Color Image Segmentation and Multi-Level Thresholding by Maximization of Conditional Entrophy // International Journal of Signal Processing. - 2006. - V. 3. - № 1. - P. 121-125.
9. Mcinerney T., Terzopoulos D. Deformable models in medical image analysis: A survey // Med Image Anal.
- 1996. - P. 91-108.
10. Cormen T.H., Leisorson C.E., Rivest R.L. Introduction to Algorithms. - MIT Press, 1990. - P. 185-191.
11. Canny J.F. A computational approach to edge detection // IEEE Trans. Pattern Analysis and Machine Intelligence. - 1986. - P. 679-698.
Дырнаев Алексей Валерьевич - Санкт-Петербургский национальный исследовательский университет
информационных технологий, механики и оптики, студент, [email protected]
УДК 535.4, 535-14
СПЕКТРАЛЬНО-ВРЕМЕННАЯ ЭВОЛЮЦИЯ ЭЛЕКТРИЧЕСКОГО ПОЛЯ ТЕРАГЕРЦОВОГО ИМПУЛЬСА ПРИ ДИФРАКЦИИ ФРАУНГОФЕРА НА ЩЕЛИ М.С. Куля, Я.В. Грачев, А.А. Городецкий, В.Г. Беспалов
Экспериментально и теоретически получены временные зависимости амплитуды электрического поля терагерцового электромагнитного импульса из нескольких осцилляций для различных углов дифракции на щели в дальней зоне. Выявлено, что с ростом угла дифракции наблюдаются существенные изменения временной формы электрического поля исследуемого импульса, а также соответствующие изменения спектра. Сравнение результатов исследования показывает хорошее согласование расчетов и экспериментальных результатов. Полученные результаты могут быть использованы в импульсной терагерцовой оптике и спектроскопии.
Ключевые слова: оптика импульсов из малого числа колебаний электрического поля, импульсное терагерцовое излучение, дифракция на щели.
Введение
Терагерцовым (ТГц) обычно называют излучение, лежащее в интервале частот 0,1-10 ТГц, которое по длинам волн находится между инфракрасным и миллиметровым/субмиллиметровым диапазонами, от 0,03 мм до 3 мм. Исследования физики и техники ТГц излучения начались достаточно давно [1, 2], однако только с развитием фемтосекундной оптики и микроэлектроники в данной области наметился значительный сдвиг, возрос интерес к фундаментальным и прикладным исследованиям, были опубликованы первые монографии [3, 4]. Были созданы методы получения и детектирования ТГц излучения с использованием
световых импульсов фемтосекундной длительности при возбуждении носителей заряда в фотополупровод-
(2)
никах и сверхпроводниках; с помощью оптического детектирования в х ' нелинейных средах; при оптическом пробое газов излучением первой и второй гармоник фемтосекундного лазера [5-7].
Одним из фундаментальных явлений оптики является дифракция света, однако исследований пространственно-временной структуры при распространении сверхкороткого широкополосного импульса, обладающего сложной амплитудно-фазовой структурой, в достаточной мере не проводилось [8-11]. Экспериментально дифракцию импульсов из нескольких осцилляций электромагнитного поля для волн в видимом и ИК-диапазонах спектра на практике достаточно сложно изучать, так как практически нет возможности прямой регистрации амплитуды и фазы, и только итерационные процедуры расчета позволяют восстановить исходный временной или пространственный профиль [12, 13].
Экспериментальное изучение дифракции сверхкороткого импульса со сложной амплитудно-фазовой структурой стало возможным с появлением техники регистрации временных форм электрического поля импульса широкополосного ТГц излучения. ТГц волновой фронт при возбуждении фемтосе-кундными импульсами, как правило, представляет собой несколько колебаний электромагнитного поля и обладает сложной амплитудно-фазовой структурой, связанной как с процессом генерации и дальнейшей дифракцией на образованной апертуре, так и с влиянием дисперсии среды [14]. Целью настоящей работы являлось получение экспериментальных результатов и сравнение расчетных зависимостей для спектрально-временных форм электрического поля волны при дифракции импульсов из малого числа колебаний на примере ТГц излучения с широкополосным спектром и длительностью несколько пикосекунд. Особенность работы заключалась в прямой регистрации временных форм дифракции предельно короткого импульса, что было затруднительно в оптике импульсов из большого числа колебаний. Данное исследование может быть полезным при интерпретации результатов измерений в импульсной ТГц оптике и спектроскопии.
Математическая модель
Рассмотрим одномерную задачу об определении распределения электрического поля Е((] дифрагированного излучения в дальней зоне плоской ТГц волны, падающей на бесконечную в направлении оси г щель шириной 2а в системе координат х, у, г. Исходное импульсное электрическое поле Е(() для терагерцовой волны записывается в виде
Е (/) =- • ехр(^), (1)
х х
где t - время; т - длительность импульса.
Комплексный спектр данного поля определяется прямым преобразованием Фурье:
да t — t2
О(у) = [-• ехр(——)• ехр(—/2лу • t)Л , (2)
—I х х2
где V - частота излучения.
Широкополосный спектр ТГц излучения (2) можно представить в виде суперпозиции монохрома, 2лу
тических компонент, каждая из которых описывается волновым числом к =-. Далее для каждой от-
с
дельной монохроматической компоненты можно записать дифракционный интеграл Фраунгофера. Так как щель бесконечна в плоскости г, то волна отклоняется только в плоскости ху, а разложение поля проводится по поперечной координате х. Таким образом, пространственно-временной спектр в дальней зоне дифракции может быть записан в следующем виде:
да t _12 а
О(V,6) = [ — • ехр(——) • ехр(—/2лу • t• [ ехр(—¡к х)ёх , (3)
х х2 у
—да ^ ^ —а
где
ку = — 8Ш(6), (4)
с
с - скорость света в вакууме; а - половина ширины щели; ку - проекция волнового вектора дифрагированной ТГц волны на оптическую ось у; 6 - угол между к у и осью у. С учетом (4) выражение (3) преобразуется в
О^, 6) =
' • а•V • (6)ЛА
81П|--51П(6)
яа•V . --51П(6)
Е(О • ехр(—/2^ • t)Л . (5)
Таким образом, получаем для каждого значения дифракционного угла 9 и частоты соответствующий комплексный спектр дифракции. Временная форма ТГц импульса может быть получена с помощью обратного преобразования Фурье от выражения (5):
да
Е^, 6) = | О(у, 6) • ехр(—/ 2Я . (6)
—да
Используя соотношения (5) и (6), можно рассчитать пространственно-временные картины дифракции импульсного ТГц излучения на щели. На рис. 1 показаны графики спектрального и временного дифракционного распределения для идеального ТГц импульса из одного колебания поля, задаваемого формулой (1). Видно, что с ростом угла дифракции 9 наблюдаются изменения временной формы им-
с
с
пульса, а также соответствующие изменения спектра. Ширина спектра сужается, а центральная частота смещается в область низких частот, при этом появляются вторичные максимумы в спектре дифракции. Во временной области с ростом угла дифракции максимумы амплитуды поля импульса смещаются в противоположных друг от друга направлениях по оси времени, что свидетельствует об искривлении исходного плоского волнового фронта терагерцовой волны при дифракции.
а б
Рис. 1. Спектральное (а) и временное (б) распределение для дифрагированного идеального однопериодного ТГц импульса Графики представлены для различных дифракционных углов: 1 - 9=0 рад;
2 - 9=0,02 рад; 3 - 9=0,03 рад; 4 - 9=0,04 рад. Ширина щели 300 мкм
Схема эксперимента
В схеме (рис. 2) лазерный пучок от фемтосекундного УЪ: KYW лазера (длина волны Х=1040 нм; длительность импульса по уровню половинной ширины х^200 фс; частота следования равна 80 МГц; средняя мощность ^ср=1,2 Вт) с помощью делителя пучка разделяется на пробный пучок и пучок накачки. Пучок накачки, пройдя через управляемую линию оптической задержки, попадает на полупроводниковый кристалл 1пАб, помещенный в постоянное магнитное поле, где происходит генерация ТГц излучения [6, 7]. Затем ТГц излучение коллимируется параболическим зеркалом диаметром 9 см, за которым находится тефлоновый фильтр, пропускающий ТГц излучение и рассеивающий мощное инфракрасное.
Рис. 2. Схема экспериментального исследования дифракции на щели ТГц импульса
Между двумя параболическими зеркалами помещается дифракционная щель из двух тонких алюминиевых пластин. Расстояние между металлическими пластинами выбиралось сопоставимым с длиной волны ТГц импульса для выполнения условия дифракции. Габаритные размеры сконструированной щели были равны апертуре сколлимированного пучка ТГц излучения, распространяющегося после первого параболического зеркала. Дифракционная щель закреплялась на микрометрической подвижке, таким
образом, обеспечивалась точность при перемещении щели поперек пучка, что было эквивалентно перемещению системы регистрации по пространственной картине дифракции.
Оптический импульс пробного пучка по мере распространения по электрооптическому кристаллу С±Те подвергается воздействию постоянного электрического поля ТГц импульса. При коллинеарном распространении в кристалле линейная поляризация пробного пучка изменяется до эллиптической. С помощью четвертьволновой пластинки линейная поляризация переводится в круговую, затем с помощью призмы Волластона пучки с разными состояниями поляризации разделяются и попадают на балансный детектор. Разность интенсивностей на диодах балансной схемы пропорциональна полю ТГц излучения. Затем сигнал с балансного детектора поступает на синхронный усилитель и далее на компьютер.
С помощью линии оптической задержки производится изменение времени пересечения ТГц излучения и пробного пучка в кристалле, и, таким образом, измеряя при различных задержках картину дву-лучепреломления, регистрируется зависимость амплитуды электрического поля ТГц импульса от времени.
Результаты эксперимента и сравнение с расчетными результатами
Для щели размером 5 мм (приблизительно 10 длин волн центральной длины волны ТГц излучения) были измерены зависимости поля £(/) для ТГц импульса в поперечных точках пространственной картины - 0 мм, 5 мм, 10 мм, 15 мм, 20 мм, 25 мм, 30 мм, что соответствовало углам дифракции 0 рад, 0,016 рад, 0,033 рад, 0,049 рад, 0,065 рад, 0,098 рад, 0,13 рад. Полученные экспериментальные зависимости £(/) имели дополнительный уровень шума вследствие того, что щель пространственно обрезала часть падающего излучения: для координаты 0 мм соотношение сигнал/шум составляло 30:1, для координаты 30 мм это значение достигало 10:1. Методами численной обработки шум частично устранялся, при вычислениях использовался фильтр, пропускающий низкие частоты и обрезающий высокие частоты, соответствующие шуму. По полученной временной форме £(/) были рассчитаны спектры дифракции ТГц импульса. Спектральная и временная формы дифрагированного импульса в центре нулевого максимума дифракции, полученные экспериментально, приведены на рис. 3. Временная зависимость исходного импульса представляет собой одно полнопериодное колебание, за которым следуют хаотические выбросы различной полярности, связанные с дисперсией воздуха и поглощением в водных парах. Из графика следует, что амплитуда отрицательного выброса на переднем фронте импульса превышает амплитуду отрицательного выброса на заднем фронте, а спектр излучения сосредоточен в области 0,1-1,2 ТГц.
Рис. 3. Поле (а) и спектр (б) терагерцовой волны в центре нулевого максимума дифракции
Представленное на рис. 3 поле £(/) является входным для численного расчета, описанного в разделе «Математическая модель». Далее приведены зависимости для дифрагированного на щели ТГц импульса в разных характерных поперечных точках пространственной картины дифракции (рис. 4). Сплошной линией представлены экспериментальные результаты, пунктирной - расчетные.
Анализируя полученные результаты, можно отметить, что с ростом угла дифракции наблюдается сужение спектра импульса и изменение его фазовых характеристик. С увеличением угла дифракции форма исходного импульса (полнопериодного колебания) претерпевает изменения, амплитуда отрицательного выброса на заднем фронте постепенно увеличивается и в крайних точках уже превышает отрицательную амплитуду на переднем фронте, что свидетельствует об изменении фазы колебания. Ширина спектра по половинному уровню монотонно уменьшается и при угле дифракции 0,13 рад составляет ~ 0,25 ТГц. С ростом угла дифракции центральная частота импульса также смещается в низкочастотную область (с 0,75 ТГц до 0,25 ТГц), у спектра дифракции появляется вторичный максимум. Таким образом, дифракционная щель ведет себя как пространственно-временной фильтр.
При сравнении экспериментальных результатов с расчетными наблюдается хорошее согласование. Однако в точке 25 мм видно, что передний фронт расчетного импульса начинает опережать передний фронт экспериментального импульса, что еще сильнее проявляется в точке 30 мм. Это говорит о возможном искривлении исходного плоского волнового фронта ТГц импульса при приближении к краю пространственной дифракционной картины.
Рис. 4. Спектры (слева) и поля (справа) дифрагированных волн для пространственных точек : а) 15 мм; б) 25 мм; в) 30 мм. Сплошная линия соответствует экспериментальным результатам, пунктирная - расчетным результатам. © - угол дифракции в радианах; х - положение щели
в поперечных координатах
Заключение
Проведен численный расчет дифракции ТГц волны на щели. Предложена схема и собрана экспериментальная установка для проведения исследования дифракции на щели однопериодного ТГц импульса. Проведен эксперимент, зарегистрированы временные формы дифракции предельно короткого ТГц электромагнитного импульса для различных углов дифракции на щели. По экспериментальным зависимостям поля от времени получены спектры дифракции. Выявлено, что с ростом угла дифракции проис-
ходит изменение временной формы ТГц импульса, при этом центральная частота смещается в более низкочастотную область, в спектре дифракции появляется вторичный максимум, что также подтверждается численными расчетами. По результатам сравнения расчета с экспериментом сделано предположение о возможном искривлении исходного плоского волнового фронта импульса при приближении к краю пространственной дифракционной картины.
Стоит также отметить, что ранее в работе [10] теоретически изучалась дифракция Фраунгофера ТГц волны в параксиальном приближении (дифракция на оптической оси), где показано изменение временной структуры импульса и смещение спектра в высокочастотную область. В настоящей работе рассматривался другой случай - непараксиальная дифракция, т.е. дифракция под большими углами к оптической оси. Этим обусловлено отличие результатов в части смещения спектра дифракции: смещение в высокочастотную область спектра при параксиальном случае и смещение в низкочастотную область с появлением вторичного максимума при непараксиальном случае. На следующем этапе исследований планируется проведение эксперимента с целью верификации искривления исходного плоского волнового фронта сверхкороткого ТГц импульса при дифракции.
Работа выполнена в рамках работы по госконтракту № 16.513.11.3070 от 19 апреля 2011 года.
Литература
1. Glagolewa-Arkadiewa A. Short electromagnetic waves of wave-length up to 82 microns // Nature. - 1924. -V. 113. - № 2844. - Р. 640.
2. Геккер И.В., Юрьев В.И. Субмиллиметровые волны. - М.-Л. - Госэнергиздат, 1961. - 60 с.
3. Zhang X.C., Xu J. Introduction to THz wave photonics. - NY: Springer, 2009. - 246 р.
4. Lee Yun-Shik. Principles of terahertz science and technology. - Springer Science+Business Media, LLC, XII. 2009. - 340 p.
5. Беспалов В.Г., Крылов В.Н., Путилин С.Э., Стаселько Д.И. Генерация излучения в дальнем ИК диапазоне спектра при фемтосекундном оптическом возбуждении полупроводника InAs в магнитном поле // Оптика и спектроскопия. - 2002. - Т. 93. - № 1. - С. 158-162.
6. Беспалов В.Г. Сверхширокополосное импульсное излучение в терагерцовой области спектра: получение и применение // Оптический журнал. - 2006. - Т. 73. - № 11. - С. 28-37.
7. Беспалов В.Г., Ефимов Ю.Н., Милоглядов Э.В. Генерация субпикосекундного терагерцового излучения в скрещенных электрическом и магнитном полях // Научно-технический вестник СПбГУ ИТМО.
- 2004. - № 15. - С. 41-44.
8. Розанов Н.Н. О дифракции предельно коротких импульсов // Оптика и спектроскопия. - 2003. - Т. 95.
- № 2. - С. 327-330.
9. Sereda L., Ferrari A., Bertolotti M. Spectral and time evolution in diffraction from a slit of polychromatic and nonstationary plane waves // J. Opt. Soc. Am. - 1996. - V. 13. - Р. 1394-1402.
10. Езерская А. А., Иванов Д.В., Беспалов В.Г., Козлов С.А. Дифракция однопериодных терагерцовых волн с гауссовым поперечным распределением // Научно-технический вестник СПбГУ ИТМО. -2011. - № 1 (71). - С. 10-16.
11. Jingzhou U., WANG Li & YANG Guozhen. Effects of spectral linewidth of ultrashort pulses on the spatiotemporal distribution of diffraction fields // Chinese Science Bulletin. - June 2001. - V. 46. - № 11. -P. 901-904.
12. Iaconis C. and I.A. Walmsley. Spectral phase Interferometry for direct electric-field reconstruction of ultrashort optical pulses // Optics Letters. - 1998. - V. 23. - Р. 792-794.
13. Петров Н., Беспалов В., Волков М. Восстановление волнового фронта при безопорной цифровой ПЗС-регистрации мультиспектральных спекл картин // Наносистемы: физика, химия, математика. -2011. - Т. 2. - № 1. - С. 82-90.
14. Hu B.B., Zhang X.-C. and D.H. Auston. Terahertz Radiation Induced by Subbandgap Femtosecond Optical Excitation of GaAs // Phys. Rev. Lett. - 1991. - V. 67. - Р. 2709.
Куля Максим Сергеевич - Санкт-Петербургский национальный исследовательский университет
информационных технологий, механики и оптики, аспирант, [email protected]
Грачев Ярослав Владимирович - Санкт-Петербургский национальный исследовательский университет
информационных технологий, механики и оптики, аспирант, [email protected]
Городецкий Андрей Александрович - Санкт-Петербургский национальный исследовательский университет
информационных технологий, механики и оптики, кандидат физ.-мат. наук, научный сотрудник, [email protected] Беспалов Виктор Георгиевич - Санкт-Петербургский национальный исследовательский университет
информационных технологий, механики и оптики, доктор физ.-мат. наук, старший научный сотрудник, [email protected]
