El número bariónico del protón no lo portarían sus quarks — la hipótesis del nudo de gluones recupera ventaja

Lo que distingue la materia de la antimateria tiene un nombre: el número bariónico. Se enseña como una propiedad de los quarks — un tercio cada uno, tres quarks, un protón. La colaboración STAR, en el colisionador RHIC de Brookhaven, publica una medición que contraría esa enseñanza: el número bariónico viaja, en las colisiones de núcleos, mucho mejor que los quarks que se supone lo portan. La explicación que mejor encaja es una vieja hipótesis que quedó en los márgenes — una topología de gluones en forma de Y, la «unión bariónica».

Fuente: science.org

En claro

Un protón son tres quarks sostenidos por pegamento: los gluones. Desde siempre se supone que la «identidad de barión» del protón — aquello que hace que cuente como una unidad de materia — está inscrita en los tres quarks, a razón de un tercio cada uno. El experimento cuenta otra historia.

Imagine dos camiones que chocan de frente. Tras el impacto, usted busca en el punto de colisión lo que pertenecía a los camiones. Los quarks son los pasajeros pesados: lanzados a toda velocidad, atraviesan y siguen su camino, muy difíciles de detener en medio. Ahora bien, en medio se encuentra mucha más «materia» de la que esos pasajeros pesados podrían depositar. Como si la identidad del camión estuviera en realidad portada por algo ligero y fácil de arrancar — un nudo en el pegamento, no los pasajeros.

Es lo que sugiere la medición. Atención: se trata de un haz de indicios convergentes que vuelve improbable la imagen de los quarks, no de una detección directa de la unión. El generador que simula la evolución completa del sistema isóbaro, UrQMD, no implementa ninguna unión y se mantiene muy por debajo del dato — lo que indica de dónde viene el problema, sin proporcionar por ello el modelo que daría cuenta de él.

Explicación técnica

1. Qué es el número bariónico, y por qué nadie sabe dónde está guardado. El número bariónico $B$B es una carga conservada: un protón vale $B = 1$B=1, un antiprotón $B = -1$B=−1. Su conservación es lo que, hasta donde sabemos, impide que el protón se desintegre y, por tanto, que la materia ordinaria se evapore. La atribución habitual — $B = 1/3$B=1/3 por quark de valencia — es una convención contable que funciona, no una medición. La cromodinámica cuántica no impone en ningún sitio que esa carga global esté localizada sobre los quarks. La alternativa propuesta en los años 1970 desplaza la carga hacia el «pegamento»: en un barión, los tres tubos de flujo de color que salen de los tres quarks deben reunirse en algún lugar para que el conjunto sea blanco de color. Ese punto de encuentro en Y es una configuración gluónica no perturbativa, eléctricamente neutra ($Q = 0$Q=0) y portadora por sí sola de $B = 1$B=1.

2. La palanca experimental: lo que se deja detener en el medio. Una colisión de núcleos ultrarrelativistas define dos extremos en rapidez, los de los haces ($y = \pm Y_{\text{beam}}$y=±Ybeam​), y un medio ($y \approx 0$y≈0). Toda la cuestión consiste en saber cuánto número bariónico logra ser frenado hasta el medio. Aquí, las dos hipótesis no dicen lo mismo, y es ahí donde se juega el mecanismo real: los quarks de valencia se llevan una gran fracción del momento del nucleón, lo que los vuelve, en los términos del artículo, «difíciles de transportar desde $Y_{\text{beam}}$Ybeam​ hasta $y \sim 0$y∼0»; la unión, en cambio, está hecha de gluones de bajo momento, luego es fácil de desprender y de depositar en el centro. Dicho de otro modo: si el número bariónico llega al medio en abundancia, es que no estaba unido a lo que va rápido.

