A Pythia8 Tune for Open Charm and Beauty Production in Fixed-Target Collisions

Muy por encima de la Tierra, partículas invisibles caen constantemente, nacidas cuando los rayos cósmicos golpean la atmósfera. Entre estas se encuentran los neutrinos, partículas fantasmales que rara vez interactúan con algo, pero que portan secretos sobre los eventos más violentos del universo. Para comprender de dónde vienen estos neutrinos, los científicos deben rastrear sus orígenes hasta el momento en que nacieron: la desintegración de partículas pesadas como los hadrones de encanto y belleza. Estas partículas pesadas se crean cuando haces de alta energía de protones o piones chocan contra un blanco, un proceso que ocurre no solo en la atmósfera superior, sino también en laboratorios especializados donde los investigadores disparan haces de partículas hacia bloques densos de materia. Predecir cuántas de estas partículas pesadas se fabrican, y exactamente cómo se mueven tras su creación, es esencial para calcular el flujo de neutrinos que alcanzan nuestros detectores. Si las predicciones fallan, el ruido de fondo en los experimentos que buscan nuevas partículas ocultas se vuelve imposible de distinguir de un descubrimiento genuino.

Durante décadas, los físicos han dependido de simulaciones por computadora para modelar estas colisiones, utilizando software complejo que actúa como un laboratorio virtual. Uno de esos programas, Pythia, ha sido la herramienta estándar durante décadas, pero fue sintonizado originalmente para coincidir con datos de colisiones de alta energía en el Gran Colisionador de Hadrones (LHC), donde las partículas chocan a energías que superan con creces lo que los experimentos de blanco fijo pueden lograr. Cuando los científicos intentaron usar esta configuración estándar para predecir qué sucede en experimentos de blanco fijo —donde los haces golpean blancos estacionarios a energías más bajas— los resultados eran a menudo erróneos. El software tendía a predecir que las partículas pesadas se moverían demasiado rápido y se producirían en mayor número de lo que realmente se observaba en los datos. Este desajuste creó un punto ciego para experimentos como SHiP, que están diseñados para cazar partículas raras y débilmente interactuantes mediante el análisis de los desechos de un enorme depósito de haces (beam dump). Sin un mapa más preciso de cómo nacen y se comportan las partículas pesadas en este rango de energía específico, la búsqueda de nueva física se ve obstaculizada por la incertidumbre.

Para solucionar esto, un equipo de investigadores de la Universidad de Gante ha desarrollado un nuevo conjunto de instrucciones para el software Pythia, adaptado específicamente para las condiciones que se encuentran en las colisiones de blanco fijo. Llaman a esta nueva configuración FTFT. El equipo no se limitó a adivinar los ajustes correctos; en su lugar, tomaron una vasta colección de mediciones del mundo real procedentes de experimentos realizados durante los últimos cuarenta años. Estas mediciones provenían de haces tanto de protones como de piones impactando contra diversos blancos metálicos, registrando la velocidad y dirección de las partículas de encanto y belleza resultantes. Al introducir estos datos históricos en sus herramientas de optimización, los investigadores ajustaron las perillas internas de la simulación hasta que las colisiones virtuales coincidieron con las reales. Se centraron en dos cosas principales: la forma de las distribuciones de partículas, como qué tan hacia adelante vuelan las partículas y cuánto momento lateral transportan, y el número total de partículas producidas.

Los resultados de este ajuste fueron sorprendentes. Los ajustes de software estándar, que funcionan bien para las colisiones de más alta energía, no lograban describir los datos de blanco fijo, prediciendo un espectro de velocidades de partículas demasiado agresivo. El nuevo ajuste FTFT corrigió esto alterando la forma en que la simulación maneja la ruptura de los quarks pesados en partículas detectables y cómo ocurren múltiples interacciones dentro de una sola colisión. Los investigadores descubrieron que las partículas pesadas en colisiones de blanco fijo se fragmentan de manera diferente a como lo hacen en el LHC, lo que requiere un proceso de ruptura más suave y gentil en la simulación. También tuvieron que ajustar cómo la simulación escala con la energía, moviendo el punto de referencia desde las energías masivas del LHC hacia el rango relevante para los haces de blanco fijo. Una vez que estos parámetros de forma quedaron fijados, el equipo aplicó un factor de escala final para hacer coincidir el número total de partículas producidas con las mediciones experimentales. Para los haces de protones, la simulación necesitaba ser potenciada por un factor de aproximadamente 2.5, mientras que para los haces de piones, el aumento era de unos 2.0. Estos factores dan cuenta de la compleja física que el modelo estándar omite a estas energías.

