Олимпиады · APhO

APhO 2022, теоретический тур

3 задачи

зад 1 Эксперимент Штерна–Герлаха

Магнетизм и индукция › Магнитное полеМагнитный дипольМагнитное поле токовРаспределения, барометрическая формула, модели атмосферысложность 6/10

Авторы задачи: Х. С. Мани (бывший директор HRI, Прайаградж) и Гаутам Датта (DAIICT, Гандинагар). Вклад Академического комитета, Группы академического развития и Международного жюри признателен.

Эксперимент Штерна–Герлаха был выполнен в 1922 году и в конечном счёте привёл к определению магнитного момента электрона. В этом эксперименте пучок атомов серебра, каждый массой \(m = 1{,}80\times10^{-25}\) кг, выходит из печи, поддерживаемой при температуре \(T = 1{,}20\times10^{3}\) К (см. рисунок 1). Считайте, что после выхода из печи все атомы имеют одинаковый импульс вдоль направления пучка (\(z\)-направление). Силой тяжести пренебрегайте.

Рисунок

Рисунок 1 — Принципиальная схема установки Штерна–Герлаха. На рисунке: Ag gas — пар серебра (Ag); Oven — печь; T = 1200K — температура печи.

A-1. Скорость атомов серебра: скорость \(v_z\) атомов серебра, вылетающих из печи, можно оценить как \(\sqrt{3k_BT/m}\) по теореме о равнораспределении. Вычислите это значение. (\(0{,}5\) балла)

Основное выражение. После выхода из печи атомы серебра движутся вдоль \(z\)-направления на расстоянии \(l_1 = 0{,}25\) м. Далее атомы серебра проходят между двумя магнитами на расстоянии \(l_2 = 0{,}5\) м. Магниты создают неоднородное магнитное поле \(B\) в \(x\)-направлении с постоянным градиентом \(dB/dx\). Считайте, что магнитный момент атома серебра направлен либо в \(+x\), либо в \(-x\), то есть \(\vec\mu_s = \pm\mu_s\hat\imath\). После прохождения магнитов атомы серебра проходят ещё расстояние \(l_3 = 0{,}25\) м до того, как ударяются об экран \(PP'\). Расстояние между двумя пучками, ударяющими экран, равно \(\Delta x\).

B-1. Выведите выражение для расстояния расщепления \(\Delta x\) на экране. (\(2\) балла)

Неоднородное магнитное поле. Эта часть посвящена установке для создания неоднородного магнитного поля (\(dB/dx \neq 0\)). Она состоит из нескольких подпунктов. Два очень длинных провода, параллельных оси \(z\), несут токи величиной \(I_0\) и расположены в точках \(A_1(0,-a,z)\) и \(A_2(0,a,z)\) (см. рисунок ниже). Ток, текущий по проводу при \(y=-a\), направлен вдоль \(-\hat k\), а ток, текущий по проводу при \(y=a\), направлен вдоль \(\hat k\). Вся система находится в среде с высокой относительной магнитной проницаемостью \(\mu_r\). Положим \(\mu = \mu_0\mu_r\). Провода изолированы, и ток в среду не утекает.

Рисунок

Рисунок 2 — Расположение для создания неоднородного магнитного поля. На рисунке: Air — воздух.

C-1. Выведите выражение для вектора магнитного поля в точке \(P_1(x,y,0)\) в плоскости \(x\)–\(y\) (см. рисунок 2). (\(1{,}5\) балла)

C-2. Рассмотрите окружность с центром \(C\) на оси \(x\) \((x_c,0)\) и радиусом, равным расстоянию \(AC\) (см. рисунок 2). Определите направление магнитного поля в точках окружности \(R\) (на оси \(x\)) и в точке \(P_0\) (\(CP_0\) параллельна оси \(y\)). (\(0{,}5\) балла)

C-3. Теперь небольшой участок среды с высокой магнитной проницаемостью между окружностями радиусов \(AC\) и \(AD\) удалён и заменён воздухом очень низкого давления (см. рисунок 2). Из соображений непрерывности можно показать, что магнитное поле в этой области зазора задаётся тем же выражением, что и при неудалённой среде. (Мы считаем это известным; доказывать это не требуется.) Запишите выражение для магнитного поля в точке \((x,0)\) в области воздушного зазора. (\(0{,}5\) балла)

Сила. Как уже говорилось, атомы серебра движутся в плоскости \((x,0,z)\), их скорости параллельны оси \(z\) и равны \(\vec v = v_z\hat k\). Напомним также, что магнитный диполь атома серебра равен \(\vec\mu_s = \pm\mu_s\hat\imath\).

D-1. Выведите выражение для модуля силы \(F_x\), действующей на атом серебра вдоль направления \(x\), через \(\mu_s\), \(I_0\), \(a\mu\) и соответствующие координаты. (\(0{,}5\) балла)

Магнитное поле и его градиент. Считаем, что эта же сила \(F_x\) действует на малом расстоянии \(l_2\) вдоль оси \(z\) (рисунок 1). Считаем также, что атомы серебра проходят через средний пункт \(P\) отрезка \(RQ\) (рисунок 2). Заданы следующие экспериментальные значения: \[\frac{\mu}{\mu_0} = 10^4; \qquad a = 0{,}60\ \text{см}; \qquad OC = 0{,}60\ \text{см}; \qquad OD = 0{,}80\ \text{см}; \qquad I_0 = 2{,}00\ \text{А}.\] Здесь \(\mu_0\) — магнитная проницаемость вакуума.

E-1. Найдите численные значения модуля магнитного поля \(B_P\) и его градиента \(dB_P/dx\) в этой средней точке в единицах СИ. (\(2{,}0\) балла)

Магнитный момент атома серебра. Если \(v_z = 500\) м\(\cdot\)с\(^{-1}\), а магнитное поле рассчитано, как указано выше, то эксперимент Штерна–Герлаха даёт расщепление \(\Delta x = 0{,}20\) см.

F-1. Найдите значение магнитного момента атома серебра \(\mu_s\) в единицах СИ. (\(1{,}5\) балла)

Разброс линии. Скорости атомов серебра могут быть не одинаковы. Пусть разброс скорости пучка составляет \(20\,\%\).

G-1. Каким будет соответствующий разброс пятна \(\delta x\) на экране? (\(0{,}5\) балла)

Погрешность магнитного момента.

H-1. Какова соответствующая погрешность \(\delta\mu_s\) при определении магнитного момента? (\(0{,}5\) балла)

▶ПодсказкаСкрыть
Подсказка

Сила Штерна–Герлаха действует на магнитный диполь только там, где есть градиент поля: \(F_x = -\mu_s\,\partial B_x/\partial x\). Смещение пучка складывается из отклонения внутри магнита и дальнейшего дрейфа с поперечной скоростью по пути до экрана.

▶ОтветСкрыть
Ответ

A-1) \(v_z = \sqrt{3k_BT/m} \approx 5{,}255\times10^{2}\) м/с;  B-1) \(\Delta x = 2\dfrac{\mu_s}{m}\dfrac{dB}{dx}\dfrac{l_2}{v_z^2}\left(\dfrac{l_2}{2}+l_3\right)\);  C-1) \(B(x,y) = \dfrac{\mu I_0 a}{\pi r_1^2 r_2^2}\left[2xy\,\hat\jmath + (x^2-y^2+a^2)\,\hat\imath\right]\), где \(r_1^2 = x^2+(y-a)^2\), \(r_2^2 = x^2+(y+a)^2\);  C-2) в точке \(R\) поле направлено вдоль \(\hat\imath\), в точке \(P_0\) — вдоль \(\hat\jmath\);  C-3) \(B(x,0) = \dfrac{\mu I_0 a}{\pi(x^2+a^2)}\hat\imath\);  D-1) \(F_x = \dfrac{2\mu_s\mu I_0}{\pi}\,\dfrac{a\,x}{(x^2+a^2)^2}\);  E-1) \(B_P \approx 0{,}16\) Тл, \(dB_P/dx \approx 17{,}3\) Тл/м;  F-1) \(\mu_s \approx 1{,}04\times10^{-23}\) Дж/Тл;  G-1) \(\delta x \approx 0{,}04\) см (пятно от \(0{,}06\) см до \(0{,}14\) см от центра);  H-1) \(\mu_s = (1{,}04 \pm 0{,}42)\times10^{-23}\) Дж/Тл.

