Эффект Ааронова-Бома и суперсимметрия (продолжение)

Это продолжение предыдущего поста. Вот его краткое содержание. Рассматривая гамильтониан частицы в поле соленоида \[ \mathrm{H}=\frac{\boldsymbol{\pi}^{2}}{2m}-\mu\sigma_{z}H_{z}, \] мы заметили, что при $g=4\mu m/e=2$ его можно записать в виде $\mathrm{H}=2Q_{1}^{2}=2Q_{2}^{2}$, где $Q_{1}=\boldsymbol{\sigma}\cdot\boldsymbol{\pi}/2\sqrt{m}$ и $Q_{2}=-i\sigma_{z}Q_{1}$ --- антикоммутирующие операторы.
Радиальное уравнение выглядит так \begin{equation} -P^{\prime\prime}\left(\rho\right)-\frac{1}{\rho}P^{\prime}\left(\rho\right)+\left(\frac{W\left(\rho\right)^{2}}{\rho^{2}}+\frac{\sigma}{\rho}W^{\prime}\left(\rho\right)\right)P\left(\rho\right)=k^{2}P\left(\rho\right)\label{eq:Pequation-2} \end{equation} Функция $W\left(\rho\right)$ зависит от распределения поля в соленоиде и обладает свойствами \[ W\left(0\right)=M,\quad W\left(\rho>R\right)=W=M-\Phi/\Phi_{0}\,. \] Благодаря суперсимметрии, нам удалось найти одно из решений этого уравнения в пренебрежении правой частью (т.е., при $\rho\ll1/k$): \begin{equation} P_{1}\left(\rho\right)=\exp\left[\sigma\int_{0}^{\rho}\frac{d\rho}{\rho}W\left(\rho\right)\right]\label{eq:sol} \end{equation} Тут следует заметить, что интеграл в экспоненте расходится на нижнем пределе если только $M\neq0$. Поэтому лучше писать, скажем, так \begin{equation} P_{1}\left(\rho\right)=\exp\left[\sigma\int_{R}^{\rho}\frac{d\rho}{\rho}W\left(\rho\right)\right]\label{eq:sol-2} \end{equation} Далее есть четыре различных случая \begin{align*} \mathrm{I}. & \sigma M\geq0\&\sigma W>0,\\ \mathrm{II}. & \sigma M\geq0\&\sigma W<0,\\ \mathrm{III}. & \sigma M<0\&\sigma W>0,\\ \mathrm{IV}. & \sigma M<0\&\sigma W<0. \end{align*} В случаях $\mathrm{I}$ и $\mathrm{II}$ функция $P_{1}\left(\rho\right)$ ведет себя при $\rho\to0$ как $\rho^{\left|M\right|}$, как и следует, конечно. При этом в случае $\mathrm{II}$ при $\rho>R$ мы получаем аномальное поведение $\rho^{-\left|W\right|}$, за которым мы, собственно, и охотились.
В случаях $\mathrm{III}$ и $\mathrm{IV}$ функция $P_{1}\left(\rho\right)$ ведет себя при $\rho\to0$ похабно --- как $\rho^{-\left|M\right|}$, что нас, конечно, не устраивает. Найдем второе решение, имеющее в нуле правильную асимптотику. Ищем его методом вариации постоянных \begin{equation} P_{2}\left(\rho\right)=C\left(\rho\right)P_{1}\left(\rho\right)\label{eq:sol-1} \end{equation}  Для функции $C\left(\rho\right)$ имеем уравнение \begin{equation} C^{\prime\prime}\left(\rho\right)+2C^{\prime}\left(\rho\right)\frac{\sigma}{\rho}W\left(\rho\right)+\frac{1}{\rho}C^{\prime}\left(\rho\right)=0\label{eq:Pequation-2-1} \end{equation} Считая выполненным условие $\sigma M<0$, получаем \begin{align*} C^{\prime}\left(\rho\right) & =\frac{1}{R}\exp\left[-\int_{R}^{\rho}\frac{d\rho}{\rho}\left(2\sigma W\left(\rho\right)+1\right)\right]=\frac{1}{\rho}\exp\left[-2\sigma\int_{R}^{\rho}\frac{d\rho}{\rho}W\left(\rho\right)\right]\\ C\left(\rho\right) & =\int_{0}^{\rho}\frac{d\rho}{\rho}\exp\left[-2\sigma\int_{R}^{\rho}\frac{d\rho}{\rho}W\left(\rho\right)\right] \end{align*} Легко определить зависимость коэффициента $C$ при $\rho\to0$: интеграл в показателе экспоненты равен $-2\sigma M\ln\rho+\mathrm{const}$, поэтому $C\left(\rho\right)\sim\rho^{-2\sigma M}=\rho^{2\left|M\right|}$. Соответственно, $P_{2}\left(\rho\right)\sim\rho^{\left|M\right|}$, как и следует.
Теперь определим поведение коэффициента $C$ при $\rho>R$: \begin{align*} C & =\int_{0}^{R}\frac{d\rho}{\rho}\exp\left[-2\sigma\int_{R}^{\rho}\frac{d\rho}{\rho}W\right]+\frac{1}{2\sigma W}\left(1-\left(\frac{\rho}{R}\right)^{-2\sigma W}\right)\\  & \sim\begin{cases} \rho^{0} & \text{при }R\ll\rho\ll1/k\text{ для случая }\mathrm{III}\\ \rho^{2\left|W\right|} & \text{при }R\ll\rho\ll1/k\text{ для случая }\mathrm{IV} \end{cases} \end{align*} Учитывая, что \[ P_{1}\left(\rho\right)\sim\rho^{\sigma W}=\begin{cases} \rho^{\left|W\right|} & \text{при }R\ll\rho\ll1/k\text{ для случая }\mathrm{III}\\ \rho^{-\left|W\right|} & \text{при }R\ll\rho\ll1/k\text{ для случая }\mathrm{IV} \end{cases}, \] получаем, что нужное нам решение $P_{2}\left(\rho\right)$ для обоих случаев $\mathrm{III}$ и $\mathrm{IV}$ при $R\ll\rho\ll1/k$ ведет себя как $\rho^{\left|W\right|}$.
Резюмируя, можно сказать, что аномальное поведение радиальной волновой функции в пределе $R\to0$ наблюдается только в случае $\mathrm{II}$, т.е., при $\sigma M\geq0\&\sigma W<0$. В случаях же $\mathrm{III}$ и $\mathrm{IV}$ состояние не зануляется действием операторов $Q_{1,2}$.

