Limites explicites de vitesse de Lieb-Robinson

22

Les limites de Lieb-Robinson décrivent comment les effets se propagent à travers un système en raison d'un hamiltonien local. Ils sont souvent décrits sous la forme où et sont des opérateurs séparés par une distance sur un réseau où le hamiltonien a des interactions locales (par exemple , du plus proche voisin) sur ce réseau, délimité par une certaine force . Les preuves de la borne de Lieb Robinson montrent généralement l'existence d'une vitesse (qui dépend de ). Ceci est souvent très utile pour délimiter des propriétés dans ces systèmes. Par exemple, il y a eu de très bons résultats ici

|[A,B(t)]|Cevtl,
ABlJv JvJ concernant le temps qu'il faut pour générer un état GHZ en utilisant un hamiltonien le plus proche voisin.

Le problème que j'ai eu est que les preuves sont suffisamment génériques qu'il est difficile d'obtenir une valeur sur ce que la serré vitesse fait est pour un système donné.

Pour être précis, imaginez une chaîne unidimensionnelle de qubits couplée par un hamiltonien où pour tout . Ici , Y_n et Z_n représentent un opérateur Pauli appliqué à un qubit n donné , et \ mathbb {I} partout ailleurs. Pouvez-vous donner une bonne limite supérieure (c'est-à-dire aussi serrée que possible) pour la vitesse de Lieb-Robinson v pour le système dans l'équation. (1)?

(1)H=n=1NBn2Zn+n=1N1Jn2(XnXn+1+YnYn+1),
JnJnXnYnZnnIv

Cette question peut être posée sous deux hypothèses différentes:

  • Le Jn et Bn sont tous fixés dans le temps
  • Le Jn et le Bn peuvent varier dans le temps.

La première est une hypothèse plus forte qui peut rendre les preuves plus faciles, tandis que la seconde est généralement incluse dans la déclaration des limites de Lieb-Robinson.


Motivation

Le calcul quantique, et plus généralement l'information quantique, revient à faire des états quantiques intéressants. À travers des travaux comme celui- ci , nous voyons que l'information prend un certain temps pour se propager d'un endroit à un autre dans un système quantique en cours d'évolution en raison d'un hamiltonien comme dans l'équation. (1), et que les états quantiques, tels que les états GHZ, ou les états avec un ordre topologique, prennent un certain temps à produire. Ce que le résultat montre actuellement est une relation d'échelle, par exemple le temps requis est Ω(N) .

Alors, disons que je viens avec un système qui fait le transfert d'information ou produit un état GHZ etc. d'une manière linéaire dans des échelles . Quelle est la qualité de ce régime en fait? Si j'ai une vitesse explicite, je peux voir à quel point le coefficient de mise à l'échelle est proche dans mon schéma par rapport à la limite inférieure.N

Si je pense qu'un jour ce que je veux voir est un protocole implémenté en laboratoire, alors je me soucie beaucoup d'optimiser ces coefficients de mise à l'échelle, pas seulement la fonctionnalité de mise à l'échelle large, car plus je peux implémenter un protocole rapidement, moins il y a de chance est pour que le bruit arrive et gâche tout.


Plus d'informations

Il y a quelques fonctionnalités intéressantes de cet hamiltonien qui, je suppose, facilitent le calcul. En particulier, l'hamiltonien a une structure de sous-espace basée sur le nombre de 1 dans la base standard (on dit qu'il préserve l'excitation) et, mieux encore, la transformation de Jordan-Wigner montre que toutes les propriétés des sous-espaces d'excitation supérieure peuvent être dérivées du sous-espace 1-excitation.N×Nh2N×2NH, où Il existe des preuves que la vitesse de Lieb-Robinson est , comme ici et ici , mais elles utilisent toutes une chaîne proche à uniformément couplée, qui a une vitesse de groupe (et je suppose que la vitesse de groupe est étroitement liée à la Vitesse de Lieb-Robinson). Cela ne prouve pas que tous les choix possibles de force de couplage ont une vitesse qui est ainsi limitée.

h=n=1NBn|nn|+n=1N1Jn(|nn+1|+|n+1n|).
v=2J2J

Je peux ajouter un peu plus à la motivation. Considérons l'évolution temporelle d'une seule excitation commençant à une extrémité de la chaîne, , et quelle est son amplitude pour arriver à l'autre extrémité de la chaîne , un court instant plus tard. Pour une première commande dans , c'est Vous pouvez voir la fonctionnalité exponentielle à laquelle vous vous attendriez en dehors du `` cône de lumière '' défini par un système Lieb-Robinson, mais plus important encore, si vous vouliez maximiser cette amplitude, vous tous les|1|Nδtδt

N|eihδt|1=δtN1(N1)!n=1N1Jn+O(δtN).
Jn=J. Ainsi, en peu de temps, le système à couplage uniforme conduit au transfert le plus rapide. En essayant de pousser plus loin, vous pouvez demander, comme un peu de fudge, quand peut-on Prendre la grande limite et utiliser la formule de Stirling sur la factorielle conduit à ce qui suggère une vitesse maximale d'environ . Fermer, mais peu rigoureux (car les termes d'ordre supérieur ne sont pas négligeables)!
tN1(N1)!n=1N1Jn1
N
etJN11,
eJ

