1. В неоднородных средах представление о распространении света вдоль лучей сохраняется, но сами лучи становятся криволинейными. Действительно, рассмотрим среду, состоящую из плоскопараллельных слоев с постоянными показателями преломления (рис. 14), меняющимися скачкообразно от слоя к слою. Световой луч, преломляясь на границах слоев, примет форму ломаной линии. Будем неограниченно увеличивать число слоев, устремляя к нулю их толщины и скачки показателя преломления. Тогда в пределе показатель преломления среды станет меняться в пространстве непрерывно, а луч перейдет в кривую с непрерывно изменяющейся касательной. Допустим теперь, что показатель преломления меняется в пространстве как угодно, но не слишком резко (см. следующий параграф). Проведем в среде поверхности равного показателя преломления. В пределах каждого достаточно малого объема эти поверхности можно считать плоскими, а среду плоскослоистой, к которой применимо рассуждение, приведенное выше.

Из изложенного следует, что геометрическую форму луча можно однозначно определить из вакона Снеллиуса путем предельного перехода. Но закон Снеллиуса получается не только в волновой, но и в корпускулярной теории Ньютона. Поэтому при определении формы светового луча можно рассуждать так, как если бы свет состоял из ньютоновых корпускул, а показатель преломления $n$ определялся формулой (3.3). Поскольку преломление определяется относительным показателем преломления, абсолютный показатель можно заменить величиной, ему пропорциональной. Ради краткости можно просто положить $n=v$. Скорость корпускулы $v$ однозначно определяется уравнением сохранения энергии, а потому б можно рассматривать как известную функцию координат. Таким образом, траектория корпускулы в потенциальном поле сил геометрически совпадает с лучом света в среРис. 14. $\partial e$, показатель преломления $n$ которой численно равен $v$. Для этого, конечно, необходимо, чтобы исходные направления этих двух кривых были одинаковы. Эта формальная аналогия между движением частицы и распространением светового луча позволяет перенести результаты, полученные в световой оптике, в электронную микроскопию, где роль световых лучей выполняют электроны, движущиеся в потенциальных электрических полях.

Воспользуемся отмеченной аналогией для вычисления радиуса кривизны $R$ светового луча. Нормальное ускорение корпускулы определяется формулой
\[
\frac{m v^{2}}{R}=F_{N}=-\frac{\partial U}{\partial N},
\]

где $F_{N}$ — составляющая действующей силы $\boldsymbol{F}$ вдоль единичного вектора главной нормали $N$, а $U$ — потенциальная энергия корпускулы. Дифференцируя вдоль $N$ уравнение энергии $1 / 2 m v^{2}+$ $+U=$ const, получаем $F_{N}=m v \partial v / \partial N$, а потому
\[
\frac{1}{R}=\frac{1}{v} \frac{\partial v}{\partial N} .
\]

Заменяя $v$ на $n$, находим выражение для кривизны луча:
\[
\frac{1}{R}=\frac{1}{n} \frac{\partial n}{\partial N}=\frac{\partial}{\partial N}(\ln n) .
\]

При этом ускорение корпускулы, а с ним и сила $\boldsymbol{F}$ не имеют составляющей вдоль бинормали $\boldsymbol{b}$ к траектории, т. е. $F_{b}=-\partial U / \partial b=$ $=m v \partial v / \partial b=0$. Отсюда $\partial v / \partial b=0$, а потому также $\partial n / \partial b=0$. Значит, вектор grad $n$ лежит в соприкасающейся плоскости светового луча. Поэтому из всех направлений, перпендикулярных к лучу, направление главной нормали $\boldsymbol{N}$ характеризуется самым быстрым изменением показателя преломления среды. Это значит, что в неоднородной среде луч изгибается в сторону наиболее быстрого изменения показателя преломления. Если среда однородна ( $n=$ const), то кривизна $1 / R$ обращается в нуль, т. е. световые лучи прямолинейны.
2. Плотность земной атмосферы, а с ней и показатель преломления убывают с высотой. Этим объясняется ряд явлений, связанных с искривлением световых лучей. К ним относится, например, астрономическая рефракция, т. е. кажущееся поднятие небесного светила из-за искривления световых лучей в земной атмосфере. Если светило стоит высоко над горизонтом, то при расчете этого явления мы не сделаем заметной ошибки, считая поверхность Земли плоской. Но так поступать нельзя, когда светило находится вблизи горизонта. В этом случае надо учитывать сферичность земной поверхности.
— Пренебрежем малыми боковыми градиентами показателя преломления воздуха $n$ и будем считать, что $n$ зависит только от высоты над земной поверхностью или, что то же самое, от расстояния $r$ от центра земного шара $O$ (рис. 15). Световой луч $A M$ от небес-