Эффект Ааронова-Бома и суперсимметрия

Рассмотрим нерелятивистскую частицу в поле соленоида. Гамильтониан Паули выглядит так:\[ \mathrm{H}=\frac{\boldsymbol{\pi}^{2}}{2m}-\mu\boldsymbol{\sigma}\cdot\mathbf{H}, \]Здесь $\boldsymbol{\pi}=\mathbf{p}-e\mathbf{A}$ --- удлиненный импульс. Поле точечного соленоида имеет вид\[ \mathbf{A}=\frac{\Phi}{2\pi}\frac{\mathbf{e}_{z}\times\mathbf{r}}{\rho^{2}}. \]Зависимость волновой функции от $z$ тривиальна и далее мы соответствующую часть гамильтониана не рассматриваем. В частности, мы считаем далее, что $\mathbf{p}=\left(p_{x},p_{y}\right)$. Переменные разделяются в координатах $\rho,\phi$ и, если искать решение в виде $\psi=e^{iM\phi}P\left(\rho\right)$, получим на функцию $P\left(\rho\right)$ уравнение\begin{equation} P^{\prime\prime}+\frac{1}{\rho}P^{\prime}+\left(k^{2}-\frac{\left(M-\Phi/\Phi_{0}\right)^{2}}{\rho^{2}}\right)P=0\label{eq:Pequation} \end{equation}Здесь $\Phi_{0}=\frac{2\pi\hbar c}{e}$ --- квант магнитного потока. Поскольку магнитное поле отсутствует во всей области $\rho>0$, это уравнение имеет одинаковый вид для частицы с любым спином. Будем считать, что $W=M-\Phi/\Phi_{0}$ не является целям числом. Решения уравнения (\ref{eq:Pequation}) очевидны --- функции Бесселя порядка $\pm\left|W\right|$. При $\rho\to0$ функция Бесселя с отрицательным индексом стремится к бесконечности, поэтому выбираем решение, имеющее приличную асимптотику в нуле: \[ P\left(\rho\right)\propto J_{\left|W\right|}\left(k\rho\right)\,. \]Здесь читателю предлагается подумать первый раз. Все правильно?

