Fully differential Higgs boson pair production at N$^3$LO with top quark mass effects

Este artículo presenta las primeras predicciones totalmente diferenciales para la producción de pares de bosones de Higgs mediante la fusión de gluones-gluones a N$^3$LO en el límite de top pesado, demostrando una reducción de tres veces en las incertidumbres de escala e incorporando efectos de la masa del quark top para proporcionar referencias teóricas altamente precisas para las búsquedas en el LHC.

Lo esencial

Imagina el universo como una máquina gigante y compleja, y el bosón de Higgs como un engranaje crucial que otorga peso a otras partículas. Los físicos quieren entender cómo funciona este engranaje chocando partículas a velocidades increíasbles en el Gran Colisionador de Hadrones (LHC). Específicamente, están tratando de ver qué sucede cuando se crean dos bosones de Higgs al mismo tiempo. Esto es como intentar atrapar dos mariposas raras y esquivas en medio de una tormenta para ver cómo interactúan.

Este artículo es un salto masivo hacia adelante en el "manual de instrucciones" (predicción teórica) sobre cómo detectar estos pares de mariposas. Aquí está el desglose en términos sencillos:

1. El Problema: Un bucle muy pesado

Para crear dos bosones de Higgs, las partículas chocan y crean un "bucle" temporal que involucra a un quark top (la partícula más pesada conocida).

• La Analogía: Imagina intentar predecir la trayectoria de una pelota rodando a través de un laberinto. El laberinto está hecho de quarks top. Debido a que el quark top es tan pesado, el laberinto es increíblemente complejo.
• La Forma Antigua: Durante años, los científicos usaron un atajo llamado "Límite del Top Pesado". Pretendían que el quark top era infinitamente pesado, lo que suavizaba el laberinto convirtiéndolo en un suelo liso y simple. Esto hacía que las matemáticas fueran más fáciles, pero no era perfectamente preciso, especialmente cuando las partículas se movían muy rápido.
• La Nueva Forma: Este artículo calcula la trayectoria a través del laberinto real (con el peso real del quark top), pero solo para el primer paso del viaje (Orden Siguiente al Líder o NLO). Sin embargo, para las partes principales y complejas del viaje, utilizan el atajo del "suelo liso", pero calculándolo con un nivel de detalle sin precedentes.

2. El Avance: Calculando a "N3LO"

El artículo reporta las primeras predicciones totalmente diferenciales a N3LO (Orden Siguiente al Siguiente al Siguiente al Líder).

• La Analogía: Piensa en calcular el clima.
  • LO (Orden Líder): "Podría llover". (Suposición muy vaga).
  • NLO: "Lloverá por la tarde". (Mejor).
  • NNLO: "Lloverá a las 3 PM con 50% de humedad". (Muy bueno).
  • N3LO: "Lloverá a las 3:04 PM, con 50% de humedad, y las gotas golpearán el suelo con un ángulo de 45 grados". (Extremadamente preciso).
• Lo que hicieron: Calcularon el "clima" de la colisión del bosón de Higgs con este nivel extremo de precisión. No solo calcularon la cantidad total de "lluvia" (sección eficaz total); calcularon exactamente dónde y cómo cae (distribuciones diferenciales), como la velocidad y el ángulo de los bosones de Higgs.

3. Los Resultados: Un enfoque más nítido

• Reducción de la Incertidumbre: Antes de este artículo, el "pronóstico" tenía un gran margen de error (como decir "podría llover entre las 10 AM y las 6 PM"). Los nuevos cálculos de N3LO reducen esa ventana significamente, reduciendo la incertidumbre aproximadamente tres veces. Ahora, la predicción es lo suficientemente precisa como para estar al "nivel de porcentaje".
• La Forma de la Tormenta: Encontraron que, si bien la cantidad total de "lluvia" no cambió mucho, la forma de la tormenta sí lo hizo. Los nuevos cálculos cambian cómo se distribuyen los bosones de Higgs en términos de su velocidad y dirección. Esto es crucial porque si el "pronóstico" (teoría) no coincide con el "clima real" (experimento), podría significar que hay nueva física escondida en los datos.

