Master Equations und Operator Sum Form

12

Ich bin eher ein Typ für Quantenoptik als ein Typ für Quanteninformation und beschäftige mich hauptsächlich mit Master-Gleichungen. Ich interessiere mich für das Operator-Summen-Formular und möchte die Fehler in diesem Formular für ein kleines Quantensystem ableiten, das ich simuliere.

Der Haken: Das Quantensystem wird durch ein externes (klassisches) Feld angetrieben, das mit einer Sinusfunktion modelliert ist, und die Dämpfungsraten sind niedrig, so dass ich keine Näherung für rotierende Wellen vornehmen kann, um diese Zeitabhängigkeit zu beseitigen. Da ich die Hauptgleichung numerisch durch Integration lösen muss und das Ergebnis jeder Integration zum Zeitpunkt nicht ausreicht, um diese Fehler herauszufinden, muss ich einige Arbeiten durchführen, um die Superoperatormatrix wiederherzustellen, die mit einer vektorisierten Dichte gearbeitet hat Matrix. dh ich füttere die Master-Gleichung mit einer vektorisierten Dichtematrix mit einer einzelnen Eingabe von 1 und der restlichen Null und baue die Matrix für eine bestimmte Zeit τ so auf . Bin ich hier auf dem richtigen Weg (Sanity Check)? Genauer gesagt, wenn v e c (tτ ist die vektorisierte (also ein Spaltenvektor) Form einer Dichtematrix mit einem einzelnen Eintrag von 1 in Position i , j bei t = 0 , die zur Zeit τ entwickelt wurde , dann eine Matrix nimm die Vektor Form der Dichtematrix von t = 0 bis t = τ ist gegeben als M = Σ i , j v e c ( ρ i j , t = 0 )vec(ρij,t=τ)i,jt=0τt=0t=τ .M=i,jvec(ρij,t=0)vec(ρij,t=τ)

Die Frage: Angesichts dieses Superoperators , der M tutM , wie kann ich KraussBetreiber für die Bediener-BetragHöhe von M , das in einer nützlichen Form ist? dh es handelt sich um ein Qubit oder Qutrit und ein anderes Qubit oder Qutrit. Ich möchte in der Lage sein, die Operator-Summe, wenn möglich, in Form von Tensor-Produkten von Spin-Matrizen auf jedem Kanal zu berechnen.Mvec(ρ0)=vec(ρτ)M

Nebenfrage: Ist eine Choi-Matrix?M

Schlussbemerkung: Ich habe Pinja die Annahme zuerkannt, da ich das von Pinja vorgeschlagene Papier verwendet habe. Ich habe selbst eine Antwort gegeben, die die Details ausfüllt.

qubyte
quelle
Was meinst du mit "fragliches System ist ein Qubit oder ein Qutrit und ein anderes Qubit oder ein Qutrit." - Was ist das "andere System"? Sprechen Sie über die Ancilla, die für die Implementierung dieses Kanals mit Unitaries + Tracing-Out erforderlich ist? Beachten Sie in diesem Fall, dass die Dimension der Ancilla bis zu D ^ 2 betragen kann, sodass Qubits nicht ausreichen.
Norbert Schuch
Nein, im Moment ist es nur ein Spielzeugmodell, das aus zwei kleinen Quantensystemen besteht, die gekoppelt sind und unterschiedliche T1- und T2-Zeiten haben. Die Antwort auf diese Frage gibt keinen Anlass zu ernster Besorgnis. Es ist eher ein Punkt von Interesse, da es hilfreich sein könnte, in Zukunft mehr darüber zu erfahren, wie dies getan werden kann.
Qubyte
Kann ich diese Frage bitte auf CS-Theorie anstatt auf Physik migrieren lassen?
Qubyte
Naja ... ich denke das wäre hier in Ordnung gewesen, aber okay.
David Z
Vielen Dank. Es tut uns leid, aber ich bin kein großer Fan von Physics.SE, und trotzdem denke ich, dass forschungsorientierte QI-Fragen hier besser passen (nachdem ich überzeugt war).
Qubyte

Antworten:

9

Ich habe in meiner Masterarbeit an einem sehr ähnlichen Problem gearbeitet, in dem ich die nicht-markovsche Dynamik eines getriebenen Qubits in einer dissipativen Umgebung untersucht habe. Mein Interesse war es zu überprüfen, ob die von mir erhaltene Master-Gleichung vollständig positiv ist, aber dies ist nur eine Seite Ihres Problems. Die Frage stellte sich als sehr nicht trivial heraus, wenn keine RWA erstellt wurde, aber ich konnte mit Ref einige Ergebnisse erzielen. [ J Mod. Opt. 54, 1695 (2007) ] und die Tatsache ausnutzen, dass das Qubit schwach an die Umgebung gekoppelt ist. Ich werde meine Trommel schlagen und auch die Ref geben. zu einem Artikel, in dem ich einige dieser Ergebnisse präsentiere, [P. Haikka und S. Maniscalco, Phys. Rev. A 81, 052103 (2010)] , finden Sie es möglicherweise nützlich.


quelle
Ah! Es stellt sich heraus, dass ich mir jetzt seit ein paar Tagen die Andersson-Zeitung ansehe. Es scheint sehr vielversprechend und gibt das konkreteste Rezept. Ich mag es, eine Methode zu haben, um Probleme zu lösen. Um ehrlich zu sein, muss ich ein wenig Zeit finden, um mich wirklich hinzusetzen und mir das anzuschauen. Es ist im Moment eher ein persönliches Projekt.
Qubyte
7

Die als Antwort auf die Quantenmechanik als Markov-Prozess gegebenen Referenzen  - insbesondere die Online-Notizen von Carlton Caves " Vollständig positive Karten, positive Karten und das Lindblad-Formular " - enthalten physikalische Ideen und mathematische Werkzeuge, die bei der Beantwortung der Frage hilfreich sind.

Ein wichtiger Punkt ist mit der spezifischen Frage verbunden, die gestellt wird: "Wie kann ich Kraus-Operatoren für das Operator-Summen-Äquivalent von , die in einer nützlichen Form vorliegen?" Für große Quantensysteme, ein generischer Superoperator M wird nicht eine algorithmisch komprimierbaren Form hat. Darüber hinaus sind Kraus-Darstellungen nicht eindeutig, und meines Wissens nach gibt es kein allgemeines und effizientes Verfahren, um Kraus-Darstellungen eines bestimmten M zu finden , die eine "nützliche Form" haben (nach welchen Kriterien auch immer) sind für ein Formular angegeben, das "nützlich" ist). Diese entscheidende Quantentrennbarkeit ist NP-schwer, was darauf hindeutet, dass kein effizienter allgemeiner Algorithmus zur Repräsentationsfindung existiert, selbst wenn MMMMM ist in seiner Gesamtheit numerisch angegeben.

