Линдбладиан

редактировать
Марковское квантовое главное уравнение для матриц плотности (смешанные состояния)

В квантовой механике Уравнение Горини – Косаковского – Сударшана – Линдблада (уравнение ГКСЛ, названное в честь, Анджей Косаковски, Джордж Сударшан и Горан Линдблад ), основное уравнение в форме Линдблада, квантовое лиувиллианское или линдбладианское является наиболее общим типом марковского и основное уравнение, описывающее (в общем, неунитарное) эволюцию матрицы плотности ρ, которая сохраняет законы квантовой механики (т. е. сохраняет след и полностью положительна для любое начальное условие).

Уравнение Шредингера является частным случаем более общего уравнения Линдблада, которое привело к некоторым предположениям, что квантовая механика может быть продуктивно расширена и расширена за счет дальнейшего применения и анализ уравнения Линдблада. Уравнение Шредингера имеет дело с векторами состояний, которые могут описывать только чистые квантовые состояния и, таким образом, являются менее общими, чем матрицы плотности, которые могут описывать смешанные состояния. тоже.

Содержание
  • 1 Мотивация
  • 2 Определение
    • 2.1 Квантовая динамическая полугруппа
    • 2.2 Свойства инвариантности
    • 2.3 Картинка Гейзенберга
  • 3 Физическое происхождение
  • 4 Примеры
  • 5 См. Также
  • 6 Ссылки
  • 7 Внешние ссылки
Мотивация

В канонической формулировке квантовой механики эволюция системы во времени определяется унитарной динамикой. Это означает, что затухания нет и фазовая когерентность сохраняется на протяжении всего процесса, и является следствием того, что учитываются все участвующие степени свободы. Однако любая реальная физическая система не является абсолютно изолированной и будет взаимодействовать со своим окружением. Это взаимодействие со степенями свободы, внешними по отношению к системе, приводит к рассеиванию энергии в окружающую среду, вызывая распад и рандомизацию фазы. Эти эффекты являются причиной того, что квантовую механику трудно наблюдать в макроскопическом масштабе. Более того, понимание взаимодействия квантовой системы с окружающей средой необходимо для понимания многих обычно наблюдаемых явлений, таких как спонтанное излучение света возбужденными атомами или работа многих квантовых технологических устройств, таких как лазер.

Некоторые математические методы были введены для изучения взаимодействия квантовой системы с окружающей средой. Одним из них является использование матрицы плотности и связанного с ней основного уравнения. Хотя в принципе этот подход к решению квантовой динамики эквивалентен картине Шредингера или картине Гейзенберга, он позволяет более легко включить несвязные процессы, которые представляют собой взаимодействия с окружающей средой. Оператор плотности обладает тем свойством, что он может представлять классическую смесь квантовых состояний, и поэтому жизненно важен для точного описания динамики так называемых открытых квантовых систем.

Определение

В более общем плане основное уравнение Линдблада для матрицы плотности ρ N-мерной системы может быть записано как (для педагогического введения вы можете обратиться к)

ρ ˙ = - я ℏ [H, ρ] + ∑ N, m = 1 N 2 - 1 hnm (A n ρ A m † - 1 2 {A m † A n, ρ}) {\ displaystyle {\ dot {\ rho}} = - {i \ over \ hbar} [H, \ rho] + \ sum _ {n, m = 1} ^ {N ^ {2} -1} h_ {nm} \ left (A_ {n} \ rho A_ {m} ^ {\ dagger} - {\ frac {1} {2}} \ left \ {A_ {m} ^ {\ dagger} A_ {n}, \ rho \ right \} \ right)}{\ displaystyle {\ dot {\ rho}} = - {i \ over \ hbar} [H, \ rho] + \ sum _ {n, m = 1} ^ {N ^ {2} -1} h_ {nm} \ left (A_ {n} \ rho A_ {m} ^ {\ dagger} - {\ frac {1} {2}} \ left \ {A_ {m} ^ {\ dagger} A_ {n}, \ rho \ right \} \ right)}

