Even-harmonic generation through nonequilibrium steady-state symmetry breaking

Resumen Técnico: Generación de Armónicos Pares mediante la Ruptura de Simetría en el Estado Estacionario Fuera del Equilibrio

Planteamiento del Problema
En sistemas cuánticos con simetría de inversión excitados por luz de polarización lineal, la generación de altos armónicos (HHG) está estrictamente restringida por la simetría para emitir únicamente armónicos impares. La observación de armónicos pares se interpreta tradicionalmente como una firma de la ruptura de la simetría de inversión, lograda típicamente mediante asimetría estructural, campos de sesgo estáticos o excitación asimétrica. Sin embargo, persiste una pregunta fundamental sobre cómo las corrientes de transporte fuera del equilibrio modifican las respuestas ópticas no lineales ultrarrápidas en sistemas centrosimétricos donde el Hamiltoniano subyacente conserva la simetría de inversión. Si bien estudios recientes sugieren que las corrientes finitas pueden inducir distribuciones electrónicas asimétricas, el mecanismo específico mediante el cual el accionamiento de fronteras disipativo en sistemas cuánticos abiertos rompe las reglas de selección para la HHG sin alterar la estructura de la red ha permanecido mayormente inexplorado.

Metodología
Los autores investigan este fenómeno utilizando una cadena de Su-Schrieffer-Heeger (SSH) unidimensional acoplada a reservorios de fuente y sumidero, formando un sistema cuántico abierto. El marco teórico se basa en la ecuación maestra de Lindblad para describir la dinámica disipativa.

• Modelo: El sistema consiste en una red dimerizada con amplitudes de salto alternantes ($v$ y $w$). La cadena está acoplada a un reservorio de ganancia (fuente) en el último sitio y a un reservorio de pérdida (sumidero) en el primer sitio con tasas $\gamma_g$ y $\gamma_l$, respectivamente.
• Dinámica: La evolución de la matriz de densidad de muchos cuerpos $\rho$ se reduce a la dinámica de la matriz de correlación de una sola partícula $C_{mn}(t) = \langle c^\dagger_m c_n \rangle$. Los autores derivan una ecuación de movimiento cerrada para $C(t)$ basada en el formalismo de Lindblad, que incluye la evolución unitaria coherente impulsada por el Hamiltoniano y términos disipativos que dan cuenta del intercambio de partículas con los reservorios.
• Estado Estacionario: El estado estacionario fuera del equilibrio (NESS) se obtiene analíticamente resolviendo la ecuación de Sylvester (Lyapunov) resultante para la matriz de correlación: $MC_{ss} + C_{ss}N = -\Gamma^+$. Este enfoque permite la determinación exacta de la matriz de densidad en estado estacionario sin necesidad de integración temporal.
• Accionamiento y Cálculo de HHG: Se aplica un campo láser espaciotemporal al sistema, localizado en el interior (bulk) para evitar el accionamiento directo de los sitios acoplados a los reservorios. El espectro de HHG se calcula a partir de la transformada de Fourier del valor esperado dependiente del tiempo del operador de corriente, $\langle \hat{J}(t) \rangle$.
• Validación: Para abordar el realismo experimental, los autores también simulan la evolución de Lindblad dependiente del tiempo partiendo de un estado inicial para observar el proceso de relajación hacia el estado estacionario, comparando la dinámica transitoria con la solución exacta del estado estacionario mediante un análisis de fidelidad.

Contribuciones Clave

1. Mecanismo de Ruptura de Simetría: El artículo demuestra que la regla de selección para la HHG en sistemas cuánticos abiertos está gobernada por la simetría del superoperador de Liouville completo y la matriz de densidad resultante en estado estacionario, más que únicamente por la simetría del Hamiltoniano.
2. Armónicos Pares Inducidos por Corriente: Se muestra que el accionamiento de fronteras disipativo genera una corriente continua (DC) finita en una cadena SSH con simetría de inversión. Este estado estacionario portador de corriente rompe la simetría de inversión de la matriz de densidad ($P\rho_{ss}P^{-1} \neq \rho_{ss}$), levantando así la supresión de los armónicos pares.
3. Correlación Directa: El estudio establece un vínculo directo entre la magnitud de la corriente de transporte y la intensidad de los armónicos pares generados. La aparición de armónicos pares coincide con el régimen donde el sistema soporta una corriente fuera del equilibrio significativa.

Resultados

• Evolución Espectral: En el límite de disipación evanescente ($\gamma \to 0$), el sistema se comporta como una cadena cerrada, y el espectro de HHG está dominado por armónicos impares, consistente con las reglas de selección convencionales. A medida que la fuerza de la disipación $\gamma$ aumenta, emergen armónicos de orden par y su intensidad crece, modificando el espectro de bajo orden.
• Dependencia de la Corriente: La corriente en estado estacionario aumenta con la fuerza del acoplamiento al reservorio. La aparición de armónicos pares coincide con el régimen donde el sistema sostiene un estado de corriente fuera del equilibrio significativo.
• Dinámica de Relajación: El análisis de la brecha de Liouville revela que, en el régimen de disipación débil, el tiempo de relajación hacia el estado estacionario diverge. Consecuentemente, la emergencia de armónicos pares depende críticamente del tiempo de relajación disponible ($\tau$). Si el sistema es sondeado antes de relajarse completamente al estado estacionario, la memoria de su estado inicial de simetría de inversión suprime los armónicos pares. Los cálculos de fidelidad confirman que la transición hacia la generación de armónicos pares está gobernada por la escala de tiempo de la relajación de Liouville.
• Distribución de Partículas: La densidad de partículas en estado estacionario transiciona de un estado similar al vacío a un estado de medio llenado a medida que $\gamma$ aumenta, siendo el régimen de generación significativa de armónicos pares correspondiente al establecimiento de un estado robusto portador de corriente.

Significancia
El artículo afirma haber descubierto un mecanismo para la generación de armónicos pares que depende únicamente de la ruptura de simetría del estado estacionario fuera del equilibrio inducida por el transporte disipativo, sin requerir ninguna modificación al Hamiltoniano subyacente o la ruptura de la simetría de inversión estructural de la red. Este hallazgo redefine la comprensión de las reglas de selección de armónicos en sistemas cuánticos abiertos, identificando la simetría de la matriz de densidad fuera del equilibrio como el factor determinante.

Los autores postulan que sus resultados proporcionan una vía para sondear corrientes de transporte mediante espectroscopía no lineal ultrarrápida en sistemas cuánticos centrosimétricos. Al establecer la HHG como una sonda sensible de los estados estacionarios fuera del equilibrio que portan corriente, el trabajo sugiere nuevas oportunidades para explorar la interacción entre transporte, disipación y ruptura de simetría en materiales cuánticos abiertos excitados. El estudio mantiene un alcance modesto, centrándose en cadenas unidimensionales no interactuantes, e identifica la extensión a sistemas interactuantes, dimensiones superiores y descripciones no markovianas (por ejemplo, mediante funciones de Green de Keldysh) como direcciones futuras necesarias.