The Tantalizing Case of the Quasar J0950+5128 -- I. Presentation of the Data and Detailed Exploration of the Binary Supermassive Black Hole Scenario

Resumen Técnico: El tentador caso del cuásar J0950+5125 – I. Presentación de los datos y exploración detallada del escenario de binaria de agujeros negros supermasivos

Problema y motivación
Se predice que las binarias de agujeros negros supermasivos (SBHB, por sus siglas en inglés) son una consecuencia natural de las fusiones de galaxias, pero la evidencia observacional inequívoca de SBHB ligadas a separaciones de sub-pársec sigue siendo esquiva. Aunque se han propuesto métodos indirectos, como la variabilidad fotométrica periódica y los cambios de velocidad radial (RV) espectroscópicos de las líneas de emisión anchas, estos métodos sufren de degeneraciones. La periodicidad fotométrica puede imitar la variabilidad estocástica de los cuásares, y los desplazamientos de las líneas espectroscópicas pueden derivarse de fenómenos no binarios, como flujos de salida (outflows), oscurecimiento por polvo o perturbaciones en la región de líneas anchas (BLR) de un solo agujero negro. El cuásar J0950+5128 (J0950), identificado en 2012 como un candidato a SBHB debido a un gran desplazamiento de velocidad en su línea ancha de H$\beta$, requiere un escrutinio adicional. Este artículo presenta 22 años de monitoreo espectroscópico para caracterizar el objeto y probar rigurosamente la hipótesis de SBHB frente a otros escenarios físicos alternativos.

Metodología
Los autores compilaron un conjunto de datos espectroscópicos que abarca 22 años (2002–2024) provenientes de los telescopios SDSS, Palomar Hale, Keck y Hobby-Eberly. El análisis consistió en:

1. Descomposición Espectral: Utilizando el código PyQSOfit, los autores modelaron y restaron la luz estelar de la galaxia anfitriona, el continuo del cuásar, la emisión de Fe II y las líneas estrechas para aislar el componente ancho de H$\beta$.
2. Medición del Desplazamiento de Velocidad: Se emplearon dos enfoques principales para medir los desplazamientos de velocidad relativos: (a) correlación cruzada de los espectros observados completos y (b) correlación cruzada de los perfiles aislados de H$\beta$ derivados de la descomposición espectral. Tanto los datos como los modelos paramétricos fueron analizados para cuantificar los errores sistemáticos, particularmente aquellos introducidos por el complejo de Fe II.
3. Simulación y Validación: Para evaluar el impacto del ensanchamiento observado del perfil de la línea H$\beta$ a lo largo del tiempo, los autores simularon series espectrales utilizando núcleos gaussianos y modelos de "respiración de la BLR". Estas simulaciones probaron si la técnica de correlación cruzada podía recuperar de manera confiable los desplazamientos de velocidad a pesar de la evolución del perfil.
4. Análisis de Curvas de Luz: Se analizaron curvas de luz ópticas (CRTS, ZTF) e infrarrojas (NEOWISE) mediante periodogramas de Lomb-Scargle y análisis de residuos para buscar variabilidad periódica consistente con el movimiento orbital binario o el realce Doppler (Doppler boosting).
5. Modelado Orbital: La curva de velocidad radial se ajustó con un modelo de órbita kepleriana utilizando el paquete radvel, incorporando un término de "jitter" (200 km s$^{-1}$) para dar cuenta de la variabilidad intrínseca del cuásar.

Resultos Clave