Para las partículas de belleza, aún más pesadas, la situación era ligeramente distinta. Debido a que la producción de belleza es tan rara y difícil de medir en detalle a estas energías, el equipo no pudo ajustar la forma de la distribución. En su lugar, simplemente aplicaron un factor de escala al número total de pares de belleza producidos. Sorprendentemente, para la belleza, la simulación estándar ya estaba bastante cerca de los datos, requiriendo solo un ajuste minúsculo de aproximadamente un 4% para los haces de protones y un 19% para los de piones. Esto sugiere que el proceso fundamental de creación de pares de belleza pesados está mejor comprendido que el proceso de creación de partículas de encanto en este régimen. El nuevo ajuste también abordó un enigma específico respecto a la carga de las partículas. En colisiones con haces de piones, existe una preferencia conocida por producir partículas de encanto negativas sobre las positivas, un fenómeno impulsado por la estructura interna del propio haz de piones. La simulación estándar luchaba por capturar esta asimetría, pero el nuevo ajuste FTFT, al incorporar un conjunto específico de reglas sobre cómo están estructurados los piones, acercó la predicción mucho más a la realidad, aunque permaneció una pequeña discrepancia.

Las implicaciones de este trabajo son inmediatas para la próxima generación de experimentos de física. El experimento SHiP, que se está construyendo actualmente para buscar partículas ocultas, depende fuertmente de predicciones precisas de los fondos de neutrinos. Las simulaciones antiguas subestimaban la producción de partículas de encanto, lo que significaba que las estimaciones del fondo probablemente eran demasiado bajas. Con el ajuste FTFT, las predicciones para el flujo de neutrinos y las tasas de fondo están ahora fundamentadas en un modelo que ha sido rigurosamente probado contra décadas de datos. Esto le otorga al equipo de SHiP una base mucho más firme para su búsqueda, permitiéndoles distinguir entre una fluctuación del fondo y una verdadera señal de nueva física con mayor confianza. Más allá de SHiP, este ajuste también es relevante para comprender el componente prompt de los neutrinos atmosféricos, que se producen cuando los rayos cósmicos impactan la atmósfera terrestre. Dado que estas cascadas atmosféricas involucran muchas colisiones secundarias a energías similares a las de los experimentos de blanco fijo, el nuevo ajuste ofrece una mejor manera de modelar las partículas pesadas que contribuyen al flujo de neutrinos medido por los telescopios.

Los investigadores advierten cuidadosamente que, si bien este ajuste es una mejora significativa, no es una respuesta final para todos los escenarios posibles. El ajuste se realizó específicamente para el rango de energía de 20 a 42 GeV, y los parámetros se ajustaron para adaptarse a los datos disponibles, que se refieren principalmente a los mesones de encanto. El comportamiento de los bariones de encanto, primos más pesados de los mesones, no estuvo restringido por los datos y sigue siendo un área de incertidumbre. Además, el ajuste asume que la producción de estas partículas escala linealmente con el tamaño del núcleo del blanco, un supuesto respaldado por los datos existentes, pero que podría probarse más a fondo. A pesar de estas limitaciones, el trabajo representa un paso crucial hacia la resolución de la brecha entre la física de colisionadores de alta energía y el régimen de blanco fijo de menor energía. Al fundamentar la simulación en la realidad de las mediciones pasadas, el equipo ha proporcionado una herramienta que convierte una fuente de incertidumbre en un predictor fiable, asegurando que la búsqueda de las partículas ocultas del universo esté guiada por el mapa más preciso disponible.