Alomtegenwoordige Kelvin-Helmholtz-instabiliteiten drijven plasmamenging aan op de zon

Gepubliceerd: 5 augustus 2026 Tijdschrift: Nature, volume 656, pagina's 595–601

Samenvatting

De interactie tussen het magnetisch veld en turbulente convectie in de fotosfeer van de zon drijft de dynamiek, evolutie en structurering van haar gemagnetiseerde atmosfeer aan. Deze interactie vindt vaak plaats op of onder de ruimtelijke resolutie van moderne waarnemingen. In dit artikel rapporteren we over waarnemingen met een hoge ruimtelijke resolutie van de zonnefotosfeer, verkregen met de eerste 4-meter klasse zonnetelescoop ter wereld, de Daniel K. Inouye Solar Telescope (DKIST) van de National Science Foundation van de VS.

Tijdreeksbeelden onthullen een veel complexer en dynamischer zonnebeeld dan voorheen is waargenomen. We identificeren alomtegenwoordige gemagnetiseerde Kelvin-Helmholtz-instabiliteiten (KHI) aan de randen van magnetische fluxconcentraties en leveren experimentele bevestiging van een langdurige theoretische voorspelling. De ontdekking van kleinschalige gemagnetiseerde Kelvin-Helmholtz-instabiliteiten in de zonnefotosfeer, die gereproduceerd kunnen worden door numerieke simulaties met een hoge resolutie, heeft verstrekkende gevolgen voor ons begrip van de creatie en dissipatie van magnetische velden die vortexbewegingen vertonen, wat kan leiden tot fluxvlechting (flux braiding).

Onze resultaten ondersteunen het beeld van losgekoppelde magnetische veldconcentraties in lagen onder het zichtbare zonneoppervlak, die verbonden zijn met monolithische fluxregio's die zichtbaar zijn als faculaire concentraties en poriën in de zonnefotosfeer. Kelvin-Helmholtz-instabiliteiten zijn een efficiënt mechanisme voor het transport van massa, energie, impuls en magnetische flux in magnetohydrodynamische systemen, en ze bieden transformatieve inzichten in processen in magnetisch actieve regio's.

Inleiding

Sinds hun oorspronkelijke formulering door Lord Kelvin en Hermann von Helmholtz zijn Kelvin-Helmholtz-instabiliteiten (KHI) waargenomen en onderzocht in vele gebieden van de fysica, waaronder vloeistofdynamica, oceanografie, heliosferische omgevingen zoals planetaire magnetosferen (bijvoorbeeld Aarde, Mars, Venus, Mercurius, Jupiter en Saturnus) en andere astrofysische systemen. In de zonneatmosfeer zijn KHI's al lang het onderwerp van zowel theoretische als observationele studies. Gegolfde grenzen die waargenomen worden in coronale massa-uitbarstingen (coronal mass ejections) worden vaak geïnterpreteerd als mogelijke tekenen van deze instabiliteit.