▶РешениеСкрыть
Решение

A-1. Используем теорему о равнораспределении. Пусть \(\overline{v^2}\) — средний квадрат скорости атомов серебра в печи при \(1200\) К. Тогда \(\dfrac{m\overline{v^2}}{2} = \dfrac{3k_BT}{2}\), откуда среднеквадратичная скорость равна \(5{,}255\times10^{2}\) м/с.

B-1. Длина \(l_1\) не входит в выражение. Ускорение атомов серебра в области \(l_2\) равно \(a = \dfrac{\mu_s}{m}\dfrac{dB}{dx}\) и направлено в \(+x\) или в \(-x\). Одновременно атомы имеют постоянную горизонтальную скорость \(v_z\). Проходя область \(l_2\) за время \(l_2/v_z\), атом получает отклонение \[\delta_1 = \frac12\,\frac{\mu_s}{m}\frac{dB}{dx}\left(\frac{l_2}{v_z}\right)^2.\] Для оставшейся части полёта атом имеет постоянную горизонтальную скорость \(v_z\) и постоянную поперечную скорость \(v_{x0} = \dfrac{\mu_s}{m}\dfrac{dB}{dx}\dfrac{l_2}{v_z}\). Из-за составляющей \(v_{x0}\) атом получает дополнительное отклонение \[\delta_2 = \frac{l_3 v_{x0}}{v_z} = \frac{\mu_s}{m}\frac{dB}{dx}\frac{l_2 l_3}{v_z^2}.\] Полное отклонение в \(+x\) направлении равно \(\delta_1+\delta_2\), а расщепление на экране в идеализированном случае вдвое больше: \(\Delta x = 2(\delta_1+\delta_2)\), откуда \[\Delta x = 2\,\frac{\mu_s}{m}\frac{dB}{dx}\,\frac{l_2}{v_z^2}\left(\frac{l_2}{2}+l_3\right).\]

C-1. Пусть \(\vec r_1 = \overrightarrow{A_1P} = x\hat\imath + (y-a)\hat\jmath\) и \(\vec r_2 = \overrightarrow{A_2P} = x\hat\imath + (y+a)\hat\jmath\). Тогда \[\vec B(x,y) = \frac{\mu I_0}{2\pi}\left[\hat k\times\frac{\vec r_1}{r_1^2} - \hat k\times\frac{\vec r_2}{r_2^2}\right] = \frac{\mu I_0}{2\pi r_1^2 r_2^2}\Big[(x\hat\jmath-(y-a)\hat\imath)(x^2+(y+a)^2) - (x\hat\jmath-(y+a)\hat\imath)(x^2+(y-a)^2)\Big],\] что после упрощения даёт \[\vec B(x,y) = \frac{\mu I_0 a}{\pi r_1^2 r_2^2}\left[2xy\,\hat\jmath + (x^2-y^2+a^2)\,\hat\imath\right].\]

C-2. Точка \(R\) имеет координаты \((x_c+\sqrt{x_c^2+a^2},\,0)\); подставляя \(y=0\), видим, что \(\hat\jmath\)-компонента обращается в нуль, поэтому \(\vec B \parallel \hat\imath\). Для точки \(P_0\) с координатами \((x_c,\,\sqrt{x_c^2+a^2})\) \(\hat\imath\)-компонента обращается в нуль, поэтому \(\vec B \parallel \hat\jmath\).

Альтернативно можно показать, что в любой точке окружности поле направлено радиально (то есть перпендикулярно окружности): уравнение окружности с центром \((x_c,0)\) и радиусом \(\sqrt{x_c^2+a^2}\) есть \(x^2-2xx_c+y^2=a^2\), а векторное произведение \(\overrightarrow{CP}\times\vec B\) пропорционально \(y\,(x^2-2xx_c+y^2-a^2)\,\hat k\), что на окружности равно нулю.

Рисунок
Рисунок

Рисунок к решению C-1 — расположение проводов и точки \(P\) (вверху) и окружности с центром \(C\) и радиусами \(AC\) и \(AD\) (внизу). На рисунке: Air — воздух; \(P_1(x,y)\) — точка наблюдения; \(P_0(x_c,y_c)\) — точка на окружности.

C-3. В воздушном зазоре поле по-прежнему задаётся выражением (2) при \(y=0\): \[\vec B(x,0) = \frac{\mu I_0 a}{\pi(x^2+a^2)}\,\hat\imath.\]

D-1. Сила, действующая на магнитный диполь вдоль \(x\): \[F_x = -\mu_s\frac{\partial B_x}{\partial x} = \mu_s\,\frac{\mu I_0}{\pi}\cdot\frac{2ax}{(x^2+a^2)^2}.\]

