За время $t$ атом $\rm ^3He$ проходит расстояние $vt$. Оценим нижнюю и верхнюю границу для $vt$ : $vt \in [5 \cdot 10^{-4}, 0.5] ~\text{нм}$. Круг радиусом $r$ пересекается с прямой длиной $vt$ если его центр лежит внутри фигуры площадью $\pi r^2 + 2rvt$. При этом нужно вычесть площадь $\pi r^2$, соответствующую центрам кругов, которые пересекают начальную точку.
Тогда вероятность $p$ того, что на пути атома встретилась нить $p=n \cdot 2 rvt \ll 1$. Значит их количество $F_1(t) = 2 rvt n N$
Рассмотрим площадь $\frac{N_0}{n}$ на которой равномерно распределены все $N_0$ нитей аэрогеля. Тогда вероятность $p_k(t)$ того, что за время $t$ атом столкнулся с $k$ нитями задается выражением: \[ p_k = \frac{N_0!}{(N_0-k)! \, k!} \left( \frac{2 rvt n}{N_0} \right)^k \left(1 - \frac{2 rvt n}{N_0} \right)^{N_0-k} = \frac{N_0!}{(N_0-k)! \, k!} \frac{\left( \frac{2 rvt n}{N_0} \right)^k}{\left(1 - \frac{2 rvt n}{N_0} \right)^{k}} \left(1 - \frac{2 rvt n}{N_0} \right)^{N_0} \] известным, как распределение Бернулли. Возьмем предел при $N_0 \to \infty$: \[ p_k(t) \to \frac{N_0^k}{\left( \frac{N_0}{2rvtn}-1 \right)^k} \frac{e^{-2rvtn}}{k!} = (2rvtn)^k \cdot \frac{e^{-2rvtn}}{k!}\] и получим известное выражение для распределения Пуассона, если ввести обозначение $2rvn = \lambda$: \[ p_k(t) = \frac{(\lambda t)^k}{k!} e^{-\lambda t} \] Среднее количество столкновений за время $t$ тогда задается выражением: \[ G(t) = \sum_1^\infty k p_k = \sum_1^\infty k \frac{ (\lambda t)^k}{k!} e^{-\lambda t} = \lambda t \, e^{-\lambda t} \sum_0^\infty \frac{(\lambda t)^{k}}{k!} = \lambda t \]
В терминах предыдущего пункта $F_2(t) = N p_0(t) = N e^{-\lambda t}$
Пусть атом налетает на нить справа на расстоянии $h$ от центра. Направление его скорости меняется на угол $\alpha=-\pi + 2 \arcsin \frac{h}{R}$ при $h>0$ и на угол $\alpha=\pi - 2 \arcsin \frac{h}{R}$ при $h<0$. Если атом налетает на нить слева, то при $h>0$ угол $\alpha = \pi -2 \arcsin \frac{h}{R} $, при $h<0$ угол $\alpha = -\pi + 2 \arcsin \frac{h}{R}$. Высота $h$ принимает равновероятно любое значение от $0$ до $R$, то есть вероятность $p_{\Delta h} = \frac{\Delta h}{R}$ - вероятность того, что высота $h$ лежит внутри какого-то отрезка длины $\Delta h$ вложенного в отрезок $[0, R]$. Изменению направления скорости на угол $\phi$ ($\phi > 0$) соответствует $h = R \sin \frac{\pi - \phi}{2} = \cos \frac{\phi}{2}$. Значит изменению направления скорости на угол $\phi + \Delta \phi$ соответствует $h - \Delta h = R \cos \frac{\phi + \Delta \phi}{2}$. Воспользуемся тем, что $\Delta \phi \ll \phi$: \[ \Delta h = R \sin \frac{\phi}{2} \cdot \frac{\Delta \phi}{2} \] Тогда вероятность $p$ того, что $\alpha \in [\phi, \phi + \Delta \phi]$: \[ p =\frac{1}{2} p_{\Delta h} = \frac{\Delta \phi}{4} \sin \frac{\phi}{2}.\] Заменой $\phi \to |\phi|$ рассматриваются случаи $\phi < 0$. В итоге $S(\phi,\Delta \phi) = \frac{N \Delta \phi}{4} \sin \frac{|\phi|}{2}$.
Из пункта А2 вероятность того, что за время $t$ произойдет один удар задается формулой $p_1(t) = 2rvn \cdot e^{-2rvn t}$. Говоря в других терминах, вероятность $p$ того, что пройдя расстояния $l$, атом ударился ровно один раз $p(l) = 2rln \cdot e^{-2rln}$. Это значит то, что вероятность $p(l,l+dl)$ того, что удар произошел в отрезке $[l, l + dl]$ равна $p(l,l+dl)=e^{-2rln} \cdot 2rn \, dl$. Тогда по определению среднего: \[ \langle \lambda \rangle = \int\limits_0^{\infty} l \, p(l,l+dl) = \frac{1}{2rn} \int\limits_0^\infty u e^{-u} du \] Интеграл $I=\int\limits_0^\infty u e^{-u} du$ возьмем по частям: \[ I = -\int\limits_0^\infty u d\left( e^{-u}\right) = -u e^{-u} \bigg|_0^\infty + \int\limits_0^\infty e^{-u} du = 1. \]
По определению среднего: \[ \langle \lambda^2 \rangle = \int\limits_0^{\infty} l^2 \, p(l,l+dl) = \frac{1}{4r^2n^2} \int\limits_0^\infty u^2 e^{-u} du. \] Полученный интеграл сводится к взятому в предыдущей части: \[ \int\limits_0^\infty u^2 e^{-u} \, du = -u^2 e^{-u} \bigg|_0^\infty + 2 \int\limits_0^\infty u e^{-u} du = 2\]
За время $t \gg \dfrac{\langle \lambda \rangle}{v}$ атом в среднем испытает $m=\dfrac{vt}{\langle \lambda \rangle}$ столкновений. То есть его положение $\vec{r}$ будет задаваться следующей суммой: \[ \vec{r} = \sum^m_{i=1} \vec{\lambda}_i\] где $\vec{\lambda}_i$ - векторы между точками соударения. В силу перестановочности операций суммирования и усреднения получим: \[ \langle \vec{r}^2 \rangle = \sum_{i=1}^m \langle \vec{\lambda}_i^2 \rangle + \sum_{i \neq j} \langle \vec{\lambda}_i \cdot \vec{\lambda}_j \rangle = m \langle \lambda^2 \rangle + 2 \sum_{i=1}^m \sum_{p=1}^{m-i} \langle \vec{\lambda}_i \cdot \vec{\lambda}_{i+p} \rangle,\] причем длины векторов $\vec{\lambda}_i$ и $\vec{\lambda}_{i+p}$ являются независимыми переменными. Таким образом нас интересуют средние вида $\langle \cos ( \vec{\lambda}_i \wedge \vec{\lambda}_{i+p} )\rangle$. Найдем их итерационно: \[ \begin{split} \langle \cos ( \vec{\lambda}_i \wedge \vec{\lambda}_{i+p} )\rangle &= \langle \cos ( \vec{\lambda}_i \wedge \vec{\lambda}_{i+p-1} + \vec{\lambda}_{i+p-1} \wedge \vec{\lambda}_{i+p} ) \rangle = \langle \cos ( \vec{\lambda}_i \wedge \vec{\lambda}_{i+p-1}) \cos(\vec{\lambda}_{i+p-1} \wedge \vec{\lambda}_{i+p} ) - \sin ( \vec{\lambda}_i \wedge \vec{\lambda}_{i+p-1}) \sin(\vec{\lambda}_{i+p-1} \wedge \vec{\lambda}_{i+p} ) \rangle =\\&= \langle \cos ( \vec{\lambda}_i \wedge \vec{\lambda}_{i+p-1}) \cos \phi \rangle - \langle \sin ( \vec{\lambda}_i \wedge \vec{\lambda}_{i+p-1}) \sin \phi \rangle \end{split} \] С точки зрения угла $\phi$ последнее столкновение никак не зависит от истории столкновений, поэтому средние можно расцепить. При этом \[ \langle \cos \phi \rangle = \int\limits_{-\pi}^{\pi} \cos u \frac{S(u, du)}{N} = 2 \int\limits_{0}^{\pi} \frac{d u}{4} \cos u \sin \frac{u}{2} = -2 \int\limits_0^{\pi} \cos^2 \frac{u}{2} \, d \left(\cos \frac{u}{2} \right) + \int\limits_0^{\pi} d\left( \cos \frac{u}{2} \right) = -\frac{1}{3} \] \[ \langle \sin \phi \rangle = \int\limits_{-\pi}^{\pi} \sin u \frac{S(u, du)}{N} = 0, \] значит \[ \langle \cos ( \vec{\lambda}_i \wedge \vec{\lambda}_{i+p} )\rangle = \left( -\frac{1}{3} \right) \langle \cos ( \vec{\lambda}_i \wedge \vec{\lambda}_{i+p-1}) \rangle = \left( -\frac{1}{3} \right)^p. \] Подставим это выражение в искомую сумму: \[ \langle \vec{r}^2 \rangle = \sum_{i=1}^m \langle \vec{\lambda}_i^2 \rangle + 2\sum_{i \neq j} \langle \vec{\lambda}_i \cdot \vec{\lambda}_j \rangle = m \langle \lambda^2 \rangle + 2 \langle \lambda \rangle^2 \sum_{i=1}^m (m-i) \cdot \left( -\frac{1}{3} \right)^i \simeq m \langle \lambda^2 \rangle + 2\langle \lambda \rangle^2 m \frac{-\frac{1}{3}}{1+\frac{1}{3}} = m \left( \langle \lambda^2 \rangle - \frac{1}{2} \langle \lambda \rangle^2 \right),\]