3. La treta de los isóbaros — anular la geometría para conservar solo la carga. El rutenio 96 y el circonio 96 tienen el mismo número másico ($A = 96$A=96) pero no el mismo número de protones ($Z = 44$Z=44 frente a $Z = 40$Z=40, es decir $\Delta Z = 4$ΔZ=4). Dos núcleos del mismo tamaño: la geometría de colisión, los volúmenes, las eficiencias de detección son prácticamente idénticos y se simplifican en una diferencia. Lo que subsiste es una diferencia de carga eléctrica conocida. El razonamiento es entonces directo: si los quarks portan a la vez la carga y el número bariónico, ambos viajan juntos y se espera $\langle B\rangle/\Delta Q \approx A/\Delta Z$⟨B⟩/ΔQ≈A/ΔZ — lo que, una vez normalizado por $\Delta Z/A$ΔZ/A, da $\langle B\rangle/\Delta Q \times \Delta Z/A = 1$⟨B⟩/ΔQ×ΔZ/A=1. Los autores se cuidan de calificar ese 1 de expectativa ingenua, válida «si el número bariónico lo portan los quarks de valencia y se ignoran todos los demás efectos». Si la carga la portan los quarks pero el número bariónico una unión neutra, este último llega en exceso y el cociente supera 1. Medición: $1{,}84 \pm 0{,}02\ (\text{stat.}) \pm 0{,}09\ (\text{syst.}) \pm 0{,}16\ (\text{feed-down})$1,84±0,02 (stat.)±0,09 (syst.)±0,16 (feed-down) en las colisiones más centrales, con un decrecimiento monótono de las centrales (0–10 %) hacia las periféricas (70–80 %). Ese decrecimiento no lo atribuyen los autores al transporte bariónico sino a un efecto de estructura nuclear, que el modelo TRENTO reproduce: la piel de neutrones más gruesa del circonio que la del rutenio reduce progresivamente la fracción de protones participantes en las colisiones Zr+Zr a medida que se va hacia lo periférico. También por esta razón los puntos periféricos no deben confrontarse directamente con los generadores — los autores lo escriben explícitamente.

4. Qué miden realmente los dobles cocientes. Las cantidades brutas son
   $Q = (N_{\pi^+} + N_{K^+} + N_p) - (N_{\pi^-} + N_{K^-} + N_{\bar p})$Q=(Nπ+​+NK+​+Np​)−(Nπ−​+NK−​+Npˉ​​)
   $B = (N_p + N_n) - (N_{\bar p} + N_{\bar n})$B=(Np​+Nn​)−(Npˉ​​+Nnˉ​)
   y la diferencia de carga entre los dos sistemas se escribe como una suma de dobles cocientes,
   $\Delta Q \approx N_\pi (R2_\pi - 1) + N_K (R2_K - 1) + N_p (R2_p - 1),$ΔQ≈Nπ​(R2π​−1)+NK​(R2K​−1)+Np​(R2p​−1),
   con, por ejemplo, $R2_\pi = (N_{\pi^+}/N_{\pi^-})_{\mathrm{Ru+Ru}} / (N_{\pi^+}/N_{\pi^-})_{\mathrm{Zr+Zr}}$R2π​=(Nπ+​/Nπ−​)Ru+Ru​/(Nπ+​/Nπ−​)Zr+Zr​. El interés metodológico reside en la manera en que se tomaron los datos, y los autores lo formulan así: los haces de rutenio y de circonio se registraron «en condiciones de acelerador y de detector casi idénticas y se analizaron según el mismo procedimiento». Un cociente de cocientes construido sobre dos conjuntos tomados en esas condiciones hace desaparecer lo esencial de lo que es común a ambos. Es lo que permite una incertidumbre estadística del orden del uno por ciento sobre una cantidad tan tenue.