E-1. Подставляя \(\mu/\mu_0 = 10^4\), \(a = 6{,}00\times10^{-3}\) м, \(OC = 6{,}00\times10^{-3}\) м, \(OD = 8{,}00\times10^{-3}\) м и \(I_0 = 2{,}00\) А. В средней точке \(P\) имеем \(y = 0\) и \[x_P = OP = \frac{(1+\sqrt2)\times 0{,}6 + 1{,}8}{2}\times10^{-2}\ \text{м} = 1{,}624\times10^{-2}\ \text{м},\] где использованы \(OD = 0{,}8\times10^{-2}\) м и \(DA = 10^{-2}\) м. Тогда \[B_x(x_P,0) = 10^4\cdot\frac{\mu_0}{\pi}\cdot\frac{I_0 a}{x_P^2+a^2} = 10^4\cdot 4\times10^{-7}\cdot\frac{2\cdot0{,}6\times10^{-2}}{(1{,}624^2+0{,}6^2)\times10^{-4}} = 0{,}16\ \text{Тл}.\] Градиент: \[\left.\frac{\partial B_x}{\partial x}\right|_{x_P} = \frac{2x_P}{x_P^2+a^2}\,B_x(x_P,0) = \frac{2\times1{,}624\times10^{-2}}{(1{,}624^2+0{,}6^2)\times10^{-4}}\times0{,}16 \approx 17{,}3\ \text{Тл/м}.\]

F-1. Из выражения для расщепления \[\mu_s = \frac{m\,\Delta x}{2\,(\partial B/\partial x)_{x_P}}\cdot\frac{1}{\dfrac{l_2}{v_z^2}\left(\dfrac{l_2}{2}+l_3\right)} = \frac{1{,}8\times10^{-25}\cdot2\times10^{-3}}{2\cdot17{,}34}\cdot10^{6} = 1{,}04\times10^{-23}\ \text{Дж/Тл}.\]

G-1. Половины расщепления \(\Delta x/2\) от центра равны \(\dfrac{\mu_s}{m}\dfrac{dB}{dx}\dfrac{l_2}{v_z^2}\left(\dfrac{l_2}{2}+l_3\right)\), то есть зависят от скорости \(v_z\). Поэтому разброс скорости приводит к разбросу расщепления: \[\delta\!\left(\frac{\Delta x}{2}\right) = \left|\frac{\partial(\Delta x/2)}{\partial v_z}\right|\delta v_z = \frac{2\,(\Delta x/2)\,\delta v_z}{v_z} = 2\,(\Delta x/2)\times0{,}2 = 0{,}04\ \text{см}.\] Таким образом, линия находится на расстоянии от центра от \(0{,}1-0{,}04 = 0{,}06\) см до \(0{,}1+0{,}04 = 0{,}14\) см.

H-1. Из предыдущего пункта следует, что расщепление меняется от \(0{,}12\) см до \(0{,}28\) см, тогда как раньше оно было \(0{,}2\) см. Зависимость расщепления от магнитного момента линейна, поэтому магнитный момент меняется от \((0{,}12/0{,}2)\) до \((0{,}28/0{,}2)\) от исходного значения. Это даёт от \(0{,}62\times10^{-23}\) до \(1{,}46\times10^{-23}\) Дж/Тл. Полный разброс составляет \(0{,}84\times10^{-23}\) Дж/Тл около среднего значения \(1{,}04\times10^{-23}\) Дж/Тл, или \[\mu_s = (1{,}04\pm0{,}42)\times10^{-23}\ \text{Дж/Тл}.\]

зад 2 Механическая модель фазовых переходов

Механика › Вращение твёрдого телаУравнение вращения, раскрутка и торможениеПлавление и кристаллизацияГрафические методысложность 8/10

Авторы задачи: Ситикантха Дас (IIT Kharagpur) и Прамендра Ранджан Сингх (Principal, Narayan College, J.P. University). Вклад Академического комитета, Группы академического развития и Международного жюри признателен.

Кольцо радиуса \(R\) имеет бусинку пренебрежимо малого размера и массы \(m\), надетую на него. Кольцо вращается вокруг своего вертикального диаметра с угловой скоростью \(\omega\), как показано на рис. 1. Вместе с этим на бусинку действует противодействующая сила \(F_f\), которая пропорциональна силе нормальной реакции \(N\) и задаётся формулой \(F_f = fkN\), где \(f\) равно \(1\) или \(-1\). Рекомендуется использовать полярные координаты \(\{r,\theta\}\). По возможности выражайте ответы через \(\omega_c = \sqrt{g/R}\), где \(g\) — модуль ускорения свободного падения.

В фазовом переходе свободная энергия \(V(\mathcal M)\) [механический аналог — потенциальная энергия] зависит от намагниченности \(\mathcal M\) следующим образом: \[V(\mathcal M) = a(T)\,\mathcal M^2 + b(T)\,\mathcal M^4,\] где \(T\) — температура, \(b(T) > 0\), а \(a(T)\) меняет знак с изменением температуры. Мы попытаемся понять явление фазового перехода с помощью описанной модели.

Примечание 1. Для кругового движения радиуса \(R\) скорость в полярных координатах равна \(\dot{\vec r} = R\dot\theta\,\hat\theta\), а ускорение равно \(\ddot{\vec r} = -R\dot\theta^2\,\hat r + R\ddot\theta\,\hat\theta\). Здесь \(\hat r\) и \(\hat\theta\) — единичные векторы радиального и касательного направлений соответственно.

Примечание 2. Направление противодействующей силы \(\vec F_f\) обозначено через \(f\). Здесь \(f=+1\), если бусинка движется в направлении против часовой стрелки (роста) \(\theta\), и \(f=-1\), если бусинка движется по часовой стрелке, например \(f=\operatorname{sgn}(\dot\theta)\), где \(\operatorname{sgn}\) равно \(+1\) или \(-1\) в зависимости от знака аргумента.

Примечание 3. Вам могут пригодиться разложения (где \(\theta\) — в радианах и \(|x|\ll1\)), полезные для некоторых частей задачи: \[\sin\theta = \theta-\frac{\theta^3}{6}+\dots,\qquad \cos\theta = 1-\frac{\theta^2}{2}+\frac{\theta^4}{24}+\dots,\qquad (1+x)^n = 1+nx+\frac{n(n-1)}{2}x^2+\frac{n(n-1)(n-2)}{6}x^3+\dots\]

Рисунок

Рисунок 1 — Бусинка на вращающемся кольце. На рисунке слово «positive» означает положительное направление (направление возрастания \(\theta\), то есть \(\hat\theta\)); \(\hat r\) — радиальный единичный вектор.

Дальше мы будем рассматривать динамику бусинки в системе отсчёта, связанной с вращающимся кольцом, и для углов в диапазоне \(-\pi/2 < \theta < \pi/2\). Диаграмма свободных сил бусинки показана на рис. 2. Пренебрегайте всеми силами, кроме показанных на ней.

Рисунок

Рисунок 2 — Диаграмма свободных сил. На рисунке: \(m\omega^2R\sin\theta\) — сила инерции; \(mg\) — сила тяжести; \(N\) — нормальная реакция; \(F_f = fkN\) — противодействующая сила.

