2025. július 30., szerda

Analitikai Modell a Kationcsere Dinamikájához Bolygóköpeny-szilikátokban, Változó P-T Körülmények Között

 

Analitikai Modell a Kationcsere Dinamikájához Bolygóköpeny-szilikátokban, Változó P-T Körülmények Között

Lengyel Ferenc

2025. július

Vezetői Összefoglaló

A köpenyásványokon belüli ioncsere alapvető folyamat, amely a Föld-típusú bolygók kémiai evolúcióját szabályozza. Olyan nagyléptékű jelenségek alapját képezi, mint a bolygódifferenciálódás és a köpenykonvekció, azáltal, hogy meghatározza a főbb elemek, például a magnézium, a vas, a kalcium és a nátrium eloszlását.1 Fontossága ellenére egy szigorú, első elveken alapuló analitikai modell, amely leírja ezen cserefolyamatok dinamikáját egy bolygóköpeny extrém nyomás-hőmérséklet (P-T) körülményei között, megoldatlan kihívás maradt. Ez a jelentés egy újszerű elméleti keretrendszert mutat be e hiányosság pótlására. Kidolgozunk egy átfogó analitikai modellt a kulcsfontosságú köpeny-szilikátásványokban zajló intrakristályos (kristályon belüli) kationcsere-dinamikára (különösen a

Mg2+↔Fe2+ és Ca2+↔Na+ cserére).

A modellt kapcsolt, nemlineáris parciális differenciálegyenletek rendszereként fogalmazzuk meg a Nernst-Planck-Poisson formalizmus alapján, amely önkonzisztens módon veszi figyelembe az ionos diffúziót, a belső diffúziós potenciálok okozta elektromigrációt és a csere-reakció termodinamikai hajtóerőit.2 A termodinamikai komponens fejlett, nemideális szilárdoldat-modelleket tartalmaz a ferromágneses szilikátokra és a földpátokra, ahol az egyensúlyi megoszlás explicit módon függ a köpeny változó P-T körülményeitől, beleértve a nagynyomású fázisátalakulásokat is.

E munka központi eredménye egy egzakt analitikai sorfejtéses megoldás levezetése erre a komplex, nemlineáris rendszerre a Homotópia Analízis Módszer (HAM) segítségével.4 Ez a hatékony matematikai technika a nemlineáris problémát megoldható lineáris egyenletek végtelen sorozatává alakítja át, egy konvergens sorfejtéses megoldást eredményezve, amelynek pontossága szisztematikusan szabályozható. Az analitikai megoldás példátlan betekintést nyújt a kémiai egyensúly beállásának karakterisztikus időskáláiba a hőmérséklet, a nyomás és a szemcseméret függvényében.

A modell alkalmazása geodinamikailag releváns forgatókönyvekre, mint például a köpenycsóvák és az alábukó lemezek, lehetővé teszi az ásványokban a köpenyen való áthaladásuk során megőrzött kémiai egyensúlyhiány számszerűsítését. Ez kvantitatív keretet teremt a köpenykőzetek „kémiai memóriájának” megértéséhez.5 Továbbá a modell rávilágít az ásványi szintű kémia és a bolygószintű dinamika közötti visszacsatolási hurokra: az ioncsere megváltoztatja az ásványok összetételét, ami módosítja a kőzetek sűrűségét és reológiáját, ezáltal befolyásolva a köpenykonvekció mintázatait és hevességét.6 Ez a munka áthidalja az ásványfizika és a geodinamika közötti szakadékot, kritikus fontosságú rácson belüli (sub-grid) fizikai modult biztosítva a bolygófejlődés következő generációs numerikus szimulációihoz, és új elméleti eszközt kínálva a köpenyből származó anyagok geokémiai megfigyeléseinek értelmezéséhez.7

1. Szekció: A Kationcsere Termodinamikai és Krisztallográfiai Alapjai Köpenyásványokban

Az ioncsere kinetikáját alapvetően a minimális Gibbs-féle szabadenergia állapotának elérésére irányuló törekvés szabályozza.9 Ezért a dinamikus modell kidolgozása előtt elengedhetetlen egy robusztus termodinamikai keretrendszer létrehozása, amely pontosan meghatározza a kationeloszlás egyensúlyi állapotát a köpenyásványokban a hőmérséklet (T), a nyomás (P) és a teljes összetétel függvényében. Ez a fejezet a kémiai megoszlás általános elveitől a

Mg2+↔Fe2+ és Ca2+↔Na+ csere-rendszerek specifikus, komplex termodinamikája felé halad a szilikátásványokban. A végső cél egy fizikailag megalapozott, P-T függő kifejezés levezetése az egyensúlyi csere-együtthatóra, a KD-re, amely a 2. Szekcióban kifejlesztett dinamikus transzportmodell termodinamikai célállapotaként szolgál majd.

1.1. Az Egyensúlyi Megoszlás és a Csere-reakciók Alapelvei

Egy kémiai fajta eloszlása két egymás mellett létező fázis között egyensúlyi állapotban a geokémia egyik sarokköve.1 Ezt az eloszlást a megoszlási együttható, D, számszerűsíti, amelyet egy i fajta koncentrációjának (C) arányaként definiálnak az

α fázisban és a β fázisban 11:

Diα/β=CiβCiα

Bár a megoszlási együttható felbecsülhetetlen értékű a nyomelemek vizsgálatában, e jelentés fókuszában a főbb, szerkezetileg hasonló kationok cseréje áll egyetlen ásványfázison belül. Az ilyen folyamatok esetében a relevánsabb paraméter a csere-együttható, amelyet gyakran KD-vel jelölnek, és amely egy csere-reakció egyensúlyi állapotát írja le.1 Az A és B kationok általános cseréjére egy ásványon belül két különböző krisztallográfiai hely, 1 és 2 között, a reakció a következőképpen írható fel 13:

A(2. hely)+B(1. hely)⇌A(1. hely)+B(2. hely)

A reakció csere-együtthatóját a tömeghatás törvénye határozza meg, a kationok móltörtjeit (X) felhasználva az egyensúlyi állapotban 15:

KD=XA2. hely⋅XB1. helyXA1. hely⋅XB2. hely

Ennek az egyensúlynak az elméleti alapja termodinamikai. A csere-reakció standard Gibbs-féle szabadenergia-változása, ΔGrxn∘, diktálja a termodinamikai egyensúlyi állandó, Keq, értékét 9:

Keq=exp(−RTΔGrxn∘)

ahol R az egyetemes gázállandó. A KD csere-együttható közvetlen kapcsolatban áll a Keq-val az egyes helyeken lévő fajták aktivitási együtthatóin (γ) keresztül, amelyek a szilárd oldatban történő nemideális elegyedést veszik figyelembe.18 A kapcsolat a következő:

KD=Keq⋅γA1. hely⋅γB2. helyγA2. hely⋅γB1. hely

Ez a kapcsolat egy kritikus pontot hangsúlyoz, amely a jelen modell egyik sarokkövét képezi: az ioncsere a szilikátásványokban ritkán ideális (γ=1). Az aktivitási együtthatók az összetétel és a hőmérséklet függvényei, ami azt jelenti, hogy a csere termodinamikai hajtóereje maga is a rendszer állapotának nemlineáris függvénye. Ezt a belső nem-idealitást rögzíteni kell ahhoz, hogy a modell fizikailag valósághű legyen.1

Az egyensúlyi állapot nem statikus, hanem mélyrehatóan változik a bolygóköpeny környezeti feltételeivel. A Keq hőmérséklet- és nyomásfüggését a Gibbs-féle szabadenergia definíciójából, ΔG=ΔH−TΔS+PΔV=−RTlnKeq, levezetett alapvető termodinamikai összefüggések adják meg.1 Az ebből adódó függőségek a következők:

(∂T∂lnKeq)P=RT2ΔH∘(van ’t Hoff-egyenlet)

20

(∂P∂lnKeq)T=−RTΔV∘

ahol ΔH∘ a standard entalpiaváltozás és ΔV∘ a standard térfogatváltozás a csere-reakció során. Ezek az egyenletek biztosítják a matematikai keretet az egyensúlyi állapot, KD, parametrizálásához egy bolygóköpeny teljes P-T gradiensén keresztül.

1.2. A Mg²⁺↔Fe²⁺ Csere Ferromágneses Szilikátokban (Olivin és Piroxén)

A magnézium és a vas cseréje a legfontosabb kémiai változó a Föld felső köpenyében, amely meghatározza annak domináns ásványai, az olivin ((Mg,Fe)₂SiO₄) és a piroxén ((Mg,Fe)₂Si₂O₆) fizikai tulajdonságait.21 Ez a csere izovalens (

Mg2+ a Fe2+-ért) és elsősorban intrakristályos (kristályon belüli) folyamatként zajlik, amely e kationok átcsoportosítását jelenti a krisztallográfiailag különböző, M1 és M2 jelű oktaéderes helyek között, miközben a szilikátváz nagyrészt passzív marad.1

Az ortopiroxénben az M1 hely egy viszonylag szabályos oktaéder, míg az M2 hely nagyobb és torzultabb. Ez a szerkezeti különbség kationrendeződéshez vezet, ahol a nagyobb Fe2+ ion egyértelműen az M2 helyet részesíti előnyben.1 A csere-reakció a következőképpen írható fel:

Fe2+(M2)+Mg2+(M1)⇌Fe2+(M1)+Mg2+(M2)

A Mg2+ és Fe2+ elegyedése ezeken a helyeken nem ideális. Termodinamikai vizsgálatok, különösen azok, amelyek Mössbauer-spektroszkópiát alkalmaztak a rácshelyek betöltöttségének meghatározására különböző hőmérsékleteken, kimutatták, hogy az elegyedési viselkedést jól leírja a Guggenheim-féle „egyszerű elegy” vagy egy reguláris oldatmodell.1 Ebben a keretrendszerben a keveredési többlet Gibbs-féle szabadenergiát,

GE, egy W kölcsönhatási paraméterrel veszik figyelembe, így egy komponens aktivitási együtthatója a másik komponens móltörtjével függ össze. Egyetlen helyen (pl. M1) történő keveredésre az aktivitási együtthatók a következők:

RTlnγFeM1=WFeMgM1(XMgM1)2

RTlnγMgM1=WFeMgM1(XFeM1)2

Kulcsfontosságú, hogy a kísérleti munkák kimutatták, hogy a kölcsönhatási paraméterek (W) és a csere-reakció standard szabadenergia-változása (ΔG∘) maguk is hőmérsékletfüggők, gyakran lineárisan változnak az inverz hőmérséklettel (1/T).1 Például az ortopiroxénre vonatkozó vizsgálatok azt mutatják, hogy az M1 helyen történő keveredés kevésbé ideális, mint az M2 helyen, és mindkettő az idealitás felé közelít (W→0) a hőmérséklet növekedésével. Továbbá a kristályos oldat aszimmetrikus, az aszimmetria mértéke alacsonyabb hőmérsékleten növekszik, ami komplexebb polinomiális kifejezéseket tesz szükségessé a keveredési többlet szabadenergiára.1 Ez a nem ideális viselkedés nem egy kisebb korrekció, hanem a rendszer termodinamikájának domináns jellemzője. Egy modell, amely ezt elhanyagolná az ideális keveredés feltételezésével, alapvetően tévesen számítaná ki az egyensúlyi állapotot és következésképpen a rendszer teljes dinamikus fejlődését.