4. Corrigiendo el atajo del "Top Pesado"

Dado que el atajo del "Top Pesado" no es perfecto cuando las partículas se mueven rápido, los autores combinaron sus cálculos de "suelo liso" ultra precisos con un cálculo más exacto del "laberinto real" para el primer paso.

• La Analogía: Imagina que tienes un mapa superdetallado de una ciudad (N3LO) pero sabes que el mapa es ligeramente erróneo respecto a la altura de los edificios. Tomas una foto de baja resolución del edificio real (NLO con masa real) y usas eso para corregir las alturas en tu mapa superdetallado.
• El Resultado: Este enfoque híbrido ofrece la imagen más precisa de la producción de pares de bosones de Higgs hasta la fecha. Encontraron que la "masa real" del quark top cambia significativamente las predicciones, especialmente para los bosones de Higgs que se mueven a altas velocidades o en direcciones específicas.

5. Por qué esto importa (según el artículo)

El artículo afirma que este nivel de precisión es esencial para los experimentos en curso en el LHC.

• El Objetivo: Los científicos buscan señales de que el potencial de Higgs (el campo de energía que otorga masa a las partículas) se comporta de manera diferente a lo que predice el Modelo Estándar.
• La Necesidad: Para encontrar estas diminutas diferencias, necesitas una "regla" (predicción teórica) que sea increíblemente precisa. Si tu regla es borrosa, no puedes distinguir si el objeto que estás midiendo es ligeramente diferente o si tu regla simplemente está mal. Este artículo proporciona una regla mucho más nítida.

En Resumen:
Este artículo es una clase magistral de precisión matemática. Toma un problema de física notoriamente difícil (dos bosones de Higgs creados a través de un bucle de quark top pesado) y lo calcula con la mayor exactitud posible disponible actualmente. Al refinar el "mapa" de cómo se comportan estas partículas, permite a los experimentalistas del LHC buscar nueva física con ojos mucho más agudos, reduciendo la "niebla" de la incertidumbre teórica que ha oscurecido la visión durante años.

Resumen técnico

Resumen Técnico: Producción de pares de bosones de Higgs totalmente diferencial a N3LO con efectos de la masa del quark top

Problema y motivación
La producción de pares de bosones de Higgs ($gg \to hh$) mediante fusión de gluones-gluones es el canal dominante para sondear el potencial del Higgs y el acoplamiento trilineal ($\lambda_{3h}$) en el Gran Colisionador de Hadrones (LHC). Aunque las correcciones de QCD de orden siguiente al líder (NLO) están disponibles con la dependencia completa de la masa del quark top, estas dejan incertidumbres de escala significativas (aprox. $\pm 13\%$) e introducen grandes incertidumbres teóricas relacionadas con el esquema de renormalización de la masa del quark top (hasta un 30% para ciertas distribuciones). Se han logrado cálculos de orden superior en el Límite de Quark Top Pesado (HTL, por sus siglas en inglés), donde $m_t \to \infty$, hasta el orden de Próximo-Siguiente-Siguiente-Orden Líder (N3LO) para secciones eficaces inclusivas y la distribución de la masa invariante di-Higgs. Sin embargo, las predicciones totalmente diferenciales a N3LO en el HTL, particularmente para otros observables cinemáticos, no estaban disponibles previamente o eran solo aproximadas. Además, la aproximación HTL deja de ser válida cuando la energía del centro de masas partónico se aproxima a $2m_t$, lo que requiere la inclusión de efectos de masa finita del quark top para una fenomenología fiable.

Metodología y marco computacional
Los autores realizan el primer cálculo de QCD a N3LO totalmente diferencial para $gg \to hh$ en el HTL, combinado con predicciones NLO que incluyen la dependencia completa de la masa del quark top ($N_{LO}^{mt}$).