DaftWullie
la source
Avez-vous calculé la meilleure limite LR à partir des preuves pour ce modèle? Comment cela se compare-t-il à la vitesse que vous citez?
Norbert Schuch
1
Ok, je concède que c'est une question d'informatique quantique, du moins la façon dont je l'interprète maintenant: "Quel est le choix de et (sous réserve de certaines contraintes) qui donne la vitesse maximale pour le transfert d'informations / état / .... " --- Est-ce la bonne interprétation? B nJnBn
Norbert Schuch
@NorbertSchuch Pas tout à fait. Je veux pouvoir dire "J'ai trouvé un ensemble de couplages qui réalise un protocole avec une certaine mise à l'échelle. Ce protocole est connu pour être contraint par les limites de Lieb-Robinson. Dans quelle mesure suis-je proche de saturer cette contrainte?" pour mesurer la vitesse de mon protocole.
DaftWullie
@DaftWullie Alors - est-ce que vous vous demandez: "Dans quelle mesure suis-je proche d'être optimal", ou "Dans quelle mesure suis-je proche d'une sorte de limite (en prenant la plus étroite possible)"?
Norbert Schuch, le
1
@ user1271772 C'est correct. B(t)=eiHtB(0)eiHt
DaftWullie

Réponses:

4

Permettez-moi d'abord de répondre à la question générale comment obtenir une vitesse de Lieb-Robinson (LR) raisonnablement serrée lorsque vous faites face à un modèle de réseau générique à interaction locale, puis je reviendrai sur le modèle 1D XY dans votre question, qui est très spécial pour être exactement soluble.


Méthode générale

La méthode pour obtenir la limite la plus stricte à ce jour (pour un modèle interactif générique à courte portée) est présentée dans Ref1 = arXiv: 1908.03997 . L'idée de base est que la norme du commutateur à temps inégauxentre des opérateurs locaux arbitraires peut être limité par la solution d'un ensemble d'équations différentielles linéaires du premier ordre vivant sur le graphe de commutativité du modèle. Le graphe de commutativité, tel qu'introduit dans la section II A de Ref1, peut être facilement tiré du modèle hamiltonien , et est conçu pour refléter les relations de commutation entre les différents opérateurs locaux présentés dans[AX(t),BY(0)]H H ω max ( i → la K ) | B n | = B > 0 | J n | = J > 0 B n , J n | J n ( t ) | J , | B n (H^H^. Dans les systèmes invariants de translation, cet ensemble d'équations différentielles peut être facilement résolu par une transformée de Fourier, et une limite supérieure de la vitesse LR peut être calculée à partir de la plus grande fréquence propre utilisant Eq. (31) de Ref1 . Dans la suite, j'appliquerai cette méthode au modèle 1D XY comme exemple pédagogique. Par souci de simplicité, je me concentrerai sur le cas indépendant du temps et invariant de la traduction , (la limite résultante ne dépend pas des signes deωmax(iκ)|Bn|=B>0|Jn|=J>0Bn,Jn). Pour le cas de traduction non invariable et dépendant du temps, vous pouvez soit résoudre l'équation différentielle numériquement (ce qui est une tâche de calcul facile pour des systèmes de milliers de sites), soit utiliser une borne supérieure globale et continuez à utiliser la méthode ci-dessous (mais cela compromet légèrement l'étanchéité par rapport à la méthode numérique).|Jn(t)|J, |Bn(t)|B

  1. Tout d'abord, nous dessinons le graphique de commutativité, comme ci-dessous. Chaque opérateur de l'hamiltonien ~ ( , , ) est représenté par un sommet, et nous deux sommets si et seulement si les opérateurs correspondants ne font pas la navette ( ou, dans le cas présent, anti-navettage). XnXn+1YnYn+1Znentrez la description de l'image ici

  2. Notez ensuite les équations différentielles Eq. (10) de Ref1 :

    γ¯˙α,n=J[γ¯α,n1(t)+γ¯α,n+1(t)]+B[γ¯3,n(t)+γ¯3,n+1(t)],  α=1,2,γ¯˙3,n=Jα=1,2[γ¯α,n1(t)+γ¯α,n(t)].

  3. Fourier transformant l'équation ci-dessus, nous avons Les fréquences propres sont . La vitesse LR est donnée par l'équation (31) de Ref1 : où

    ddt(γ¯1,kγ¯2,kγ¯3,k)=(2Jcosk0B(1+eik)02JcoskB(1+eik)J(1+eik)J(1+eik)0)(γ¯1,kγ¯2,kγ¯3,k).
    2Jcosk,Jcosk±(Jcosk)2+2BJ(1+cosk) v LR min κ > 0 ω max ( i κ )
    vLRminκ>0ωmax(iκ)κ=ZB2JJ,
    Zyminκ>0coshκ+cosh2κ+4y(1+coshκ)κ.

