A Bogoliubov-ratio framework for quantum-information diagnostics of time-dependent two-mode Boson Hamiltonian

Resumen Técnico: Un Marco de la Relación de Bogoliubov para el Diagnóstico de Información Cuántica

Planteamiento del Problema
La dinámica de los sistemas bosónicos de dos modos dependientes del tiempo, como aquellos que describen las perturbaciones cosmológicas primordiales o los amplificadores paramétricos ópticos no degenerados, se caracterizan convencionalmente mediante una amplitud de compresión (squeezing) $r_k(\eta)$ y una fase de compresión $\phi_k(\eta)$, o equivalentemente, mediante coeficientes de Bogoliubov $\alpha_k(\eta)$ y $\beta_k(\eta)$. Si bien estas parametrizaciones son completas, su implementación directa en sistemas generales dependientes del tiempo suele ser engorrosa. Las variables $r_k$ y $\phi_k$ satisfacen ecuaciones diferenciales no lineales acopladas, y la fase $\phi_k$ puede experimentar oscilaciones rápidas. Además, la evaluación de magnitudes de la teoría de la información (como la pureza y la entropía) requiere típicamente reconstruir y diagonalizar la matriz de densidad reducida para cada modelo o protocolo de excitación específico. Este enfoque de modelo por modelo oscurece la estructura de información teórica universal inherente a los estados de par de dos modos normalizados y complica las comparaciones entre diferentes realizaciones físicas (por ejemplo, cosmológica frente a óptica).

Metodología
Los autores proponen un marco compacto y unificado basado en la relación de Bogoliubov definida como $\lambda_k(\eta) \equiv \beta_k(\eta)/\alpha_k(\eta)$, sujeta a la condición de normalización $|\alpha_k|^2 - |\beta_k|^2 = 1$ y al dominio de convergencia $|\lambda_k| < 1$.

1. Representación del Estado: Los autores establecen que un estado comprimido de dos modos generalizado y normalizado puede representarse exactamente en la base de números pareados $|n, n\rangle$ utilizando una única variable compleja $\lambda_k$:
   $$|\psi_k(\eta)\rangle = \sqrt{1 - |\lambda_k(\eta)|^2} \sum_{n=0}^{\infty} \lambda_k^n(\eta) |n_k, n_{-k}\rangle.$$
   Esta formulación absorbe las deformaciones a nivel de estado (derivadas de coeficientes de Lanczos generalizados o dinámicas de fondo específicas) directamente en $\lambda_k$.

2. Reducción Dinámica: Para un Hamiltoniano cuadrático de dos modos dependiente del tiempo de forma general:
   $$\hat{H}_k(\eta) = \omega_k(\eta) (\hat{a}^\dagger_k \hat{a}_k + \hat{a}^\dagger_{-k} \hat{a}_{-k} + 1) + i g_k(\eta) (\hat{a}^\dagger_k \hat{a}^\dagger_{-k} - \hat{a}_k \hat{a}_{-k}),$$
   la ecuación de Schrödinger se reduce a una única ecuación de tipo Riccati cerrada para $\lambda_k$:
   $$\frac{d\lambda_k}{d\eta} = g_k(\eta)(1 - \lambda_k^2) - 2i\omega_k(\eta)\lambda_k.$$
   Aquí, $\omega_k$ gobierna la rotación de fase libre y $g_k$ controla la producción de pares.

3. Diagnósticos Universales: Al trazar el modo compañero ($-k$), la fase de $\lambda_k$ desaparece. El espectro de la matriz de densidad reducida de un solo modo está determinado enteramente por la magnitud al cuadrado $q_k(\eta) = |\lambda_k(\eta)|^2$:
   $$p_n(\eta) = [1 - q_k(\eta)] q_k^n(\eta).$$
   Consecuentemente, las métricas estándar de la información teórica —pureza ($\mu_k$), entropía lineal ($S_{L,k}$), entropía de Rényi-2 ($S_{2,k}$) y entropía de von Neumann ($S_{vN,k}$)— se derivan como funciones explícitas e independientes del modelo de $q_k$.

