PT-symmetry breaking phase transitions in an LMG dimer

Resumo Técnico: Transições de Fase de Quebra de Simetria PT em um Dímero LMG

Enunciado do Problema
Sistemas quânticos abertos, impulsionados tanto por campos externos quanto por dissipação, exibem uma fenomenologia mais rica do que seus equivalentes em equilíbrio, incluindo a estabilização de fases impossíveis em estados fundamentais. Embora pesquisas anteriores tenham explorado transições de fase dissipativas em modelos como o Dicke aberto e o modelo Lipkin-Meshkov-Glick (LMG) dissipativo, e investigado separadamente transições de quebra de simetria Paridade-Tempo (PT), a interação entre essas distintas classes de transições fora do equilíbrio permanece subexplorada. Especificamente, há uma necessidade de compreender como um sistema que combina um mecanismo de acoplamento PT-simétrico com uma interação LMG não linear e de quebra de simetria se comporta. Os autores visam construir um modelo que integre esses dois mecanismos para investigar a competição, a estrutura de fases e os regimes dinâmicos emergentes resultantes.

Metodologia
Os autores introduzem um modelo consistindo em um par de modelos LMG acoplados: um dissipativo e um impulsionado, ligados via um Hamiltoniano PT-simétrico. O sistema é descrito por uma equação mestre que governa a matriz densidade $\rho$ de dois spins coletivos (rotulados como A e B) de comprimento $S$.

• Construção do Modelo: O Hamiltoniano total $H$ combina os Hamiltonianos LMG para cada spin ($H_A^{LMG} \otimes I_B + I_A \otimes H_B^{LMG}$) com um termo de acoplamento PT-simétrico $H_{PT}$. A dinâmica inclui superoperadores de Lindblad representando ganho em um spin e perda no outro.
• Análise de Campo Médio: Para analisar o sistema no limite termodinâmico ($S \to \infty$), os autores derivam um conjunto de seis equações diferenciais não lineares acopladas para os componentes de spin normalizados ($X, Y, Z$). Estas equações são exatas no limite termodinâmico para esta classe de modelos.
• Investigação Numérica: Os autores resolvem numericamente as equações de campo médio para mapear o diagrama de fases do estado estacionário, identificando pontos fixos, ciclos limite e regiões caóticas. A estabilidade é determinada via os autovalores da matriz Jacobiana.
• Caracterização Dinâmica: Para distinguir entre dinâmica regular e caótica, os autores computam o expoente de Lyapunov máximo médio sobre condições iniciais simétricas.
• Simulação Quântica Exata: Para validar as previsões de campo médio, os autores realizam a diagonalização numérica exata da equação mestre para tamanhos de spin finitos ($S$), trabalhando em um espaço de Liouville de tamanho $(2S+1)^4$. Eles também computam funções Wigner de spin para comparar as representações de espaço de fase do estado quântico com as trajetórias de campo médio.

Contribuições Principais e Resultados

1. Diagrama de Fases Complexo: A competição entre a transição de PT (impulsionada pelo acoplamento inter-spin $\lambda$) e a transição ferromagnética LMG (impulsionada pela não linearidade $J$) gera um diagrama de fases com múltiplas fases de estado estacionário não presentes nos modelos constituintes.

   • Fase Normal: Em baixos valores de $\lambda$ e $|J|$, o sistema exibe um único ponto fixo estável correspondendo a uma magnetização estagiada (dominada pela dissipação).
   • Fases LMG: À medida que $|J|$ aumenta, o sistema entra em regiões com dois pontos fixos estáveis, correspondendo a estados ferromagnéticos de quebra de simetria (ordenação spin-A ou spin-B).
   • Regime PT/Caótico: Em grandes valores de $\lambda$, o sistema perde pontos fixos estáveis. A teoria de campo médio prevê a emergência de ciclos limite e, em regiões de parâmetros específicas, dinâmica caótica caracterizada por expoentes de Lyapunov máximos positivos.

2. Regimes Dinâmicos Emergentes: O artigo identifica que a interação das duas transições cria comportamentos dinâmicos ausentes nos modelos individuais. Notavelmente, existem regiões onde a análise de campo médio prevê pontos fixos estáveis, mas o expoente de Lyapunov é positivo, indicando dinâmica oscilatória complexa dependente das condições iniciais (bacias de atração).

3. Validação via Escalonamento de Tamanho Finito: Os autores demonstram que as características previstas pela teoria de campo médio sobrevivem no modelo quântico completo.

   • Resultados de diagonalização exata para $S$ finito coincidem qualitativamente com as fronteiras de fase de campo médio.
   • Efeitos de tamanho finito suavizam transições abruptas (por exemplo, a transição PT em $\lambda = g/2$) e garantem um operador densidade de estado estacionário único, mesmo quando a teoria de campo médio prevê ciclos limite.
   • Funções Wigner de spin confirmam que as distribuições de probabilidade quântica alinham-se com os pontos fixos de campo médio (na fase LMG) ou distribuições difusas correspondentes a ciclos limite (na fase PT).

Significância e Alegações
O artigo alega demonstrar que o acoplamento de sistemas com distintos mecanismos de quebra de simetria dissipativa pode gerar fases dinâmicas emergentes, incluindo caos e ciclos limite complexos, que estão ausentes de qualquer um dos modelos constituintes. Os autores propõem este trabalho como um arcabouço geral para projetar fases em sistemas impulsionados-dissipativos através do acoplamento de modelos com fases conhecidas por meio de interações que suportam distintas transições não de equilíbrio.

Os autores afirmam que seu modelo é acessível utilizando tecnologias experimentais atuais, citando realizações existentes de modelos do tipo LMG em qudits supercondutores, Dysprosium ultracold, Cavidade QED e Circuit QED, combinadas com desenvolvimentos em dissipação engenheirada. Eles sugerem que tais técnicas poderiam ser relevantes para tecnologias quânticas que requerem estados estacionários ou pontos críticos engenheirados, como computação quântica dissipativa, sensoriamento quântico e termodinâmica quântica. O trabalho destaca o potencial para descobrir novas classes de universalidade e comportamentos de escalonamento de tamanho finito, embora os autores notem que determinar essas propriedades exigiria abordagens numéricas mais especializadas.