An approximate formula for the entropy of the negative binomial distribution

Resumen Técnico: Una Fórmula Aproximada para la Entropía de la Distribución Binomial Negativa

Planteamiento del Problema
La Distribución Binomial Negativa (DBN) sirve como la parametrización estándar para las multiplicidades de partículas cargadas en la física de partículas de altas energías, aplicable a sistemas que van desde la dispersión inelástica profunda ($e^-p$) y colisiones protón-protón ($pp$) hasta interacciones de iones pesados. Desarrollos teóricos recientes han establecido un vínculo potencial entre la entropía de entrelazamiento del estado inicial y la entropía de Shannon del estado final de estas distribuciones. En consecuencia, calcular la entropía de Shannon de la DBN es de gran interés para la comunidad. Sin embargo, no existe una expresión analítica cerrada, simple y exacta para esta entropía. Si bien son posibles evaluaciones numéricas utilizando representaciones integrales (específicamente la Ecuación (23) en la Ref. [20]), estas requieren una integración numérica computacionalmente intensiva, lo que limita su utilidad para obtener una comprensión analítica rápida o realizar escaneos amplios de parámetros.

Metodología
El artículo aborda la falta de una solución de forma cerrada derivando una fórmula analítica aproximada para la entropía de la DBN. El enfoque procede de la siguiente manera:

1. Reformulación de la Integral: Los autores comienzan con la representación integral conocida de la entropía de la DBN (Ecuación 1). Aíslan el componente integral, $I(k, \langle n \rangle)$, que depende de los parámetros de la DBN: la multiplicidad media $\langle n \rangle$ y el parámetro de dispersión $k$.
2. Simplificación Analítica: Utilizando la identidad $(1-z)^{k-1}-1 = \ln(1-z) \int_0^{k-1} (1-z)^t dt$, los autores separan la integral en un término que puede evaluarse analíticamente ($\langle n \rangle \ln k$) y una integral restante, $J(k, \langle n \rangle)$, que aún carece de forma cerrada.
3. Transformación de Variables: Para facilitar la aproximación, se aplica un cambio de variable $u = -\ln(1-z)$, transformando el dominio de integración de $[0, 1]$ a $[0, \infty)$.
4. Aproximación Exponencial: El núcleo de la metodología consiste en aproximar el término complejo entre corchetes dentro de la integral transformada. Los autores reemplazan el término $\left(1 + \frac{\langle n \rangle}{k}(1-e^{-u})\right)^{-k} - 1$ con una única función exponencial de saturación de la forma $-D(1 - e^{-\lambda u})$.
5. Derivación de la Forma Cerrada: Esta sustitución permite evaluar la integral restante analíticamente en términos de la función Gamma ($\Gamma$). Los parámetros $D$ y $\lambda$ se definen explícitamente en función de $\langle n \rangle$ y $k$ (Ecuación 12).

Contribuciones y Resultados Clave
La contribución principal es la derivación de una fórmula aproximada de forma cerrada para la entropía de Shannon de la DBN, presentada como la Ecuación (13):

$$S_{NBD} \approx (\langle n \rangle + k) \ln(\langle n \rangle + k) - (\langle n \rangle \ln \langle n \rangle + k \ln k) - D \ln \left[ \frac{\Gamma(k + \lambda)}{\Gamma(k)\Gamma(1 + \lambda)} \right]$$

donde $D = 1 - (1 + \langle n \rangle/k)^{-k}$ y $\lambda = \langle n \rangle / D$.

La precisión de esta aproximación se validó comparando la entropía aproximada ($S_{approx}$) con la entropía "exacta" ($S_{exact}$) obtenida mediante integración numérica de la fórmula original. Se encontró que la diferencia relativa, $\delta S = (S_{approx} - S_{exact}) / S_{exact}$, es:

• Típicamente alrededor del 10% dentro del dominio de parámetros estándar de la DBN.
• Dentro del 20% para valores extremos de $\langle n \rangle$ y $k$ (específicamente en casos de sobredispersión).

Se proporciona un escaneo detallado del espacio de parámetros en la Figura 1 del artículo, ilustrando la distribución de estos errores.

Significado y Afirmaciones
El artículo afirma modestamente que, si bien una determinación precisa de la entropía de la DBN sigue requiriendo la implementación numérica de la integral original, la fórmula derivada ofrece una alternativa valiosa cuando la precisión de $\sim10\%$ (o hasta $\sim20\%$ en extremos) es satisfactoria. La fórmula proporciona una expresión de forma cerrada que involucra únicamente funciones estándar (logaritmos y funciones Gamma), evitando así la necesidad de integración numérica. Esto facilita una manipulación analítica más sencilla y una evaluación más rápida en contextos donde la entropía de la DBN es un observable central, particularmente en estudios que vinculan las multiplicidades de partículas con la entropía de entrelazamiento. Los autores señalan explícitamente que no se conocía ninguna solución de forma cerrada previamente en la literatura.