1. Descomposición HTL: El cálculo en el HTL se organiza utilizando un lagrangiano efectivo donde los bucles de quarks top se integran en interacciones de contacto locales. La amplitud al cuadrado se descompone en tres clases basadas en el número de inserciones de vértices efectivos:

   • Clase-a: Dos vértices efectivos (dominante, $O(\alpha_s^2)$ en LO).
   • Clase-b: Tres vértices efectivos ($O(\alpha_s^3)$ en LO).
   • Clase-c: Cuatro vértices efectivos ($O(\alpha_s^4)$ en LO).
     Para alcanzar la precisión N3LO ($O(\alpha_s^5)$), los autores computan las correcciones N3LO para la Clase-a, las correcciones NNLO para la Clase-b, y las correcciones NLO para la Clase-c.

2. Técnicas computacionales:

   • Clase-a: Calculada utilizando el marco NNLOJET. Los autores relacionan la sección eficaz di-Higgs con la sección eficaz de un solo Higgs (conocida a N3LO) mediante un mapeo cinemático. Emplean el método de $q_T$-slicing para manejar las divergencias infrarrojas, separando el espacio de fase en una región de baja $q_T$ (computada mediante resúmen de $q_T$ en la Teoría de Colineales Blandos o SCET) y una región de alta $q_T$ (computada como $hh + \text{jet}$ a NNLO).
   • Clase-b: Las correcciones NNLO se computan utilizando una combinación de NNLOJET (para la parte de resúmen de baja $q_T$) y MadGraph5_aMC@NLO (para la parte de radiación real de alta $q_T$), utilizando el esquema de sustracción FKS.
   • Clase-c: Computada enteramente a NLO usando MadGraph5_aMC@NLO con sustracción FKS.
   • Validación: La implementación se valida contra resultados inclusivos conocidos de iHixs2 y resultados diferenciales de MadGraph5_aMC@NLO a órdenes inferiores, mostrando un acuerdo perfecto. Se verifica la independencia de los resultados respecto al parámetro de slicing $p_T^{\text{veto}}$, confirmando la cancelación de las correcciones de potencia.

3. Incorporación de efectos de la masa del quark top:
   Para dar cuenta de los efectos de la masa finita del quark top, los autores combinan los resultados N3LO del HTL con los resultados NLO con dependencia completa de $m_t$ (computados usando el generador ggxy). Evalúan tres esquemas de combinación:

   • Aditivo ($N_{kLO} \oplus N_{LO}^{mt}$): Suma las correcciones de orden superior del HTL al resultado NLO de $m_t$ completo.
   • Mejorado por Born ($N_{kLO}^{B-i} \oplus N_{LO}^{mt}$): Repondera las correcciones HTL mediante la razón de la sección eficaz LO de $m_t$ completo frente a la sección eficaz LO del HTL.
   • Multiplicativo ($N_{kLO} \otimes N_{LO}^{mt}$): Repondera las predicciones HTL mediante la razón de NLO de $m_t$ completo frente a NLO del HTL.
     El estudio identifica el esquema multiplicativo como el más estable perturbativamente para las distribuciones diferenciales.

Resultados clave
Los cálculos se presentan para colisiones protón-protón a $\sqrt{s} = 14$ TeV bajo cortes fiduciales realistas ($p_{T,h1} > 30$ GeV, $p_{T,h2} > 20$ GeV, $|y_h| < 2.4$).

1. Secciones eficaces fiduciales (HTL):

   • Las correcciones de QCD N3LO aumentan la sección eficaz fiducial NNLO en aproximadamente un 3.8%.
   • Las incertidumbres de escala se reducen significativamente: de $\pm 18\%$ en NLO a $\pm 5.2\%$ en NNLO, y hasta $+1.4\%/-2.9\%$ en N3LO. Esto representa una reducción de un factor de aproximadamente tres comparado con NNLO.
   • Las contribuciones de la Clase-a dominan ($\sim 103\%$ del total), mientras que la Clase-b y la Clase-c son insignificantes ($\sim -3.6\%$ y $\sim 0.03\%$, respectivamente).