Remarque: Cette limite diverge lorsque , tandis que la vitesse de propagation des informations physiques reste limitée. Nous pouvons nous débarrasser de ce problème en utilisant la méthode de la Sec. VI de Ref1 . Le résultat est dans cette limite, où est défini comme la solution de l'équation .B/Jv LR4 X 0 J X y x a r c s i n h ( x ) = vLR4X0JXyxarcsinh(x)=x2+1+y


Limites de vitesse pour certains modèles classiques

La méthode ci-dessus est complètement générale. Si vous souhaitez en savoir plus, j'ai répertorié les limites de vitesse pour certains modèles classiques dans le tableau suivant, obtenues de manière similaire. Notez que la vitesse LR est limitée par la plus petite de toutes les expressions répertoriées (donc dans différentes régions de paramètres, différentes expressions doivent être utilisées). La fonction est définie comme la plus grande racine deTous les paramètres sont supposés positifs (il suffit de prendre la valeur absolue pour les cas négatifs).vLRF(Jx,Jy,Jz)x3(JxJy+JxJz+JyJz)x2JxJyJz=0.

ModelvLRd-dimensional TFIM2X0dJh=3.02dJhH^=JmnXmXn+hnZn4Xd1ddJ8.93dJ4X0dh=6.04dhd-dimensional Fermi-Hubbard2X3U4dJdJH^=Jmn,s=↑,(am,san,s+H.c.)8Xd1ddJ17.9dJ  +UnananananZU/JJ (d=1)1D Heisenberg XYZ4X0F(Jx,Jy,Jz)H^=n(JxXnXn+1+JyYnYn+1+JzZnZn+1)34.6max{Jx,Jy}

Quant à la qualité de ces bornes, je n'ai pas enquêté en général, mais pour le TFIM 1D au point critique , la solution exacte donne , tandis que la borne ci-dessus donne . De même, au point de FH et au point de Heisenberg XYZ, les limites ci-dessus sont toutes plus grandes que la solution exacte d'un facteur . [En fait, à ces points spéciaux, les deux derniers sont équivalents à des chaînes découplées de TFIM, comme on peut le juger directement à partir de leur graphe de commutativité.]J=hvLR=2J2X0J3.02JU=0Jx=Jy,Jz=0X01.50888


Limite plus serrée pour 1D XY par mappage sur des fermions libres

Parlons maintenant davantage du modèle 1D XY. Comme vous l'avez remarqué, il est exactement résoluble en mappant sur des fermions libres: Pour le général vous devez résoudre le problème des fermions libres numériquement, mais permettez-moi de mentionner deux cas spéciaux qui sont analytiquement traitables.

H^=nBn(anan1/2)+nJn(anan+1+H.c.).
Bn(t),Jn(t)

  1. Bn(t)=B,Jn(t)=J sont fixes et invariants en translation. La solution exacte est alors où est la fonction de Bessel d'ordre. La vitesse LR est donc .

    an(t)=12πππa~kei2Jtcoskeikxdk=mJ|nm|(2Jt)am(0),
    J|nm|(2Jt)|nm|vLRXY=2J

  2. Bn,Jn sont fixes dans le temps mais sont complètement aléatoires (trouble éteint). Ensuite, en raison de la localisation à plusieurs corps (ou de la localisation d'Anderson dans l'image du fermion), les informations ne se propagent pas dans ce système, donc . Plus rigoureusement, dans arXiv: quant-ph / 0703209 , la borne suivante est prouvée pour les cas désordonnés: avec un cône lumineux logarithmique décélérant .vLR=0

    [AX(t),BY(0)]const. t edXY/ξ,
    dXY=ξlnt

Lagrenge
la source
Dois-je déduire de ce que vous dites que pour chaque modèle (y compris ceux sans invariance de translation) avec , que la vitesse est ? XY|Jn|JvLRXY2J
DaftWullie
@DaftWullie Non, vous ne pouvez utiliser qu'une limite supérieure globale pour les paramètres de la méthode générale, car la méthode générale donne toujours une limite strictement non décroissante dans la valeur absolue de tout coefficient. La borne est obtenue à partir de la solution exacte de fermion libre, dans laquelle vous ne pouvez pas utiliser une borne supérieure globale pour les paramètres et devez résoudre au cas par cas. Si le est invariant par rapport à la traduction, alors vous pouvez définir dans la méthode générale puisque le terme commute avec , et obtenez . B n ( t ) B = 0 B H v LR2 X 0 J = 3,02 J2JBn(t)B=0BH^vLR2X0J=3.02J
Lagrenge
@DaftWullie Cher DaftWullie, si vous pensez que quelque chose manque encore dans ma réponse, ou si un point n'est toujours pas clair, faites-le moi savoir.
Lagrenge
la réponse semble potentiellement utile. Je n'ai pas encore eu le temps de regarder votre article (cela peut prendre quelques semaines). En supposant que je comprends tout OK, c'est le point que j'accepterai votre réponse.
DaftWullie