5. El eslabón débil, nombrado por los autores: los neutrones. $B$B contiene los neutrones, que STAR no mide directamente. Su rendimiento se estima, apoyándose en el cociente deuterón/antideuterón, mediante
   $N_n^{\text{pri}} = N_{\bar p}^{\text{pri}}\sqrt{N_d/N_{\bar d}}.$Nnpri​=Npˉ​pri​Nd​/Ndˉ​​.
   A ello se suma la desintegración de los hiperones, estimada a partir de sus rendimientos y de razones de ramificación conocidas, que contribuye en torno al 30 % al rendimiento total de neutrones. Es exactamente por eso que el resultado lleva una tercera barra de error, separada y dominante: $\pm 0{,}16$±0,16 de «feed-down», es decir, más que las incertidumbres estadística y sistemática reunidas.

6. La prueba fotonuclear, más limpia porque el fotón no aporta nada. En los eventos $\gamma+\mathrm{Au}$γ+Au, uno de los dos núcleos solo interactúa por su campo electromagnético: el fotón no porta ningún número bariónico. Todo exceso de bariones viene por tanto de un solo lado, lo que suprime la contaminación del compañero. Se ajusta entonces la densidad de protones netos en rapidez mediante una exponencial en la distancia al haz,
   $f(y) \propto \exp(-\alpha_B \Delta y), \qquad \Delta y = Y_{\text{beam}} - y,$f(y)∝exp(−αB​Δy),Δy=Ybeam​−y,
   con $Y_{\text{beam}} = 4{,}06$Ybeam​=4,06. El parámetro $\alpha_B$αB​ indica a qué velocidad se agota el transporte del número bariónico cuando uno se aleja del haz: cuanto menor es, más lejos viaja la carga. La teoría de Regge aplicada a la unión predice $0{,}42 < \alpha_B < 1$0,42<αB​<1; la medición da $\alpha_B = 1{,}04 \pm 0{,}22$αB​=1,04±0,22. Hay que leer esa cifra por lo que es: el valor central se sitúa ligeramente por encima de la cota alta predicha, y la compatibilidad con el intervalo solo se sostiene por la incertidumbre de $\pm 0{,}22$±0,22.

7. El control que hace defendible la medición. Una pendiente siempre puede ser un artefacto de selección de eventos. Los autores miden por ello la misma pendiente para los antiprotones: $\alpha_{\bar B} = 0{,}02 \pm 0{,}05$αBˉ​=0,02±0,05, es decir, plana. Los antiprotones se crean por pares en la colisión y no transportan el número bariónico de los haces; si no muestran ninguna dependencia en rapidez mientras que los protones netos sí la muestran, es que la dependencia observada es efectivamente una propiedad del transporte bariónico y no del dispositivo. Es un control experimental, y no un argumento teórico, lo que aquí vuelve improbable el artefacto. Por último, en las colisiones Au+Au, la pendiente media vale $0{,}64 \pm 0{,}05$0,64±0,05 y no depende de la centralidad — coherente con el valor fotonuclear dentro de 1,7 σ — allí donde la imagen de los quarks haría esperar una variación con el número de dispersiones múltiples.

Por qué funcionó

El experimento funciona porque transforma una cuestión de localización — dónde está guardada una carga — en una cuestión de transporte, que sí se mide. El contraste cuantitativo es nítido: allí donde la imagen de los quarks de valencia conduce ingenuamente a $\langle B\rangle/\Delta Q \times \Delta Z/A = 1$⟨B⟩/ΔQ×ΔZ/A=1, el dato muestra 1,84, es decir, casi el doble, con una incertidumbre total del orden de 0,2 sumando las tres contribuciones. El generador que simula la evolución completa del sistema isóbaro, UrQMD, no implementa ninguna unión y predice un valor de 0,5 a 0,7 — «significativamente por debajo» del dato experimental, escriben los autores.