2. Distribuciones diferenciales (HTL):

   • Masa invariante ($m_{hh}$): Las correcciones N3LO son relativamente planas ($K \approx 1.04$) excepto cerca del umbral ($m_{hh} \to 2m_h$), donde las incertidumbres de escala siguen siendo grandes.
   • Momento transversal ($p_T$): Los espectros de $p_T$ de los bosones de Higgs líder y sub-líder muestran modificaciones de forma significativas a N3LO. En la región de baja $p_T$, las predicciones de orden fijo son inestables (correcciones negativas, grandes incertidumbres), lo que indica la necesidad de resúmen. En la región de alta $p_T$ ($>100$ GeV), las correcciones N3LO modifican la forma de las distribuciones NNLO.
   • Ángulo azimutal: La distribución $|\phi_{h1} - \phi_{h2}|/\pi$ muestra una excelente convergencia entre NNLO y N3LO, con correcciones que alcanzan el $\sim 10\%$ en bins específicos.
   • Rapidez: Las distribuciones de rapidez muestran factores $K$ planos y las bandas de N3LO se encuentran dentro de las bandas de incertidumbre de NNLO.

3. Efectos de la masa finita del quark top:

   • La combinación del HTL N3LO con el NLO de $m_t$ completo (vía el esquema multiplicativo) revela que los efectos de la masa finita del quark top son esenciales incluso a masas invariantes bajas ($m_{hh} \approx 250$ GeV).
   • Comparado con el HTL, los resultados de $m_t$ completo muestran un aumento de un factor de dos en la región del umbral y una supresión en la cola de alta masa.
   • Los espectros de $p_T$ se suavizan en las colas cuando se incluye la dependencia completa de la masa.
   • La distribución del ángulo azimutal se suprime en más del $10\%$ en la región de retroceso (back-to-back) y por un factor de $\sim 1.7$ en la región colineal.
   • Las incertidumbres de escala para las predicciones combinadas son generalmente pequeñas, aunque la dependencia de escala de $N_{LO}^{mt}$ difiere ligeramente del caso HTL.

Significado y afirmaciones
El artículo afirma proporcionar las primeras predicciones totalmente diferenciales a N3LO para la producción de pares de bosones de Higgs en el HTL. Al combinar estas con efectos de masa de quark top a NLO, los autores presentan "una de las predicciones diferenciales a nivel de partones más precisas hasta la fecha" para las búsquedas en el LHC.

Las afirmaciones clave sobre el impacto de este trabajo incluyen:

• Reducción de la incertidumbre: Las correcciones de QCD a N3LO reducen las incertidumbres de escala de las predicciones fiduciales y diferenciales de NNLO en aproximadamente un factor de tres, llevando la incertidumbre teórica al nivel del percentil en el HTL.
• Modificaciones de forma: Los resultados demuestran que las correcciones N3LO no son meramente un factor de normalización; alteran significativamente las formas de las distribuciones de momento transversal y del ángulo azimutal, resaltando la necesidad de incluir correcciones N3LO para una fenomenología de precisión.
• Necesidad de la masa del top: El estudio confirma que los efectos de la masa finita del quark top no pueden ignorarse incluso a las masas invariantes más bajas y deben incluirse al orden más alto disponible (NLO en este contexto) para evitar distorsiones significativas en las distribuciones diferenciales.

Los autores señalan modestamente que, si bien las incertidumbres de escala están bajo control, otras incertidumbres teóricas permanecen, incluyendo las correcciones de masa finita de $m_t$ omitidas en NNLO y más allá, las incertidumbres del esquema de masa del quark top (OS vs. $\overline{MS}$), las correcciones electrodébiles de orden superior omitidas y las contribuciones de los bucles de quarks bottom. Sugieren que el trabajo futuro extendiendo la dependencia de masa completa a NNLO podría reducir aún más estas incertidumbres.