• Desplazamiento de Velocidad Monotónico: La línea ancha de H$\beta$ exhibe un desplazamiento de velocidad sistemático y monotónico de aproximadamente 3500 km s$^{-1}$ durante la línea de base de 22 años, cambiando de signo de negativo (desplazamiento al azul) a positivo (desplazamiento al rojo).
• Ensanchamiento del Perfil: El ancho de la mitad del máximo (FWHM) de la línea ancha de H$\beta$ ha aumentado por un factor de $\sim$3 durante el periodo de monitoreo, una tendencia que es inconsistente con los modelos simples de "respiración de la BLR", los cuales predicen un declive de flujo más fuerte de lo observado.
• Fiabilidad de la Medición: Las simulaciones confirmaron que la correlación cruzada puede recuperar de manera confiable los desplazamientos de velocidad incluso cuando los perfiles se ensanchan. Los autores determinaron que la correlación cruzada de los modelos aislados de H$\beta$ proporciona los desplazamientos más fiables, minimizando los errores sistemáticos de la contaminación de Fe II.
• Ausencia de Periodicidad Fotométrica: El análisis de la curva de luz no reveló señales periódicas significativas. La variabilidad es consistente con el comportamiento estocástico de los cuásares, y no se encontró evidencia de realce Doppler (las variaciones de magnitud predichas de $\sim$0.07 mag están por debajo de los umbrales de detección).
• Ajuste Orbital: El ajuste de la curva de velocidad radial a un modelo kepleriano arroja un periodo orbital de $33^{+7}_{-5}$ años (en el marco observado) y una excentricidad de $0.65^{+0.13}_{-0.13}$. El ajuste implica un límite inferior de la masa total de los agujeros negros de $\sim 10^7 M_\odot$ y un semieje mayor de $\sim 10^{-2}$ pc.

Interpretación y Exclusión de Alternativas
Los autores evaluaron varias interpretaciones para los desplazamientos observados:

• Oscurecimiento por Polvo: Se descartó porque la escala de tiempo para que una nube de polvo cruce la BLR requeriría una masa de agujero negro ($\gtrsim 10^{10} M_\odot$) que es físicamente implausible, y la evolución del perfil de la línea observado no coincide con las firmas de oscurecimiento esperadas.
• Flujos de Salida (Outflows): Se descartaron ya que no pueden explicar el perfil de línea simétrico observado en 2015 ni el cambio de velocidad suave y monotónico.
• Agujero Negro con Retroceso (Recoiling Black Hole): Se descartó debido a la discrepancia entre la curva de velocidad observada (que muestra aceleración constante y cambio de signo) y la dinámica esperada de un agujero negro expulsado por retroceso, así como la improbabilidad de observar tales velocidades de retroceso tan altas.
• BLR Perturbada: Aunque una BLR única y perturbada, de tipo disco, sigue siendo una explicación física viable, los autores señalan que la tendencia monotónica específica y el cambio de signo en J0950 son distintos de los emisores de doble pico conocidos y requieren un modelado adicional para confirmarse.

Significancia y Conclusiones
El artículo concluye que J0950 es un candidato convincente a SBHB. La interpretación binaria es el único escenario para el cual un modelo físico (órbita kepleriana) puede compararse inmediatamente con los datos, arrojando un ajuste consistente. Sin embargo, los autores enfatizan que los datos disponibles cubren solo un segmento de la supuesta órbita, y el escenario de una BLR única y perturbada sigue siendo viable a la espera de un modelado adicional.

La significancia de este trabajo radica en la caracterización rigurosa de la variabilidad espectroscópica de J0950 sobre una línea de base de varias décadas, la demostración de que el ensanchamiento del perfil no invalida las mediciones de desplazamiento de velocidad, y el establecimiento de J0950 como un objetivo para el monitoreo continuo. Si J0950 es efectivamente una binaria, sus parámetros orbitales actuales sugieren que podría ser detectable por las Matrices de Temporización de Púlsares (PTAs) en la banda de ondas gravitacionales de nanohertz, con una señal de fusión potencialmente detectable por futuros interferómetros espaciales como LISA. Los autores subrayan que el monitoreo espectroscópico continuo es esencial para observar el giro previsto en la curva de velocidad radial y para distinguir definitivamente entre el escenario de binaria y el de la BLR perturbada.