Разберемся. Понятно, что бесконечно тонкий соленоид означает предел соленоида конечного радиуса $R$ при $R\to0$. Или даже не предел, а просто ситуацию, когда $kR\ll1$. Рассмотрим эту ситуацию подробно. Для соленоида конечного радиуса имеем \[ \mathbf{A}=\frac{\Phi}{2\pi}\frac{\mathbf{e}_{z}\times\mathbf{r}}{\rho^{2}}f\left(\rho/R\right), \] где функция $f\left(x\right)$ стремится к нулю при $x\to0$ и равна единице при $x\geq 1$. Например, для случая однородного поля внутри соленоида эта функция равна $f\left(x\right)=\min\left(1,x^{2}\right)$. В общем случае магнитное поле внутри соленоида выражается через производную этой функции:\[ H\left(\rho\right)=\frac{\Phi f^{\prime}\left(\rho/R\right)}{2\pi\rho R} \] Соответственно, уравнение на $P\left(\rho\right)$ имеет вид\begin{equation} P^{\prime\prime}+\frac{1}{\rho}P^{\prime}+\left(k_{\perp}^{2}-\frac{\left(M-\frac{\Phi}{\Phi_{0}}f\left(\rho/R\right)\right)^{2}}{\rho^{2}}+\sigma\frac{g}{2}\frac{\Phi}{\Phi_{0}}\frac{f^{\prime}\left(\rho/R\right)}{\rho R}\right)P=0\label{eq:Pequation1} \end{equation}Здесь $\sigma=\pm1$ для спина вверх/вниз, $g=4m\mu/e$. Если частица бесспиновая, можно положить $g=0$. Поскольку $kR\ll1$, анализ проводится абсолютно аналогично анализу медленных частиц в § 132 3-его тома Ландау&Лифшица. На расстояниях больше радиуса соленоида, но много меньше $1/k$ уравнение имеет вид\begin{equation} P^{\prime\prime}+\frac{1}{\rho}P^{\prime}-\frac{W^{2}}{\rho^{2}}P=0\label{eq:Pequation2} \end{equation}Общее решение --- комбинация степеней, которую мы запишем для удобства как\begin{equation} c_{1}\left(\rho/R\right)^{\left|W\right|}+c_{2}\left(\rho/R\right)^{-\left|W\right|}\label{eq:rhollR} \end{equation}Коэффициенты $c_{1,2}$ не зависят от $k$ и определяются граничными условиями в нуле. Из соображений размерности сразу очевидно, что, с точностью до общей нормировки, они не зависят и от $R$. Волновая функция при произвольных $\rho>R$ является линейной комбинацией $J_{\pm\left|W\right|}$ с коэффициентами, которые, используя асимптотику функций Бесселя, легко выразить через $c_{1,2}$:\[ P\left(\rho\right)=c_{1}\Gamma\left(\left|W\right|+1\right)\left(kR/2\right)^{-\left|W\right|}J_{\left|W\right|}\left(k\rho\right)+c_{2}\Gamma\left(-\left|W\right|+1\right)\left(kR/2\right)^{\left|W\right|}J_{-\left|W\right|}\left(k\rho\right) \]Ну вот теперь в этой формуле мы можем перейти к пределу $R\to0$. В этом пределе первый член доминирует над вторым и, да, мы получаем\[ P\left(\rho\right)\propto J_{\left|W\right|}\left(k\rho\right)\,. \]Здесь читателю предлагается подумать второй раз. Все правильно?