Contribuciones Clave

• Reducción de Variables: El marco reduce la descripción de la dinámica de compresión de dos modos de dos variables reales acopladas ($r_k, \phi_k$) a una única variable compleja ($\lambda_k$) gobernada por una ecuación diferencial de primer orden.
• Fórmulas Independientes del Modelo: Los autores proporcionan fórmulas explícitas para la pureza y las entropías del estado reducido que dependen únicamente de $|\lambda_k|^2$. Esto elimina la necesidad de reconstruir o diagonalizar la matriz de densidad reducida para cada nuevo Hamiltoniano dependiente del tiempo.
• Separación de la Dinámica y el Diagnóstico: El marco factoriza claramente el problema en dos pasos:
  1. Dinámica dependiente del modelo: Resolver la ecuación de Riccati para $\lambda_k$ utilizando entradas específicas de $\omega_k(\eta)$ y $g_k(\eta)$.
  2. Diagnósticos universales: Convertir $|\lambda_k|^2$ en medidas de información utilizando fórmulas fijas y universales.

Resultados
El marco se demuestra en dos configuraciones no estacionarias distintas:

1. Perturbaciones Cosmológicas Primordiales:

   • Los autores aplican el marco a modos con frecuencia efectiva $\omega_k = \sqrt{k^2 + \alpha/\eta^2}$ y acoplamiento $g_k = -1/\eta$.
   • Las soluciones numéricas muestran que la rotación de fase inducida por el fondo y el ablandamiento de la frecuencia (controlado por el parámetro $\nu$) regulan el crecimiento de la compresión.
   • A medida que la frecuencia efectiva se ablanda (mayor $\nu$), la competencia entre la rotación de fase y la producción de pares cambia, lo que conduce a un $|\lambda_k|$ mayor, una pureza de un solo modo reducida y mayores entropías.
   • Los resultados aclaran cómo la dinámica de fondo suprime la acumulación de compresión coherente en comparación con un benchmark resonante ideal.

2. Amplificador Paramétrico No Degenerado de Pulso con Chirp Óptico:

   • El marco se aplica a un sistema con un envolvente de bomba gaussiana y un chirp de frecuencia dependiente del tiempo (desintonización).
   • Los resultados capturan el inicio retrasado de la producción de correlación (localizado cerca del centro del pulso) y la supresión de la acumulación de compresión debido al chirp.
   • A diferencia del límite de acoplamiento constante donde las entropías crecen indefinidamente, la duración finita de la bomba y el chirp de frecuencia conducen a valores de saturación finitos en la pureza y la entropía en tiempos tardíos.
   • La recuperación numérica del límite estándar $\tanh(\tau)$ para el acoplamiento resonante constante valida las convenciones de signo y la normalización de la propuesta ecuación de Riccati.

Significado y Reivindicaciones
El artículo sostiene que esta formulación de la relación de Bogoliubov ofrece una herramienta de diagnóstico eficiente y estandarizada para una amplia clase de sistemas bosónicos cuadráticos impulsados paramétricamente. Su principal significado radica en:

• Unificación: Aísla el contenido de información universal de la producción de pares de dos modos de los mecanismos físicos específicos (expansión cosmológica frente a bombeo óptico) responsables de la evolución.
• Eficiencia Computacional: Al evitar la reconstrucción de matrices de densidad reducidas para cada modelo, agiliza el cálculo de las métricas de información cuántica.
• Perspectiva Física: Proporciona un mecanismo claro para entender cómo la rotación de fase inducida por el fondo y el ablandamiento de la frecuencia regulan la mixtura del estado y el crecimiento de la compresión.

Los autores señalan que el marco es aplicable a cualquier sistema donde el estado evolucionado mantenga una representación de Bogoliubov de par normalizado generada por un Hamiltoniano cuadrático hermítico. Expresan explícitamente que las extensiones a la dinámica de sistemas genuinamente abiertos (por ejemplo, evolución de Lindblad) o interacciones no gaussianas débiles están fuera del alcance del presente trabajo.