Фазовая проблема и синтез оптических полей

Автор: Абрамочкин Е.Г., Волостников В.Г.

Журнал: Компьютерная оптика @computer-optics

Рубрика: Методы и элементы компьютерной оптики

Статья в выпуске: 10-11, 1992 года.

Бесплатный доступ

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

Короткий адрес: https://sciup.org/14058258

IDR: 14058258

Текст научной статьи Фазовая проблема и синтез оптических полей

•♦      1           .     1 [ dl где (rocj)^-^-(vixvy)             ay ” ay dx J *

При этом в изолированных нулях интенсивности имеет место соотношение

При распространении поля в зоне Френеля изолированные нули образуются парами с противоположными знаками (rot j)t в них. Численные эксперименты показали, что, во-первых, проблема восстановления полей по интенсивности связана с вихревыми свойствами векторного поля j, во-вторых, информация о знаке и модуле (rot j) в изолированных нулях интенсивности определяет в основном структуру поля и обеспечивает устойчивое восстановление векторного поля j.

К настоящему времени в общем случае не решен вопрос определения знаков (rot j)f в изолированных нулях. Используя аппарат теории функций комплексной переменной, удается лишь разделить нули по парам с противоположными знаками (rot j) .

Для задачи синтеза эта неопределенность означает, что могут существовать близкие (или равные) по интенсивности поля, обладающие существенно различной структурой, и пренебрежение соотношениями типа ( 1) при постановке задачи может привести к тому, что соотношение (1) все же реализуется, но при этом результат синтеза может отличаться от исходно заданного.

Например, соотношение (1) «запрещает» отдельные варианты изменения интенсивности при распространении поля: покажем, что невырожденный минимум интенсивности не может стать невырожденным нулем (см. рис. 16. Невырожденность нуля здесь понимается в смысле невырожденности для функции F(x, у), то есть F(x, у) = О, но VF(х, у) * 0). Рассмотрим область

Плоскость (х.у)

Рис. 1. Невозможный (б) и возможный (в, г) варианты реализации невырожденного нуля из локального минимума интенсивности (а) при распространении поля вдоль 1

Q = |(x, y)gR2, I(x, y) < IQ = constj-, содержащую точку минимума (семейство таких областей всегда существует). Используя соотношение (1), получим

jlVydr=lJvpdr=^ = О

5Q        SQ Q

Так как (rot j)i - вещественная непрерывная по х, у, 1 функция, то существует кривая в области Q, на которой (rot j) равен нулю и которая в •*                 1                             1

силу соотношения (rot =  —^- (V I х V ф)[ проходит через точку минимума. Пусть при распространении поля при некотором 1 интенсивность в минимуме достигает нулевого значения. В силу непрерывности (rot j)1 линия нулевого ротора должна проходить через этот глобальный минимум. С другой стороны, в невырожденном нуле интенсивности (rot j)j в соответствии с (2) имеет конечное ненулевое значение. Следовательно, это вырожденный нуль (рис. 1в). При малом возмущении он распадается на пару невырожденных нулей с конечными значениями ротора, разделенных нулевой линией (rot j) . Поэтому ротор в этих нулях имеет противоположные знаки (рис. 1г).

В качестве примера учета вихревой составляющей при синтезе световых полей рассмотрим задачу фокусировки лазерного излучения в кольцо радиуса

г на расстоянии f .

о г                 о

В параксиальном приближении известным решением, полученным по методу стационарной фазы [2], является фокусатор «сферическая линза + аксикон>:

<ро(г)

, 2.

krкгг

----+----

2ff оо

Легко видеть.

что данное решение в приближении геометрической

оптики реализует отображение

( £, Л) -> (х, у) с нулевым якобианом:

y=71+f0

Л --4-

го

Л

- г sin а

f

1 О

f

0

/ „2, 2 v С +л

0

Соответствующее дифракционное поле Фокальной плоскости имеет вид:

обладает нулевым ротором и в

FO(P)-

kEo Г ikp2

---exp ikf +---- if          ° 2f о                     о

r j Jo

kp ---- r dr,

где

Eq - амплитуда светового поля в плоскости (£,

л) ;

Р - полярный радиус;

а - радиус фокусатора.

График зависимости IQ(p) - |FQ(p)|2 от полярного радиуса р приведен на рис. 2 (п-0). Появление локального максимума интенсивности в центре обусловлено различием между геометрооптическим приближением и его волновой реализацией. Из вышесказанного ясно, что устойчивая нулевая интенсивность в центре может существовать лишь при нулевом роторе.

Простейшими решениями

фокусировки в кольцо с ненулевым ротором являются решения вида кг2 кгго

Ф (г, а)---+------+па,

"           2f f о о где п - целое число, а - arg(£+i,T)).

При

F (р, п

этом дифракционное поле Ffi(p,

в) примет вид

кЕо            ikp2

в) - ----exp ikf +---- +ine + if          ° 2f о                    о

3nin

1кгг а _____О

J е f° rj_

о

кр

г dr. (8)

Особенность при г - О не играет существенной область (г, а) € [с, а]х[О, 2л], где с<<а.

Графики зависимости 1п(р) ■ lFn(P) |2 от Р Для п ■

роли: можно взять

1, 2, 3 показаны на

рис. 2.

Рис. 2. Зависимость интенсивности от полярного радиуса для фаз, вычисляемых по формуле (7), при различных п

Если обратиться к аналогии между поведением жидкости и распространением энергии светового поля, то решения (7) и (4) соотносятся так же, как два способа создать воронку на поверхности жидкости: вращением и центрально-симметричным давлением на нее.

Фокусатор (7) реализует взаимно однозначное (при г*0) отображение с ненулевым якобианом:

nf x-r cos 9---sin е, кг nf                                      (9)

у-г sin е+----cos 6. о , кг

Однозначность отображения (9) означает, что фокусатор (7) не производит фокусировки в смысле метода стационарной фазы. Более того, положение максимума интенсивности г^ приходится на область 2,2 n f

Го+--2~2"^Т0 обусловлено тем» что решения (7) являются волновыми и не могут быть получены в рамках геометрической оптики.

Соответствующая дискретная задача при конечном числе N градаций фазы ф^(г, а) дает осцилляции интенсивности в области кольца. Например, при N=6 осцилляции не превышают 13%. На рисунках 3, 4 показаны сечения поверхностей интенсивности для шестиградационной фазы при п - 0, 1 и распределение 1} (р, 6) на окружности, где достигается максимальная интенсивность.

Решения типа (7) легко обобщаются на ряд задач фокусировки в кривую (например, эллипс). Авторами разработаны несколько подходов к решению задач фокусировки в произвольные область и кривую.

б а

Рис. 3. Сечения профилей интенсивности при шести градациях фазы в безвихревом (а) и вихревом (б) случаях

Рис. 4. Зависимость интенсивности от полярного угла на кольце максимальной интенсивности

Статья научная