Нет, не совсем. Это рассуждение годится, когда $c_{1}\neq0$. А вот если $g=2$ может получиться как раз $c_{1}=0$! Чтобы это доказать, обозначим $W\left(\rho\right)=M-\frac{\Phi}{\Phi_{0}}f\left(\rho/R\right)$. Тогда во всей области $\rho\ll 1/k$ мы имеем \begin{equation} P^{\prime\prime}+\frac{1}{\rho}P^{\prime}-\left(\frac{W^{2}}{\rho^{2}}+\frac{\sigma}{\rho}W^{\prime}\right)P=0\label{eq:Pequation3} \end{equation}Подстановкой легко проверить, что у этого уравнения есть решение\begin{equation} P\left(\rho\right)=\exp\left[\sigma\int_{0}^{\rho}\frac{d\rho}{\rho}W\left(\rho\right)\right]\label{eq:sol} \end{equation}Учитывая, что при $\rho>R$ величина $W\left(\rho\right)=W$ не меняется, получаем\[ P\left(\rho>R\right)=C\left(\frac{\rho}{R}\right)^{\sigma W}, \]где $C=\exp\left[\sigma\int_{0}^{R}\frac{d\rho}{\rho}W\left(\rho\right)\right]$. Таким образом, когда $\mathrm{sgn}W=-\sigma$, мы получаем чудо: для любой функции $f\left(x\right)$ коэффициент $c_{1}$ в (\ref{eq:rhollR}) строго зануляется. Конечно, хотелось бы понять, откуда такое чудо взялось. Вспоминаем, что при $g=2$ гамильтониан Паули можно записать в виде\[ \mathrm{H}=\frac{\left[\boldsymbol{\sigma}\cdot\boldsymbol{\pi}\right]^{2}}{2m}\,. \]Введем операторы \[ Q_{1}=\boldsymbol{\sigma}\cdot\boldsymbol{\pi}/2\sqrt{m},\quad Q_{2}=-i\sigma_{3}Q_{1}\,. \]Напомню, что жирные символами мы обозначаем двумерные вектора в плоскости $xy$. Легко сообразить, что\[ \left\{ Q_{i},Q_{j}\right\} =\delta_{ij}\mathrm{H}. \] А что это за соотношение? Правильно, это самый простейший пример суперсимметричной алгебры. Нам, правда, хотелось бы сразу увидеть суперсимметрию в радиальном гамильтониане. Но тут есть проблемы. Дело в том, что хотя сам гамильтониан коммутирует с $l_{z}=-i\partial_{\varphi}$, операторы $Q_{1,2}$ --- уже нет. Поэтому определить операторы сразу для радиального уравнения не удается. Тем не менее, выбрасывая $k^{2}$ из уравнения, мы, фактически, получаем уравнение на состояние с энергией, равной нулю. Поскольку гамильтониан выражается через квадрат эрмитова оператора, это состояние, если существует, является основным и может быть найдено из уравнения $Q_{1}\psi=0$ или $Q_{2}\psi=0$. Если мы будем искать решение в виде собственной функции операторов $l_{z}$ и $\sigma_{z}$, как $\psi=e^{iM\varphi}\phi_{\sigma}P\left(\rho\right)$, то как раз получим (\ref{eq:sol}). Итак, подытожим.
Почти для любого $\mu$ наша прескрипция выбирать в качестве решения радиального уравнения функцию $J_{\left|W\right|}\left(kr\right)$ является правильной. Исключение составляет случай, когда $\mu=\mu_{B}=e/2m$ (или $\mu=-\mu_{B}$). Причиной такого необычного исключения можно считать суперсимметрию гамильтониана для этого случая. Для этих случаев мы должны выбирать $J_{-\sigma\left|W\right|}\left(kr\right)$. Замечательно этот исключительный случай как раз соответствует нерелятивистскому электрону.

---

P.S. Ну и, да, здесь читателю тоже предлагается подумать, не обманул ли я где-нибудь его.
P.P.S. Исправил некоторые опечатки (ну и ошибки, чего уж там). Но это еще не конец истории, продолжение следует.

Разложение на простые дроби (профессиональное)