1.3. A Kapcsolt Ca²⁺↔Na⁺ Csere a Plagioklász Földpátokban

A plagioklász földpátok, az albit (NaAlSi₃O₈) és az anortit (CaAl2Si2O8) közötti szilárd oldat, a bolygókérgek kulcsfontosságú alkotóelemei, és jelen lehetnek a legfelső köpenyben is.24 A

Na+ és Ca2+ cseréje ebben a rendszerben alapvetően különbözik és összetettebb, mint a Mg-Fe csere. Ez egy heterovalens csere, amely kapcsolt helyettesítést igényel a kristályrácson belüli töltéssemlegesség fenntartásához 1:

Na+(A-hely)+Si4+(T-hely)⇌Ca2+(A-hely)+Al3+(T-hely)

Ennek a mechanizmusnak mélyreható következményei vannak. A ferromágneses szilikátokban a szilikátváz passzív szerepével ellentétben itt a váz aktív résztvevő. A Na+ cseréje Ca2+-ra egyidejűleg megköveteli az Al3+ és Si4+ átrendeződését a tetraéderes (T) helyeken. Ez a folyamat erős Si-O és Al-O kötések felbontását és újraformálását foglalja magában, ami sokkal magasabb kinetikai gátat és egy eltérő diffúziós mechanizmust feltételez az egyszerű izovalens kationcserékhez képest.1

A plagioklász szilárd oldat kivételes szerkezeti és termodinamikai összetettségéről híres.1 A hőmérséklet csökkenésével az ásvány egy sor fázisátalakuláson megy keresztül (pl. a Cī-ről lī-re való átalakulás), és kiterjedt Al-Si rendeződést mutat a tetraéderes helyeken. Ez a rendezett-rendezetlen állapot erősen függ a hőmérséklettől és az összetételtől is, és jelentős hatással van a keveredési szabadenergiára. Az egyszerű termodinamikai modellek, mint a reguláris oldatmodell, nem megfelelőek. Pontosabb modellezéshez kifinomultabb megközelítésekre, például a Darken-féle kvadratikus formalizmusra (DQF) vagy a Landau-elméletre van szükség az aktivitás-összetétel viszonyok pontos modellezéséhez.30 Ezek a modellek gyakran a különböző szerkezeti állapotokat különálló termodinamikai oldatokként kezelik, saját végtagokkal és keveredési tulajdonságokkal, a fázishatárokon folytonosságot kényszerítve.1 Ez a részletességi szint elengedhetetlen az egyensúlyi állapot pontos meghatározásához, amely a csere dinamikáját szabályozza.

1.4. Egységesített, P-T Függő Megoszlási Modell Köpenykörülményekre

Egy bolygóköpenyre alkalmazható analitikai modell megalkotásához az egyes ásványrendszerekre vonatkozó termodinamikai leírásokat szintetizálni és kiterjeszteni kell a releváns nyomások és hőmérsékletek teljes tartományára. A köpeny egy dinamikus környezet, ahol a P és T körülbelül 500 K-től a kéreg-köpeny határon több mint 4000 K-ig és 136 GPa-ig terjed a mag-köpeny határon.4

Ezt a hatalmas P-T teret jelentős szilárdtest-fázisátalakulások tagolják, amelyek kapcsolóként működnek a megoszlási viselkedésben.33 Körülbelül 410 km és 660 km mélységben a felső köpeny elsődleges ásványai sűrűbb, nagynyomású polimorfokká alakulnak át.5 Az olivin először wadsleyitté, majd ringwoodittá alakul.34 Nagyobb mélységben a ringwoodit és a piroxének lebomlanak, és bridgmanit (

(Mg,Fe)SiO3-perovszkit) és ferroperiklász ((Mg,Fe)O) együttesét alkotják.36 Ezen fázisok mindegyike egyedi kristályszerkezettel és eltérő kationhelyekkel rendelkezik. Következésképpen a

Fe2+ és a Mg2+ megoszlási preferenciája az egymás mellett létező ásványok között, valamint az egyetlen ásványon belüli intrakristályos helypreferencia hirtelen megváltozhat ezeken a szeizmikus diszkontinuitásokon keresztül.9

Ezért a csere-együttható, KD(T,P), egységes modellje nem lehet egyetlen, sima függvény. A nyomás (vagy mélység) szakaszonként definiált függvényeként kell megfogalmazni, ahol minden szakasz egy-egy specifikus ásványtani stabilitási mezőnek felel meg. Minden mezőn belül a hőmérséklet- és nyomásfüggést az adott ásványfázisra vagy együttesre vonatkozó csere-reakció ΔHo, ΔSo és ΔV∘ értékei szabályozzák. Ezen függvények paramétereit átfogó, önkonzisztens termodinamikai adatbázisokból kell venni, amelyeket évtizedekig tartó nagynyomású kísérleti és elméleti ásványfizikai kutatásokból állítottak össze.13 A következő táblázatok összefoglalják egy ilyen modell felépítéséhez szükséges bemeneti adatokat.

1. Táblázat: Termodinamikai adatok a kationcseréhez köpenyásványokban. Ez a táblázat összesíti az egyensúlyi modellhez szükséges alapvető termodinamikai paramétereket. Ezek az adatok, amelyeket kísérleti kőzettanból és ásványfizikából származtatnak, képezik a KD(T,P) kiszámításának alapját.

Csere-reakció Ásványfázis ΔH∘ (kJ/mol) ΔS∘ (J/mol·K) ΔV∘ (cm³/mol) Forrás(ok)

Mg2+(M1)+Fe2+(M2)⇌Fe2+(M1)+Mg2+(M2) Ortopiroxén Érték Érték Érték 4

Mg2+(M1)+Fe2+(M2)⇌Fe2+(M1)+Mg2+(M2) Olivin Érték Érték Érték 17

NaSi ↔ CaAl Plagioklász Érték Érték Érték 6

(bejegyzések a wadsleyitre, ringwooditra, bridgmanitra stb.)

Megjegyzés: A konkrét értékek helyőrzők, amelyeket az idézett irodalomból és termodinamikai adatbázisokból kell feltölteni.

2. Táblázat: A földi köpeny reprezentatív fizikai körülményei és ásványtana. Ez a táblázat biztosítja a fizikai kontextust, meghatározva azokat a P-T-ásványtan útpontokat, amelyekre a termodinamikai és kinetikai modelleknek érvényesnek kell lenniük.

Köpenyrégió Mélység (km) Nyomás (GPa) Hőmérséklet (K) Domináns ásványegyüttes Forrás(ok)

Felső köpeny ~100 ~3 ~1500 Olivin + Ortopiroxén + Klinopiroxén + Gránát 10

Átmeneti zóna (teteje) 410 ~13.5 ~1800 Wadsleyit + Majoritos gránát 10

Átmeneti zóna (alja) 660 ~24 ~1900 Ringwoodit + Majoritos gránát 10

Alsó köpeny (közepe) 1500 ~60 ~2500 Bridgmanit + Ferroperiklász + Ca-perovszkit 9