Рис. 15. ного светила будет лежать в вертикальной плоскости, проходящей через это светило и глаз наблюдателя. На основании $(4.1)$
\[
\frac{1}{R}=\frac{d \alpha}{d s}=\frac{\partial}{\partial N}(\ln n)=\frac{d(\ln n)}{d r} \frac{\partial r}{\partial N}=-\frac{d(\ln n)}{d r} \sin \beta,
\]

где $\alpha$ — так называемое зенитное расстояние, т. е. переменный угол между вертикалью места $M Z$ и касательной к лучу, а $s$ — длина луча, отсчитываемая от глаза набйюдаля $M$, Қак видно из рис. 15 , .
\[
\frac{d \alpha}{d s}=\frac{d \alpha}{d r} \frac{d r}{d s}=\frac{d \alpha}{d r} \cos \beta,
\]
a потому
\[
\frac{1}{\operatorname{tg} \beta} \frac{d \alpha}{d r}=-\frac{d}{d r}(\ln n) .
\]

Далее, $r d \gamma=d \mathrm{~s} \sin \beta$, откуда
\[
\frac{1}{\operatorname{tg} \beta} \frac{d \gamma}{d r}=\frac{1}{r} \frac{d s}{d r} \cos \beta=\frac{1}{r} .
\]

Вычитая это равенство из (4.2) и приннмая во внимание, что $\alpha-\gamma \approx \beta$, получим
\[
\frac{d \beta}{\operatorname{tg} \beta}=-d \ln n-\frac{d r}{r}=-d \ln (n r) .
\]

Интегрирование этого уравнения дает
\[
n r \sin \beta=n_{0} r_{0} \sin \alpha_{0},
\]

где нуликом обозначены величины $n, r, \alpha$ в точке $M$ (в которой $\alpha \equiv \beta$ ), Вычислив отсюда $\operatorname{tg} \beta$ и подставив его значение в (4.2), найдем
\[
\alpha_{\infty}-\alpha_{0}=-\bar{n}_{0} r_{0} \sin \alpha_{0} \int_{r_{0}}^{\infty} \frac{d \ln n}{d r} \frac{d r}{\sqrt{n^{2} r^{2}-n_{0}^{2} r_{0}^{2} \sin ^{2} \alpha_{0}}} .
\]

Здесь $\alpha_{0}$ — видимое зенитное расстояние светила в точке $M$, а $\alpha_{\infty}$ — угол между асимптотой луча и вертикалью в той же точке. Разность $\alpha_{\infty}-\alpha_{0}$ называется рефракцией. Для ее вычисления и служит формула (4.4). Зная зависимость плотности воздуха $\rho$ от высоты над земной поверхностью, находят $n$ по формуле $(n-1) / \rho=$ const, после чего производят численное интегрирование в формуле (4.4).

Если светило находится не слишком близко от горизонта, то кривизной земной поверхности можно пренебречь. Тогда рефракция не зависит от характера изменения показателя преломления с высотой и может быть вычислена по формуле
\[
\sin \alpha_{\infty}=n_{0} \sin \alpha_{0}
\]
(см. задачу 3 к § 2). Ошибка, даваемая этой формулой, быстро убывает с уменьшением $\alpha$. Уже при $\alpha_{0}=80^{\circ}$ она не превосходит $3 \%$. При $\alpha_{0}<75^{\circ}$ формула (4.5) дает достаточную точность.

Если светило находится на горизонте ( $\alpha=90^{\circ}$ ), то средняя рефракция при $10^{\circ} \mathrm{C}$ и 760 мм рт. ст. составляет $35^{\prime} 24^{\prime \prime}$. Она быстро уменьшается по мере поднятия светила над горизонтом. Уже при $\alpha_{10}=89^{\circ}$ средняя рефракция уменьшается до $24^{\prime} 37^{\prime \prime}$. Этим объясняется сплюснутая форма Солнца при восходе и заходе.

С рефракцией связано некоторое удлинение дня. Для средних широт бі составляет в среднем 3-4 минуты.

Рефракция зависит от длины волны. Поэтому при заходе Солнца сначала должны исчезать красные и желтые лучи, а оставшийся сегмент солнечного диска должен окрашиваться на одну-две секунды в зеленый или даже синий цвет. При восходе Солнца, наоборот, должна сначала появляться кратковременная зеленая вспышка. Это явление «зеленого луча» наблюдается на море, да и то крайне редко, так как необходимы исключительно спокойная атмосфера и ясная погода.
3. Иногда вблизи земной поверхности из-за сильного нагревания или охлаждения возникают большие градиенты показателя преломления воздуха. Тогда шарообразность Земли можно не учитывать. Если градиент $n$ направлен вертикально, то можно воспользоваться формулой (4.3), полагая в ней $\beta=\alpha, r=r_{0}$. В результате получится
\[
n \sin \alpha=n_{0} \sin \alpha_{0}=\text { const. }
\]

