Even-harmonic generation through nonequilibrium steady-state symmetry breaking

Resumo Técnico: Geração de Harmônicos Pares através da Quebra de Simetria em Estado Estacionário de Não Equilíbrio

Definição do Problema
Em sistemas quânticos com simetria de inversão submetidos a luz polarizada linearmente, a geração de harmônicos de alta ordem (HHG) é estritamente limitada pela simetria para emitir apenas harmônicos ímpares. A observação de harmônicos pares é tradicionalmente interpretada como uma assinatura de quebra de simetria de inversão, tipicamente alcançada através de assimetria estrutural, campos de polarização estática ou condução assimétrica. No entanto, uma questão fundamental permanece sobre como correntes de transporte fora do equilíbrio modificam as respostas ópticas não lineares ultrarrápidas em sistemas centrosimétricos onde o Hamiltoniano subjacente retém a simetria de inversão. Embora estudos recentes sugiram que correntes finitas podem induzir distribuições eletrônicas assimétricas, o mecanismo específico pelo qual a condução dissipativa de fronteira em sistemas quânticos abertos quebra as regras de seleção para HHG, sem alterar a estrutura da rede, permaneceu amplamente inexplorado.

Metodologia
Os autores investigam este fenômeno utilizando uma cadeia Su-Schrieffer-Heeger (SSH) unidimensional acoplada a reservatórios de fonte e sumidouro, formando um sistema quântico aberto. O arcabouço teórico baseia-se na equação mestra de Lindblad para descrever a dinâmica dissipativa.

• Modelo: O sistema consiste em uma rede dimerizada com amplitudes de salto alternadas ($v$ e $w$). A cadeia é acoplada a um reservatório de ganho (fonte) no último síto e a um reservatório de perda (sumidouro) no primeiro síto com taxas $\gamma_g$ e $\gamma_l$, respectivamente.
• Dinâmica: A evolução da matriz de densidade de muitos corpos $\rho$ é reduzida à dinâmica da matriz de correlação de partícula única $C_{mn}(t) = \langle c^\dagger_m c_n \rangle$. Os autores derivam uma equação de movimento fechada para $C(t)$ baseada no formalismo de Lindblad, que inclui a evolução unitária coerente impulsionada pelo Hamiltoniano e termos dissipativos que contabilizam a troca de partículas com os reservatórios.
• Estado Estacionário: O estado estacionário de não equilíbrio (NESS) é obtido analiticamente resolvendo a equação de Sylvester (Lyapunov) resultante para a matriz de correlação: $MC_{ss} + C_{ss}N = -\Gamma^+$. Esta abordagem permite a determinação exata da matriz de densidade do estado estacionário sem integração temporal.
• Condução e Cálculo de HHG: Um campo laser espaço-temporal é aplicado ao sistema, localizado no interior (bulk) para evitar o acionamento direto dos sítos acoplados ao reservatório. O espectro de HHG é calculado a partir da transformada de Fourier do valor esperado do operador de corrente no tempo, $\langle \hat{J}(t) \rangle$.
• Validação: Para abordar o realismo experimental, os autores também simulam a evolução de Lindblad dependente do tempo partindo de um estado inicial para observar o processo de relaxação em direção ao estado estacionário, comparando a dinâmica transiente com a solução exata do estado estacionário via análise de fidelidade.

Principais Contribuições

1. Mecanismo de Quebra de Simetria: O artigo demonstra que a regra de seleção para HHG em sistemas quânticos abertos é governada pela simetria do superoperador de Liouville completo e pela matriz de densidade de estado estacionário resultante, e não apenas pela simetria do Hamiltoniano.
2. Harmônicos Pares Induzidos por Corrente: Demonstra-se que a condução dissipativa de fronteira gera uma corrente contínua (DC) finita em uma cadeia SSH com inversão de simetria. Este estado de condução de corrente quebra a simetria de inversão da matriz de densidade ($P\rho_{ss}P^{-1} \neq \rho_{ss}$), levantando assim a supressão de harmônicos pares.
3. Correlação Direta: O estudo estabelece um vínculo direto entre a magnitude da corrente de transporte e a intensidade dos harmônicos pares gerados. A emergência de harmônicos pares não se deve à distorção estrutural da rede, mas é uma consequência direta da quebra de simetria do estado estacionário de não equilíbrio.

Resultados

• Evolução Espectral: No limite de dissipação evanescente ($\gamma \to 0$), o sistema comporta-se como uma cadeia fechada, e o espectro de HHG é dominado por harmônicos ímpares, consistente com as regras de seleção convencionais. À medida que a força da dissipação $\gamma$ aumenta, harmônicos de ordem par emergem e sua intensidade cresce, modificando o espectro de baixa ordem.
• Dependência da Corrente: A corrente de estado estacionário aumenta com a força do acoplamento do reservatório. O aparecimento de harmônicos pares coincide com o regime onde o sistema suporta uma corrente de não equilíbrio significativa.
• Dinâmica de Relaxação: A análise do gap de Liouville revela que, no regime de baixa dissipação, o tempo de relaxação para o estado estacionário diverge. Consequentemente, a emergência de harmônicos pares depende criticamente do tempo de relaxação disponível ($\tau$). Se o sistema for sondado antes de relaxar totalmente para o estado estacionário, a memória do estado inicial de simetria de inversão suprime os harmônicos pares. Cálculos de fidelidade confirmam que o cruzamento para a geração de harmônicos pares é governado pela escala de tempo da relaxação de Liouville.
• Distribuição de Partículas: A densidade de partículas no estado estacionário transita de um estado do tipo vácuo para valores muito baixos de $\gamma$ para um estado de preenchimento de meia densidade em altos valores de $\gamma$, sendo que o regime de geração significativa de harmônicos pares corresponde ao estabelecimento de um estado robusto de condução de corrente.

Significância
O artigo afirma descobrir um mecanismo para a geração de harmônicos pares que depende exclusivamente da quebra de simetria do estado estacionário de não equilíbrio induzida pelo transporte dissipativo, sem requerer qualquer modificação ao Hamiltoniano subjacente ou quebra de simetria de inversão estrutural da rede. Esta descoberta redefine a compreensão das regras de seleção de harmônicos em sistemas quânticos abertos, identificando a simetria da matriz de densidade de não equilíbrio como o fator determinante.

Os autores postulam que seus resultados fornecem uma rota para sondar correntes de transporte através de espectroscopia não linear ultrarrápida em sistemas quânticos centrosimétricos. Ao estabelecer a HHG como uma sonda sensível de estados estacionários de não equilíbrio que conduzem corrente, o trabalho sugere novas oportunidades para explorar a interdependência entre transporte, dissipação e quebra de simetria em materiais quânticos abertos sob condução. O estudo mantém um escopo modesto, focando em cadeias unidimensionais não interagentes, e identifica a extensão para sistemas interagentes, dimensões superiores e descrições não markovianas (por exemplo, via funções de Green de Keldysh) como direções futuras necessárias.