TensorMixedStates: a Julia library for simulating pure and mixed quantum states using matrix product states

Resumo Técnico de "TensorMixedStates: uma biblioteca Julia para simular estados quânticos puros e mistos usando estados de produto matricial"

Enunciação do Problema
A simulação de sistemas quânticos de muitos corpos abertos é um desafio crítico na física moderna, impulsionado pelo progresso experimental em átomos frios, íons presos e circuitos supercondutores. Diferentemente de sistemas fechados, sistemas abertos interagem com seu ambiente, levando a ruído, dissipação e decoerência, que são frequentemente descritos pela equação mestra de Lindblad. Embora métodos de rede de tensores, especificamente Estados de Produto Matricial (MPS), tenham se tornado o padrão para simular estados puros (usando bibliotecas como ITensor ou TenPy), o software para simular estados mistos (matrizes densidade) dentro do framework MPS é significativamente mais limitado. Soluções existentes frequentemente carecem de flexibilidade, são restritas a espaços de Hilbert locais específicos (por exemplo, apenas qubits) ou não se integram perfeitamente com algoritmos de última geração. Além disso, representar estados mistos exige tratar matrizes densidade como vetores em um espaço de Hilbert duplicado, introduzindo complexidades nas ações de operadores (por exemplo, $O\rho$, $O\rho O^\dagger$, comutadores) que frameworks de estado puro padrão não suportam nativamente.

Metodologia
Os autores introduzem o TensorMixedStates (TMS), uma biblioteca Julia construída sobre a biblioteca ITensor. O TMS foi projetado para simular tanto estados quânticos puros quanto mistos usando representações MPS.

• Arquitetura Central: O TMS aproveita a manipulação robusta de tensores de baixo nível e os algoritmos de alto desempenho (DMRG, TDVP) do ITensor. No entanto, reconhecendo que o framework nativo do ITensor é otimizado para estados puros, o TMS implementa um novo framework de estado e operador do zero para acomodar estados mistos.
• Representação de Estado Misto: A biblioteca representa a matriz densidade $\rho$ como um estado vetorizado $|\rho\rangle\rangle$ em um espaço maior de dimensão $d^2$. Isso permite o uso de técnicas MPS onde a matriz densidade é decomposta em um produto de tensores de 4 índices (ou tensores de 3 índices na base vetorizada).
• Framework de Operadores: Uma escolha de design chave é o desenvolvimento de um framework de operadores flexível capaz de lidar com as quatro maneiras distintas pelas quais um operador pode atuar sobre uma matriz densidade: $O\rho$, $O\rho O^\dagger$, $[O, \rho]$ e $\{O^\dagger O, \rho\}$. Isso inclui suporte a dissipadores de Lindblad e canais quânticos gerais (operadores de Kraus).
• Algoritmos:
  • Evolução Temporal: O TMS suporta evolução temporal contínua via equação de Lindblad usando dois algoritmos primários: TDVP (Princípio Variacional Dependente do Tempo) e aproximações WI/WII (aproximações MPO do mapa exponencial). Os métodos WI/WII são notados como complementares ao TDVP, particularmente para evoluções não hermitianas.
  • Estados Estacionários: A biblioteca calcula estados estacionários encontrando o estado fundamental do operador hermitiano $L^\dagger L$ usando DMRG, evitando a necessidade de integração de longo tempo.
  • Evolução Discreta: Suporta a aplicação de portas não unitárias e canais quânticos.
• Interface: O TMS oferece uma interface dupla: uma API de baixo nível para controle granular e uma interface de alto nível (runTMS) que permite aos usuários definir fluxos de trabalho de simulação complexos (criação de estado, evolução, medição) em poucas linhas de código.

Contribuições Principais

1. Simulação Generalizada de Estados Mistos: O TMS estende as capacidades MPS além de qubits para espaços de Hilbert locais gerais, incluindo bósons, férmions, spins e modelos $t-J$, abordando uma lacuna no ecossistema de software para sistemas quânticos abertos.
2. Framework de Operadores Flexível: Fornece um conjunto abrangente de ferramentas para construir Hamiltonianos, dissipadores de Lindblad e canais quânticos gerais, lidando com a complexidade algébrica das operações de estado misto.
3. Suíte Algorítmica: A biblioteca integra TDVP e aproximações WI/WII para evolução temporal e DMRG para cálculos de estado estacionário especificamente adaptados ao formalismo de Lindbladiano.
4. Design Amigável ao Usuário: A interface de alto nível simplifica a configuração de simulações sofisticadas, enquanto a arquitetura subjacente Julia/ITensor garante multithreading e contrações de tensores eficientes.
5. Validação: A biblioteca é validada contra soluções analíticas exatas para vários modelos quase-livres (cadeias fermiônicas e bosônicas com desfaseamento ou injeção, cadeias de spin XX com dissipação de fronteira) e resultados conhecidos para decoerência de estados de grafo e circuitos quânticos ruidosos.

Resultados
O artigo demonstra as capacidades da biblioteca através de seis exemplos:

• Cadeia Fermiônica com Desfaseamento: Simulações da densidade de férmions e Entropia de Entrelaçamento do Espaço de Operadores (OSEE) mostram excelente concordância com resultados assintóticos exatos para várias intensidades de dissipação.
• Cadeia de Spin XX com Dissipação de Fronteira: A evolução temporal da magnetização e da corrente de spin coincide com resultados da literatura anterior.
• Bósons Livres com Fonte Localizada: A biblioteca modela com sucesso a dinâmica de uma cadeia bosônica com injeção de partículas, capturando perfis de densidade e crescimento do número total de partículas, embora o custo computacional aumente com a dimensão do espaço de Hilbert local.
• Férmions Livres com Fonte Localizada: Comparações entre os algoritmos TDVP e WI/WII mostram que ambos oferecem precisão similar para este modelo, embora o TDVP requira significativamente mais tempo de CPU.
• Decoerência de um Estado de Grafo Completo: A biblioteca lida com tamanhos de sistema grandes (até $N=512$) de forma eficiente, pois o baixo entrelaçamento do estado permite dimensões de ligação pequenas.
• Circuito Quântico Ruidoso: A biblioteca simula um circuito quântico em formato de parede de tijolos com portas unitárias e ruído de despolarização, rastreando o crescimento e subsequente decaimento da OSEE devido ao ruído.

Significado e Alegações
Os autores posicionam o TMS como uma ferramenta necessária para preencher a lacuna em software para simular sistemas quânticos abertos com MPS. Eles afirmam que a biblioteca oferece uma "interface amigável ao usuário" permitindo simulações sofisticadas em "poucas linhas de código", mantendo ao mesmo tempo o acesso a "algoritmos de última geração" via ITensor. O artigo enfatiza que a biblioteca não implementa desenovelamento de trajetórias, mas foca na evolução direta da matriz densidade. O significado reside em sua capacidade de lidar com espaços de Hilbert locais gerais (diferentemente de seu predecessor, Lindbladmpo) e sua integração com o ecossistema moderno Julia. Os autores notam modestamente que, embora a biblioteca seja poderosa, o custo computacional para sistemas com grandes dimensões locais (como bósons) permanece alto, e trabalhos futuros focarão na incorporação de números quânticos conservados (QNS) e compressão automática de MPO para melhorar ainda mais a eficiência.