где H - (эрмитова ) гамильтонова часть, а {A m} {\ displaystyle \ {A_ {m} \}}{\ displaystyle \ {A_ {m} \}} - произвольный ортонормированный базис операторов Гильберта-Шмидта на гильбертовом пространстве системы с ограничением, что A N пропорционально тождественному оператору. Наше соглашение подразумевает, что другие A m не имеют следов, и обратите внимание, что суммирование выполняется только до N - 1, таким образом исключая единственную базисную матрицу с ненулевым следом. Матрица коэффициентов h вместе с гамильтонианом определяет динамику системы. Матрица h должна быть положительно полуопределенной, чтобы гарантировать, что уравнение сохраняет след и полностью положительно. антикоммутатор определяется как {a, b} = a b + b a. {\ displaystyle \ {a, b \} = ab + ba.}\ { a, b \} = ab + ba.

Если все h mn равны нулю, то это сводится к квантовому уравнению Лиувилля для замкнутого система, ρ ˙ знак равно - (я / ℏ) [H, ρ] {\ displaystyle {\ dot {\ rho}} = - (i / \ HBAR) [H, \ rho]}{\ displaystyle {\ dot {\ rho}} = - (i / \ hbar) [H, \ rho]} . Это также известно как уравнение фон Неймана и является квантовым аналогом классического уравнения Лиувилля.

. Поскольку матрица h положительно полуопределенная, ее можно диагонализовать с помощью унитарного преобразование u:

u † hu = [γ 1 0 ⋯ 0 0 γ 2 ⋯ 0 ⋮ ⋮ ⋱ ⋮ 0 0 ⋯ γ N 2 - 1] {\ displaystyle u ^ {\ dagger} hu = {\ begin {bmatrix} \ gamma _ {1} 0 \ cdots 0 \\ 0 \ gamma _ {2} \ cdots 0 \\\ vdots \ vdots \ ddots \ vdots \\ 0 0 \ cdots \ gamma _ { N ^ {2} -1} \ end {bmatrix}}}u ^ {\ dagger} hu = {\ begin {bmatrix} \ gamma _ {1} 0 \ cdots 0 \\ 0 \ gamma _ {2} \ cdots 0 \\\ vdots \ vdots \ ddots \ vdots \\ 0 0 \ cdots \ gamma _ {{N ^ {2} -1}} \ end {bmatrix}}

где собственные значения γ i неотрицательны. Если мы определим базис другого ортонормированного оператора

L i = ∑ j = 1 N 2 - 1 uji A j {\ displaystyle L_ {i} = \ sum _ {j = 1} ^ {N ^ {2} -1} u_ {ji} A_ {j}}{\ displaystyle L_ {i} = \ sum _ {j = 1} ^ {N ^ {2} -1} u_ {ji} A_ {j} }

мы можем переписать уравнение Линдблада в диагональной форме

ρ ˙ = - i ℏ [H, ρ] + ∑ i = 1 N 2 - 1 γ i (L i ρ L i † - 1 2 {L i † L i, ρ}). {\ Displaystyle {\ точка {\ rho}} = - {я \ над \ hbar} [H, \ rho] + \ сумма _ {я = 1} ^ {N ^ {2} -1} \ gamma _ {я } \ left (L_ {i} \ rho L_ {i} ^ {\ dagger} - {\ frac {1} {2}} \ left \ {L_ {i} ^ {\ dagger} L_ {i}, \ rho \ right \} \ right).}{\ displaystyle {\ dot {\ rho}} = - {i \ над \ hbar} [H, \ rho] + \ sum _ {i = 1} ^ {N ^ {2} -1} \ gamma _ {i} \ left (L_ {i} \ rho L_ {i} ^ { \ dagger} - {\ frac {1} {2}} \ left \ {L_ {i} ^ {\ dagger} L_ {i}, \ rho \ right \} \ right).}

Новые операторы L i обычно называют операторами Линдблада или операторами перехода системы.

Квантовая динамическая полугруппа

Карты, сгенерированные линдбладианом за разное время, вместе называются квантовой динамической полугруппой - семейством квантовых динамических отображений ϕ t {\ displaystyle \ phi _ {t}}\ phi _ {t} на пространстве матриц плотности, индексированных одним временным параметром t ≥ 0 {\ displaystyle t \ geq 0}t \ geq 0 , которые подчиняются полугруппе свойству

ϕ s (ϕ t (ρ)) = ϕ t + s (ρ), t, s ≥ 0. { \ Displaystyle \ phi _ {s} (\ phi _ {t} (\ rho)) = \ phi _ {t + s} (\ rho), \ qquad t, s \ geq 0.}\ phi_s (\ phi_t (\ rho)) = \ phi_ {t + s} ( \ rho), \ qquad t, s \ ge 0.

Уравнение Линдблада можно получить следующим образом:

L (ρ) = lim Δ t → 0 ϕ Δ t (ρ) - ϕ 0 (ρ) Δ t {\ displaystyle {\ mathcal {L}} (\ rho) = \ mathrm {lim } _ {\ Delta t \ to 0} {\ frac {\ phi _ {\ Delta t} (\ rho) - \ phi _ {0} (\ rho)} {\ Delta t}}}\ mathcal {L} (\ rho) = \ mathrm {lim} _ {\ Delta t \ to 0} \ frac {\ phi _ {\ Delta t} (\ rho) - \ phi_0 (\ rho)} {\ Delta t}

который, по линейности ϕ t {\ displaystyle \ phi _ {t}}\ phi _ {t} , является линейным супероператором. Полугруппа восстанавливается как

ϕ t + s (ρ) = e L s ϕ t (ρ). {\ displaystyle \ phi _ {t + s} (\ rho) = e ^ {{\ mathcal {L}} s} \ phi _ {t} (\ rho).}\ phi_ {t + s} (\ rho) = e ^ {\ mathcal {L} s} \ phi_t ( \ rho).

Свойства инвариантности

Уравнение Линдблада инвариантно относительно любого унитарного преобразования v операторов и констант Линдблада,

γ i L i → γ i ′ L i ′ = ∑ j = 1 N 2 - 1 vij γ j L j, {\ displaystyle { \ sqrt {\ gamma _ {i}}} L_ {i} \ to {\ sqrt {\ gamma _ {i} '}} L_ {i}' = \ sum _ {j = 1} ^ {N ^ {2 } -1} v_ {ij} {\ sqrt {\ gamma _ {j}}} L_ {j},}{\displaystyle {\sqrt {\gamma _{i}}}L_{i}\to {\sqrt {\gamma _{i}'}}L_{i}'=\sum _{j=1}^{N^{2}-1}v_{ij}{\sqrt {\gamma _{j}}}L_{j},}

, а также при неоднородном преобразовании

L i → L i ′ = L i + ai I, {\ Displaystyle L_ {i} \ к L_ {i} '= L_ {i} + a_ {i} I,}{\displaystyle L_{i}\to L_{i}'=L_{i}+a_{i}I,}
H → H' = H + 1 2 я ∑ j = 1 N 2 - 1 γ j (aj ∗ L j - aj L j †) + b I, {\ displaystyle H \ to H '= H + {\ frac {1} {2i}} \ sum _ {j = 1} ^ {N ^ {2 } -1} \ gamma _ {j} \ left (a_ {j} ^ {*} L_ {j} -a_ {j} L_ {j} ^ {\ dagger} \ right) + bI,}{\displaystyle H\to H'=H+{\frac {1}{2i}}\sum _{j=1}^{N^{2}-1}\gamma _{j}\left(a_{j}^{*}L_{j}-a_{j}L_{j}^{\dagger }\right)+bI,}

где a i - комплексные числа, а b - действительное число. Однако первое преобразование нарушает ортонормированность операторов L i (если все γ i не равны), а второе преобразование уничтожает бесследовательность. Следовательно, с точностью до вырождения среди γ i, L i диагональной формы уравнения Линдблада однозначно определяются динамикой, если мы требуем, чтобы они были ортонормированными и бесследными..

Картинка Гейзенберга

Эволюция типа Линдблада матрицы плотности на картине Шредингера может быть эквивалентно описана на картинке Гейзенберга, используя следующее (диагонализованное) уравнение движения для каждой квантовой наблюдаемой X:

X ˙ = i ℏ [H, X] + ∑ i = 1 N 2 - 1 γ i (L i † XL i - 1 2 {L i † L i, X}). {\ displaystyle {\ dot {X}} = {\ frac {i} {\ hbar}} [H, X] + \ sum _ {i = 1} ^ {N ^ {2} -1} \ gamma _ { i} \ left (L_ {i} ^ {\ dagger} XL_ {i} - {\ frac {1} {2}} \ left \ {L_ {i} ^ {\ dagger} L_ {i}, X \ right \} \ right).}{\ displaystyle {\ dot {X}} = {\ frac {i} {\ hbar} } [H, X] + \ sum _ {i = 1} ^ {N ^ {2} -1} \ gamma _ {i} \ left (L_ {i} ^ {\ dagger} XL_ {i} - {\ гидроразрыв {1} {2}} \ left \ {L_ {i} ^ {\ dagger} L_ {i}, X \ right \} \ right).}

Подобное уравнение описывает временную эволюцию ожидаемых значений наблюдаемых, задаваемых теоремой Эренфеста. В соответствии со свойством сохранения следа уравнения Линдблада изображения Шредингера уравнение изображения Гейзенберга является унитальным, т.е. оно сохраняет тождественный оператор.

Физическое происхождение

Основное уравнение Линдблада описывает эволюцию различных типов открытых квантовых систем, например система, слабо связанная с марковским резервуаром. Обратите внимание, что H, фигурирующая в уравнении, не обязательно совпадает с гамильтонианом затравочной системы, но также может включать эффективную унитарную динамику, возникающую в результате взаимодействия системы и окружающей среды.

Эвристический вывод, например, в примечаниях Прескилла, начинается с более общей формы открытой квантовой системы и преобразует ее в форму Линдблада, делая марковское предположение и расширяясь за короткое время. Более физически мотивированная стандартная трактовка охватывает три общих типа выводов линдбладиана, начиная с гамильтониана, действующего как на систему, так и на окружающую среду: предел слабой связи (подробно описанный ниже), приближение низкой плотности и предел сингулярной связи. Каждый из них основан на определенных физических предположениях, касающихся, например, корреляционных функций окружающей среды. Например, при выводе предела слабой связи обычно предполагается, что (а) корреляции системы с окружающей средой развиваются медленно, (б) возбуждения среды, вызванные распадом системы, и (в) члены, которые быстро осциллируют. по сравнению с интересующей системой шкалой времени можно пренебречь. Эти три приближения называются борновской, марковской и вращающейся волной соответственно.

Вывод предела слабой связи предполагает квантовую систему с конечным числом степеней свободы, соединенную с термостатом, содержащим бесконечное количество степеней свободы. свобода. Система и ванна обладают гамильтонианом, записанным в терминах операторов, действующих только на соответствующем подпространстве полного гильбертова пространства. Эти гамильтонианы управляют внутренней динамикой несвязанной системы и ванны. Существует третий гамильтониан, который содержит произведения операторов системы и ванны, таким образом связывая систему и ванну. Самая общая форма этого гамильтониана:

H = HS + HB + HBS {\ displaystyle H = H_ {S} + H_ {B} + H_ {BS} \,}H = H_S + H_B + H_ {BS} \,

Динамика всей системы может описывается уравнением движения Лиувилля, χ ˙ = - i [H, χ] {\ displaystyle {\ dot {\ chi}} = - i [H, \ chi]}\ dot {\ chi} = - i [H, \ chi] . Это уравнение, содержащее бесконечное количество степеней свободы, невозможно решить аналитически, за исключением очень частных случаев. Более того, при определенных приближениях нет необходимости рассматривать степени свободы ванны, и эффективное главное уравнение может быть получено в терминах матрицы плотности системы, ρ = tr B ⁡ χ {\ displaystyle \ rho = \ operatorname {tr} _ {B} \ chi}\ rho = \ operatorname {tr} _B \ chi . Проблему легче проанализировать, перейдя к картине взаимодействия, определяемой унитарным преобразованием M ~ = U 0 MU 0 † {\ displaystyle {\ tilde {M}} = U_ {0} MU_ {0} ^ {\ dagger}}{\ displaystyle {\ tilde {M}} = U_ {0} MU_ {0} ^ {\ dagger}} , где M {\ displaystyle M}M - произвольный оператор, а U 0 = ei (HS + HB) t {\ displaystyle U_ {0} = e ^ {i (H_ {S} + H_ {B}) t}}{\ displaystyle U_ {0} = e ^ {i (H_ {S} + H_ {B}) t}} . Также обратите внимание, что U (t, t 0) {\ displaystyle U (t, t_ {0})}{\ displaystyle U (t, t_ {0})} является полным унитарным оператором всей системы. Несложно подтвердить, что уравнение Лиувилля принимает вид

χ ~ ˙ = - i [H ~ BS, χ ~] {\ displaystyle {\ dot {\ tilde {\ chi}}} = - i [{\ tilde { H}} _ {BS}, {\ tilde {\ chi}}] \,}\ точка {\ тильда {\ чи}} = - я [\ тильда {H} _ {BS}, \ тильда {\ хи}] \,

где гамильтониан H ~ BS = ei (HS + HB) t HBS e - i (HS + HB) t {\ displaystyle {\ tilde {H}} _ {BS} = e ^ {i (H_ {S} + H_ {B}) t} H_ {BS} e ^ {- i (H_ {S} + H_ {B }) t}}\ tilde {H} _ {BS} = e ^ {я (H_S + H_B) t} H_ {BS} e ^ {- i (H_S + H_B) t} явно зависит от времени. Кроме того, согласно изображению взаимодействия, χ ~ = UBS (t, t 0) χ UBS † (t, t 0) {\ displaystyle {\ tilde {\ chi}} = U_ {BS} (t, t_ {0}) \ чи U_ {BS} ^ {\ dagger} (t, t_ {0})}{\ displaystyle {\ тильда {\ chi}} = U_ {BS} (t, t_ {0}) \ chi U_ {BS} ^ {\ dagger} (t, t_ {0})} , где UBS = U 0 † U (t, t 0) {\ displaystyle U_ {BS} = U_ {0} ^ {\ dagger} U (t, t_ {0})}{\ displaystyle U_ {BS} = U_ {0} ^ {\ dagger} U (t, t_ {0})} . Это уравнение можно интегрировать напрямую, чтобы получить

χ ~ (t) = χ ~ (0) - i ∫ 0 tdt ′ [H ~ BS (t ′), χ ~ (t ′)] {\ displaystyle {\ tilde {\ chi}} (t) = {\ tilde {\ chi}} (0) -i \ int _ {0} ^ {t} dt '[{\ tilde {H}} _ {BS} (t'), {\ tilde {\ chi}} (t ')]} \tilde{\chi}(t)=\tilde{\chi}(0) -i\int^t_0 dt' [\tilde{H}_{BS}(t'),\tilde{\chi}(t')]

Это неявное уравнение для χ ~ {\ displaystyle {\ tilde {\ chi}}}\ tilde {\ chi} можно снова подставить в уравнение Лиувилля для получения точного дифференциально-интегрального уравнения

χ ~ ˙ = - i [H ~ BS, χ ~ (0)] - ∫ 0 tdt ′ [H ~ BS (t), [H ~ BS (t ′), Χ ~ (t ′)]] {\ Displaystyle {\ dot {\ tilde {\ chi}}} = - я [{\ tilde {H}} _ {BS}, {\ tilde {\ chi}} (0)] - \ int _ {0} ^ {t} dt '[{\ tilde {H}} _ {BS} (t), [{\ tilde {H}} _ {BS} (t'), {\ tilde {\ chi}} (t ')]]} \dot{\tilde{\chi}}=-i[\tilde{H}_{BS},\tilde{\chi}(0)] - \int^t_0 dt' [\tilde{H}_{BS}(t),[\tilde{H}_{BS}(t'),\tilde{\chi}(t')]]

Мы продолжаем вывод, предполагая, что взаимодействие инициируется в t = 0 {\ displaystyle t = 0}t = 0 , и в то время нет корреляции между системой и ванной. Это означает, что начальное условие факторизуемо как χ (0) = ρ (0) R 0 {\ displaystyle \ chi (0) = \ rho (0) R_ {0}}\ chi (0) = \ rho (0) R_0 , где R 0 {\ displaystyle R_ {0}}R_ {0} изначально является оператором плотности ванны.

Трассировка по степеням свободы ванны, tr R ⁡ χ ~ = ρ ~ {\ displaystyle \ operatorname {tr} _ {R} {\ tilde {\ chi}} = {\ tilde { \ rho}}}\ operatorname {tr} _R \ тильда {\ chi} = \ тильда {\ rho} вышеупомянутого дифференциально-интегрального уравнения дает

ρ ~ ˙ = - ∫ 0 tdt ′ tr R ⁡ {[H ~ BS (t), [H ~ BS (t ′), Χ ~ (t ′)]]} {\ displaystyle {\ dot {\ tilde {\ rho}}} = - \ int _ {0} ^ {t} dt '\ operatorname {tr} _ {R} \ {[{\ tilde {H}} _ {BS} (t), [{\ tilde {H}} _ {BS} (t '), {\ tilde {\ chi}} (t')]] \ }} \dot{\tilde{\rho}}= - \int^t_0 dt' \operatorname{tr}_R\{[\tilde{H}_{BS}(t),[\tilde{H}_{BS}(t'),\tilde{\chi}(t')]]\}

Это уравнение является точным для временной динамики матрицы плотности системы, но требует полного знания динамики степеней свободы ванны. Упрощающее допущение, называемое приближением Борна, основывается на большом размере ванны и относительной слабости связи, то есть связь системы с ванной не должна существенно изменять собственные состояния ванны. В этом случае полная матрица плотности может быть разложена на все времена как χ ~ (t) = ρ ~ (t) R 0 {\ displaystyle {\ tilde {\ chi}} (t) = {\ tilde {\ rho }} (t) R_ {0}}\ тильда {\ чи} (t) = \ тильда {\ rho} (t) R_0 . Основное уравнение принимает вид

ρ ~ ˙ = - ∫ 0 tdt ′ tr R ⁡ {[H ~ BS (t), [H ~ BS (t ′), ρ ~ (t ′) R 0]]} {\ displaystyle {\ dot {\ tilde {\ rho}}} = - \ int _ {0} ^ {t} dt '\ operatorname {tr} _ {R} \ {[{\ tilde {H}} _ {BS} (t), [{\ tilde {H}} _ {BS} (t '), {\ tilde {\ rho}} (t') R_ {0}]] \}} \dot{\tilde{\rho}}= - \int^t_0 dt' \operatorname{tr}_R\{[\tilde{H}_{BS}(t),[\tilde{H}_{BS}(t'),\tilde{\rho}(t')R_0]]\}

Теперь уравнение явное в системе степеней свободы, но это очень сложно решить. Последним предположением является приближение Борна-Маркова, согласно которому производная по времени матрицы плотности зависит только от ее текущего состояния, а не от ее прошлого. Это предположение справедливо при быстрой динамике ванны, когда корреляции внутри ванны теряются очень быстро, и сводится к замене ρ (t ') → ρ (t) {\ displaystyle \ rho (t') \ rightarrow \ rho ( t)} \rho(t')\rightarrow \rho(t)в правой части уравнения.

ρ ~ ˙ = - ∫ 0 tdt ′ tr R ⁡ {[H ~ BS (t), [H ~ BS (t ′), ρ ~ (t) R 0]]} {\ displaystyle {\ dot { \ tilde {\ rho}}} = - \ int _ {0} ^ {t} dt '\ operatorname {tr} _ {R} \ {[{\ tilde {H}} _ {BS} (t), [ {\ tilde {H}} _ {BS} (t '), {\ tilde {\ rho}} (t) R_ {0}]] \}} \dot{\tilde{\rho}}= - \int^t_0 dt' \operatorname{tr}_R\{[\tilde{H}_{BS}(t),[\tilde{H}_{BS}(t'),\tilde{\rho}(t)R_0]]\}

