On the Theory of Bulk Viscosity of Cold Plasmas and Thermodynamics of Alkali-Noble Gas Cocktails

Dit artikel leidt een uitdrukking af voor de volumebinding in koude plasma's, waarbij wordt aangetoond dat deze veel groter kan zijn dan de schuifviscositeit en dat de Mandelstam-Leontovich-benadering exact is, met toepassingen voor akoestische verwarming van de zonne-atmosfeer en thermodynamica van alkali-edelgas-mengsels.

De kern

De Geluidsdempende "Kleefkracht" van Sterrenstof: Een Simpele Uitleg

Stel je voor dat je een glas champagne proost. Je hoort een helder, rinkelend geluid. Nu stel je je voor dat je datzelfde glas vult met een mengsel van zout en water, of met een heel speciaal soort "sterrensoep" (een plasma zoals in de zon). Als je daarop zou tikken, zou het geluid misschien niet zo helder klinken, maar juist dof worden, alsof het geluid wordt opgeslokt door een onzichtbare spons.

Wetenschappers Albert Varonov en Todor Mishonov hebben in dit artikel ontdekt waarom dat gebeurt in koude plasma's (zoals in de buitenste atmosfeer van de zon) en hoe ze dit kunnen berekenen.

Hier zijn de belangrijkste punten, vertaald naar alledaagse taal:

1. De Twee Soorten "Viskeuze" Krachten

In de natuurkunde kennen we twee soorten "stroperigheid" (viscositeit) in vloeistoffen en gassen:

• Scherende viscositeit (Shear viscosity): Dit is als het wrijven van je handen tegen elkaar. Het is de weerstand die je voelt als je een vloeistof probeert te schuiven of te roeren. Dit kennen we allemaal.
• Bulkviscositeit (Bulk viscosity): Dit is de weerstand die optreedt als je iets in- en uitademt (samendrukt en uitbreidt). Denk aan een spons die je in elkaar knijpt en weer laat los. Als de spons niet snel genoeg water kan opzuigen of loslaten, voelt dat als een soort interne weerstand.

Het grote geheim: In de meeste gassen (zoals de lucht die we inademen) is de eerste soort veel belangrijker. Maar in de "koude plasma's" van de zon (waar atomen worden gesplitst in elektronen en ionen) is de bulkviscositeit duizenden tot miljoenen keren sterker dan de scherende viscositeit! Het is alsof je probeert een auto te duwen, maar de wielen zijn vastgeplakt aan de grond door honing.

2. De "Kleefkracht" van de Zon

Waar komt deze enorme weerstand vandaan?
Stel je voor dat je een groep mensen (atomen) in een kamer hebt. Als je de kamer snel kleiner maakt (geluidsgolf), moeten de mensen snel van de ene stoel naar de andere springen (ionisatie: atomen verliezen een elektron).

• Als je de kamer heel langzaam kleiner maakt, springen ze rustig over. Alles is in evenwicht.
• Als je de kamer heel snel kleiner maakt, blijven ze achter op hun stoel. Ze hebben tijd nodig om te reageren.

Die vertraging – het feit dat de atomen niet direct kunnen "reageren" op de drukverandering – creëert die enorme weerstand (bulkviscositeit). In de zon zorgt dit ervoor dat geluidsgolven hun energie kwijtraken aan het "reorganiseren" van de atomen. Die energie wordt omgezet in warmte.

3. De Zon als een Geluidsoven

De auteurs gebruiken dit om een oud mysterie op te lossen: Waarom is de atmosfeer van de zon zo heet?
De zon is aan de buitenkant (fotosfeer) ongeveer 5500 graden, maar hoger in de atmosfeer (chromosfeer) stijgt de temperatuur plotseling naar tienduizenden graden. Hoe kan dat?
Het antwoord in dit artikel: Geluid.
De zon "zingt" voortdurend. Geluidsgolven reizen omhoog. In de koude lagen van de zon (waar de temperatuur lager is dan de energie die nodig is om atomen te splitsen), botsen deze golven tegen de "ionisatie-muur". Door de enorme bulkviscositeit worden deze geluidsgolven gedempt en verandert hun bewegingsenergie in warmte. Het is alsof de zon zichzelf verwarmt met een onzichtbare geluidsoven.

4. De "Wannier-berg" en de Wiskunde