Zonne- en stellaire magneto-convectie produceren complexe plasmastromen in de fotosfeer die evolueren van coherente laminaire bewegingen naar kleinschalige, instabiele structuren. Theoretische studies geven aan dat deze stromingen in de zonnefotosfeer snelheids-shear-lagen (velocity shear layers) rond magnetische fluxconcentraties (MFC's) kunnen genereren en de KHI kunnen triggeren. Het bestaan van KHI's in de fotosferische lagen van de zon is in het verleden niet waargenomen, waarschijnlijk omdat de karakteristieke ruimtelijke schaal niet gemakkelijk toegankelijk is voor telescopen met een apertuur kleiner dan 2 meter.

Waarnemingen

De waarnemingen die hier worden gerapporteerd, zijn gemaakt met de Daniel K. Inouye Solar Telescope (DKIST) bij een golflengte van 416 nm in een magnetisch actieve regio nabij een zonnevlek (Fig. 1). Bij deze golflengte wordt het zichtbare zonneoppervlak gedomineerd door convectieve bewegingen (granules) vermengd met kleinschalige MFC's. Dit tafereel bevat verschillende grotere coherente structuren (poriën).

De interfaces tussen de magnetische kenmerken en de granulatie in deze beelden vertonen een structuur die zeer verschilt van de relatief gladde, wazige interfaces die gewoonlijk worden gezien in beelden met een lagere resolutie. DKIST onthult dat deze interfaces bijna volledig zijn samengesteld uit vortex-achtige structuren en striaties (strepen), wat het zonnebeeld drastisch verandert ten opzichte van eerdere waarnemingen.

De tijdreeks van beelden onthult de alomtegenwoordige vorming en evolutie van vortices aan de interface van magnetische elementen, waarbij de kleinste vortices zich dicht bij de theoretische diffractielimiet van DKIST bevinden (ongeveer 19 km op de zon bij deze golflengte). Een analyse van 47 vortices in het waargenomen veld toont aan dat de karakteristieke ruimtelijke golflengte — de ruimtelijke afstand tussen vortices — 65 km bedroeg, met formaten variërend van 25 km tot 170 km.

De groeisnelheid $\gamma_{\mathrm{im}}$, gedefinieerd als de exponentiële factor waarmee de vervorming van de magnetische interface in de loop van de tijd toeneemt, heeft waarden tussen $0,014\text{ s}^{-1}$ en $0,054\text{ s}^{-1}$ voor prominente voorbeelden in het gezichtsveld (FoV). De gemeten schijnbare snelheden (de horizontale voortplantingssnelheid van vortices rond de magnetische elementen) varieerden tussen $0,67\text{ km s}^{-1}$ en $3,0\text{ km s}^{-1}$.

Numerieke Simulaties

We hebben de waargenomen structuren vergeleken met numerieke stralings-magnetohydrodynamische (MHD) simulaties van de zonneatmosfeer, uitgevoerd met de MURaM-code. De simulaties modelleren een plage-regio met een netto flux, beginnend vanuit een aanvankelijk verticaal veld dat gelijk is aan dat van de waargenomen scène, berekend met een ruimtelijke discretisatiegrootte van 3,2 km.

De numerieke simulaties vertonen structuren die overeenkomen met die in de waarnemingen. De prominente vorm van de vortices in zowel de observationele als synthetische data komt overeen met de morfologische structuur van de KHI.

Dynamica van de instabiliteit

De simulaties tonen een belangrijk ingrediënt dat niet door de waarnemingen is gemeten: de horizontale plasmasnelheid. Deze loopt parallel aan de scherpe interface van het sterke magnetische veld en vertoont daar een sterk gradiënt, waardoor de optimale omgeving voor het ontstaan van KHI wordt gecreëerd.

Magnetische velden kunnen een stabiliserende invloed hebben op de KHI: de component parallel aan de shear-stroom neigt de instabiliteit te onderdrukken, terwijl de loodrechte component geen stabilisatie biedt. De sterke velden in de waargenomen en gesimuleerde plage-regio's zijn aan het zonneoppervlak bijna verticaal en dus hoofdzakelijk loodrecht op de horizontale shear-stroom. De lineaire theorie voorspelt dat de groei van de instabiliteit wordt gemaximaliseerd wanneer de Kelvin-Helmholtz-golfvector is uitgelijnd met de stroomrichting.

Analyse en Resultaten

Vergelijking tussen simulaties en waarnemingen

We hebben de instabiliteitsparameters in de MURaM-simulaties onderzocht met dezelfde technieken als bij de waargenomen beelden. De resultaten tonen aan dat het histogram piekte bij een karakteristieke instabiliteitsgolflengte van 49 km. De groeisnelheid voor verschillende opvallende voorbeelden in de synthetische beelden lag in het bereik van $0,027\text{ s}^{-1}$ tot $0,059\text{ s}^{-1}$. De instabiliteitsgolflengte en groeisnelheid in de synthetische beelden zijn vergelijkbaar met die in de waarnemingen, evenals de gemeten schijnbare vortexsnelheden ($1,6\text{ km s}^{-1}$ tot $2,8\text{ km s}^{-1}$).

Fysische mechanismen

De KHI-structuren treden snel in het niet-lineaire regime. Ze vertonen volledig ontwikkelde vortices samen met snelle vortex-fusie en de vorming van secundaire, kleinschaligere KHI's op de grotere vortices. Dit toont aan dat de KHI een kritische bijdrage levert aan de ontwikkeling van turbulentie in deze grenslaag.

De KHI's kunnen worden beschouwd als een bron van turbulente diffusiviteit $\eta_{\text{tur}} \approx 1/3 u \lambda$. Met een golflengte $\lambda \approx 65\text{ km}$ en gemiddelde shear-snelheid $u \approx 3\text{ km s}^{-1}$ resulteert dit in een geschatte turbulente diffusiviteit van $0,65 \times 10^{12}\text{ cm}^2\text{ s}^{-1}$ in regio's met een sterk magnetisch veld. Deze diffusiviteit leidt tot menging tussen gemagnetiseerd en grotendeels ongemagnetiseerd plasma.

Magnetische stratificatie en striaties

In de MURaM-simulaties ontwikkelen KHI-vortices zich in horizontale vlakken over een verticaal dieptebereik $z \in \sim [100, -400]\text{ km}$. In de diepere lagen ($z < 0\text{ km}$) is de magnetische stabilisatie zwakker, waardoor de KHI leidt tot fragmentatie en het monolithische magnetische element opsplitst in verschillende kleinere structuren. Deze magnetische stratificatie vertoont gelijkenissen met het model van Parker voor zonnevlekken.

Daarnaast tonen zowel de waarnemingen als de synthetische kaarten fijne donkere striaties die geworteld zijn in de KHI-vortices. We tonen aan dat deze variaties voortkomen uit de KHI, die de grenzen van de MFC's golvend maakt.

Implicaties voor de Zonneatmosfeer

Magnetische structuren zoals spicules, prominenties en coronale lussen zijn geworteld in fotosferische MFC's. De KHI kan, door magnetische elementen continu te vervormen, MHD-golven exciteren en energie dissiperen via vortices die vervolgens kunnen uiteenvallen in turbulentie.

De voetpuntbewegingen van magnetische structuren in de fotosfeer zijn de primaire bron van vrije magnetische energie, die de meeste zonneactiviteit aandrijft (van nano-flares tot flares, jets en coronale massa-uitbarstingen). Een centraal theoretisch concept hiervoor is magnetische fluxvlechting (flux braiding). De waargenomen vortices gegenereerd door KHI zouden een efficiënte bron kunnen zijn van kleinschalige, alomtegenwoordige magnetische fluxvlechting.

Methoden

Waarnemingen en datareductie

De waarnemingen zijn gedaan op 14 april 2025 tussen 21:38 UT en 21:41 UT met DKIST. De focus lag op een regio met poriën nabij een actieve regio (NOAA 14060). De gebruikte golflengte was 416 nm met een smalle banddoorlaat van 0,5 nm. De camera behaalde een pixelgrootte van $\approx 6\text{ km px}^{-1}$ op het zonneoppervlak. Beelden werden gereconstrueerd met multi-frame blind deconvolution, wat resulteerde in een ruimtelijke resolutie van 19 km.

Numerieke simulaties en spectrale synthese

Simulaties zijn uitgevoerd met de MURaM-code in een domein van $6,144 \times 6,144 \times 2,048\text{ Mm}^3$. De simulatie startte met een verticaal magnetisch veld gebaseerd op HMI-magnetogrammen. De grid-schaal was 3,2 km. Spectrale synthese werd uitgevoerd met de Rybicki–Hummer code in lokaal thermodynamisch evenwicht (LTE).

Lineaire MHD-theorie voor KHI-groei

De lineaire theorie stelt dat de instabiliteit het snelst groeit wanneer de shear-stroom loodrecht op het magnetische veld staat. Voor een slab-geometrie met een eindige breedte van de shear-laag $d$ is de dispersierelatie:

$${{\rm{e}}}^{-2\kappa }=\left[1-\frac{\kappa {(\nu +1)}^{2}}{(\nu +1)+1/2{\epsilon }\kappa {(\nu +1)}^{2}}\right]\cdot \left[1+\frac{\kappa {(\nu -1)}^{2}}{(\nu -1)+1/2{\epsilon }\kappa {(\nu -1)}^{2}}\right]$$

Waarbij:

• $\nu = \gamma / k U_0$
• $\kappa = kd$
• $\gamma$ is de frequentie
• $k$ is het golfgetal
• $\rho_0$ is de dichtheid in de shear-laag
• $\epsilon$ is de dichtheidscontrastfactor

Schatten van KHI-parameters

De groeisnelheid $\gamma{\mathrm{im}}$ werd gemeten uit de temporele verandering in de amplitude $A(t)$ van de instabiliteit tijdens de lineaire fase: $A(t)=A({t}{0})\,\exp ({\gamma }_{{\rm{i}}{\rm{m}}}t)$. De fasesnelheid werd bepaald met tijd-afstand-diagrammen.

Referenties

1. Karpen, J. T., et al. (1993). Astrophys. J. 403, 769–784.
2. Kolesnikov, F., et al. (2004). Astron. Astrophys. 420, 737–746.
3. Thomson, W. (1871). Philos. Mag. 42, 362–377.
4. Helmholtz, H. L. F. (1868). Monatsberichte K. Preuss. Akad. Wiss. Berl. 23, 215–228.
5. Drazin, P. G. & Reid, W. H. (1981). Hydrodynamic Stability. Cambridge Univ. Press.

... (volledige lijst beschikbaar in het originele document)