Limites de velocidade explícitos de Lieb-Robinson

22

Os limites de Lieb-Robinson descrevem como os efeitos são propagados através de um sistema devido a um hamiltoniano local. Eles são frequentemente descritos no formato que e são operadores que são separados por uma distância em uma rede onde o Hamiltoniano tem (por exemplo, vizinho mais próximo) interações locais nessa estrutura, delimitadas por alguma força . As provas do limite de Lieb Robinson mostram tipicamente a existência de uma velocidade (que depende de ). Isso geralmente é realmente útil para delimitar propriedades nesses sistemas. Por exemplo, houve alguns resultados muito bons aqui

|[A,B(t)]|Cevtl,
ABlJvJ quanto tempo leva para gerar um estado GHZ usando um Hamiltoniano do vizinho mais próximo.

O problema que tive é que as provas são suficientemente genéricas e é difícil obter um valor apertado sobre qual é realmente a velocidade de um determinado sistema.

Para ser específico, imagine uma cadeia unidimensional de qubits acoplada por um Hamiltoniano que para todos os . Aqui , e representam um operador de Pauli sendo aplicada a um determinado qubit e em qualquer outro lugar. Você pode dar um limite superior bom (ou seja, o mais rígido possível) para a velocidade de Lieb-Robinson para o sistema na Eq. (1)

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

Esta pergunta pode ser feita sob duas suposições diferentes:

  • O e o são todos fixados no tempoJnBn
  • O e o podem variar no tempo.JnBn

O primeiro é um pressuposto mais forte que pode facilitar as provas, enquanto o segundo geralmente é incluído na declaração dos limites de Lieb-Robinson.


Motivação

A computação quântica e, em geral, a informação quântica, resume-se a criar estados quânticos interessantes. Através de trabalhos como esse , vemos que a informação leva um certo tempo para se propagar de um lugar para outro em um sistema quântico em evolução devido a um hamiltoniano como na Eq. (1), e que estados quânticos, como estados GHZ ou estados com uma ordem topológica, levam um certo tempo para serem produzidos. O que o resultado mostra atualmente é uma relação de escala, por exemplo, o tempo necessário é .Ω(N)

Então, digamos que eu vir para cima com um esquema que faz transferência de informação, ou produz um estado etc. GHZ de uma forma que as escalas linearmente em . Quão bom é esse esquema, na verdade? Se eu tiver uma velocidade explícita, posso ver como o coeficiente de escala está muito próximo do meu esquema em comparação com o limite inferior.N

Se eu acho que um dia o que eu quero ver é um protocolo implementado no laboratório, então eu me preocupo muito em otimizar esses coeficientes de escala, não apenas a ampla funcionalidade de escala, porque quanto mais rápido eu puder implementar um protocolo, menos chances haverá. é para o barulho aparecer e estragar tudo.


Outras informações

Existem algumas características interessantes deste Hamiltoniano que, presumo, facilitam o cálculo. Em particular, o hamiltoniano tem uma estrutura de subespaço baseada no número de 1s na base padrão (diz-se preservar a excitação) e, melhor ainda, a transformação de Jordan-Wigner mostra que todas as propriedades de subespaços de excitação mais altos podem ser derivadas do subespaço de 1 excitação.N×Nh2N×2NH, em que Há alguma evidência de que a velocidade de Lieb-Robinson é , como aqui e aqui , mas todas elas usam uma cadeia próxima a uniformemente acoplada, que tem uma velocidade de grupo (e presumo que a velocidade de grupo esteja intimamente conectada à Velocidade de Lieb-Robinson). Isso não prova que todas as opções possíveis de força de acoplamento têm uma velocidade que é tão limitada.

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

Eu posso acrescentar um pouco mais à motivação. Considere a evolução no tempo de uma única excitação começando em uma extremidade da cadeia, , e qual é a sua amplitude para chegar à outra extremidade da cadeia , pouco tempo mais tarde. Para a primeira ordem em , este é Você pode ver a funcionalidade exponencial que você esperaria estar fora do 'cone de luz' definido por um sistema Lieb-Robinson, mas mais importante, se você quisesse maximizar essa amplitude, todo o|1|Nδtδt