Сегодня вспомнил, общаясь с дочкой, как в школе учат раскладывать на простые дроби. Для правильных дробей с square-free знаменателями от одной переменной рецепт формулируется так. Записываем сумму простых дробей с неопределенными коэффициентами $A,\ B,\ C$: \begin{equation*}\frac{x^2-8x-3}{(x-{\color{red}1})(x-{\color{green}3})(x-{\color{blue}6})}=\frac{A}{x-{\color{red}1}}+\frac{B}{x-{\color{green}3}}+\frac{C}{x-{\color{blue}6}}\end{equation*} Затем приводим к общему знаменателю, собираем в числителе коэффициенты при разных степенях $x$ и получаем линейную систему на $A,\ B,\ C$: \begin{align*}
A+B+C&=1,\\-9 A-7 B-4 C&=-8,\\18 A+6 B+3 C&=-3.
\end{align*}Решая эту  систему, наконец, находим \begin{equation*}
\frac{x^2-8x-3}{(x-{\color{red}1})(x-{\color{green}3})(x-{\color{blue}6})}=-\frac{1}{x-{\color{red}1}}+\frac{3}{x-{\color{green}3}}-\frac{1}{x-{\color{blue}6}}
\end{equation*} И да, я, вроде, тоже так когда-то делал. Пока не знал ТФКП.
А теперь я делаю так:\begin{multline*}
\frac{x^2-8x-3}{(x-{\color{red}1})(x-{\color{green}3})(x-{\color{blue}6})}=\frac{{\color{red}1}^2-8\cdot{\color{red}1}-3}{(x-{\color{red}1})({\color{red}1}-{\color{green}3})({\color{red}1}-{\color{blue}6})}\\
+\frac{{\color{green}3}^2-8\cdot{\color{green}3}-3}{({\color{green}3}-{\color{red}1})(x-{\color{green}3})({\color{green}3}-{\color{blue}6})}
+\frac{{\color{blue}6}^2-8\cdot{\color{blue}6}-3}{({\color{blue}6}-{\color{red}1})({\color{blue}6}-{\color{green}3})(x-{\color{blue}6})}\\
 =-\frac{1}{x-{\color{red}1}}+\frac{3}{x-{\color{green}3}}-\frac{1}{x-{\color{blue}6}}
\end{multline*} Если в знаменателе есть степени, нужно, естественно, итерировать. Кстати, рецепт обобщается и на случай нескольких переменных. Научите ребенка — и он сможет при случае щегольнуть быстрым счетом в школе.

Про идеалы

Идеалов в реальном мире не бывает, но в математике они существуют еще как. Идеалы в  математике бывают разные, но здесь я расскажу про идеалы коммутативных колец.
Вообще, кольцо --- это множество, замкнутое относительно двух операций со свойствами сложения и умножения. Эти операции поэтому часто так и называют: сложение и умножение. Правда, коммутативность умножения не требуется, но нас интересует как раз кольцо с коммутативным умножением --- коммутативное кольцо. Ну еще хотят чтобы был нуль, и был противоположный элемент. Существование противоположного элемента позволяет определить в этом множестве и операцию вычитания по формуле $a-b\equiv a+(-b)$, где $-b$ --- элемент, противоположный $b$. Заметим, что существование единичного элемента и обратного элемента для ненулевого $a$ не требуется, и в этом принципиальное отличие кольца от поля.
Говоря неформально, в кольце деление даже на ненулевой элемент может быть невыполнимо. И, естественно, встает вопрос, а когда же деление все-таки выполнимо? Возьмем элемент кольца $a$ и спросим, какие элементы $b$ на него делятся, для каких $b$ разрешимо уравнение \begin{equation}b=ca.\label{eq1}\end{equation} Множество допустимых значений $b$ является примером идеала.

Возьмем самый простой пример коммутативного кольца — множество целых чисел $\mathbb{Z}$ с обычным сложением и умножением. Тогда, например, множество чисел, кратных числу $a$ является идеалом.

Можно расширить определение идеала так: пусть заданы несколько элементов кольца $a_1,\ldots,a_n$. Тогда идеал определим как множество значений $b$ (из кольца, конечно), для которых уравнение\begin{equation}b=c_1a_1+\ldots+c_na_n\label{eq2}\end{equation} разрешимо относительно $c_1,\ldots,c_n$ (в кольце, конечно). Про такие идеалы говорят, что они порождаются элементами $a_1,\ldots,a_n$. Для кольца целых чисел наше обобщение \eqref{eq2} не добавляет ничего нового. Чтобы это понять, заметим, что, во-первых, достаточно рассматривать идеалы, порожденные натуральными числами. Во-вторых, идеал, порожденный натуральными числами $a_1,\ldots,a_n$, совпадает с идеалом, порожденным  $\mathrm{GCD}(a_1,\ldots,a_n)$. Таким образом, в кольце целых чисел не только каждому натуральному числу соответствует идеал, но и наоборот, что открывает возможность переформулировать утверждения о натуральных числах в утверждения об идеалах. Вот, например, идеал, соответствующий простому числу (по понятным причинам называется простым идеалом), --- это собственный идеал (собственный $\equiv$ не совпадающий со всем кольцом), обладающий следующим свойством: для любого его элемента, имеющего вид произведения, по-крайней мере один из сомножителей также принадлежит идеалу. Теперь, когда у нас есть такое определение, мы можем его использовать для идеалов любых коммутативных колец, и, естественно, надеяться, что у этой и подобных конструкций есть простой смысл и большая польза.

