CC BY
38
A.B. Кобытев, E.B. Курмпшев, H.H. Сисакян
ФОКАЛЬНЫЕ РАСПРЕДЕЛЕНИЯ ПОЛЕЙ ПРИ НАКЛОННОМ ПАДЕНИИ ИЗЛУЧЕНИЯ НА ГЕОМЕТРООПТИЧЕСКИЕ ФОКУСАТОРЫ. ФОКУСИРОВКА ПЛОСКОЙ ВОЛНЫ В КОЛЬЦО
Плоские оптические элементы (ПОЭ), получающие все более широкое распространение в различных областях физического эксперимента и прогрессивных технологий
, как правило, используются в отражательном режиме при наклонном падении излучения, что продиктовано очевидными технологическими требованиями построения оптических систем, а также необходимостью охлаждения таких элементов, предназначенных для фокусировки и коррекции излучения мощных технологических лазеров.
Наиболее развитыми как в математическом плане, так и с точки зрения приложений являются математические модели синтеза ПОЗ, базирующиеся на приближении геометрической оптики [2-5]. Однако, как показывает анализ области применимости указанных алгоритмов и практики использования ПОЭ, в определенных областях параметров фокусаторов требуется более точное знание распределений интенсивности
в фокальных областях с учетом волновой природы излучения. Задача формирования волновых полей с помощью оптических элементов для случая наклонного падения излучения имеет, как правило, более низкую симметрию в сравнении со случаем нормального падения волны. Этот факт значительно усложняет или делает невозможным аналитические расчеты (даже в зонах Френеля и Фраунгофера) , обычно использующие сведение двукратного интеграла Кирхгофа к однократным и, возможно, последующее разложение в ряды по функциям Бесселя для осесимметричного случая или в ряды по другим функциям в зависимости от симметрии задачи [6-8], Поэтому в данной работе проведено численное исследование фокальных распределений полей, формируемых геометрооптическими элементами, фокусирующими наклонно падающую плоскую волну в кольцо с равномерным распределением интенсивности вдоль кольца.
Поставленная задача решалась с помощью соответствующей программы расчета двумерного интеграла Кирхгофа
k pikr
u(P) -_*_// <p(E,tj)£ • dgdn, ...
2txí £ U)
где
k=2n/\ - волновое число;
P(x,y,z) - точка наблюдения с координатами (x,y,z);
ф(5,П) = A(g,Ti)e^"* _ комплексное поле в плоскости z=+0 с вещественной
амплитудой A(g,n) и вещественной фазой х(Е,П);
г2 = (х-£)3 + (у - т))2 + z2 и интегрирование ведется по плоской области Е (выходная плоскость ПОЭ) ,
При выполнении условия kr » 1 интеграл (1) является решением первой краевой задачи уравнения Гельмгольца Ди + к2и = 0 с граничным условием Кирхгофа [7,8]. Отметим, что используемая нами программа расчета интеграла (1), созданная на базе метода Симпсона , позволяет рассчитывать волновые поля от ПОЭ произвольной (не только осесимметричной или прямоугольной) формы в широком диапазоне волновых параметров D = Xz/тса2 10-э, где а - характерный размер апертуры [9] .
Итак, пусть ПОЭ располагается в плоскости ху (рис. 1). Плоская однородная волна падает на этот элемент под углом а к оси oz. Требуется сфокусировать такую волну в фокальное кольцо с равномерным распределением интенсивности вдоль кольца. Легко показать, что распределение фазы на выходной плоскости ПОЭ, решающее поставленную задачу в приближении геометрической оптики, дается выражением
x(E,n) = - Ali + n-sin а)2 + (Po - Л2-ros2 а + £2)2, (2)
f - фокусное расстояние ОЭ (расстояние от центра 0 до фокальной плоскости);
р0 - радиус фокального кольца;
а - угол падения.
ПОЭ имеет форму эллипса с большой полуосью а = R/cos а и малой полуосью b=R,где R - радиус круглого фокусатора для случая нормального падения излучения -b<£<b; -а < n < а .
С целью исследования описанного геометрооптического фокусатора мы провели численные расчеты волнового поля с помощью интеграла (1), используя в качестве выходной амплитудно-фазовой функции <р(Е,г|) = причем амплитуда была
постоянной А = const, а фаза определялась выражением (2). Для всех расчетов длина волны А. = 1.06 • Ю-9 м, радиус круглого фокусатора R = Ю-2 м, радиус фокаль-
y(n)
X(£j
Рис. 1. Геометрия задачи
ного кольца р0 = 2-10-3 м. Изменяемыми параметрами задачи были угол падения
а = 0°; 30°{ 45°; 75° и фокусное расстояние f = 3-Ю-2 м; 5-Ю~2 м; 10~1 м;
3-Ю"1 м; 5-Ю"1 м; 1 м.
Введем величину относительной интенсивности излучения l(f,p; а,3) =
I(f,p; a,&)/I(f,p0; а,0°) в окрестности фокального кольца х'2 + z'2 = р2, рас-
о
положенного в плоскости x'z' на расстоянии f от центра О оптического элемента. Здесь a - угол падения; Р - полярный угол, отсчитываемый от оси o'z' в плоскости x'z' (см. рисунки 1, 2). В силу симметрии задачи относительно плоскости yzCy'z'), определяемой, кроме всего прочего, выбором расчетной сетки, интенсивность I(f,p; a,3) вычислялась при значениях полярного угла 0 = 0°; 45°; 90°> 135°; 180°, а в случае необходимости и при других значениях (3. Отметим, что дополнительные контрольные расчеты для углов 3 = 225°; 270°; 315° полностью подтвердили указанную симметрию.
Рис. 2. Фокальное кольцо
Неравномерность распределения интенсивности вдоль фокального кольца описывается величиной 1(£,рв; а,в) при заданных значениях параметров £, рс, а и различных 3. В таблице приведены результаты расчетов I для различных значений
Распределение интенсивности поля вдоль фокального кольца
^Nf (м) о 3. 3 • 1 о-2 5-1 О"2 кг1 3-Ю"1 5-Ю"1 1
0° 1 0.999 1 0.999 0. 999 1 0.999 1.001 0.997 0.999 1 0.999 1 0.999 0.999 1 1 0.999 1 1 1 0.999 1 0.999 1 1 1 0.999 1 1
30° 1 0.958 1 . 0А6 1 ,И»6 1.203 1 1 . 0 2 i» 1 .078 1 . 1 <И 1 .179 1 1.012 1 . 0^7 1 .069 » .095 1 1.012 1 .038 1 . 0<»7 1 . 041» 1 1 . 006 1.017 1 . 023 1 . 025 1 1 . 004 1.013 1 .022 1 .026
45° 1 1 . 06^4 1 .2Ь2 1 . ^ 7 8 1 .622 1 1 .оз 1.129 1 .257 1.312 1 1 .021 1 . 06 1.121 1 . 1 1 1.012 1 . 06 1 . 065 1 .057 1 1.012 1 . 03^ 1 . 01(3 1 .047 1 1 . 006 1 .023 1 .038 1 . 045
75° 1 1 .088 1 . 1 2.1(32 3.576 1 1 .053 1 .2^5 1 .529 1 Л1Ь 1 0.999 1 .072 1 .159 1 . 1 69 1 1.001 1 . 099 1 .098 1 .093 1 1.017 1 .073 1 . 1 36 1.16
f и а в точках фокального кольца, определяемых значениями полярного угла 0=0° 45°; 90°; 135°; 180°.
В случае наклонного падения излучения на ПОЭ данные таблицы обнаруживают увеличение относительной интенсивности I при изменении полярного угла. Отмеченная регулярная неравномерность в распределении интенсивности вдоль фокального кольца при фиксированном а уменьшается с ростом фокусного расстояния f, а при заданном f она значительно возрастает с увеличением а. Причину наблюдаемого эффекта нетрудно понять из следующих соображений.
Рассмотрим ход лучей от фокусаторов в соответствующие точки фокального коль ца в плоскости х = 0 (см. рис. 1). В качестве наиболее контрастного примера геометрические построения приведены на рис. 3 для о. = 75° и f ■ 30, 50 мм. Для других f и а рисунки и рассуждения аналогичны.
Рис. 3. Положение фокального кольца относительно фокусатора в плоскости у=0 при а=75° и £=30, 50 мм
Согласно геометpoonтичес кому алгоритму каждая точка фокального кольца формируется соответствующим лучевым отрезком сэ, так, например, реперная точка кольца Р с полярными координатами р0 = 2 мм и ß = 0° формируется отрезком (0, а = R/cos а) оси оу, a отрезок (-a = -R/cos а, 0) отображается в точку фокального кольца Р' с координатами р0 = 2 мм и ß = 180°. Из рис. 3 и аналогичных ему рисунков для других f и а легко видеть, что при увеличении а из-за увеличения большой полуоси а = R/cos а апертуры фокусатора и в связи с приближением фокального кольца к ПОЭ происходит уширение формирующих пучков (на рис. 3 эти пучки ограничены лучами (-аР') и (оР'), (аР) и (оР)) , приводящее к нарушению условия параксиальности пучков (особенно для случаев малых f = 30, 50 мм), в рамках которого геометрооптическое решение задачи является корректным. Отмеченная непараксиальность приводит к замедлению сходимости результатов расчетов при измельчении сетки, то есть с увеличением а установление происходит на более мелких сетках.
Другим следствием непараксиальности пучков является увеличение регулярной неравномерности распределения интенсивности вдоль геометрооптического кольца при росте а. Рассматривая участки лучевых отрезков фокусатора как источники вторичных сферических волн, амплитуда которых затухает по закону А~1~1, легко оценить отношение интенсивностей, например, в точках фокального кольца Р и Р' по приближенной формуле
I(P') /KP) я [1 (Р) /1(Р')]2, (3)
где
1(Р) - расстояние от середины отрезка (о,а) до точки Р;
1(Р') -расстояние от середины отрезка (-а,о) до точки Р'. Оценки по формуле (3) для различных значений f и а не только качественны, но и по порядку величин верно отражают установленную в результате численных расчетов регулярную неравномерность распределения интенсивности вдоль кольца, объясняя тем самым ее происхождение и отличие результатов расчета по формуле (1) от гео-метрооптических расчетов фокусатора [2], предсказывающих равномерность интенсивности вдоль кольца.
Практический интерес представляет ширина фокального кольца, то есть р -распределение интенсивности при заданных значениях параметров f, а, В. В данной работе относительные интенсивности I(f,p; a,ß) = I(f,p; a,ß)/I(f,р0; a,ß = 0o) при различных фокусных расстояниях f = 30, 50, 100, 300, 500, 1000 мм и углах падения а = 0°, 30°, 45°, 75° вычислялись в фокальной плоскости x'z' вдоль оси
о1 Z 1 .
Расчет распределений интенсивности при а = 0° и f = 30, 50, 100, 300, 500, 1000 мм дает следующие результаты, которые представлены в виде характерных графиков на рис. А. Отметим, что без потери существенных особенностей поведения кривых на графиках представлено значительно меньшее число расчетных точек, отмеченных квадратиками. Ввиду общей структурной схожести графиков для фокусных расстояний f = 30, 50, 100 мм эти распределения представлены одним графиком с f = 50 мм (рис. ^а); для всех указанных фокусных расстояний четко выделяется основной пик распределения в окрестности точки р * pQ. Ширина фокального кольца, определяемая полушириной главного пика распределения Кр) в окрестности точки
Р = Р0, для случая f = 30 мм составляет гг^Х, для f = 50 мм - =г7Х, для f = 100мм я13Х.
focus -so angle - 0° wave - 0.0106
19 1С 21 22 2}
FOCI'S - t0° ANCLE - 0* HAVE - 0.0)06
• a
• *
•г* -jo -is *to
Рис. it. Радиальное распределение интенсивности поля в окрестности геометрического кольца при нормальном падении излучения: a) f = 50 мм; б) f = 300 мм; в) f «= 500 мм; г) f = Ю00 мм
Для указанных выше фокусных расстояний величины вторичных максимумов в окрестности основного максимума о = D0 не превышает 20t от величины основного пика, имея тенденцию к росту с увеличением f. Расчеты для фокусных расстояний f = 300, 500, 1000 мм проводились не только в окрестности точки о = 00, но и по всему диаметру фокального кольца. Ширина фокального кольца для f = 300 мм со-ставляет г для = 500 мм -г75Х, для f = 1 000 мм - г150А (см. рис. Ь б , в , г).
Таким образом, с увеличением f наблюдается сильное дифракционное размытие геоме трооптичес ко го кольца, которое сопровождается следующими эффектами:
- происходит значительное изменение формы главного пика (рис. Ur)¡
- максимум главного пика смещается от точки р = рд в сторону центра фокального кольца (рис. Ьв,г);
- происходит рост максимума в центре фокального кольца с последующим разделением этого пика на два относительно центральной точки х' = z' = 0 фокального кол ь ца .
Исследование зависимости формы главного пика р - распределения интенсивности фокального кольца от различных углов падения а и фокусных расстояний f проводилось следующим образом. На основе формул (1) и (2) были рассчитаны распределения относительных интенсиеностей I(f,p; a,ß=0°) для фокусных расстояний f = 30, 50, 100 мм и углов падения излучения а = 0°, 30°, 45°, 75°. В случае f = 100 мм график р - распределения I (р) для всего диапазона углов падения а = 0°, 30°, 45°, 75° практически не меняется и имеет вид, аналогичный рис. 4а, причем полуширина главного пика составляет примерно 13Х. Для f = 30, 50 мм общая картина зависимости соответствующих р - распределений I(f,p; a,ß=0°) от углов а похожа на случай с f = 100 мм - графики указанных р - распределений в окрестности точки р = р0 остаются практически неизменными и имеют вид рис. 4а с полуширинами главных пиков 7\ и 4А. соответственно. Однако следует отметить, что с уменьшением f и увеличением а вследствие нарушения параксиальности формирующих пучков сходимость результатов расчетов при увеличении числа узлов расчетной сетки происходит на значительно более мелких сетках так, что для случая f = 30 мм и а = 75° при разумном времени счета на ЭВМ типа VAX р-распределение не установилось.
С целью определения влияния граничных условий Кирхгофа [8] на полученные выше результаты для случая f = 50 мм и а = 0°, 30°, 45°, 75° были получены распределения интенсивности I(f,p; а,0) с учетом аподизации П0Э. Модуляция амплитуды отраженной волны с фазой (2) на выходной плоскости фокусатора описывалась выражением [10]
А(Б,П) = А0 exp (In + al)N), (4)
где
а = R/cos а;
b = R - большая и малая полуоси эллипса;
Tmin " минимальнь,й уровень отражения падающей плоской волны с амплитудой А0 на границе апертуры;
N - параметр "жесткости" аподизации.
Как показали расчеты, учет аподизации при т = 10_3 и N = 50 приводит к
шт
изменениям распределений относительной интенсивности не более чем на от соответствующих распределений без учета аподизации. Таким образом, в данной задаче разрывные граничные условия Кирхгофа оказывают сравнительно небольшое влияние на результаты расчетов. Следует отметить, что использование аподизации П0Э в форме, отличной от (4), может привести к выравниванию распределения интенсивности вдоль фокального кольца.
В заключение кратко суммируем результаты работы.
1. Наблюдается сравнительно быстрая сходимость результатов расчетов (до уровня Д1 = 11(N) - I (N') I/1(N) < е = 10~3) при переходе от сетки N х (N/cos а) узлов к более мелкой сетке N' х (N'/cos а) . Сходимость зависит как от f, так и от а: для f > ТОО мм, как правило, достаточно сетки с N = 256(5'2); для f £ 50 мм необходимы сетки с N > 512(700).
2. В случае нормального падения излучения в диапазоне f = 30-1000 мм наблюдается :
- значительное дифракционное размытие основного пика р - распределения интенсивности в окрестности геометрооптического фокального кольца при f > 100 мм;
- наличие и изменение пика р - распределения интенсивности в центре фокального кольца при изменении £;
- смещение максимума основного пика р - распределения к центру фокального кольца (при £ = 1000 мм,Др = -4Л).
3. Для исследованных в работе геометрооптических фокусаторов, рассчитанных под различные углы падения излучения а, с увеличением а при данном £ наблюдается рост регулярной неравномерности распределения интенсивности вдоль фокального кольца. Указанная неравномерность при фиксированном угле а убывает с ростом фокусного расстояния
Ь. В данной задаче разрывное граничное условие Кирхгофа оказывает незначительное влияние на результаты расчетов интеграла Кирхгофа.
ЛИТЕРАТУРА
1. С и с а к я н И.Н., С о й ф е р В.А. Компьютерная оптика. Достижения
и проблемы. - В сб.: Компьютерная оптика. М.: МЦНТИ, 1987, вып. 1, с. 5-19.
2. Гончарский A.B. Математические модели в задачах синтеза плоских оптических элементов. - В сб.: Компьютерная оптика. М.: МЦНТИ, 1987, вып. 1, с. 19-31 .
3. Воронцов М.А., Матвеев А.Н., Сивоконь В.П. К расчету фокусаторов лазерного излучения в дифракционном приближении. - В сб.: Компьютерная оптика. М.: МЦНТИ, 1987 , вып. 1, с. 7W8.
U. Гончарский A.B., Данилов В.А., Попов В.В. и др. Решение обратной задачи фокусировки лазерного излучения в произвольную кривую // Доклады АН СССР, 1.983 , т. 273 , tf 3, с. 605-608.
5. Гончарский A.B., Данилов В.А., Попов В.В. и др. Квантовая электроника, 1986, т. 13, с. 660.
6. К о с о б у р j Т.П., Степанов Н.С. Особенности дифракции при косом падении света на транспарант // Оптика и спектроскопия, 1986, т. 60, вып. 3, с. 588-592.
7. Б о р н М., Вольф Э. Основы оптики. М.: Наука, 1973. - 720 с.
8. Виноградова М.Б., Р у д е н к о 0.8., Сухорукое А.П. Теория волн. М.: Наука, 1 979. - 38*» с.
9- Кобытев A.B., Курмышее Е.В., С и с а к я н И.Н. Численное исследование пространственной структуры когерентных волновых полей оптического диапазона. - В сб.: Компьютерная оптика. М.: МЦНТИ, 1988, вып. 3, с. 57-70.
10. Цветков А.Д., Потапова Н.И., Седов Б.М. Влияние параметров супергауссовых аподизирующих диафрагм на распределение интенсивности излучения в дальней зоне // Квантовая электроника, 1988, т. 15, N* 5, с. 10itl-10i)5.