N|eihδt|1=δtN1(N1)!n=1N1Jn+O(δtN).
Jn=J. Assim, em momentos curtos, o sistema uniformemente acoplado leva à transferência mais rápida. Tentando forçar isso ainda mais, você pode perguntar, como se fosse um fudge, quando pode Pegar o limite grande e usar a fórmula de Stirling no fatorial leva a que sugere uma velocidade máxima de . Fechar, mas dificilmente rigoroso (como os termos de ordem superior não são desprezíveis)!
tN1(N1)!n=1N1Jn1
N
etJN11,
eJ

DaftWullie
fonte
Você calculou o melhor limite de LR das provas para esse modelo? Como ele se compara à velocidade que você cita?
Norbert Schuch
1
Ok, admito que é uma questão de computação quântica, pelo menos da maneira como a interpreto agora: "Qual é a escolha de e (sujeito a algumas restrições) que produz a velocidade máxima para a transferência de informações / estado / .... " --- Essa é a interpretação correta? JnBn
Norbert Schuch
@NorbertSchuch Não é bem assim. Quero poder dizer: "Eu criei um conjunto de acoplamentos que atingem um protocolo com uma certa escala. Esse protocolo é conhecido por ser limitado pelos limites de Lieb-Robinson. Quão perto estou de saturar essa restrição?" como uma medida de quão rápido é o meu protocolo.
DaftWullie #
@DaftWullie Então - você pergunta: "Quão perto estou de ser o ideal" ou "Quão perto estou de algum tipo de limite (assumindo o mais rígido possível)"?
Norbert Schuch
1
@ user1271772 Isso está correto. B(t)=eiHtB(0)eiHt
iHt

Respostas:

4

Deixe-me primeiro responder à pergunta geral sobre como obter uma velocidade Lieb-Robinson (LR) razoavelmente reduzida quando você estiver enfrentando um modelo de treliça de interação local genérica e depois voltarei ao modelo 1D XY na sua pergunta, o que é muito especial para ser exatamente solucionável.


Método Geral

O método para obter o limite mais restrito até a data (para um modelo de interação genérico de curto alcance) é introduzido em Ref1 = arXiv: 1908.03997 . A idéia básica é que a norma do comutador de tempo desigualentre operadores locais arbitrários pode ser delimitado pela solução para um conjunto de equações diferenciais lineares de primeira ordem que vivem no gráfico de comutatividade do modelo. O gráfico de comutatividade, conforme apresentado na Seção II A da Ref1, pode ser facilmente desenhado a partir do modelo Hamiltoniano e foi projetado para refletir as relações de comutação entre diferentes operadores locais apresentados em[AX(t),BY(0)]H^H^. Em sistemas invariantes de conversão, esse conjunto de equações diferenciais pode ser facilmente resolvido por uma transformada de Fourier, e um limite superior da velocidade LR pode ser calculado a partir da maior frequência de autenticidade usando Eq. (31) da Ref . A seguir, aplicarei esse método ao modelo 1D XY como um exemplo pedagógico. Por simplicidade, vou me concentrar no caso invariante e independente de tempo , (o limite resultante não depende dos sinais deωmax(iκ)|Bn|=B>0|Jn|=J>0Bn,Jn) Para o caso não invariável e dependente do tempo da tradução, é possível resolver numericamente a equação diferencial (que é uma tarefa computacional fácil para sistemas de milhares de sites) ou usar um limite superior geral e continue usando o método abaixo (mas isso compromete levemente a tensão em comparação com o método numérico).|Jn(t)|J, |Bn(t)|B

  1. Primeiro, desenhamos o gráfico da comutatividade, como abaixo. Cada operador no Hamiltoniano ~ ( , , ) é representado por um vértice e vinculamos dois vértices se e somente se os operadores correspondentes não comutarem ( ou, no caso atual, anti-deslocamento). XnXn+1YnYn+1Zninsira a descrição da imagem aqui

  2. Em seguida, anote as equações diferenciais Eq. (10) da 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 transformando a equação acima, temos As freqüências próprias são . A velocidade LR é dada pela Eq. (31) da Ref1 : em que

    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)
    vLRminκ>0ωmax(iκ)κ=ZB2JJ,
    Zyminκ>0coshκ+cosh2κ+4y(1+coshκ)κ.

