Optimizing Quantum Photonic Integrated Circuits using Differentiable Tensor Networks

Resumo Técnico: Otimização de Circuitos Integrados Fotônicos Quânticos usando Redes de Tensores Diferenciáveis

Declaração do Problema
Avanços recentes na fotônica integrada demonstraram grandes não linearidades ópticas em dispositivos semicondutores (especificamente heteroestruturas planares baseadas em GaAs) via forte acoplamento luz-matéria excitônico. Essas capacidades permitem o processamento quântico no regime de baixa ocupação de fótons, onde estatísticas de fótons não clássicas podem superar contrapartes clássicas em computação e sensoriamento. No entanto, o projeto de Circuitos Integrados Fotônicos Quânticos (qPICs) para esses regimes permanece desafiador. Diferentemente de circuitos lineares clássicos, os qPICs envolvem portas unitárias não lineares e componentes estocásticos não unitários (perdas) que quebram estatísticas de Poisson coerentes. As estruturas de otimização existentes para PICs clássicos não capturam as estatísticas quânticas relevantes exigidas para qPICs. Além disso, descrever estados quânticos fotônicos tipicamente requer um espaço de Hilbert que cresce exponencialmente, tornando a otimização sistemática difícil.

Metodologia
Os autores propõem uma estrutura de otimização baseada em gradiente para qPICs que alavanca redes de tensores diferenciáveis, especificamente a representação de Estado de Produto Matricial (MPS). A metodologia integra três componentes principais:

1. Caracterização de Portas: O circuito físico é modelado com base em simulações de campo 2D de amostras de GaAs gravadas. O qPIC consiste em camadas empilhadas de portas de acoplamento não linear de dois modos. Cada porta é definida por um Hamiltoniano $H^{(2)}(J, U)$ contendo uma taxa de tunelamento $J$ (acoplamento linear) e uma não linearidade intra-guia de onda $U$ (tipo Kerr). As portas são parametrizadas por valores adimensionais $(J\Delta t, U\Delta t)$, derivados de parâmetros experimentais realistas incluindo a velocidade de grupo de polaritons e frações de excitons.
2. Simulação MPS Diferenciável: O estado de fótons bosônico multimodo é representado como um MPS, onde cada modo é modelado por um espaço de Fock local truncado. Esta abordagem é válida no regime de baixa ocupação de fótons com emaranhamento intermodo modesto. A contração do circuito envolve a aplicação de portas unitárias aos tensores MPS. Crucialmente, os autores implementam diferenciação automática de contrações de tensores (usando PyTorch) para rastrear gradientes em relação aos parâmetros de acoplamento $J_{l,d}$. Para lidar com instabilidades numéricas na Decomposição em Valores Singulares (SVD) durante a retropropagação, eles empregam truncamento explícito e regularização de valores singulares quase degenerados.
3. Amostragem Estocástica de Perdas: Para contabilizar perdas fotônicas e excitônicas (mecanismos de dissipação dominantes), o método incorpora amostragem de trajetórias quânticas Monte Carlo (MC). Os tensores MPS são aumentados com um índice de lote para representar um conjunto de trajetórias. Portas estocásticas não unitárias são aplicadas para simular eventos de perda de fótons (operadores de salto) ou passos de tempo não unitários (eventos de não clique). Os valores esperados de observáveis são reconstruídos a partir deste conjunto, mantendo a diferenciabilidade das operações de tensor necessárias para a otimização baseada em gradiente.

Principais Contribuições
O artigo contribui com uma estrutura abrangente para o projeto e otimização de qPICs:

• Estrutura de Projeto: Uma descrição e parametrização consolidadas de portas de circuito baseadas em simulações de campo 2D de amostras de GaAs.
• Estrutura de Tensores Diferenciáveis: Desenvolvimento de um repositório de código personalizado para simulação e otimização quântica fotônica. Esta estrutura suporta exclusivamente simulação quântica bosônica, otimização dentro da variedade de portas não lineares de Kerr, inclusão de ruído fotônico via amostragem estocástica e tarefas específicas de otimização fotônica.
• Otimização de Geração de Estados: Demonstração da otimização de qPICs para preparação de estados quânticos, especificamente geração de estados de gato de Schrödinger de circuitos profundos e geração de fótons únicos ruidosa.
• Otimização de Sensoriamento: Demonstração da otimização de qPICs para leitura de sensoriamento de fase quântica, maximizando a sensibilidade a pequenos deslocamentos de fase em luz de entrada coerente.

Resultados
Os autores validaram o método através de dois casos de uso primários:

1. Geração de Estados Quânticos:

   • Geração de Estados de Gato: O método otimizou com sucesso circuitos para gerar estados de gato de Schrödinger ímpares. Os resultados indicam que, para geração de modo único, interferências sintonizadas e não linearidade de Kerr acumulada usando apenas três guias de onda são suficientes para alcançar alta fidelidade (>90% de fidelidade condicionada). Aumentar o número de guias de onda (por exemplo, para 11) não melhorou significativamente a fidelidade ou a robustez contra erros de fabricação em comparação com o caso de 3 guias de onda. A não linearidade acumulada ($U \cdot T_{circ}$) foi identificada como o fator primário para o sucesso.
   • Geração de Fótons Únicos: Em um ambiente ruidoso, o método otimizou para forte antiagrupamento ($g^{(2)}(0) < 1$) e eficiência de fóton único. Os resultados mostraram que a não linearidade do circuito é o fator dominante para gerar estatísticas antiagrupadas. Embora aumentar a profundidade do circuito ($D$) oferecesse melhorias marginais na probabilidade de geração, não alterou significativamente o grau de antiagrupamento.

2. Sensoriamento de Fase Quântica:

   • A otimização visou maximizar a Informação de Fisher (FI) clássica da distribuição do número de fótons na saída para detectar pequenos deslocamentos de fase na entrada.
   • O estudo encontrou que qPICs profundos com não linearidades modestas ($U\Delta t > 0$) podem superar os limites de sensibilidade da interferometria fotônica linear padrão (e detecção homódina) no regime de baixo número de fótons.
   • Embora interferências lineares ofereçam benefícios significativos sobre a detecção homódina padrão, a inclusão de não linearidade introduz estatísticas de fótons não gaussianas, o que aumenta ainda mais a Informação de Fisher além do que as estatísticas gaussianas sozinhas preveriam.

Significado e Alegações
O artigo alega fornecer uma estrutura prática e escalável para o projeto de circuitos fotônicos quânticos no regime de baixo número de fótons, um domínio que anteriormente carecia de ferramentas de otimização sistemática. Ao combinar representações MPS com programação diferenciável e amostragem estocástica, os autores demonstram que:

• A otimização baseada em gradiente é viável para circuitos fotônicos quânticos complexos e dissipativos.
• A não linearidade acumulada é um recurso crítico para gerar estados não clássicos e melhorar capacidades de sensoriamento, frequentemente mais do que o mero número de guias de onda ou a profundidade do circuito.
• Estatísticas não gaussianas decorrentes de não linearidades de Kerr oferecem uma vantagem genuína para sensoriamento quântico, permitindo que circuitos se aproximem ou excedam limites de detecção clássica em regimes de baixa energia.

Os autores concluem que, embora sua abordagem atual seja limitada pelo truncamento do espaço de Fock e pela suposição de emaranhamento modesto, ela fornece uma linha de base robusta para tarefas de processamento quântico de curto prazo. Eles sugerem que extensões futuras poderiam incluir a otimização simultânea da preparação de estados e da leitura, incorporando robustez explícita contra erros de fabricação e explorando métodos de simulação de variáveis contínuas para ocupações de fótons mais altas. O código para a estrutura desenvolvida está disponível publicamente.