Om dit te berekenen, kijken de auteurs naar hoe elektronen atomen raken. Ze gebruiken een wiskundig concept dat de "Wannier-berg" heet.

• Vergelijking: Stel je voor dat een elektron een berg beklimt om een atoom te raken. Op de top van die berg (de drempel) gebeurt er iets speciaals: twee elektronen bewegen als spiegelbeelden van elkaar.
• Dankzij deze specifieke wiskundige regel kunnen de auteurs een heel nauwkeurige formule maken. Ze ontdekten dat een oude theorie uit de jaren '30 (van Mandelstam en Leontovich) in dit geval perfect werkt. Het is alsof ze een oude, vergeten sleutel vonden die precies in het slot van dit complexe probleem past.

5. Wat betekent dit voor ons?

• Voor de sterrenkunde: Het bevestigt dat geluidsgolven een cruciale rol spelen bij het verwarmen van de zon. We moeten deze "bulkviscositeit" meenemen in onze modellen, anders begrijpen we de zon niet goed.
• Voor de toekomst: De auteurs suggereren dat we dit in het lab kunnen nabootsen met mengsels van edelgassen en alkali-metalen (zoals natrium en neon). Als we daar geluid door sturen, kunnen we zien hoe het opwarmt. Dit zou kunnen helpen bij het begrijpen van hoe we plasma's kunnen beheersen, misschien zelfs voor toekomstige energiebronnen of raketten.

Samenvatting in één zin:

De zon is warm omdat geluidsgolven erin vastlopen in een soort "moleculaire modder" (bulkviscositeit), waardoor de atomen moe worden en de energie afgeven als warmte; en de auteurs hebben de perfecte wiskundige formule gevonden om dit proces te beschrijven.

Technische samenvatting

Probleemstelling

De auteurs adresseren een fundamenteel probleem in de plasmafysica en de astrofysica: het ontbreken van een eerste-principes berekening voor de bulkviscositeit (volumeviscositeit, $\zeta$) in koud plasma. Hoewel bulkviscositeit in veel systemen (zoals de zonnestraling, oceanen en deeltjesversnellers) dominant kan zijn ten opzichte van de schuifviscositeit ($\eta$), wordt deze vaak verwaarloosd of benaderd met onnauwkeurige modellen.

Specifiek richten de auteurs zich op:

1. Het verklaren van het mechanisme voor akoestische verwarming van de zonnechromosfeer.
2. Het afleiden van een expliciete formule voor bulkviscositeit in koud plasma, gedefinieerd als plasma waar de temperatuur $T$ veel lager is dan de eerste ionisatiepotentiaal $I_a$ ($T \ll I_a$).
3. Het valideren of de klassieke Mandelstam-Leontovich (ML) benadering (die uitgaat van één relaxatietijd) geldig is voor deze systemen.

Methodologie

De auteurs hanteren een combinatie van kinetische theorie, thermodynamica en numerieke simulatie:

1. Kinetische Vergelijking: Ze lossen de kinetische vergelijking op voor ionisatie- en recombinatieprocessen in een koud plasma. Ze gebruiken de Saha-vergelijking voor evenwicht en lineairiseren deze voor kleine oscillaties in dichtheid en temperatuur.
2. Wannier-Nabijheid: Voor de ionisatiesnelheid ($\beta$) gebruiken ze het Wannier-model voor de ionisatiewerkingsdoorsnede vlak boven de drempelwaarde. Dit is een analytische oplossing die geldt voor $s$-elektronen en zeer nauwkeurig is voor koud plasma.
3. Thermodynamica: Ze berekenen de thermodynamische variabelen (enthalpie, interne energie, warmtecapaciteit $C_p$ en $C_v$) voor een mengsel van alkali- en edelgassen (geïllustreerd met een H-He mengsel voor de zon).
4. Lineaire Respons: Ze analyseren de reactie van het plasma op harmonische drukoscillaties. Dit leidt tot een gegeneraliseerde complexe polytropische index $\hat{\gamma}(\omega)$, die de dissipatie door ionisatie/recombinatie beschrijft.
5. Numerieke Implementatie: Ze lossen een stelsel lineaire vergelijkingen op voor een cocktail van 6 elementen (H, He, C, Mg, Si, Fe) gebaseerd op het AL08-model van de zonnechromosfeer.