Вот, например, рассмотрим полиномиальное кольцо ---  множество полиномов от $n$ переменных $x_1,\ldots,x_n$ с комплексными, скажем, коэффициентами. Полиномы, как и целые числа, можно складывать и умножать.  Обозначать это множество будем так: $\mathbb{C}[x_1,\ldots,x_n]$. В этом кольце наше обобщение \eqref{eq2} уже не сводится к \eqref{eq1}. Для кольца целых чисел каждый идеал у нас соответствовал натуральному числу, а для кольца полиномов идеалу можно придать гораздо более интересный смысл. А именно, идеалу, порожденному полиномами $p_1,\ldots,p_n$, сопоставим алгебраическое многообразие в $R^n$, заданное полиномиальными уравнениями \begin{align}p_1=0\nonumber\\\vdots\label{eq3}\\p_n=0\nonumber\end{align}Только вот это соответствие неоднозначное. В частности, если мы любой полином в системе \eqref{eq3} заменим на его степень, многообразие, очевидно не поменяется, а идеал, вообще говоря, поменяется. Чтобы соответствие было однозначным, нам нужно рассматривать не все идеалы, а то, что называется радикал-идеалы. Радикал-идеал $\mathcal{I}$ --- это идеал, для которого из условия $a^n\in \mathcal{I}$ следует $a\in \mathcal{I}$. Так вот соответствие между радикал-идеалами и алгебраическими многообразиями однозначно. Это означает, что, исследуя радикал-идеалы, можно исследовать алгебраические многообразия. Кстати, радикал-идеалы в кольце $\mathbb{Z}$ соответствуют натуральным числам, в разложении которых на простые нет степеней. Поэтому радикал-идеалы называются еще полупростыми. Очевидно, что простые идеалы являются также полупростыми, а обратное не всегда верно.

Пример приводимого многообразия,
заданного уравнением
$(x^2+y^2-4)(x^2-y)=0$

Чтобы понять, чему соответствуют простые идеалы, нам нужно ввести понятие неприводимых алгебраических многообразий. Неприводимые алгебраические многообразия ---  алгебраические многообразия, которые нельзя представить в виде объединения нескольких алгебраических многообразий. Так вот, неприводимые алгебраические многообразия играют роль простых чисел --- им соответствуют простые идеалы.

Красота идеальной математики на этом не кончается, но вот пост пора закруглять. Могу только порекомендовать для самообразования книжку Кокс, Литтл, О'Ши Идеалы, многообразия и алгоритмы. 

P.S. Все, можно статьи писать прямо в блоггере. Последний MathJax поддерживает нумерацию уравнений + ссылки через \ref и \eqref. Чтобы узнать как это включить, смотреть в коде поста скрипт с  src="...MathJax.js"

Булевы алгебры, проекторы и квантовые игры

