Method to study $CP$ violation in $B_s^0\to K_S K^\pm \pi^\mp$ decays

Este artículo propone y demuestra la viabilidad de un nuevo método de análisis de diagramas de Dalitz dependiente del tiempo de decaimiento con marcado simultáneo para decaimientos $B_s^0\to K_S K^\pm \pi^\mp$ para medir con precisión la fase débil $\phi_s^{\rm eff}$ y buscar nuevas fuentes de violación de $CP$, una técnica implementada en el paquete Laura++ y validada mediante pseudoexperimentos.

Lo esencial

Imagina que el universo es una pista de baile gigante y compleja donde las partículas diminutas llamadas mesones B son los bailarines. Los físicos quieren entender las reglas de este baile para ver si coinciden con el "Modelo Estándar" (el libro de reglas actual de la física) o si existen nuevos movimientos secretos (Nueva Física) que aún no hemos descubierto.

Este artículo trata sobre un movimiento de baile específico y complicado: la desintegración del mesón $B^0_s$ (descomponiéndose en tres otras partículas: $K^0_S$, $K$ y $\pi$).

Aquí tienes el desglose de lo que hicieron los autores, utilizando analogías sencillas:

1. El misterio: Bailarines "diestros" vs. "zurdos"

En el mundo de las partículas, existe un concepto llamado violación de CP. Piensa en esto como si fuera así: Si filmas la desintegración de una partícula y luego reproduces la película hacia atrás (o miras su imagen especular), ¿se ve exactamente igual?

• Modelo Estándar: Por lo general, sí. El baile se ve igual hacia adelante que hacia atrás.
• El objetivo: Los autores buscan un baile donde la versión hacia adelante se vea diferente de la versión hacia atrás. Encontrar esta diferencia es la "prueba irrefutable" de una nueva física.

2. El desafío: Una rutina de baile de dos etapas

Normalmente, los físicos estudian un tipo de baile a la vez. Pero esta desintegración específica es especial porque puede ocurrir de dos maneras diferentes (dos estados finales) que son imágenes especulares entre sí.

• El problema: Para capturar la "violación de CP", no puedes limitarte a observar a un solo bailarín. Tienes que observar a ambos bailarines simultáneamente y comparar sus movimientos fotograma a fotograma a lo largo del tiempo.
• La analogía: Imagina intentar detectar una diferencia sutil entre un tenista diestro y uno zurdo. Si solo observas al zurdo, no puedes saber si está haciendo algo extraño. Tienes que observar a ambos jugadores en la cancha al mismo tiempo, comparando cada golpe, cada paso y cuánto tiempo permanecen en la cancha.

3. El método: El mapa de "Dalitz Plot"

Los autores proponen una nueva forma de analizar estos datos llamada análisis de Dalitz-plot.

• El mapa: Imagina un mapa de una ciudad donde cada punto representa una forma diferente en la que las partículas podrían salir disparadas.
• El factor tiempo: Esto no es solo un mapa estático; es una película. Los autores están creando un método para observar cómo los bailarines se mueven a través de este mapa a lo largo del tiempo.
• La etiqueta: Para que esto funcione, necesitan saber qué bailarín comenzó como la versión "diestra" y cuál comenzó como la versión "zurda". Esto se llama "etiquetado de sabor" (flavor tagging). Es como ponerle un sombrero rojo a un bailarín y un sombrero azul al otro al inicio de la carrera.

4. El experimento: La pista de baile "falsa"

Dado que aún no tienen suficientes datos reales del experimento LHCb para realizar este análisis tan complejo, construyeron una simulación (llamada "pseudoexperimentos").

• La simulación: Crearon un programa de computadora que generó 500 conjuntos de datos "falsos", pretendiendo ser los datos reales que el LHCb recolectará en el futuro (específicamente de los Runs 1, 2 y 3).
• La prueba: Alimentaron estos conjuntos de datos falsos con su nuevo método de análisis para ver si el método podía encontrar con éxito las señales ocultas de "violación de CP" que habían plantado en el código.

5. Los resultados: ¡Funciona!

El artículo afirma que su nuevo método es factible.

• Éxito: Cuando ejecutaron sus experimentos "falsos", el método recuperó con éxito los parámetros ocultos. Pudo distinguir entre los dos estilos de baile y medir la "diferencia de fase débil" (el ángulo de la violación de CP) con buena precisión.
• Precisión: Encontraron que con los datos que el LHCb está recolectando actualmente (Runs 1–3), pueden medir este ángulo con mucha exactitud. Si esperan a tener incluso más datos en el futuro (Runs 4–6), la precisión mejorará aún más.
• La herramienta: Ya han integrado este método en un paquete de software llamado Laura++, que otros científicos pueden utilizar.

