УДК 535.317
П.А. Носов, И.И. Пахомов, А.Ф. Ширанков
СОСТОЯНИЕ И ПЕРСПЕКТИВЫ РАЗВИТИЯ МЕТОДОВ РАСЧЕТА ПРЕОБРАЗОВАНИЯ ЛАЗЕРНОГО ИЗЛУЧЕНИЯ ОПТИЧЕСКИМИ СИСТЕМАМИ
Рассмотрены существующие методы расчета преобразования оптическими системами лазерного излучения и определены тенденции и перспективы их развития.
E-mail: pan@bmstu.ru
Ключевые слова, лазер, лазерный пучок, оптическая система, лазерно-
оптическая система, формирование лазерного излучения.
Появление лазеров (оптических квантовых генераторов) способствовало становлению и интенсивному развитию новых научных направлений. Одним из этих направлений является лазерная оптика — раздел оптики, изучающий вопросы формирования лазерного излучения оптическими резонаторами и преобразования его оптическими системами, а также вопросы синтеза таких систем, обеспечивающих формирование лазерного пучка с требуемыми параметрами. В настоящее время это направление вполне сформировалось и продолжает развиваться, способствуя улучшению характеристик лазеров и расширению областей их применения.
Важную роль в теории лазерной оптики играют гауссовы пучки, которые формируют устойчивые резонаторы. В работах [1,2] установлено, что для гауссовых пучков при их распространении в свободном пространстве характер распределения амплитуды поля остается неизменным и лишь меняется масштаб этого распределения, определяющий поперечный размер пучка. Кроме того, на удалении z поперечного сечения гауссова пучка от перетяжки он имеет сферический волновой фронт, радиус кривизны которого определяется выражением
1 + (±42
V z
где zк — параметр конфокальности пучка — продольное расстояние от перетяжки пучка, на котором размер пучка увеличивается в \/2 раз.
Анализ выражения (1) показывает, что радиус кривизны волнового фронта гауссова пучка изменяется от 2zк до бесконечности [1,2]. Радиус кривизны волнового фронта равен бесконечности в сечении перетяжки и в дальней зоне при z ^ <х>. Поэтому формируемые устойчивыми лазерными резонаторами гауссовы пучки обладают свойствами, не присущими излучению классических источников (тепловых, люминесцентных).
Яф(г)
(1)
Пространственные параметры эрмито- и лагерро-гауссовых пучков удовлетворяют инварианту [2]. Для случая малой угловой расходимости пучка инвариант имеет вид
Jmn Ьртпвтп СОП81 М , (2)
П
где 2Нртп, 2втп — диаметр перетяжки и угловая расходимость пучка моды ТЕМтп соответственно; М2 — коэффициент распространения пучка [3]; А — длина волны лазерного излучения.
Параметр конфокальности и размер перетяжки пучка связаны сле-
Пк2
" ' "р тп
дующим соотношением гк = •
В большинстве практических случаев параметры сформированного резонатором лазерного пучка не соответствуют требуемым для решения конкретной задачи. Поэтому выходящий из резонатора пучок преобразовывают с помощью внешней (вне резонатора) оптической системы — лазерно-оптической системы (ЛОС).
В современной формулировке задачу преобразования лазерного пучка оптической системой обозначим следующим образом: по заданным параметрам и характеристикам пучка на входе в оптическую систему и ее конструктивным параметрам (г, й, п) требуется определить распределение комплексной амплитуды поля ф'(х', у', гПА) в плоскости анализа гПА и пространственные параметры преобразованного пучка с учетом совокупности реальных факторов, имеющих на практике существенное значение (дифракционные эффекты, тепловые эффекты, самофокусировка излучения и др.).
Эта задача может быть решена как строго, так и приближенно. Для получения приближенного решения вводят различные допущения. Так, рассматривая компоненты ЛОС идеальными (безаберрационными), при отсутствии дифракции получено решение в параксиальном приближении [4, 5]. Для уточнения этого решения, усложняют модель за счет учета только аберраций ЛОС или одновременно аберраций ЛОС и дифракционных эффектов. Для этого используют методы геометрической и волновой оптики.
Решение задачи преобразования лазерного излучения оптической системой с учетом дифракционных эффектов и аберраций ЛОС можно получить в рамках волновой оптики. Для этого используют скалярную и векторную теории дифракции [1,5-7]. Метод скалярной теории дифракции применяют, когда необходимо знать точную структуру электромагнитного поля лазерного пучка, преобразованного оптической системой, либо когда существенную роль играет диафрагмирование пучка, что не может быть учтено методом геометрической оптики. Векторная теория дифракции азируется на уравнениях Максвелла, т. е.
на наиболее полном описании электромагнитного поля, и используется в том случае, когда допущения остальных методов (в том числе и скалярной теории дифракции) оказываются недопустимыми.
Перейдем к рассмотрению существующих методов расчета преобразования лазерного излучения оптическими системами.
Волновые методы описания преобразования гауссова пучка идеальной оптической системой. Лазерная оптическая система считается идеальной, если она не вносит дополнительных искажений в форму волнового фронта, т. е. изменяет фазу пучка по квадратичному закону [1,5]. Так, при прохождении лазерного пучка через тонкую линзу с задним фокусным расстоянием комплексная амплитуда выходного поля имеет вид
Ф'(£> П) = > П) exP
-f^ 2 + п2)
где ф, ф' — распределения комплексной амплитуды поля на входе и выходе линзы соответственно; £, п — поперечные координаты на линзе.
Поскольку для всех мод зависимость изменения радиуса кривизны волнового фронта (1) и параметр конфокальности гк пучка одинаковы, то исследование преобразования многомодового лазерного пучка идеальной ЛОС ограничивают рассмотрением основной моды [1,2,8]. Пространственными параметрами гауссова пучка с длиной волны излучения А являются: размер перетяжки Нр (по уровню 1 /е2 от осевой интенсивности); параметр конфокальности гк; угловая расходимость 2000 (по уровню 1/е2 от осевой интенсивности). Кроме того, положение перетяжки входного пучка задается величиной гр (ее удаление относительно переднего фокуса ЛОС) или отрезком вр (ее удаление относительно первой поверхности ЛОС).
В настоящее время для определения параметров гауссова пучка на выходе идеальной ЛОС используют следующие методы: 1) метод вариансов [2,8]; 2) ABCD-метод [1,2,7]; 3) метод сопряженных плоскостей [1,2]. В методе вариансов и ABCD-методе параметры гауссова пучка в произвольном сечении описываются комплексным числом: вариансом V и комплексным параметром кривизны Когельника 1/ц_. В рассматриваемом сечении пучка комплексная величина определяется размером пучка и радиусом кривизны волнового фронта. При распространении гауссова пучка в свободном пространстве и при его преобразовании оптической системой, варианс и комплексный параметр кривизны изменяются в соответствии с соотношениями, приведенными в работах [2,7,8]. Преимуществом использования этих методов является компактность записи аналитических выражений в комплексной форме.
В работе [4] получены выражения для параметров преобразованного гауссова пучка методом сопряженных плоскостей. Этот метод основан на принципе подобия оптических полей в оптически сопряженных плоскостях [2,5]: в параксиальном приближении при отсутствии диафрагмирования амплитуды электромагнитных полей любого вида, если рассматривать их в оптически сопряженных плоскостях, подобны друг другу, причем коэффициент подобия равен линейному увеличению в. Обозначим через атп (х, у) распределение амплитуды поля моды ТЕМтп во входной плоскости Оху. Для распределения амплитуды поля в выходной плоскости О'х'у', оптически сопряженной с Оху, можно записать
т. е. поверхности уровня амплитуды (или интенсивности) являются изображением поверхностей уровня амплитуд (интенсивностей) поля в опорной плоскости, оптически с ней сопряженной. Доказательство этого принципа приведено в работе [5] при волновом рассмотрении задачи преобразования распределения поля излучения тонкой линзой.
В работах [1,2,4,7,8] установлено, что гауссов пучок после преобразования безаберрационным оптическим элементом сохраняет гауссово распределение амплитуды поля, а его параметры определяются выражениями
Здесь zK, z'K — параметры конфокальности исходного и преобразованного оптическим элементом пучка; zp, z'p — отрезки, определяющие положение перетяжки исходного и преобразованного пучка относительно переднего F и заднего F' фокусов элемента; 2hp, 2h'p и 2900, 2в'0о — диаметры перетяжки пучка и угловая расходимость до и после преобразования соответственно. При расчете продольного увеличения преобразующего элемента а знак «+» берется для линзы, а знак «—» — для сферического зеркала. Правило знаков для zp и z'p: вправо от фокусов — положительные знаки, влево — отрицательные [2].
Выражения (3) являются обобщением формулы Ньютона [8, 9] и показывают, что уже в параксиальной области принципиально отличаются законы преобразования параметров лазерного и классического (тепловых, люминесцентных и др.) излучений. Так как преобразование параметров гауссова пучка описывается формулами лазерной
(3)
оптики, то перетяжку преобразуемого пучка, имеющую плоский волновой фронт и гауссово распределение амплитуды, нельзя рассматривать как предмет для оптической системы. В классической оптике, если предмет расположен в передней фокальной плоскости, то его изображение находится в бесконечности. В лазерной оптике, если перетяжку входного пучка совместить с передней фокальной плоскостью, то перетяжка выходного пучка будет располагаться в задней фокальной плоскости.
Соотношения (3) также позволяют рассчитывать параметры гауссова пучка на выходе многокомпонентной ЛОС с реальными толщинами. В этом случае в них f' есть эквивалентное фокусное расстояние ЛОС, а отрезки zp и z'p отсчитываются от соответствующих фокусов эквивалентной ЛОС.
В то же время пространственные параметры лазерного пучка на выходе идеальной многокомпонентной ЛОС можно найти, используя хорошо известные методы вычислительной оптики [1,2,8]. Так, расходимость выходного пучка моды TEMmn
, 0' _ hmnpn
0mn |f | '
где 2hmnFn — диаметр моды пучка в передней фокальной плоскости последнего n-го компонента ЛОС; fn — его переднее фокусное расстояние. Применяя формулу линейного увеличения [8], определяем в обратном ходе плоскость, оптически сопряженную с передней фокальной плоскостью n-го компонента, и находим 2hFn. Далее, используя инвариант лазерного пучка (2), определяем размер перетяжки и параметр конфокальности выходного пучка.
Поскольку на выходе идеальной ЛОС для всех поперечных мод положение сечения перетяжки и параметр конфокальности одинаковы [1,2], то характер амплитудно-фазового распределения сохраняется. Поэтому распределение поля в произвольной плоскости анализа, отсчитываемой от сечения перетяжки выходного гауссова пучка, описывается известным выражением [1,2], в котором параметры пучка определяются по формулам (3).
Распределение поля гауссова пучка с учетом аберрационных и дифракционных искажений оптической системы. В соответствии со скалярной теорией дифракции распределение комплексной амплитуды поля ф'(х', y', znA) на выходе ЛОС можно определить используя формулу Рэлея — Зоммерфельда [7]:
ф'(х', y', ZnA) = - Щ №, п) cos 7exp Гкr'2) dE> (4)
где п) — распределение комплексной амплитуды поля на выходе оптической системы; 7 — угол между осью г и направлением от точки (£, п) к точке анализа (х', у', гПА); к = 2п/Л — волновое число; г12 — расстояние между точкой (£, п) на выходе оптической системы и точкой анализа (х', у', гПА). Интегрирование ведется по апертуре ЛОС £.
В работах [2,4] разложение волновой аберрации ЛОС от поперечных координат (£, п) на выходе ЛОС представлено в виде
Ъ2
W(£, n) = -f О + аЪ + ßh + шЪ),
(5)
где к2 = £2 + п2; /' — заднее фокусное расстояние ЛОС; а, в, ш — коэффициенты, определяющие аберрации ЛОС 3-го, 5-го и 7-го порядков.
В частном случае тонкой линзы коэффициенты разложения волновой аберрации ЛОС рассчитывают по формулам:
1 Р1 - РЬ2. 5 р\ - рт2
а = iPi-4 ß
ш
4 Р1 - р2 8 Р1 - р2 64 Р1 - р2'
где р 1,2 = 1 /г 1,2 — кривизна преломляющей поверхности линзы; г 1,2 — радиусы кривизны преломляющих поверхностей линзы.
Волновая аберрация тонкого двухлинзового склеенного объектива описывается выражением (5), в котором коэффициенты а, в, ш имеют следующий вид:
Ф
а
Ф; ß = L ß
Фп
aq = -
1 PAq - р3+1
4 pq - pq+1 '
ßq =
q. q Ф '
ш
Е
шп
*1
Ф
1p5 - р5+1 8 Pq - Pq+1 '
Шп =
7 7
5 р7 — P7+1
64 Pq - Pq+1
1,2,
где Ф = 1//' = Ф1 +Ф2 — оптическая сила тонкого объектива; Фд — оптическая сила д-го компонента объектива (Ф1 = (щ - 1) (р1 - Р2), Ф2 = (щ2 - 1) (рз - р4)); щ1,2 — показатель преломления марки стекла компонента объектива.
В фазовом множителе подынтегрального выражения дифракционного интеграла (4) расстояние между произвольными точками на линзе и в плоскости анализа представлено в виде степенного ряда с учетом четырех членов разложения. В частном случае фокусировки излучения оно имеет вид [2,4]
Г12 =
1+1(x - о2+(y1 - n)2 1( +_L( Ъ2| _
2
8
16
128
Ъ2
-А
q
4
5
z
ПА
z
z
z
ПА
ПА
ПА
ПА
Рассматривая преобразование гауссова пучка и считая в дифракционном интеграле (4) cos 7 ~ 1, r12 ~ ^ПА (в знаменателе подынтегрального выражения), получаем в результате вычисления дифракционного интеграла [2,4] распределение комплексной амплитуды поля в произвольной точке анализа после тонкой оптической системы:
eikznA
ф'(х', у', ^па) = Л—-exp
глгпА
к
2zn
х'2 + у/2
X
х exp — (а£2 + bx£)] exp — (arf + byn)] x
x exp<j -i f
^h8mSixT8(p) + hmax (ߣ + 2^ShLx)
+ hmax ( «S + 3ßZhmax + 7^S^ax ) T4(p)
d^dn, (6)
где
^ кЛ 1
2 - г~2 И + Z"
1 1 h2
a = _______hmax
= Ro f' f'
+ 1 f'
az = a + -
4 ^ПА
ПА
15
16
• кХ , . ку'
Ьх i , Ьу i ;
ZnA ZnA
7
aZ + TF ßZ hmax + 8
B
kh0'
ßz
ß 1 f' ; + 5 f'
ß - 8Uz = w + 64 ZT-
64 z7
В этих выражениях Ао и Ко — амплитуда поля и радиус кривизны волнового фронта гауссова пучка на оси на входе линзы; Л0 — размер входного пучка на линзе (по уровню 1 /е2 от осевой интенсивности); 2Л-тах — световой диаметр линзы; Т8 (р = к/ктах), Тб(р), Т4(р) — полиномы Чебышева, наименее уклоняющиеся от нуля на отрезке [ —1, +1] [2,10]; коэффициенты аЕ, вЕ и характеризуют аберрации оптической системы и свободного пространства.
При записи подынтегрального выражения (6) члены Л4, Л6, Л8 были выражены через полиномы Чебышева и сгруппированы по Т8(р), ^Мр^ Т4 (Р) •
Так как величина А определяется коэффициентами аЕ, вЕ и ше, которые, в свою очередь, зависят от гПА, то А есть функция от гПА, т. е. А = А(^па).
Пренебрегая членами с полиномами Чебышева 4-го, 6-го и 8-го порядков и расширяя пределы интегрирования до после инте-
грирования получаем распределение поля с гауссовым распределением амплитуды поля в любом поперечном сечении. Этот пучок назван аберрационным гауссовым пучком. Его распределение комплексной
а
2
амплитуды поля определяется выражением [2,4]
ф'(х', у', ¿па) =
л ho
Ä0 щА>exp
" x'2 + y'2'
[ h'2(znA)J exp
i kznA -
n x'2 + y'2 X И'ф (znA)
Выражение (7) позволяет найти зависимости и параметры аберрационного гауссова пучка [2,4]:
размер пучка в плоскости анализа по уровню интенсивности 1 /е2
h'(znA) = h^(znAÄ + 1)2 + (z„AB )2 =
X
nhoVÄ2 + B2 \
1 +
znA + Ä/ (Ä2 + B2) B/ (Ä2 + B2)
радиус кривизны волнового фронта в точке (0,0, гПА) ( ) = 4а(А2 + Б2) + 2^ + 1
КфЫ = гПА(А2 + Б2) + А ' угловая расходимость (по уровню интенсивности 1 /е2)
2^0 о = 2^ А*, + Б2,
(8)
(9)
(10)
где А
2.
R0- f
h
2
тах
« + "[бßhmax + 8^ax^ J
у ■ ~тах ■ 8
положение перетяжки выходного пучка в'р относительно последней поверхности оптической системы и размер пучка Н'р в этом сечении (по уровню интенсивности 1 /е2) соответственно
sp = -
ÄW
h'
X
1
[ÄM] + B2 P [Ä (sp)}2 + B 2
(11)
Согласно выражениям (7)-(11), распределение поля, зависимости и параметры аберрационного гауссова пучка определяются аберрациями оптической системы и свободного пространства (они вошли в величину А). Поэтому данный пучок гораздо ближе по форме к реальному пучку.
Используя формулы (10), (11), получаем выражения для инварианта и параметра М2 аберрационного гауссова пучка:
Ja6
П
\
Äl + B2 .
[Ä (s'P)}2 + B2'
M
аб
\
Äl + B2
[Ä (s'P)}2 + B2
(12)
Исследование дифракционных искажений поля лазерного пучка
2
выполнено в работе [2], где показано, что ошибка определения размера пучка зависит от отношения светового диаметра линзы и диаметра огибающей пучка на ней (2ктах/2к0). Приведенные в этой работе результаты исследований позволяют сделать следующий вывод: размер перетяжки увеличивается в 1, 32 раза при 2ктах/2к0 = 1,0; в 1,08 раза при 2ктах/2к0 = 1,5 ив 1,008 раза при 2ктах/2к0 = 2,0.
Геометрооптическое описание преобразования гауссова пучка оптической системой. При геометрическом подходе используют метод лучевых пакетов [2,11].
Совместив центр системы координат с центром перетяжки гауссова пучка, для распределения относительной интенсивности поля в произвольной точке наблюдения (х, у, г) получим
2" гк (х2 + у2)1 . (13)
Ioo (x, y, z) = exp
Л г2 + г2
В цилиндрических координатах (р, р, г) поверхность равной относительной интенсивности имеет вид
22
Р ' 1. (14)
Ро
2 z2 zк
Здесь р2 = х2 + у2, р0 = х0 + у0, причем р = \/х2 + у2 — расстояние от точки наблюдения в пучке до оси пучка, а р0 — размер пучка в сечении перетяжки (г = 0), соответствующий тому же уровню плотности мощности, что и точка наблюдения (х, у, г).
Зависимость (14) описывает однополостный гиперболоид вращения. Поэтому лазерный пучок с распределением относительной интенсивности поля (13) можно рассматривать как совокупность вложенных гиперболоидов вращения вида (14) с одинаковым параметром конфокальности гк и разными размерами пучка в сечении перетяжки р0, соответствующими различным уровням энергии (рис. 1).
Рис. 1. Линии равной относительной интенсивности
Для однополостного гиперболоида (14) выражения для направляющих векторов V1 и У2 его прямолинейных образующих имеют вид [11]
V" =Ш± у); Ш - у); 4 (>5)
Лучевым пакетом называется семейство всех лучей с направляющими векторами V! или У2, принадлежащих одному однополостно-му гиперболоиду (рис. 2). Поскольку лучевой пакет описывает часть пучка только на одном уровне плотности мощности (одном значении относительной интенсивности), то весь пучок представляется совокупностью лучевых пакетов.
Сечение перетяжки
Рис. 2. Лучевой пакет
Используя известные алгоритмы расчета хода луча через оптическую систему [8,9], определяют параметры каждого лучевого пакета на выходе ЛОС.
Несмотря на то, что метод лучевого пакета базируется на геометрической теории, и в нем для описания лазерных полей используются лучи, он достаточно строго обоснован теоретически и на практике дает адекватные с волновой теорией результаты [12-15].
Заключение. Рассмотренные в статье методы расчета преобразования лазерного излучения оптическими системами позволяют сделать следующие выводы: 1) разработаны методы расчета распределения поля и параметров лазерного пучка на выходе идеальной ЛОС; 2) для комплексной амплитуды поля на выходе реальной ЛОС в случае основной моды получено приближенное аналитическое выражение,
а для многомодового пучка разработана методика численными методами; 3) разработанные методы (аналитические и расчетные) позволяют описать искажение основной моды лазерного пучка и проводить синтез ЛОС с малыми искажениями выходного поля.
В связи с этим представляется актуальной разработка методов расчета комплексной амплитуды поля на выходе оптической системы для многомодового лазерного излучения и отдельно для высших поперечных мод. Это обусловлено их широким применением в задачах лазерной технологии при использовании мощных лазеров (волоконных, дисковых, газовых и др.). Кроме того, весьма важно анализировать лучевую стойкость элементов оптической системы и ее термоаберрации.
СПИСОК ЛИТЕРАТУРЫ
1. Пахомов И.И., Рожков О.В., Рождествин В.Н. Оптико-электронные квантовые приборы. — М.: Радио и связь, 1982. — 456 с.
2. Пахомов И.И., Цибуля А.Б. Расчет оптических систем лазерных приборов. — М.: Радио и связь, 1986. — 152 с.
3. ГОСТ Р ИСО 11146-1-2008. Лазеры и лазерные установки (системы). Методы измерений ширин, углов расходимости и коэффициентов распространения лазерных пучков. Ч. 1. Стигматические (гомоцентрические) и слабоастигматические пучки. Введен 01.01.2010. — М.: Стандартинформ, 2010. — 13 с.
4. Пахомов И.И. Расчет преобразования лазерного пучка в оптических системах. — М.: МВТУ, 1984. — 54 с.
5. Маркузе Д. Оптические волноводы. — М.: Мир, 1974. — 457 с.
6. Борн М., Вольф Э. Основы оптики. — М.: Наука, 1970. — 856 с.
7. Быков В.П., Силичев О.О. Лазерные резонаторы. — М.: Физматлит, 2004. — 320 с.
8. Прикладная оптика / М.И. Апенко и др. — М.: Машиностроение, 1992. — 480 с.
9. Русинов М.М., Грамматин А.П., Иванов П.Д. Вычислительная оптика: справочник. — М.: ЛКИ, 2008. — 423 с.
10. Пашковский С. Вычислительные применения многочленов и рядов Чебы-шева. — М.: Наука, 1983. — 384 с.
11. Аберрационный синтез оптических систем, предназначенных для преобразования лазерных пучков / П.А. Носов, В.Ю. Павлов, И.И. Пахомов, А.Ф. Ширанков // Оптический журнал. — 2011. — Т. 78; № 9. — С. 34-44.
12. Цибуля А.Б., Чертов В.Г., Шерешев А.Б. Пространственная структура лазерных пучков и геометрическая оптика // Оптико-механическая промышленность. — 1977. — № 10. — С. 66-72.
13. Быков В.П., Вайнштейн Л.А. Геометрическая оптика открытых резонаторов // ЖЭТФ. — 1964. — Т. 47; № 8. — С. 508-517.
14. Кравцов Ю.А. Об одной модификации метода геометрической оптики // Изв. вузов. Радиофизика. — 1965. — Т. 7; № 4. — С. 664-673.
15. Кравцов Ю.А. Модификация метода геометрической оптики для волны, просачивающейся через каустику // Изв. вузов. Радиофизика. — 1965. — Т. 8; № 4. — С. 659-667.
Статья поступила в редакцию 04.10.2012