Global Coronal Equilibria with Solar Wind Outflow II -- Optimizing the Outflow Model

Probleemstelling

Traditionele modellen voor het extrapoleren van het magnetische veld van de zonnekorona, zoals het Potential Field Source Surface (PFSS) model, hebben aanzienlijke tekortkomingen:

1. Het "Open Flux Problem": PFSS-modellen onderschatten de totale Open Solar Flux (OSF) – de totale onbevestigde magnetische flux die de zon verlaat – met ongeveer 45% ten opzichte van in-situ metingen in de interplanetaire ruimte.
2. Afhankelijkheid van de Bronoppervlakte: De resultaten van PFSS-modellen zijn sterk afhankelijk van de willekeurig gekozen hoogte van het "bronoppervlak" (meestal $r_{ss} = 2.5 R_\odot$). Als deze hoogte verandert, verandert de topologie van het magnetische veld drastisch, wat fysisch onrealistisch is.
3. Onrealistische Structuren: PFSS-modellen voorspellen vaak dat magnetische stromers (streamers) zich uitstrekken tot aan het bronoppervlak, wat niet overeenkomt met waarnemingen tijdens zonsverduisteringen.

In eerdere werken (Paper I) introduceerden de auteurs "Outflow Fields" als alternatief, maar deze paper (Paper II) heeft tot doel deze verbeteringen te kwantificeren en de modelparameters te optimaliseren.

Methodologie

De auteurs gebruiken een magnetofriktie (MF) model dat een evenwichtszustand zoekt voor het magnetische veld, inclusief een term voor de zonnewind, zonder de noodzaak voor kostbare tijdsafhankelijke simulaties.

1. Het Outflow Model:

   • In tegenstelling tot PFSS ($\mathbf{B} = \nabla \Phi$), laat het Outflow-model meer radiale variatie toe: $\mathbf{B} = f(r)\nabla \Phi$.
   • Het model lost de vergelijking op: $\left( \frac{(\nabla \times \mathbf{B}) \times \mathbf{B}}{\nu_0 |\mathbf{B}|^2} + v(r)\mathbf{e}_r \right) \times \mathbf{B} = 0$, waarbij $v(r)$ het zonnewindprofiel is en $\nu_0$ de magnetofriktieconstante.
   • In plaats van een vaste Parker-windoplossing, wordt een dimensieloze snelheid $V(r)$ gebruikt die empirisch wordt geoptimaliseerd.

2. Data Voorbereiding:

   • Magnetogramdata van HMI (SDO) en MDI (SOHO) worden verwerkt (downsampled, gladgestreken met Legendre-polynomen en gefilterd) om discontinuïteiten te elimineren.
   • Er wordt gebruik gemaakt van in-situ OSF-data van de Wind en ACE-satellieten (gecorrigeerd door Frost et al., 2022) als "ground truth".

3. Optimalisatie via Evolutie-algoritme:

   • Om de parameters te optimaliseren, vergelijken de auteurs de hoeken van de magnetische veldlijnen in hun modellen met waarnemingen van 12 zonsverduisteringen (2006-2024).
   • Een CMA-ES (Covariance Matrix Adaptation Evolution Strategy) algoritme wordt gebruikt om het zonnewindprofiel $V(r)$ te vinden dat de gemiddelde afwijking van de veldlijnen ten opzichte van de radiale richting minimaliseert, in vergelijking met de "ground truth" uit verduisteringsfoto's.
   • Er worden synthetische verduisteringsafbeeldingen gegenereerd om de topologie visueel te vergelijken.

Belangrijkste Bijdragen

1. Onafhankelijkheid van de Bronoppervlakte: Het paper toont aan dat bij Outflow Fields de OSF en de veldtopologie veel minder afhankelijk zijn van de gekozen hoogte van het bronoppervlak ($r_{ss}$) dan bij PFSS. De zonnewind bepaalt de structuur, niet de hoogtegrens.
2. Empirische Optimalisatie: Voor het eerst wordt een methode gepresenteerd om de zonnewindprofielen te optimaliseren op basis van de topologie van de korona tijdens verduisteringen, in plaats van alleen op basis van OSF-metingen.
3. Open Source Software: De auteurs introduceren het Python-pakket 'outflowpy', dat ontworpen is om compatibel te zijn met het bestaande 'pfsspy'-pakket, waardoor de gemeenschap direct gebruik kan maken van deze verbeterde modellen.

Resultaten

• Verbetering van de OSF-schatting:
  • PFSS onderschat de OSF met ongeveer 45% (bij $r_{ss} = 2.5 R_\odot$).
  • Een geoptimaliseerd Outflow-model vermindert deze discrepantie tot ongeveer 24%.
  • Hoewel het model de OSF nog niet perfect voorspelt, is dit een aanzienlijke verbetering. Het paper stelt dat het volledig matchen van de OSF vaak leidt tot onrealistisch radiale veldlijnen, wat de topologie verstoort.
• Topologische Overeenkomst:
  • Geoptimaliseerde Outflow Fields matchen de hoeken van de magnetische veldlijnen (gemeten tijdens verduisteringen) veel nauwkeuriger dan PFSS-modellen.
  • PFSS-modellen kunnen de waargenomen niet-radiale veldlijnen niet reproduceren, ongeacht de gekozen hoogte van het bronoppervlak.
  • De geoptimaliseerde modellen tonen stromers die zich uitstrekken tot 1-3 $R_\odot$, wat overeenkomt met waarnemingen, terwijl PFSS-stromers vaak tot aan de bronoppervlakte (bijv. 5 $R_\odot$) doorgaan.
• Correlatie: De correlatiecoëfficiënt tussen de voorspelde OSF en in-situ metingen stijgt van 0.748 (PFSS) naar 0.864 (geoptimaliseerd Outflow).

Betekenis en Conclusie

Deze paper bevestigt dat Outflow Fields een superieur alternatief zijn voor PFSS-modellen voor het modelleren van het globale magnetische veld van de zon. Ze bieden een betere balans tussen rekenkundige efficiëntie en fysische realisme.

• Praktische Toepassing: Het model kan direct worden gebruikt in plaats van PFSS voor ruimteweersvoorspellingen en als initiële conditie voor complexere MHD-simulaties (zoals aangetoond door Li et al., 2025).
• Beperkingen: Het model onderschat de OSF nog steeds enigszins. De auteurs suggereren dat dit kan worden veroorzaakt door niet-potentiële structuren in de lagere korona (zoals flux ropes en CME's) of door kalibratieproblemen in de magnetogramdata (HMI vs. MDI).
• Toekomst: Hoewel het huidige model alleen afhankelijk is van de hoogte, zou latere afhankelijkheid van de zonnewindsnelheid (bijv. op breedtegraad) de nauwkeurigheid verder kunnen verbeteren.

Kortom, het paper levert een robuust, geoptimaliseerd model dat de "Open Flux Problem" aanzienlijk vermindert en de topologie van de korona realistischer weergeeft dan eerdere benaderingen.