Матричное представление светоделителя для численного расчета выходных данных на основе заданного количества фотонов (состояние Фока) на входе

Вопрос: Можно ли выразить влияние простого светоделителя на 50% на количество фотонов с помощью матриц, чтобы результат можно было вычислить с помощью матричных вычислений, а не путем ручной подстановки уравнений?

Чтобы объяснить проблему, рассмотрим светоделитель 50% и определим: а 1 , 2 + = операторы создания для ввода и б 1 , 2 + = операторы создания для вывода. Для сбалансированного 50%-ного светоделителя мы можем написать, как светоделитель сопоставляет ввод с операторами создания вывода:

а 1 + 1 2 ( б 1 + + б 2 + )
а 2 + 1 2 ( б 1 + б 2 + )

Затем мы можем вычислить распределение количества фотонов на выходе для заданных состояний числа фотонов на входе, записав состояния Фока в терминах операторов рождения с последующей подстановкой в ​​приведенные выше уравнения.

Например, рассмотрим 2 фотона, падающих на светоделитель в режиме MODE1, и 1 фотон, падающий на вход MODE2. Сначала мы запишем их в виде числовых состояний фотонов, а затем выразим в терминах операторов рождения, действующих на вакуум.

| ψ я н "=" | 2 1 | 1 2 "=" 1 2 а 1 + 2 а 2 + | 0 1 | 0 2

Мы можем найти выходное состояние, переписав а 1 + & а 2 + с точки зрения б 1 + & б 2 + (используя приведенные выше уравнения):

| ψ о ты т "=" 1 2 2 ( б 1 + + б 2 + ) 2 ( б 1 + б 2 + ) | 0 1 | 0 2

который затем можно оценить, раскрыв скобки и оценив влияние операторов рождения на вакуум (т.е. б 1 + | 0 1 "=" | 1 1 и б 1 + 2 | 0 1 "=" 2 | 2 1 и т. д.). Результат этого примера таков:

| ψ о ты т "=" 1 4 ( 6 | 3 1 | 0 2 + 2 | 2 1 | 1 2 2 | 1 1 | 2 2 6 | 0 1 | 3 2 )

Описанный выше метод работает, но требует подстановки алгебраического уравнения, что не является математической операцией. Мой вопрос заключается в том, как мы можем сформулировать проблему так, чтобы входные данные можно было вводить в виде вектора/матрицы, а выходные данные (например, в виде вектора, показывающего коэффициенты каждого возможного выходного состояния числа фотонов) вычислялись с помощью матричных операций, которые хорошо подходит для исполнения на компьютере. Я обычно видел эффект светоделителя, выраженный в виде матрицы: например