Размышляя о проекторах, вот, что я придумал! Рассмотрим следующие функции\[n(x)=1-x,\]\[a(x,y)=xy,\]\[o(x,y)=x+y-xy\]На множестве {0,1} эти функции дают обычные логические not,and и or, то есть задают простейшую булеву алгебру.
То же самое можно сказать более вычурно: Функции $n,a,o$ задают булеву алгебру  на множестве корней уравнения $x^2-x=0$.
Кажется, что это утверждение не несет ничего нового, ведь уравнение как раз и имеет два корня --- $0$ и $1$. Однако тут мы вспоминаем, что это уравнение является также определяющим уравнением для проекторов. Легко проверить, что если мы возьмем множество коммутирующих проекторов, замкнутое относительно операций $n$ и $a$, то приведенные выше функции определяют булеву алгебру. Далее операции $n,a,o$ будем обозначать как в логике:$\lnot,\land,\lor$.
Кто-то скажет "что толку в такой булевой алгебре". Попробую объяснить. Вспомним, что в квантовой механике коммутирующие эрмитовы операторы соответствуют одновременно измеримым величинам. Более того, если, скажем, есть у нас проектор $P_0=|\psi_0\rangle\langle\psi_0|$ на какое-нибудь состояние, то вероятность обнаружения этого состояния в состоянии $\psi$ равна среднему от проектора \[\langle\psi|P_0|\psi\rangle\] Ага, теперь все становится понятно. Пусть у нас есть запутанное состояние двух спинов. Например, такое\[|\psi\rangle=\frac1{\sqrt{2}}(|\uparrow\rangle|\downarrow\rangle-|\downarrow\rangle|\uparrow\rangle).\]  И мы, скажем, хотим определить вероятность измерения, которое покажет, что первый спин торчит вдоль оси $z$ или второй спин торчит вдоль оси $x$ (или то и другое). Очень просто,эта вероятность равна:\[\langle\psi|P_1^z\lor P_2^x|\psi\rangle=\frac34,\] где $P_1^z=\frac12(1+\sigma_1^z)$, $P_2^x=\frac12(1+\sigma_2^x)$. То есть наше определение измерения напрямую транслировалось в ответ, даже задумываться не пришлось. Это, конечно, простейший пример, но идея понятна. Например, операторы, соответствующие особо запутанным измерениям можно, перед тем как усреднять, упрощать по правилам логики.

Заметим, между прочим, что определенные функции $n,a,o$ имеют простую классическую вероятностную интерпретацию: если $x$ и $y$ --- вероятности независимых событий, то $n,a,o$, соответственно, вероятности того, что первое не произойдет, что произойдут оба и что произойдет хотя бы одно. Результаты измерения спинов в запутанном состоянии, вообще говоря, не являются независимыми, но для проекторов действуют все те же простые правила.

Дальше воображение рисует кучу приложений, например, в теории квантовых игр или в квантовых компьютерах.

Может кто это уже придумал все, а? Уж больно красиво.

Клин

На прошлой неделе один обед протупил над следующим "парадоксом".

Рассмотрим матрицу \[
\Lambda_{\mathbf{p}}=\frac12\left(1+\frac{\varepsilon+\beta m}{p^2}\boldsymbol{\alpha}\mathbf{p}\right).
\] В голове не укладывались 4 утверждения, касающиеся этой матрицы:

1. $\Lambda_{\mathbf{p}}$ --- проектор, т.е. $\Lambda_{\mathbf{p}}^2=\Lambda_{\mathbf{p}}$. 
2.  $\Lambda_{\mathbf{p}}u_{\mathbf{p}}=u_{\mathbf{p}}$,  т.е., это проектор на пространство решений уравнения $[\varepsilon-\boldsymbol{\alpha}\mathbf{p}-\beta m]u_{\mathbf{p}}=0$.
3. $\Lambda_{\mathbf{p}}+\Lambda_{-\mathbf{p}}=1$, откуда следует, что $\Lambda_{\mathbf{p}}\Lambda_{-\mathbf{p}}=0$.
4. $u^{\dagger}_{\mathbf{p}}u_{\mathbf{-p}}\neq0$.

Только написав на доске понял в чем дело.

Docked Cells в Mathematica