Um Fortschritte zu erzielen, kann es hilfreich sein, heuristische Fragen zu stellen: "Was ist das Besondere an meinem Superoperator? Kann ich eine Reihe von Lindbladian-Generatoren dafür ausstellen, die nützliche Symmetrieeigenschaften haben und / oder kompatible Druckflüsse auf dem Hilbert-Zustandsraum erzeugen ? Sind diese Lindbladian-Eigenschaften mit einer natürlichen Hilbert-Basis verbunden, in der eine spärliche, faktorisierte oder anderweitig algorithmisch komprimierbare Darstellung hat? "M

Wenn Fragen wie diese durch "Drehen einer algorithmischen Kurbel" effizient beantwortet werden könnten, wäre die Quantenphysik ein weitaus weniger interessantes Thema! :)

John Sidles
quelle
Dies ist so ziemlich das, was ich gehofft hatte, was nicht der Fall war, aber ich dachte es wäre. Leider weist das System nur dann eine ausnutzbare Symmetrie auf, wenn nur dephasiert wird, ohne dass die Bevölkerung entvölkert wird. Es gibt eine sehr ansprechende Form der Lindblad-Meistergleichung, die Terme, die nicht von Krauss stammen, zu einem nicht-hermitischen Hamiltoninan zusammenfasst. Für den Fall, dass das Hamilton keine Zeitabhängigkeit aufweist, kann eine Basis gewählt werden, die natürlich den Zerfall ausdrückt wie die übrigen Krauss-Begriffe. Ordentlich, aber keine Hilfe für mich.
Qubyte
Eine der Referenzen in Caves 'Notizen ist Wolf and Cirac Dividing Quantum Channels (arXiv: math-ph / 0611057), die ich ohne die geringste Garantie empfehle, die (vielen und subtilen) quanteninformatischen Probleme, die in diesem Artikel behandelt werden, persönlich erfasst zu haben! :)
Schön, das schaue ich mir an. Eine interessante Sache, die ich möglicherweise über die nicht angetriebene Version des obigen Systems hätte bemerken müssen, ist, dass die Zeitunabhängigkeit bedeutet, dass Sie direkt mit Matrixexponentiation finden können (nicht effizient, aber dieses System ist klein genug). Sie können auch konstruierenM mit einigen erratenen Krauss-Operatoren, und ein paar Gleichungen später ergeben die Zuordnung zwischen den beiden, wodurch die Fehlerraten auf den verschiedenen Kanälen extrahiert werden können. M
Qubyte
6

