Explicit bounds for Dickman's function

Resumen Técnico: Límites Explícitos para la Función de Dickman

Planteamiento del Problema
La función de Dickman, $\rho(u)$, se define como la solución continua de la ecuación diferencial de retardo $u\rho'(u) + \rho(u-1) = 0$ con la condición inicial $\rho(u) = 1$ para $0 \le u \le 1$. Representa la densidad natural de los enteros $n$ cuyo mayor factor primo no supera a $n^{1/u}$. Si bien el comportamiento asintótico de $\rho(u)$ cuando $u \to \infty$ está bien comprendido mediante el trabajo de de Bruijn y Alladi, estos resultados suelen involucrar términos de error implícitos (por ejemplo, $O(1/u)$). El principal desafío abordado en este artículo es establecer límites superiores e inferiores numéricamente explícitos para $\rho(u)$ que permitan una evaluación aproximada sin necesidad de resolver numéricamente la ecuación diferencial de retardo.

Metodología
El autor emplea el método del punto de silla aplicado a la representación de la integral de Laplace inversa de $\rho(u)$. El enfoque refina la técnica de prueba encontrada en el trabajo de Tenenbaum, el cual estableció el resultado de Alladi con un error relativo de $O(u^{-1})$. Para lograr los límites explícitos más ajustados requeridos aquí, la metodología implica:

1. Seguimiento Preciso del Error: Mantener el registro de todos los factores constantes en lugar de absorberlos en la notación de orden grande (big-O).
2. Integración Refinada: Dividir el rango de integración en intervalos adicionales para aislar las contribuciones de manera más precisa.
3. Expansión de Orden Superior: Utilizar un polinomio de Taylor de orden siete (en lugar de tres) cerca del punto de silla para capturar términos de orden superior en la expansión del integrando.
4. Lemas Explícitos: Derivar desigualdades específicas para las derivadas de la función integral $I(s) = \int_0^s \frac{e^z-1}{z} dz$ y acotar las contribuciones de la cola de la integral utilizando lemas relacionados con la función gamma incompleta.

La función $\tilde{\rho}(u)$, que sirve como la aproximación principal, se define como:
$$\tilde{\rho}(u) := \frac{\exp\{\gamma - u\xi + I(\xi)\}}{\sqrt{2\pi I''(\xi)}}$$
donde $\xi = \xi(u)$ es la única solución positiva de $e^\xi = 1 + u\xi$, y $\gamma$ es la constante de Euler.

Contribuciones Clave y Resultados
El artículo establece dos teoremas principales que proporcionan límites explícitos y estimaciones para $\rho(u)$:

• Teorema 1 (Límites Explícitos): Para $u > 1$, la función $\rho(u)$ está acotada por:
  $$\tilde{\rho}(u) \left(1 - \frac{1}{12u}\right) < \rho(u) < \tilde{\rho}(u) \left(1 - \frac{1}{12u(1 + \frac{1}{\log u})}\right)$$
  Esto proporciona un rango conveniente para $\rho(u)$ utilizando únicamente funciones elementales y la función W de Lambert (para resolver $\xi$).

• Teorema 2 (Estimación Refinada con Término de Error): Para $u > 1$, la función puede expresarse como:
  $$\rho(u) = \tilde{\rho}(u) \left(1 + \alpha(u) + \frac{\theta(u)}{u^2}\right)$$
  donde $\alpha(u)$ es una expresión explícita que involucra las derivadas de $I(\xi)$:
  $$\alpha(u) := \frac{1}{8}\frac{I^{(4)}(\xi)}{(I''(\xi))^2} - \frac{5}{24}\frac{(I^{(3)}(\xi))^2}{(I''(\xi))^3}$$
  Crucialmente, el término de error está acotado explícitamente: $|\theta(u)| < 0.005$ para todo $u \ge 5$. Esto resulta en un error relativo de menos de $0.005/u^2$.

El artículo también proporciona un método para el cálculo numérico de $\rho(u)$ para verificación:

• Para $1 \le u \le 50$, utiliza el algoritmo de polinomio de Taylor de Marsaglia, Zaman y Marsaglia.
• Para valores de $u$ mayores, evalúa la integral de Laplace inversa numéricamente, observando que el integrando converge a la función de densidad normal estándar.

Significancia y Reivindicaciones
El artículo afirma que estos resultados permiten la "evaluación aproximada de $\rho(u)$ sin la necesidad de resolver numéricamente la ecuación diferencial de retardo". Al proporcionar constantes explícitas y límites de error, el trabajo hace que las estimaciones asintóticas de de Bruijn y Alladi sean "enteramente explícitas".

El autor señala que, si bien otros trabajos (Smida, Xuan) proporcionan estimaciones con magnitudes de error menores, carecen de explicitud numérica. Los límites explícitos derivados aquí son particularmente útiles para aplicaciones que involucran $\Psi(x, y)$, el conteo de enteros hasta $x$ cuyo mayor factor primo es $\le y$. Específicamente, combinados con resultados condicionales de Gorodetsky (asumiendo la Hipótesis de Riemann), estos límites producen límites inferiores condicionales para $\Psi(x, y)$ que son "fáciles de evaluar numéricamente".

El artículo valida sus resultados comparándolos con los valores obtenidos por van de Lune y Wattel, quienes utilizaron la regla del trapecio para aproximar $\rho(u)$ con al menos cinco dígitos significativos para $u \le 1000$. Los valores derivados de los Teoremas 1 y 2 son consistentes con estos referentes numéricos establecidos.