Оценка дифракционного размытия фокальной линии геометрооптических фокусаторов
Автор: Голуб М.А., Казанский Н.Л., Сисакян И.Н., Сойфер В.А., Харитонов С.И.
Журнал: Компьютерная оптика @computer-optics
Рубрика: Численные методы компьютерной оптики
Статья в выпуске: 5, 1989 года.
Бесплатный доступ
Решается задача дифракционного расчета интенсивности поля вблизи фокальной линии, фокусатора когерентного излучения. Приводятся примеры расчетов волновых полей в случае фокусировки в отрезок с равномерным распределением энергии из плоского и гауссова пучка. Приводятся оценки для дифракционной ширины линии и энергетической эффективности фокусировки.
Короткий адрес: https://sciup.org/14058176
IDR: 14058176
Текст научной статьи Оценка дифракционного размытия фокальной линии геометрооптических фокусаторов
6 данной работе решается задача дифракционного расчета интенсивности поля вблизи фокальной линии фокусатора когерентного излучения. Метод основан на вычислении интеграла Кирхгофа в параксиальном приближении в системе координат, связанной со слоями на фокусаторе, уравнения которых находятся в процессе реше ния обратной задачи. Приведены примеры расчетов световых полей в случае фокуси ровки плоского гауссова пучка в отрезок с равномерным распределением интенсивности. Получены простые выражения для дифракционной ширины фокального отрезка и энергетической эффективности фокусатора.
Пусть рассматривается задача фокусировки плоского пучка когерентного излучения с длиной волны X в плоскую фокальную линию F с параметрическими уравнени ям и х=х0(Е), Ее[о, ь], (1)
где 1(E)* линейная плотность; Х=(х, у) - декартовы координаты в плоскости фокусировки; Е - натуральный параметр кривой.
Фокальная плоскость отстоит на расстояние f от плоскости фокусатора (это расстояние будем называть фокусным) .
Определим координаты в окрестности фокальной линии. Для этого введем линейное продолжение F линии F касательными в концевых точках согласно уравнениям:
х=х0(Е)+х;е(Е-Ё), ё -
где
Е» если Е€[0, L] ;
О, если Е<0;
L, если E>L,
При g£[0, L] величина I g-g I есть расстояние от точки х до соответствующей граничной точки прямой, измеренное по касательной.
Далее введем в плоскости фокусировки криволинейные координаты g, п, связанные с х, у соотношениями:
У=Уо^>+хоЕ' п; (3)
х=х0(Е)-у;£п.
Координата п» выбранная таким образом, характеризует отклонение точки х от кривой F, отсчитанное по нормали.
В геометрооптическом решении обратной задачи фокусировки [l-З] каждая точка g фокальной линии формируется соответствующим слоем, имеющим в параксиальном приближении уравнение (u-x0(g))x'(g)=c(g)f0, (4)
где функция c(g) находится из энергетических соотношений.
Результатом решения обратной задачи фокусировки является фазовая функция фокусатора кФ [2], определяемая уравнением 1 u (u-x^2 U-Uo)xc
Ф(и) = т- / Xo(g(u))dl----—г----- , го u0 го где . 2п к = у— - волновое число; u=(u, v) - декартовы координаты в плоскости фокусатора.
В дифракционном решении прямой задачи фокусировки будем считать, что поле в каждой точке в окрестности фокальной линии формируется лишь той частью излучения, которая проходит вблизи соответствующего слоя фокусатора. Отдельно для каждой точки Е линии F построим систему декартовых координат в плоскости фокуса тора так, чтобы ось В совпадала с соответствующим слоем U=(c(g)fo+xo(g)x'oS(g))x;(g)+x;g(g)xB); v=(c(g)f0+x0(g)x;E(g))y'(g)+(x;E(g)-в), (5)
где ххв=ха-у0, (х-В)=хВ+уа.
Расчет поля w(x, f0) в фокальной плоскости проводим с помощью интеграла Кирхгофа в параксиальном приближении [^]
w(x, f0)=eikf° 2nif~ Vw0(u, 0)exp(1кФ)d3u, (6)
где
Ф = ^j-(x-u)2+T. (7)
Анализируя фазу Ф, можно убедиться, что множество точек стационарной фазы образует прямую линию, совпадающую со слоем на фокусаторе.
Сделаем в (6) замену переменных (u, v)—(a, В) и разложим Ф(а, В) в ряд Тей лора, ограничившись квадратичными членами, а функцию wQ аппроксимируем ее зна чением на соответствующем слое (то есть при а=0) . Далее вычисляется интеграл по В; получаем
w(x, f0)=e1K -^— /dB(w0 (и (0, В), 0)ехр(^(хо 2пг 1 о

G (В)
G_(B)’
^П *
где
-
Y0=^о^*а=о' Yi ~ Эа^°1а=0,
4-^l , E(x> . Л/е^ах,
-
1 2 la'Ia-O n
x - фаза, не зависящая от a и g, a=G1(В) и a=G3(B) уравнения границ фокусатора, которые находятся из уравнений границ в переменных (u, v) с использованием соотношений (4) и (5).
Если в (8) 0,(0)=», Ga(B)=-“, получим асимптотическую Формулу
w(x, f0) =
^V^ ^д|
w(u(0, 0) ,O)exp(^|p)
В частности, формула (9) справедлива всюду, если фокальная кривая замкнута и гладкая, а также справедлива для окрестности внутренних точек, достаточно удаленных от концов произвольной фокальной кривой. Интенсивность поля рассчитывается по формулам (8) и (9) путем вычисления квадрата модуля.
В качестве примера рассмотрим фокусировку в фокальный отрезок прямой с равномерной линейной плотностью (рисунок). Отрезок предполагается расположенным вдоль оси и симметрично относительно начала координат, а фокусатор - прямоугольным размера 2а х 2Ь. Для него Е=х, п=У, а слои представляют прямые, параллельные оси V.

Распределение интенсивности от фокусатора гауссова пучка в отрезок (А—10,6 мк, а-10 мм, а-12,5 мм, L-15 мм)
Для фокусируемого пучка, имеющего плоский волновой фронт с равномерным рас пределением интенсивности
Ф(и, v) = - «1(1 L ) - ,
2to 23 2t0
а поле в фокальной области определяется согласно (8), выражения w (Ей) = f Ь sine (^т^) х [sign (а - хЕ ((7—) ) + nLfo nf0 3 L До
+ sign (а + ^)xE(ka_(g + |) з} .
/И"о I 1 X>0
До = v т—— ; Д = g - •= sign x = < 0 x=0
-
-1 x<0 .
Соотношение (10) позволяет представить данный фокусатор как две скрещенные линзы с фокусами f0/(l - ^-) и fQ по осям и и v соответственно. Заметим, что по ле (11), рассчитанное по формуле (8) вдоль оси g, описывается формулами дифракции Френеля сферического сходящегося пучка на щели шириной 2а, а вдоль оси п -формулами дифракции Фраунгофера на щели шириной 2Ь.
При L—0 мы получаем поле от линзы с апертурой 2а х 2Ь и фокусным расстоянием fQ.
При —д = |Д|>>Д0, то есть для внутренних точек, лежащих вдали от концов отрезка, формула (11) асимптотически переходит в
w(gn) = wosinc(^p), (12)
/я v т> описывает усредненное от оси значение комплексной ампли- туды, а sine-функция описывает дифракционное уширение До фокальной кривой.
При Д0«Ь для уровня 6 по интенсивности получим
2nf° Г
До = -k^ arsine ve, где символом arsine обозначена функция обратная к sine (х). Формула (12) позволяет оценить дифракционную эффективность фокусатора ед, то есть долю энергии, попадающую в прямоугольник с шириной Дд и длиной L eg = ~(SI (2пх) -nO^I x=arsincV§.
Для фокусируемого гауссова пучка с плоским волновым фронтом и распределением интенсивности
1=1 ехр(- —^—),
2а3
согласно (8) и (9), выражение для амплитуды имеет вид:
е 2rrf L w(gn) = ' eix Zta \/2nerf (^y)k |
) - w2 при Д»Д0; (13) x (signAE ( (y-) 3) - ^Mwn при А°Дог |


Lsxp (-где I2f Lexp(- ---) , \ / 2aa . A_ a r V koVSnerf (^5) wn = / exp(i£nv) exp(- - |
— )f 1 2a3 Д=Е - j; —)dv. |
В случае, когда размеры фокусатора больше, чем характерный размер пучка о
Wy принимает вид 2 wn=2o\HT ехр(- ——) , 452
где 3 = ^ , совпадающий с известной формулой для фокусировки гауссова пучка линзой.
При этом 2 f 1 /
Де = kS?(- ine) =;
eg = erf (v^Tn©) .