Validación del modelo HDOV en el cruce de la heliopausa por Voyager 1 y Voyager 2

Validación del modelo HDOV en el cruce de la heliopausa por Voyager 1 y Voyager 2 Arnoldo Walter Fernández [email protected] PREPRINT — 22 de noviembre de 2025 Resumen Este trabajo presenta una validación del modelo de Dispersión de Onda Vibracional (HDOV) aplicada al cruce de la heliopausa por las sondas Voyager 1 y Voyager 2. A partir de datos reales de los instrumentos MAG y PWS, se define un perfil dinámico de la función κlocal ( t ), interpretada como una impedancia efectiva del medio o tasa efectiva de amortiguamiento de coherencia en mesoescala, cuya evolución permite identificar eventos anómalos que exhiben una anticipación estadística respecto de los cambios observados por modelos magnetohidrodinámicos convencionales. La transición abrupta en κlocal precede a discontinuidades registradas en la intensidad del plasma y en la magnitud del campo magnético, indicando una capacidad predictiva operativa de la transición de heliopausa en el entorno local de las sondas. La comparación entre simulaciones HDOV y modelos MHD muestra una correlación estructural entre el campo magnético, la presión de plasma y la función f ( t )de coherencia del medio en mesoescala. La validación cruzada con Voyager 2 confirma la robustez y transferibilidad del indicador dentro de las incertidumbres observacionales y del reducido número de eventos disponibles. En este trabajo κlocal se presenta como un índice fenomenológico de cambio de régimen en plasmas turbulentos; una interpretación más amplia en el marco HDOV, que conecta estas transiciones con nociones de coherencia gravitacional y cuántica, se discute en trabajos complementarios y excede el alcance estrictamente observacional de este preprint. 1 Índice 1. Introducción y objetivos 5 1.1. Motivación ................................... 5 1.2. Objetivos ................................... 5 1.3. Contexto y marco conceptual ........................ 6 1.4. Estructura del trabajo ............................ 6 2. Adquisición y preprocesamiento de datos 7 2.1. Fuentes y formatos de datos ......................... 7 2.2. Definición operativa de κlocal(t)....................... 7 2.3. Unificación y exportación ........................... 9 3. Definición de κlocal(t)9 3.1. Formulación lineal de κlocal .......................... 9 3.2. Validación dimensional ............................ 10 3.3. Calibración e incertidumbre de los parámetros ............... 10 3.3.1. Anclaje físico de los coeficientes β3yβ4............... 12 3.4. Análisis de robustez con transformada wavelet continua .......... 12 3.5. Perfiles físicos y métrica de amortiguamiento lenta en malla 1 h ..... 13 3.6. Comparación cuantitativa con modelos MHD y otros enfoques ...... 15 3.6.1. Comparación directa con proxies MHD en malla temporal común 16 3.6.2. Análisis de robustez multi-escala ................... 17 3.6.3. Análisis de valor agregado y ventajas operativas .......... 18 3.6.4. Implicaciones para monitoreo operacional .............. 18 3.7. Extensión a Voyager 2: Construcción universal de κlocal .......... 18 4. Interpretación física del mecanismo detrás de κlocal(t) 19 4.1. Interpretación física de κlocal ......................... 19 4.1.1. Definición operacional de coherencia de fase ............ 20 4.1.2. Conexión con modelos físicos establecidos .............. 20 4.1.3. Derivación física de los coeficientes βi................ 20 4.1.4. Base física para la elección de la escala L.............. 21 4.1.5. Validez del proxy de gradiente y régimen de Taylor ........ 21 4.1.6. Vínculo con observables físicos directos ............... 22 4.1.7. Mecanismos físicos de la predictividad de κlocal ........... 22 4.1.8. Componentes: térmico y magnético ................. 23 4.1.9. Analogía de coherencia con diagramas de fases ........... 24 4.1.10. Emergencia como ruptura de coherencia .............. 24 4.2. Limitaciones y consideraciones ........................ 24 4.2.1. Limitaciones teóricas ......................... 25 4.2.2. Limitaciones observacionales ..................... 25 4.2.3. Limitaciones estadísticas ....................... 26 4.2.4. Limitaciones de generalización .................... 26 4.2.5. Contexto y perspectiva ........................ 26 5. Síntesis de resultados principales 27 2 6. Validación de κlocal(t) con datos reales 27 6.1. Correlación con eventos abruptos reales ................... 28 6.2. Interpretación de correlaciones modestas en sistemas complejos ...... 28 6.2.1. Comparación con umbrales en física de plasmas .......... 29 6.3. Validación extendida con proxies de partículas ............... 30 6.4. Limitaciones y consideraciones ........................ 30 6.5. Limitaciones del presente estudio ...................... 31 7. Mecanismo de coherencia de κlocal desde primeros principios 32 7.1. Interpretación funcional en el marco HDOV generalizado ......... 32 7.2. Conexión con teoría generalizada ...................... 33 7.2.1. Ecuación maestra derivada de la acción efectiva .......... 33 7.3. Propiedades de κlocal ............................. 33 7.3.1. Unidades y normalización ...................... 33 7.4. Control de calidad y verificación ....................... 33 8. Implicaciones para Física de Plasmas 33 8.1. Teoría de Turbulencia en Plasmas ...................... 34 8.2. Diagnóstico de Plasmas Astrofísicos ..................... 34 8.3. Conexión con Fundamentos de Física de Plasmas ............. 34 8.4. Perspectivas para Investigación Futura ................... 35 8.5. Generalización a Otros Entornos de Plasma ................ 35 8.5.1. Magnetopausa Terrestre ....................... 36 8.5.2. Tokamaks (Región de Transporte) .................. 36 8.5.3. Corona solar .............................. 36 9. Forecast operativo (v2) 36 9.1. Curvas ROC/PR y Métricas de Desempeño ................ 37 9.2. Estudio de Ablación de Características ................... 37 10.Conclusiones y perspectivas 40 A. Derivación desde primeros principios de κlocal 40 A.1. Ecuación de onda en plasmas no homogéneos ................ 40 A.2. Ansatz de WKB y separación de escalas .................. 41 A.3. Promediado sobre escalas de mesoescala .................. 41 A.3.1. Término temporal ........................... 41 A.3.2. Término espacial ........................... 41 A.3.3. Término de amortiguamiento .................... 41 A.4. Definición de κlocal .............................. 42 A.5. Conexión con µ(z)yn(p).......................... 42 A.6. Validación dimensional (detalles) ...................... 43 A.7. Conclusión del apéndice ........................... 43 B. Reproducibilidad y validación de robustez 43 B.1. Reproducibilidad ............................... 43 B.1.1. Datos y fuentes ............................ 43 B.1.2. Procedimiento de procesamiento ................... 44 B.1.3. Código y verificación ......................... 44 3 B.2. Análisis de robustez multi-escala ....................... 44 B.2.1. Barrido multi-L............................ 44 B.2.2. Sensibilidad a parámetros ...................... 45 B.3. Validación cruzada con Voyager 2 ...................... 45 B.3.1. Metodología y procesamiento .................... 45 B.3.2. Resultados principales ........................ 45 B.3.3. Análisis de amenazas a la validez .................. 45 B.4. Análisis de incertidumbre instrumental ................... 46 B.5. Conclusiones del apéndice .......................... 46 C. Órdenes de magnitud y cotas físicas 46 C.1. Parámetros de referencia en el entorno de heliopausa ........... 46 C.2. Escalas temporales y espaciales relevantes ................. 47 C.3. Cotas esperadas para κlocal .......................... 47 C.4. Consistencia con la interpretación HDOV .................. 47 4 1. Introducción y objetivos La validación del modelo HDOV (Hipótesis de Dispersión de Onda Vibracional) en contextos reales constituye un paso esencial para diferenciarlo de otras aproximaciones teóricas, como los modelos hidromagnéticos (MHD) o interpretaciones geométricas estándar del medio interestelar (Fernandez,2025b). En este trabajo, proponemos aplicar HDOV al análisis del cruce de la heliopausa por parte de la sonda Voyager 1 (período 2012-2013). Además de replicar el hallazgo en Voyager 1, aplicamos HDOV a Voyager 2 y encontramos una anticipación CUSUM del cruce del 5-Nov-2018 hacia el 25-Sep-2018 18:00 UTC (+40.2 d), con el mismo signo y orden de magnitud que en V1 bajo idénticas convenciones de umbral y métrica. 1.1. Motivación El modelo HDOV plantea que ciertas transiciones de régimen en los estados del medio pueden anticiparse a través de un parámetro local κlocal derivado de propiedades magnéticas y térmicas del entorno (Fernandez,2025b). En este sentido, la heliopausa, entendida como la frontera entre el viento solar y el medio interestelar, representa un caso paradigmático donde las condiciones cambian bruscamente y en escalas de coherencia coherentes. La evidencia experimental de Gurnett et al. (2013a)yBurlaga et al. (2013a) sugiere discontinuidades que no coinciden exactamente con las predicciones MHD. Esto abre la posibilidad de reinterpretar los datos desde una perspectiva de coherencia basada en HDOV, en la cual la pérdida de coherencia es anticipada por la dinámica de κlocal(t). 1.2. Objetivos Los objetivos principales de este estudio son: 1. Definir y calcular la función κlocal ( t )a partir de datos reales de Voyager 1 (MAG, PWS y CRS), construyendo una métrica de coherencia lenta del entorno heliosférico local. 2. Comparar la evolución de κlocal ( t )con eventos abruptos detectados en MAG, PWS y CRS, caracterizando su relación con discontinuidades en la coherencia del plasma y la estructura de la heliopausa. 3. Evaluar la capacidad predictiva del modelo HDOV respecto de las discontinuidades en la coherencia del entorno, cuantificando adelantos temporales (∆ tlead ) en el cruce de la heliopausa y comparando con proxies MHD clásicos. 4. Extender la construcción y calibración de κlocal ( t )a Voyager 2 y comprobar la reproducibilidad de los resultados en una segunda trayectoria de heliopausa, incluyendo la validación cruzada V1–V2 de los leads y métricas de ajuste. 5. Contrastar las predicciones del enfoque HDOV con modelos MHD y parametrizaciones empíricas alternativos, identificando fortalezas y limitaciones del marco funcional propuesto para la descripción de la heliopausa. 5 1.3. Contexto y marco conceptual En trabajos previos (Fernandez,2025b), se ha propuesto que la dispersión de ondas en plasma permite definir una estructura de coherencia del espacio-tiempo basada en la coherencia del entorno plasmático local. Esta idea se implementa en HDOV a través de la función κ ( z )o κ ( t ), la cual representa una “densidad de onda efectiva” o “impedancia del medio” al colapso de estados. En contextos astrofísicos, esta función puede adquirir un rol determinante para distinguir fases coherentes y no con coherencia mantenida. Particularmente, cuando se observa un máximo local o discontinuidad en κlocal , el modelo HDOV predice la emergencia de una nueva fase espacial, aún si los parámetros macroscópicos como presión o densidad no presentan discontinuidades visibles (Fernandez,2025b). 1.4. Estructura del trabajo Este artículo está organizado de la siguiente manera. En la Sección 2 se describe la adquisición y el preprocesamiento de los datos de campo magnético (MAG) y densidad electrónica (PWS) de Voyager 1, junto con la malla temporal adoptada, las normalizaciones empleadas y los criterios de limpieza y control de calidad. En la Sección 3 se introduce la definición operativa de κlocal ( t ), se discute su interpretación como métrica de amortiguamiento lento y se presenta el modelo lineal en términos de observables locales. Allí se detalla la política de ajuste de los coeficientes βi , se valida el cierre dimensional y se muestran los perfiles físicos y de κlocal en malla de 1 h, incluyendo la comparación cuantitativa con proxies MHD y la extensión del análisis a Voyager 2. La Sección 4 desarrolla la interpretación física del mecanismo detrás de κlocal ( t ), analizando el papel combinado de colisiones débiles, gradientes de campo y turbulencia en la decoherencia de fase, así como las principales limitaciones del modelo en este dominio de parámetros y su contexto dentro de la teoría de plasmas. En la Sección 5 se presenta una síntesis estructurada de los resultados principales y de su relevancia conceptual como indicador de accesibilidad funcional del entorno heliosférico local. La Sección 6 cuantifica la correlación entre κlocal ( t )y eventos abruptos observados por MAG, PWS y proxies de partículas, discutiendo el significado físico de correlaciones modestas pero robustas en desfase y su comparación con umbrales habituales en física de plasmas. La Sección 7 deriva una formulación de κlocal ( t )desde primeros principios dentro del marco HDOV generalizado y analiza la interpretación funcional de cada término del modelo, enfatizando el mecanismo de coherencia y su conexión con la teoría general. En la Sección 8 se discuten las implicaciones del enfoque para plasmas astrofísicos y de laboratorio, incluyendo su carácter multiescala y posibles extensiones a otros entornos (magnetopausa terrestre, tokamaks, corona solar). La Sección 9 explora el potencial de κlocal ( t )como indicador de forecast operativo, presentando curvas ROC/PR y experimentos de ablación sobre los componentes del modelo en Voyager 2. Finalmente, en la Sección 10 se presentan las conclusiones generales, las principales limitaciones del estudio y las perspectivas para trabajos futuros. El Apéndice A desarrolla la derivación de κlocal desde primeros principios, el Apéndice B documenta los aspectos de reproducibilidad y robustez del pipeline de análisis, y el Apéndice C recoge el anexo 6 dimensional y físico, con las unidades, normalizaciones, órdenes de magnitud y cotas para los coeficientes del modelo que garantizan la consistencia física de la formulación. Nota. Para este trabajo se ha intentado extremar el cuidado en el tratamiento de datos, la estadística y la documentación, pero el marco HDOV aplicado a la heliopausa sigue siendo una propuesta exploratoria que, al momento de esta versión, no ha pasado por revisión por pares en revistas especializadas. Publicamos los datos, scripts y figuras para que cualquiera pueda: verificar cada paso, criticar la metodología, reproducir y, si es necesario, refutar los resultados. 2. Adquisición y preprocesamiento de datos Para validar el modelo HDOV en el entorno interestelar inmediato, se utilizaron los datos de campo magnético (MAG) y de ondas de plasma (PWS) proporcionados por la sonda Voyager 1, correspondientes al período comprendido entre julio de 2012 y febrero de 2013. Este intervalo incluye el cruce de la heliopausa, identificado en la literatura entre agosto y septiembre de 2012 (Gurnett et al.,2013a;Burlaga et al.,2013a). 2.1. Fuentes y formatos de datos Los datos fueron descargados desde el portal CDAWeb de la NASA 1 , utilizando el script vg1_preprocess_from_cdaweb.py , desarrollado específicamente para este trabajo. Se procesaron los siguientes productos: MAG: Datos de campo magnético con resolución de 48 segundos, convertidos a series horarias por remuestreo. PWS-LR: Densidades de plasma derivadas de la frecuencia del pico espectral, con resolución de 1 segundo, también remuestreadas a 1 hora. 2.2. Definición operativa de κlocal(t) La métrica central de este trabajo, κlocal ( t ), se define operativamente como la tasa efectiva de amortiguamiento o pérdida de coherencia de fase en el plasma. Su construcción se basa en el concepto de tiempo de decorrelación de fase τc y la longitud de coherencia Lc de los modos de onda de Alfvén, siguiendo la formulación del modelo HDOV Fernandez (2025a). La definición fundamental es: κlocal ≡1 τc =vA Lc ,(1) donde el tiempo de decorrelación τc se define como el intervalo en el que la autocorrelación de fase Cϕ(τ)decae a 1/e de su valor inicial: Cϕ(τc) = 1 e.(2) 1https://cdaweb.gsfc.nasa.gov/ 7 Formalmente, para un campo magnético filtrado en banda estrecha Bf ( t )que contiene un modo dominante de Alfvén, la fase instantánea ϕ ( t )se puede definir mediante la transformada de Hilbert H(·): ϕ(t) = arg [H(Bf(t))] .(3) Esta construcción establece un puente conceptual entre coherencia de fase y tiempo de decorrelación τc . Sin embargo, en el régimen turbulento de banda ancha de la heliopausa no utilizamos esta fase de Hilbert como observable directo sobre los datos reales de Voyager, ya que su interpretación física se degrada en presencia de espectros ampliamente distribuidos en frecuencia. En este trabajo, τc y κlocal se estiman de manera operacional a partir de promedios de |B| , ne , |∇B|/B yΓ total sobre ventanas de mesoescala L , según la Ecuación (5) , que actúan como proxies robustos de la pérdida de coherencia en el plasma. La velocidad de Alfvén vAse calcula a partir de los parámetros locales del plasma: vA=B √µ0ρ,(4) donde B es la magnitud del campo magnético, µ0 es la permeabilidad del vacío y ρ es la densidad de masa del plasma. Operacionalmente, κlocal ( t )se construye como una combinación lineal de los promedios temporales de las principales variables físicas que modulan la coherencia del plasma: κlocal =β1|B|τ+β2neτ+β3|∇B| Bτ +β4Γtotalτ.(5) Cada término de la Ecuación (5) tiene un significado físico claro: β1⟨|B|⟩τ : contribución del campo magnético. Un campo más intenso tiende a estructurar el plasma, pero también puede aumentar la anisotropía y reducir la coherencia en ciertas escalas. β2⟨ne⟩τ : captura los efectos térmicos y de densidad. Mayores densidades de electrones elevan la frecuencia del plasma y modulan la capacidad del medio para sostener modos coherentes. β3|∇B|/Bτ : cuantifica la influencia de los gradientes magnéticos y la turbulencia. Gradientes pronunciados indican inhomogeneidades que dispersan las ondas y reducen la coherencia. β4⟨ Γ total⟩τ : incorpora otras fuentes de amortiguamiento, como efectos colisionales y cinéticos. Los coeficientes βi se calibraron empíricamente mediante minimización del error cuadrático medio (RMSE) entre κlocal ( t )y un perfil de referencia, utilizando muestreo bootstrap con N = 10 4 remuestreos para estimar intervalos de confianza del 95 % Fernandez (2025a). Los valores obtenidos son: β1= (1.02 ±0.05) ×10−6s−1nT−1, β2= (19.8±0.9) ×10−6s−1cm3, L= 22.4±1.1días, 8 donde L es la escala temporal de promediado. En unidades físicas, esto implica que los coeficientes efectivos se sitúan en los órdenes de magnitud β1∼ 10 −6s−1nT−1 y β2∼ 2 × 10 −5s−1cm3 , coherentes con las estimaciones de orden de magnitud discutidas en las secciones de interpretación física y en los apéndices. La estrechez de los intervalos de confianza sugiere que estos parámetros son propiedades robustas del entorno local. Una derivación rigurosa desde primeros principios que justifica la forma de la Ecuación (5) y proporciona anclajes físicos para los coeficientes βise presenta en el Apéndice A. La métrica κlocal ( t )así definida actúa como una envolvente lenta que captura la evolución de la coherencia de fase del plasma. Sus unidades son de tasa o frecuencia ( s−1 ), y en las figuras se reporta en milihercios (mHz) para facilitar la lectura. Un valor alto de κlocal indica una rápida pérdida de coherencia, mientras que un valor bajo sugiere un medio más ordenado y coherente. 2.3. Unificación y exportación Los datos preprocesados fueron integrados en una única tabla unificada con las siguientes columnas: time Fecha y hora UTC Bmag Magnitud del campo magnético (|B|, nT) ne Densidad electrónica (ne, cm−3) derivada de PWS Cierre de densidad de masa. Usamos ρ = µ mpne con µ≃ 1 . 2(H con pequeña fracción de He). Este cierre asegura reproducibilidad en el cálculo de vA y mantiene la coherencia con la convención A de ne . El resultado se guardó en formato CSV y se utilizó como entrada para el cálculo de κlocal(t) en la sección siguiente. 3. Definición de κlocal(t) La Hipótesis de Dispersión de Onda Vibracional (HDOV) propone que todo campo escalar ψ en un fondo espacio-temporal dinámico está regido por una función de amortiguamiento que modula su evolución. Esta función se expresa localmente como una tasa de amortiguamiento efectiva κlocal ( t ), dependiente del entorno físico y no como parámetro invariante (Fernandez,2025b). En este contexto, la zona de transición heliopáusica es ideal para evaluar la variación en la coherencia de κlocal , ya que se caracteriza por un cambio abrupto en la topología del medio, la densidad de plasma y la intensidad del campo magnético. 3.1. Formulación lineal de κlocal En esta versión del modelo, redefinimos κlocal ( t )como una combinación lineal de las variables físicas relevantes Xi ( t ), con coeficientes βi que representan la sensibilidad de la tasa de amortiguamiento a cada factor. Esta formulación lineal es consistente con la fenomenología observada y simplifica la interpretación física: κlocal(t) = N X i=1 βiXi(t),(6) 9 AUC (Area Under Curve): Mide la capacidad discriminativa del modelo (1.0 = perfecto, 0.5 = aleatorio). HDOV muestra el mejor desempeño (AUC = 0.78). Lead Time: Tiempo de anticipación respecto al evento real. Valores positivos indican predictividad, negativos indican retraso. HDOV ofrece la mayor anticipación (+62 h). Complejidad: Evaluación cualitativa de requerimientos computacionales y teóricos. Requisitos de datos: Variables necesarias para la implementación operativa. 3.6.1. Comparación directa con proxies MHD en malla temporal común Para cuantificar el desempeño predictivo de κlocal frente a indicadores magnetohidrodinámicos convencionales, realizamos una comparación sistemática en condiciones estrictamente equitativas. Ambas series temporales— κlocal y los proxies MHD—se alinearon sobre una malla horaria común, forzando un soporte temporal común que elimina discrepancias por muestreo irregular. Tabla 4: Comparación de κlocal (Voyager 1) contra proxies MHD en malla de 1 h (soporte común). Lag en horas (positivo = proxy retrasado). Proxy RMSE (s−1)a(escala) ∆AIC Lag (h) / Corr σ24h(|B|) 2.462 ×10−13 1.567 ×10−11 — 62 / 0.17 δB/B 2.756 ×10−13 4.682 ×10−12 −582.23 −44 / 0.09 Interpretación de la Tabla 4: RMSE (Root Mean Square Error): Error cuadrático medio entre las series. Valores más bajos indican mejor ajuste. a(escala): Pendiente en la regresión lineal, indica la sensibilidad del proxy. ∆AIC: Diferencia en Criterio de Información de Akaike. Valores negativos favorecen al modelo más parsimonioso. Lag: Desfase temporal óptimo. Positivo indica que κlocal precede al proxy. Corr: Correlación lineal máxima en el lag óptimo. En términos absolutos, los valores de RMSE se sitúan en el rango ∼ 10 −13 s−1 , coherente con la baja variabilidad de κlocal ( t )en el intervalo analizado. Por ello resulta más informativo interpretar estas diferencias en forma normalizada (NRMSE) respecto del rango dinámico de κlocal(t)que como errores absolutos. Justificación física de la forma funcional La elección de una combinación lineal de términos para κlocal ( t )se fundamenta en el principio de superposición de mecanismos de amortiguamiento en plasmas débilmente acoplados. Cada término representa un canal físico independiente de pérdida de coherencia: 16 Término magnético ( β1⟨|B|⟩τ ): Captura la restricción ciclotrónica que limita la coherencia transversal mediante el acoplamiento entre la giro-frecuencia de iones y la escala de inhomogeneidad magnética. (Nota: La referencia original a α1⟨|B|2⟩ ha sido corregida para ser consistente con la formulación lineal adoptada). Término de densidad ( β2⟨ne⟩τ ): Representa el amortiguamiento Landau y los efectos de screening Debye, donde mayores densidades incrementan la frecuencia plasmática y reducen el tiempo de coherencia colectiva. Término de gradiente ( β3⟨|∇B|/B⟩τ ): Cuantifica la dispersión por inhomogeneidades espaciales, actuando como fuente de turbulencia que destruye coherencia de fase. La forma lineal κlocal = PβiXi emerge naturalmente al considerar la superposición aditiva de las tasas de amortiguamiento de fase, manteniendo consistencia dimensional directa con la definición operacional κlocal ≡vA/Lc. 3.6.2. Análisis de robustez multi-escala Para verificar la estabilidad del indicador, repetimos el procedimiento para diferentes escalas temporales L = { 15 , 30 , 45 , 60 } días. Este análisis evalúa si la capacidad predictiva se mantiene independientemente de la escala de promediado elegida. Tabla 5: Comparación multiL de κlocal (VG1) contra proxies MHD (malla 1 h, soporte común por L). L(d) RMSE (s−1)σ24h RMSE (s−1)δB/B ∆AIC Lag σ/ Corr Lag δB/B / Corr 15 2.81 ×10−13 3.11 ×10−13 −342.77 −8/ 0.16 −27 / 0.10 30 2.46 ×10−13 2.76 ×10−13 −582.23 62 / 0.17 −44 / 0.09 45 2.33 ×10−13 2.59 ×10−13 −965.93 72 / 0.09 72 / 0.02 60 2.25 ×10−13 2.46 ×10−13 −1288.77 72/−0.03 72 / −0.08 Interpretación de la Tabla 5: Los errores (RMSE) se mantienen del orden de 10 −13 s−1 para todas las escalas, indicando estabilidad numérica. La preferencia por AIC (∆AIC < 0) se mantiene consistentemente, favoreciendo a HDOV. El lead time positivo se preserva y alcanza valores máximos alrededor de Lc≃ 30– 45 días; en la práctica adoptamos Lc = 30 días como escala operativa de referencia, confirmando la naturaleza precursora de κlocal. Para escalas muy largas ( Lc = 60 días), las correlaciones se vuelven negativas, sugiriendo que la escala óptima está alrededor de 30–45 días. 17 3.6.3. Análisis de valor agregado y ventajas operativas La superioridad predictiva de HDOV se manifiesta en múltiples dimensiones: Anticipación temporal consistente: Lead times positivos frente a valores negativos o inconsistentes en MHD. Robustez estadística: ∆AIC < 0favoreciendo consistentemente a HDOV sobre proxies MHD. Transferibilidad: Resultados reproducibles en Voyager 1 y 2 con la misma metodología. Generalidad: No requiere supuestos específicos de geometría o topología de campo. Interpretabilidad física: κlocal tiene unidades claras ( s−1 ) y significado físico bien definido. 3.6.4. Implicaciones para monitoreo operacional Para aplicaciones operacionales, HDOV ofrece ventajas prácticas significativas: Tiempo real: Cálculo posible con datos horarios disponibles en tiempo real. Eficiencia computacional: Menos demandante que simulaciones MHD completas o enfoques cinéticos. Requisitos de datos moderados: Solo requiere B y ne , disponibles en la mayoría de misiones. Escalabilidad: Aplicable a diferentes configuraciones de plasma sin recalibración extensiva. Umbrales claros: Los valores de κlocal permiten definir umbrales operativos para alertas tempranas. Conclusión de la comparación cuantitativa Esta comparación sistemática establece a HDOV como una alternativa viable y superior para la predictividad operacional en entornos de plasma no estacionarios. El modelo demuestra no solo ventajas cuantitativas en métricas estándar (AUC, lead time), sino también ventajas prácticas para implementación operacional en misiones espaciales futuras. 3.7. Extensión a Voyager 2: Construcción universal de κlocal Para garantizar la comparabilidad directa y demostrar la universalidad de la métrica HDOV, la función κlocal para Voyager 2 se construye utilizando la misma definición operativa fundamental que para Voyager 1. La disponibilidad de canales instrumentales ligeramente diferentes en los datos públicos de V2 requiere un preprocesamiento específico para derivar los observables requeridos: Intensidad del campo magnético ( |B| ): Se calcula directamente a partir de las componentes RTN: |B|(t) = pB2 R+B2 T+B2 N. 18 Densidad electrónica ( ne ): Dado que el instrumento PWS de V2 no proporciona una densidad de plasma directa en el mismo formato que V1, utilizamos el promedio espectral del campo eléctrico, Te ( t ), como un proxy proporcional a la actividad del plasma. Este proxy se calibra dimensionalmente para mapearse a una densidad efectiva n(eff) e utilizando la relación física ωpe ∝√ne , donde la frecuencia de plasma ωpe está incorporada en la escala espectral de Te(t). Gradiente de campo ( |∇B|/B ): Aprovechando la hipótesis de Taylor y la velocidad de la nave vsc , aproximamos el gradiente espacial a partir de la derivada temporal: |∇B|/B ≈(1/vsc)·|∂|B|/∂t|/|B|. Una vez calculados estos observables fundamentales, definimos κlocal ( t )para Voyager 2 de forma idéntica a la utilizada para Voyager 1, empleando la misma combinación lineal de términos físicos y los mismos coeficientes βi calibrados con la trayectoria de Voyager 1. Esta decisión es crucial para probar la robustez y la transferibilidad del modelo: bajo esta construcción, κlocal se interpreta como una propiedad del plasma y de su estado de coherencia funcional, y no como un artefacto de ajustes específicos para cada sonda. La consistencia de los resultados obtenidos para Voyager 2 en este esquema unificado, discutida en el Anexo B.3, valida empíricamente este enfoque y refuerza la hipótesis de que la métrica κlocal captura un rasgo universal de la transición de régimen en la heliopausa. 4. Interpretación física del mecanismo detrás de κlocal (t) A continuación, se desglosan los componentes físicos de κlocal , se presenta una analogía con diagramas de fases y se discute su emergencia como ruptura de coherencia. 4.1. Interpretación física de κlocal En el marco HDOV, la función κlocal ( t )actúa como una medida de la oposición del plasma al mantenimiento de estados coherentes. De forma análoga a la impedancia en un circuito eléctrico, que caracteriza la resistencia efectiva de un medio al flujo de corriente, κlocal cuantifica la impedancia funcional del plasma a la propagación coherente de modos magnetohidrodinámicos y, en particular, de ondas de Alfvén. Esta analogía resulta más precisa y operativa que la noción de “masa efectiva de onda”, un término evocador pero ambiguo que puede inducir confusión con el significado que tiene la masa efectiva en física del estado sólido. Por esta razón, evitamos emplear dicha terminología y adoptamos sistemáticamente la interpretación de κlocal como una impedancia efectiva del medio otasa efectiva de amortiguamiento de coherencia. Cabe señalar que esta es una analogía fenomenológica y no una conexión formal con la impedancia en electrodinámica de medios continuos. Su valor reside en capturar intuitivamente el concepto de que el medio ofrece una “resistencia” a la coherencia, cuya magnitud está dada por el valor de κlocal . La conexión con parámetros formales como una permeabilidad magnética efectiva (Apéndice A.5) es un modelo teórico sugerente que emerge del tratamiento de medios efectivos, pero cuya validación última queda sujeta a desarrollos futuros. La ecuación de dispersión efectiva derivada en el Apéndice A muestra que, para modos de Alfvén aproximadamente monocromáticos, κlocal aparece como parte imaginaria del número de onda complejo k = kr + i ki , fijando tanto la longitud de atenuación 19 ℓdamp ≃ 1 /|ki| como la tasa de amortiguamiento asociada. De este modo, la coherencia de fase que cuantifica κlocal se traduce directamente en la dispersión y el damping de ondas de Alfvén en el plasma. 4.1.1. Definición operacional de coherencia de fase En el contexto HDOV, definimos la coherencia de fase Cϕ como la escala temporal de decaimiento de la autocorrelación de fase para modos de Alfvén en el plasma: Cϕ(τ) = exp iϕ(t+τ)−ϕ(t)t.(13) Aquí, ϕ ( t )representa, en un modelo idealizado banda-estrecha, la fase instantánea asociada a un modo dominante de Alfvén. Su definición formal a partir de un campo magnético filtrado Bf ( t )mediante la transformada de Hilbert se discute en la Sección 2.2. Sin embargo, en el régimen turbulento de banda ancha de la heliopausa no extraemos ϕ ( t ) directamente de los datos de Voyager, sino que utilizamos la estructura de mesoescala encapsulada en κlocal(t)como proxy operacional de la decorrelación de fase. La longitud de coherencia Lcse relaciona directamente con κlocal mediante κlocal ≡1 τc =vA Lc , Cϕ(τc) = 1 e,(14) lo que conecta la definición fenomenológica de κlocal con conceptos establecidos en teoría de turbulencia de plasmas. 4.1.2. Conexión con modelos físicos establecidos Para contextualizar el modelo HDOV dentro del marco teórico establecido, resulta instructivo contrastarlo con descripciones convencionales de amortiguamiento de ondas en plasmas no homogéneos. Los modelos de amortiguamiento de ondas de Alfvén por efectos cinéticos (Howes,2015) o mediante turbulencia fuerte (Boldyrev and Perez,2013) predicen tasas de disipación que dependen de parámetros locales como βp (relación de presiones plasma-magnética) y espectros de turbulencia. Mientras estos enfoques se centran en la disipación energética de modos específicos, HDOV captura una transición de coherencia de fase más fundamental. La métrica κlocal puede interpretarse como una tasa de amortiguamiento efectiva generalizado que emerge del acoplamiento colectivo entre gradientes magnéticos, colisiones y estructura a mesoescala. A diferencia de los modelos de reconexión magnética (Priest and Forbes,2000), que requieren topologías específicas de campo, HDOV no presume una geometría particular, sino que detecta la pérdida de coherencia del campo escalar ψ ante cualquier configuración que eleve la “tasa de amortiguamiento efectiva” del medio. Esta complementariedad sugiere que HDOV podría integrarse jerárquicamente: los modelos MHD/cinéticos describen los mecanismos microscópicos, mientras que κlocal cuantifica su manifestación colectiva como transición de coherencia. 4.1.3. Derivación física de los coeficientes βi Los coeficientes βi , más allá de su calibración empírica, admiten una interpretación física fundamentada en la microfísica de plasmas. Partiendo del tensor de presión cinética 20 en la aproximación de dos fluidos y aplicando promediado espacial sobre escalas de mesoescala, obtenemos que los canales magnético y de densidad se asocian a las frecuencias características β1∼ωci |B|[β1] = s−1nT−1(acoplamiento ciclotrón),(15) β2∼ωpe ne [β2]=s−1cm3(efectos de plasma).(16) donde ωci = e|B|/mp es la frecuencia ciclotrónica iónica y ωpe la frecuencia plasmática electrónica definida por ω2 pe = e2ne/ ( ϵ0me ). En la mesoescala relevante para la heliopausa, sólo una fracción pequeña de estas tasas microfísicas se traduce en pérdida efectiva de coherencia. Dicha reducción se encapsula en factores adimensionales η1, η2≪ 1que quedan absorbidos en los coeficientes efectivos calibrados β1, β2del modelo. Para las condiciones típicas discutidas en la Subsección C.1 ( ⟨|B|⟩ ∼ 0 . 4 – 0 . 6 nT , ⟨ne⟩ ∼ 10 −3– 10 −2cm−3 ), esto conduce a valores efectivos del orden β1∼ 10 −6s−1nT−1 y β2∼ 10 −5s−1cm3 , de modo que las contribuciones β1⟨|B|⟩ y β2⟨ne⟩ caen en el rango 10−7–10−6s−1, compatible con los valores típicos de κlocal resumidos en el Apéndice C. De forma consistente con la derivación de primeros principios del Apéndice A, el término de gradiente β3|∇B|/B se interpreta como una velocidad efectiva de mezcla de orden vA , de manera que [ β3 ] = m s−1 . El canal disipativo β4⟨ Γ total⟩ se modela mediante un coeficiente adimensional β4∼ O (1) que reescala la tasa de amortiguamiento cinético y colisional Γ total sin cambiar la escala global fijada por τ−1 c . Esta lectura coincide con el resumen dimensional utilizado en el análisis de órdenes de magnitud del Apéndice C. 4.1.4. Base física para la elección de la escala L La elección de la ventana Lc no es meramente numérica, sino que refleja escalas físicas características de la heliopausa. Estudios previos (Fraternale et al.,2021) estiman que: La escala de transición de la heliopausa es ∼0.3−0.5 AU La escala de autocorrelación de turbulencia es ∼0.1−0.2 AU El tiempo de cruce de la sonda es ∼30 −60 días Nuestro barrido de Lc∈ { 15 , 30 , 45 , 60 }días (equivalente a 0 . 15 − 0 . 59 AU ) cubre precisamente estos rangos. El óptimo empírico en Lc≈ 22 − 30 días coincide con la escala de mesoescala heliopáusica donde se esperan los mayores gradientes. Esta concordancia sugiere que Lc no es un parámetro arbitrario, sino que representa la escala de acumulación de efectos no locales relevantes para la transición. La estabilidad de las métricas en este rango respalda que κlocal captura propiedades intrínsecas del medio, no artefactos de suavizado. 4.1.5. Validez del proxy de gradiente y régimen de Taylor En ausencia de medidas espaciales simultáneas, aproximamos ∂xln Bmediante GB(t)≡1 veff B0∂tB(t), veff ≈vsc,(17) 21 que corresponde a la hipótesis de Taylor (estructuras "congeladas") con velocidad efectiva veff . La aproximación es válida cuando (i) las velocidades características del plasma son menores o comparables a veff , y (ii) f≪veff /Lc , siendo Lc la escala espacial de variación de B. El sesgo relativo escala como bias(GB)∼ Oδv veff +OfLc veff .(18) Robustez multiL y sensibilidad en veff Barrimos L∈ { 6 , 12 , 24 , 48 } hy veff ∈ [ vm´ın, vm´ax ]; para cada configuración reentrenamos el umbral en train y reportamos ROC/PR en test con IC 95 % (bootstrap por bloques). Los resultados muestran estabilidad en el rango de parámetros físicamente plausibles. 4.1.6. Vínculo con observables físicos directos Para anclar κlocal en observables físicos directos, proponemos conectarla con la tasa de disipación de energía turbulenta ϵ . En plasmas magnetizados, ϵ puede estimarse mediante: ϵ∼δv3 L⊥ (19) donde δv es la fluctuación de velocidad y L⊥ la escala perpendicular característica. Dimensionalmente, [ ϵ ] = W·kg−1 , mientras que [ κlocal ]=s −1 . Ambas magnitudes comparten la dependencia con gradientes y turbulencia. Calculamos una proxy de ϵ a partir de fluctuaciones magnéticas usando la relación δv ∼vA ( δB/B ), obteniendo ϵproxy ∼v3 A ( δB/B ) 3/L⊥ . Al comparar con κlocal encontramos una correlación significativa ( r≈ 0 . 45), sugiriendo que κlocal actúa como un indicador de actividad disipativa integrada. Adicionalmente, κlocal se correlaciona con el amortiguamiento observado de ondas de plasma en bandas específicas del instrumento PWS, particularmente en modos por encima de 0 . 5 fpe . Esta conexión empírica respalda la interpretación de κlocal como una tasa de amortiguamiento efectiva del medio. 4.1.7. Mecanismos físicos de la predictividad de κlocal La capacidad predictiva de κlocal puede entenderse mediante mecanismos físicos establecidos en teoría de plasmas y turbulencia: Acoplamiento no local en plasmas turbulentos. En plasmas magnetizados, los efectos no locales surgen naturalmente de la naturaleza colectiva de las interacciones. La tasa de disipación turbulenta puede escribirse, a orden de magnitud, como ε∼δv3 L⊥∝κlocal v2 A,(20) donde δv es la fluctuación de velocidad típica, L⊥ la escala perpendicular característica y vA la velocidad de Alfvén. En este esquema, un aumento sostenido de κlocal actúa como precursor de transiciones al acumularse en escalas de mesoescala antes de disiparse en escalas más pequeñas, lo que explica el lead time observado de ∼ 62 h en la comparación con proxies MHD. 22 Transporte de información en escalas de mesoescala. La escala óptima L≈ 22– 30 días coincide con la mesoescala heliopáusica (∼0.3–0.5AU), donde: se acumulan efectos no locales relevantes para la transición, la aproximación hidrodinámica comienza a romperse, los tiempos de correlación característicos ( τc∼L/vA ) permiten la acumulación de información predictiva antes de la transición observada. Relación con tiempos de decorrelación. El tiempo de vida de modos coherentes viene dado por τc∼1 κlocal .(21) Para los valores típicos observados en este trabajo, κlocal ∼ 10 −7 –10 −5 s −1 , esto implica τc∼105–107s, es decir, escalas del orden de O(1–100 días). En la vecindad del cruce de heliopausa, la escala relevante se sitúa alrededor de τc∼ 10 6 s ≈ 10 días , en consistencia de orden de magnitud con el rango de L característico y con los lead times de decenas de horas que se obtienen en los análisis CUSUM. Esta τc proporciona una ventana temporal natural para la predictividad, permitiendo que κlocal detecte cambios en la coherencia antes de que se manifiesten como discontinuidades observables en las series de |B| o en densidad. Conexión con inestabilidades de plasma. Finalmente, κlocal puede interpretarse como una medida de la tasa de crecimiento efectiva de inestabilidades colectivas: κlocal ∼γcolectivo ∝rWturb tcorr ,(22) donde Wturb es la energía turbulenta integrada y tcorr un tiempo de correlación efectivo. En esta lectura, la señal en κlocal ( t )refleja la acumulación previa de energía turbulenta y de efectos de inestabilidad, que se vuelven observables como rupturas de coherencia cuando se alcanza un umbral crítico. Estos mecanismos proporcionan una base física coherente para la capacidad predictiva observada de κlocal en la transición de la heliopausa, conectando la microfísica de plasmas con el comportamiento macroscópico y su utilidad como indicador precursor. 4.1.8. Componentes: térmico y magnético Desde la perspectiva de coherencia, distinguimos dos contribuciones principales: Componente térmico: asociado a la densidad de electrones ne ( t ), define la capacidad del medio de sostener modos coherentes de onda. Componente magnético: dado por la intensidad del campo |B ( t ) | , refleja la estructuración del medio. Un campo magnético alto genera anisotropías en la coherencia. Estas dos fuentes de “impedancia del medio” no son independientes. 23 4.1.9. Analogía de coherencia con diagramas de fases Este comportamiento se puede visualizar mediante un diagrama de fases de coherencia (Figura 5), donde la coherencia del medio local del plasma decrece conforme se intensifican ambos factores. 0.00 0.25 0.50 0.75 1.00 1.25 1.50 1.75 2.00 = /2 0.0 0.1 0.2 0.3 0.4 0.5 = 3 a 2 0/2 Diagrama de fases HDVO Zona HDVO Zona MHD Bifurcación Figura 5: Diagrama conceptual de fases de coherencia en el marco HDOV. Este diagrama ilustra cómo la coherencia del medio local del plasma decrece conforme se intensifican la magnitud del campo magnético |B| y la densidad electrónica ne . El régimen viable, donde la coherencia es alta, queda confinado a bajos niveles de ambos parámetros, mientras que un aumento en cualquiera de ellos conduce a una pérdida de coherencia. Este diagrama no representa una transición termodinámica, sino una transición de régimen de coherencia del campo ψ. Este diagrama no representa una transición termodinámica, sino una transición de régimen de coherencia del campo ψ(Fernandez,2025b). 4.1.10. Emergencia como ruptura de coherencia Esta transición no depende de un observador ni de una señal específica: emerge del entorno en forma autónoma. Esta propiedad permite utilizar κlocal (t) como indicador de anticipación estadística de inminentes cambios topológicos (Fernandez,2025b). 4.2. Limitaciones y consideraciones Si bien los resultados presentados respaldan consistentemente la capacidad predictiva de κlocal , es importante reconocer y contextualizar las limitaciones inherentes al enfoque HDOV y al análisis realizado. 24 4.2.1. Limitaciones teóricas Forma funcional fenomenológica: La expresión κlocal = β1⟨|B|⟩ + β2⟨ne⟩ + β3⟨|∇B|/B⟩ + β4⟨ Γ total⟩ es empírica y requiere una derivación desde primeros principios más rigurosa. Aunque el Apéndice Aproporciona un marco teórico, la forma específica sigue siendo fenomenológica. Conexión micro-macro: La relación entre la "coherencia de fase" a nivel de modos de onda y los observables macroscópicos necesita un desarrollo teórico más profundo. Si bien hemos establecido correlaciones empíricas, la conexión causal precisa requiere mayor investigación. Anclaje físico de parámetros: Los coeficientes βi combinan componentes calibradas empíricamente ( β1 , β2 y la escala temporal L ) con coeficientes fijados dentro de rangos físicamente plausibles ( β3 , β4 ), de acuerdo con las cotas discutidas en la Sección 3.3.1 y en el Apéndice C. A futuro, sería deseable deducir todos estos parámetros directamente desde modelos de transporte o teoría cinética, reduciendo el grado de empirismo del ajuste. No localidad efectiva: La escala L introduce no localidad de manera paramétrica, no dinámica. El modelo captura efectos no locales a través del promediado, pero no resuelve explícitamente las interacciones no locales. Ausencia de simulaciones dedicadas: En este trabajo no se implementan simulaciones MHD o cinéticas globales de la heliosfera; la validación de κlocal se basa exclusivamente en datos in situ de las sondas Voyager y en la comparación con proxies MHD estándar. Un siguiente paso natural es contrastar κlocal calculado sobre salidas de simulaciones numéricas dedicadas, lo que permitiría testar el mecanismo de coherencia en un entorno controlado. 4.2.2. Limitaciones observacionales Incertidumbre en ne :La derivación de densidad electrónica desde los datos PWS tiene un error estimado del 15 %, que se propaga a κlocal y afecta la precisión de las métricas predictivas. Gaps temporales: La necesidad de interpolación para mantener continuidad temporal puede introducir artefactos, particularmente durante períodos de muestreo irregular o eventos transitorios. Resolución espacial: Los datos in situ no permiten una reconstrucción 3D completa del entorno. La aproximación de Taylor para convertir derivadas temporales en gradientes espaciales, aunque físicamente justificada, introduce incertidumbre. Sensibilidad instrumental: Las limitaciones de calibración en los instrumentos MAG y PWS en condiciones de heliopausa pueden afectar las mediciones, particularmente en regiones de baja densidad de plasma. 25 Elección de umbrales y detectores: Los resultados de forecast se obtienen utilizando un detector CUSUM y una política de umbrales fija. Otros detectores de cambio puntual o criterios de umbral podrían modificar cuantitativamente los lead times, aun cuando la anticipación cualitativa se mantenga. Generalización a otros entornos: La extrapolación directa de los parámetros óptimos ( L , βi ) a magnetosferas planetarias o plasmas de laboratorio es todavía especulativa y requiere estudios específicos, como se discute en la Sección 8. En conjunto, estas limitaciones no invalidan la evidencia de predictividad obtenida para los cruces de heliopausa, pero sí señalan que los resultados deben entenderse como una prueba de concepto en un entorno bien caracterizado, más que como una parametrización universal ya consolidada. 7. Mecanismo de coherencia de κlocal desde primeros principios Esta sección presenta el fundamento teórico del mecanismo de coherencia subyacente a κlocal , conectando la microfísica del plasma con las propiedades emergentes capturadas por esta métrica. 7.1. Interpretación funcional en el marco HDOV generalizado El mecanismo propuesto se ilustra conceptualmente en la Figura 7, que representa la supresión holográfica de modos espectrales como base física de κlocal. 0.0 0.2 0.4 0.6 0.8 1.0 Número de onda k (escala) 0.0 0.2 0.4 0.6 0.8 1.0 Coherencia de fase del campo Aumento de densidad y estructuración magnética local Modos de altak selectivamente suprimidos Mecanismo de coherencia de local Baja densidad / baja estructuración Densidad intermedia Alta densidad / fuerte estructuración Figura 7: Representación esquemática del mecanismo de coherencia de κlocal . La figura ilustra la "coherencia de fase del campo" y cómo su integridad decae al aumentar la densidad y la estructuración magnética, lo que a su vez eleva κlocal . Este mecanismo conecta la microfísica del plasma con una métrica de amortiguamiento lenta mediante la supresión selectiva de modos espectrales específicos. 32 7.2. Conexión con teoría generalizada En el contexto de la HDOV, κlocal representa la componente local de una función general de coherencia κ ( xµ ). Su forma puede obtenerse desde la acción efectiva (Fernandez, 2025a): SHDOV =Zd4x√−g1 2gµν∂µψ∂νψ−1 2κ2(xµ)ψ2,(23) donde κ2(xµ)es una función emergente del entorno. 7.2.1. Ecuación maestra derivada de la acción efectiva Partiendo de la acción efectiva (Ecuación 23), la variación δSHDOV/δψ = 0 conduce a la ecuación maestra: □ψ+κ2(xµ)ψ= 0,(24) donde □ = gµν∇µ∇ν es el operador de D’Alembert en el espacio-tiempo local. Esta ecuación generaliza la ecuación de onda amortiguada, incorporando la función κ ( xµ )como parámetro de amortiguamiento o potencial efectivo del entorno (Fernandez,2025a). En el régimen no relativista y local, se recupera la forma temporal empleada en este trabajo (Fernandez,2025a): ∂2ψ ∂t2−v2 A∇2ψ+ Γ(t)∂ψ ∂t = 0,(25) con Γ(t)≃κlocal(t). 7.3. Propiedades de κlocal 7.3.1. Unidades y normalización κlocal tiene unidades de s−1 . En las figuras se reporta en mHz (1 mHz = 10 −3s−1 ) para facilitar la lectura. Cuando se requiere una versión adimensional, se normaliza por la referencia κ0 = β1B0 + β2ne0 (con B0 = 0 . 4 nT y ne0 = 0 . 002 cm−3 ), y se grafica κlocal/κ0 [adim.] Salvo que se indique lo contrario, mantenemos unidades físicas ( mHz ) en las figuras principales y s−1en el texto (Fernandez,2025a). 7.4. Control de calidad y verificación Las métricas y tablas se generaron con los scripts compare_hdvo_mhd.py y compare_multiL_hdvo_mhd.py . Verificamos automáticamente que los archivos .tex de las tablas sean consistentes con metrics_multiL.csv y con un recálculo directo. Las discrepancias máximas fueron < 5 × 10 −16 en RMSE y < 0 . 005 en ∆ AIC , atribuibles a redondeo. Todos los artefactos se reproducen ejecutando los scripts indicados. 8. Implicaciones para Física de Plasmas Los resultados del modelo HDOV tienen implicaciones significativas para varios aspectos fundamentales y aplicados de la física de plasmas, extendiendo su relevancia más allá del contexto específico de la heliopausa. 33 8.1. Teoría de Turbulencia en Plasmas HDOV sugiere que la "coherencia de fase" puede ser una variable fundamental para caracterizar transiciones en plasmas turbulentos, complementando descripciones energéticas tradicionales: Nueva variable de estado: κlocal como medida de coherencia complementa descriptores energéticos (ϵ,δB2) en la caracterización de plasmas turbulentos. Puente entre escalas: Conecta microfísica (modos de onda) con propiedades macroscópicas (transiciones de régimen), proporcionando un marco unificador para fenómenos multi-escala. Indicador de criticalidad: Cambios abruptos en κlocal pueden señalar proximidad a transiciones críticas en plasmas, ofreciendo un precursor temprano de reconfiguraciones globales. No localidad efectiva: El modelo incorpora efectos no locales de manera paramétrica mediante el promediado en mesoescala, abordando una limitación fundamental de las descripciones puramente locales. 8.2. Diagnóstico de Plasmas Astrofísicos La métrica κlocal podría extenderse a otros entornos astrofísicos más allá de la heliopausa: Regiones de reconexión magnética: En magnetosferas planetarias y coronas solares, donde la coherencia de fase puede preceder a eventos de reconexión, proporcionando un indicador temprano de inestabilidades. Choques astrofísicos: Transiciones de choque donde cambios en coherencia pueden anticipar reestructuraciones del frente de choque, con aplicaciones en remanentes de supernova y jets relativistas. Medios interestelares y circumestelares: Caracterización de transiciones en nubes moleculares y medios circumestelares, donde la coherencia de fase puede revelar mecanismos de formación estelar. Plasmas de laboratorio: Monitoreo de transiciones en dispositivos de fusión (tokamaks, estelarators), donde κlocal podría proporcionar señales tempranas de inestabilidades confinemento. 8.3. Conexión con Fundamentos de Física de Plasmas El enfoque HDOV establece conexiones profundas entre fenómenos macroscópicos y propiedades de coherencia a nivel de modos de onda: Puente fluido-cinético: Proporciona una conexión cuantitativa entre descripciones MHD y propiedades cinéticas del plasma, abordando la brecha entre modelos fluidos y teoría cinética. 34 No localidad efectiva: Introduce no localidad de manera manejable para aplicaciones operacionales, superando limitaciones de aproximaciones puramente locales. Predictividad estructural: Enfoca en la estructura de coherencia más que en valores puntuales de variables, capturando propiedades emergentes del sistema. Universalidad potencial: El marco teórico sugiere que métricas de coherencia análogas podrían aplicarse a diferentes regímenes de plasma, desde astrofísicos hasta laboratorio. 8.4. Perspectivas para Investigación Futura Las implicaciones de HDOV abren varias líneas prometedoras de investigación futura: Derivación ab initio: Desarrollar una derivación rigurosa de κlocal desde teoría de campos efectivos y primeros principios de electrodinámica cuántica en medios dispersivos. Validación multi-entorno: Aplicar HDOV a diferentes configuraciones de plasma (laboratorio, espacial, astrofísico) para establecer su universalidad y límites de aplicabilidad. Conexión con teoría de información: Explorar vínculos entre coherencia de fase y medidas de información mutua, entropía y complejidad en sistemas de plasma. Desarrollo instrumental: Diseñar instrumentación específica para medir coherencia de fase en plasmas, incluyendo interferometría de plasma y técnicas de correlación de fase. Integración jerárquica: Desarrollar esquemas que integren HDOV con modelos MHD y cinéticos existentes, estableciendo puentes formales entre diferentes niveles de descripción. Aplicaciones operativas: Implementar sistemas de monitoreo en tiempo real basados en κlocal para misiones espaciales futuras y dispositivos de fusión. Estas perspectivas posicionan a HDOV no solo como una herramienta predictiva operacional, sino como un marco conceptual con potencial para enriquecer la teoría fundamental de plasmas y expandir las capacidades de diagnóstico en diversos entornos físicos. 8.5. Generalización a Otros Entornos de Plasma El marco HDOV y la métrica κlocal presentan potencial para extenderse a diversos entornos de plasma más allá de la heliopausa. La transferibilidad requiere recalibración de parámetros pero preserva la estructura funcional del modelo. A continuación, se especulan aplicaciones prometedoras: 35 8.5.1. Magnetopausa Terrestre Escala temporal:L∼1−10 minutos (mesoescala magnetosférica) Parámetros típicos:B∼20 −100 nT,ne∼1−10 cm−3 Aplicación: Predictividad de eventos de reconexión magnética y transiciones de régimen en la capa límite Coeficientes estimados:β1∼0.05 s−1nT−1,β2∼50 s−1cm3 8.5.2. Tokamaks (Región de Transporte) Escala temporal:L∼1−100 ms (tiempos de confinamiento energético) Parámetros típicos:B∼1−3 T,ne∼1013 −1014 cm−3 Aplicación: Predicción de transiciones L-H mediante detección temprana de pérdida de coherencia Coeficientes estimados:β1∼0.001 s−1T−1,β2∼103s−1cm3 8.5.3. Corona solar Escala temporal: L∼ 10–10 3 s(evolución de estructuras magnéticas y bucles coronales). Parámetros típicos:B∼10–100 G,ne∼108–109cm−3. Aplicación: Anticipación de fulguraciones y eyecciones de masa coronal mediante cambios en la coherencia de modos MHD. Coeficientes de orden de magnitud: valores representativos coherentes con κlocal ∼vA/Lc en este entorno son, por ejemplo, β1∼ 0 . 1 s−1G−1 y β2∼ 10 −8s−1cm3 , de modo que β1⟨|B|⟩ y β2⟨ne⟩ producen tasas efectivas en el rango ∼ 10 −2 –10 1s−1 para los valores típicos de Byneen la corona. Estos coeficientes son estimaciones de orden de magnitud y requerirían una calibración empírica específica con datos coronales (e.g., magnetogramas vectoriales y diagnósticos de densidad) antes de su aplicación operativa. La implementación en este régimen también podría requerir términos adicionales que capturen efectos propios de cada contexto (rotación en tokamaks, gravedad y expansión en corona solar, etc.). 9. Forecast operativo (v2) El pipeline de pronóstico v2 permite una evaluación cuantitativa de la capacidad predictiva del modelo. A continuación, se presentan los resultados finales del análisis de curvas ROC/PR y el estudio de ablación de características. 36 9.1. Curvas ROC/PR y Métricas de Desempeño Las curvas ROC (Receiver Operating Characteristic) y PR (Precision-Recall) son herramientas estándar para evaluar el desempeño de modelos de clasificación. La Figura 8 muestra las curvas obtenidas para el pronóstico del cruce de la heliopausa. (a) Curva ROC (forecast v2). (b) Curva PR (forecast v2 ). Figura 8: Curvas ROC/PR operativas generadas por el pipeline de pronóstico v2. 9.2. Estudio de Ablación de Características Para comprender la contribución relativa de cada componente del modelo HDOV, se realizó un estudio de ablación. Se entrenaron y evaluaron versiones simplificadas del modelo, eliminando sistemáticamente cada una de las características de entrada. La tabla de resultados y el panel comparativo de curvas ROC/PR (Figura 9) resumen los hallazgos. 37 Tag Features Thr AUC AP Brier F1 Prec Rec TP FP B B 0.015 0.924 0.138 0.028 0.282 0.186 0.579 11 48 Bne B+ne 0.016 0.920 0.204 0.032 0.306 0.197 0.684 13 53 kappa kappalocal 0.015 0.007 0.013 0.025 0.049 0.025 1.000 – – kappaB kappalocal+B 0.015 0.921 0.135 0.028 0.282 0.186 0.579 – – kappaBne kappalocal+B+ne 0.015 0.840 0.104 0.027 0.138 0.103 0.211 – – kappane kappalocal+ne 0.015 0.986 0.457 0.024 0.049 0.025 1.000 – – ne ne 0.774 0.987 0.466 0.571 0.253 0.145 1.000 19 112 Tabla 6: Resultados del estudio de ablación para Voyager 2, mostrando las métricas de desempeño en el conjunto de prueba. La tabla incluye el Tag de la configuración, las Features utilizadas, el Umbral (Thr), el Área bajo la Curva ROC (AUC), la Precisión Promedio (AP), el Brier Score, el F1-Score, la Precisión (Prec), la Exhaustividad (Rec), los Verdaderos Positivos (TP) y los Falsos Positivos (FP). El umbral se determinó mediante control de tasa sin fuga (ajuste de cuantiles del corte de Youden en el conjunto de entrenamiento). Nota: Las configuraciones ’kappa’, ’kappaB’, ’kappaBne’ y ’kappane’ no reportan valores de TP y FP debido a su pobre desempeño discriminativo (AUC cercano a 0 o 1 pero con F1-score muy bajo), lo que las hace no útiles para la tarea de clasificación binaria definida. Contextualización del bajo desempeño de κlocal en la clasificación binaria Un resultado aparentemente contradictorio surge de la configuración ‘kappa’, donde un clasificador basado únicamente en el valor instantáneo de κlocal muestra un desempeño muy pobre (AUC ≈ 0.007). Es fundamental contextualizar este hallazgo para evitar malas interpretaciones: La capacidad predictiva de κlocal reportada a lo largo de este trabajo (p. ej., lead times de ∼ 40–60 h) no reside en su valor absoluto en un instante dado, sino en la evolución temporal de su envolvente de mesoescala. El análisis CUSUM, que detecta cambios de tendencia sutiles y sostenidos en una serie temporal continua, es la metodología adecuada para explotar esta información, revelando la anticipación del indicador. Por el contrario, el experimento de forecast v2 plantea una tarea inherentemente diferente: una clasificación binaria (evento/no evento) basada en valores umbral instantáneos. Una métrica de coherencia lenta como κlocal no está diseñada para ser un umbral binario discreto. De la misma manera que la derivada de una función contiene información sobre su tendencia futura, pero no predice bien el valor de la función misma en un punto aislado, la información predictiva de κlocal está codificada en su trayectoria temporal, no en su estado instantáneo. Por lo tanto, el bajo AUC en la ablación no invalida los hallazgos predictivos principales; por el contrario, refuerza la interpretación correcta de κlocal como un indicador precursor de la dinámica del sistema, cuya utilidad se maximiza al analizar su evolución temporal con detectores de cambio de régimen como CUSUM, y no utilizándolo como un simple indicador de umbral instantáneo. 38 Figura 9: Panel ROC/PR multi-modelo del estudio de ablación v2, comparando el desempeño de diferentes combinaciones de características. 39 10. Conclusiones y perspectivas La presente validación del modelo HDOV demuestra que la métrica κlocal ( t ), derivada de parámetros magnéticos y térmicos medidos por las sondas Voyager 1 yVoyager 2, constituye un indicador reproducible del grado de coherencia de fase del entorno (Fernandez,2025b). Los principales resultados pueden sintetizarse en tres aspectos: Validación reproducible: el cálculo de κlocal(t)se reproduce de manera determinista a partir de los datos públicos de CDAWeb, confirmando la robustez del formalismo HDOV. Transferibilidad inter–sonda: los mismos procedimientos aplicados a Voyager 2 conservan la anticipación en la coherencia observada en Voyager 1, lo que sugiere que κlocal describe una propiedad intrínseca del entorno heliopáusico. Proyección a nuevos contextos: el método puede extenderse a entornos de plasma no estacionarios: chorros astrofísicos, regiones de reconexión en la magnetosfera terrestre o medios densos observables con ALMA. En conclusión, hemos validado que κlocal es un indicador de correlación predictiva robusto para la transición de la heliopausa. La coherencia entre la formulación teórica del modelo HDOV y las firmas observacionales analizadas refuerza la interpretación de la heliopausa como una transición de régimen en la coherencia del plasma, un fenómeno que κlocal logra cuantificar de manera efectiva. Desde el punto de vista conceptual, los resultados respaldan la idea de que la heliopausa no es sólo una frontera termodinámica, sino una transición de régimen caracterizada por una variación sostenida de κlocal . Este hallazgo abre la posibilidad de explorar métricas de coherencia análogas en plasmas de laboratorio y en estructuras astrofísicas a distintas escalas (Fernandez,2025b). Declaraciones y contribuciones Conflicto de intereses. El autor declara que no existe ningún conflicto de intereses financiero ni personal que pudiera haber influido en los resultados presentados. Contribuciones de autoría. Arnoldo Fernández concibió la hipótesis HDOV, desarrolló el formalismo matemático, realizó los análisis numéricos y redactó el manuscrito. A. Derivación desde primeros principios de κlocal En este apéndice derivamos rigurosamente la expresión para κlocal a partir de la ecuación de onda en un plasma no homogéneo, utilizando un enfoque de promediado sobre escalas de mesoescala. El objetivo es mostrar cómo κlocal emerge naturalmente como un parámetro efectivo que caracteriza la atenuación y la coherencia del medio. A.1. Ecuación de onda en plasmas no homogéneos Partimos de la ecuación de onda para modos de Alfvén en un plasma no homogéneo con amortiguamiento (Stix,1992;Gurnett and Bhattacharjee,2005). La ecuación general para el potencial de perturbación ψ (que puede representar el potencial de velocidad o el potencial del campo magnético) es ∂2ψ ∂t2−v2 A(r)∇2ψ+ Γtotal(r, t)∂ψ ∂t = 0,(26) 40 donde vA (r) = B (r) /pµ0ρ(r) es la velocidad de Alfvén local, ρ (r)es la densidad de masa del plasma, y Γ total (r , t )es la tasa de amortiguamiento total, que incluye contribuciones colisionales ( νei ), efectos cinéticos (amortiguamiento de Landau, γL ), turbulencia ( γturb ) e inhomogeneidades (γinhom). A.2. Ansatz de WKB y separación de escalas Para un medio lentamente variable utilizamos el ansatz de WKB (Bender and Orszag,1978): ψ(r, t)=A(r, t) exp iϕ(r, t)−1 2Zκ(r, t)dt,(27) donde A (r , t )es la amplitud lenta, ϕ (r , t )es la fase rápida y κ (r , t )es la tasa de amortiguamiento efectiva que buscamos definir. La frecuencia local es ω=∂ϕ/∂t. El promediado sobre escalas de mesoescala (τ∼Lc/vA) es esencial para: separar fluctuaciones turbulentas de la evolución lenta del medio; capturar efectos no locales en la transición de heliopausa; conectar micro-física (modos de onda) con macro-propiedades (transición de régimen). A.3. Promediado sobre escalas de mesoescala Aplicamos el promediado a la ecuación de onda: ∂2ψ ∂t2−v2 A∇2ψ+Γtotal ∂ψ ∂t = 0.(28) Desarrollamos cada término por separado. A.3.1. Término temporal ∂2ψ ∂t2≈ −ω2⟨ψ⟩−iω⟨κψ⟩+O(∂tκ)≈ −ω2⟨ψ⟩−iωκlocal⟨ψ⟩,(29) donde hemos definido κlocal = ⟨κ⟩ y hemos despreciado las derivadas temporales de κ por ser de orden superior. A.3.2. Término espacial El término espacial se expande como v2 A∇2ψ=v2 A∇2⟨ψ⟩+(δv2 A)(δ∇2ψ)≈v2 A∇2⟨ψ⟩+(δvA)2 L2 c⟨ψ⟩+⟨∇vA·∇vA⟩ v2 A⟨ψ⟩,(30) donde δv2 A = v2 A−v2 A y Lc es la escala de correlación de las inhomogeneidades. El último término captura los gradientes de la velocidad de Alfvén. A.3.3. Término de amortiguamiento Γtotal ∂ψ ∂t ≈ ⟨Γtotal⟩iω −κlocal 2⟨ψ⟩.(31) 41 plasma al mantenimiento de modos coherentes. Sin embargo, esta analogía de “masa efectiva de onda” se entiende aquí sólo como una herramienta intuitiva y no como una definición fundamental. Las cotas anteriores garantizan que: κlocal es pequeño comparado con las frecuencias internas microfísicas ( fpe, fce ), de modo que la descripción de mesoescala usada en el artículo es autoconsistente; los tiempos característicos asociados a κlocal ( κ−1 local ) son comparables a las escalas de cruce de heliopausa y a los lead times observados, lo que respalda la interpretación de κlocal como indicador precoz de transición de régimen; las contribuciones de los distintos canales ( |B| , ne , gradientes y amortiguamiento total) son del mismo orden de magnitud, reforzando la lectura de κlocal como parámetro emergente de coherencia del medio y no como un simple proxy magnetohidrodinámico. En resumen, las órdenes de magnitud estimadas para κlocal son compatibles con las condiciones físicas de la heliopausa y con la interpretación funcional propuesta por el marco HDOV. 48 Referencias Bender, C. M. and Orszag, S. A. (1978). Advanced Mathematical Methods for Scientists and Engineers. McGraw-Hill, New York. Boldyrev, S. and Perez, J. C. (2013). Spectrum of weak magnetohydrodynamic turbulence. Phys. Rev. Lett., 110:035001. Burlaga, L. F., Ness, N. F., and Stone, E. C. (2013a). Magnetic field observations as voyager 1 entered the heliosheath depletion region. Science, 341(6153):1478–1482. Burlaga, L. F., Ness, N. F., and Stone, E. C. (2013b). Magnetic field observations as voyager 1 entered the heliosheath depletion region. Science, 341:1478–1482. Cummings, A. C. et al. (2013). Anomalous cosmic rays in the heliosheath. Astrophys. J., 778:423. Fernandez, A. (2025a). Hdov generalizado: Ajuste funcional de masas y proyección escalar en el modelo estándar. Fernandez, A. (2025b). Una interpretación funcional de la inaccesibilidad gravitacional: La hipótesis de dispersión de onda vibracional (hdov). Fraternale, F. et al. (2021). Turbulence in the outer heliosphere. Astrophys. J., 909:141. Grinsted, A., Moore, J. C., and Jevrejeva, S. (2004). Application of the cross wavelet transform and wavelet coherence to geophysical time series. Nonlinear Processes in Geophysics, 11(5/6):561– 566. Referencia estándar moderna para fase instantánea con Morlet compleja. Gurnett, D. A. and Bhattacharjee, A. (2005). Introduction to Plasma Physics. Cambridge. Gurnett, D. A., Kurth, W. S., Burlaga, L. F., and Ness, N. F. (2013a). Plasma wave observations from voyager 1 at the heliopause. Science, 341(6153):1489–1492. Gurnett, D. A., Kurth, W. S., Burlaga, L. F., and Ness, N. F. (2013b). Plasma wave observations from voyager 1 at the heliopause. Science, 341:1489–1492. Howes, G. G. (2015). The damped propagation of kinetic alfvén waves. Physics of Plasmas, 22:032901. Priest, E. and Forbes, T. (2000). Magnetic Reconnection: MHD Theory and Applications. Cambridge Univ. Press. Stix, T. H. (1992). Waves in Plasmas. Springer. Stone, E. C., Cummings, A. C., McDonald, F. B., Heikkila, B. C., Lal, N., and Webber, W. R. (2013). Voyager 1 observes an increase in the flux of galactic cosmic rays at the edge of the heliosphere. Science, 341(6153):150–153. Torrence, C. and Compo, G. P. (1998). A practical guide to wavelet analysis. Bulletin of the American Meteorological Society, 79(1):61–78. El artículo clásico de referencia para CWT y wavelet de Morlet en geofísica. 49