Б "=" 1 2 ( 1 1 1 1 )
[ссылка: https://physics.stackexchange.com/questions/330888/how-to-write-the-output-state-of-a-beam-splitter]

так что потенциально можно было бы написать:

( б 1 + б 2 + ) "=" 1 2 ( 1 1 1 1 ) ( а 1 + а 2 + )

но это выражает отношения только в терминах одиночных операторов создания, и неясно, как сформулировать матричную задачу для поддержки произвольных входных данных (т. е. когда каждый вход может содержать несколько фотонов).

Возможные решения : Из моего исследования я отмечаю:

Все предложения или идеи очень ценятся, спасибо.

Конечно, это возможно, просто проработайте это, заменив действие операторов создания/уничтожения, как вы обсуждали выше, состояниями Фока. Отвечает ли это на ваш вопрос? Или вам нужно простое закрытое выражение для такого преобразования? (Обоснованное предположение: может быть не так уж сложно угадать и доказать это с помощью индукции или чего-то подобного.)
Спасибо за предложение. Рад слышать, что это возможно, но не могли бы вы показать пример, потому что я не совсем понимаю, как бы я это сделал. Я ищу матричное выражение/подход, чтобы я мог закодировать его на общем языке программирования, а затем получить распределение числа выходных фотонов для любого распределения числа входных фотонов. например, если я ввожу входные данные как 2 фотона в РЕЖИМ1 и 1 фотон в РЕЖИМ2, как в моем примере, я бы хотел, чтобы выходные данные позволяли мне видеть коэффициенты каждого возможного состояния номера выходного фотона (чтобы я мог возвести их в квадрат, чтобы найти вероятность возникновения ). Спасибо
Что ж, для случая не более одного фотона в каждом падающем луче вы сами разобрались выше! (Извините, проверил еще раз: не совсем то, что вы сделали, но вы можете сразу изменить то, что вы сделали для этого случая.)
Я пытался решить это, но не могу получить ни одного выражения для вывода. В частности, светоделитель представляет собой матрицу 2 x 2 и действует на вектор (2 x 1), в результате чего получается вектор (2 x 1). Это означает, что операторы b1 и b2 разделены. Однако известно, что некоторые члены выходных данных требуют ОБА b1 и b2 вместе, например, |1>|2>, так как же выходной вектор (2 x 1) может содержать эту информацию? Поскольку существует 4 возможных состояния вывода (|3>|0>, |2>|1>, |1>|2> и |0>|3>, мне нужно как-то расширить размер моей проблемы? Извините, если Я упускаю что-то очевидное? Спасибо.
Действительно, вы хотите расширить, по крайней мере, ваше выходное пространство до полного пространства, которое вы можете охватить 3 фотонами, то есть 4 состояния, которые вы упомянули. (Вы также можете расширить входное пространство до этого измерения, но это не обязательно.) Каждый фотон может двигаться в любом направлении, поэтому априори вы можете получить любое количество фотонов на любом выходе (если нет деструктивных помех, таких как в ХОМ.)
Спасибо за помощь. Исходя из этого и первого комментария, могу ли я сформировать матрицу M так, чтобы: MX = Y, где Y — выход, а X — входное состояние. В моем примере я мог бы позволить X = (0, 1, 0, 0), чтобы указать состояние входа |2>|1>, а выход Y должен быть: 1/4(√6, √2, -√2, -√6). M должна быть матрицей 4 x 4, которая включает в себя действие преобразования состояний Фока в оператор создания, а также операцию светоделителя (например, матрицу B 2 x 2 в моем исходном посте, но, возможно, расположенную по диагонали матрицы 4 x 4). ). Это правильно? Или же (лучшее) решение включает в себя внешние продукты для расширения размерности?
Я думаю, это имеет смысл. Я не думаю, что вы сможете написать его с помощью простого внешнего произведения из-за коммутационных отношений (в основном комбинаторного фактора sqrt (n)). Но я уверен, что есть хорошая закрытая форма (с каким-то комбинаторным объектом, биномиальным коэффициентом и т.п.).

Ответы (1)

Я не совсем уверен, что вы имеете в виду, но я думаю, что это легко сделать.

Что вам нужно помнить, так это то, что ваш светоделитель действует на продукты операторов создания как преобразование продукта. Таким образом, если ваш светоделитель представлен матрицей

(1) U "=" ( U 11 U 12 U 21 U 22 )
затем:
а я Дж а Дж U Дж я , (2) а 1 а 1 U 11 + а 2 U 21 , а 2 а 2 U 12 + а 2 U 22
и, таким образом, на состоянии продукта
| в "=" 1 2 а 1 а 1 а 2 | 0
вы получаете вывод
(3) | вне "=" 1 2 я Дж к а я а Дж а к U я 1 U Дж 1 U к 2 .
Обратите внимание, что ничто не ограничивает (1) до 2 × 2 матрицы. В целом индексы я Дж к в (3) выйдет за пределы режимов, поэтому в 3 -канальный интерферометр, от которого они будут работать 1 к 3 например.

В вашем конкретном примере коэффициент правильно нормализованного состояния вывода состояния 1 3 ! | 3 1 | 0 находится путем ограничения суммы на все комбинации ( я Дж к ) которые умножают коэффициент ( а 1 ) 3 . вот это явно ( я Дж к ) "=" ( 1 , 1 , 1 ) поэтому коэффициент 1 2 U 11 U 11 U 12 "=" 1 4 . Таким образом, коэффициент | 3 1 | 0 2 "=" 6 / 4 .

Коэффициент правильно нормированного состояния 1 2 | 2 1 | 1 2 будет суммой слагаемых с ( я Дж к ) "=" ( 1 , 1 , 2 ) , ( 1 , 2 , 1 ) или ( 2 , 1 , 1 ) :

1 2 ( U 11 U 11 U 22 + U 11 U 21 U 12 + U 21 U 11 U 12 ) "=" 1 2 1 2 2 ( 1 + ( 1 ) + ( 1 ) ) "=" 1 4
поэтому коэффициент | 2 1 | 1 2 является 2 4 . У меня есть разница в знаках с вами, вероятно, из-за того, как я размещаю свои индексы в (2): ваше определение преобразования может иметь индексы ( я Дж ) перевернулся с моим способом написания вещей. Если вы перевернете индексы, вы получите 2 4 .

Обратите внимание, что в моих обозначениях операторы создания отображаются в векторы

а 1 | 0 | 1 "=" ( 1 0 ) , а 2 | 0 | 2 "=" ( 0 1 ) ,
так что
U | 1 "=" | 1 1 | U | 1 + | 2 2 | U | 1 , "=" | 1 U 11 + | 2 U 21 , "=" а 1 | 0 U 11 + а 2 | 0 U 21
так что действительно (2) выглядит правильно.

Есть еще...

Вы можете еще больше упростить задачу, используя представление Швингера (или карту Жордана )

так что состояние углового момента

| Дж м ( а 1 ) Дж + м ( а 2 ) Дж м ( Дж + м ) ! ( Дж м ) ! | 0 .
Таким образом, ваше состояние ввода отображается в состояние углового момента
| ψ я н | 3 2 , 1 2 .
Тогда ваше светоделительное преобразование соответствует (см. Yurke, B., McCall, SL и Klauder, JR, 1986. SU (2) и SU (1, 1) интерферометры. Physical Review A, 33(6), p.4033.) ротация U "=" р г ( 0 ) р у ( π / 2 ) р г ( 0 ) чтобы потом можно было одним махом написать
U | ψ я н р у ( π / 2 ) | 3 2 , 1 2 , "=" м | 3 2 , м Д м , 1 / 2 3 / 2 ( 0 , π / 2 , 0 )
где Д м м Дж ( α , β , γ ) является D-матрицей Вигнера .

Тогда легко получить результат в терминах состояний углового момента:

| 3 2 , 3 2 ( 3 2 2 ) + | 3 2 , 3 2 ( 1 2 2 ) | 3 2 , 1 2 + 1 2 2 | 3 2 , 1 2 + | 3 2 , 3 2 ( 3 2 2 )
а затем преобразовать обратно в состояния числа занятий.

Спасибо за очень полезный ответ. Первый подход, кажется, работает хорошо (второй, использующий Dmatrix, тоже хорош, но мне нужно больше читать, чтобы понять отображение состояний углового момента). Мой единственный оставшийся вопрос: как мы можем использовать этот подход с входными данными в состояниях суперпозиции. Давайте рассмотрим пример, когда входной фотон находится в состоянии суперпозиции для обоих режимов ввода: |INPUT> = 1/sqrt(2) * (|0>|1> +|1>|0>). В конечном итоге я хочу масштабировать этот анализ, чтобы учитывать больше режимов (2 временных интервала + 2 входных пути), поэтому я думаю, что могу просто использовать больше индексов для большего количества режимов, как вы предлагаете. Спасибо
@Qconfused9102 Лучшее место для чтения о представителе Schwinger. Лекции Гордона Бейма по квантовой механике. Да, раздача линейных комбинаций произвольных состояний становится бременем, если вы используете метод умножения, и у вас есть действие на каждую часть комбо отдельно, после чего следите за коэффициентом комбо. С этой точки зрения Д -function трюк помогает, так как пример в вашем комментарии будет похож на вращение 1 2 ( | 1 / 2 , 1 / 2 > + | 1 / 2 , 1 / 2 > ) .
еще раз спасибо за это. Ради интереса, как вы думаете, может ли этот подход работать и с когерентными состояниями? то есть замена операторов создания операторами смещения... либо в исходном подходе, либо в подходе Джордана-Швингера?
Этот трюк Швингера в основном для С U ( н ) поэтому, если вы генерируете когерентное состояние углового момента, тогда да. В случае стандартного когерентного состояния можно расширить 2-модовое когерентное состояние до С U ( 2 ) подпространство в основном с использованием карты Джордана, т.е. взять состояние | н 1 н 2 в вашем 2-режимном когерентном состоянии и думайте об этом как о состоянии углового момента, но в остальном я не на 100% то, что вы делаете, «заменяя операторы создания операторами смещения». Предположительно, вы могли бы заменить вращения смещениями...
... и тогда вы получите Д -функции для группы трансляций (думаю, это функции Бесселя).