A1  ^0.45 Постройте точно ход лучей в плоско-выпуклых и в плоско-вогнутой линзе. Радиусы кривизны поверхностей $R=5 см$, показатель преломления стекла $n=1.5$. Укажите точку $F'$ схождения (расхождения) параллельного параксиального пучка.

Заметим, что радиус линзы составляет 4 см, толщина линзы – 2 см, а падающие лучи нарисованы через каждый сантиметр от оптической оси линзы.
Для плоско-выпуклой (первая) и плоско-вогнутой (третья) линзы преломление происходит только на одной поверхности. Положение фокуса для первой и третей линзы отмечаем используя формулу тонкой линзы:
$$F = \frac{R}{n-1} = 10 см.$$
Это расстояние отсчитывается от изогнутой поверхности.
Для выпукло-плоской (вторая) линзы преломление происходит сразу на двух поверхностях. Чтобы найти положение фокуса, учтём, что такая линза является толстой. После первого преломления лучи внутри стекла будут собираться на расстоянии $\frac{n}{n-1}R=15 см$ от изогнутой поверхности, то есть на расстоянии 13 см от плоской поверхности. При преломлении на плоской поверхности изображение окажется в $n$ раз ближе, то есть на расстоянии 8.67 см. Итого, получаем фокус на расстоянии 10.67 см от изогнутой поверхности.
Чтобы построить преломление лучей, используем транспортир, линейку и закон Снелиуса:
$n_1 \sin \alpha = n_2 \sin \beta,$
где $\alpha$ – угол падения, $\beta$ – угол преломления, $n_1$ и $n_2$ – коэффициенты преломления сред по разные стороны границы.

A2  ^0.30 Измерьте величину максимальной сферической аберрации $\delta s'$ для всех построений. Запишите результаты (в мм) в лист ответов.

По результатам построений измеряем максимальную величину аберрации. Для выпукло-плоской эффект аберраций оказывается меньше. Ниже приведены результаты с точностью до десятых мм, полученные путём расчётов.

B1  ^0.25 Чему равно фокусное расстояние $f'$ для тонкой плоско-выпуклой линзы?