A-1. Запишите уравнения движения для радиальной \(F_r\) и касательной \(F_\theta\) компонент силы, действующей на бусинку. Считайте, что \(\theta\) возрастает против часовой стрелки. (\(0{,}5\) балла)

Для пунктов B-1–B-9 считайте \(k = 0\).

B-1. Укажите связь между углом (углами) равновесия \(\theta_0\) и величинами \(\omega\), \(\omega_c\). (\(1{,}0\) балла)

B-2. Качественно начертите график зависимости \(\theta_0\) (по оси \(y\)) от \(\omega/\omega_c\) (по оси \(x\)). (\(0{,}5\) балла)

B-3. Качественно начертите график величины силы нормальной реакции бусинки в зависимости от \(\omega/\omega_c\) при устойчивом равновесии. (\(0{,}5\) балла)

B-4. Определим потенциальную энергию, соответствующую касательной силе \(F_\theta\): \[F_\theta = -\frac{1}{R}\frac{d}{d\theta}V(\theta),\] с нулём потенциальной энергии при \(\theta=0\). Если \(V(\theta)\) записана в виде \(P + Q\cos\theta + S\sin^2\theta\), найдите \(P\), \(Q\), \(S\). (\(1{,}0\) балла)

B-5. Для малых \(\theta\) можно записать \(V(\theta) = a(\omega)\,\theta^2 + b(\omega)\,\theta^4\). Найдите коэффициенты \(a(\omega)\) и \(b(\omega)\). (\(1{,}0\) балла)

B-6. Постройте типичные графики \(V(\theta)\) в зависимости от \(\theta\) для значений \(\omega/\omega_c\), чуть меньших \(1{,}0\) (например, \(0{,}9\)), и для большого \(\omega/\omega_c\), например \(5{,}0\). Требуются только качественные наброски, без подробных расчётов. (\(1{,}0\) балла)

B-7. Теория Ландау фазовых переходов второго рода позволяет разделить простые магнитные системы на две фазы. При температурах \(T\), больших критической температуры \(T_c\), система парамагнитна. При \(T<T_c\) система ферромагнитна, и намагниченность \(\mathcal M\) задаётся формулой \[\mathcal M(T) = \mathcal M_0\,(1-T/T_c)^{1/2}, \qquad T<T_c.\] Обозначим показатель \(1/2\) через \(\beta\). Сравните это поведение с задачей о бусинке. Какие аналоги \(\mathcal M\), \(T_c\), \(T/T_c\) в нашем случае? Чему равно эквивалентное значение \(\beta\) в нашем случае? (\(1{,}0\) балла)

B-8. Определите угловую частоту колебаний \(\Omega_0\) бусинки при отклонении от «равновесного» положения \(\theta_0\). Замечание: для малых колебаний \[\Omega_0 = \frac{1}{R}\sqrt{\frac{V''(\theta_0)}{m}}.\] (\(1{,}0\) балла)

B-9. Качественно начертите график зависимости \(\Omega_0\) от \(\omega\). (\(1{,}0\) балла)

Для пунктов C-1–C-2 считайте \(k \neq 0\).

C-1. Считайте \(f = 1\) и выразите \(k = \tan\alpha\). Условие для углов равновесия \(\theta_0\) можно записать в виде \[\left(\frac{\omega}{\omega_c}\right)^2 = \frac{\tan(x)}{\sin(y)}.\] Найдите \(x\) и \(y\). (\(1{,}0\) балла)

C-2. Дано, что \(f=1\) и \(k = 0{,}05\). Найдите углы равновесия \(\theta_0\), если они существуют, для следующих случаев:

  1. \(\omega/\omega_c = 0{,}50\);

  2. \(\omega/\omega_c = 0{,}70\).

(\(0{,}5\) балла)

▶ПодсказкаСкрыть
Подсказка

Для B-пунктов достаточно записать тангенциальную компоненту уравнения движения при \(\ddot\theta = 0\) и найти нули правой части. Потенциал получается интегрированием \(F_\theta\); разложение по малым \(\theta\) дает коэффициенты \(a\) и \(b\).

▶ОтветСкрыть
Ответ

A-1) \(m R\dot\theta^2 = N - mg\cos\theta - mR\omega^2\sin^2\theta\);  \(mR\ddot\theta = mR\omega^2\sin\theta\cos\theta - mg\sin\theta - \operatorname{sgn}(\dot\theta)\,kN\);  B-1) \(\theta_0 = 0\) при любом \(\omega\); при \(\omega>\omega_c\) также \(\cos\theta_0 = \omega_c^2/\omega^2\), \(\theta_0 = \pm\arccos(\omega_c^2/\omega^2)\);  B-2) при \(\omega<\omega_c\) \(\theta_0=0\); при \(\omega>\omega_c\) появляются две ветви \(\pm\theta_0\), стремящиеся к \(\pm\pi/2\) при \(\omega\to\infty\);  B-3) \(N = mg\) при \(\omega<\omega_c\), \(N = m\omega^2R\) при \(\omega>\omega_c\);  B-4) \(P = mgR\), \(Q = -mgR\), \(S = -\dfrac{\omega^2}{2\omega_c^2}\,mgR\);  B-5) \(a(\omega) = \dfrac{mgR}{2}\left(1-\dfrac{\omega^2}{\omega_c^2}\right)\), \(b(\omega) = \dfrac{mgR}{6}\left(\dfrac{\omega^2}{\omega_c^2}-\dfrac14\right)\);  B-6) см. графики в решении;  B-7) \(\mathcal M \to \theta_0\), \(T_c \to 1/\omega_c^2\), \(T/T_c \to \omega_c^2/\omega^2\), \(\beta = 1/2\);  B-8) \(\Omega_0 = \sqrt{\omega_c^2-\omega^2}\) при \(\omega<\omega_c\), \(\Omega_0 = \omega\sqrt{1-\omega_c^4/\omega^4}\) при \(\omega>\omega_c\);  B-9) см. график в решении;  C-1) \(x = \theta_0+\alpha\), \(y = \theta_0\);  C-2) (1) \(\theta_0 \approx -0{,}07\) рад; (2) \(\theta_0 \approx -0{,}1\) рад.

▶РешениеСкрыть
Решение

A-1. Радиальная компонента: \[mR\dot\theta^2 = N - mg\cos\theta - mR\sin^2\theta\,\omega^2 \qquad (1)\] Касательная компонента: \[mR\ddot\theta = mR\sin\theta\cos\theta\,\omega^2 - mg\sin\theta - \operatorname{sgn}(\dot\theta)\,kN \qquad (2)\] или, в обозначениях \(f\): \[mR\ddot\theta = mR\sin\theta\cos\theta\,\omega^2 - mg\sin\theta - fkN, \quad f=1. \qquad (3)\]