Все, кто программировал в Математике, знают, что одна из самых полезных операций называется Quit Kernel. Удивительно, что для нее нет ни клавиатурного сокращения, ни кнопки. Нужно лезть в меню Evaluation>Quit Kernel>Local, а затем еще жать Ok в появившемся окошке с предупреждением. Есть немного более простой способ, конечно --- ввести Quit и нажать Shift-Enter. Но если ядро занято другими вычислениями, этот способ не работает. Есть и третий способ: создать кнопку командой  Button["Quit", Quit[]] и жать на нее в свое удовольствие. Проблема в том, что кнопочка эта постепенно уплывает из зоны видимости. Можно, конечно, сгенерировать палитру (palette) и для кого-то это будет самый удобный вариант. Понятно, что для настоящего художника в палитре должна быть не только черная краска, так что в палитру можно оформить несколько ваших любимых кнопок.
Однако, в том, что палитра хранится отдельно от файла есть как достоинства, так и недостатки. Достоинства, понятно, в том, что одну палитру можно использовать при работе в любом файле. Недостаток --- в том, что взяв с собой файл на другой компьютер, нужно не забыть взять и установить туда и палитру. Выходом из этого концептуального тупика может стать Docked Cell. По сути, это та же палитра, только определенного вида. Важно, что хранится она  прямо в файле. Вот пример создания Docked Cell-а с тремя кнопками:
SetOptions[InputNotebook[], DockedCells -> {Cell[BoxData[ToBoxes[
      Row[{Button["Quit", Quit[]],
        Button["Abort", FrontEndTokenExecute["FindEvaluatingCell"];
         FrontEndTokenExecute["EvaluatorAbort"]],
        Button["Evaluate", FrontEndTokenExecute["EvaluateCells"];
         SelectionMove[InputNotebook[], Next, Cell]]}]
      ]],
    Background -> RGBColor[0.9, 0.9, 0.9],
    CellMargins -> 0,
    CellFrame -> {{0, 0}, {1, 0}}]}]
Жмем Shift-Enter и вуаля. Кому-то может показаться, что вбивать такую команду в каждом новом файле утомительно. Согласен, но тут можно и палитру сделать:
CreatePalette[Button["Install tools",SetOptions[InputNotebook[], DockedCells -> {Cell[BoxData[ToBoxes[
      Row[{Button["Quit", Quit[]],
        Button["Abort", FrontEndTokenExecute["FindEvaluatingCell"];
         FrontEndTokenExecute["EvaluatorAbort"]],
        Button["Evaluate", FrontEndTokenExecute["EvaluateCells"];
         SelectionMove[InputNotebook[], Next, Cell]]}]
      ]],
    Background -> RGBColor[0.9, 0.9, 0.9],
    CellMargins -> 0,
    CellFrame -> {{0, 0}, {1, 0}}]}]]]
А еще лучше --- сделать такую палитру
CreatePalette[{Row[{Button["Quit", Quit[]],
    Button["Abort", FrontEndTokenExecute["FindEvaluatingCell"];
     FrontEndTokenExecute["EvaluatorAbort"]],
    Button["Evaluate", FrontEndTokenExecute["EvaluateCells"];
     SelectionMove[InputNotebook[], Next, Cell]]}],
  Button["embed tools",
   SetOptions[InputNotebook[],
    DockedCells -> Most@First@NotebookGet[ButtonNotebook[]]],
   ContentPadding -> False, FrameMargins -> 0]}]

Еще идефикс

Все, конечно, знают, что риманова поверхность функции \[\sqrt{P_{2g+2}(z)},\] где $P_{2g+2}(z)$ --- полином $2g+2$ порядка, в общем случае гомеоморфна сфере с $g$ ручками. В частности, риманова поверхность функции $f(z)=\sqrt{z^4-1}$ должна быть гомеоморфна тору. Я вот тоже это недавно узнал. Ну, я когда-то знал это раньше, но забыл, да, да. Полезно потренировать воображение и сообразить, как непрерывной деформацией превратить тор в дважды накрытую сферу Римана с четырьмя точками ветвления. Я такими вещами обычно занимаюсь, когда не могу заснуть.
Проблема в том, что когда картинка понятна, очень хочется ее нарисовать, и сопротивляться этому желанию у меня пока не получается. Иначе картинка будет висеть в бэкграунде и тормозить весь мыслительный процесс. Вот только в данном случае нужно не нарисовать, а анимировать. К счастью, Wolfram Mathematica заботится о нас, толкает свой CDF формат (как альтернативу PDF-у, что-ли?). Ладно, попробуем. Ниже результат.
Когда ползунок в самой правой позиции, имеем дважды накрытую сферу Римана, на которой однозначно определена наша функция $f(z)$. Четыре конца красных линий соответствуют корням нашего полинома $\pm1,\pm i$. Сами линии обозначают один из возможных способов выбрать разрезы.

---

---

P.S. Главное страдание --- выдумать подходящую функцию, интерполирующую между двумя крайними положениями. Первоначальная простая, как три копейки, мысль --- интерполировать как $(1-t)*тор+t*сфера2$ --- дает убогий результат.