Если предположить, что гамильтониан взаимодействия имеет вид

HBS = ∑ α я Γ i {\ displaystyle H_ {BS} = \ sum \ alpha _ {i} \ Gamma _ {i}}H_ {BS} = \ sum \ alpha_i \ Gamma_i

для системных операторов α i {\ displaystyle \ alpha _ { i}}\ alpha _ {i} и операторы ванны Γ i {\ displaystyle \ Gamma _ {i}}\ Gamma_i , основное уравнение принимает вид

ρ ~ ˙ = - ∑ ∫ 0 tdt ′ Тр р ⁡ {[α я (t) Γ я (t), [α j (t ′) Γ j (t ′), ρ ~ (t) R 0]]} {\ displaystyle {\ dot {\ tilde {\ rho}}} = - \ sum \ int _ {0} ^ {t} dt '\ operatorname {tr} _ {R} \ {[\ alpha _ {i} (t) \ Gamma _ {i} ( t), [\ alpha _ {j} (t ') \ Gamma _ {j} (t'), {\ tilde {\ rho}} (t) R_ {0}]] \}} \dot{\tilde{\rho}}= - \sum\int^t_0 dt' \operatorname{tr}_R\{[\alpha_i(t) \Gamma_i(t),[\alpha_j(t') \Gamma_j(t'),\tilde{\rho}(t)R_0]]\}

который может раскладывается как

ρ ~ ˙ = - ∑ ∫ 0 tdt ′ (α i (t) α j (t ′) ρ - α j (t ′) ρ (t) α i (t)) ⟨Γ i ( t) Γ j (t ′)⟩ + (ρ (t) α J (T ′) α я (T) - α I (T) ρ (t) α J (t ′)) ⟨Γ j (t ′) Γ i (t)⟩ {\ Displaystyle {\ dot {\ tilde) {\ rho}}} = - \ sum \ int _ {0} ^ {t} dt '\ left (\ alpha _ {i} (t) \ alpha _ {j} (t') \ rho - \ alpha _ {j} (t ') \ rho (t) \ alpha _ {i} (t) \ right) \ langle \ Gamma _ {i} (t) \ Gamma _ {j} (t') \ rangle + \ left (\ rho (t) \ alpha _ {j} (t ') \ alpha _ {i} (t) - \ alpha _ {i} (t) \ rho (t) \ alpha _ {j} (t') \ right) \ langle \ Gamma _ {j} (t ') \ Gamma _ {i} (t) \ rangle}\dot{\tilde{\rho}}=- \sum\int^t_0 dt' \left(\alpha_i(t)\alpha_j(t')\rho-\alpha_j(t')\rho(t)\alpha_i(t)\right)\langle\Gamma_i(t)\Gamma_j(t')\rangle + \left(\rho(t)\alpha_j(t')\alpha_i(t)-\alpha_i(t)\rho(t)\alpha_j(t')\right)\langle\Gamma_j(t')\Gamma_i(t)\rangle

Ожидаемые значения ⟨Γ i Γ j⟩ = tr ⁡ {Γ i Γ j R 0} {\ displaystyle \ langle \ Gamma _ {i} \ Gamma _ {j} \ rangle = \ operatorname {tr} \ {\ Gamma _ {i} \ Gamma _ {j} R_ {0} \}}\ langle \ Gamma_i \ Gamma_j \ rangle = \ oper atorname {tr} \ {\ Gamma_i \ Gamma_jR_0 \} относятся к степеням свободы ванны. Предполагая быстрое затухание этих корреляций (в идеале ⟨Γ i (t) Γ j (t ′)⟩ ∝ δ (t - t ′) {\ displaystyle \ langle \ Gamma _ {i} (t) \ Gamma _ {j} (t ') \ rangle \ propto \ delta (t-t')}{\displaystyle \langle \Gamma _{i}(t)\Gamma _{j}(t')\rangle \propto \delta (t-t')}) достигается указанная выше форма супероператора Lindblad L.

Примеры

Для одного F {\ displaystyle F}Fи без единой эволюции супероператор Lindblad , действующий на матрица плотности ρ {\ displaystyle \ rho}\ rho , равно

D [F] (ρ) = F ρ F † - 1 2 (F † F ρ + ρ F † F) {\ displaystyle {\ mathcal {D}} [F] (\ rho) = F \ rho F ^ {\ dagger} - {\ frac {1} {2}} \ left (F ^ {\ dagger} F \ rho + \ rho F ^ {\ dagger} F \ right)}{\ displaystyle {\ mathcal {D}} [F ] (\ rho) = F \ rho F ^ {\ dagger} - {\ frac {1} {2}} \ left (F ^ {\ dagger} F \ rho + \ rho F ^ {\ dagger} F \ right)}

Такой член регулярно встречается в уравнении Линдблада, используемом в квантовой оптике, где он может выражать поглощение или излучение фотоны из резервуара. Если кто-то хочет иметь и поглощение, и излучение, для каждого потребуется оператор перехода. Это приводит к наиболее распространенному уравнению Линдблада, описывающему затухание квантового гармонического осциллятора (представляющего, например, резонатор Фабри – Перо ), соединенного с термостатом, с операторы скачка

F 1 = a, γ 1 = γ 2 (n ¯ + 1), F 2 = a †, γ 2 = γ 2 n ¯. {\ displaystyle {\ begin {align} F_ {1} = a, \ gamma _ {1} = {\ tfrac {\ gamma} {2}} \ left ({\ overline {n}} + 1 \ справа), \\ F_ {2} = a ^ {\ dagger}, \ gamma _ {2} = {\ tfrac {\ gamma} {2}} {\ overline {n}}. \ end {выровнено }}}{\ displaystyle {\ begin {align} F_ {1} = a, \ gamma _ { 1} = {\ tfrac {\ gamma} {2}} \ left ({\ overline {n}} + 1 \ right), \\ F_ {2} = a ^ {\ dagger}, \ gamma _ {2} = {\ tfrac {\ gamma} {2}} {\ overline {n}}. \ End {align}}}

Здесь n ¯ {\ displaystyle {\ overline {n}}}\ overline {n } - это среднее число возбуждений в резервуаре, подавляющее осциллятор, а γ - скорость затухания. Если мы также добавим дополнительную унитарную эволюцию, генерируемую квантовым гармоническим осциллятором с частотой ω c {\ displaystyle \ omega _ {c}}{\ displaystyle \ omega _ {c}} , мы получим

ρ ˙ = - i [ω ca † a, ρ] + D [F 1] (ρ) + D [F 2] (ρ). {\ displaystyle {\ dot {\ rho}} = - я [\ omega _ {c} a ^ {\ dagger} a, \ rho] + {\ mathcal {D}} [F_ {1}] (\ rho) + {\ mathcal {D}} [F_ {2}] (\ rho).}{\ displaystyle {\ dot {\ rho}} = - я [\ omega _ {c} a ^ {\ dagger} a, \ rho] + {\ mathcal {D}} [F_ {1}] (\ rho) + {\ mathcal {D}} [F_ {2}] (\ rho).}

Дополнительные операторы Линдблада могут быть включены для моделирования различных форм дефазировки и колебательной релаксации. Эти методы были включены в методы распространения матрицы плотности на основе сетки.

См. Также
Литература
  • Хрушинский, Дариуш; Паскацио, Саверио. «Краткая история уравнения ГКЛС». arXiv : 1710.05993.
  • Косаковский, А. (1972). «О квантовой статистической механике негамильтоновых систем». Rep. Math. Phys. 3 (4): 247. Bibcode : 1972RpMP.... 3..247K. doi : 10.1016 / 0034-4877 (72) 90010-9.
  • Белавин, A.A.; Зельдович, Б. Я.; Переломов, А.М.; Попов, В. (1969). «Релаксация квантовых систем с эквидистантными спектрами». ЖЭТФ. 29 : 145. Bibcode : 1969JETP... 29..145B.
  • Линдблад, Г. (1976). «О генераторах квантовых динамических полугрупп». Commun. Математика. Phys. 48 (2): 119. Bibcode : 1976CMaPh..48..119L. doi : 10.1007 / BF01608499.
  • Gorini, V.; Косаковский, А.; Сударшан, E.C.G. (1976). «Вполне положительные полугруппы N-уровневых систем». J. Math. Phys. 17 (5): 821. Bibcode : 1976JMP.... 17..821G. doi : 10.1063 / 1.522979.
  • Banks, T.; Сасскинд, Л.; Пескин, М.Е. (1984). «Трудности эволюции чистых состояний в смешанные». Ядерная физика Б. 244 : 125–134. Bibcode : 1984NuPhB.244..125B. doi : 10.1016 / 0550-3213 (84) 90184-6.
  • Аккарди, Луиджи; Лу, Юнь Ган; Волович, И. (2002). Квантовая теория и ее стохастический предел. Нью-Йорк: Springer Verlag. ISBN 978-3-5404-1928-0.
  • Алики, Роберт. «Приглашение в квантовые динамические полугруппы». arXiv : Quant-ph / 0205188.
  • Алики, Роберт; Ленди, Карл (1987). Квантовые динамические полугруппы и приложения. Берлин: Springer Verlag. ISBN 978-0-3871-8276-6.
  • Атталь, Стефан; Джой, Ален; Пилле, Клод-Ален (2006). Открытые квантовые системы II: марковский подход. Springer. ISBN 978-3-5403-0992-5.
  • Gardiner, C.W.; Золлер, Питер (2010). Квантовый шум. Серия Спрингера в синергетике (3-е изд.). Берлин Гейдельберг: Springer-Verlag. ISBN 978-3-642-06094-6.
  • Ingarden, Roman S.; Косаковский, А.; Охя, М. (1997). Информационная динамика и открытые системы: классический и квантовый подход. Нью-Йорк: Springer Verlag. ISBN 978-0-7923-4473-5.
  • Линдблад, Г. (1983). Неравновесная энтропия и необратимость. Дордрехт: Дельта Рейдел. ISBN 1-4020-0320-X.; Comm. Математика. Phys. 48 (1976), 119-130. онлайн
  • Тарасов, Василий Евгеньевич (2008). Квантовая механика негамильтоновых и диссипативных систем. Амстердам, Бостон, Лондон, Нью-Йорк: Elsevier Science. ISBN 978-0-0805-5971-1.
  • Перл, П. (2012). «Простой вывод уравнения Линдблада». European Journal of Physics, 33 (4), 805.
Внешние ссылки
Последняя правка сделана 2021-05-27 10:16:03
Содержание доступно по лицензии CC BY-SA 3.0 (если не указано иное).
Обратная связь: support@alphapedia.ru
Соглашение
О проекте