B-1. Полагаем \(k=0\) в касательной компоненте силы: \[mR\ddot\theta = mR\sin\theta\cos\theta\,\omega^2 - mg\sin\theta \qquad (4)\] Для равновесия \(\ddot\theta_0 = 0\). Тогда \(\theta_0 = 0\) — угол равновесия при всех \(\omega\). Остальные значения находятся из \[\cos\theta_0 = \frac{g}{\omega^2R} = \frac{\omega_c^2}{\omega^2} \qquad (5)\] \[\theta_0 = \pm\left|\cos^{-1}\frac{\omega_c^2}{\omega^2}\right| \qquad (6)\] со значениями \(\theta_0\) в интервале от \(-\pi/2\) до \(\pi/2\). Знак \(\pm\) означает два равноправных положения: бусинка может подняться с любой стороны от оси. При \(\omega<\omega_c\) из (5) следует \(\cos\theta_0>1\), что физически невозможно; поэтому \(\theta_0=0\) при \(\omega<\omega_c\).

B-2. Эскиз \(\theta_0\) как функции \(\omega/\omega_c\) приведён на рисунке ниже: при \(\omega<\omega_c\) \(\theta_0=0\); при \(\omega>\omega_c\) есть две ветви \(\pm\theta_0\).

Рисунок

B-3. Эскиз величины нормальной реакции в устойчивом равновесии приведён ниже. При \(\omega<\omega_c\) имеем \(N = mg\); при \(\omega>\omega_c\) имеем \(N = m\omega^2R\) (это следует из (1) с учётом (5)).

Рисунок

B-4. Интегрируя (4) с \(V(0)=0\): \[-R\int_0^\theta F_\theta\,d\theta = \int_0^V dV = V.\] Подставляя \(F_\theta\) из (4), левая часть равна \[-R\int_0^\theta F_\theta\,d\theta = -\frac{m\omega^2R^2}{2}\int_0^\theta\sin2\theta\,d\theta + mgR\int_0^\theta\sin\theta\,d\theta = \frac{m\omega^2R^2}{4}\big(\cos2\theta-1\big) - mgR\big(\cos\theta-1\big),\] а с учётом \(\cos2\theta-1 = -2\sin^2\theta\) и \(\omega_c^2 = g/R\) получаем, откуда \[V(\theta) = mgR\left[(1-\cos\theta) - \frac{\omega^2}{2\omega_c^2}\sin^2\theta\right]. \qquad (9)\] Отсюда \(P = mgR\), \(Q = -mgR\), \(S = -\dfrac{\omega^2}{2\omega_c^2}\,mgR\). Тот же результат получается дифференцированием \(V = P+Q\cos\theta+S\sin^2\theta\) и сравнением с \(F_\theta\).

B-5. Используя разложения, сохраняя члены до \(\theta^4\): \[V(\theta) \approx mgR\left[\frac{\theta^2}{2}-\frac{\theta^4}{24} - \frac{\omega^2}{2\omega_c^2}\left(\theta-\frac{\theta^3}{6}\right)^2\right] = \frac{mgR}{2}\left(1-\frac{\omega^2}{\omega_c^2}\right)\theta^2 + \frac{mgR}{6}\left(\frac{\omega^2}{\omega_c^2}-\frac14\right)\theta^4.\] Отсюда \(a(\omega)\) и \(b(\omega)\) как в ответе. B-6. Качественные графики \(V(\theta)\) приведены ниже. При \(\omega<\omega_c\) потенциал имеет один минимум при \(\theta=0\) (красная кривая). При \(\omega>\omega_c\) в точке \(\theta=0\) — максимум, а минимумы — при \(\pm\theta_0\) (синяя кривая). Ниже приведены графики при \(\theta\in[-\pi/2,\pi/2]\) и \(\theta\in[-\pi,\pi]\). На рисунке: \(\omega<\omega_c\) — красная кривая, \(\omega>\omega_c\) — синяя кривая; \(-\theta_0\), \(\theta_0\) — положения минимумов.

Рисунок
Рисунок

B-7. Решение 1. При \(\omega\to\omega_c^+\) \(\theta_0\) близко к нулю. Разлагая косинус в (5), получаем \(1-\dfrac{\theta_0^2}{2} = \dfrac{\omega_c^2}{\omega^2}\), откуда \[\theta_0 = \pm\sqrt2\left[1-\frac{\omega_c^2}{\omega^2}\right]^{1/2}.\] При \(\omega\to\infty\) имеем \(\theta_0\to\pm\pi/2\). Аналогия с фазовым переходом: роль \(\mathcal M\) играет \(\theta_0\), а температура обратно пропорциональна \(\omega\). Итого: \[\mathcal M \to \theta_0, \qquad T_c \to \frac{1}{\omega_c^2}, \qquad \frac{T}{T_c} \to \frac{\omega_c^2}{\omega^2}.\] Эквивалентное значение \(\beta\) для бусинки равно \(1/2\).

Решение 2. При \(\omega>\omega_c\) \(\cos\theta_0 = \omega_c^2/\omega^2\). Используя \(\sin^2\theta_0 = 1-\cos^2\theta_0\), получаем \(\sin\theta_0 = \left(1-\omega_c^4/\omega^4\right)^{1/2}\). Аналогия: \(\mathcal M \to \sin\theta_0\) (или \(\theta_0\) при малом угле), \(T_c \to 1/\omega_c^4\), \(T/T_c \to \omega_c^4/\omega^4\). Эквивалентное значение \(\beta\) равно \(1/2\).

Замечание: показатель затухания намагниченности \(1/2\) как в нашем случае, так и в теории Ландау, но в эксперименте и в более сложных теориях показатель равен \(1/3\).

B-8. Частота колебаний бусинки около «равновесного» положения \(\theta_0\): \[\Omega_0 = \frac{1}{R}\sqrt{\frac{V''(\theta_0)}{m}}.\] Вторая производная потенциала: \[V''(\theta) = mgR\cos\theta\left[1-\frac{\omega^2}{\omega_c^2}\cos\theta\right] + mgR\frac{\omega^2}{\omega_c^2}\sin^2\theta.\] Для \(\theta_0 = \pm\cos^{-1}(\omega_c^2/\omega^2)\) при \(\omega>\omega_c\): \[V''(\theta_0) = mgR\,\frac{\omega^2}{\omega_c^2}\left(1-\frac{\omega_c^4}{\omega^4}\right) > 0.\] Для \(\omega<\omega_c\) \(\theta_0 = 0\), и из последнего соотношения получаем \[\Omega_0 = \left(\omega_c^2-\omega^2\right)^{1/2}.\] Аналогично для \(\omega>\omega_c\): \[\Omega_0 = \omega\left(1-\frac{\omega_c^4}{\omega^4}\right)^{1/2}.\]