6. Por qué es importante

• Nueva Física: Si los datos reales (cuando lleguen) muestran un resultado diferente al que predice el Modelo Estándar, significa que hay una nueva física escondida en las sombras.
• Un plano: Este artículo no solo estudia una desintegración; proporciona un plano (una receta) sobre cómo estudiar cualquier desintegración de partículas compleja que tenga dos resultados de imagen especular.

Resumen

Piensa en este artículo como un manual de entrenamiento para un juego de detectives muy difícil. Los autores inventaron una nueva forma de comparar dos bailes de partículas que son imágenes especulares a lo largo del tiempo. Probaron su método en una simulación por computadora y demostraron que funciona. Ahora, están listos para aplicar este método a los datos reales que provienen del Gran Colisionador de Hadrones para ver si al universo le quedan movimientos secretos que rompan las reglas por descubrir.

Resumen técnico

Resumen Técnico: Método para estudiar la violación de CP en decaimientos $B^0_s \to K^0_S K^\pm \pi^\mp$

Problema y Motivación
El estudio de la violación de CP en decaimientos de mesones $B$ sin carga es esencial para probar el Modelo Estándar (SM) y sondear la nueva física (NP). El canal de decaimiento $B^0_s \to K^0_S K^\pm \pi^\mp$ es de particular interés. En el SM, cuando interviene un resonancia $K^{*0}$ intermedia, el decaimiento está dominado por un único diagrama de penguin de bucle $t$, lo que implica la ausencia de violación de CP directa y de violación de CP inducida por mezcla (ya que la diferencia de fase débil se anula). Sin embargo, las contribuciones de diagramas de bucle $c$ y $u$ de orden superior o de NP podrían modificar esta expectativa. Por el contrario, cuando intervienen resonancias $K^{*\pm}$, contribuyen diagramas de nivel árbol, permitiendo potencialmente la violación de CP y la sensibilidad al ángulo CKM $\gamma$.

Para acceder al conjunto completo de observables de violación de CP en estos decaimientos, se requiere un análisis de Dalitz-plot dependiente del tiempo de decaimiento y con marcado de sabor (flavor-tagged). Aunque se han realizado tales análisis para estados finales auto-conjugados únicos (por ejemplo, $B^0 \to K^0_S \pi^+ \pi^-$) y para decaimientos de cuatro cuerpos, nunca se ha realizado un análisis simultáneo de los dos estados finales de carga conjugada ($K^0_S K^+ \pi^-$ y $K^0_S K^- \pi^+$) con dependencia del tiempo de decaimiento. Los datos actuales de LHCb de los Runs 1 y 2 carecen de la estadística para un análisis con marcado de sabor, pero se espera que el conjunto de datos del Run 3 (y ejecuciones futuras) proporcione el rendimiento suficiente.

Metodología
Los autores proponen un formalismo y un método para realizar un análisis de Dalitz-plot simultáneo, con marcado de sabor y dependiente del tiempo de decaimiento para los dos estados finales $f \equiv K^0_S K^+ \pi^-$ y $\bar{f} \equiv K^0_S K^- \pi^+$.

1. Formalismo: Las anchuras de decaimiento diferencial para los mesones $B^0_s$ y $\bar{B}^0_s$ producidos inicialmente que decaen a $f$ y $\bar{f}$ se expresan como funciones del tiempo de decaimiento $t$ y de las variables de fase-espacio $\Phi_3$. Estas ecuaciones incorporan los parámetros de mezcla $B^0_s$-$\bar{B}^0_s$ ($\Delta m_s$, $\Delta \Gamma_s$, $q/p$) y las amplitudes de decaimiento ($A_f, \bar{A}_f, A_{\bar{f}}, \bar{A}_{\bar{f}}$). La violación de CP se manifiesta como diferencias en las magnitudes de las amplitudes (CPV directa) y en la interferencia entre la mezcla y el decaimiento (CPV inducida por mezcla), caracterizada por una diferencia de fase débil efectiva $\phi^{\text{eff}}_s$.
2. Modelo de Amplitud: Las amplitudes de decaimiento se construyen utilizando el formalismo isobárico, sumando las contribuciones de resonancias intermedias. Para el estudio de viabilidad, se adopta un modelo simplificado que incluye únicamente las resonancias $K^*(892)$ ($K^{*0}$, $K^{*+}$, $K^{*-}$). Los coeficientes complejos ($c_i$) se descomponen en partes que conservan la CP ($c$) y partes que violan la CP ($\Delta c$).
3. Configuración de Pseudoexperimentos: El método se implementa en el paquete de análisis de Dalitz-plot Laura++. Los pseudoexperimentos se generan basados en el conjunto de datos LHCb Run 1–3 (aproximadamente 130,000 de rendimiento total, ~6,500 de rendimiento con marcado efectivo).
4. Estrategia de Ajuste: Se realiza un ajuste de máxima verosimilitud no parametrizado (unbinned maximum-likelihood fit) sobre conjuntos de 500 pseudoexperimentos por escenario. Para resolver las ambigüedades del ajuste, se aplican restricciones específicas:
   • El coeficiente de amplitud de referencia $c_-$ se fija en la unidad.
   • La parte imaginaria de $\Delta c_-$ se establece en cero.
   • Se aplica una restricción a la fase relativa entre $A_{\bar{f}}$ y $\bar{A}_{\bar{f}}$.
   • Las soluciones para $\phi^{\text{eff}}_s$ se restringen dentro de $\pm \pi/3$ del valor de entrada para resolver la ambigüedad de $\pi$.
   • Se realizan 100 ajustes con puntos de partida aleatorios por pseudoexperimento para asegurar la convergencia al mínimo global.