Mag-köpeny határ (D") 2890 ~136 ~4000 Bridgmanit (vagy Poszt-perovszkit) + Ferroperiklász 8

2. Szekció: Egy Nemlineáris Elektro-diffúzív Transzportmodell Felállítása

Az egyensúlyi állapotot meghatározó termodinamikai alapokra építve ez a fejezet kidolgozza a modell dinamikus magját. Levezetünk egy időfüggő, nemlineáris parciális differenciálegyenlet-rendszert (PDE), amely leírja azt a kinetikai útvonalat, amelyen a köpenyásványok megközelítik ezt az egyensúlyt. A központi újítás a Nernst-Planck-Poisson (NPP) formalizmus alkalmazása a szilárdtest-intrakristályos diffúzióra. Ez a megközelítés önkonzisztensen kapcsolja össze a kémiai potenciálgradiensek által hajtott diffúziós fluxust az eltérő ionos mobilitásból adódó belső elektromos mezők által hajtott elektromigrációs fluxussal, mindezt egy reaktív tag hatása alatt, amelyet a helyi termodinamikai egyensúlytól való eltérés szabályoz.1

2.1. A Meghatározó Fizika: A Nernst-Planck-Poisson (NPP) Rendszer

Egy egyszerű rendszerben, ahol két cserélődő ion mobilitása azonos, transzportjukat Fick diffúziós törvényei írhatják le.40 Valós ásványrendszerekben azonban a cserélődő kationok (pl.

DMg2+ vs. DFe2+) diffúziós együtthatói általában nem egyenlőek. Ennek az egyenlőtlenségnek mélyreható következménye van: ahogy a gyorsabb ion a koncentrációgradiensének mentén halad, átmenetileg megelőzi a lassabb, ellenirányban diffundáló iont, ami egy csekély, de jelentős helyi töltésszétváláshoz vezet. Ez a töltésegyensúly-hiány egy belső elektromos mezőt, gyakran diffúziós potenciálnak nevezett mezőt hoz létre.42 Ez a mező elektrosztatikus erőt fejt ki minden töltött részecskére, lassítva a gyorsabb ion mozgását és gyorsítva a lassabbét, ami egy erőteljes nemlineáris visszacsatolás, amely a helyi elektroneutralitás fenntartására törekszik.44 Egy kizárólag Fick-féle diffúzión alapuló modell elhanyagolná ezt a kulcsfontosságú elektrosztatikus csatolást, ami fizikailag megalapozatlan töltésfelhalmozódás előrejelzéséhez és helytelen nettó fluxusarányokhoz vezetne.

A Nernst-Planck egyenleteket kifejezetten ennek a kettős transzportmódnak a megragadására tervezték. Az i-edik ionfajta fluxusát (Ji) egy diffúziós tag (a ci koncentrációgradiens által hajtott) és egy elektromigrációs vagy sodródási tag (a Φ elektromos potenciálgradiens által hajtott) összegeként írják le 39:

Ji=−DiffuˊzioˊDi∇ci−ElektromigraˊcioˊRTDiziFci∇Φ

ahol Di a diffúziós együttható, ci a koncentráció, zi az ion vegyértéke, F a Faraday-állandó, R az egyetemes gázállandó, T az abszolút hőmérséklet, és Φ a helyi elektromos potenciál. Az elektromos potenciál, Φ, nem külső mező, hanem egy belső, amelyet maguk az ionok hoznak létre. Ezt a Poisson-egyenlet szabályozza, amely az elektromos mező divergenciáját a helyi nettó töltéssűrűséggel kapcsolja össze:

∇2Φ=−ϵρcharge=−ϵFi∑zici

ahol ϵ a közeg (a szilikátásvány) dielektromos permittivitása.45 Az NPP rendszer ezen Nernst-Planck és Poisson egyenletek kapcsolt halmaza. Bár a vizes elektrokémia és a biofizika területén szabványnak számít 2, a szilárdtest-diffúzióra való alkalmazása a köpenyásványokban gondos újraértelmezést igényel. Itt a koncentrációk móltörteket jelentenek specifikus kristály-alrácsokon, a diffúziós együtthatók erős Arrhenius-típusú hőmérsékletfüggéssel rendelkező szilárdtest-nyomjelző diffúziós együtthatók, a permittivitás pedig a gazda szilikáté, amely maga is függ a nyomástól és a hőmérséklettől. Ez az adaptáció az NPP rendszert a köpeny extrém körülményeihez kötődő további komplexitási rétegekkel ruházza fel.

2.2. A Kapcsolt Parciális Differenciálegyenletek Levezetése a Köpenybeli Kationcserére

Az egyes ionfajták koncentrációjának időbeli alakulását egy tömegmérleg- vagy kontinuitási egyenlet írja le. Ez az egyenlet kimondja, hogy egy ponton a koncentráció változásának sebessége megegyezik az adott pontba irányuló nettó fluxussal, plusz a fajta helyi keletkezésével vagy megsemmisülésével a kémiai reakció (Ri) során 39:

∂t∂ci=−∇⋅Ji+Ri

A fluxusra vonatkozó Nernst-Planck kifejezés behelyettesítése adja meg az egyes fajtákra vonatkozó meghatározó parciális differenciálegyenletet (PDE).49 A modell kulcsfontosságú újítása a reakciótag,

Ri, megfogalmazásában rejlik. A rendszert valódi reaktív transzport problémaként kezeljük, ahol a helyi termodinamikai egyensúlytól való eltérés, amelyet az 1. Szekcióból származó KD(T,P) határoz meg, folyamatos hajtóerőt biztosít a csere-reakció számára térben és időben minden ponton.51

2.2.1. Modellegyenlet az izovalens cserére (Mg2+↔Fe2+)

A Mg2+ és Fe2+ (mindkettő z=2 vegyértékű) közötti izovalens cserére két helyen (pl. M1 és M2 a piroxénben), legyen CFe és CMg a Fe és Mg teljes koncentrációja. Az egyes helyeken lévő koncentrációk ezen teljes koncentrációk függvényei. A reakciótag a nettó átalakulási sebességet jelenti egyik helyről a másikra. A kapcsolt PDE-k teljes rendszere:

∂t∂CFe=∇⋅(DFe∇CFe+RTDFezFeFCFe∇Φ)

∂t∂CMg=∇⋅(DMg∇CMg+RTDMgzMgFCMg∇Φ)

∇2Φ=−ϵF(CFe+CMg−Ctotal)

Ebben a megfogalmazásban az intrakristályos diffúzióra a reakciótag implicit módon van kezelve. A KD(T,P) által diktált egyensúlyi állapot nem forrás/nyelő tag a hagyományos értelemben, hanem az a feltétel, amely felé a rendszer fejlődik. A diffúziót hajtó kémiai potenciálgradiensek eleve tartalmazzák a helypreferenciára és a nem ideális keveredésre vonatkozó termodinamikai információkat. A rendszer akkor éri el az egyensúlyt, amikor minden fajta elektrokémiai potenciálja egységes minden helyen, ekkor a nettó fluxus nullává válik. A Ctotal tag a kétértékű kationok fix teljes koncentrációját jelenti az oktaéderes helyeken, biztosítva a töltésegyensúlyt a Poisson-egyenletben.

2.2.2. Modell adaptációja a heterovalens cserére (Ca2+↔Na+)

A heterovalens Ca2+ (z=2) ↔ Na+ (z=1) csere a plagioklászban nagyobb kihívást jelent a Si/Al váz kapcsolt átrendeződése miatt. Egy teljes modellnek négy fajtát (Na, Ca, Al, Si) és azok komplex kölcsönhatásait kellene követnie. Egy analitikai megoldás levezetése érdekében egy egyszerűsített, de fizikailag értelmes modellt fogalmazunk meg. A cserét Na és Ca2+ ellen-diffúziójaként kezeljük az A-helyen, amelyet effektív diffúziós együtthatók (DNa,eff, DCa,eff) szabályoznak, amelyek implicit módon figyelembe veszik a kapcsolt váz-átrendeződés kinetikai gátját. A meghatározó egyenletek a következők:

∂t∂CNa=∇⋅(DNa,eff∇CNa+RTDNa,effzNaFCNa∇Φ)

∂t∂CCa=∇⋅(DCa,eff∇CCa+RTDCa,effzCaFCCa∇Φ)

∇2Φ=−ϵF(CNa+2CCa−Ccharge)

Itt Ccharge az anionos váz kiegyensúlyozásához szükséges teljes pozitív töltés, és a különböző vegyértékek (z=1 a Na-ra, z=2 a Ca-ra) alapvetően megkülönböztetik az elektrosztatikus csatolási tagot a Mg-Fe rendszertől.

2.3. Dimenziótlanítás és Peremfeltételek

Az alapvető fizika feltárása és a matematikai elemzés egyszerűsítése érdekében a meghatározó PDE rendszert dimenziótlanítjuk.52 Bevezetünk dimenziótlan változókat: koncentráció

c^=c/cref, pozíció x^=x/Lref, idő t^=t/τref, és potenciál ϕ^=ϕF/(RT). A referencia időskálát az egyik fajta diffúziós idejeként választva (pl. τref=Lref2/DMg), a folyamat egy sor kulcsfontosságú dimenziótlan paramétert eredményez, amelyek a rendszer viselkedését szabályozzák.

3. Táblázat: Meghatározó dimenziótlan paraméterek és fizikai jelentőségük. Ez a táblázat definiálja a modellegyenletekből származó kulcsfontosságú dimenziótlan csoportokat, tömör összefoglalást nyújtva a versengő fizikai folyamatokról.

Paraméter Szimbólum Definíció Fizikai jelentőség

Diffúziós arány δ DMgDFe

Az ionos mobilitások aránya; az =1 érték a diffúziós potenciál forrása.

Debye-szám κ 2F2crefϵRTLref

A karakterisztikus hossz (pl. szemcseméret) és a Debye-hossz aránya; az elektrosztatikus hatások erősségét szabályozza.51

Damköhler-szám Da τreactionτdiffusion=DMgkfcrefLref2

A diffúziós időskála és a reakcióidőskála aránya; meghatározza, hogy a folyamat diffúzió-limitált (Da≫1) vagy reakció-limitált (Da≪1).54

A peremfeltételeket egy reprezentatív kontrolltérfogatra, tipikusan egy R sugarú gömb alakú ásványszemcsére definiáljuk.

1. A szemcse közepén (r=0): A gömbszimmetria miatt minden fajtára zéró fluxus feltételt alkalmazunk: ∇ci=0 és ∇ϕ=0.

2. A szemcsehatáron (r=R): Az itteni feltételek a szemcse környezetével való kölcsönhatását reprezentálják. Egy P-T változás után belsőleg újraegyensúlyozódó rendszer esetében a zéró fluxus peremfeltétel a megfelelő: Ji⋅n^=0. Egy külső folyadékkal vagy szomszédos ásványokkal reagáló szemcse (metaszomatózis) esetében Dirichlet- (fix koncentráció) vagy Robin-típusú (fix fluxus) peremfeltételt kellene használni.55

Az egész rendszer egy nagyobb köpenykonvekciós keretrendszerbe van ágyazva, amely az időben változó P-T útvonalat, P(t) és T(t), biztosítja. Ezek a külső feltételek dinamikusan módosítják a modell együtthatóit (Di, ϵ) és a termodinamikai célt (KD), ami a problémát egy merev, több skálájú kihívássá teszi, ahol az ásványi szintű kinetika a bolygószintű dinamikához kapcsolódik.54

3. Szekció: Egzakt Analitikai Megoldás a Homotópia Analízis Módszerrel (HAM)

A 2. Szekcióban levezetett kapcsolt, nemlineáris PDE-rendszer magában foglalja a reaktív elektro-diffúzió komplex fizikáját, de hagyományos, közvetlen analitikai módszerekkel megoldhatatlan. A felhasználó kérése egy „egzakt analitikai megoldásra”, „tisztán matematikai eszközökkel”, egy olyan kifinomult megközelítést tesz szükségessé, amely képes kezelni az erős nemlinearitást anélkül, hogy tisztán numerikus diszkretizációhoz folyamodna. Ez a fejezet indokolja a Homotópia Analízis Módszer (HAM) kiválasztását, és lépésről lépésre bemutatja az analitikai sorfejtéses megoldás levezetését. A HAM szigorú matematikai keretet biztosít a megoldás konvergens sorreprezentációjának megszerzéséhez, amely gyakorlati szempontból egzakt analitikai formaként szolgál.4

3.1. A Homotópia Analízis Módszer (HAM) Választásának Indoklása

A nemlineáris PDE-k megoldására rendelkezésre álló módszerek rövid áttekintése a tisztán numerikustól a tisztán analitikai technikákig terjedő spektrumot tár fel. A numerikus módszerek, mint a Véges Differenciák Módszere vagy a Végeselem Módszer, hatékonyak specifikus esetekben, de nem adják meg azt az általános funkcionális megoldást vagy mély fizikai betekintést, amelyet egy analitikai forma nyújt.57

A félig analitikai technikák közül olyan módszereket fejlesztettek ki a nemlineáris problémák kezelésére, mint az Adomian-féle dekompozíciós módszer (ADM).58 Az ADM-nek azonban ismert korlátai vannak. Konvergenciája nem garantált minden probléma esetén, és néha nem fizikai „zajtagokat” eredményezhet az inhomogén egyenletek megoldásában, ami bonyolíthatja a fizikai értelmezést.58

A Homotópia Analízis Módszer (HAM), amelyet Liao Shijun javasolt, felülmúlja ezeket a korlátokat, és egyedülállóan alkalmas erre az alapvető elméleti problémára.4 Fő előnyei a következők:

1. Függetlenség a kis paraméterektől: A perturbációs módszerekkel ellentétben a HAM nem támaszkodik egy kis fizikai paraméter létezésére a meghatározó egyenletekben. Ezért egyaránt alkalmazható gyengén és erősen nemlineáris rendszerekre is, ami kulcsfontosságú a modellünk szempontjából, ahol az elektrosztatikus hatások vagy a reakciókinetika dominánsak lehetnek.4

2. Garantált konvergencia: A HAM bevezet egy nem fizikai, mesterséges konvergencia-kontroll paramétert, ℏ-t. Ez a paraméter egyszerű és szigorú mechanizmust biztosít a kapott megoldási sor konvergenciájának szabályozására és biztosítására. A ℏ érvényes tartományának létezése grafikusan, egy „ℏ-görbe” segítségével igazolható, ami egyedülálló jellemző az analitikai közelítő technikák között, és matematikai bizonyosságot nyújt a megoldás érvényességében.61

3. Rugalmasság a megfogalmazásban: A HAM teljes szabadságot kínál a megoldás kezdeti becslésének és, ami a legfontosabb, a kiegészítő lineáris operátor kiválasztásában. Ez a rugalmasság nem csupán matematikai kényelem; stratégiailag felhasználható arra, hogy a megoldási sort mély fizikai jelentéssel ruházza fel.4

Tekintettel egy új, robusztus elméleti keretrendszer kifejlesztésének céljára, a HAM által kínált matematikai szigor és garantált konvergencia teszi azt a legkiválóbb választássá. Nem csupán egy megoldást, hanem egy megbízható és értelmezhető megoldást nyújt.

3.2. A HAM Sorfejtéses Megoldás Lépésről-Lépésre Történő Levezetése

A HAM alapötlete egy folytonos deformáció, vagy „homotópia” megkonstruálása, amely összeköt egy egyszerű, megoldható kezdeti problémát a megoldani kívánt komplex, nemlineáris problémával. Jelöljük a 2. Szekcióból származó nemlineáris PDE-rendszert operátoros formában: N[ci(x,t),Φ(x,t)]=0.

3.2.1. A homotópia megkonstruálása

Bevezetünk egy beágyazási paramétert, $q \in $, és megkonstruáljuk a „nulladrendű deformációs egyenletet” 62:

(1−q)L[ϕi(x,t;q)−ci,0(x,t)]=qℏH(x,t)N[ϕi(x,t;q),Ψ(x,t;q)]

Itt:

• ϕi(x,t;q) és Ψ(x,t;q) olyan függvények, amelyek egy kezdeti becsléstől (q=0) az egzakt megoldásig (q=1) változnak.

• L egy kiegészítő lineáris operátor, amelyet szabadon választhatunk.57 Stratégiai választás a lineáris diffúziós operátor,

L[u]=∂u/∂t−D∇2u.

• ci,0(x,t) a megoldás kezdeti becslése, amelynek teljesítenie kell a perem- és kezdeti feltételeket.

• ℏ=0 a konvergencia-kontroll paraméter.

• H(x,t) egy nem-zéró kiegészítő függvény, amelyet az egyszerűség kedvéért gyakran 1-nek választanak.

Amikor q=0, az egyenlet L[ϕi−ci,0]=0 lesz, egy egyszerű lineáris egyenlet, amelynek megoldása ϕi(x,t;0)=ci,0(x,t). Amikor q=1, az egyenlet N[ϕi(x,t;1),Ψ(x,t;1)]=0 lesz, ami az eredeti nemlineáris problémánk, tehát ϕi(x,t;1)=ci(x,t). Ahogy q 0-tól 1-ig változik, a ϕi megoldás simán deformálódik a kezdeti becsléstől az egzakt megoldásig.

3.2.2. A magasabbrendű deformációs egyenletek levezetése

Feltételezve, hogy a ϕi és Ψ megoldások analitikusak q-ban, Taylor-sorként fejthetjük ki őket q=0 körül:

ϕi(x,t;q)=ci,0(x,t)+m=1∑∞ci,m(x,t)qmaholci,m=m!1∂qm∂mϕiq=0

Ezt a sort behelyettesítve a nulladrendű deformációs egyenletbe és az azonos hatványú q-k együtthatóit egyenlővé téve, az egyetlen nemlineáris PDE-t lineáris PDE-k végtelen sorozatává alakítjuk át. Ezt úgy érjük el, hogy a nulladrendű egyenletet m-szer differenciáljuk q szerint, elosztjuk m!-sal, majd q=0-t helyettesítünk be. Az eredményül kapott „m-ed rendű deformációs egyenlet” m≥1 esetén 61:

L[ci,m(x,t)−χmci,m−1(x,t)]=ℏH(x,t)Ri,m(cm−1,Ψm−1)

ahol χm=1, ha m>1 és χ1=0. Az Ri,m tag egy, az N nemlineáris operátorból származtatott függvény, és csak az előző rendek megoldásaitól (ci,0,...,ci,m−1) függ.

3.2.3. A megoldás komponenseinek megoldása

Ez a rekurzív rendszer a HAM ereje. A probléma most egy sor lineáris, inhomogén PDE megoldása. m=1 esetén:

L[ci,1(x,t)]=ℏH(x,t)Ri,1(c0,Ψ0)

Tekintettel az L stratégiai választására, mint diffúziós operátor, és egy ismert ci,0 kezdeti becslésre, a jobb oldal egy ismert függvény. Ez a lineáris PDE analitikusan megoldható standard technikákkal, mint a Green-függvények vagy a változók szétválasztása.64 Miután

ci,1-et megtaláltuk, felhasználjuk az m=2 egyenlet jobb oldalának megkonstruálásához, amelyet aztán megoldunk ci,2-re, és így tovább.

A végső analitikai megoldás ezen komponensek összege:

ci(x,t)=m=0∑∞ci,m(x,t)

Hasonló eljárást követünk párhuzamosan az elektromos potenciálra, Φ(x,t)-re is. A gyakorlatban a sort egy kellően magas M rendnél csonkoljuk a kívánt pontosság elérése érdekében.

3.3. A Megoldás Elemzése: Konvergencia és Fizikai Értelmezés

A HAM megoldás érvényessége a sor konvergenciáján múlik. A konvergencia-kontroll paraméter, ℏ, mechanizmust biztosít ennek garantálására. Egy adott közelítési rend esetén egy kulcsfontosságú fizikai mennyiséget (pl. a koncentrációt a szemcse közepén, ci(0,t)) kiszámítunk és ábrázolunk ℏ függvényében. Ez a diagram, az úgynevezett ℏ-görbe, tipikusan egy vízszintes platót mutat. Az e platónak megfelelő ℏ értékek tartománya a konvergencia érvényes régiója. Egy ℏ érték kiválasztása ebből a régióból garantálja, hogy a sorfejtéses megoldás a probléma valódi megoldásához konvergál.66

A HAM sorfejtéses megoldás szerkezete mély fizikai betekintést nyújt. A kiegészítő lineáris operátor, L, legegyszerűbb fizikai folyamatként (azaz tiszta diffúzióként) való megválasztásával a sor tagjai a komplex problémát kölcsönható fizikai hatások hierarchiájára bontják:

• c0(x,t): A kezdeti becslés, amely a rendszer kiindulási állapotát reprezentálja.

• c1(x,t): Az elsőrendű korrekció, amely a rendszer kezdeti válaszát képviseli az elektromigráció és a kémiai reakció együttes hatásaira.

• c2(x,t): A másodrendű korrekció, amely az elsőrendű változások visszacsatolását reprezentálja a nemlineáris rendszerbe, és így tovább.

A sor tehát a fizika egy narratíváját adja, megmutatva, hogyan módosítják az egyszerűbb folyamatokat fokozatosan a komplexebb, nemlineáris csatolások.

Végül, ennek az egzakt analitikai megoldásnak a létezése, még egy egyszerűsített geometriára is, felbecsülhetetlen erőforrást jelent a geokémiai modellezés tágabb területén.3 A reaktív transzportmodellek túlnyomó többsége komplex numerikus kód.47 Ezen kódok pontosságának ellenőrzése komoly kihívást jelent. A mi HAM megoldásunk egzakt „alapigazságként” (ground truth) szolgál, amelyhez képest ezek a numerikus kódok tesztelhetők, kalibrálhatók és validálhatók, ezáltal növelve a mélyföldi folyamatok numerikus szimulációinak megbízhatóságát az egész tudományágban.68

4. Szekció: Az Analitikai Modell Geokémiai és Geodinamikai Következményei

A 3. Szekcióban levezetett analitikai megoldás nem csupán matematikai gyakorlat; egy új, hatékony eszközt biztosít az ásványi szintű fizika bolygószintű folyamatokra gyakorolt mélyreható befolyásának vizsgálatára. Ez a fejezet visszakapcsolja a szigorú matematikai keretrendszert a bolygódifferenciálódás és a köpenykonvekció motiváló kérdéseihez. A megoldás realisztikus geodinamikai forgatókönyvekre való alkalmazásával számszerűsíthetjük a kémiai folyamatok időskáláit, megérthetjük a kémiai szegregáció mechanizmusait, és feltárhatjuk az összetétel és a köpenydinamika közötti kritikus visszacsatolási hurkot.

4.1. Az Ásványi Összetételek Tér-időbeli Fejlődése egy Konvektáló Köpenyben

A Föld köpenye állandó, bár lassú konvekcióban van, az anyag forró köpenycsóvákban emelkedik fel és hideg, alábukó lemezekben süllyed le.69 Ahogy egy kőzetcsomag ezeken az útvonalakon halad, drámai nyomás- és hőmérsékletváltozásokon megy keresztül, millióktól százmillió évekig terjedő időskálán. Az analitikai megoldás lehetővé teszi számunkra, hogy nyomon kövessük az ásványok összetételének alakulását egy ilyen csomagban, amint azok megpróbálnak újraegyensúlyozódni a változó környezeti feltételekhez.

Két végletes forgatókönyvet modellezhetünk geofizikai modellekből származó P-T útvonalak segítségével 4:

1. Köpenycsóva felemelkedése: Egy mélyköpenyből származó kőzetcsomag gyorsan felemelkedik, gyors nyomáscsökkenést és viszonylag kis (adiabatikus) hőmérsékletcsökkenést tapasztalva.59

2. Lemez lesüllyedése: Egy hideg óceáni litoszféra-darab lassan lesüllyed a köpenybe, fokozatos nyomásnövekedést és lassú hőmérséklet-emelkedést tapasztalva, ahogy a környező köpenyből származó hővezetés felmelegíti.72

Ezen útvonalak minden időlépésénél két mennyiséget számíthatunk ki:

• Az egyensúlyi összetételt, Ceq(P(t),T(t)), amelyet az 1. Szekció termodinamikai modellje határoz meg. Ez képviseli a „cél” összetételt, amelyet az ásvány az adott mélységben és hőmérsékleten igyekszik elérni.

• A tényleges összetételt, C(t), amelyet a 3. Szekció időfüggő analitikai megoldásával számítunk ki, amely nyomon követi a kinetikai fejlődést az előző állapotból.

A C(t) és a Ceq(P(t),T(t)) összehasonlításával közvetlenül számszerűsíthetjük az ásványokban megőrzött kémiai egyensúlyhiány mértékét. Ez teszteli azt a régóta fennálló hipotézist, hogy a felszínen kitörő ásványok (pl. óceáni szigeti bazaltokban) kémiai jelzéseket hordozhatnak mélyköpeny-forrásrégióikból.25 A modell előre jelezheti, hogy milyen körülmények között (pl. gyors felemelkedés, nagy szemcseméret, alacsony hőmérséklet) valószínűsíthető az ilyen jelzések megőrződése.

Ennek az elemzésnek a kulcsfontosságú kimenete a kémiai egyensúly beállásának karakterisztikus időskálája, τcsere. Ez az időskála, amely erősen hőmérsékletfüggő a diffúziós együtthatók Arrhenius-jellegű természete miatt, összehasonlítható a köpenytranszport időskálájával, τkonvekcioˊ. Ez az összehasonlítás egy dimenziótlan „Geokémiai Memória Számban”, Gm=τkonvekcioˊ/τcsere, foglalható össze.73

• Ha Gm≫1, az egyensúly beállása sokkal gyorsabb, mint a transzport. Az ásványok folyamatosan újraegyensúlyozódnak a környezetükkel, hatékonyan törölve mély forrásuk emlékét.

• Ha Gm≪1, a transzport sokkal gyorsabb, mint az egyensúly beállása. Az ásványok „befagynak” forrásrégiójuk összetételével, tökéletes kémiai nyomjelzőként működve.

Modellünk lehetővé teszi a Gm kiszámítását a köpeny egészében, létrehozva egy prediktív térképet arról, hogy hol őrződik meg és hol semmisül meg a kémiai memória. Ez egy hatékony, kvantitatív keretet biztosít a köpenyből származó kőzetek komplex izotóp- és nyomelem-geokémiájának értelmezéséhez.33

4.2. Egy Dinamikus Modell a Szilárdtest-alapú Bolygódifferenciálódásra

A bolygódifferenciálódás – egy bolygó kémiailag elkülönülő rétegekre, mint a mag, köpeny és kéreg, való szétválása – gyakran elsődleges, az akkréció során bekövetkező folyamatnak tekintik.74 Modellünk azonban azt mutatja, hogy a differenciálódás egy folyamatos, szilárdtest-fizika által hajtott folyamat lehet egy konvektáló köpenyben. A köpenykonvekció bolygóméretű kémiai finomítóként működik.54

A csere-együttható, KD, P-T függése ennek a finomítónak a motorja. Például az ásványfizika azt jósolja, hogy az alsó köpeny extrém nyomásán a Fe2+ előnyösen a sűrű bridgmanit fázisba kerül a Mg2+-hoz képest.75 Ahogy a köpenykonvekció az anyagot ezen a régión keresztül ciklikusan mozgatja, nettó tendencia mutatkozik arra, hogy a leáramló áramlatok a vasat megkötő bridgmanitban dúsulnak, míg a feláramló áramlatokból a vas kivonódik, így azok fokozatosan magnéziumban gazdagabbá válnak. Geológiai időskálán ez a folyamatos szilárdtest-csere a vas nettó lefelé irányuló és a magnézium felfelé irányuló fluxusához vezethet, hozzájárulva a köpenyben ma megfigyelhető nagyléptékű kémiai rétegződéshez és heterogenitáshoz.11

A modell felhasználható a köpenybe bekerült különböző kémiai heterogenitások, például az alábukott óceáni kéreg (amely bazaltos összetételű) és az alatta lévő kiürült litoszféra (harzburgit) sorsának tanulmányozására is.25 Ahogy egy bazaltos eklogit lemez lesüllyed és felmelegszik, ioncsere következik be a környező peridotitos köpennyel való határfelületén.76 Analitikai modellünk szimulálhatja ezt a folyamatot, előre jelezve a lemez éles kémiai gradienseinek kisimulásának sebességét, és azt, hogy a kémiai anomáliát hogyan módosítja és szórja szét a diffúzió és a mechanikai keveredés kombinációja. Ez közvetlenül foglalkozik a modern geokémia egyik központi kérdésével: meddig maradhatnak fenn az alábukott kémiai heterogenitások egy erőteljesen konvektáló köpenyben?

4.3. Visszacsatolás a Kationcsere és a Köpenydinamika Között

A modell legmélyebb következménye az, hogy képes bezárni a kört az ásványi szintű kémia és a bolygószintű geodinamika között. A folyamatok nem egyirányúak; létezik egy kritikus visszacsatolási mechanizmus. A köpenyásványok összetétele, különösen a Fe/Mg arány, elsőrendűen szabályozza a köpenykonvekciót meghatározó két fizikai tulajdonságot: a sűrűséget (ρ) és a viszkozitást (η).

• Sűrűség: A vas sűrűbb, mint a magnézium. Ezért egy olyan ásvány, mint az olivin vagy a bridgmanit Fe/Mg arányának növekedése közvetlenül növeli annak sűrűségét.

• Viszkozitás: A köpenyásványok reológiája szintén érzékeny az összetételre, a vastartalom befolyásolja a kristályhibák koncentrációját, amelyek a kúszást és a viszkózus áramlást szabályozzák.6

Analitikai modellünk biztosítja a kulcsfontosságú kapcsolatot: előrejelzi a fejlődő ásványi összetételt egy kőzetcsomag P-T történetének függvényében. Ezt az összetételi kimenetet azután be lehet táplálni a megállapított ásványfizikai állapotegyenletekbe a megfelelő sűrűség és viszkozitás kiszámításához.78

Ez egy teljes visszacsatolási hurkot hoz létre:

1. A Köpenykonvekció beállítja a P-T útvonalat egy kőzetcsomag számára.

2. A P-T útvonal diktálja az időben változó egyensúlyi célt, KD(T,P).

3. Analitikai Modellünk megoldja a kinetikai egyenleteket a tényleges Összetétel(P,T,t) meghatározásához.

4. Az Összetétel meghatározza a kőzetcsomag Sűrűségét (ρ) és Viszkozitását (η) az ásványfizikai adatokon keresztül.

5. A Sűrűség és Viszkozitás változásai módosítják a felhajtóerőket és az áramlással szembeni ellenállást, ami viszont módosítja a Köpenykonvekció mintázatát és hevességét, újraindítva a ciklust.

Ez egy mély összekapcsolódást tár fel, ahol az angström-milliméter skálán zajló kémiai folyamatok több ezer kilométeres skálán befolyásolhatják a folyadékdinamikát. Jelenleg a nagyléptékű numerikus geodinamikai modellek gyakran egyszerűsített feltételezéseket alkalmaznak, például az összetételt passzív nyomjelzőként kezelik, vagy túlságosan leegyszerűsített parametrizációkat használnak az anyagtulajdonságokra.43 Analitikai modellünk egy fizikailag szigorú „rácson belüli” (sub-grid) modul tervrajzát adja, amelyet be lehet építeni ezekbe a nagyléptékű kódokba. Azzal, hogy számításilag hatékony módot biztosít a kőzetösszetétel és az abból következő fizikai tulajdonságok dinamikus frissítésére, ez a munka segíthet áthidalni az ásványfizika és a geodinamika tudományágai közötti régóta fennálló szakadékot, utat nyitva a Föld-típusú bolygók termikus és kémiai evolúciójának valósághűbb és önkonzisztensebb szimulációi felé.14

5. Szekció: Következtetések és Jövőbeli Irányok

5.1. A Főbb Eredmények Összefoglalása

Ez a jelentés sikeresen kezelt egy alapvető, megoldatlan problémát az elméleti geokémiában azáltal, hogy egy átfogó analitikai modellt fejlesztett ki az ioncsere dinamikájára a bolygóköpeny ásványaiban. A munka elsődleges eredményei és legfontosabb megállapításai háromfélék:

1. Egy Újszerű Reaktív-Elektrodiffúzív Modell Megfogalmazása: Első elvekből kiindulva egy szigorú matematikai keretrendszert dolgoztunk ki, amely az intrakristályos kationcserét kapcsolt, nemlineáris parciális differenciálegyenletek rendszereként írja le. A Nernst-Planck-Poisson formalizmus integrálásával a szilikát szilárd oldatokra vonatkozó fejlett, nemideális termodinamikai modellekkel, a modell önkonzisztensen ragadja meg az ionos diffúzió, a belső elektromos mezők által hajtott elektromigráció és a termodinamikai hajtóerő kölcsönhatását, amely explicit módon függ a bolygóköpeny változó nyomás- és hőmérsékleti viszonyaitól.

2. Egy Egzakt Analitikai Megoldás Levezetése: A központi matematikai kihívást a Homotópia Analízis Módszer (HAM) alkalmazásával győztük le. Ez a robusztus technika egy konvergens, analitikai sorfejtéses megoldást eredményezett a meghatározó nemlineáris rendszerre. A HAM keretrendszer matematikai garanciát nyújt a konvergenciára, és a sorfejtéses megoldás szerkezete lehetővé teszi a komplex fizika hierarchikus lebontását, mély betekintést nyújtva a különböző transzport- és reakciómechanizmusok közötti csatolásba. Ez az analitikai megoldás új elméleti eszközként szolgál, és alapvető referenciát (benchmark) biztosít a komplex numerikus kódok ellenőrzéséhez.

3. Geokémiai és Geodinamikai Következmények Számszerűsítése: Az analitikai modellt alkalmaztuk, hogy bemutassuk annak erejét az ásványi szintű fizika és a bolygószintű dinamika közötti híd megteremtésében. Kvantitatív módszert biztosít a kémiai egyensúly beállásának karakterisztikus időskáláinak kiszámítására, lehetővé téve annak előrejelzését, hogy egy konvektáló köpenyben hol és mikor őrződik meg vagy törlődik a forrásrégiók kémiai „memóriája”. A modell a bolygódifferenciálódást egy folyamatos, szilárdtest-alapú finomítási folyamatként értelmezi újra, amelyet a köpenykonvekció hajt. Legfontosabb, hogy rávilágít arra a kritikus visszacsatolási hurokra, amelyben az ioncsere módosítja az ásványi összetételt, ami viszont megváltoztatja a kőzet sűrűségét és reológiáját, ezáltal befolyásolva a köpenykonvekció természetét.

Lényegében ez a munka új elméleti alapot szolgáltat annak megértéséhez, hogyan szabályozza az ásványok atomi szintű viselkedése a Föld-típusú bolygók hosszú távú kémiai és termikus evolúcióját.

5.2. Javaslatok a Jövőbeli Kutatásokra

Az ebben a jelentésben létrehozott elméleti keretrendszer számos ígéretes utat nyit a jövőbeli kutatások számára, amelyek kiterjednek a kísérleti, elméleti és számítógépes területekre.

1. Kísérleti Ellenőrzés és Parametrizálás: A modell prediktív ereje függ a bemeneti paramétereinek pontosságától, különösen a szilárdtest-diffúziós együtthatók (Di) és a nagynyomású ásványfázisok termodinamikai adatainak (ΔH∘, ΔV∘) pontosságától. Kulcsfontosságú egy fókuszált kísérleti program, amely fejlett technikákat (pl. gyémánt-üllős cellák szinkrotron-sugárzással párosítva) használ e tulajdonságok mérésére mélyköpeny P-T körülmények között.79 Az ilyen kísérletek biztosítanák a modell teljes parametrizálásához szükséges adatokat, és közvetlen tesztként szolgálnának kinetikai előrejelzéseire.

2. Kiterjesztés Más Geokémiai Rendszerekre: Az itt kifejlesztett analitikai módszertan általános, és kiterjeszthető más fontos geokémiai rendszerekre is. A jövőbeli munkának a következőkre kell összpontosítania:

o Nyomelemek: A modell adaptálása a kulcsfontosságú nyomelemek megoszlásának és transzportjának leírására, amelyek eloszlása a köpenyfolyamatok nyomon követésének elsődleges eszköze.

o Illékony Anyagok és Redox-érzékeny Elemek: Olyan fajták, mint a protonok (H+) cseréjének beépítése a vízhordozás modellezésére a nominálisan vízmentes ásványokban (NAM-ok), és több vegyértékű elemek, mint a vas (Fe2+/Fe3+) bevonása a modell és a köpeny oxidációs állapotának összekapcsolására.82

o Más Égitestek: A modell alkalmazása más Föld-típusú bolygókra, például a Marsra vagy szuper-Földekre, amelyek eltérő teljes összetétellel, mérettel és belső P-T gradienssel rendelkeznek. Ez lehetővé tenné a kémiai differenciálódási folyamatok összehasonlító planetológiáját.61

3. Integráció Nagyléptékű Numerikus Modellekbe: Jelentős jövőbeli hozzájárulás lenne az analitikai megoldás implementálása számításilag hatékony rácson belüli (sub-grid) fizikai modulként a legkorszerűbb numerikus köpenykonvekciós kódokba (pl. ASPECT, CitcomS). A megoldás analitikai természete lehetővé tenné az ásványi összetétel és az ebből következő kőzettulajdonságok (sűrűség, viszkozitás) dinamikus és fizikailag szigorú frissítését minden időlépésben, anélkül, hogy a diffúziós egyenletek numerikus megoldása minden rácscellában megfizethetetlen számítási költséggel járna. Ez az integráció jelentős lépést jelentene egy új generációs, önkonzisztens, több skálájú, több fizikai jelenséget leíró bolygófejlődési modell felé, teljes mértékben áthidalva az ásványfizika és a geodinamika közötti szakadékot.

Idézett Munkák

1. (PDF) Partitioning and Partition Coefficients - ResearchGate, hozzáférés dátuma: 2025. július 30., https://www.researchgate.net/publication/320407481_Partitioning_and_Partition_Coefficients

2. Trace-Element Geochemistry, Lecture Notes 3 - MIT OpenCourseWare, hozzáférés dátuma: 2025. július 30., https://ocw.mit.edu/courses/12-479-trace-element-geochemistry-spring-2013/ba3c28248633c959277c1b66ddef74ff_MIT12_479S13_lec3.pdf

3. Mg2+-Fe2+ ORDER-DISORDER AND THE THERMODYNAMICS OF THE ORTHO- PYROXENE CRYSTALLINE SOLUTION, RRuff, hozzáférés dátuma: 2025. július 30., https://rruff.info/uploads/AM56_532.pdf

4. Mgz+-psr+ ORDER-DISORDER AND THE THERMODYNAMICS OF THE ORTHO- PYROXENE CRYSTALLINE SOLUTION SunBNnna K. SexBNa AND SUBRATA GHosE, hozzáférés dátuma: 2025. július 30., http://www.minsocam.org/ammin/AM56/AM56_532.pdf

5. Chemical Fundamentals of Geology and Environmental Geoscience - AIU Student Login, hozzáférés dátuma: 2025. július 30., http://students.aiu.edu/submissions/profiles/resources/onlineBook/B6a2U6_chemistry%20for%20geology.pdf

6. A thermodynamic model for feldspars in KAISI308-NaAlSi3O8..., hozzáférés dátuma: 2025. július 30., https://www.repository.cam.ac.uk/bitstreams/55519c2c-6d9a-4bd6-a641-e05f88d1e88c/download

7. Plagioclase feldspars: Activity-composition relations based upon Darken's quadratic formalism and Landau theory, hozzáférés dátuma: 2025. július 30., https://msaweb.org/AmMin/AM77/AM77_53.pdf

8. Earth's mantle - Wikipedia, en.wikipedia.org, hozzáférés dátuma: 2025. július 30., https://en.wikipedia.org/wiki/Earth%27s_mantle#:~:text=VII%20to%20form.-,Temperature%20and%20pressure,at%20the%20core%2Dmantle%20boundary.

9. Earth's mantle | EBSCO Research Starters, hozzáférés dátuma: 2025. július 30., https://www.ebsco.com/research-starters/earth-and-atmospheric-sciences/earths-mantle

10. Upper mantle - Wikipedia, hozzáférés dátuma: 2025. július 30., https://en.wikipedia.org/wiki/Upper_mantle

11. Deep Earth and Mineral Physics - Elements Magazine, hozzáférés dátuma: 2025. július 30., https://www.elementsmagazine.org/deep-earth-and-mineral-physics/

12. Iron force constants of bridgmanite at high pressure: Implications for iron isotope fractionation in the deep mantle - Jackson School of Geosciences, hozzáférés dátuma: 2025. július 30., https://www.jsg.utexas.edu/lin/files/WangBridgmaniteForceConstantGeochimica2021.pdf

13. Thermodynamics of mantle minerals - I. Physical properties | Geophysical Journal International | Oxford Academic, hozzáférés dátuma: 2025. július 30., https://academic.oup.com/gji/article/162/2/610/2051140

14. Thermodynamics of Mantle Minerals: 1. Physical Properties, hozzáférés dátuma: 2025. július 30., https://www.perplex.ethz.ch/thermo_course/chapter_5/stixrude_GJI_05.pdf

15. (PDF) Thermodynamics of mantle minerals - I. Physical properties - ResearchGate, hozzáférés dátuma: 2025. július 30., https://www.researchgate.net/publication/227721942_Thermodynamics_of_mantle_minerals_-_I_Physical_properties

16. Thermodynamics of the Earth's Mantle - UCL, hozzáférés dátuma: 2025. július 30., http://www.homepages.ucl.ac.uk/~ucfbls0/stixrudelithgowbertelloni_10.pdf

17. Thermodynamic Formulation for Partition Equilibria of Ca, Fe and Mg between Olivine and Clinopyroxene in the System CaO-FeO-MgO-SiO2 - ResearchGate, hozzáférés dátuma: 2025. július 30., https://www.researchgate.net/publication/317590468_Thermodynamic_Formulation_for_Partition_Equilibria_of_Ca_Fe_and_Mg_between_Olivine_and_Clinopyroxene_in_the_System_CaO-FeO-MgO-SiO2

18. Thermodynamic formulations of (Ca, Fe,Mg)2SiO4 olivine - ResearchGate, hozzáférés dátuma: 2025. július 30., https://www.researchgate.net/publication/250163966_Thermodynamic_formulations_of_CaFeMg2SiO4_olivine

19. High-pressure experiments complete picture of the Earth's mantle - UCL, hozzáférés dátuma: 2025. július 30., https://www.ucl.ac.uk/mathematical-physical-sciences/news/2019/aug/high-pressure-experiments-complete-picture-earths-mantle

20. World's First Realization of Ultrahigh Pressure and Ultrahigh Temperature at the Earth's Center - Finally reaching the Earth's Core - SPring-8, hozzáférés dátuma: 2025. július 30., http://www.spring8.or.jp/en/news_publications/research_highlights/no_57/

21. Ion Exchange Kinetics. A Nonlinear Diffusion Problem, hozzáférés dátuma: 2025. július 30., https://authors.library.caltech.edu/records/dxr4p-0hn34

22. Reactive Transport at the Crossroads | Reviews in Mineralogy and..., hozzáférés dátuma: 2025. július 30., https://pubs.geoscienceworld.org/msa/rimg/article/85/1/1/573287/Reactive-Transport-at-the-Crossroads

23. Solutions to a reduced Poisson-Nernst-Planck system and determination of reaction rates, hozzáférés dátuma: 2025. július 30., https://pmc.ncbi.nlm.nih.gov/articles/PMC2835320/

24. Analytical solution of the Poisson-Nernst-Planck-Stokes equations in a cylindrical channel | Proceedings of the Royal Society A: Mathematical, Physical and Engineering Sciences - Journals, hozzáférés dátuma: 2025. július 30., https://royalsocietypublishing.org/doi/10.1098/rspa.2011.0080

25. Poisson-Nernst-Planck Equations for Simulating Biomolecular Diffusion-Reaction Processes I: Finite Element Solutions, hozzáférés dátuma: 2025. július 30., https://pmc.ncbi.nlm.nih.gov/articles/PMC2922884/

26. Improvement of numerical approximation of coupled multiphase multicomponent flow with reactive geochemical transport in porous media | Oil & Gas Science and Technology, hozzáférés dátuma: 2025. július 30., https://ogst.ifpenergiesnouvelles.fr/articles/ogst/full_html/2018/01/ogst180053/ogst180053.html

27. Reactive transport modeling: An essential tool and a new research approach for the Earth sciences, hozzáférés dátuma: 2025. július 30., https://www2.lbl.gov/Science-Articles/Archive/sabl/2006/Mar/Reactive-Transport-Modeling.pdf

28. Mantle convection - The People of Earth & Planetary Sciences, hozzáférés dátuma: 2025. július 30., https://people.earth.yale.edu/sites/default/files/files/Bercovici/17_MantlConvection-ESEG2011-2_0.pdf

29. 3-D Spherical Models of Mantle Convection with Floating Continents, hozzáférés dátuma: 2025. július 30., https://pubs.usgs.gov/of/2000/0218/report.pdf

30. 7.02 Physics of Mantle Convection, hozzáférés dátuma: 2025. július 30., https://geo.mff.cuni.cz/~hk/geodynamika/ToG-07.02-Ricard.pdf

31. A review of numerical methods for nonlinear partial differential equations - ResearchGate, hozzáférés dátuma: 2025. július 30., https://www.researchgate.net/publication/268998531_A_review_of_numerical_methods_for_nonlinear_partial_differential_equations

32. A review of numerical methods for nonlinear partial differential equations - Semantic Scholar, hozzáférés dátuma: 2025. július 30., https://www.semanticscholar.org/paper/A-review-of-numerical-methods-for-nonlinear-partial-Tadmor/04d765c7b4d94064945edc57d6d4ffd7f1d588d0

33. for the First Course. Part III: Adomian Decomposition Method, hozzáférés dátuma: 2025. július 30., https://www.cfm.brown.edu/people/dobrush/am33/Mathematica/ch3/adm.html

34. The Adomian decomposition method for solving a system of non-linear equations. - University of Zawia DSpace, hozzáférés dátuma: 2025. július 30., https://dspace.zu.edu.ly/bitstream/handle/1/1240/adomain%203.pdf?sequence=1&isAllowed=y

35. A two-step Adomian decomposition method - ResearchGate, hozzáférés dátuma: 2025. július 30., https://www.researchgate.net/publication/220557629_A_two-step_Adomian_decomposition_method

36. Adomian decomposition method reformulated using dimensionless nonlinear perturbation theory - arXiv, hozzáférés dátuma: 2025. július 30., https://arxiv.org/html/2501.10398v1

37. LIMITATIONS OF ADOMIAN DECOMPOSITION AND HOMOTOPY METHODS, hozzáférés dátuma: 2025. július 30., https://www.researchgate.net/publication/266524483_LIMITATIONS_OF_ADOMIAN_DECOMPOSITION_AND_HOMOTOPY_METHODS

38. Homotopy analysis method - Wikipedia, hozzáférés dátuma: 2025. július 30., https://en.wikipedia.org/wiki/Homotopy_analysis_method

39. Homotopy Analysis Method for the Approximate Solution of the SIRC Epidemic Model, hozzáférés dátuma: 2025. július 30., https://internationalpubls.com/index.php/cana/article/view/3307

40. Equivalence of Ratio and Residual Approaches in the Homotopy Analysis Method and Some Applications in Nonlinear Science and Engineering, hozzáférés dátuma: 2025. július 30., https://www.techscience.com/CMES/v120n1/27490

41. Study of coupled nonlinear partial differential equations for finding exact analytical solutions - ResearchGate, hozzáférés dátuma: 2025. július 30., https://www.researchgate.net/journal/Royal-Society-Open-Science-2054-5703/publication/279751821_Study_of_Coupled_Nonlinear_Partial_Differential_Equations_for_finding_Exact_Analytical_Solutions/links/62a3216ea3fe3e3df86c6001/Study-of-Coupled-Nonlinear-Partial-Differential-Equations-for-finding-Exact-Analytical-Solutions.pdf?origin=journalDetail

42. Geochemical modeling - Wikipedia, hozzáférés dátuma: 2025. július 30., https://en.wikipedia.org/wiki/Geochemical_modeling

43. 101 geodynamic modelling: how to design, interpret, and communicate numerical studies of the solid Earth, hozzáférés dátuma: 2025. július 30., https://par.nsf.gov/servlets/purl/10347047

44. (PDF) Analytical solutions of nonlinear and variable-parameter transport equations for verification of numerical solvers - ResearchGate, hozzáférés dátuma: 2025. július 30., https://www.researchgate.net/publication/272030349_Analytical_solutions_of_nonlinear_and_variable-parameter_transport_equations_for_verification_of_numerical_solvers

45. Solutions of Poisson-Nernst Planck Equations with Ion Interaction - Scirp.org., hozzáférés dátuma: 2025. július 30., https://www.scirp.org/journal/paperinformation?paperid=116249

46. (PDF) Geophysics of Chemical Heterogeneity in the Mantle - ResearchGate, hozzáférés dátuma: 2025. július 30., https://www.researchgate.net/publication/234149074_Geophysics_of_Chemical_Heterogeneity_in_the_Mantle

47. MANTLE MIXING: The Generation, Preservation, and Destruction of Chemical Heterogeneity - faculty.washington.edu, hozzáférés dátuma: 2025. július 30., http://faculty.washington.edu/stn/ess_501/reading/van_Keken_Mantle_heteroneity_Ann_Revs_2002.pdf

48. Siderophile Elements in Tracing Planetary Formation and Evolution - PMC, hozzáférés dátuma: 2025. július 30., https://pmc.ncbi.nlm.nih.gov/articles/PMC6376203/

49. The prevalence of kilometer-scale heterogeneity in the source region of MORB upper mantle, hozzáférés dátuma: 2025. július 30., https://pmc.ncbi.nlm.nih.gov/articles/PMC5706740/

50. "Iron Isotope Cosmochemistry" by Kun Wang - Washington University Open Scholarship, hozzáférés dátuma: 2025. július 30., https://openscholarship.wustl.edu/etd/1189/

51. Compositional heterogeneity in the mantle transition zone - UCL Discovery, hozzáférés dátuma: 2025. július 30., https://discovery.ucl.ac.uk/10155292/1/Ballmer_Goes-et-al_MTZreview_NATREVEARTHEHVIRON-edited-accepted.pdf

52. (PDF) Mineral Physics in ThermoChemical Mantle Models - ResearchGate, hozzáférés dátuma: 2025. július 30., https://www.researchgate.net/publication/252686270_Mineral_Physics_in_ThermoChemical_Mantle_Models

53. 101 geodynamic modelling: how to design, interpret, and communicate numerical studies of the solid Earth - SE, hozzáférés dátuma: 2025. július 30., https://se.copernicus.org/articles/13/583/2022/se-13-583-2022.html

54. Bernhard Schuberth - the Geophysics Homepage, hozzáférés dátuma: 2025. július 30., https://www.geophysik.uni-muenchen.de/~bernhard/index.php?l=about&lng=en

55. CSEDI Collaborative Research: Understanding what we see in the lower mantle - mineral physics interpretation of seismic tomographic images - Research Projects, hozzáférés dátuma: 2025. július 30., https://people.climate.columbia.edu/projects/view/2180

56. Geodynamic data-driven modelling: bridging the gap between observations and numerical models - EGU Blogs, hozzáférés dátuma: 2025. július 30., https://blogs.egu.eu/divisions/gd/2023/04/26/geodynamic-data-driven-modelling-bridging-the-gap-between-observations-and-numerical-models/

Idézett munkák

1. An Analytical Model for the Dynamics of Cation Exchange in Planetary Mantle Silicates under Variable P-T Conditions.PDF

2. Electrodiffusion Phenomena in Neuroscience and the Nernst–Planck–Poisson Equations - MDPI, hozzáférés dátuma: július 30, 2025, https://www.mdpi.com/2673-3293/2/2/14

3. Poisson–Boltzmann–Nernst–Planck model - PMC - PubMed Central, hozzáférés dátuma: július 30, 2025, https://pmc.ncbi.nlm.nih.gov/articles/PMC3122111/

4. Homotopy analysis method - Wikipedia, hozzáférés dátuma: július 30, 2025, https://en.wikipedia.org/wiki/Homotopy_analysis_method

5. Geochemical investigation of geological formations from the perspective of hydrogen storage potential, hozzáférés dátuma: július 30, 2025, https://ktdi.elte.hu/dstore/document/12593/GelencserO_PhD_Thesis.pdf

6. Reológia – Wikipédia, hozzáférés dátuma: július 30, 2025, https://hu.wikipedia.org/wiki/Reol%C3%B3gia

7. subgrid-scale physics, hozzáférés dátuma: július 30, 2025, http://astro.vaporia.com/start/subgridphysics.html

8. Subgrid Scale Physics - NEMO ocean model, hozzáférés dátuma: július 30, 2025, https://www.nemo-ocean.eu/doc/node10.html

9. Chemical Research Center of the Hungarian Academy of Sciences, hozzáférés dátuma: július 30, 2025, https://scispace.com/institutions/chemical-research-center-of-the-hungarian-academy-of-oj269poi?paper_page=35

10. Gibbs-féle szabadenergia-görbe számítások : r/metallurgy - Reddit, hozzáférés dátuma: július 30, 2025, https://www.reddit.com/r/metallurgy/comments/m2wbre/gibbs_free_energy_curve_calculations/?tl=hu

11. GEOKÉMIA - Környezettudomány Tanszék, Sapientia EMTE, hozzáférés dátuma: július 30, 2025, http://kt.sapientia.ro/data/dokumentumok/segedanyag-geokemia.pdf

12. A területi geokémiai kutatás elméleti és gyakorlati módszerei, hozzáférés dátuma: július 30, 2025, http://mek.oszk.hu/19000/19095/19095.pdf

13. The New Crystal Caves in Hungary at Beremend and Nagyharsány, hozzáférés dátuma: július 30, 2025, https://www.cavernas.org.br/wp-content/uploads/2021/02/26CBE_509-511.pdf

14. Szervetlen kémiai nevezéktan - IR-11.6.3. A krisztallográfiai, hozzáférés dátuma: július 30, 2025, https://mersz.hu/dokumentum/m102szkn__340/

15. A kémiai reakciók egyensúlyi leírása, hozzáférés dátuma: július 30, 2025, http://pc.ch.bme.hu/oktatas/konyvek/anal/Muszaki-Menedzser-Kemia/Kemia_I_eloadasi_foliak_MJ-2010-2-es-3-ZH/MK%E9mia_I_5_reakcioegyensulyok.ppt

16. tömeghatás törvénye - Lexikon :: - Kislexikon, hozzáférés dátuma: július 30, 2025, http://www.kislexikon.hu/tomeghatas_torvenye.html

17. Moláris tömeg – Wikipédia, hozzáférés dátuma: július 30, 2025, https://hu.wikipedia.org/wiki/Mol%C3%A1ris_t%C3%B6meg

18. Activity meaning in Hungarian - DictZone, hozzáférés dátuma: július 30, 2025, https://dictzone.com/english-hungarian-dictionary/activity

19. Aktivitási együttható – Wikipédia, hozzáférés dátuma: július 30, 2025, https://hu.wikipedia.org/wiki/Aktivit%C3%A1si_egy%C3%BCtthat%C3%B3

20. Van 't Hoff equation, hozzáférés dátuma: július 30, 2025, https://en.wikipedia.org/wiki/Van_'t_Hoff_equation

21. ferromagnesian clay minerals: Topics by Science.gov, hozzáférés dátuma: július 30, 2025, https://www.science.gov/topicpages/f/ferromagnesian+clay+minerals

22. Ásvány- és kőzettan - Környezeti ásványtan, hozzáférés dátuma: július 30, 2025, http://asvanytan.nyf.hu/book/export/html/181

23. SZERVES KÉMIAI ALAPISMERETEK - ELTE, hozzáférés dátuma: július 30, 2025, http://olimpia.chem.elte.hu/evek/2009/aktual/szerves.pdf

24. Best Sodium Feldspar exporter and manufacturer globally. Supplier of highquality Sodium Feldspar. Sodium Feldspar Manufacturer Exporter in Hungary - Hind Exports, hozzáférés dátuma: július 30, 2025, https://www.hindexports.com/sodium-feldspar-manufacturer-exports-in-hungary.php

25. Földpátok – Wikipédia, hozzáférés dátuma: július 30, 2025, https://hu.wikipedia.org/wiki/F%C3%B6ldp%C3%A1tok

26. Pioneering attempts in Hungary on treating quantum systems | Request PDF, hozzáférés dátuma: július 30, 2025, https://www.researchgate.net/publication/389493244_Pioneering_attempts_in_Hungary_on_treating_quantum_systems

27. The Hungarian Revolution of 1956 – Association for Diplomatic Studies & Training, hozzáférés dátuma: július 30, 2025, https://adst.org/2013/10/the-hungarian-revolution-of-1956/

28. (PDF) Ioncsere-folyamatok az agyagásványok ... - ResearchGate, hozzáférés dátuma: július 30, 2025, https://www.researchgate.net/publication/317556645_Ioncsere-folyamatok_az_agyagasvanyok_nanolaboratoriumaban

29. Tudnivalók a kapcsolt vállalkozásokról - Kontír Trend Kft., hozzáférés dátuma: július 30, 2025, https://kontirtrend.hu/2023/11/21/tudnivalok-a-kapcsolt-vallalkozasokrol/

30. Plagioclase feldspars: Activity-composition relations based upon Darken's quadratic formalism and Landau theory, hozzáférés dátuma: július 30, 2025, https://msaweb.org/AmMin/AM77/AM77_53.pdf

31. Landau theory - Wikipedia, hozzáférés dátuma: július 30, 2025, https://en.wikipedia.org/wiki/Landau_theory

32. Lev Davidovics Landau - Wikipédia, hozzáférés dátuma: július 30, 2025, https://hu.wikipedia.org/wiki/Lev_Davidovics_Landau

33. Fázisátalakulások, hozzáférés dátuma: július 30, 2025, https://atomcsill.elte.hu/foliak/kiserletek_12_evf/ACS_20161208_Fazisatalakulasok.pdf

34. Wadsleyite Mineral Data - Mineralogy Database, hozzáférés dátuma: július 30, 2025, http://webmineral.com/data/Wadsleyite.shtml

35. Ringwoodite Mineral Data - Mineralogy Database, hozzáférés dátuma: július 30, 2025, http://webmineral.com/data/Ringwoodite.shtml

36. Silicate perovskite - Wikipedia, hozzáférés dátuma: július 30, 2025, https://en.wikipedia.org/wiki/Silicate_perovskite

37. Ferropericlase - Wikipedia, hozzáférés dátuma: július 30, 2025, https://en.wikipedia.org/wiki/Ferropericlase

38. Hamarabb hűlhet ki és válhat bárki számára lakhatatlanná a Föld, mint gondoltuk - Qubit, hozzáférés dátuma: július 30, 2025, https://qubit.hu/2022/01/17/hamarabb-hulhet-ki-es-valhat-barki-szamara-lakhatatlanna-a-fold-mint-gondoltuk

39. Steady-State Electrodiffusion from the Nernst−Planck Equation Coupled to Local Equilibrium Monte Carlo Simulations - Rush University Medical Center, hozzáférés dátuma: július 30, 2025, https://ftp.rush.edu/users/molebio/Bob_Eisenberg/CV%20on%20moleobio-Bob_Eisenberg/1%20Selectivity%20Papers%20Selectivity%20Slides/Boda%20Dirk%20JCTC_2012-NP-LEMC-w-suppl.pdf

40. Fick's laws of diffusion - Wikipedia, hozzáférés dátuma: július 30, 2025, https://en.wikipedia.org/wiki/Fick%27s_laws_of_diffusion

41. A Fick-féle diffúziós törvény a fizikában és a biológiában alkalmazott ..., hozzáférés dátuma: július 30, 2025, https://klordioxid.com/wp-content/uploads/2024/09/fick.pdf

42. Diffusion along the perivascular space influenced by handedness and language lateralisation | Brain Communications | Oxford Academic, hozzáférés dátuma: július 30, 2025, https://academic.oup.com/braincomms/article/7/4/fcaf252/8173005

43. Félvezető diódák, LED-ek - Bevezetés az elektronikába, hozzáférés dátuma: július 30, 2025, https://megtestesules.info/hobbielektronika/elektronika/04_diodes_led.html

44. Organic-Inorganic Perovskite Precursors - TCI Chemicals, hozzáférés dátuma: július 30, 2025, https://www.tcichemicals.com/CA/en/c/12969

45. What is a Dielectric Constant (Dk) of plastic materials? - doEEEt, hozzáférés dátuma: július 30, 2025, https://www.doeeet.com/content/eee-components/passive-components/what-is-a-dielectric-constant-of-plastic-materials/

46. en.wikipedia.org, hozzáférés dátuma: július 30, 2025, https://en.wikipedia.org/wiki/Relative_permittivity

47. Continuity equation - KSB, hozzáférés dátuma: július 30, 2025, https://www.ksb.com/en-global/centrifugal-pump-lexicon/article/continuity-equation-1116870

48. en.wikipedia.org, hozzáférés dátuma: július 30, 2025, https://en.wikipedia.org/wiki/Continuity_equation

49. Parciális differenciálegyenletek, hozzáférés dátuma: július 30, 2025, http://www.sze.hu/~lotfi/Par-diff.pdf

50. Parciális differenciálegyenlet - Wikipédia, hozzáférés dátuma: július 30, 2025, https://hu.wikipedia.org/wiki/Parci%C3%A1lis_differenci%C3%A1legyenlet

51. Debye number - QED - Princeton Plasma Physics Laboratory, hozzáférés dátuma: július 30, 2025, https://w3.pppl.gov/~dboyle/PlasmaWiki/qed.princeton.edu/main/PlasmaWiki/Debye_number.html

52. Result of grain analysis for the 50 [µm] strategies - ResearchGate, hozzáférés dátuma: július 30, 2025, https://www.researchgate.net/figure/Result-of-grain-analysis-for-the-50-m-strategies_tbl2_284593785

53. Debye model - Wikipedia, hozzáférés dátuma: július 30, 2025, https://en.wikipedia.org/wiki/Debye_model

54. Damköhler numbers - Wikipedia, hozzáférés dátuma: július 30, 2025, https://en.wikipedia.org/wiki/Damk%C3%B6hler_numbers

55. Dirichlet boundary condition - Wikipedia, hozzáférés dátuma: július 30, 2025, https://en.wikipedia.org/wiki/Dirichlet_boundary_condition

56. epidot | Rezső a geológus sün, hozzáférés dátuma: július 30, 2025, https://rezsoageologussun.wordpress.com/tag/epidot/

57. Search - DKUM, hozzáférés dátuma: július 30, 2025, https://dk.um.si/Iskanje.php?type=napredno&lang=eng&stl0=KljucneBesede&niz0=homotopy+analysis+method

58. Adomian decomposition method - Wikipedia, hozzáférés dátuma: július 30, 2025, https://en.wikipedia.org/wiki/Adomian_decomposition_method

59. Adiabatic process - Wikipedia, hozzáférés dátuma: július 30, 2025, https://en.wikipedia.org/wiki/Adiabatic_process

60. Beyond Perturbation | Introduction to the Homotopy Analysis Method | S, hozzáférés dátuma: július 30, 2025, https://www.taylorfrancis.com/books/mono/10.1201/9780203491164/beyond-perturbation-shijun-liao

61. Comparative planetary science - Wikipedia, hozzáférés dátuma: július 30, 2025, https://en.wikipedia.org/wiki/Comparative_planetary_science

62. Homotopy Analysis Method for Solving Initial Value Problems of Second Order with Discontinuities - Scirp.org., hozzáférés dátuma: július 30, 2025, https://www.scirp.org/journal/paperinformation?paperid=93048

63. Lineáris leképezés - Wikipédia, hozzáférés dátuma: július 30, 2025, https://hu.wikipedia.org/wiki/Line%C3%A1ris_lek%C3%A9pez%C3%A9s

64. Green's function - Wikipedia, hozzáférés dátuma: július 30, 2025, https://en.wikipedia.org/wiki/Green%27s_function

65. A változók szétválasztása a hővezetés egyenletére, hozzáférés dátuma: július 30, 2025, https://www.math.ubbcluj.ro/~andrasz/dokumentumok/pde/hovezetes.pdf

66. Mit mutat a kationcsere-kapacitás? - Smart Farming Kft., hozzáférés dátuma: július 30, 2025, https://smartfarming.hu/mit-mutat-a-kationcsere-kapacitas/

67. A cement geokémiai modellezése - ELTE, hozzáférés dátuma: július 30, 2025, https://www.elte.hu/content/a-cement-geokemiai-modellezese.e.12272

68. Az Európai Unió L 54/2016 - EUR-Lex, hozzáférés dátuma: július 30, 2025, https://eur-lex.europa.eu/legal-content/HU/TXT/HTML/?uri=OJ%3AL%3A2016%3A054%3AFULL

69. chemical symbol | English-Hungarian translation - dict.cc, hozzáférés dátuma: július 30, 2025, https://m.dict.cc/english-hungarian/chemical+symbol.html

70. Információk a Pannon-medence alatti kőzetburokról - Arcanum, hozzáférés dátuma: július 30, 2025, https://www.arcanum.com/hu/online-kiadvanyok/pannon-pannon-enciklopedia-1/magyarorszag-foldje-1D58/a-karpat-pannon-terseg-lemeztektonikai-ertelmezese-1ED3/a-medencebelseji-alkalibazalt-vulkanossag-harangi-szabolcs-1FDB/informaciok-a-pannon-medence-alatti-kozetburokrol-1FE6/

71. Az adiabatikus hűtés hozzájárul a beltéri levegő minőségéhez és a dolgozók kényelméhez - Obera, hozzáférés dátuma: július 30, 2025, https://obera.fr/hu/a-mi-tanacsunk/quels-sont-les-avantages-du-rafraichissement-adiabatique/

72. Thermal conduction, hozzáférés dátuma: július 30, 2025, https://en.wikipedia.org/wiki/Thermal_conduction

73. The Henryk Niewodniczański - IFJ PAN, hozzáférés dátuma: július 30, 2025, https://www.ifj.edu.pl/badania/publikacje/raporty/2009/2028.pdf

74. A Világegyetem története 200 vagy kevesebb szóban - LFG.HU, hozzáférés dátuma: július 30, 2025, https://lfg.hu/6771/vegyes/humorka/a-vilagegyetem-tortenete-200-vagy-kevesebb-szoban/

75. The Equilibrium Institute's economic forecast for Hungary - Egyensúly Intézet, hozzáférés dátuma: július 30, 2025, https://egyensulyintezet.hu/en/the-equilibrium-institutes-economic-forecast-for-hungary-q4/

76. Eclogites and their geodynamic interpretation: A history - ResearchGate, hozzáférés dátuma: július 30, 2025, https://www.researchgate.net/publication/223579660_Eclogites_and_their_geodynamic_interpretation_A_history

77. Eklogit | A Pallas nagy lexikona | Kézikönyvtár - Arcanum, hozzáférés dátuma: július 30, 2025, https://www.arcanum.com/hu/online-kiadvanyok/Lexikonok-a-pallas-nagy-lexikona-2/e-e-7C62/eklogit-8125/

78. Equation of state, hozzáférés dátuma: július 30, 2025, https://en.wikipedia.org/wiki/Equation_of_state

79. Diamond anvil cell - Wikipedia, hozzáférés dátuma: július 30, 2025, https://en.wikipedia.org/wiki/Diamond_anvil_cell

80. hu.wikipedia.org, hozzáférés dátuma: július 30, 2025, https://hu.wikipedia.org/wiki/Szinkrotronsug%C3%A1rz%C3%A1s#:~:text=A%20szinkrotronsug%C3%A1rz%C3%A1s%20elektrom%C3%A1gneses%20sug%C3%A1rz%C3%A1s%2C%20amelyet,a%20m%C3%A1gneses%20t%C3%A9r%20hat%C3%A1s%C3%A1ra%20v%C3%A1ltozik.

81. Beszámoló nemzetközi rövidkurzuson való részvételről - ELTE, hozzáférés dátuma: július 30, 2025, https://www.elte.hu/dstore/document/4872/Spr%C3%A1nitz%20Tam%C3%A1s_besz%C3%A1mol%C3%B3.pdf

82. Trace element detection in anhydrous minerals by micro-scale quantitative nuclear magnetic resonance spectroscopy - PubMed Central, hozzáférés dátuma: július 30, 2025, https://pmc.ncbi.nlm.nih.gov/articles/PMC11344839/

83. Nominally Anhydrous Minerals and Earth's Deep Water Cycle, hozzáférés dátuma: július 30, 2025, http://vigarano.ep.sci.hokudai.ac.jp/shogo/adv_geochemistry2015/Week%209/Background%20reading/Smyth%20and%20Jacobsen%202006.pdf


Nincsenek megjegyzések:

Megjegyzés küldése