B-9. График \(\Omega_0\) в зависимости от \(\omega\) приведён ниже. При \(\omega<\omega_c\) (красная кривая) \(\Omega_0\) убывает от \(\omega_c\) до \(0\); при \(\omega>\omega_c\) (синяя кривая) \(\Omega_0\) растёт от \(0\) и при больших \(\omega\) стремится к прямой \(\Omega_0=\omega\) (пунктир).

Рисунок

C-1. Подставляем выражение для нормальной реакции (1) в угловое уравнение (3) и получаем \[mR\ddot\theta = mR\sin\theta\cos\theta\,\omega^2 - mg\sin\theta - fk\left(mg\cos\theta + mR\sin^2\theta\,\omega^2 + mR\dot\theta^2\right).\] Используя \(\omega_c^2 = g/R\) и перегруппировывая члены, имеем \[\ddot\theta = \omega_c^2\left[\left(\sin\theta\right)\left(\cos\theta-fk\sin\theta\right)\left(\frac{\omega}{\omega_c}\right)^2 - \sin\theta - fk\cos\theta - fk\left(\frac{\dot\theta}{\omega_c}\right)^2\right].\] В равновесии \(\dot\theta=0\), \(\ddot\theta=0\), а \(f=\operatorname{sgn}(\dot\theta)=\pm1\) зависит от того, как было достигнуто равновесие. Для \(\theta_0\in(-\pi/2,\pi/2)\) получаем \[\sin\theta_0\left(\cos\theta_0-fk\sin\theta_0\right)\left(\frac{\omega}{\omega_c}\right)^2 = \sin\theta_0 + fk\cos\theta_0.\] При \(f=1\) и \(k=\tan\alpha\) это записывается в виде \[\left(\frac{\omega}{\omega_c}\right)^2 = \frac{\sin\theta_0+\tan\alpha\cos\theta_0}{\sin\theta_0\left(\cos\theta_0-\tan\alpha\sin\theta_0\right)} = \frac{\tan(\theta_0+\alpha)}{\sin\theta_0}.\] Таким образом, \(x = \theta_0+\alpha\), \(y = \theta_0\).

C-2. Так как \(k = 0{,}05\) мало, \(\alpha\approx k\), \(\sin\theta_0\approx\theta_0\), \(\tan(\theta_0+\alpha)\approx\theta_0+\alpha\), поэтому \[\left(\frac{\omega}{\omega_c}\right)^2 \approx 1 + \frac{k}{\theta_0}.\] Получаем:

  1. при \(\omega/\omega_c = 0{,}50\): \(\theta_0 \approx -0{,}07\) рад;

  2. при \(\omega/\omega_c = 0{,}70\): \(\theta_0 \approx -0{,}1\) рад.

Кривая перестаёт быть симметричной относительно нуля.

зад 3 Максвелл, Рэлей и Эверест

Волновая оптика и электромагнитные волны › Электромагнитные волныЭлектромагнитные волны, излучениеВынужденные колебания, резонансОценки, размерность и подобиесложность 7/10

Авторы задачи: Амитабх Вирмани (CMI, Ченнаи) и А. С. Бияни (в отставке, Государственный колледж естественных наук имени Нагарджуны, Райпур). Вклад Академического комитета, Группы академического развития и Международного жюри признателен.

Лорд Рэлей посетил Дарджилинг (Индия) в 1897 году. Глядя на Эверест с расстояния 170 км, он вспомнил вопрос Джеймса Максвелла, заданный 26 годами ранее, об ослаблении света в воздухе и видимости далёких вершин. Через два года, в 1899 году, лорд Рэлей опубликовал известную научную работу, в которой исследовал эту задачу. В данной задаче мы попытаемся восстановить, хотя бы частично, современную версию рассуждений лорда Рэлея.

Колебания электронного облака. Моделируем типичную нейтральную молекулу воздуха неподвижным положительным зарядом \(q\), окружённым сферическим равномерно заряженным облаком массы \(m\), радиуса \(r\) и заряда \(-q\). Собственная угловая частота колебаний молекулы равна \(\omega_0\). На молекулу падает свет, и отрицательное облако колеблется, сохраняя сферическую форму, с угловой частотой \(\omega\): \[y = y_0\cos(\omega t), \qquad (1)\] под действием электрического поля световой волны \[\vec E(t) = E_0\cos(\omega t)\,\hat\jmath. \qquad (2)\] Здесь \(y\) — расстояние между неподвижным положительным зарядом и центром отрицательного облака молекулы.

A-1. Запишите уравнение движения для \(y\). Примите, что вклад электрического поля в модуль ускорения равен \(E(t)\,q/m\). (\(0{,}5\) балла)

A-2. На основании приведённой информации решите уравнение для \(y\). Найдите амплитуду. (\(0{,}5\) балла)

A-3. Найдите модуль \(p(t)\) дипольного момента молекулы воздуха как функцию времени при \(\omega\ll\omega_0\). (\(0{,}5\) балла)

A-4. Получите выражение для \(\omega_0\) через \(q\), \(m\) и \(r\). (\(0{,}5\) балла)

Излучаемая мощность. Синусоидально изменяющийся во времени диполь излучает электромагнитное излучение. Мощность излучения зависит от амплитуды дипольного момента \(p_0 = qy_0\), электрической постоянной вакуума \(\epsilon_0\), скорости света \(c\) и частоты колебаний \(\omega\).

B-1. Используя размерный анализ, выразите среднюю излучаемую мощность \(s\) через эти величины. (\(1\) балл)

B-2. Примите коэффициент пропорциональности равным \(1/(12\pi)\). Выразите мощность излучения \(s\) через \(E_0\), \(\omega_0\), \(\omega\) и связанные величины. Считайте \(\omega\ll\omega_0\). (\(0{,}2\) балла)

Затухание интенсивности \(I(x)\). Напомним, что интенсивность электромагнитной волны равна \[I = \frac12 c\,\epsilon_0E_0^2. \qquad (3)\] Интенсивность уменьшается вдоль пути света, потому что в единице объёма теряется мощность \[S = n_0\,s, \qquad (4)\] где \(n_0\) — число молекул в единице объёма.

C-1. Составьте дифференциальное уравнение для интенсивности \(I(x)\) как функции расстояния \(x\). (\(1\) балл)

C-2. Получите выражение для интенсивности \(I(x)\) как функции \(x\) через характерный масштаб длины \(L\), связанный с падением интенсивности. Начальная интенсивность равна \(I_0\). (\(0{,}5\) балла)

C-3. Считайте \(m\) массой электрона (обычно один электрон составляет облако заряда) и примите \[n_0 = 2{,}54\times10^{25}\ \text{м}^{-3}, \qquad \omega_0 = 1{,}25\times10^{16}\ \text{рад/с}, \qquad \omega = 3{,}25\times10^{15}\ \text{рад/с}.\] Найдите численное значение \(L\) в километрах. (\(0{,}3\) балла)