Для тонкой линзы
$$\frac{1}{f'} = (n-1) \left(\frac{1}{R_1}+\frac{1}{R_2} \right),$$
в нашем случае $R_1 = \infty$, $R_2 = R$. Отсюда получаем ответ.

B2  ^1.50 Для данной линзы найдите зависимость величины продольной сферической аберрации $\delta s'$ от полуширины пучка $h$. Ответ выразите через $R, n$ и $h$.
Запишите приближенное выражение для $\delta s'(h)$, разложив ответ в ряд Тейлора по степеням $h$ до $h^2$.
Постройте схематический график полученной зависимости.

Согласно геометрии (рис. B2) и закону Снеллиуса:
$$n \sin \alpha = \sin (\alpha+\beta).$$

Выразим углы через отрезки и показатель преломления:
$$\alpha = \arcsin \frac{h}{R}, \quad \alpha+\beta = \arcsin \frac{nh}{R}.$$

Расстояние $l$ вдоль оптической оси от точки преломления до точки пересечения с оптической осью:
$$l = h \, \mathrm{ctg} \beta = h \, \mathrm{ctg} \left( \mathrm{arcsin} \left(\frac{n h}{R} \right)- \mathrm{arcsin} \left(\frac{h}{R} \right) \right).$$

Преломление происходит в точке, смещенной на
$$\delta x = R-\sqrt{R^2-h^2}.$$

Выразим $\delta s'$:
$$\delta s' = l -\delta x -f'.$$

Подставим $l$ и $\delta x$ найденные ранее, подставим $f'$ из B1, получим ответ.

B3  ^0.50 Найдите величину поперечной сферической аберрации $\delta g'$. Ответ выразите через $\delta s', f'$ и $h$.

Из подобия треугольников (рис. B2) получаем
$$\frac{h}{l}=\frac{\delta g'}{|\delta s '|}.$$
Считая аберрации малыми $(l=f')$, выражаем ответ.

B4  ^1.50 Разобьем кружок рассеяния кольцами радиусов $0.2\delta g', 0.4\delta g', 0.6\delta g', 0.8\delta g'$.
С какой доли (в %) поперечной площади падающего пучка собираются лучи в центральный круг (радиуса $0.2\delta g'$)?
С какой доли (в %) поперечной площади падающего пучка собираются лучи в последующие кольца (между $0.2\delta g'$ и $0.4\delta g'$ и так далее соответственно)?
Какую долю площади кружка рассеяния занимает каждое кольцо? Сведите все расчеты в таблицу в листе ответов.
Постройте схематический график зависимости интенсивности света от радиальной координаты в кружке рассеяния.

Считаем доли площади каждого кольца:
$$0 \div 0.2\delta g' : \quad 0.2^2-0^2 = 0.04 = 4%;$$
$$0.2\delta g' \div 0.4\delta g' : \quad 0.4^2-0.2^2 = 0.12 = 12%;$$
$$0.4\delta g' \div 0.6\delta g' : \quad 0.6^2-0.4^2 = 0.20 = 20%;$$
$$0.6\delta g' \div 0.8\delta g' : \quad 0.8^2-0.6^2 = 0.28 = 28%;$$
$$0.8\delta g' \div 1.0\delta g' : \quad 1.0^2-0.8^2 = 0.36 = 36%;$$

Из пункта B2 получаем $|\delta s'| \sim h^2$. Учитывая пункт B3 получаем, что $\delta g' \sim h^3.$

Считаем доли площади колец в падающем пучке:
$$0 \div 0.2\delta g' : \quad 0.2^{2/3}-0^{2/3} = 0.34 = 34%;$$
$$0.2\delta g' \div 0.4\delta g' : \quad 0.4^{2/3}-0.2^{2/3} = 0.20 = 20%;$$
$$0.4\delta g' \div 0.6\delta g' : \quad 0.6^{2/3}-0.4^{2/3} = 0.17 = 17%;$$
$$0.6\delta g' \div 0.8\delta g' : \quad 0.8^{2/3}-0.6^{2/3} = 0.15 = 15%;$$
$$0.8\delta g' \div 1.0\delta g' : \quad 1.0^{2/3}-0.8^{2/3} = 0.14 = 14%;$$

Если разделить долю площади в падающем кольце на долю площади кольца в фокальной плоскости, то получим величину пропорциональную интенсивности. Получается качественно график $I\sim r^{-4/3}$.

C1  ^0.30 Используя теорему синусов для треугольников $MSC$ и $MS'C$, а также закон преломления Снеллиуса, получите инвариант вида $n\cdot f(r, s, p)=n'\cdot f(r, s', p')=\mathrm{inv}$.

Теорема синусов для треугольника $MSC$ (см. рис. 2 из условия):
$$\frac{\sin \theta}{p} = \frac{\sin \angle SMC}{r-s}.$$
Теорема синусов для треугольника $MS'C$:
$$\frac{\sin \theta}{p'} = \frac{\sin \angle S'MC}{r-s'}.$$
Закон преломления:
$$n \sin \angle SMC = n' \sin \angle S'MC.$$

Выразим из теоремы синусов $\sin \angle SMC$ и $\sin \angle S'MC$, подставим в закон преломления, получим инвариант

C2  ^0.10 Воспользовавшись теоремой косинусов в треугольнике $MSC$, перейдите в инварианте от величины $p$ к углу $\theta$.
Считая угол $\theta$ малым, разложите множители в виде корней в ряд Тейлора до первого слагаемого с ненулевой степенью $\theta$.
Разделите полученное приближенное значение на «$-r$» и получите приведенное значение инварианта.

Теорема косинусов
$$p^2 = r^2 + (r-s)^2 - 2 r (r-s) \cos \theta.$$

Используя разложение косинуса и полный квадрат
$$p^2 = (r - (r-s))^2 + 2 r(r-s) \cdot \frac{ \theta^2 }{ 2}.$$

Подставляем в инвариант, получаем:
$$\cfrac{n(r-s)}{\sqrt{s^2+r(r-s)\theta^2}}=\cfrac{n(r-s)}{s} \left( 1 - \frac{r(r-s)}{2s^2} \right).$$

Разделим на «$-r$», получим ответ.

C3  ^0.10 Записав инвариант из пункта C2 для параксиальных лучей до и после преломления, покажите инвариантность величины $Q(s_0)$.

Положим $\theta = 0$ в инварианте пункта C2.

C4  ^0.10 Левая часть инварианта, полученного в С2, является функцией $f$ от $s$ и $\theta$. Найдите значение $f(s,\theta)-f(s_0,\theta)$. Считайте, что $\delta s$ и $\theta$ являются малыми параметрами.

Инвариант является функцией от $s$ и от $\theta$:
$$f(s,\theta) = n\left(\cfrac{1}{r}-\cfrac{1}{s}\right)\left(1-\cfrac{r(r-s)\theta^2}{2s^2}\right).$$
Находим значение разницы, учитываем что $\theta$ одинаковый:
$$f(s,\theta)-f(s_0,\theta) = \frac{\partial f}{\partial s}ds.$$

Так как $\theta$ тоже малая величина, то
$$\frac{\partial f}{\partial s} = n \frac{1}{s^2}$$

C5  ^0.20 Используя результаты предыдущих пунктов, выразив малый угол $\theta$ через радиус кривизны поверхности и малую высоту $h$, получите искомое соотношение:
$$
\cfrac{n'\delta s'}{s'^2_0}-\cfrac{n\delta s}{s^2_0}=-\cfrac{h^2}{2} Q^2(s_0)\left(\cfrac{1}{n's'_0}-\cfrac{1}{ns_0}\right).
$$

Выразим угол $\theta$ через $h$:
$$\theta = \frac{h}{r}.$$

Представим выражение через $Q(s)$:
$$\cfrac{r(r-s)\theta^2}{2s^2}=\cfrac{(r-s)h^2}{2s^2 r}=n\left(\frac{1}{s}-\frac{1}{r} \right) \frac{h^2}{2ns}=Q(s) \frac{h^2}{2ns}.$$

Следующие величины инвариантны
$$f(s, \theta) = f(s', \theta),$$
$$f(s_0, \theta) + \frac{n \delta s_0}{s_0^2} = f(s'_0, \theta) + \frac{n' \delta s'}{{s'}_0^2},$$
$$ n\left(\cfrac{1}{r}-\cfrac{1}{s_0}\right)\left(1-\cfrac{r(r-s_0)\theta^2}{2s_0^2}\right) + \frac{n \delta s}{s_0^2} = n\left(\cfrac{1}{r}-\cfrac{1}{s'_0}\right)\left(1-\cfrac{r(r-s'_0)\theta^2}{2{s'}_0^2}\right) + \frac{n' \delta s'}{{s'}_0^2},$$
$$ Q(s_0) - \frac{h^2}{2ns_0}Q^2(s_0) + \frac{n \delta s}{s_0^2} = Q(s'_0) - \frac{h^2}{2n's'_0}Q^2(s'_0) + \frac{n' \delta s'}{{s'}_0^2}$$

Пользуемся тем, что $Q(s_0)=Q(s'_0)$, получаем выражение
$$
\cfrac{n'\delta s'}{s'^2_0}-\cfrac{n\delta s}{s^2_0}=-\cfrac{h^2}{2} Q^2(s_0)\left(\cfrac{1}{n's'_0}-\cfrac{1}{ns_0}\right).
$$

C6  ^0.20 Покажите, что аберрации предмета и изображения после $k$-й преломляющей поверхности, связаны соотношением:
$$
n_k'\, \delta s_k' \, \alpha_k'^2 - n_1\, \delta s_1\, \alpha_1^2=-\cfrac{1}{2}\sum_{i=1}^{k} h_i^4\, Q^2_i(s)\, \left(\cfrac{1}{n'_{i}s'_{i}}-\cfrac{1}{n_is_i}\right).
$$
Иными словами, для вычисления суммарной аберрации нужно по всем поверхностям сложить правые части доказанного соотношения, умноженного на $h_i^2$.

Перейдём от $s_i$ к углам в левой части формулы из C5:
$$\alpha_i = \frac{h_i}{s_i}, \quad \alpha '=\frac{h_i}{s_i'}.$$

При прохождении луча между преломляющими поверхностями меняется расстояние от луча до оптической оси, однако, сохраняется угол наклона:
$$\alpha_i'=\alpha_{i+1},$$
$$n'_i = n_{i+1},$$
$$\delta s'_i = \delta s_{i+1}.$$

Складывая такие выражения
$$
n_i'\, \delta s_i' \, \alpha_i'^2 - n_i\, \delta s_i\, \alpha_i^2=-\cfrac{h_i^4}{2} Q_i^2(s_i)\left(\cfrac{1}{n'_is'_i}-\cfrac{1}{n_i s_i}\right).
$$
в левой части остаётся только два слагаемых: одно соответствует вхождению лучей в оптическую систему, другое – вхождению из неё. Тем самым можно вычислять аберрации для сложных систем.

D1  ^2.00 Найдите сферическую аберрацию $\delta s_2'$ тонкой собирающей линзы для падающего параллельного пучка света. Примените результат пункта C6 для двух преломляющих поверхностей и формулу тонкой линзы, оставив зависимость только от радиуса кривизны первой поверхности (куда падает исходный пучок света). Ответ выразите через показатель преломления $n$ и обратные величины $\rho=1/r$ и $\varphi=1/f'$.

Будем использовать соглашение о знаках как в части C. Тогда
$$s_1 = \infty, \quad s'_2 = -f.$$
Положение изображения после первого преломления:
$$s'_1=s_2 = r \frac{n}{n-1}.$$
Формула тонкой линзы позволяет найти радиус кривизны $r_2$ второй поверхности при заданном фокусе $f$ и радиусе первой поверхности $r$:
$$\frac{1}{f} = (n-1)\left( -\frac{1}{r} + \frac{1}{r_2} \right).$$
Так как линза тонкая, то $h_1=h_2=h$.
Левая часть формулы C6:
$$n_2'\, \delta s_2' \, \alpha_2'^2 - n_1\, \delta s_1\, \alpha_1^2 = \delta s_2' \, \left( \frac{h}{f} \right)^2.$$
$$Q_1=\frac{1}{r} = \rho.$$
$$ Q_2 = \frac{1}{r_2} - \frac{1}{-f} = \frac{n}{n-1}\varphi + \rho.$$

Собираем всё вместе:
$$\delta s_2' = -\frac{h^2}{2 \varphi^2} \left( Q_1^2\frac{n-1}{n^2}\rho + Q_2^2 \left(-\varphi -\frac{n-1}{n^2}\rho \right) \right)$$

получаем ответ.

D2  ^0.50 Найдите значение $\rho^*$, при котором сферическая аберрация $\delta s_2'$ минимальна. Ответ выразите через $n$ и $\varphi$.

Выражение в квадратных скобках в D1 представляет собой параболу относительно $\rho$, ветви вверх. Минимальная аберрация (по модулю) будет соответствовать вершине. Отсюда получаем значение для $\rho^*$.

Заметим, что данные рассуждения верны только лишь тогда, когда парабола не пересекает горизонтальную ось. Ведь если она пересекает, то вообще говоря можно получить нулевую аберрацию. Предположим, что пересекает. Тогда дискриминант должен быть неотрицательным:
$$D=\frac{1-4 n}{(n-1)^2}\ge0,$$
$$n\le1/4.$$
Но для собирающей линзы $n>1$, поэтому, в действительности парабола не пересекает горизонтальную ось в наших условиях.

D3  ^0.50 Рассчитайте значение радиусов кривизны поверхностей линзы. Нарисуйте (качественно) вашу рассчитанную линзу в листе ответов.

D4  ^0.50 Схематично изобразите график зависимости модуля продольной сферической аберрации $\delta s_2'$ для линз, приведенных на рисунке 4.

E1  ^1.50 В какую сторону и на какое расстояние нужно сместить плоскость наблюдения от фокальной плоскости, чтобы кружок рассеяния был минимального размера?

Используется факт из части B, что аберрация квадратична по удалённости луча от оптической оси $h$:
$$\delta s = \frac{h^2}{\varepsilon}, \quad \delta s' = \frac{H^2}{\varepsilon},$$
где $\varepsilon$ – некий коэффициент (размерности длины). Здесь $H$ – радиус линзы (максимально удалённый луч), $h$ – произвольный луч. Аналогично, $\delta s'$ – максимальная аберрация, $\delta s$ – точка пересечения оптической оси лучем $h$. Если что тут $\delta s$ и $\delta s'$ положительные, отчитываются от фокуса в сторону линзы.
Расстояние $r$ от оси до точки луча на расстоянии $x$ от фокуса, который до преломления идёт на расстоянии $h$ от оси можно найти из подобия (рис. E1):
$$
\frac{r}{x-\delta s} = \frac{h}{f}.
$$
Получается что $r$ есть функция плоскости наблюдения $x$ и луча $h$:
$$r(x, h) =\frac{1}{f}\left (xh - \frac{h^3}{\varepsilon}\right).$$
Находим, при каком $h$ при некотором фиксированном $x$ будет максимум
$$h_1 = \sqrt{\frac{\varepsilon}{3}}x^{1/2}.$$
Максимальный радиус пучка в зависимости от плоскости наблюдения
$$r_{\mathrm{max}}(x) = \frac{2 \varepsilon^{1/2} x^{3/2}}{3^{3/2}f}.$$
Это получается некоторая огибающая, которую мы получаем, смотря на лучи, которые не пересекают оптическую ось. Чем ближе к фокусу, тем меньше размер кружка, который дают эти лучи. Однако, ещё есть лучи, которые пересекают оптическую ось и при приближении к фокусу увеличивают размер кружка.
Зависимость расстояния от оси до точки луча для крайнего противоположного луча от координаты
$$y(x) = \frac{H}{f}(\delta s' - x).$$
Минимальный кружок будет там, где эти эффекты равны $r_{\mathrm{max}}(x) = y(x)$:
$$\frac{2 \varepsilon^{1/2} x^{3/2}}{3^{3/2}f} = \frac{H}{f}(\delta s' - x),$$
$$\frac{2 }{3^{3/2}} \left( \frac{x}{\delta s'} \right)^{3/2} = \left(1 - \frac{x}{\delta s'} \right).$$
Решая это уравнение, получаем ответ
$$x = \frac{3}{4} \delta s'.$$

E2  ^1.50 Постройте схематический график зависимости интенсивности света от радиальной координаты в минимальном кружке рассеяния. Для этого выполните шаги, аналогичные пункту B4.

Ответим в начале на вопрос: как радиус окружности ($r$) в плоскости E1 соотносится с радиусом окружности ($h$) в пучке, падающем на линзу?
Из прошлого пункта оказывается, что они связаны следующим образом:
$$r = h |3-4h^2|,$$
здесь $r$ – отнормировано на максимальный радиус пятна $ \left( \frac{1}{4} \frac{\delta s' H}{f} \right)$ в рассматриваемой плоскости, $h$ – на радиус пучка $H$.
В общем случае одному радиусу в плоскости наблюдения соответсвует целых три радиуса в падающем пучке. Доля площади кружка оказывается такой же, как и в прошлый раз. Чтобы найти долю пучка, нужно сложить долю площади трех колец, то есть величины вида $h_2^2-h_1^2$. Для нахождения $h_1$ и $h_2$ каждого кольца нужно решить уравнение третей степени.