Nota: Esse limite diverge quando , enquanto a velocidade de propagação da informação física permanece finita. Podemos nos livrar desse problema usando o método no segundo. VI da Ref1 . O resultado é neste limite, onde é definido como a solução para a equação .B/JvLR4X0JXyxarcsinh(x)=x2+1+y


Limites de velocidade para alguns modelos clássicos

O método acima é completamente geral. Caso você esteja interessado em mais, listei os limites de velocidade para alguns modelos clássicos na tabela a seguir, obtidos de maneira semelhante. Observe que a velocidade LR é limitada pela menor de todas as expressões listadas (portanto, em diferentes regiões de parâmetros, expressões diferentes devem ser usadas). A função é definida como a maior raiz deTodos os parâmetros são considerados positivos (basta assumir o valor absoluto para os casos negativos).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}

Quanto à qualidade desses limites, não investiguei em geral, mas para o 1D TFIM no ponto crítico , a solução exata fornece , enquanto o limite acima fornece . Da mesma forma, no ponto de FH e ponto de Heisenberg XYZ, o limite acima é todo maior que a solução exata por um fator de . [Na verdade, nesses pontos especiais, as duas últimas são equivalentes a cadeias dissociadas de TFIM, como pode ser julgado diretamente pelo seu gráfico de comutatividade.]J=hvLR=2J2X0J3.02JU=0Jx=Jy,Jz=0X01.50888


Limite mais apertado para 1D XY, mapeando para férmions livres

Agora vamos falar mais sobre o modelo 1D XY. Como você notou, é exatamente solucionável ao mapear para liberar férmions: Para geral é necessário resolver numericamente o problema de fermion livre, mas deixe-me mencionar dois casos especiais analiticamente tratáveis.

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

  1. Bn(t)=B,Jn(t)=J são fixos e a tradução é invariável. Então a solução exata é que é a função de Bessel da ordem. Portanto, a velocidade da LR é .

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

  2. Bn,Jn são fixados no tempo, mas são completamente aleatórios (desordem extinta). Então, devido à localização de muitos corpos (ou localização de Anderson na figura do férmion), as informações não se propagam neste sistema, portanto . Mais rigorosamente, em arXiv: quant-ph / 0703209 , o seguinte limite é comprovado para casos desordenados: com um cone de luz logarítmico em desaceleração .vLR=0

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

Lagrenge
fonte
Devo inferir do que você diz que para todo modelo (incluindo aqueles sem invariância de conversão) com , que a velocidade é ? XY|Jn|JvLRXY2J
DaftWullie 25/08
@DaftWullie Não, você só pode usar um limite superior geral para os parâmetros no método geral, pois o método geral sempre fornece um limite estritamente não decrescente no valor absoluto de qualquer coeficiente. O limite é obtido da solução exata de férmion livre, na qual você não pode usar um limite superior geral para parâmetros e precisa resolver caso a caso. Se for invariável na conversão, você poderá definir no método geral, pois o termo comuta com e obtém . 2JBn(t)B=0BH^vLR2X0J=3.02J
Lagrenge em
@DaftWullie Caro DaftWullie, se você acha que ainda falta alguma coisa na minha resposta ou se algum ponto ainda não está claro, entre em contato.
Lagrenge
a resposta parece potencialmente útil. Ainda não tive tempo de analisar seu trabalho (pode demorar algumas semanas). Supondo que eu entenda tudo bem, esse é o ponto em que aceitarei sua resposta.
DaftWullie 27/08