АБЕРРАЦИИ ВЫСШИХ ПОРЯДКОВ ОБЪЕМНОГО ГОЛОГРАММНОГО ОПТИЧЕСКОГО ЭЛЕМЕНТА
Юрий Цыдыпович Батомункуев
Сибирская государственная геодезическая академия, 630108, г. Новосибирск, ул. Плахотного,10, кандидат технических наук, доцент, е-mail: [email protected]
Представлены выражения коэффициентов аберраций седьмого и девятого порядков объемного голограммного оптического элемента. Показана необходимость учета аберраций девятого порядка для голограммных оптических элементов с относительным отверстием больше 1 : 1. Отмечается, что изменением температуры могут быть скорректированы отдельные типы аберраций.
Ключевые слова: голограммный оптический элемент, аберрации высших порядков, деформационные аберрации.
HIGHER-ORDER ABERRATIONS OF VOLUME HOLOGRAPHIC OPTICAL ELEMENT
Yuri Ts. Batomunkuev
Siberian State Academy of Geodesy, 10 Plakhotnogo St., Novosibirsk, 630108, Russian Federation, associate professor, е-mail: [email protected]
The formulas for coefficients of the seventh and ninth order aberrations for volume holographic optical element are presented. It is shown that aberrations of the ninth order for holographic optical element with relative aperture greater then 1:1 are necessary taken into account. It is noted that same types of aberrations can be corrected by change of temperature.
Key words: holographic optical element, higher-order aberrations, deformation aberrations.
Известно, что фокусное расстояние голограммного оптического элемента (ГОЭ) очень быстро уменьшается при увеличении рабочей длины волны. При этом возрастает величина относительного отверстия ГОЭ, что приводит к росту аберраций высших порядков. Также известно, что для светосильного тонкого ГОЭ возможно исправление сферических аберраций третьего, пятого и седьмого порядков [1-3]. Выражения коэффициентов аберраций высших порядков тонкого ГОЭ представлены в [4,5]. Однако при расчете аберраций объемного ГОЭ требуется учитывать его толщину и селективные свойства [6,7]. Целью работы является расчет аберраций высших порядков объемного ГОЭ с относительным отверстием 1:1 и более.
На рис. 1а представлена рабочая схема пропускающего объемного ГОЭ, где Pc(xc,yc,zc) - точка предмета, Pi(xi,yi,zi) - соответствующая ей точка изображения, M(x,y,z) - произвольная точка объемного ГОЭ. Ось Оz декартовой системы координат на рисунке 1 является оптической осью и направлена перпендикулярно поверхностям ГОЭ. Плоскость координат хОу проходит через центр ГОЭ. На рис. 1б представлена схема записи объемного пропускающего ГОЭ опорной и объектной сферическими волнами, где Pr(0,0,zr) - точечный
источник опорной (референтной) волны, Ро(0,0^о) - точечный источник объектной волны, М'(х',у',^) - точка объема ГОЭ при записи, соответствующая точке М(х,у^) при использовании. Обозначим расстояния в рабочей схеме как: РСМ = 1с, МРІ = 1і, а расстояния, указанные в схеме записи, как: РГМ' = 1г, М'Ро = 1о, где х'=х/тх, у'=у/ту, z'=z/mz (тх, ту, т2 - коэффициенты усадки или набухания объемного ГОЭ вдоль соответствующих осей декартовой системы координат).
Рис. 1. Рабочая схема (а) и схема записи (б) пропускающего объемного ГОЭ
Для объемного осевого ГОЭ характеристическая функция У(х,у,т) записывается в виде [4]
У(х,у^) = 1с - Тс ± (11 - т) - к^(х,у^)Хс/Пс, (1)
где т(х,у,т) = по(1г - тг ± (1о - то))/Хо - количество периодов пространственного распределения модуляции показателя преломления, коэффициента поглощения или толщины объемного ГОЭ; п и пс - начальное и среднее конечное значения показателя преломления; Хо - длина волны лазера, используемого при записи; Хс - рабочая длина волны; к - порядок дифракции.
В силу пространственной и угловой селективности объемных ГОЭ практический интерес представляет случай, когда поперечные размеры (диаметр) объемного ГОЭ намного меньше расстояний тс, z1 до плоскостей предмета и изображения, но больше размеров предмета и изображения. Кроме этого считаем, что при записи отклонение от оптической оси источников опорной и объектной волн меньше поперечных размеров объемного ГОЭ, а размеры ГОЭ много меньше расстояний тг, то от нее до источников опорной и объектной волн. Также считаем, что толщина ГОЭ много меньше расстояний тг, т0. Выражения для расстояний 1с могут быть разложены в ряд по величинам малости хс/гс, ус/гс, х/гс, у/гс, т/гс, а выражения для расстояния 1г - по величинам малости хг/гг, уг/гг, х/гг, у/гг, т/гг. Аналогичным образом, разлагая выражения для расстояний 11, 10 и подставляя в характеристическую функцию (1), получаем У(х,у,т) = Vlt(x,y) + УЦх,у^) + Уз1<х,у) + Узу(х,у,т) +
+ У51<Х,у) + V5у(х,у,т) + У^(х,у) + У7у(х,у,т) + У*(х,у) + У^(х,у,т) + ..., (2)
где Vlt(x,y) = х2р200/2 + у2р020/2 - хБш - уБою - члены разложения, характеризующие аберрации первого порядка тонкого осевого ГОЭ;
У1у(х,у,т) = (х2/2)(2р201 + 2^202) + ^/2)^021 + т2Б022) - хгБю1 - х^Б^ -утБ011 - ут Б012 - члены разложения, характеризующие объемные аберрации первого порядка;
У3,(х,у) = - (x4/8)S4оо - (х2у2/4^220 - (y4/8)S040 + (Х3/2)С300 + (х2у/2)С210 +
+ (у2х/2)С120 + (у3/2)С030 - (х2/2)А200 - хуАш - (у2/2)Л02„ -
- (х2/4)(Е200 - (у2/4)Е020 - (x/2)D100 - (у/2)Д)10 (3)
- члены разложения, характеризующие аберрации третьего порядка тонкого осевого ГОЭ;
Уз,(х,у,т) = - (x4/8)(3zS401 + 6т^402) - (x2y2/4)(3zS221 + 6т%22) -
- (y4/8)(3zSо4l + 6z2Sо42) + (х3/2)(3тСз01 + 6т2Сз02) + (х2у/2)(3тС211 + 6т2С212) +
+ (ху2/2)(3тС121 + 6т2С122) + (у3/2)(3тС031 + 6т2С(т) -
- (х2/2)(3тА201 + 6т2 А202) - ху(3тАш + 6т2Аш) - (у2/2)(3тА<ш + 6т2 А022) -
- (Х2/4)(3тЕ201 + 62^202) - (у2/4)(3тЕ021 + 672Еш) +
+ (х/2Х3г^01 + 6z2DЮ2) + (y/2)(3zDоll + 6т2БШ2) (4)
У5,(х,у), У7,(х,у), У9,(х,у) - члены разложения, характеризующие аберрации пятого, седьмого и девятого порядков тонкого осевого ГОЭ;
У5у(х,у,т), У7у(х,у,т), У9у(х,у,т) - члены разложения, характеризующие объемные аберрации пятого, седьмого и девятого порядков.
В отличие от [3] в выражении (2) характеристическая функция У(х,у,т) представлена в виде суммы членов разложения, соответствующих разным порядкам аберраций, причем в явном виде выделены типы объемных аберраций. По сравнению с характеристической функцией тонкого ГОЭ, в (2) учтены толщина, изменения толщины (усадка), изменение размеров и изменение показателя преломления объемного ГОЭ. По аналогии с известными аберрационными коэффициентами тонкого ГОЭ аберрационные коэффициенты в У3,(х,у), У3у(х,у,т) соответствуют следующим типам аберраций третьего порядка [3]: - сферической аберрации, Сук - коме, А^ - астигматизму, Еук -
кривизне поверхности изображения, Dijk - дисторсии, где индексы ij,к совпадают с соответствующими индексами аберрационных коэффициентов. Индексы ¡¿,к также совпадают соответственно с показателями степени переменных х, у, т в У3,(х,у), У3у(х,у,т). Аналитические выражения членов разложения У5,(х,у), У5у(х,у,т), У7,(х,у), У7у(х,у,т), У«и(х,у), У9у(х,у,т) и их коэффициентов аберраций приведены в приложении. Выражения для коэффициентов при индексе к = 0 совпадают с известными выражениями коэффициентов аберраций высших порядков тонкого ГОЭ [3-5].
На рис. 2 представлены расчетные значения первых сферических аберраций седьмого (пунктирная линия) и девятого (сплошная линия) порядков объемного ГОЭ с относительным отверстием примерно 1:1 на рабочей длине волны 532 нм. Радиус сферической аберрации девятого порядка превышает аберрацию седьмого порядка. Таким образом, при расчете объемных ГОЭ с относительным отверстием 1:1 и более требуется учитывать аберрации девятого порядка.
В выражениях коэффициентов аберраций в явном виде можно выделить составляющие
монохроматической и хроматической аберраций и аберраций, возникающих из-за изменения среднего показателя преломления и размеров объемного ГОЭ, в том числе
фотоиндуцированных и
термоиндуцированных аберраций. Так, средний показатель преломления nc объемного ГОЭ может быть представлен в виде [8]
nc = По + Anf + Ant + An¿, где По - начальное значение показателя преломления, Anf и Ant -фотоиндуцированное и термоиндуцированное изменения показателя преломления, Ana - изменение показателя преломления из-за деформации объемного ГОЭ.
Изменение размеров голограммного элемента вдоль каждой из осей системы координат может быть представлено в виде коэффициента [8] m = 1 + Amf + Amt + Am¿,
где Amf, Amt - относительные фотоиндуцированное и термоиндуцированное изменения размеров, Am¿ - относительное изменение размеров из-за деформации объемного ГОЭ.
Известно, что в отличие от тонкого ГОЭ при расчете аберраций объемного ГОЭ необходимо учитывать его селективность [3]. В первом приближении координаты (х,у) крайних от оптической оси точек падения лучей, для которых выполняется условие локальной селективности, в зависимости от величины коэффициентов объемных монохроматических аберраций и толщины голографической среды определяются в непараксиальной области для объемного осевого ГОЭ из формул
х = (SAA3nzS40i/8 + 3nz2S402/4))1'4, у = (§yU(3nzS04i/8 + 3nz2S042/4))1/4. (5)
При оценочных расчетах локальные величины ^х и ^у можно принять равными приблизительно п.
Из характеристической функции V(x,y,z) объемного ГОЭ могут быть определены геометрические аберрации Aх, Ay точек изображения по известным формулам
Aх ~ Zi(dV(x,y,z)/dx),
Ay ~ Z^SV^y^/Sy).
Величины аберраций Aх, Ay рассчитываются с учетом (5). При этом получаем, что дисторсия не зависит от поперечных размеров и селективности объемного ГОЭ. Отметим что, селективность объемного ГОЭ приводит к уменьшению всех аберраций по сравнению с соответствующими аберрациями тонкого ГОЭ из-за уменьшения объема ГОЭ, в котором еще выполняется
Рис. 2. Сферическая аберрация седьмого и девятого порядков
Лу,
0,3 + мм
Дх, мм
условие дифракции Брэгга. Аналогичным образом могут быть получены выражения для геометрических аберраций пятого порядка и более высших порядков.
Выражения, описывающие геометрические хроматические аберрации, то есть отклонение координат лучей в плоскости изображения на величины Дх^, Az^ от параксиальной точки изображения могут быть получены вычислением частных производных по длине волны Хс из выражений для Дх и Az
Дхх = (0(Дх)/0Хс)ДХс,
Azx = (Э^)/ЭХс)ДХс, (6)
где ДХс - изменение рабочей длины волны объемного осевого ГОЭ от заданного значения. Отметим, что при рассмотрении аберраций первого и третьего порядков объемного осевого ГОЭ от рабочей длины волны Хс в явном виде зависят только коэффициенты Fijk и Sijk. Из (6), следует, что хроматические аберрации объемного ГОЭ складываются из хроматических аберраций тонкого ГОЭ и объемной хроматической аберрации, зависящей от толщины ГОЭ. Объемная хроматическая аберрация зависит от толщины ГОЭ в первом приближении обратно пропорционально.
Выражения, описывающие геометрические деформационные аберрации Дх^ Дуd объемного осевого ГОЭ, могут быть получены из Дх и Ду вычислением частных производных по напряжению а Дхd = (Э(Дх)/Эа)Да,
Дуd = (0(ДуУЭа)Да, (7)
где Да - изменение механического напряжения в объеме ГОЭ.
Из (7) получаем следующие выражения для определения изменения координат Дх^ Ду^ Дzd луча в плоскости изображения из-за изменения на Да механического напряжения объемного осевого ГОЭ:
1) при деформационной аберрации положения
Л Л
Дzd = (дгУЭа)Да = (кц/тх )(1/zr ± 1/го)(2Дтх^тх + Дщ/n^ ,
Л 5 Л
Дzd = (дгУЭа)Да = (кц/ту )(1/zr ± 1^о)(2Дтууту + Дщ/п^,
2) при деформационной аберрации увеличения Дxd = (B100 + zBioi + z^B^^Zd,
Дyd = (B010 + zBoil + Z2Bol2)ДZd,
3) при деформационной дефокусировке
Дхd = (0F2oo/da + z0F2oi/5a + z2ЭF202/да)xZiДа + (F200 + zF20i + z^oi^z^
Дyd = (0Fo2o/da + zdFo2i/da + z2ЭFo22/да)yZiДа + (F200 + zFo2i + z2Fo22)yДzd, где 0F2OO/0a = k^(1/zr ± 1/zо)(2Дmx,d + Дп^п^/Да;
0F2O1/0a = k^(1/zr2 ± 1/zо2)(2дmx,d + Дmz,d + Д^^уДа;
0F2O2/0a = k^(1/zr 3 ± 1/zо3)(2Дmx,d + 2ДmZ,d + Д^^уДа;
0FO2O/0a = k^(1/zr ± 1^о)(2Дту^ + Д^^уДа;
0FO21/0a = k^(l/zr2 ± 1^о2)(2Дту^ + Дmz,d + Дщ/^УДа;
5FO22/0a = k^(1/zr 3 ± 1^о3)(2Дту4 + 2ДmZ,d + Д^^уДа;
4) при деформационной сферической аберрации третьего порядка Дхd = (^400/5а + 3zЭS401/Эа + 6z2ЭS4o2/Эа)x3ziДа/2 +
+ (S400 + 3zS4oi + 6z2S402)x3Дzd/2 +
+ (^220/Эа + 3z0S221/0а + 6z2ЭS222/Эа)xy2ziДа/2 + (S220 + 3zS221 +
+ 6z2S222)xy2Azd/2,
Ayd = (0SO4O/ög + 3zöSO41/0G + 6z20SO42/0Q)y3ziAa/2 +
+ (S400 + 3zSo4i + 6z2So42)y3Azd/2 +
+ (öS22o/öa + 3z0S22i/ÖG + 6z20S222/ÖG)yx2ziAö/2 + (S220 + 3zS22i +
+ 6z2S222)yx2Azd/2,
-5 -5
где 0S4OO/ög = k^(1/zr ± 1/z0 )(4Amxd + And/no)/Aa; öS401/ög = k^(1/zr4 ± 1/z04)(4Amx,d + Amzd + And/no)/Aa; öS402/ög = k^(l/zr5 ± 1/z05)(4Amx,d + 2Amzd + And/no)/Aa; öS040/ög = k^(l/zr3 ± 1/z03)(4Amy=d + And/no)/Aa; öS041/ög = k^(l/zr4 ± 1/z04)(4Amy;d + Amzd + And/no)/Aa; öS042/ög = k^(l/zr5 ± 1/z05){4Amy;d + 2Amzd + And/no)/Aa;
5) при деформационной коме третьего порядка
Axd = 3(C300 + 3zC301 + 6z2C302) x2Azd/2 + (C210 + 3zC211 + 6z2C212)xyAzd +
+ (C120 + 3zCi2i + 6z2Ci22)y2Azd/2,
Ayd = (C210 + 3zC2ii + 6z2C2i2)x2Azd/2 + (C120 + 3zCm + 6z2Ci22)xyAzd +
+ 3(Co3o + 3zCo3i + 6z2Co32)y2 Azd/2;
6) при деформационном астигматизме третьего порядка
л л
Axd = (A200 + 3zA20i + 6z A202)xAzd + (АПо + 3zAm + 6z Am)yAzd,
Л Л
Ayd = (A110 + 3zAin + 6z Aii2)xAzd + (A020 + 3zAo2i + 6z Ao22)yAzd;
7) при деформационной кривизне поверxности изображения третьего порядка
Axd = (E200 + 3zE20i + 6z2E202)xAzd/2,
Ayd = (E200 + 3zE20i + 6z2E202)yAzd/2;
8) при деформационной дисторсии третьего порядка
Axd = (D100 + 3zD101 + 6z2D102)Azd/2,
Ayd = (D100 + 3zD101 + 6z2D102)Azd/2.
Отметим, что величина дисторсии не зависит от размеров объемного осевого ГОЭ и определяется величиной деформационной аберрации положения Azd.
Важность полyчения аналитическиx выражений деформационной аберрации объемного ГОЭ заключается в том, что они позволяют получить количественные величины этиx аберраций, сравнить иx с заданными допустимыми погрешностями. Аналитические выражения для расчета термоиндуцированные и фотоиндуцированные аберраций могут быть получены из выражений для деформационные аберраций заменой в ниx деформационные изменений показателя преломления и коэффициента усадки на соответствующие фотоиндуцированные и термоиндуцированные изменения.
Полученные в явном виде выражения коэффициентов аберрации позволяют выявить новые возможности компенсации фотоиндуцированные, термоиндуцированные и деформационные аберраций и установить более оптимальные условия использования объемные ГОЭ. Так, например, фотоиндуцированная, термоиндуцированная и деформационная сферические аберрации взаимно компенсируются при изменении температуры на величину AT = [±(1/zC+j+k-1 - 1/zi j-1)/(^Ao + A^AoXl/zj-1 - 1^+к-1) + 1 -i(Amx,f + Amx,d) - j(Amy,f +Amy,d) - k(Amz,f + Amz,d) - An^/n -
-Ancd/no]/(iax + + jay + kaz + Эпс/ПдЭТ + Эп^/п^Т), (8)
где ax, ay, az - коэффициенты температурного расширения вдоль осей Ox, Oy, Oz; ЭпсупоЭТ - фотоиндуцированный термооптический коэффициент объемного осевого ГОЭ; 0пс^/по0Т - деформационный термооптический коэффициент объемного осевого ГОЭ; Amxf, Amxd Amyf, Amyd, Amzf, Amzd -фотоиндуцированное и деформационное изменения коэффициента усадки относительно осей Ох, Оу, Oz; An^f, Anc,d - фотоиндуцированное и деформационное изменения показателя преломления.
Формула (8) может быть использована и локально при неоднородном изменении показателя преломления и усадки ГОЭ. В этом случае требуется формирование соответствующего пространственного поля температур AT(x,y,z) в объеме ГОЭ.
Таким образом, в работе получены выражения коэффициентов объемных аберраций седьмого и девятого порядков объемного голограммного оптического элемента. Показана необходимость учета аберраций девятого порядка для голограммных оптических элементов с относительным отверстием больше 1:1. В явном виде выделены деформационные аберрации объемного голограммного оптического элемента. Получено, выражение, для вычисления температуры, при которой могут быть скорректированы отдельные типы аберраций.
БИБЛИОГРАФИЧЕСКИЙ СПИСОК
1. Meier R.W. Magnification and third-order aberration in holography // JOSA.1965.V. 55. P. 987-992.
2. Latta J.N. Fifth-Order Hologram Aberrations // Appl. Optics. - 1971. - V. 10, N 3. - P. 666-667.
3. Ган М.А. Теория и методы расчета голограммных и киноформных оптических элементов. Л., ГОИ. 1984. 140 с.
4. Mehta P.S., Rao K., Syam Sunder, Hradaynath R. Higher order aberrations in
holographic lenses / Appl. Opt.1982.V. 21, № 24. P. 4553-4558.
5. Rebordao J.M. General form for aberration coefficients in holography // J. Opt. Soc. Am. A. 1984. V.1. N. 7. P. 788-790.
6. Forshaw M.R.B. The imaging properties and aberrations of thick transmission holograms // Opt. Acta.1973.V. 20, № 9. P. 669-686.
7. Батомункуев Ю.Ц., Сандер Е.А., Шойдин С.А. Аберрации объемных голограмм // Тезисы Всесоюзного семинара ’’Автоматизация проектирования оптических систем“(26 - 28 октября 1988 г. Москва).- М., 1989.-C.101-112.
8. Батомункуев Ю.Ц. Особенности аберраций рентгеновского голограммного элемента / Сб. трудов VI Международной конференции «ГОЛОГРАФИЯ ЭКСПО-2009». -Киев, 2009. - С.151-154.
Приложение
Члены разложения пятого порядка характеристической функции объемного ГОЭ V5t(x,у) = (S600X6 + 3S420X4y2 + 3S240X2y4 + So60y6)/16 + 3(S400X4 + 2 S220X2y2 + So40y4)/16 -
- 3(C500X5 + 2C320X3 y2 + C140Xy4 + C410X4y + 2C230X2 y3 + C050y5)/8 -
- 3(C300X3 + C210X2y + C120Xy2 + C030y3)/4 + 3^200X2 + 2Аlloxy + А020У2)/4 +
+ 3(E200X2 + E020y2)/16 - 3(D100X + D010y)/8 +
+ 3(P 400X4 + P220X2y2 + 2P310X3y + 2PB0Xy3 + P040y4)/4 -(Z300X3 + 3Z210X2y + 3Z120Xy2 + Z030y3)/2,
V5v (x,y,z) = 5z(S60ix6+ 3S42ix4y2 + 3S24ix2y4 + So6iy6)/16 +
+ 15z(S40ix4 + 2S22ix2y2 + So4iy4)/16 -
- i5z(C50iX5 + 2C32iX3 y2 + Ci4ixy4 + C4iix4y + 2C23iX2 y3 + Co5iy5)/8 -
- i5z(C30ix3 + C2iix2y + Ci2ixy2 + Co3iy3)/4 +
+ i5z(Ä20ix2 + 2Aiiixy + Ä02iy2)/4 + i5z(E20ix2 + E20iy2)/i6 -
- i5z(Di0ix + D0iiy)/8 +
+ i5z(P40ix4 + P22ix2y2 + 2P3iix3y + 2Pi3ixy3 + P04iy4)/4 +
5z(Z30ix3 + 3Z2iix2y) + 3Zi2ixy2 + Z03iy3)/2,
где коэффициенты аберраций пятого порядка объемного осевого ГОЭ:
Sijk = 1/z/+k ± 1/zi5+k - (k^/mximyjmzk)(i/zr5+k ± 1/zo5+k ) (ijk = 600; 60i; 060; 061; 420; 42i, 240, 24i),
Sijk = (xc2 + yc2)/zc5+k ± (x2 + yi2)/zi5+k (ijk = 400; 401; 040; 041; 220; 221),
Cjk = xc(xc2 + yc2)(5j2/zc5+k ± xi(xi2 + yi2)(5j2/zi5+k (ijk = 500; 501; 320; 321; 140; 141; 300;
301; 210; 211),
Cjk = yc(xc2 + yc2)(5 j2/zc5+k ± yi(xi2 + yi2)(5-l-j)/2lzi5+k (ijk = 050; 051; 230; 231; 410; 411; 030;
031; 120; 121),
Ajk = xciycj(xc2 + yc2)/zc5+k ± Xiiyj(xi2 + yi2)/zi5+k (ijk = 200; 201; 020; 021; 110; 111),
Ejk = (xc2 + yc2)2/zc5+k ± (Xi2 + yi2)2/zi5+k (ijk = 200; 201; 020; 021),
Djk = xciycj(xc2 + yc2)2/zc5+k ± Xiiy/(xi2 + yi2)2/zi5+k (ijk = 100; 101; 010; 011),
Pjk = xci/2ycj/2/zc5+k ± Xii/2y//2/zi5+k (ijk = 400; 401; 040; 041),
Pjk = xcyc/zc5+k ± Xiyi/zi5+k (ijk = 310; 311; 130; 131),
Pjk = (xc2 + yc2)/zc +k ± (Xi2 + yi2)/zi5+k (ijk = 220; 221),
Zjjk = xciycj(xc2 + yc2)4"i"j/zc5+k ± Xiiyij(xi2 + yi2)4-jzi5+k (ijk = 300; 301; 210; 211; 120; 121; 030;
031).
Члены разложения седьмого порядка характеристической функции объемного ГОЭ V7t(x,y) = - 5(S800 x8+ 4S620x6y2 + 6S440x4y4 + 4S260x2y6 + S080y8)/128 -
- 5(S600x6 + 3S420x4y2 + 3S240x2y4 + S060y6)/32 - 15(S400x4 + 2 S220x2y2 + S040y4)/64 -
+ 5(C700x7 + 3C520x5y2 + 3C340x3y4 + Ci60xy6 + C610x6y + 3C430x4y3 + 3C250x2y5 + C070y7)/16 +
+ 15(C500x5 + C410x4y + 2C320x3y2 + 2C230x2y3 + Ci40xy4 + C050y5)/16 +
+ 15(C300x3 + C210x2y + Ci20xy2 + C030y3)/16 -
- 15(Ä200x2 + 2Aii0xy + Ä020y2)/16 - 5(E200x2 + E020y2)/32 + 5(Di00x + D0i0y)/16 -
- 15(P600x6 + 2P5i0x5y + P420x4y2 + 2P 420x4y2 + 4P330x3y3 + 2P240x2y4 + P240x2y4 + 2Pi50xy5 +
+ P060y6) /16 + 5(P500x5 + 3P4i0x4 y + P320x y2 + P230x2y3 + 3Pi40xy4 + P050y5)/4 -
- 15(P400x4 + 2P3i0x3y + P220x2y2 + 2Pi30xy3 + P040y4)/8 -
- 5(Z400x4 + 4Z3i0x3 y + 6Z220x y2 + 4Zi30xy3 + Z040y4)/8 +
+ 5(Z300x3 + 3Z210x2 y + 3Zi20xy2 + Z030y3)/4,
V7v(x,y,z) = - 35z(S801 x8 + 4S62ix6y2 + 6S44ix4y4 + 4S26ix2y6 + S08iy8)/128 -
- 35z(S60ix6 + 3S42ix4y2 + 3S24ix2y4 + S06iy6)/32 - 105z(S40ix4 + 2S22ix2y2 + S04iy4)/64 -
+ 35z(C70ix7 + 3C52ix5y2 + 3C34ix3y4 + Ci6ixy6 + C61ix6y + 3C43ix4y3 + 3C25ix2y5 + C07iy7)/16 +
+ 105z(C50ix5 + C4iix4y + 2C32ix3y2 + 2C23ix2y3 + Ci4ixy4 + C05iy5)/16 +
+ 105z(C30ix3 + C2iix2y + Ci2ixy2 + C03iy3)/16 -
- 105z(Ä20ix2 + 2Aiiixy + Ä02iy2)/16 - 35z(E20ix2 + E02iy2)/32 + 35z(Dmx + D0iiy)/16 -
- 105z(P60ix6 + 2P5iix5y + P 42ix4y2 + 2P 42ix4y2 + 4P33ix3y3 + 2P 24ix2y4 + P 24ix2y4 + 2Pi5ixy5 +
+ P06iy6) /16 + 35z(P50ix5 + 3P4iix4 y + P32ix y2 + P23ix2y3 + 3Pi4ixy4 + P05iy5)/4 -
- 105z(P40ix4 + 2P3iix3y + P 22ix2y2 + 2Pi3ixy3 + P04iy4)/8 -
- 35z(Z40ix4 + 4Z3iix3 y + 6Z22ix y2 + 4Zi3ixy3 + Z04iy4)/8 +
+ 35z(Z30ix3 + 3Z2iix2 y + 3Zi2ixy2 + Z03iy3)/4,
где коэффициенты аберраций седьмого порядка объемного осевого ГОЭ:
8ук = ШС+к ± 1/^7+к - (кц/тхгшу7т/)(1/2г7+к ± 1^7+к ) (ук = 800; 801; 080; 081; 620; 621, 260,
261, 440, 441),
с / 2 I 2\ / 1+к _|_ / 2 I 2\/ 7+к
Ьгук = (Хс + Ус )/^о ± + Уi )/Zi
241),
2 2 2 7+к 2 2 2 7+к
8ук = (Хс + Ус ) /Zc ± (Хi + у ) /^
Сук = Хс(хс2 + yc2)(7-г-у)/2/zc7+k ± ^(х2 + у/2)(7-2у%7+к
161; 500; 501; 320; 321; 140; 141; 300; 301; 120; 121),
Сук = Ус(Хс2 + Ус2)(7-^с7+к ± Уi(Хi2 + у/*™
611; 050; 051; 230; 231; 410; 411; 030; 031; 210; 211),
л 2 У í 2 | 2\2 / 1+к I / у / 2 I 2\2 / 1+к
Лук = Хс Ус (Хс + Ус ) /Zc ± Xi у/(х/ + у/ ) /Zi
т-. / 2 , 2ч3 / 7+к , / 2 , 2Ч3, 7+к
Еук = (Хс + Ус ) ^с ± (Xi + Уi ) /^
тч 2 /V 2 | 2\3 / 7+к I 2 // 2 I 2\3 / 7+к
Дук = Хс Ус (Хс + Ус ) ^с ± Xi у/(х/ + yi ) /Zi
^ 2 У/ 2 , 2\ 4-2-2 / 7+к I 2 2/ 2 , 2\ 4-2-2 / 7+к
Яук = Хс УсУ(Хс + Ус ) у^с ± Xi Уi/(xi + Уi ) ^
130;131; 040; 041; 300; 301; 210; 211; 120; 121; 030; 031),
тл 2-4 2 / 7+к | 2-4 2 / 7+к
Ргук = Хс Ус^с ± Xi У//^
т« 2 У-4 / 7+к , 2 У-4 / 7+к
Рук = ХсУс2 /Zc ± XiУi/ /Zi
т> / 7+к , , 7+к
Рук = ХсУс^с ± XiУi/Zi
2 2 7+к 2 2 7+к
Р/7к = (2Хс + Ус )/Zc ± (2х/ + у/
2 2 7+к 2 2 7+к
Ргук = (Хс + 2ус )/Zc ± (Х/ + 3у/ )/Zi
тл 2-3 2 / 7+к I 2-3 2 / 7+к
Р ук = Хс Ус^с ± X/ у//^
тл 2 1~3 I 1+к I 2 )-3 / 1+к
Рук = ХсУс2 /Zc ± X/ У!2 /^
Р 27'к = (3Хсус2 + Хc3)/Zc7+k ± (3Xiyi2 + Х/3)/^7+к
Рг7к = (3Хс Ус + Ус3)^с7+к ± (3Xi у/ + У/3)/^7+к тл 2-2 2/ 2 I 2\ / 7+к I 2-2 2/ 2 I 2\ / 7+к
Р^к = Хс Ус2(Хс + Ус )^с ± X/ у/(х/ + у/ )/Zi
тл 2 2 I 2\ / 7+к I 2 у-2/ 2 I 2\ / 7+к
Рук = Хс Ус2 (Хс + Ус)^с ± X/ у/ (Х/ + у/ )/Zi
2 2 2 7+к 2 2 2 7+к
Р^к = (Хс + Ус ) /Zc ± (X/ + у/ ) /Zi
Члены разложения девятого порядка характеристической функции объемного ГОЭ У*(х,у) = + 1(Б1000Х10+ 5Б820Х8у + 10 8640Х6у4 + 108460Х4у6 + 58280Х2у8 + Б0100у10)/256 + 35(8800 х8 + 48620Х6у2 + 6 8440Х4у4 + 48260Х2у6 + 8080у8)/256 +
+ 35 (8600Х6 + 3 8420Х4у2 + 38240Х2у4 + 8060у6)/128 + 35(8400 х4 + 28220 Х2у2 +
+ 8040у4)/128 - 35(С900Х9 + 4С720Х7у2 + 6С540Х5у4 + 4С360Х3у6 + С180ху8 + С810Х8у +
+ 4С630Х6у3 + 6С450Х4у5 + 4С270Х2у7 + С090у9)/128 - 35(С100Х7 + С610Х у + 3С520Х5у2 + 3С430Х4у3 +
+3Сэ40Х3у4 + 3С250Х2у5 + С160Ху6 + С070у7)/32 - 105(С500Х5 + С410Х4у +
+ 2С320Х у2 + 2С230Х у3 + С140ХУ4 + С050У )/64 - 35(С300Х3 + С2юх2у + Сшху2 +
+ С030у3)/32 + 35(Л200Х2 + 2Лшху + Л020у2)/32 + 35(Е200Х2 + Е020у2)/256 -
- 35(ДШ0Х + Д010у)/128 + 35(Р800Х8 + 2Рц0Х7у + Р620Х6у2 + 6Р530Х у3 + 3Р440Х4у4 +
+ 6Р 350Х3у5 + Р260Х2у6 + 2Р170Ху7 + Р080у8)/32 + 105(Р600Х6 + 2Р510Х5у + Р420Х4у2 +
+ 4Р330Х3у3 + Р240Х2у4 + 2Р150ХУ5 + Р060у6)/32 + 35(0600Х6 + 40510Х5у + 0420Х4у2 +
+ 4QззoХ3y3 + 0240Х2у4 + 4Ql50Хy5 + 0060у6)/16 + 105(Р400Х4 + 2Р310Х3у + Р220Х2у2 +
+ 2РШХу + Р040у4)/32 - 35(Р100Х7 + 3Р610Х6у + Р520Х у2 + Р430Х4у3 + Р340Х3у4 +
+ Р250Х2у5 + 3Р160Ху6 + Р070у7)/16 - 35(Р500Х5 + 3Р410Х4у + Р320Х3у2 + Р230Х2у3 +
+ 3Р 140ху4 + Р050у5)/8 - 7(2500Х5 + 5Ь410Х4 у + 102320 Х3у2 + 102230Х2у3 + 52шху4 +
+ 2050у5)/8 + 35(2400Х4 + 423юх3у + 62220Х2у2 + 42шху3 + 2040у4)/16 - 35(2300Х3 +
+ 32210Х2у + 32шху2 + 2030у3)/16,
У9v(x,y,z) = 63(81001Х10+ 5 8 821Х8у2 + 10 8641Х6у4 + 108461Х4у6 + 58281Х2у8 +
+ 80101у10^/256 + 35(8801 х8 + 48621Х6у2 + 6 8441Х4у4 + 48261Х2у6 + 8081у8^/256 +
+ 35 (8601Х6 + 38421Х4у2 + 3 8241Х2у4 + 8061у6^/128 + 35(8401 х + 28221 х2у2 +
+ 8o4ly4)z/128 - 35(С901Х9 + 4С721Х7у2 + 6С541Х5у4 + 4С361Х3у6 + С181ху8 + С811Х8у +
+ 4С631Х6у3 + 6С451Х4у5 + 4С271Х2у7 + С091у9^/128 - 35(С701Х7 + С611Х6у +3С521Х5у2 +
+ 3С431Х4у3 + 3С341Х3у4 + 3С251Х2у5 + С161ХУ6 + С071у7>/32 - 105(С501Х5 + С411Х4у +
(ук = 600; 601);
(ук = 060; 061);
(ук = 510; 511; 330; 331; 150; 151); (ук = 420; 421);
(ук = 240; 241);
(ук = 500; 501; 410; 411);
(ук = 050; 051; 140; 141);
(ук = 320; 321);
(ук = 230; 231);
(ук = 400; 401; 310; 311);
(ук = 040; 041; 130; 131);
(ук = 420; 421; 240; 241).
(ук = 600; 601; 060; 061; 420; 421, 240,
(ук = 400; 401; 040; 041; 220; 221),
(ук = 700; 701; 520; 521; 340; 341; 160;
(ук = 070; 071; 250; 251; 430; 431; 610;
(ук = 200; 201; 020; 021; 110; 111),
(ук = 200; 201; 020; 021),
(ук = 100; 101; 010; 011),
(ук = 400; 401; 310; 311; 220; 221;
+ 2Сз21Х3у2 + 2С231Х2у3 + Ci4ixy4 + Co5iy5)z/64 - 35(Сз01Х3 + С211Х2у + С121ХУ2 +
+ Co31y3)z/32 + 35(Ä201X2 + 2Лшху + Ao21y2)z/32 + 35(E201X2 + Eo21y2)z/256 -- 35(D101X + Dony)z/128 + 35(Pgo1X8 + 2P7Ux7y + P621X6y2 + 6P531X5y3 + 3P441X4y4 + + 6P 351X3y5 + P 261X2y6 + 2P 171Xy7 + Po81y8)z/32 + 105(P601X6 + 2P511X5y + P 421X4y2 + + 4P331X3y3 + P 241X2y4 + 2P151Xy5 + Po61y6)z/32 + 35(Q601X6 + 4Q5UX5y + Q421X4y2 +
+
4Q331X3y3 + Q241X2y4 + 4Q151Xy5 + Qo61y6)z/16 + 105(P401X4 + 2P311X3y + P221X2y + 2P131Xy + Po41y4)z/32 - 35(P701X7 + 3P6nX6y + P521X y2 + P431X4y3 + P341X3y4 +
+ P 251X2y5 + 3P161Xy6 + Po71y7)z/16 - 35(P501X5 + 3P411X4y + P321X3y2 + P231X2y3 +
+ 3P141Xy4 + P051y5)z/8 - 7(Z501X5 + 5L411x4 y + 10Z321 X3y2 + 10Z231X2y3 + 5Z141Xy4 + + Zo51y5)z/8 + 35(Z401X4 + 4Z3nX3y + 6Z221X2y2 + 4Z131Xy3 + Zo41y4)z/16 - 35(Z301X3 + + 3Z211X2y + 3Z121Xy2 + Zo31y3)z/16,
24
22
+
где коэффициенты аберраций девятого порядка объемного осевого ГОЭ:
Sy* = 1/zc9+k ± 1/z9+k - (k^/mX'my7mzk)(1/zr9+k ± 1/zo9+k ) (ijk = 1000; 1001; 0100; 0101; 820; 821, 280, 281, 640, 641, 460, 461),
Sijk = (Xc2 + yc2)(10'ij/2/zc9+k ± (Xi2 + y2)(10jv+k (ijk = 800; 801; 080; 081; 620; 621, 260,
261, 440, 441),
Cjk = Xc(Xc2 + yc2)(9 j2/zc9+k ± Xi(Xi2 + yi2)(9-ij2/zi9+k (ijk = 900; 901; 720; 721; 540; 541; 360;
361; 180; 181; 700; 701; 520; 521; 340; 341; 160; 161; 500; 501; 320; 321; 140; 141; 300; 301; 120; 121),
Pi-6 i / 9+k , i-6 i / '
ijk = Xc yc/zc ± Xi yiVzi
tj i j-61 9+k 1 i j-61 9+k
Pijk = Xc yc /zc ± Xiy/ /zi
Pijk = Xcyc/zc9+k ± Xiyi/zi9+k 171);
2 2 9+k 2 2 9+k
Pijk = (3Xc + yc )/zc ± (3xi + yi J/zi
2 2 9+k 2 2 9+k
Pijk = (Xc + 3yc )/zc ± (Xi + 3yi )/zi
t> / 2 I 2\ / 9+^i/ 2 _|_ 2\ / 9+k
Pijk = (Xc + yc )/zc ± (Xi + yi )/zi
tj i-4 j / 9+k , i-4 j / 9+k
Pijk = Xc Ус/zc ± Xi y//zi
i j-4 9+k i j-4 9+k
Pijk = Xcycj /zc ± Xiy/ /zi
Pijk = (3Xcyc2 + 2Xc3)/zc9+k ± (3xiyi2 + 2Xi3)/zi9+k
Pijk = (3Xc2yc + 2yc3)/zc9+k ± (3xi2yi + 2yi3)/zi9+k
Pijk = (6Xc2yc + yc3)/zc9+k ± (6Xi2yi + yi3)/zi9+k
Pijk = (6Xcyc2 + Xc3)/zc9+k ± (6Xiyi2 + Xi3)/zi9+k
Pijk = Xc"4ycj(Xc2 + yc2)/zc9+k ± Xii-V(Xi2 + yi2)/zi'
Pijk = Xciycj-4(Xc2 + yc2)/zc9+k ± Xiiyij-4(Xi2 + yi2)/zi'
9+k
9+k
Cj = yc(Xc2 + yc2f-1-j)/2/zc9+k ± yi(Xi2 + yi2)(9-jV+k (ijk = 090; 091; 270; 271; 450; 451; 630;
631; 810;811; 070; 071; 250; 251; 430; 431; 610; 611; 050; 051; 230; 231; 410; 411; 030; 031; 210; 211),
Aijk = Xciycj(Xc2 + yc2)3/zc9+k ± Xiiy/(Xi2 + yi2)3/zi9+k (ijk = 200; 201; 020; 021; 110; 111),
Ejk = (Xc2 + yc2)3/zc9+k ± (Xi2 + yi2)3/z9+k (ijk = 200; 201; 020; 021),
Dijk = Xciycj (Xc2 + yc2)4/zc9+k ± Xiiy/(Xi2 + yi2)4/zi9+k (ijk = 100; 101; 010; 011),
Zjjk = XcVcW + yc2)5-i-j/zc9+k ± Xiiyij(Xi2 + yi2)5-i-j/zi9+k (ijk = 500; 501; 410; 411; 320; 321; 230;
231; 140; 141; 050; 051; 400; 401; 310; 311; 220; 221; 130; 131; 040; 041; 300; 301; 210; 211; 120; 121; 030; 031),
(ijk = 800; 801);
(ijk = 080; 081);
(ijk = 710; 711; 530; 531; 350; 351; 170;
(ijk = 620; 621);
(ijk = 260; 261);
(ijk = 440; 441);
(ijk = 700; 701; 610; 611);
(ijk = 070; 071; 160; 161);
(ijk = 520; 521);
(ijk = 250; 251);
(ijk = 430; 431);
(ijk = 340; 341);
(ijk = 600; 601; 510; 511);
(ijk = 060; 061; 150; 151);
(ijk = 420; 421);
(ijk = 240; 241);
9+k
9+k
c yc (Xc + yc )/zc
Pijk = (2Xc2 + yc2)(Xc2 + yc2)/zc9+k ± (2x2 + yi2)(Xi2 + yi2)/zi9+k Pijk = (Xc2 + 2yc2)(Xc2 + yc2)/zc9+k ± (Xi2 + 2yi2)(Xi2 + yi2)/z9+k Pijk = Q ijk = Xcyc(Xc2 + yc2)/zc9+k ± Xiyi(Xi2 + yi2)/zi9+k (ijk = 330; 331);
Qijk = Xci"2ycj/zc9+k ± Xii-2yij/zi9+k (ijk = 600; 601; 510; 511);
Qijk = Xciycj-2/zc9+k ± Xiiy/-2/zi9+k (ijk = 060; 061; 150; 151);
Qijk = (6Xc2yc2 + Xc4)/zc9+k ± (6xi2yi2 + Xi4)/zi9+k (ijk = 420; 421);
Qijk = (6Xc2yc2 + yc4)/zc9+k ± (6xi2yi 2+ yi4)/zi9+k (ijk = 240; 241);
Pijk = xci'2ycj(xc2 + yc2)/zc9+k і Xii-2yiJ(Xi2 + yi2)/Zi9+k (ijk = 500; 501; 410; 411; 400; 401; 310;
311);
Pijk = xciycj-2(xc2 + yc2)/zc9+k і XiiyiJ-2(Xi2 + yi2)/Zi9+k (ijk = 050; 051; 140; 141;040; 041; 130;
131);
Pijk = (Xc3 + 3xcyc2)(xc2 + yc2)/zc9+k і (xi3 + 3xiyi2)(xi2 + yi2)/z9+k (ijk = 320; 321);
Pijk = (Xc3 + 3xc2yc)(xc2 + yc2)/zc9+k і (xi3 + 3xi2yi)(xi2 + yi2)/zi9+k (ijk = 230; 231);
Pijk = (xc2 + yc2) /Zc +k і (xi2 + yi2)3/zi9+k (ijk = 220; 22l).
© Ю.Ц. Батомункуев, 2012