Conviene, sin embargo, señalar con claridad la distancia entre lo que la medición establece y lo que uno se sentiría tentado de concluir. En primer lugar, la comparación con los generadores exige una precaución de lectura. Los autores sí llevan a la misma figura que sus puntos de medición cinco cálculos de modelos — UrQMD, TRENTO, PYTHIA 8.3, PYTHIA 8.3 con reconexión de color en modo 2, y HERWIG 7.2 — de los cuales el cuarto es precisamente el que fabrica dinámicamente uniones. Pero los valores de PYTHIA y de HERWIG — alrededor de 0,5 a 0,6 para los ajustes por defecto, $0{,}99 \pm 0{,}03$0,99±0,03 para el modo con uniones — están calculados para el caso límite de colisiones protón-protón a 200 GeV, y no para el sistema isóbaro en sí. Ahora bien, los únicos puntos isóbaros cuyo número de nucleones participantes se aproxima a ese límite son los puntos periféricos, y es precisamente sobre ellos sobre los que los autores advierten: a causa del efecto de piel de neutrones, «no deberían compararse directamente con las predicciones de HERWIG y PYTHIA». No existe, por tanto, una confrontación cuantitativa limpia entre un generador con uniones y el 1,84 de las colisiones centrales: la comparación se muestra, pero no es explotable como una sustracción.

En segundo lugar, la unión tampoco escapa a la crítica sobre el otro observable: los dos ajustes de PYTHIA 8.3, incluido el que produce uniones, sobrestiman significativamente la pendiente medida en $\gamma+\mathrm{Au}$γ+Au — «quizá debido a la ausencia de ciertos ingredientes físicos en el modelo», señalan los autores. La unión gana una comparación a dos frente a la imagen de los quarks; no por ello logra un ajuste.

En tercer lugar, los autores no anuncian ninguna significancia estadística en σ para rechazar la imagen de los quarks: su conclusión es que estos resultados, respaldados por los datos Au+Au publicados anteriormente por STAR a varias energías, «desfavorecen la imagen de los quarks de valencia». Desfavorecer no es excluir, y una preferencia de modelo no es una detección. Por último, el objeto mismo escapa al cálculo: una unión es no perturbativa, luego imposible de tratar con las herramientas habituales de la QCD perturbativa, lo que deja la predicción teórica en un intervalo amplio más que en un valor. A ello se suman limitaciones de muestra reconocidas: los eventos fotonucleares seleccionados están sesgados hacia las altas multiplicidades, y los autores indican que $\alpha_B$αB​ podría depender de la cinemática de colisión sin que las mediciones actuales permitan dirimirlo.

Cadena causal

Confinamiento del color: un barión debe ser blanco → los tres tubos de flujo que salen de los quarks deben reunirse en un punto → esa topología en Y se predice como objeto propio por Artru (1975) y luego Rossi y Veneziano (1977) → la unión, hecha de gluones de bajo momento, se separa de los quarks rápidos durante una colisión → es frenada hasta la mitad de rapidez mientras los quarks de valencia siguen su camino → el número bariónico llega al centro en exceso respecto de la carga eléctrica → medido sobre isóbaros Ru/Zr, ese cociente vale 1,84 allí donde la expectativa ingenua de los quarks de valencia da 1 → medido en fotonuclear, donde el fotón no aporta ningún barión, la pendiente de transporte es compatible con la ventana prevista para la unión, superando su valor central ligeramente la cota alta → la imagen de los quarks de valencia queda desfavorecida por tres observables independientes → pero los generadores con uniones solo están calculados para el caso límite protón-protón, y sobrestiman la pendiente fotonuclear → el colisionador electrón-ión en construcción en Brookhaven heredará la cuestión.

Anécdota

Los cuatro mil millones de colisiones de rutenio y circonio — dos mil millones por especie — no se tomaron para esa física. En 2018, RHIC alternó esos dos haces — a veces varias veces el mismo día — para rastrear un fenómeno completamente distinto: el efecto magnético quiral. El razonamiento era que el rutenio, con cuatro protones más, genera un campo magnético más intenso que el circonio y debería por tanto producir una separación de cargas más marcada. Para protegerse de su propio deseo de encontrar, los físicos de STAR montaron un análisis a ciegas en tres etapas, con cinco grupos independientes trabajando sin saber de qué especie provenían los datos. El veredicto, emitido en 2021, fue negativo: ninguna señal de efecto magnético quiral conforme a las expectativas. Es ese conjunto de datos, constituido para una pregunta y que no dio la respuesta esperada, el que sirve hoy para pesar una hipótesis de medio siglo sobre la naturaleza del protón.

