A Bayesian framework for posterior bounds on self-replicating probes in the Solar System under observational incompleteness

Resumen Técnico: Un marco bayesiano para los límites posteriores de sondas autorreplicantes en el Sistema Solar

Planteamiento del problema
A pesar de décadas de búsquedas sistemáticas, no se ha reportado evidencia robusta de artefactos extraterrestres (sondas autorreplicantes) dentro del Sistema Solar. Mientras que la paradoja de Fermi ha sido analizada extensamente con respecto a la prevalencia de civilizaciones en la Galaxia, existe una falta de una inferencia unificada y de extremo a extremo que traduzca cuantitativamente la ausencia de detecciones en regiones específicas del Sistema Solar en una restricción falsable sobre el número esperado ($\lambda$) de sondas residentes. Los marcos existentes se centran o bien en los parámetros de la ecuación de Drake para civilizaciones galácticas (por ejemplo, Sandberg et al., 2018) o discuten cualitativamente los objetivos de búsqueda sin integrar datos heterogéneos de completitud de sondeos en una única distribución posterior.

Metodología
Los autores construyen un marco de inferencia bayesiana que combina un prior paramétrico con una verosimilitud de Poisson derivada de resultados observacionales nulos.

1. Construcción del Prior: El modelo adopta un prior log-uniforme sobre $\lambda$, derivado de un producto de factores tipo ecuación de Drake (tasa de formación estelar, fracción de planetas habitables, probabilidad de inteligencia, etc.) extendido para incluir parámetros específicos de las sondas ($f_{probe}$ y el número en estado estacionario de sondas por sistema $\langle n_{probe/sys} \rangle$). Esto resulta en un prior de cola pesada con un límite superior del 95% de $\lambda_{prior} \approx 230$.
2. Función de Verosimilitud: La verosimilitud se modela como un evento de detección cero de un solo ensayo de Poisson: $L(\lambda) \propto e^{-\lambda p_{det}}$. Aquí, $p_{det}$ representa la probabilidad efectiva de que una sola sonda, si estuviera presente, habría sido detectada por los sondeos existentes.
3. Descomposición de la Probabilidad de Detección ($p_{det}$): El marco descompone $p_{det}$ a través de siete regiones distintas del Sistema Solar: Objetos Cercanos a la Tierra (NEOs), co-orbitales de la Tierra, puntos de Lagrange Tierra-Luna, la superficie lunar, el cinturón principal de asteroides, el cinturón de Kuiper y el Sistema Solar exterior/misceláneo.
   • Para cada región, se construye una curva de completitud dependiente del tamaño $p^{(i)}(d)$ mediante la inversión o interpolación de datos de eficiencia de sondeo publicados (por ejemplo, Harris & Chodas para NEOs, Gregg & Wiegert para puntos de Lagrange).
   • Se aplica un escenario de peso de residencia ($w_i$) para ponderar la contribución de cada región al $p_{det}$ total. Se prueban tres escenarios: "Default" (dominado por el Sistema Solar exterior), "Diffuse" (distribución de volumen uniforme) y "Lurker" (concentrado en co-orbitales de la Tierra).
4. Implementación: La inferencia se ejecuta mediante muestreo de Monte Carlo ($N = 2 \times 10^5$) utilizando un paquete de Python (seta_bayesian). El código incluye una suite de pruebas rigurosa que verifica los límites del prior, la monotonicidad de la completitud y la convergencia.

Resultados Clave
El estudio reporta los límites superiores posteriores del 95% para $\lambda$ para tres tamaños de sonda canónicos (10 m, 100 m y 1 km) bajo los tres escenarios de residencia, basados en los datos disponibles hasta 2025.

• Restricciones Actuales (2025):
  • Sondas de 10 m: Bajo el modelo de residencia por defecto, el límite superior del 95% es $\lambda < 14.3$. El límite es laxo porque la probabilidad de detección efectiva es baja ($p_{det} \approx 0.0019$), impulsada por la suposición de que la mayoría de las sondas residen en el Sistema Solar exterior donde la completitud del sondeo es insignificante para este tamaño.
  • Sondas de 100 m: El límite se estrecha a $\lambda < 3.16$.
  • Sondas de 1 km: El límite se estrecha significativamente a $\lambda < 0.68$. A esta escala, el posterior está "saturado por el prior" porque los sondeos existentes (NEO, Cinturón Principal) ya están cerca de la completitud ( $>90\%$) para objetos de este tamaño; los datos añaden poca información más allá de la cola superior del prior.
• Sensibilidad: El escenario "Lurker" (sondas en co-orbitales de la Tierra) produce los límites más estrictos (por ejemplo, $\lambda < 2.27$ para 10 m), mientras que el escenario "Diffuse" produce resultados intermedios. El escenario "Default" produce los límites más laxos debido al alto peso asignado a las regiones exteriores insuficientemente sondeadas.
• Pronósticos Futuros (2025–2035): El marco proyecta el impacto de las próximas misiones (LSST, NEO Surveyor, Telescopio Roman).
  • Para sondas de 10 m, se prevé que el límite de residencia por defecto se comprima por un factor de 7.6× (estrechando el límite de ~14.3 a ~1.9) para 2035. Esta compresión es altamente sensible a la realización de los objetivos de completitud por debajo de los 100 m (intervalo: 4.8×–10.0×).
  • Para sondas de 1 km, el límite se mueve solo un ~14–16%, confirmando que los sondeos a escala de km ya están en el régimen dominado por el prior.

Significancia y Reivindicaciones
Los autores declaran explícitamente que la contribución principal no es una probabilidad definitiva de que las sondas existan, sino un marco reproducible y basado en datos que:

1. Cuantifica los Resultados Nulos: Convierte los datos de sondeo heterogéneos de no detección en una restricción posterior rigurosa sobre la abundancia de sondas.
2. Audita los Supuestos: Mediante el uso de priors fijados por YAML y una procedencia explícita para cada entrada de completitud, el marco hace que los supuestos del prior y las fuentes de datos sean transparentes y auditables.
3. Identifica el Canal de Descubrimiento Crítico: El análisis revela que la capacidad de la próxima década para restringir la existencia de sondas depende casi enteramente del descubrimiento de sub-100 m (específicamente el régimen de 10 m). Los sondeos que apuntan a objetos de escala de km ofrecen rendimientos decrecientes porque el prior ya satura el límite en esos tamaños.
4. Define la Estructura de Regímenes: El marco distingue entre un régimen "saturado por el prior" (objetos grandes, donde los datos añaden poco) y un régimen "restringido por los datos" (objetos pequeños, donde la mejora en la completitud del sondeo estrecha directamente el límite).

El artículo concluye que, si bien los datos actuales no pueden descartar la existencia de sondas pequeñas, el marco proporciona un límite superior falsable que se estrechará significativamente si los próximos sondeos logran sus objetivos de completitud para objetos menores a 100 metros.