Contribuciones Clave

• Marco de Análisis Novedoso: Este artículo describe, por primera vez, el método para realizar un análisis de Dalitz-plot simultáneo, con marcado de sabor y dependiente del tiempo de decaimiento sobre dos estados finales de carga conjugada.
• Implementación: El método se implementa en el paquete Laura++, estableciendo un marco de referencia para esta clase de mediciones.
• Estudio de Sensibilidad: Los autores evalúan la sensibilidad a la fase débil efectiva $\phi^{\text{eff}}_s$ y a las fases de violación de CP en los coeficientes de amplitud a través de varias configuraciones de modelo (variando magnitudes, fases e introduciendo violación de CP en las amplitudes de $K^*$).

Resultados

• Viabilidad: Los pseudoexperimentos demuestran que los parámetros de entrada, incluyendo $\phi^{\text{eff}}_s$ y los coeficientes de amplitud, pueden determinarse de manera fiable.
• Precisión en $\phi^{\text{eff}}_s$: Para el conjunto de datos LHCb Run 1–3, se puede alcanzar una precisión estadística de aproximadamente 26 mrad en $\phi^{\text{eff}}_s$. Esto es comparable a las mediciones actuales de LHCb/CMS usando $B^0_s \to J/\psi \phi$ (22 mrad) y significativamente mejor que $B^0_s \to \phi \phi$ (75 mrad).
• Precisión en las Fases de la Amplitud: La sensibilidad a las fases de los componentes de la amplitud de violación de CP ($\arg(\Delta c_0)$, $\arg(\Delta \bar{c}_0)$) es menor, con incertidumbres alrededor de 200 mrad. Esto se atribuye al limitado solapamiento de fase-espacio entre las amplitudes de referencia de $K^*$ cargadas y las amplitudes de $K^*$ neutras.
• Proyecciones Futuras: Escalar las estadísticas a los conjuntos de datos esperados de LHCb Run 4 (60 fb$^{-1}$) y Run 6 (300 fb$^{-1}$), la incertidumbre proyectada en $\phi^{\text{eff}}_s$ mejora a 17 mrad y 8.6 mrad, respectivamente, siguiendo el escalado de raíz cuadrada inversa con la luminosidad.
• Dependencia del Modelo: Se muestra que la precisión en $\phi^{\text{eff}}_s$ depende del modelo de amplitud subyacente. En el límite de casi dos cuerpos, la sensibilidad disminuye cuando las diferencias de fase fuerte se aproximan a $\pm \pi/2$. Sin embargo, los autores señalan que un modelo más realista que incluya resonancias adicionales (por ejemplo, $K^*_0(1430)$) probablemente reduciría esta dependencia.

Significancia y Reivindicaciones
El artículo afirma que el método propuesto es viable y proporciona un camino para medir la violación de CP en los decaimientos $B^0_s \to K^0_S K^\pm \pi^\mp$ con el conjunto de datos de LHCb Run 3. Este estudio establece una base necesaria para futuros análisis de datos, que se espera que produzcan mejoras adicionales en la precisión con el Run 4 y posteriores. Los autores enfatizan que, si bien este estudio actual utiliza un modelo simplificado y descuida efectos experimentales (fondo, resolución, imperfecciones de marcado), la sensibilidad demostrada sugiere que se puede alcanzar un nivel de precisión interesante. El método también se señala como extensible a otros decaimientos multibody con múltiples estados finales. Se destaca el decaimiento contraparte $B^0 \to K^0_S K^\pm \pi^\mp$ como un objetivo futuro donde la simetría U-spin podría relacionar las amplitudes entre los sistemas $B^0_s$ y $B^0$ para controlar las incertidumbres hadrónicas.