Ich denke, was Sie vielleicht suchen, ist Folgendes: Die Real Density Matrix . Es gibt Ihnen ein Rezept für die Konvertierung zwischen verschiedenen Superoperator-Darstellungen (einschließlich der Verwendung einer Tensorproduktbasis von Paulis). Ein detailliertes Experiment zur Quantenprozess-Tomographie unter Verwendung der Ergebnisse finden Sie hier: Quantenprozess-Tomographie der Quanten-Fourier-Transformation . Ganz allgemein hat Havel hier auch Algorithmen abgeleitet, um in minimale Kraus-Darstellungen umzuwandeln: Verfahren zum Umwandeln zwischen Lindblad-, Kraus- und Matrix-Darstellungen von quantendynamischen Halbgruppen .

vec(ρ)ρvec(|ij|)=|i|jvec(|ij|)=|j|icol(ρ)Mrowvec(ρ0)=vec(ρt)Mcolcol(ρ0)=col(ρt)

C=i,j(1|ij|)Mrow(|ij|1),
C=i,j(|ij|1)Mcol(1|ij|).
{|ij||kl|}
Chris Ferrie
quelle
Das ist interessant, es könnte genau das sein, wonach ich suche ...
Qubyte
Ich habe gerade Ihren Zusatz gesehen. Danke, das ist sehr nützlich. Ich habe ursprünglich deine Version von vec genommen, aber jetzt verwende ich die gestapelten Spalten. Danke Wikipedia dafür. Vielleicht sollte ich Ihre Notation zur Klarheit übernehmen.
Qubyte
4

Wie Pinja bemerkte, wurde in einem Aufsatz von Andersson et al. ( arXiv ) ( DOI ) war besonders nützlich. Das Papier geht sehr ins Detail, und ich habe mich heute endlich hingesetzt, um es mir genauer anzusehen. Als Beispiel für ein Problem habe ich zwei Qubits mit einer Austauschinteraktion ausgewählt, um dies zu überprüfen. Dies ist eine minimale Version dessen, was ich in Betracht ziehe. Zunächst ist die Hauptgleichung gegeben durch

ρ˙=Λ(ρ).

Das Verfahren erfordert, dass Basisoperatoren des Systems ausgewählt werden. Es ist zweckmäßig, diese in Form der Pauli-Matrizen für zwei Qubits anzugeben, aber für einen Qutrit würde man die Gell-Mann-Matrizen verwenden. Definierenσich=1,σx,σy,σz Dieses System hat für jedes Qubit eine Basis aufgebaut aus deren Tensorprodukten mit einem Faktor von 1/2 zur Normalisierung ergeben sich 16 Operatoren Gich z.B G5=Gxx=(σxσx)/2. Sticking with Hermitian operators keeps things neat as well, since some daggers can be neglected.

A special matrix is now composed called L, which is related to the master equation.

Ln,m=Tr[GnΛ(Gm)].

If we are dealing with the master equation as a matrix acting on a vectorised density operator as discussed in the question, then this can be expressed as

Ln,m=vec(Gn)Λvec(Gm),

which allows L to be derived in a single matrix equation, but that's getting a little off topic.

In the sample case I considered, L is does not contain time varying terms, so it may be exponentiated to get a new matrix F, which is related to the solution of the master equation ϕ

F(t)=exp(Lt).

F can be used to get a Choi matrix S, which is exactly what I need. At this point, a basis needs to be chosen for the future Krauss operators. I'm quite happy with the Pauli operators so I'll stick with those for this next equation,

Sa,b=n,mFm,nTr[GnGaGsGb].

Finally, the wonderful part.

ρt=ϕn,m(ρ0,t)=Sn,m(t)Gnρ0Gm

As you can see, S is a matrix of weights for a sum of superoperators in a useful basis that I can select. This has been referred to as the process matrix (arXiv)(DOI) which is unique to a process in a given basis. In the sample case, in which the master equation has no time dependent terms on the RHS, the solution can be directly verified by representing Λ in matrix form and exponentiating it to get ϕ(t)=exp(Λt).

This works in the time independent case for quits and qutrits as expected. I need to check that this works in the case of time dependence.

qubyte
quelle