Допустим ради определенности, что световой луч распространяется вверх под малым углом к горизонту, а показатель преломдения $n=n$ (z) убывает с ғысотой $z$. Может случиться, что на некоторой высоте угол $\alpha$ обратится в $90^{\circ}$. Тогда касательная к лучу станет горизонтальной. Высота $z=h$, на которой это может произойти, определится из соотношения
\[
n(h)=\vec{n}_{0} \sin \alpha_{0} .
\]

Так как показатель преломления возрастает вниз, то, достигнув высоты $z=h$, луч в дальнейшем должен загнуться книзу. На указанной высоте происходит нечто аналогичное полному отражению (рис. 16). Такое явление может возникнуть при сильно аномальном распределении плотности воздуха по высоте и является причиной различного рода миражей, наблюдаемых в атмосфере.
Обычно наблюдается вер $x$ ний или нижний мираж. При верхнем мираже, помимо самих предметюз, видны их изображения, расположенные сверху; при нижнем мираже изображение получается ниже самого предмета. Нижний мираж наблюдается в нустынях й в степях в теплое время года, когда прилегающий к земной поверхности слой воздуха сильно нагрет, а его плотность и показатель преломления быстро возрастают с высотой. Из көждой точки предмета в глаз наблюдателя всегда попадают прямые лучи, $t, e$, лучи, не испытавшие полного отражения в воздухе; им соответствует обычное — прямое — изображение предмета. Но при больших градиентах показателя преломления могут также попасть лучи, испытавшие полное отражение. Они дают обратное изображение предмета, как в зеркале (рис. 17). При этом лучи, выходящие из различных точек предмета, претерпевают полное отражение на несколько разных высотах; точка $M$, в которой отражается луч, выходящий из $A$, лежит несколько выше точки $N$, в которой отражается луч, выходящий из $B$. Наблюдатель видит два изображения: прямое и обратное. Создается иллюзия водной поверхности, в которой, как в зеркале, видно изображение неба.

Аналогично объясняется и верхний мираж. Он наблюдается зимой в холодных странах, когда вблизи земной поверхности образуется холодный слой воздуха, в котором показатель преломления быстро убывает с высотой. В горах, хотя и очень редко, наблюдается боковой мираж, связанный с изменением показателя преломления воздуха в боковом направРис, 17. лении. Вообще, в зависимости от характера распределения показателя преломления воздуха, мираж проявляется в весьма разнообразных и часто причудливых формах.
З А д А ч и
1. Показать, что в пренебрежении кривизной земной поверхности справедлива формула
\[
\alpha_{\infty}-\alpha_{0}=\left(n_{0}-1\right) \operatorname{tg} \alpha_{0} .
\]
2. Определить, при каких градиентах температуры $T$ воздуха возможен нижний мираж.

Р ешение. Показатель преломления, а с ним и плотность воздуха $\rho$ должны увеличиваться с высотой $h$ (т. е. $d \rho / d h>0$ ), Воспользуемся уравнением -состояния $\rho=\mu P /(R T)$, где $P$ — давление, $\mu$ — относительная молекулярная масса воздуха, а $R$ — универсальная газовая постоянная, Из него находим
\[
\frac{1}{\rho} \frac{d \rho}{d h}=\frac{1}{P} \frac{d P}{d h}-\frac{1}{T} \frac{d T}{d h} .
\]

При механическом равновесии воздуха $d P / d h=-\rho g$, где $g-$ ускорение свободного падения, В результате получаем
\[
\frac{d T}{d h}<\frac{\mu g}{R} .
\]

Используя соотношение $c_{P}-c_{v}=R / \mu$, этому условию можно придать вид
\[
\frac{d T}{d h}<-\frac{g}{c_{P}-c_{v}} \approx-0,025 \mathrm{~K} / \mathrm{M},
\]

где $c_{p}$ и $c_{v}$ — удельные теплоемкости воздуха при постоянном давлении и постоянном объеме. Такое распределение температур конвективно неустойчиво, так как лля конвективной устойчивости необходимо
(сM, T. II
\[
\frac{d T}{d h} \geqslant-\frac{g}{c_{P}} \approx-0,01 \mathrm{~K} / \mathrm{M}
\]