Высота гор, видимая наблюдателем. На рисунке точка \(P\) обозначает Дарджилинг, горную станцию в восточных Гималаях на высоте \(h = 2042\) м над уровнем моря. Линия \(BS\) обозначает Эверест, находящийся на расстоянии \(d = 170\) км от Дарджилинга и имеющий высоту \(H = 8848\) м. Другой пик, Канченджанга (на рисунке не показан), находится на расстоянии \(75\) км от Дарджилинга и имеет высоту \(8586\) м.

Рисунок

Рисунок 1 — Большой круг, на котором лежат гора \(BS\) высотой \(H\) и наблюдатель \(P\) на высоте \(h\). Рисунок не в масштабе. На рисунке: \(O\) — центр Земли; \(R\) — радиус Земли; \(\theta\) — угол; \(d\) — расстояние; \(h\) — высота; \(H\) — высота горы.

D-1. Получите выражение и численные значения вертикальной высоты \(H'\) этих гор, видимых наблюдателем из Дарджилинга, через перечисленные величины. Считайте, что наблюдатель не видит ниже местного горизонта. Нарисуйте подходящую схему. Радиус Земли \(R\) примите равным \(6378\) км. (\(2\) балла)

E-1. Примите интенсивность света от Канченджанги в Дарджилинге за опорное значение. Какова интенсивность света от Эвереста относительно Канченджанги, если смотреть из Дарджилинга? В задаче не учитываем изменение плотности молекул воздуха с высотой. Если интенсивность составляет \(5\,\%\) или более от опорного значения, гора считается видимой. Будет ли виден Эверест из Дарджилинга? (\(1\) балл)

Длина ослабления \(L_p\) из-за загрязнения аэрозолями. Выше мы вычислили характерную длину \(L\), связанную с падением интенсивности из-за рассеяния на молекулах воздуха. Теперь нас интересует характерная длина \(L_p\), связанная с падением интенсивности из-за рассеяния на аэрозольных частицах (загрязнении). \(L_p\) зависит от концентрации частиц \(n\) и площади поперечного сечения \(\pi r^2\) аэрозольной частицы радиуса \(r\). Получите эту зависимость, используя физические соображения и размерный анализ. Безразмерную константу примите равной \(1/8\). При умеренном загрязнении средняя плотность аэрозолей в Дарджилинге равна \(\rho_p = 5\) мкг/м\(^3\), а их средний радиус равен \(500\) нанометров. Какова длина \(L_p\)? Плотность отдельной аэрозольной частицы \(\rho = 3\) г/см\(^3\). Напомним, \(1\) мкг \(=10^{-9}\) кг, \(1\) нм \(=10^{-9}\) м.

F-1. Найдите \(L_p\). (\(1\) балл)

G-1. Оцените относительную интенсивность Канченджанги и Эвереста при указанном уровне загрязнения относительно опорного значения. Какие из этих пиков, если вообще какие-либо, будут видны из Дарджилинга? Считайте, что загрязнение однородно по всему пути света. (\(1\) балл)

▶ПодсказкаСкрыть
Подсказка

Ослабление складывается из двух независимых вкладов: рассеяния на молекулах воздуха (длина \(L\)) и рассеяния на аэрозолях (длина \(L_p\)). В экспоненте они входят как сумма \(d/L + d/L_p\), а геометрия учитывается через кривизну Земли.

▶ОтветСкрыть
Ответ

A-1) \(\ddot y = -\omega_0^2y - \dfrac{qE_0}{m}\cos\omega t\);  A-2) \(y_0 = \dfrac{qE_0/m}{\omega^2-\omega_0^2}\);  A-3) \(p(t) = qy(t) \approx \dfrac{q^2E_0}{m\omega_0^2}\cos\omega t\);  A-4) \(\omega_0 = \dfrac{q}{\sqrt{4\pi\epsilon_0 m r^3}}\);  B-1) \(s = k\,\dfrac{p_0^2\omega^4}{\epsilon_0c^3}\);  B-2) \(s = \dfrac{1}{12\pi}\,\dfrac{q^4E_0^2}{m^2\epsilon_0c^3}\,\dfrac{\omega^4}{\omega_0^4}\);  C-1) \(-\dfrac{dI}{dx} = n_0 s\);  C-2) \(I(x) = I_0e^{-x/L}\), \(L = \dfrac{6\pi\epsilon_0^2m^2c^4}{n_0q^4}\left(\dfrac{\omega_0}{\omega}\right)^4\);  C-3) \(L\approx130\) км;  D-1) \(H' = R+h-\dfrac{(R+h)^2}{R+H}\sec\theta \approx R+h-\dfrac{(R+h)^2}{R+H}\left(1+\dfrac{d^2}{2R^2}\right)\); для Канченджанги \(H'\approx 6096\) м, для Эвереста \(H'\approx4534\) м;  E-1) \(I_{\text{Эверест}}/I_{\text{Канченджанга}}\approx0{,}093\); Эверест виден;  F-1) \(L_p = \dfrac{1}{8n\pi r^2} = \dfrac{r\rho}{6\rho_p}\approx50\) км;  G-1) \(I_K/I_{\text{ref}}\approx0{,}22\) — Канченджанга видна; \(I_E/I_{\text{ref}}\approx0{,}003\) — Эверест не виден.

▶РешениеСкрыть
Решение

A-1. \(E(t)\) — электрическое поле в месте молекулы. Уравнение движения заряда без \(\vec E(t)\) есть \(\ddot y = -\omega_0^2y\); при вынужденных колебаниях \[\ddot y = -\omega_0^2y - \frac{qE_0}{m}\cos\omega t. \qquad (2)\]

A-2. Подставляя \(y = y_0\cos\omega t\) в (2), получаем \[-\omega^2y_0 = -\omega_0^2y_0 - \frac{qE_0}{m}, \qquad (3)\] откуда амплитуда колебаний \[y_0 = \frac{qE_0/m}{\omega^2-\omega_0^2}. \qquad (4)\]

A-3. Поскольку \(y\) — расстояние между зарядами, модуль дипольного момента равен \[p(t) = qy(t) \approx \frac{q^2E_0}{m\omega_0^2}\cos\omega t. \qquad (5)\]

A-4. Атом моделируем как неподвижный положительный точечный заряд \(q\), окружённый сферическим отрицательным облаком полного заряда \(-q\), радиуса \(r\) и массы \(m\). При смещении облака на малое расстояние \(y\) электростатическая сила со стороны центрального положительного заряда равна \[\vec F_{\text{эл}} = m\ddot y\,\hat\jmath = -\frac{q^2}{4\pi\epsilon_0r^3}\,y\,\hat\jmath = -m\omega_0^2y\,\hat\jmath, \qquad (6)\] откуда собственная частота колебаний \[\omega_0 = \frac{q}{\sqrt{4\pi\epsilon_0mr^3}}. \qquad (7)\]

Рисунок

Рисунок к решению A-4 — модель атома: облако заряда \(-q\) радиуса \(r\) смещено на \(y\) относительно заряда \(+q\). На рисунке: \(-q\), \(+q\) — заряды; \(r\) — радиус облака; \(y\) — смещение; \(\hat\jmath\), \(\hat\imath\) — единичные векторы.

B-1. Размерности: \([s] = \text{кг}\cdot\text{м}^2\cdot\text{с}^{-3}\), \([p_0] = \text{Кл}\cdot\text{м}\), \([\omega] = \text{с}^{-1}\), \([\epsilon_0] = \text{Кл}^2\cdot\text{Н}^{-1}\cdot\text{м}^{-2} = \text{Кл}^2\cdot\text{кг}^{-1}\cdot\text{м}^{-3}\cdot\text{с}^2\), \([c] = \text{м}\cdot\text{с}^{-1}\). Пусть \(s = p_0^{\alpha}\omega^{\beta}\epsilon_0^{\gamma}c^{\delta}\). Приравнивая показатели, получаем \(\alpha=2\), \(\beta=4\), \(\gamma=-1\), \(\delta=-3\), то есть \[s = k\,\frac{p_0^2\omega^4}{\epsilon_0c^3}. \qquad (8)\]

B-2. Подставляя \(p_0 = qy_0\) и \(y_0 \approx q E_0/(m\omega_0^2)\) при \(\omega\ll\omega_0\), получаем \[s = \frac{1}{12\pi}\,\frac{q^2y_0^2\omega^4}{\epsilon_0c^3} = \frac{1}{12\pi}\,\frac{q^4E_0^2}{m^2\epsilon_0c^3}\,\frac{\omega^4}{\omega_0^4}. \qquad (9)\]

C-1. Рассмотрим горизонтальный столб атмосферы сечения \(A\) и длины \(\Delta x\). Пусть падающая интенсивность равна \(I(x)\), а прошедшая — \(I(x+\Delta x)\). Падение интенсивности обусловлено рассеянием света на молекулах воздуха. Если концентрация молекул равна \(n_0\), то полная излучаемая мощность в единице объёма равна \(n_0 s\). Тогда \[I(x)A - I(x+\Delta x)A = n_0 s\,(A\Delta x), \qquad (10)\] откуда \[-\frac{dI}{dx} = n_0 s. \qquad (11)\]

C-2. Так как \(s\propto E_0^2\) и \(I\propto E_0^2\), имеем \[-\frac{dI}{dx} = \frac{I}{L}, \qquad (12)\] где \[L = \frac{6\pi\epsilon_0^2m^2c^4}{n_0q^4}\left(\frac{\omega_0}{\omega}\right)^4. \qquad (13)\] Решение дифференциального уравнения: \(I(x) = I_0e^{-x/L}\).

C-3. Подставляя числа, получаем \(L\approx130\) км.

D-1. В треугольнике \(OPB'\) имеем \(OB' = OP\sec\theta = (R+h)\sec\theta\). Так как \(\angle OSP = \angle B'SS'\) и \(\angle SOP = \angle SB'S'\), треугольники \(OSP\) и \(B'SS'\) подобны. Отсюда \[\frac{B'S'}{OP} = \frac{B'S}{OS} = \frac{OS-OB'}{OS} = 1-\frac{OB'}{OS}.\] С учётом \(B'S' = H'\), \(OP = R+h\), \(OS = R+H\): \[\frac{H'}{R+h} = 1 - \frac{(R+h)\sec\theta}{R+H}, \qquad H' = R+h-\frac{(R+h)^2}{R+H}\sec\theta. \qquad (14)\] Так как \(\cos\theta\approx1-\theta^2/2\) и \(\theta = d/R\), получаем \[H' \simeq R+h-\frac{(R+h)^2}{R+H}\left(1+\frac{d^2}{2R^2}\right). \qquad (15)\] Численно: \(H' = 6096\) м для Канченджанги и \(H' = 4534\) м для Эвереста. Альтернативная запись: \(H' = H-h-\dfrac{d^2}{2R}\).

Рисунок

Рисунок к решению D-1 — большой круг с горой \(BS\) высотой \(H\), наблюдатель \(P\) на высоте \(h\), точка \(B'\) — видимая вершина, \(S'\) — проекция. Рисунок не в масштабе.

E-1. В точке наблюдения интенсивность, которая была бы без ослабления, равна \(I_0 = P_0/(4\pi d^2)\), так что \[I = \frac{P_0}{4\pi d^2}\exp\left[-\frac{d}{L}\right]. \qquad (16)\] Относительная интенсивность Эвереста по отношению к Канченджанге: \[\frac{I_{\text{Эверест}}}{I_{\text{Канченджанга}}} = \frac{d_{\text{Канч}}^2}{d_{\text{Эверест}}^2}\exp\left[-\frac{d_{\text{Эверест}}-d_{\text{Канч}}}{L}\right] = 0{,}093.\] Так как \(0{,}093 > 0{,}05\), Эверест виден из Дарджилинга.

F-1. Длина ослабления из-за аэрозолей \[L_p = \frac{1}{8n\pi r^2}. \qquad (17)\] Здесь \(n = \rho_p/m\) и \(m = \dfrac{4\pi}{3}r^3\rho\). Отсюда \[L_p = \frac{r\rho}{6\rho_p}.\] Подставляя \(r = 500\times10^{-9}\) м, \(\rho = 3\times10^{3}\) кг/м\(^3\), \(\rho_p = 5\times10^{-9}\) кг/м\(^3\), получаем \(L_p = 50\) км.

G-1. Новое выражение для ослабления интенсивности: \[I = \frac{P_0}{4\pi d^2}\exp\left[-\frac{d}{L}-\frac{d}{L_p}\right].\] Для Канченджанги \[\frac{I_K}{I_{\text{ref}}} = \exp\left[-\frac{d_K}{L_p}\right] = 0{,}22.\] Интенсивность снижается до \(22\,\%\) от опорного значения, поэтому Канченджанга будет видна из Дарджилинга. Для Эвереста \[\frac{I_E}{I_{\text{ref}}} = 0{,}093\cdot\exp\left[-\frac{d_E}{L_p}\right] = 0{,}093\times0{,}033 = 0{,}003.\] Интенсивность снижается до \(0{,}3\,\%\) от опорного значения, поэтому Эверест из Дарджилинга видно не будет.

Источник: Азиатская физическая олимпиада (APhO) 2022, теоретический тур — условия (1) · условия (2) · условия (3) · решения (1) · решения (2) · решения (3). Задачи перенабраны на botay.pro, решения — по официальным.