An Increase in the Galactic Planet Host Fraction Fails to Reproduce the Galactic Height Trend in Planet Occurrence

Resumen Técnico: Un aumento en la fracción de planetas anfitriones galácticos no logra reproducir la tendencia de altura galáctica en la ocurrencia de planetas

Planteamiento del problema
Resultados observacionales recientes, específicamente Zink et al. (2023) (en adelante Z23), han identificado una fuerte tendencia donde la tasa de ocurrencia de planetas pequeños (Super-Tierras y Sub-Neptunos) alrededor de enanas FGK disminuye significativamente con el aumento de la amplitud de oscilación vertical máxima ($Z_{max}$) desde el plano medio galáctico. Aunque se sabe que la metalicidad estelar se correlaciona con la formación planetaria, el gradiente de metalicidad galáctica por sí solo es insuficiente para explicar la pronunciación de este declive observado. La cuestión central abordada en este estudio es si un aumento en la fracción de anfitriones de planetas ($f$) dependiente del tiempo a lo largo de la historia de la Vía Láctea podría reproducir la tendencia de ocurrencia-$Z_{max}$ observada. Los autores investigan si un "impulso" en la eficiencia de la formación planetaria en algún momento del pasado de la Galaxia, combinado con el calentamiento cinemático de las órbitas estelares a lo largo del tiempo, puede explicar el déficit de planetas a grandes alturas galácticas.

Metodología
Los autores desarrollaron un nuevo código de síntesis de población de sistemas planetarios, psps, para simular poblaciones planetarias sintéticas alrededor de una muestra de enanas FGK. La metodología consistió en los siguientes pasos:

1. Muestras Estelares: Se construyeron dos muestras estelares distintas para probar la robustez de los resultados frente a las incertidumbres de edad:

   • Muestra Real: Una muestra conjunta de Kepler y K2 (22,184 estrellas) extraída del catálogo de Hardegree-Ullman et al. (2025), cotejada con edades de isocronas de Berger et al. (2020b, 2025). Los valores de $Z_{max}$ se calcularon utilizando la cinemática de Gaia DR3 y un modelo de potencial galáctico mediante el paquete Gala.
   • Muestra Sintética: Una muestra de síntesis de población TRILEGAL (49,109 estrellas) ajustada para asemejarse a los campos de Kepler y K2. Esta muestra proporcionó un escenario "optimista" con edades estelares arbitrariamente precisas y modelos de calentamiento cinemático conocidos, lo que permitió a los autores aislar los efectos de la precisión de la edad de los modelos físicos.

2. Modelos de Ocurrencia Dependientes del Tiempo: El estudio probó dos familias de modelos para la evolución de la fracción de estrellas que albergan planetas ($f$) como función de la edad estelar (tiempo de retroceso):

   • Funciones de Escalón: Un salto repentino de una fracción inicial $f_1$ a una fracción actual $f_2$ en un tiempo de retroceso $t$ específico.
   • Funciones por Partes: Un aumento lineal gradual de $f_1$ a $f_2$ comenzando en el tiempo $t$.
   • Todos los modelos fueron restringidos para producir una fracción de anfitriones de planetas actual de $\sim30-35\%$, consistente con las restricciones observacionales.

3. Síntesis de Población y Observación: Utilizando psps, los autores generaron sistemas planetarios basados en la $f$ prescrita para cada estrella. Se asignaron arquitecturas a los sistemas (sistemas multiplanetarios compactos vs. planetas únicos) y parámetros orbitales (excentricidad, inclinación) basados en restricciones demográficas existentes (ej. Lam & Ballard 2024). Las poblaciones sintéticas fueron "observadas" utilizando métodos de eficiencia de detección inversa (IDEM) que tuvieron en cuenta la probabilidad geométrica de tránsito y las funciones de sensibilidad específicas de las misiones Kepler y K2.