Legado y datos actuales

Lo que está en juego va más allá de la contabilidad interna del protón. Si el número bariónico es una propiedad topológica del campo de gluones más que una etiqueta portada por partículas, entonces los escenarios que buscan explicar por qué el universo contiene materia y casi nada de antimateria — la bariogénesis — trabajan sobre un objeto cuya naturaleza no es la que se suponía. A estas alturas, hay que ser claro sobre lo que no está establecido: esta medición no dice nada sobre la violación del número bariónico, no proporciona ninguna tasa y no tiene implicación cuantificada para la desintegración del protón. Lo adquirido es más modesto y más sólido: un mecanismo de transporte privilegiado, medido en dos sistemas de colisión independientes, más una discrepancia con los generadores sin unión que los autores aportan al expediente.

La mirada del investigador — preguntas abiertas

Esta sección corresponde a la interpretación y no a los resultados del artículo. Tres experimentos dirimirían la cuestión. El barrido en energía, en primer lugar: si $\alpha_B$αB​ es efectivamente una propiedad de la unión y no un accidente de cinemática, debe permanecer estable al cambiar $\sqrt{s_{NN}}$sNN​​ — los propios autores señalan que sus datos no permiten verificarlo. La prueba de portabilidad, después: la misma medición sobre otros pares isóbaros, con un $\Delta Z$ΔZ distinto, indicaría si 1,84 es una constante del mecanismo o depende del par elegido. La restricción teórica, por último: al ser la unión no perturbativa, solo la QCD en la red puede en principio proporcionar una predicción cuantificada de su contribución al número bariónico. Mientras ese cálculo falte — y mientras los generadores con uniones solo estén calculados para el caso límite protón-protón, fuera del régimen en que se mide el 1,84 — no se podrá excluir que un efecto de transporte ordinario, no modelizado, explique una parte de la discrepancia. Una observación de método para terminar, que es mía y no de los autores: mostrar la incertidumbre de feed-down como una tercera barra de error separada, en lugar de fundirla en la sistemática, es una elección de transparencia — hace visible que el límite dominante de esta medición es la estimación de los neutrones.

Fuentes

Referencias verificadas durante la auditoría de fact-checking del 14 de agosto de 2026 — texto íntegro de la versión depositada en arXiv consultado (Methods, Results y leyendas de figuras), metadatos de publicación verificados por separado. Las predicciones de generadores citadas aquí se resitúan en el sistema de colisión para el que los autores las calcularon.

• STAR Collaboration, «Tracking the baryon number with nuclear collisions», Science 393, n.º 6812, p. 727–731, 13 de agosto de 2026. DOI: 10.1126/science.ads5962 — artículo evaluado por pares.
• Versión íntegra depositada: arXiv:2408.15441, enviada el 27 de agosto de 2024. arxiv.org/abs/2408.15441 — de esta versión provienen los valores numéricos detallados, las ecuaciones y las limitaciones citadas aquí.

Referencias de fondo

Estas referencias sustentan los recordatorios de contexto y no provienen del estudio comentado.

• X. Artru, «Classical string phenomenology» (1975) y G. C. Rossi, G. Veneziano (1977) — propuestas de origen de la unión bariónica, citadas como tales en las referencias del estudio.
• STAR Collaboration, «Methods for a blind analysis of isobar data collected by the STAR collaboration», arXiv:1911.00596 — protocolo de análisis a ciegas del programa isóbaro, para la anécdota sobre el origen de los datos de 2018.
• Brookhaven National Laboratory, comunicado «Results from search for "Chiral Magnetic Effect" at RHIC» (2021) — resultado negativo de la búsqueda para la que se habían programado las colisiones isóbaras.