среда, 25 января 2012 г.

Скрипт для задач на построение


Тестовая версия скрипта для задач на построение с помощью циркуля и линейки. Инструменты:
point
ставит свободные точки и пересечения кривых. Чтобы поставить точку на кривой нужно два раза по ней кликнуть. Чтобы поставить точки пересечения двух кривых нужно кликнуть на одну (она выделится), а затем на вторую.
line
функционал линейки. Чтобы провести линию через две точки, нужно на них кликнуть по очереди.
circle
функционал циркуля. Чтобы провести окружность с центром в заданной точке заданного радиуса, нужно кликнуть сначала на точку центра, а затем на пару точек, определяющих радиус.
delete
стирательная резинка.

пятница, 20 января 2012 г.

Фактор Зоммерфельда-Гамова-Сахарова

Фактор Зоммерфельда-Гамова-Сахарова является, как известно, множителем в сечении или вероятности, учитывающим кулоновское взаимодействие конечных (или начальных) частиц при малых скоростях. Важность этого множителя обусловлена тем, что даже при малом взаимодействии он может существенно отличаться от единицы. Дело в том, что теория возмущений строится формально как степенной ряд по параметру, характеризующему силу взаимодействия. Например, в квантовой электродинамике (КЭД) этот параметр равен $\alpha=e^2/\hbar c\approx 1/137.036$ (Далее работаем в системе $\hbar=c=1$). Наивно можно полагать, что поправки на взаимодействие подавлены по этому параметру (в КЭД --- на два порядка). Однако есть определенный класс поправок, которые оказываются пропорциональны не просто $\alpha^n$, а $(\alpha/v)^n$, где $v$ --- относительная скорость частиц. Если частицы релятивистские, т.е. $v\sim 1$, такие поправки малы, как и поправки, не содержащие скорость в знаменателе, но когда $v\ll 1$ такие члены могут стопроцентно изменить величину амплитуды.
К счастью, поправки такого вида можно вычислить с помощью обычной нерелятивистской квантовой механики. Пусть для определенности пара медленных частиц появляется в конечном состоянии. Тогда амплитуда $M$ вероятности такого процесса антилинейна по волновой функции пары. Обозначая за $\psi(\mathbf{r})$ координатную часть волновой функции этой пары, зависящую от относительного положения частиц, мы можем записать\[ M=\int d\mathbf{r}\psi^*(\mathbf r)K(\mathbf{r}), \] где $K(\mathbf r)$ --- функция, зависящая от процесса. В пренебрежении взаимодействием $\psi(\mathbf{r})$, имеет вид \[ \psi_0(\mathbf{r})=e^{i\mathbf{p r}}. \] Кулоновское взаимодействие $-\alpha/r$, конечно, меняет волновую функцию на (см. ЛЛ III,136.9):\[ \psi_С(\mathbf r)=e^{{\pi\alpha\over 2v}}\Gamma\left(1+i{\alpha\over v}\right)e^{i\mathbf{pr}}F\left(-i{\alpha\over v},1, -i(pr+\mathbf{pr})\right). \] Таким образом, фактор Зоммерфельда-Гамова-Сахарова равен\[ S_{\mathrm{SGS}}=\left|\frac{\int d\mathbf{r}\psi_C^*(\mathbf r)K(\mathbf{r})}{\int d\mathbf{r}\psi_0^*(\mathbf r)K(\mathbf{r})}\right|^2 \] Польза от этой формулы появляется в том случае, когда у нас есть возможность посчитать этот фактор не из процесс-зависящего определения, а универсально. Это можно сделать, когда характерный размер спадания ядра $K(\mathbf r)$ мал по сравнению с характерным размером изменения волновой функции. Типичный пример --- рождение пары на пороге. Тогда ядро $K$ имеет размер порядка комптоновской длины волны рождающихся частиц $\hbar/mc$, что действительно мало по сравнению с $\hbar/mv$ (мы рассматриваем нерелятивизм). Тогда мы выносим плавно меняющуюся $\psi$ из-под интеграла и получаем \[ S_{\mathrm{SGS}}=\left|\frac{\psi_С^*(0)}{\psi_0^*(0)}\right|^2={2\pi\alpha/v\over 1-\exp(-2\pi\alpha/v)} \] На первый взгляд кажется, что это универсальное выражение справедливо всегда независимо от специфики функции $K(\mathbf{r})$, при условии, конечно, что ее размер мал по сравнению с $\hbar/mv$. Но на самом деле тут есть одна засада: все это годится, если $\int d\mathbf r K(\mathbf r)$ не равен нулю. Это не так, если в амплитуду дают вклад только ненулевые орбитальные моменты. В этом случае интеграл от $K(\mathbf r)$ зануляется за счет интегрирования по углам и мы получаем неопределенность типа 0/0. В этом случае нам нужно начать раскладывать $\psi(\mathbf r)$ и оставить первый неисчезающий после интегрирования член разложения. Можно вместо этого сразу использовать парциальные разложения\[ \psi_0(\mathbf{r})=\sum(2l+1) i^{l}P_{l}(\mathbf{n}_r\mathbf{n}_p)R_{pl}^0(r) \,,\quad \psi_C(\mathbf{r})=\sum(2l+1) i^{l}P_{l}(\mathbf{n}_r\mathbf{n}_p)R_{pl}^C(r) \]Тогда фактор Зоммерфельда-Гамова-Сахарова для углового момента $l$ выражается как предел отношения радиальных волновых функций\[ S_{\mathrm{SGS}}(l)=\lim_{r\to 0} \left|R_{pl}^C(r)\over R_{pl}^0(r)\right|^2. \]Радиальные волновые функции имеют не очень простой вид но нам нужны только разложения при маленьких $r$:\[ R_{pl}^{0}(r)\approx\frac{\left(pr\right)^{l}}{\left(2l+1\right)!!}=\frac{\left(2pr\right)^{l}l!}{\left(2l+1\right)!},\quad R_{pl}^C(r)=e^{-\frac{\pi\alpha}{2v}}\Gamma\left(l+1+i\frac{\alpha}{v}\right)\frac{(2pr)^{l}}{(2l+1)!} \]Получаем окончательно\[ S_{\mathrm{SGS}}(l)=e^{\frac{\pi\alpha}{v}}\left|\Gamma\left(l+1+i\frac{\alpha}{v}\right)\right|^2={2\pi\alpha/v\over 1-\exp(-2\pi\alpha/v)}\prod_{s=1}^{l}\left(1+{\alpha^2\over v^2s^2}\right) \]Заметим, что при очень маленьких $v$ фактор $S_{\mathrm{SGS}}(l)$ ведет себя как $v^{-2l-1}$. Если вспомнить множитель $(pr)^l$ в радиальной функции и множитель $p$ из фазового объема родившихся частиц, то можно сделать вывод о конечности порогового сечения для любого $l$ в случае притяжения.
P.S. И кстати, теперь я совершенно не понимаю, зачем нужно этот множитель обобщать на релятивистский случай.