4. Análisis: Las tasas de ocurrencia simuladas se agruparon por $Z_{max}$ (100–1000 pc) y se ajustaron con una ley de potencia ($y \propto Z_{max}^\tau$) para determinar la pendiente ($\tau$) y la normalización ($\eta$). Estos resultados se compararon con la tendencia observada de Z23 ($\tau \approx -0.28$).

Resultados Clave

• Insuficiencia de los Aumentos en la Fracción de Anfitriones: El estudio encuentra que aumentar la fracción de anfitriones de planetas ($f$) a lo largo del tiempo, ya sea de forma instantánea o gradual, es generalmente insuficiente para reproducir la pendiente de la tendencia de ocurrencia-$Z_{max}$ observada.
• Restricciones del Modelo:
  • Para la muestra real de Kepler-K2, la gran mayoría de los modelos produjeron pendientes ($\tau$) demasiado tenues (más cercanas a 0) en comparación con el $\tau \approx -0.28$ observado. Los modelos que más se acercaron a coincidir con la tendencia requirieron parámetros extremos y físicamente implausibles (ej. $f_1 \approx 1\%$ o $f_2 \approx 100\%$) y tiempos de escalón ($t$) restringidos a una ventana estrecha entre 2.2 y 3.2 Gyr atrás.
  • Para la muestra sintética TRILEGAL (con edades precisas), los resultados fueron similares. Si bien un aumento tardío e instantáneo en $f$ podría coincidir con la pendiente de Z23, nuevamente requirió un ajuste hacia los extremos teóricos ($f_1=1\%$ o $f_2=100\%$). Esto sugiere que incluso con una precisión de edad perfecta, un simple aumento en $f$ no puede impulsar por sí solo la tendencia observada.
• Análisis Solo de Kepler: Al analizar únicamente el campo de Kepler (que tiene una distribución de edad diferente y carece de los datos de K2), la pendiente de ocurrencia-$Z_{max}$ derivada fue similar a la de la muestra conjunta. Un aumento de tipo escalón en $f$ coincidió ligeramente mejor con esta tendencia de solo Kepler, pero la dificultad fundamental de reproducir la pendiente persistió.
• Precisión de la Edad: La comparación entre la muestra real y la sintética indica que la incapacidad de los modelos para coincidir con la tendencia no se debe primordialmente a la falta de precisión en las mediciones de la edad estelar. Incluso con edades idealizadas, los modelos de "aumento gradual" fueron descartados, y los modelos de "escalón" requirieron parámetros irreales.

Significancia y Reivindicaciones
El artículo concluye que la fuerte correlación negativa observada entre la ocurrencia de planetas pequeños y la altura galáctica no puede explicarse únicamente por un aumento dependiente del tiempo en la fracción de estrellas que albergan planetas. Los autores afirman que, si bien un impulso tardío en $f$ (específicamente una función de escalón ocurrida hace 2–3 Gyr) puede reproducir marginalmente la tendencia, esto requiere un ajuste hacia los límites físicos del modelo (ej. casi el 100% de las estrellas albergan planetas hoy en día).

Los autores argumentan que esta discrepancia implica que otros factores deben estar en juego. Sugieren que los modelos futuros podrían lograr un mejor ajuste:

1. Diferenciando entre Super-Tierras y Sub-Neptunos (como encontró Z23, la tendencia es impulsada principalmente por los Sub-Neptunos).
2. Permitiendo que la fracción de sistemas multiplanetarios compactos "intactos" evolucione con el tiempo, en lugar de mantenerla constante.
3. Investigando mecanismos físicos que mejoren la formación planetaria reciente o que desestabilicen progresivamente los sistemas más antiguos (ej. mediante encuentros estelares cercanos, mareas galácticas o cambios en la turbulencia del medio interestelar).

El estudio enfatiza que la imagen canónica del enriquecimiento metálico local como único vínculo entre el contexto galáctico y la demografía planetaria es insuficiente en la era de Gaia. Los resultados resaltan la necesidad de conectar la historia dinámica a escala galáctica (fusiones, llegadas de satélites) con los procesos locales que afectan los discos protoplanetarios y la estabilidad de los sistemas para explicar plenamente la demografía de los exoplanetas.