Belangrijkste Bijdragen en Resultaten

1. Dominantie van Bulkviscositeit

Het meest opvallende resultaat is dat de bulkviscositeit $\zeta$ in koud plasma vele ordes van grootte groter kan zijn dan de schuifviscositeit $\eta$.

• De verhouding $P_{\zeta/\eta} = \zeta_0 / \eta$ bereikt waarden tot $10^{11}$ in de zonnechromosfeer.
• Dit betekent dat voor geluidsgolven in dit regime de demping en verwarming bijna uitsluitend worden bepaald door bulkviscositeit en niet door schuifviscositeit of warmtegeleiding.

2. Validatie van de Mandelstam-Leontovich (ML) Benadering

De auteurs tonen aan dat de ML-benadering, die de frequentieafhankelijkheid van de viscositeit beschrijft met één relaxatietijd $\tau$, in de praktijk exact is voor koud plasma.

• De berekende complexe polytropische index $\hat{\gamma}(\omega)$ valt perfect samen met de ML-curve (een cirkelvormige Cole-Cole plot).
• Voor zuiver waterstofplasma is de ML-benadering wiskundig exact; voor mengsels is de afwijking verwaarloosbaar klein (< 2%).
• Dit vereenvoudigt de modellering aanzienlijk, omdat men niet de volledige kinetische vergelijking hoeft op te lossen voor elke frequentie, maar slechts de relaxatietijd $\tau$ en de limietwaarden $\gamma_0$ en $\gamma_\infty$ nodig heeft.

3. Expliciete Formules voor Relaxatietijd en Viscositeit

Ze leiden expliciete analytische formules af voor:

• De relaxatietijd $\tau$, die afhankelijk is van de ionisatiegraad $\alpha$ en de temperatuur.
• De statische bulkviscositeit $\zeta_0 = \tau p (\gamma_\infty - \gamma_0)$.
• De frequentieafhankelijke demping van geluidsgolven.

4. Akoestische Verwarming van de Zon

De resultaten hebben directe implicaties voor de verwarming van de zonnechromosfeer:

• Bij lage frequenties ($f < f_c$) wordt de demping van longitudinale geluidsgolven gedomineerd door bulkviscositeit.
• Er bestaat een kritieke frequentie $f_c$ (ongeveer 35 mHz in hun voorbeeld) waarbij bulkviscositeit en schuifviscositeit gelijk zijn. Boven deze frequentie domineert schuifviscositeit.
• De auteurs concluderen dat akoestische golven met een significant spectrum onder $f_c$ een cruciale rol spelen bij het verwarmen van de chromosfeer via ionisatie-recombinatie dissipatie.

Significantie en Toepassingen

1. Astrofysica: Het onderzoek biedt een robuust theoretisch kader om het "verwarmingsprobleem" van de zonnechromosfeer op te lossen. Het toont aan dat bulkviscositeit een onmisbaar mechanisme is in modellen voor akoestische verwarming, wat eerder vaak werd genegeerd.
2. Laboratoriumplasma's: De theorie is niet beperkt tot de zon. De afgeleide formules zijn toepasbaar op laboratoriumplasma's van alkali-noble gas mengsels (bijv. Na-Ne). Dit biedt een route om de theorie experimenteel te verifiëren door akoestische verwarming in dergelijke mengsels te meten.
3. Methodologische Doorbraak: Het is een van de weinige gevallen waarbij een kinetische coëfficiënt (bulkviscositeit) volledig "ab initio" kan worden berekend voor een complex chemisch systeem (ionisatie/recombinatie) zonder empirische aanpassingen.
4. Vereenvoudiging van Modellen: De bevestiging dat de ML-benadering exact is, stelt onderzoekers in staat om complexe hydrodynamische simulaties van plasma's te vereenvoudigen door gebruik te maken van de Drude-achtige frequentieafhankelijkheid voor bulkviscositeit.

Conclusie:
Varonov en Mishonov hebben succesvol een brug geslagen tussen microscopische kinetische processen (ionisatie) en macroscopische hydrodynamische eigenschappen (bulkviscositeit). Hun werk bewijst dat bulkviscositeit de dominante dissipatiemechanisme is in koud, gedeeltelijk geïoniseerd plasma en dat dit mechanisme essentieel is voor het begrijpen van de energietransportprocessen in de zonneatmosfeer.