A Kialakuló Kristályosodási Frontok Analitikus Dinamikája Bolygói Magmaóceánokban

A Kialakuló Kristályosodási Frontok Analitikus Dinamikája Bolygói Magmaóceánokban

Lengyel Ferenc

2025. július

Vezetői Összefoglaló

A bolygói magmaóceánok megszilárdulása a kőzetbolygók kialakulásának és fejlődésének alapvető korszakát jelenti. Ez a folyamat az elsődleges, bolygóméretű differenciációs esemény, amely létrehozza a kéreg, a köpeny és a légkör kezdeti állapotát, ezáltal meghatározva minden későbbi geológiai és potenciálisan biológiai evolúció peremfeltételeit. A folyamat dinamikája azonban továbbra is intenzív elméleti és kísérleti kutatások tárgyát képezi. Jelen jelentés egy központi és megoldatlan elméleti kihívással foglalkozik: a hűlő magmaóceánon áthaladó kristályosodási front dinamikájának analitikus megoldásainak levezetésével. E jelenség kvantitatív, első elveken alapuló modellje elengedhetetlen a bolygói kéregképződés, a köpeny rétegződésének és az illóanyagok kigázosodási történetének előrejelző modelljeinek fejlesztéséhez.

A vizsgálat magját a fizikai rendszer többkomponensű, egydimenziós, mozgó határfelülettel rendelkező reakció-diffúzió problémaként való megfogalmazása képezi. Ez a keretrendszer, amely a klasszikus Stefan-probléma általánosítása, megragadja a hőtranszport, a kémiai diffúzió és a fázisváltozás összekapcsolt fizikájának lényegét. Az irányító rendszer csatolt parciális differenciálegyenletekből áll, amelyek a szilárd és folyékony fázisok hőmérsékleti mezejét, valamint a folyékony olvadékban lévő több kémiai komponens koncentrációmezejét írják le. A kritikus komplexitás a mozgó határfelületen jelentkező nemlineáris csatolásból fakad: a határfelületi hőmérséklet a likvidusz-hőmérséklet, amely a helyi olvadékösszetétel (beleértve az illóanyagokat is) függvénye, a határfelület sebességét pedig csatolt energia- és tömegmérleg-feltételek (a Stefan-feltétel és a kémiai megoszlási feltételek) szabályozzák, amelyek a fronton lévő hőmérsékleti és kémiai gradiensektől függenek.

Felismerve, hogy a teljes nemlineáris rendszer matematikailag kezelhetetlen egy általános, zárt alakú megoldás számára, ez a jelentés egy sor jól bevált analitikus technikát alkalmaz egy közelítő, de fizikailag mély betekintést nyújtó megoldás levezetésére. Az elsődleges módszer a hasonlósági transzformáció, amely az irányító parciális differenciálegyenleteket csatolt, nemlineáris közönséges differenciálegyenletek rendszerévé redukálja. Ezt a transzformált rendszert ezután analitikusan, reguláris perturbációs sorfejtéssel oldjuk meg a Stefan-szám, egy dimenzió nélküli paraméter szerint, amely a látens és a szenzibilis hő arányát képviseli, és amely a szilikátos rendszerek esetében jellemzően kicsi.

A levezetett analitikus megoldás megadja a kristályosodási front helyzetét az idő függvényében, s(t), valamint a hőmérséklet- és kémiai koncentrációprofilokat a szilárd és folyékony köpenyben. A legfontosabb megállapítások a következők:

1. A front sebessége és annak függőségei: A modell explicit kifejezést ad a kristályosodási front sebességére. Bemutatja, hogy a megszilárdult köpeny vastagsága arányos az idő négyzetgyökével, ami egy jellegzetes diffúzív skálatörvény. Kimutattuk, hogy a növekedés üteme a bolygó hűlési sebességének (amelyet a termikus peremfeltételek határoznak meg), a magmaóceán kezdeti tömbösszetételének, az ásvány-olvadék megoszlási hányadosok által leírt kémiai frakcionáció hatékonyságának és az oldott illóanyagok koncentrációjának komplex függvénye.

2. Az illóanyag-telítettség mint nemlineáris kiváltó ok: A modell kvantitatívan leírja az inkompatibilis elemek és illóanyagok feldúsulásának folyamatát egy kémiai határrétegben a haladó front előtt. Előrejelzi, hogy ez a feldúsulás az olvadékot az illóanyag-telítettség állapotába hozhatja. A telítődéskor az illóanyagok kiválása éles változást okoz a likvidusz-hőmérsékletben, ami viszont a kristályosodási front nemlineáris gyorsulásához vezet. Ez egy robusztus fizikai mechanizmust biztosít a gyors, késői stádiumú kigázosodási eseményekhez, amelyek kritikusak a bolygó ősi légkörének kialakulásában.

3. Kémiai differenciáció és fejlett rezervoárok kialakulása: A kémiai koncentrációprofilokra adott megoldás közvetlen analitikus modellt szolgáltat a magmás differenciációra a frakcionált kristályosodás révén. Előrejelzi a maradék olvadék kémiai evolúcióját a megszilárdulás előrehaladtával, kvantitatív alapot kínálva a magasan fejlett, inkompatibilis elemekben gazdag litológiák, mint például a Holdon található KREEP (kálium, ritkaföldfémek, foszfor) rezervoár kialakulásának megértéséhez.

4. A kristályosodás-redox-atmoszféra visszacsatolási hurok: A jelentés rávilágít egy kritikus visszacsatolási mechanizmusra, amely összekapcsolja a köpeny megszilárdulását a légkör fejlődésével. Az olyan ásványok kristályosodása, mint az olivin és a piroxén, előnyben részesítve vonja ki a vas(II)-iont (Fe2+) az olvadékból, fokozatosan oxidálva a maradék folyadékot. A magmaóceán felszínének ez a változó redoxállapota diktálja a kigázosodott illóanyagok kémiai formáját (pl. a CO/H2 és a CO2/H2O arányát). A keletkező légkör összetétele pedig visszahatva szabályozza annak sugárzási tulajdonságait, és így magának a bolygónak a hűlési sebességét. A modell keretet biztosít ennek a visszacsatolásnak a számszerűsítésére, bemutatva, hogy a hűlési sebesség nem független paraméter, hanem endogén módon kapcsolódik a kristályosodási folyamathoz.

A jelentésben kidolgozott analitikus keretrendszer jelentős előrelépést kínál a korai bolygófejlődés kvantitatív modellezésében. Új, fizikailag megalapozott eszközt nyújt az elsődleges kéreg kialakulásának időskáláinak és a köpenyből történő illóanyag-felszabadulás történetének előrejelzésére. Azáltal, hogy a mélybelső dinamikáját összekapcsolja a légkör összetételével, ezek a megállapítások mélyreható következményekkel bírnak a bolygók lakhatóságának kezdeti feltételeinek és az olyan bolygók, mint a Föld, a Vénusz és a Mars eltérő fejlődési útjainak megértésében.

1. Szekció: A Magmaóceánok Megszilárdulásának Elméleti Keretrendszere

A kőzetbolygók keletkezése és későbbi evolúciója rendkívüli energiájú és összetettségű folyamat. Ezen égitestek legkorábbi történetének központi eleme a magmaóceán fogalma – egy globális vagy közel globális olvadt szilikátkőzet-réteg. E magmaóceán lehűlése és megszilárdulása vitathatatlanul a legjelentősebb differenciációs esemény egy bolygó életében, amely megalapozza az összes későbbi geológiai tevékenységet. E folyamat dinamikájának megértése ezért alapvető a bolygótudomány számára. Ez a szekció bemutatja a magmaóceán-kristályosodás tudományos kontextusát, matematikailag kezelhető mozgó határfelületi problémaként keretezve azt, és formálisan meghatározza e vizsgálat célkitűzéseit.

1.1 A Magmaóceánok Szerepe a Bolygófejlődésben

A belső Naprendszer kőzetbolygóinak kialakulása erőszakos és energetikus folyamat volt, amelyet planetezimálok akkréciója és nagyméretű ütközések jellemeztek. Két elsődleges energiaforrás biztosította, hogy ezek a születőben lévő égitestek egy vagy több olyan fázison menjenek keresztül, amelyben köpenyük részben vagy teljesen megolvadt: a korán kialakuló égitestekben a rövid felezési idejű radionuklidok, mint a ^{26}Al és a ^{60}Fe bomlása, valamint a később kialakuló, nagyobb bolygók esetében az akkréciós becsapódások során felszabaduló hatalmas kinetikus energia. Az óriásbecsapódásokat, mint amilyen a Föld Holdját feltételezhetően létrehozó esemény volt, a bolygóképződés gyakori következményének tekintik, és elegendőek lettek volna egy mély magmaóceán létrehozásához. Következésképpen a magmaóceánokat a kőzetbolygók életciklusának általános szakaszának tartják.

E hatalmas magmatestek megszilárdulása az a kohó, amelyben egy bolygó elsődleges szerkezete kialakul. Ez a bolygói differenciáció első és legmélyrehatóbb epizódja, amely meghatározza az elemek és illóanyagok kezdeti eloszlását a belső és a felszín között. Ahogy a magmaóceán hűl és kristályosodik, egy frakcionált kristályosodásnak nevezett folyamat játszódik le, amelynek során különböző ásványok válnak ki az olvadékból különböző hőmérsékleten és nyomáson. E kristályok leülepedése vagy felúszása megváltoztatja a megmaradó folyékony magma összetételét, ami egy kémiailag rétegzett köpeny és egy összetételében eltérő elsődleges kéreg kialakulásához vezet. Például a holdi felföldek, amelyek túlnyomórészt anortozit kőzetből állnak, széles körben elfogadott nézet szerint az alacsony sűrűségű plagioklász kristályok felúszásából (flotációjából) jöttek létre a holdi magmaóceánban.

Ezzel egyidejűleg ez a folyamat diktálja a bolygó első légkörének kialakulását. Az illékony elemek és vegyületek (pl. víz, szén-dioxid, nitrogén), amelyek kezdetben oldott állapotban voltak az olvadt köpenyben, felszabadulnak – vagyis "kigázosodnak" –, ahogy az olvadék kristályosodik és oldhatóságuk csökken. Ennek az ősi légkörnek az összetétele és tömege tehát a magmaóceán evolúciójának közvetlen következménye, ami pedig mélyreható hatással van a bolygó korai klímájára és a lakhatóság lehetőségére. Ennek a korai differenciációnak az öröksége geológiai időskálákon is fennmaradhat; egyes kutatók szerint a ma megfigyelhető mélyköpeny-struktúrák, mint például a Föld köpenyének alján található Nagy, Alacsony Nyírásihullám-sebességű Tartományok (LLSVP-k), a földi magmaóceán utolsó maradékai lehetnek. Így a magmaóceán megszilárdulásának kvantitatív megértése nem csupán egy bolygó múltjának vizsgálata, hanem szükséges lépés a mai állapotának magyarázatához is.

1.2 A Kristályosodási Front mint Mozgó Határfelületi Probléma

A magmaóceán megszilárdulása nem egy pillanatszerű, az egész bolygóra kiterjedő esemény, amely tökéletes egyensúlyi körülmények között zajlik. Inkább egy dinamikus folyamat, amelyet a bolygó felszínéről az űrbe irányuló folyamatos hőveszteség hajt. Ahogy a bolygó energiát sugároz ki, a magma lehűl, és egy kristályosodási front – egy határfelület, amely elválasztja az újonnan képződött szilárd kumulátumhalmazt a megmaradt folyékony magmától – halad át a köpenyen. Egy ilyen terjedő fázisváltó határfelület vizsgálata a matematikai fizika klasszikus témája, amelyet mozgó határfelületi problémaként ismerünk. A leghíresebb példa a Stefan-probléma, amely a jég vízzé olvadását írja le. A magmaóceán megszilárdulása a Stefan-probléma egy rendkívül összetett, többkomponensű, geokémiai analógjaként fogható fel.

Ennek a kristályosodási frontnak a dinamikáját fizikai és kémiai folyamatok szorosan összekapcsolt kölcsönhatása szabályozza. A hőt mind a szilárd, mind a folyékony fázison keresztül kell szállítani vezetéssel és konvekcióval. Ezzel egyidejűleg a kémiai komponensek diffundálnak az olvadékban a haladó front felé vagy attól távolodva. A frakcionált kristályosodás folyamata, amely során bizonyos ásványok előnyben részesítve épülnek be a szilárd fázisba, kémiai "reakcióként" vagy nyelőként működik a határfelületen, folyamatosan megváltoztatva a szomszédos folyadék összetételét.

Ez a rendszer egy kritikus nemlineáris csatolást mutat, amely összetettségének elsődleges forrása és az egyszerű megoldás fő akadálya. A kristályosodási front hőmérséklete nem egy rögzített érték, hanem az olvadék likvidusz-hőmérséklete határozza meg. A likvidusz-hőmérséklet pedig érzékeny függvénye a helyi nyomásnak és, ami döntő fontosságú, az olvadék kémiai összetételének, beleértve az oldott illóanyagok koncentrációját is. Ezért a hőmérsékleti mező határozza meg a kristályosodás sebességét, de maga a kristályosodási folyamat megváltoztatja a kémiai mezőt, ami aztán visszahatva módosítja a határfelületen lévő termikus peremfeltételt. Ez a visszacsatolási hurok egységes kezelést tesz szükségessé. A termikus problémát (hőtranszport) és a kémiai problémát (anyagtranszport) nem lehet egymástól függetlenül megoldani; azokat csatolt parciális differenciálegyenletek rendszerként kell megoldani. A probléma ezen alapvető természete, mint csatolt, nemlineáris transzportjelenségek problémája, irányítja a jelentés teljes matematikai megközelítését.

1.3 A Probléma Formális Meghatározása

Tekintettel a magmaóceán megszilárdulásának alapvető fontosságára, a kristályosodási front dinamikájának kvantitatív, előrejelző modellje kiemelt tudományos értékkel bír. Míg a nagyléptékű numerikus szimulációk hatékony eszközök, számításigényesek lehetnek, és elfedhetik a paraméterek közötti alapvető fizikai kapcsolatokat. Egy analitikus megoldás, még egy egyszerűsített rendszerre is, mély betekintést nyújthat a skálatörvényekbe és a különböző fizikai folyamatok kölcsönhatásába.

A jelentés központi célkitűzése egy olyan zárt alakú analitikus megoldás levezetése és elemzése, amely leírja egy kristályosodási front egydimenziós terjedését, amelyet a s(t) pozíciója jelöl, egy félig végtelen bolygói magmaóceánba. A modellnek kellően kifinomultnak kell lennie ahhoz, hogy megragadja a folyamat kulcsfontosságú fizikáját. Konkrétan a modellnek figyelembe kell vennie a következőket:

1. Hőtranszport: A hő diffúziója mind a szilárd kumulátum köpenyben, mind a folyékony magmaóceánban.

2. Anyagtranszport: A fő szilikátalkotó elemek és a kulcsfontosságú illékony komponensek többkomponensű kémiai diffúziója a folyékony fázisban.

3. Fázisváltási reakció: A frakcionált kristályosodás folyamata a határfelületen, amely nyelőként működik a kémiai komponensek számára és forrásként a látens hő számára.

4. Bolygói hűlés: A bolygó hűlési sebességének hatása, amelyet egy termikus peremfeltétel képvisel, amely az egész folyamatot hajtja.

5. Illóanyag-hatások: Az illóanyag-koncentráció explicit hatása a fázisátalakulás termodinamikájára, különösen a likvidusz-hőmérséklet csökkenésére és a telítettség által kiváltott kristályosodás lehetőségére.

A kívánt megoldásnak a rendszer evolúciójának önkonzisztens leírását kell nyújtania, megadva a front helyzetét s(t), a hőmérsékletprofilokat a szilárd és folyékony fázisban, T_s(z,t) és T_l(z,t), valamint a releváns komponensek kémiai koncentrációprofiljait az olvadékban, C_i(z,t). Egy ilyen megoldás hatékony elméleti eszközt jelentene az elsődleges kéregképződés, a köpeny-differenciáció és a kőzetbolygókon történő illóanyag-kigázosodás történetének kvantitatív előrejelzéseinek javítására.

2. Szekció: Egy Többkomponensű Kristályosodó Rendszer Irányító Egyenletei

A kristályosodási front dinamikájára vonatkozó analitikus megoldás levezetéséhez először az 1. szekcióban leírt fizikai és geokémiai elveket kell szigorú matematikai keretrendszerbe ültetni. Ez a szekció fejleszti ki az irányító parciális differenciálegyenletek, határfelületi feltételek és peremfeltételek teljes rendszerét, amelyek leírják egy megszilárduló magmaóceán csatolt termikus és kémiai evolúcióját.

2.1 A Tartomány és a Koordináta-rendszer

Egy egyszerűsített egydimenziós geometriát veszünk figyelembe, amely megragadja a függőleges hő- és anyagtranszport lényeges fizikáját. Ez egy gyakori és hatékony kiindulópont a bolygóméretű folyamatok modellezéséhez, ahol a radiális gradiensek dominálnak. A magmaóceán megszilárdulása általában alulról felfelé halad, ami a nyomásfüggő olvadási görbe (szolidusz) meredekebb lejtésének következménye a magma adiabatikus hőmérsékleti profiljához képest. Ahogy az egész köpeny hűl, először a legnagyobb nyomáson, azaz az alján kezd kristályosodni. Ezért a tartományunkat egy olyan koordináta-rendszerben definiáljuk, ahol z=0 a mag-köpeny határt (CMB) jelöli, és a z felfelé növekszik a bolygó felszíne felé. A rendszert bármely t időpontban két régióra osztja a mozgó kristályosodási front a z=s(t) pozícióban:

• Szilárd fázis (Kumulátum köpeny): A 0 \le z < s(t) régió a megszilárdult ásványegyüttesből áll.

• Folyékony fázis (Magmaóceán): A z > s(t) régió az olvadt szilikátmagmából áll.

Az s(t) függvény a megszilárdult köpeny vastagságát jelenti, és a megoldás ismeretlen része. A feladat a hőmérséklet- és összetételmezők megoldása ezeken a tartományokon belül, miközben egyidejűleg meghatározzuk a s(t) határfelület evolúcióját is.

2.2 A Hőmérsékleti Mező (Hőtranszport)

A hőmérséklet evolúcióját mind a szilárd, mind a folyékony fázisban a hődiffúziós egyenlet szabályozza. Bár a magmaóceánok erőteljesen konvektívek, ami közel adiabatikus hőmérsékleti profilt eredményez a folyadék nagy részében, a diffúzió dominálja a hőtranszportot a határfelületek melletti termikus határrétegekben. Modellünk ezekre a diffúzív folyamatokra összpontosít, amelyek végső soron sebességmeghatározóak a határfelületi dinamikában. Az irányító egyenletek a szilárd, T_s, és a folyékony, T_l, fázis hőmérsékletére a következők:

• Szilárd fázis (0 \le z < s(t)): \frac{\partial T_s}{\partial t} = \kappa_s \frac{\partial^2 T_s}{\partial z^2} (2.1. egyenlet)

• Folyékony fázis (z > s(t)): \frac{\partial T_l}{\partial t} = \kappa_l \frac{\partial^2 T_l}{\partial z^2} (2.2. egyenlet)

Itt \kappa_s és \kappa_l a szilárd és folyékony köpeny hőmérséklet-vezetési tényezője (termikus diffuzivitása). A termikus diffuzivitás definíciója \kappa = k/(\rho c_p), ahol k a hővezető képesség, \rho a sűrűség, és c_p a fajhőkapacitás állandó nyomáson. Geológiai alkalmazásokban a termikus diffuzivitás tipikus értéke 10^{-6} \, \text{m}^2/\text{s} nagyságrendű. Ezek az egyenletek leírják, hogyan változik a hőmérséklet egy térbeli pontban a pontba be- vagy onnan kiáramló nettó hőáram következtében.

2.3 A Kémiai Mező (Anyagtranszport)

A magmaóceán egy összetett kémiai oldat. Az analitikus kezelhetőség érdekében ezt a valóságot egyszerűsítenünk kell, miközben megőrizzük a magmás differenciáció lényeges vonásait. Míg egy teljes modell a köpeny pirolitjára releváns teljes FCMAS (FeO-CaO-MgO-Al_2O_3-SiO_2) rendszert is figyelembe vehetné, a problémát egy általános n-komponensű rendszerre fogalmazzuk meg. A végső levezetéshez ezt specializálhatjuk egy ternér rendszerre (pl. olivin, piroxén és plagioklász komponenseket reprezentálva) és egy illékony komponensre.

Az egyes i kémiai komponens C_i koncentrációjának evolúcióját a folyékony fázisban egy reakció-diffúzió egyenlet írja le. A folyadék belsejében, távol a kristályosodó határfelülettől, feltételezzük, hogy nincsenek kémiai reakciók. Az elsődleges "reakció" – a kristályosodás – kizárólag a mozgó határfelületen, z=s(t)-nél történik. Ezért a reakciótag az irányító egyenletben a folyadék belsejében nulla. A koncentráció változása bármely ponton a folyadékban kizárólag a kémiai diffúziónak köszönhető:

• Folyékony fázis (z > s(t)): $$ \frac{\partial C_i}{\partial t} = D_i \frac{\partial^2 C_i}{\partial z^2} \quad \text{ahol } i=1,2,...,n $$ (2.3. egyenlet)

Itt D_i az i komponens kémiai diffúziós együtthatója (diffuzivitása) a szilikátolvadékban. A szilikátolvadékokban a kémiai diffuzivitások jellemzően sokkal kisebbek, mint a termikus diffuzivitások, gyakran 10^{-12} \, \text{m}^2/\text{s} nagyságrendűek. Ez a különbség a transzportsebességekben a geokémiai rendszerek egyik kulcsfontosságú jellemzője, ami a termikus határrétegeknél jóval vékonyabb kémiai határrétegek kialakulásához vezet.

2.4 Határfelületi Feltételek z=s(t)-nél (Az Általánosított Stefan-feltételek)

A fázisváltozás fizikája teljes egészében a mozgó határfelületen, z=s(t)-nél alkalmazott feltételekben rejlik. Ezek a feltételek összekapcsolják a termikus és kémiai mezőket, és szabályozzák a front mozgását.

2.4.1 Termodinamikai Egyensúly

A határfelületen a szilárd és folyékony fázisok helyi termodinamikai egyensúlyban vannak. Ez azt jelenti, hogy a hőmérséklet folytonos a határfelületen, és egyenlő az olvadék likvidusz-hőmérsékletével az adott nyomáson és összetételben.

T_s(s(t),t) = T_l(s(t),t) = T_{liq}(P, C_1, C_2,..., C_n, C_v) (2.4. egyenlet)

Ez a probléma legkritikusabb nemlineáris csatolási feltétele. A T_{liq} likvidusz-hőmérséklet a nyomás (P), a fő szilikátkomponensek koncentrációinak (C_i) és az oldott illóanyagok koncentrációjának (C_v) komplex függvénye. Ez a feltétel kimondja, hogy amint az olvadék összetétele a határfelületen a frakcionált kristályosodás miatt megváltozik, magának a kristályosodási hőmérsékletnek is ennek megfelelően kell változnia.

2.4.2 Az Illóanyag-telítettség Hatása

Az illóanyagok, mint a H_2O és a CO_2, döntő szerepet játszanak a magmás folyamatokban. Jellemzően "inkompatibilis" elemek, ami azt jelenti, hogy inkább az olvadékban maradnak, mintsem beépülnének a közönséges köpenyásványok kristályszerkezetébe. Az oldott illóanyagok jelenléte általában csökkenti a szilikátolvadék likvidusz-hőmérsékletét. Ezt a függést explicit módon modellezhetjük a likvidusz-függvényben, például egy lineáris közelítéssel alacsony koncentrációkra: T_{liq} \approx T_{liq,dry} - \gamma C_v, ahol \gamma egy állandó.

Kulcsfontosságú jelenség következik be, amikor egy illóanyag, C_v, koncentrációja az olvadékban eléri a telítettségi határát, C_{sat}(P,T). Ekkor egy különálló, illóanyagokban gazdag fluid fázis (gáz vagy szuperkritikus fluidum) válik ki az olvadékból. Ez a kigázosodási folyamat drámaian megnövelheti a megmaradó szilikátolvadék likvidusz-hőmérsékletét, potenciálisan egy gyors kristályosodási rohamot váltva ki. Ez a hatás a modell azon képességének sarokköve, hogy magyarázatot adjon a nemlineáris dinamikára, és beépítésre kerül a T_{liq} funkcionális formájába.

2.4.3 Hőenergia-mérleg (Stefan-feltétel)

A határfelület mozgását a hőelvonás hajtja, ami magában foglalja mind a folyadék hűlését (szenzibilis hő), mind a fázisváltozás során felszabaduló hőt (látens hő). A Stefan-feltétel az energiamegmaradás tétele a határfelületen. A sebesség, amellyel a fúzió látens hője, L, felszabadul, egyenlő kell legyen a határfelülettől távozó nettó hőárammal.

$$ \rho L \frac{ds}{dt} = \left. k_s \frac{\partial T_s}{\partial z} \right|{z=s(t)^{-}} - \left. k_l \frac{\partial T_l}{\partial z} \right|{z=s(t)^{+}} $$ (2.5. egyenlet)

Itt \rho a sűrűség (az egyszerűség kedvéért mindkét fázisra azonosnak feltételezve), L a fúzió látens hője tömegegységenként, k_s és k_l pedig a szilárd és folyékony fázis hővezető képessége. A ds/dt tag a kristályosodási front sebessége. Ez az egyenlet közvetlenül összekapcsolja a front sebességét a két oldalán lévő hőmérséklet-gradiensekkel.

2.4.4 Tömegmérleg (Kémiai Megoszlási Feltétel)

Ahogy a front halad előre, kémiai komponenseket épít be a folyadékból a szilárd fázisba. Ennek a beépülésnek a hatékonyságát a K_{D,i} ásvány-olvadék megoszlási hányados írja le, amelyet az i komponens szilárd fázisbeli (C_{i,s}) és folyadékfázisbeli (C_{i,l}) koncentrációjának arányaként definiálunk egyensúlyban: K_{D,i} = C_{i,s}/C_{i,l}.

Minden i komponensre a tömegmegmaradás a határfelületen megköveteli, hogy az a sebesség, amellyel az adott komponens tömege "elvész" a folyadékból, egyensúlyban legyen az adott komponens diffúzív áramával a tömbmagmából a határfelület felé. Ez egy sor kémiai ugrási feltételt eredményez:

$$ (C_{i,l} - C_{i,s}) \frac{ds}{dt} = -D_i \left. \frac{\partial C_{i,l}}{\partial z} \right|_{z=s(t)^{+}} $$ (2.6. egyenlet)

A megoszlási hányados definícióját behelyettesítve ez a következővé válik:

$$ C_{i,l}(s(t), t) (1-K_{D,i}) \frac{ds}{dt} = -D_i \left. \frac{\partial C_{i,l}}{\partial z} \right|_{z=s(t)^{+}} \quad \text{ahol } i=1,2,...,n $$ (2.7. egyenlet)

Ez egy kulcsfontosságú egyenletrendszer. Egy "inkompatibilis" elem esetén (K_{D,i} < 1), a (1-K_{D,i}) tag pozitív, ami azt jelenti, hogy a komponenst a szilárd fázis visszautasítja. Ehhez pozitív koncentrációgradiensre (\partial C_{i,l}/\partial z > 0) van szükség a határfelületen, hogy a visszautasított komponenst eldiffundáltassa az olvadékba, ami annak feldúsulásához vezet egy határrétegben. Ezzel szemben egy "kompatibilis" elem esetén (K_{D,i} > 1) a gradiensnek negatívnak kell lennie, ami elszegényedéshez vezet. Ezek az egyenletek összekapcsolják a front sebességét, ds/dt-t, az olvadékban lévő összes kémiai komponens koncentrációgradiensével.

2.5 Kezdeti és Távoltéri Peremfeltételek

Az PDE-rendszer egyértelmű megoldásának eléréséhez meg kell határoznunk a rendszer állapotát a folyamat kezdetén (kezdeti feltételek) és a tartomány határain (peremfeltételek).

Kezdeti feltételek (t=0-nál):

• A megszilárdulás a köpeny alján kezdődik: s(0) = 0 (2.8. egyenlet)

• A magmaóceán kezdetben egyenletes, likvidusz feletti T_0 hőmérsékleten van: T_l(z,0) = T_0 \quad \text{ahol } z>0 (2.9. egyenlet)

• A magmaóceán kezdetben kémiailag homogén, C_{i,0} összetételű: C_i(z,0) = C_{i,0} \quad \text{ahol } z>0 (2.10. egyenlet)

Peremfeltételek (t>0-ra):

• A mag-köpeny határon (z=0) a hőmérsékletet a mag hűlése szabályozza. Az egyszerűség kedvéért feltételezhetünk egy rögzített T_{CMB} hőmérsékletet, vagy egy lassan változó időfüggvényt, T_{CMB}(t): T_s(0,t) = T_{CMB}(t) (2.11. egyenlet)

• A távoltérben, mélyen a magmaóceánban (z \rightarrow \infty), a hőmérséklet és az összetétel a kezdeti, zavartalan értékén marad: \lim_{z\to\infty} T_l(z,t) = T_0 (2.12. egyenlet) \lim_{z\to\infty} C_i(z,t) = C_{i,0} (2.13. egyenlet)

A (2.1)–(2.13) egyenletrendszer teljes matematikai leírást ad az idealizált, egydimenziós magmaóceán-megszilárdulási problémára.

2.6 Termokémiai Paraméterek a Magmaóceán Modellezéséhez

Az irányító egyenletek megoldása számos fizikai és kémiai paraméter értékétől függ. Az alábbi táblázat egy szilikát magmaóceánra vonatkozó reprezentatív értékeket gyűjt össze, amelyeket kísérleti petrológiai és geofizikai modellezésből származtattunk, és amelyek az analitikus megoldást egy reális bolygói kontextusba helyezik.

1. Táblázat: Reprezentatív Termokémiai Paraméterek a Magmaóceán Modellezéséhez

Paraméter Szimbólum Érték/Tartomány Egység Forrás

Termikus tulajdonságok

Termikus diffuzivitás \kappa 1 \times 10^{-6} \text{m}^2/\text{s}

Fúzió látens hője L 4 \times 10^5 \text{J/kg}

Fajhőkapacitás c_p 1000-1200 \text{J/(kg}\cdot\text{K)}

Folyadék sűrűsége \rho_l 3500 \text{kg/m}^3

Kémiai tulajdonságok

Kémiai diffuzivitás (olvadék) D_i 10^{-12} - 10^{-10} \text{m}^2/\text{s}

Olvadék viszkozitása (alacsony nyomáson) \eta 0.01-10 \text{Pa}\cdot\text{s}

Megoszlási hányadosok (K_D = C_{szilárd}/C_{folyadék})

Olivin (Fe-Mg) K_{D, \text{Fe-Mg}}^{\text{ol}} 0.32 dimenzió nélküli

Gránát (Fe-Mg) K_{D, \text{Fe-Mg}}^{\text{grt}} 0.48 dimenzió nélküli

Bridgmanit (Ca) K_{D, \text{Ca}}^{\text{brg}} ~0.3 dimenzió nélküli

Illóanyag tulajdonságok

H_2O oldhatóság (bazalt, 1 GPa) C_{sat, H_2O} 5-7 tömeg%

CO_2 oldhatóság (bazalt, 1 GPa) C_{sat, CO_2} 0.1-0.5 tömeg%

Ez a táblázat biztosítja a modell parametrizálásához szükséges bemeneti adatokat. A termikus és kémiai diffuzivitások közötti jelentős (hat nagyságrendű) különbség azonnal azt sugallja, hogy a kristályosodási fronton lévő kémiai határrétegek sokkal élesebbek és vékonyabbak lesznek, mint a termikus határrétegek. A megoszlási hányadosok értékei határozzák meg a kémiai frakcionáció irányát és nagyságát, ahol az egynél kisebb értékek (mint például az Fe-Mg az olivinben) azt jelzik, hogy az elem inkompatibilis és a maradék olvadékban fog feldúsulni.

3. Szekció: Az Analitikus Megoldás Levezetése

A 2. szekcióban megfogalmazott, csatolt, nemlineáris parciális differenciálegyenletekből álló, mozgó határfelülettel rendelkező rendszer komoly matematikai kihívást jelent. Az általános esetre nem ismert közvetlen, egzakt analitikus megoldás. Ez gyakori jellemzője az ilyen egyenletekkel leírt összetett fizikai rendszereknek, beleértve a Stefan-problémát és a nemlineáris reakció-diffúzió rendszereket is. A továbbhaladás útja tehát nem egy nyers erővel történő megoldási kísérlet, hanem egy sor jól bevált analitikus technika alkalmazása, amelyek a problémát kezelhető formára egyszerűsítik, miközben megőrzik annak lényeges fizikai viselkedését. Ez a szekció részletezi ezt a többlépcsős levezetést, kezdve a dimenziómentesítéssel, folytatva a hasonlósági transzformációval, és egy perturbációs megoldásban csúcsosodva ki.

3.1 Az Egzakt Megoldás Kihívása és az Egyszerűsítés Indoklása

A közvetlen analitikus megoldás elsődleges akadályai a problémában rejlő nemlinearitások. Két fő forrásuk van:

1. A mozgó határfelület: A probléma tartománya maga is az idő függvénye, s(t). Ez azt jelenti, hogy a peremfeltételeket egy olyan helyen alkalmazzuk, amely a megoldás része, ami a szabad vagy mozgó határfelületi problémák meghatározó jellemzője, és alapvetően nehezebbé teszi őket a rögzített tartományú problémáknál.

2. Az összetételfüggő határfelületi hőmérséklet: A határfelületen érvényes termodinamikai feltétel, T(s(t),t) = T_{liq}(C_1, C_2,...), közvetlenül és nemlineáris módon kapcsolja össze a termikus és kémiai mezőket. A kristályosodás hőmérséklete magának a kristályosodási folyamatnak a kimenetelétől függ.

E kihívások fényében a "zárt alakú analitikus megoldás" keresése stratégiai megközelítést igényel. A cél egy tökéletes, mindent átfogó képlet megtalálása helyett egy olyan matematikailag átlátható modell levezetése, amely megvilágítja a mögöttes fizikát és a skálatörvényeket. A diffúzió által dominált, ilyen jellegű problémákra a leghatékonyabb technika a hasonlósági transzformáció, amely azon az elgondoláson alapul, hogy a rendszer térbeli szerkezete idővel önmagához hasonló módon fejlődik. Ez a módszer képezi analitikus stratégiánk sarokkövét. A probléma egy új koordináta-rendszerbe való transzformálásával a parciális differenciálegyenletekről a közönséges differenciálegyenletekre csökkenthetjük a komplexitást, amelyeket lényegesen könnyebb megoldani.

3.2 Dimenziómentesítés

Bármely szisztematikus elemzés első lépése az irányító egyenletek dimenziómentesítése. Ez a folyamat a problémát dimenzió nélküli változókban fogalmazza újra, aminek két fő előnye van. Először is, egyszerűsíti az egyenleteket a független fizikai paraméterek számának csökkentésével. Másodszor, feltárja azokat az alapvető dimenzió nélküli csoportokat, amelyek a rendszer viselkedését szabályozzák, mély fizikai betekintést nyújtva a versenyző folyamatokba.

A következő dimenzió nélküli változókat vezetjük be, ahol a tildével (~) jelölt változók dimenzió nélküliek:

• Tér: \tilde{z} = z/L_{ref}

• Idő: \tilde{t} = t/t_{ref}

• Hőmérséklet: \tilde{T} = (T - T_{ref,1}) / (T_{ref,2} - T_{ref,1})

• Koncentráció: \tilde{C}_i = (C_i - C_{i,ref,1}) / (C_{i,ref,2} - C_{i,ref,1})

A megfelelő referencia-skálák kiválasztása kulcsfontosságú. Egy természetes időskála a termikus diffúziós idő egy karakterisztikus L_{ref} hossz-skálán, tehát t_{ref} = L_{ref}^2/\kappa_l. A hőmérsékleti skálát a kezdeti magma-hőmérséklet és a CMB-hőmérséklet alapján választhatjuk meg, pl. T_{ref,1} = T_{CMB} és T_{ref,2} = T_0. Ezeket behelyettesítve az irányító egyenletekbe (2.1-2.13) és átrendezve egy dimenzió nélküli rendszert kapunk, ahol a viselkedést kevés dimenzió nélküli paraméter szabályozza.

3.2.1 Kulcsfontosságú Dimenzió Nélküli Csoportok

A dimenziómentesítés folyamata feltárja a következő kritikus dimenzió nélküli számokat, amelyek a magmaóceán megszilárdulásának dinamikáját szabályozzák.

2. Táblázat: Kulcsfontosságú Dimenzió Nélküli Csoportok a Magmaóceán Megszilárdulásában

Dimenzió Nélküli Csoport Szimbólum Képlet Fizikai Jelentés

Stefan-szám Ste \frac{c_p(T_{liq}-T_{CMB})}{L} A fázisváltozás során felszabaduló látens hő és a szilárd anyag hőmérsékletének megváltoztatásához szükséges szenzibilis hő aránya. A kristályosodás termikus hatását szabályozza. A kis Ste érték azt jelenti, hogy a látens hő csekély perturbációt jelent a teljes hőmérsékleti költségvetésben.

Lewis-szám Le_i \frac{\kappa_l}{D_i} A termikus diffuzivitás és a kémiai diffuzivitás aránya. A termikus és kémiai határrétegek relatív vastagságát szabályozza. Mivel a magmák esetében Le_i \gg 1, a kémiai határrétegek sokkal vékonyabbak, mint a termikusak.

Megoszlási szám \Pi_i 1-K_{D,i} Az i komponens kémiai frakcionációjának hatékonyságát jelenti. \Pi_i \approx 1 a magasan inkompatibilis elemeknél; \Pi_i < 0 a kompatibilis elemeknél.

Illóanyag-befolyás paraméter \Psi \frac{\gamma C_{sat}}{T_0 - T_{CMB}} Egy újonnan definiált paraméter, amely az illóanyag-csökkenés miatti maximális lehetséges likvidusz-hőmérséklet változást reprezentálja, a rendszer teljes hőmérséklet-különbségével skálázva. Számszerűsíti az illóanyagok potenciális hatását a rendszer termodinamikájára.

Diffuzivitás arány \tilde{\kappa} \frac{\kappa_s}{\kappa_l} A szilárd és folyékony fázis termikus diffuzivitásának aránya. Jellemzően 1-hez közeli.

Hővezető-képesség arány \tilde{k} \frac{k_s}{k_l} A szilárd és folyékony fázis hővezető-képességének aránya.

Ezek a dimenzió nélküli csoportok hatékony lencsét biztosítanak a probléma vizsgálatához. Például bármely két magmaóceán dinamikája, függetlenül azok specifikus bolygói kontextusától, azonos lesz, ha ezeknek a dimenzió nélküli számoknak az értékei megegyeznek. Analitikus megoldásunkat ezekben az alapvető paraméterekben fogjuk kifejezni.

3.3 A Hasonlósági Transzformáció

A diffúzió által vezérelt mozgó határfelületi problémák gyakran megengednek egy speciális megoldásosztályt, az úgynevezett hasonlósági megoldásokat. Az alapelv az, hogy a hőmérséklet és a koncentráció térbeli profiljai, megfelelően skálázva, idővel ugyanazt az alakot tartják meg. A megoldás bármely t időpontban úgy néz ki, mint a megoldás egy másik időpontban lévő nyújtott vagy összenyomott változata. Ezt a viselkedést egyetlen hasonlósági változó, \eta, bevezetésével ragadjuk meg, amely egyesíti a teret és az időt.

Az egydimenziós diffúziós egyenlethez a helyes hasonlósági változó : \eta = \frac{z}{2\sqrt{\kappa_l t}} (3.1. egyenlet)

Ezután olyan megoldást keresünk, ahol a dimenzió nélküli hőmérséklet \tilde{T} és koncentrációk \tilde{C}_i csak \eta függvényei, nem pedig z és t független függvényei: \tilde{T}(z,t) = \Theta(\eta) (3.2. egyenlet) \tilde{C}_i(z,t) = \Gamma_i(\eta) (3.3. egyenlet)

Ahhoz, hogy ez a transzformáció működjön, a mozgó határfelületnek, s(t)-nek, is követnie kell ugyanezt a skálázást. Ez azt jelenti, hogy a front pozíciójának a következő alakúnak kell lennie: s(t) = 2\lambda\sqrt{\kappa_l t} (3.4. egyenlet) ahol \lambda egy dimenzió nélküli állandó, amely a front normalizált sebességét képviseli. A \lambda értéke ismeretlen, és a megoldás részeként kell meghatározni.

A láncszabály alkalmazásával a differenciálásra (pl. \partial/\partial t = -(\eta/2t) \cdot d/d\eta és \partial/\partial z = (1/2\sqrt{\kappa_l t}) \cdot d/d\eta), az eredeti parciális differenciálegyenlet-rendszert (PDE) (z,t)-ben átalakíthatjuk egy csatolt, nemlineáris közönséges differenciálegyenlet-rendszerré (ODE) az egyetlen \eta változóban.

Például a hőegyenlet a folyékony fázisban (2.2. egyenlet) a következővé válik: -2\eta \frac{d\Theta_l}{d\eta} = \frac{d^2\Theta_l}{d\eta^2} (3.5. egyenlet)

És a kémiai diffúziós egyenlet (2.3. egyenlet) a következővé válik: -2\eta \frac{d\Gamma_i}{d\eta} = \frac{1}{Le_i} \frac{d^2\Gamma_i}{d\eta^2} (3.6. egyenlet)

A határfelületi és peremfeltételek hasonlóképpen átalakulnak a \Theta(\eta) és \Gamma_i(\eta) függvényekre és azok deriváltjaira vonatkozó feltételekké \eta=\lambda-nál és \eta=0, \infty-nél. Az eredmény egy rendkívül nemlineáris, kétpontos peremérték-probléma egy ODE-rendszerre. Bár még mindig összetett, ez jelentős egyszerűsítés az eredeti PDE-rendszerhez képest.

3.4 A Transzformált ODE-rendszer Aszimptotikus Megoldása

Az ODE-rendszer a határfelületi feltételek miatt továbbra is nemlineáris. Az analitikus továbblépés érdekében egy reguláris perturbációs módszert, más néven aszimptotikus sorfejtést alkalmazunk. Ez a technika akkor alkalmazható, ha a probléma egyik dimenzió nélküli paramétere nagyon kicsi. Számos geológiai rendszerben, így a magmaóceánokban is, a Stefan-szám (Ste) kicsi, mivel a fúzió látens hője elhanyagolható a köpeny több száz fokos lehűlésével járó, nagy mennyiségű szenzibilis hőhöz képest.

Ezért a Ste-t egy kis \epsilon paraméterként kezeljük, és a megoldást a Ste hatványsoraként keressük. Az ismeretlen \Theta(\eta), \Gamma_i(\eta) függvényeket és az ismeretlen \lambda állandót a következőképpen fejtjük ki: \Theta(\eta) = \Theta^{(0)}(\eta) + \text{Ste} \cdot \Theta^{(1)}(\eta) + O(\text{Ste}^2) (3.7. egyenlet) \Gamma_i(\eta) = \Gamma_i^{(0)}(\eta) + \text{Ste} \cdot \Gamma_i^{(1)}(\eta) + O(\text{Ste}^2) (3.8. egyenlet) \lambda = \lambda^{(0)} + \text{Ste} \cdot \lambda^{(1)} + O(\text{Ste}^2) (3.9. egyenlet)

Ezeket a sorfejtéseket behelyettesítjük a teljes ODE-rendszerbe és a peremfeltételekbe. Ezután összegyűjtjük az azonos Ste hatványú tagokat, és egymástól függetlenül nullával tesszük őket egyenlővé. Ez az eljárás az egyetlen, nehéz nemlineáris problémát egyszerűbb, lineáris peremérték-problémák hierarchiájára bontja, amelyeket egymás után meg lehet oldani.

3.4.1 Nulladrendű Megoldás (O(\text{Ste}^0))

A Ste-t nem tartalmazó tagok összegyűjtése adja a nulladrendű problémát. Ez fizikailag egy nulla látens hővel (L=0) rendelkező rendszernek felel meg. Az irányító ODE-k lineárisak, és általános megoldásaik a hibafüggvény (erf) és a komplementer hibafüggvény (erfc) segítségével fejezhetők ki, amelyek az egyszerű diffúziós egyenlet alapvető megoldásai.

A nulladrendű hőmérséklet általános megoldása a folyadékban, \Theta_l^{(0)}(\eta), például a következő alakú lesz: \Theta_l^{(0)}(\eta) = A_0 + B_0 \cdot \text{erf}(\eta) (3.10. egyenlet) ahol A_0 és B_0 integrálási állandók. Hasonló formák léteznek a szilárd fázis hőmérsékletére és a kémiai koncentrációkra (a hibafüggvény argumentumai Le_i-vel skálázva).

A peremfeltételek és a határfelületi feltételek alkalmazása \eta=\lambda^{(0)}-nál egy transzcendens algebrai egyenletrendszert eredményez az ismeretlen integrálási állandókra és \lambda^{(0)}-ra. A T_{liq} likvidusz-hőmérséklet függvény ezen a ponton lép be. Még nulladrendben is a probléma csatolt és nemlineáris, ha a T_{liq} függ az összetételtől. A további analitikus haladás érdekében a likvidusz-függvényt linearizáljuk a távoltéri összetétel, C_{i,0} körül. Ez egy lineáris algebrai egyenletrendszert eredményez, amelyet meg lehet oldani, hogy explicit kifejezéseket találjunk \lambda^{(0)}-ra és az integrálási állandókra. Ez adja a front sebességének és a termikus/kémiai profiloknak a vezető rendű közelítését.

3.4.2 Elsőrendű Megoldás (O(\text{Ste}^1))

Ezután összegyűjtjük az összes \text{Ste}^1-gyel arányos tagot. Ez egy új lineáris ODE-rendszert ad az elsőrendű korrekciós tagokra (\Theta^{(1)}, \Gamma_i^{(1)}). A kulcsfontosságú különbség az, hogy ezek az egyenletek most inhomogének. Az ezen egyenletekben szereplő forrástagok a most megtalált nulladrendű megoldásoktól függenek.

Például az elsőrendű hőenergia-mérleg (Stefan-feltétel) most már expliciten tartalmazni fogja a látens hő tagot, és ezt fogjuk használni a front sebességének elsőrendű korrekciója, \lambda^{(1)} megoldására. A megoldási folyamat hasonló: megoldjuk a lineáris inhomogén ODE-ket és alkalmazzuk a peremfeltételeket az új integrálási állandók és \lambda^{(1)} meghatározásához.

A végső analitikus megoldás, amely a Stefan-számban első rendig helyes, a nullad- és elsőrendű megoldások összege: s(t) \approx 2(\lambda^{(0)} + \text{Ste} \cdot \lambda^{(1)})\sqrt{\kappa_l t} (3.11. egyenlet)

Ez a kifejezés képviseli a problémafelvetésben keresett "zárt alakú analitikus megoldást". Ez egy közelítő megoldás, de olyan, amelyet első elvekből vezettünk le, és amely explicit, analitikus kapcsolatokat biztosít a front sebessége és a rendszer összes kulcsfontosságú fizikai paramétere között, amint azt a dimenzió nélküli csoportok megragadják. Az algebrai rendszerek és megoldásaik részletes levezetése hosszadalmas, és az A és B függelékekben található. A következő szekciók ennek a megoldásnak a fizikai jelentését és következményeit elemzik.

4. Szekció: A Kristályosodási Front Dinamikájának Elemzése

A 3. szekcióban levezetett matematikai megoldás, bár absztrakt, rengeteg információt tartalmaz egy megszilárduló magmaóceán fizikai viselkedéséről. Ez a szekció értelmezi ezt a megoldást a kristályosodási front dinamikájának elemzéséhez, fókuszálva annak sebességére, az illóanyag-telítettség kritikus szerepére, a magma ebből következő kémiai differenciációjára, valamint a kristályosodás és a légkör fejlődése közötti átfogó visszacsatolásra.

4.1 A Kristályosodási Front Sebessége

Az analitikus levezetés elsődleges eredménye a kristályosodási front helyzetének kifejezése az idő függvényében: s(t) = 2\lambda\sqrt{\kappa_l t} (4.1. egyenlet)

Ez az egyenlet felfedi, hogy a megszilárdult köpeny vastagsága arányos az idő négyzetgyökével. Ez egy klasszikus skálatörvény a diffúzió által szabályozott folyamatokra, ahol az előrehaladás sebessége idővel lelassul, ahogy a front távolabb kerül a hőnyelőtől (a magtól), és a termikus/kémiai gradiensek ellaposodnak.

A lényegi fizika a \lambda dimenzió nélküli állandóban van beágyazva, amelyet a perturbációs analízis során határoztunk meg. A \lambda-ra levezetett kifejezés (konkrétan a vezető rendű \lambda^{(0)} tag) felfedi, hogyan függ a front sebessége a rendszer kulcsfontosságú paramétereitől. Sematikusan a megoldás azt mutatja, hogy \lambda nagyobb (azaz a kristályosodás gyorsabb), amikor:

• A hőmérsékleti gradiens meredekebb: A sebesség arányos a forró magmaóceán (T_0) és a hűvösebb határfelület (T_{CMB}) közötti kezdeti hőmérséklet-különbséggel. A nagyobb hőmérséklet-különbség nagyobb hőáramot hajt, ami gyorsabb hűléshez és kristályosodáshoz vezet. Ez képviseli a megszilárdulás elsődleges hajtóerejét.

• A kezdeti összetétel kedvez a korai kristályosodásnak: A likvidusz-hőmérséklet, T_{liq}, függ a kezdeti összetételtől, C_{i,0}. Azok az összetételek, amelyeknek magasabb a likvidusz-hőmérséklete, egy adott termikus állapotban korábban és gyorsabban kezdenek kristályosodni.

• A frakcionáció kevésbé hatékony az inkompatibilis elemekre: Az inkompatibilis elemek (alacsony K_{D,i}-vel) felhalmozódása a front előtti határrétegben csökkenti a helyi likvidusz-hőmérsékletet, ami a kohászatban "konstitucionális túlhűtésként" ismert jelenség. Ez fékként hat a kristályosodásra. A modell azt mutatja, hogy \lambda csökken, ahogy az inkompatibilis elemek koncentrációja növekszik, számszerűsítve ezt az önkorlátozó hatást.

• Az illóanyag-tartalom alacsonyabb: Az oldott illóanyagok jelenléte szintén csökkenti a likvidusz-hőmérsékletet. A modell azt mutatja, hogy a magasabb kezdeti illóanyag-tartalom, C_{v,0}, kisebb \lambda-t és így lassabb kristályosodást eredményez, mivel a magma alacsonyabb hőmérsékleten is folyékony maradhat.

A \lambda analitikus formája lehetővé teszi a kvantitatív érzékenységi elemzést, pontosan megmutatva, hogy mennyivel változik a kristályosodási sebesség például a köpeny kezdeti víztartalmának vagy a bolygó hűlési sebességének változására.

4.2 Az Illóanyag-telítettség Szerepe mint Kiváltó Mechanizmus

A modell hatékony keretet biztosít a nemlineáris viselkedések vizsgálatára, különösen az illóanyag-telítettség hatásaira. A koncentrációprofilokra, \Gamma_i(\eta), adott analitikus megoldás előrejelzi az inkompatibilis elemek – beleértve az olyan illóanyagokat, mint a víz – felhalmozódását a haladó front előtti diffúzív határrétegben.

Egy illóanyag koncentrációja a határfelületen, C_v(s(t),t), idővel növekedni fog a kristályosodás előrehaladtával. Ez egy kritikus küszöbjelenséghez vezet. Egy bizonyos t_{sat} időpontban az illóanyag koncentrációja a határfelületen elérheti a telítettségi határát az adott nyomásra és hőmérsékletre, C_{sat}(P,T). Ez az esemény alapvetően megváltoztatja a termodinamikát a határfelületen. Ahogy az illékony komponens elkezd kiválni az olvadékból egy különálló fluid fázis képzésére, a megmaradó, most már illóanyagban szegény szilikátolvadék likvidusz-hőmérséklete élesen megnőhet.

Ezt a folyamatot analitikusan modellezhetjük a keretrendszerünkön belül.

1. A telítettség előrejelzése: A C_v(s(t),t) megoldását használjuk a t_{sat} időpont meghatározására, amelynél C_v(s(t_{sat}),t_{sat}) = C_{sat}.

2. A peremfeltétel megváltozása: t > t_{sat} esetén a határfelületi feltételben (2.4. egyenlet) használt T_{liq} likvidusz-hőmérséklet függvényt módosítjuk, hogy tükrözze a kigázosodott olvadék új, magasabb likviduszát. Az olvadék a határfelületen most a telítettségi feltételnél van pufferelve.

3. Új frontsebesség: A problémát újra megoldjuk ezzel az új peremfeltétellel. Az eredmény egy új, nagyobb érték a frontsebesség paraméterére, \lambda_{new}.

Ez a sorozat azt mutatja, hogy a kristályosodási front, egy potenciálisan hosszú, lassú, egyenletes növekedési periódus után, hirtelen, nemlineáris gyorsuláson mehet keresztül az illóanyag-telítettség elérésekor. Ez egy robusztus, első elveken alapuló fizikai mechanizmust biztosít az olyan késői stádiumú eseményekhez, amelyeket például a robbanásos vulkáni kitörések kiváltó okaként javasolnak. A gyorsulással járó gyors kristályosodás nagy mennyiségű illóanyagot taszítana ki a megszilárduló rétegből, potenciálisan hatalmas kigázosodási eseményeket hajtva, amelyek döntő fontosságúak egy bolygó ősi légkörének kialakulásában. Ez az átmenet egy lassú, diffúzió által vezérelt rezsimből egy gyors, telítettség által hajtott rezsimbe a modell egyik kulcsfontosságú prediktív eredménye.

4.3 A Maradék Olvadék Kémiai Evolúciója és a Magmás Differenciáció

A koncentrációprofilokra, \Gamma_i(\eta), adott analitikus megoldás a frakcionált kristályosodás és annak következményeinek közvetlen, kvantitatív modelljét nyújtja. A megoldás expliciten mutatja egy kémiai határréteg kialakulását a front előtt, ahol az inkompatibilis elemek (K_{D,i} < 1) feldúsulnak, a kompatibilis elemek (K_{D,i} > 1) pedig elszegényednek a távoltéri magmaösszetételhez, C_{i,0}-hoz képest. Ez a folyamat a magmás differenciáció lényege.

A megmaradó folyékony magma tömbösszetételének evolúcióját követve, ahogy a megszilárdulási front s(t) halad előre, előre jelezhetjük az egész rendszer kémiai pályáját. Ez közvetlenül analóg a reaktív olvadéktranszport modellezésével, ahol az olvadék összetétele fejlődik, miközben kölcsönhatásba lép egy szilárd mátrixszal.

Különösen fontos alkalmazás a magasan fejlett, geokémiailag dúsult rezervoárok kialakulásának modellezése. Ahogy a magmaóceán a kristályosodásának utolsó stádiumaihoz közelít (pl. >90% szilárd), a megmaradt kis mennyiségű olvadék rendkívül feldúsul minden olyan inkompatibilis elemben, amelyet szisztematikusan kizártak a kristályosodó szilárd anyagokból. Ide tartoznak olyan elemek, mint a kálium (K), a ritkaföldfémek (REE), a foszfor (P), az urán (U) és a tórium (Th), valamint az illóanyagok.

Modellünk kvantitatív előrejelzést adhat ennek a végső maradék olvadéknak az összetételére. A megszilárdult frakciót, s(t)/D_{MO} (ahol D_{MO} a magmaóceán teljes mélysége), egy 0.95-höz hasonló értékre állítva kiszámíthatjuk a megmaradó 5% folyadék megfelelő tömbösszetételét. Ez hatékony elméleti eszközt biztosít olyan rezervoárok eredetének megértéséhez, mint a holdi KREEP, egy geokémiai komponens, amelyről ismert, hogy ezekben az elemekben magasan dúsult, és feltételezhetően a holdi magmaóceán utolsó maradékait képviseli. A modell tesztelheti, hogy egy adott kezdeti holdi összetételből származó egyszerű frakcionált kristályosodási folyamat elegendő-e egy KREEP-szerű maradék folyadék létrehozásához.

4.4 A Kristályosodás-Redox-Atmoszféra Visszacsatolási Hurok

A probléma holisztikus szemléletéből egy kulcsfontosságú felismerés bontakozik ki: a bolygó hűlési sebessége nem egy független, külsőleg rákényszerített paraméter. Ehelyett endogén módon kapcsolódik magához a kristályosodási folyamathoz egy összetett visszacsatolási hurkon keresztül, amely magában foglalja a magma redoxállapotát és a kigázosodott légkör összetételét.

Ennek a visszacsatolásnak a sorrendje a következő:

1. Kristályosodás és Redox Evolúció: A fő köpenyásványok, különösen az olivin és a piroxén kristályosodása magában foglalja a vas megoszlását a vas(II) (Fe^{2+}) és vas(III) (Fe^{3+}) oxidációs állapotai között. Kísérleti bizonyítékok azt mutatják, hogy a Fe^{2+} előnyben részesítve épül be a kristályszerkezetbe a Fe^{3+}-tal szemben. Ennek eredményeként a frakcionált kristályosodás szisztematikusan eltávolítja a redukált vasat az olvadékból, aminek következtében a maradék folyékony magma Fe^{3+}/\Sigma Fe aránya idővel növekszik. A magmaóceán tehát fokozatosan oxidálódik, ahogy megszilárdul.

2. Redox Kontroll a Légköri Speciáción: A magma oxidációs állapota a felszínen beállítja a felette lévő légkör oxigénfugacitását (fO_2). Ez pedig szabályozza a kigázosodott illóanyagok kémiai speciációját. Egy redukált magmaóceán (alacsony fO_2) redukált légkört fog kigázosítani, amely gazdag olyan fajokban, mint a hidrogén (H_2) és a szén-monoxid (CO). Ahogy a magmaóceán oxidálódik, egyre oxidáltabb légkört kezd kigázosítani, amelyet a vízgőz (H_2O) és a szén-dioxid (CO_2) dominál.

3. Légköri Kontroll a Bolygói Hűlésen: A légkör összetétele mélyreható hatással van annak sugárzási tulajdonságaira és hőcsapda-képességére (üvegházhatás). A H_2-t és CO-t tartalmazó redukált légkörök általában kevésbé hatékonyak a hő űrbe történő kisugárzásában, mint az erős üvegházhatású gázokat, H_2O-t és CO_2-t tartalmazó oxidált légkörök. Ezért, ahogy a légkör átmenetet képez a redukáltról az oxidáltra, a bolygó hűtési képessége növekszik.

Ez a teljes visszacsatolási hurok beépíthető analitikus keretrendszerünkbe. A bolygó hűlési sebességét végső soron a felszíni hőáram határozza meg, amely a termikus probléma távoltéri peremfeltételeként szolgál. Ezt a felszíni hőáramot parametrizálhatjuk a megmaradó folyékony olvadék tömbösszetételének függvényeként, amely meghatározza a légkör összetételét. Például: \text{Felszíni hőáram} = F(fO_2(\text{Tömb } C_{Fe^{3+}/\Sigma Fe})) (4.2. egyenlet)

Azzal, hogy a termikus peremfeltételt a rendszer fejlődő állapotának függvényévé tesszük, a modell teljesen önkonzisztenssé válik. A hűlési sebesség már nem egy tetszőleges bemeneti adat, hanem a csatolt köpeny-légkör evolúciójának egy emergens tulajdonsága. Ez lehetővé teszi a modell számára, hogy megragadja az egyik leg alapvetőbb kölcsönhatást, amely egy bolygó átmenetét szabályozza egy olvadt gömbből egy stabil légkörrel rendelkező szilárd testté.

5. Szekció: Geokémiai és Bolygótudományi Következtetések

A kristályosodási front dinamikájának analitikus modellje, miután levezetésre és elemzésre került, most már alkalmazható a bolygótudomány kulcsfontosságú kérdéseinek megválaszolására. A mozgó határfelület mikroszkopikus fizikájának és egy bolygó makroszkopikus evolúciójának összekapcsolásával a modell kvantitatív előrejelzéseket kínál az ősi kéreg kialakulására, az illóanyagok kigázosodására és a hosszú távú bolygófejlődés kezdeti feltételeire. Ez a szekció ezeket a széles körű következményeket vizsgálja, összehasonlítva a modell előrejelzéseit a bolygói megfigyelésekkel, és kiemelve a jövőbeli kutatások irányait.

5.1 Előrejelzések az Ősi Kéreg Kialakulására

A magmaóceán teljes megszilárdulásához szükséges idő és a keletkező elsődleges kéreg végső vastagsága alapvető bolygói tulajdonságok. Modellünk közvetlen módszert kínál e mennyiségek becslésére. A frontsebesség, ds/dt, integrálásával t=0-tól addig, amíg a front eléri a magmaóceán tetejét (vagy amíg a megszilárdult frakció el nem ér egy kritikus küszöböt, pl. ~99%), kiszámíthatjuk a teljes megszilárdulási időskálát. A front végső pozíciója, s(t_{final}), megadja a köpeny kumulátumhalmazának vastagságát, amely az elsődleges kéreg alatt helyezkedik el.

Maga a kéreg jellege a felső magmaóceánban kristályosodó ásványok felhajtóerejétől függ. A modell előrejelzése a maradék olvadék kémiai evolúciójáról lehetővé teszi számunkra, hogy meghatározzuk, mikor válnak az alacsony sűrűségű ásványok, mint a plagioklász, a domináns kristályosodó fázissá.

• Alkalmazás a Holdra: A Hold a magmaóceánból képződött elsődleges kéreg kanonikus példáját nyújtja. A holdi felföldek anortozitból állnak, egy olyan kőzetből, amelyet a plagioklász ásvány dominál. Ezt a plagioklász-flotációs hipotézis magyarázza, amely szerint az alacsony sűrűségű plagioklász kristályok a holdi magmaóceán (LMO) tetejére úsztak, hogy vastag kérget képezzenek. Modellünket e megfigyeléssel szemben tesztelhetjük. Hold-specifikus paramétereket használva a kezdeti összetételre és a nyomás-hőmérséklet viszonyokra, szimulálhatjuk az LMO frakcionált kristályosodási sorozatát. Egy sikeres modellnek azt kellene előre jeleznie, hogy a mélységi mafikus ásványok (olivin, piroxén) kezdeti kristályosodási fázisa után a maradék olvadék annyira feldúsul alumíniumban és kalciumban, hogy a plagioklász válik a likvidusz fázissá, és alacsony sűrűsége miatt flotációs kérget képez. A modell előre jelezheti a termelt plagioklász teljes mennyiségét, lehetővé téve a közvetlen összehasonlítást a megfigyelt, ~30-40 km vastag holdi kéreggel. Fényt deríthet az eltérésekre is, például arra, hogy egyes numerikus modellek miért küzdenek a modern, GRAIL-korszakbeli kéregvastagság-mérések reprodukálásával, talán rávilágítva a plagioklász-termelés érzékenységére a kezdeti illóanyag-tartalomra vagy a hűlési sebességre, olyan tényezőkre, amelyeket modellünk expliciten kezel.

• Alkalmazás a Marsra: A bizonyítékok arra utalnak, hogy a Mars is átélt egy magmaóceán fázist, amely egy elsődleges kérget hozott volna létre. A marsi kéreg-dichotómia – az idősebb, vastagabb déli felföldek és a fiatalabb, vékonyabb északi alföldek közötti éles különbség – gyökerei ebben a korai korszakban keresendők. Modellünk első rendű becslést adhat az ősi marsi kéreg vastagságára és átlagos összetételére, egy alulról felfelé történő megszilárdulási forgatókönyvet feltételezve. Ez egy alapot biztosít, amelyhez képest a későbbi folyamatok, mint például az óriásbecsapódások vagy a köpeny-felborulás hatásai értékelhetők.

• Alkalmazás a Földre: Bár a Föld ősi kérge már rég elpusztult és újrahasznosult a több milliárd éves lemeztektonika által, kezdeti állapotának megértése kulcsfontosságú a hadei geodinamika modelljei számára. A Föld magmaóceánjának megszilárdulása határozta meg a köpeny kezdeti termikus és összetételi állapotát, befolyásolva, hogy a lemeztektonika mikor és hogyan kezdődhetett el. Modellünk előre jelezheti a magmaóceán fázis időtartamát a Földön, egy kritikus paramétert, amely meghatározza a bolygó kezdeti termikus állapotát és azt az időtartamot, amíg egy szilárd litoszféra kialakulhatott. A modell előrejelzi egy sűrű, vasban gazdag és hidroszus kumulátumréteg kialakulását is a megszilárduló köpeny tetején, amely gravitációsan instabil lett volna és hajlamos lett volna a felborulásra, ami a korai köpenydinamika kulcsfontosságú eseménye.

5.2 Az Illóanyag-kigázosodás és a Légkör Kialakulásának Modellezése

Egy bolygó másodlagos légkörének kialakulása elválaszthatatlanul kapcsolódik köpenyének megszilárdulásához. Az analitikus modell ennek a folyamatnak a kvantitatív történetét nyújtja. Bármely illékony komponens áramlása a magmaóceán felszínéről a légkörbe kiszámítható a modell által előrejelzett koncentrációgradiensekből. Ezt az áramlást a megszilárdulás időtartama alatt integrálva meghatározhatjuk az egyes illóanyagok (pl. H_2O, CO_2, N_2) teljes tömegét, amely a köpenyből a légkörbe kerül.

Még hatékonyabban, a kristályosodás-redox-atmoszféra visszacsatolási hurok (4.4. szekció) beépítésével a modell szimulálhatja a köpeny és a légkör együttes evolúcióját. A kristályosodás előrehaladtával a modell követi az olvadék változó redoxállapotát. Ez pedig minden lépésben meghatározza a kigázosodott illóanyagok speciációját. A modell így előre jelezheti a légköri összetétel időbeli alakulását, kezdve egy korai, redukált, H_2-ben és CO-ban domináló légkörrel, és áttérve egy oxidáltabb, H_2O-ban és CO_2-ben gazdag légkörre, ahogy az alatta lévő magmaóceán oxidálódik.

Ez közvetlen, fizikai kapcsolatot biztosít a megszilárdult köpeny százalékos aránya, valamint a légkör nyomása és összetétele között. Például a modell előre jelezheti, hogy hány földi óceánnyi víz gázosodna ki, mire a köpeny 50%-ban megszilárdul. Előre jelezheti a CO_2-dominált légkörre való átmenet időzítését is, ami a korai Föld és Vénusz evolúciójának kritikus lépése. Ezek az előrejelzések kulcsfontosságú bemeneti adatokat szolgáltatnak a korai Föld klímamodelljeihez és azoknak a körülményeknek a megértéséhez, amelyek kedvezőek lehettek vagy nem az élet eredetéhez.

5.3 A Modell Korlátai és a Jövőbeli Kutatások Irányai

Bár a jelentésben levezetett analitikus megoldás hatékony betekintést nyújt, egy sor egyszerűsítő feltételezésen alapul. E korlátok elismerése kulcsfontosságú az eredmények helyes értelmezéséhez és a jövőbeli kutatások ígéretes irányainak azonosításához.

• Egydimenziós geometria: A modell egydimenziós és laterális homogenitást feltételez. A valós magmaóceánok laterális hőmérséklet- és összetétel-változásokat tapasztaltak volna, különösen a heterogén becsapódások vagy a nagyléptékű konvektív mintázatok kialakulása miatt.

• A konvekció elhanyagolása a tömbmagmában: A bolygói magmaóceánok várhatóan erősen turbulensek, a Rayleigh-számok potenciálisan meghaladhatják a 10^{30}-t. Ez a konvekció hatékonyan keverné a magmaóceán nagy részét, a modellben elemzett diffúzív határrétegeket nagyon vékony régiókra korlátozva a kristályosodási front mellett. Megoldásunk ezért egy diffúzió-dominált végpontot képvisel. A valós dinamika a tömbfolyadékban zajló turbulens transzport és a határrétegekben zajló diffúzív transzport komplex kölcsönhatását foglalja magában. Modellünk előnye, hogy feloldja a fizikát ezen a kritikus határrétegen belül, ahol a fázisváltozás ténylegesen bekövetkezik.

• Egyszerűsített kémia és termodinamika: Egy egyszerűsített ternér rendszer és linearizált likvidusz-függvények használata szükséges kompromisszum az analitikus kezelhetőség érdekében. A valós magmák összetett, többfázisú, többkomponensű keverékek, rendkívül nemlineáris termodinamikai tulajdonságokkal. A modell megragadja a frakcionált kristályosodás elsőrendű hatásait, de nem veszi figyelembe a bonyolultabb jelenségeket, mint például az ásványreakciókat vagy a folyadék-immiszcibilitást.

• Éles határfelület feltételezése: A modell egy éles, végtelenül vékony határfelületet feltételez a szilárd és a folyékony fázis között. A valóságban valószínűleg létezik egy "kásás zóna" ("mushy zone") – egy részlegesen olvadt régió, amely szilárd kristályok hálózatát tartalmazza. A reológia és a transzporttulajdonságok ebben a kásás zónában összetettek, és befolyásolnák az általános dinamikát.

E korlátok ellenére az analitikus modell alapvető célt szolgál. Világos, fizikailag átlátható alapot biztosít a működő alapvető folyamatok megértéséhez. A jövőbeli kutatásoknak erre az analitikus alapra kell építeniük. A határfelületi dinamikára levezetett megoldás hatékony és számításilag takarékos szubgríd modellként használható a köpenykonvekció nagyléptékű numerikus szimulációin belül. Ez lehetővé tenné a nagyléptékű modellek számára, hogy beépítsék a kristályosodási front alapvető mikrofizikáját anélkül, hogy azt megfizethetetlen számítási költséggel kellene feloldaniuk. A jövőbeli elméleti munka kiterjeszthetné az analitikus megközelítést, például perturbációs módszerekkel bevonva a kásás zóna elsőrendű hatásait, vagy összekapcsolva az 1D modellt a légköri szökés egy parametrizált modelljével. Az egyszerűsített analitikus elmélet és a komplex numerikus szimuláció közötti szakadék áthidalásával tovább finomíthatjuk a bolygók születésének és fejlődésének megértését.

Függelékek

A. Függelék: A Transzformált ODE-rendszer Részletes Levezetése

Ez a függelék részletezi az irányító parciális differenciálegyenletek (PDE-k) és peremfeltételek transzformációját a (z,t) koordináta-rendszerből a \eta = z/(2\sqrt{\kappa_l t}) hasonlósági koordinátára.

1. Deriváltak transzformációja: A láncszabályt használva a t és z szerinti parciális deriváltak a következőképpen transzformálódnak: $$ \frac{\partial}{\partial t} = \frac{\partial\eta}{\partial t} \frac{d}{d\eta} = -\frac{z}{4t\sqrt{\kappa_l t}} \frac{d}{d\eta} = -\frac{\eta}{2t} \frac{d}{d\eta} $$ $$ \frac{\partial}{\partial z} = \frac{\partial\eta}{\partial z} \frac{d}{d\eta} = \frac{1}{2\sqrt{\kappa_l t}} \frac{d}{d\eta} $$ $$ \frac{\partial^2}{\partial z^2} = \frac{1}{4\kappa_l t} \frac{d^2}{d\eta^2} $$

2. Az irányító egyenletek transzformációja: Ezeket behelyettesítve az irányító egyenletekbe (2.1-2.3):

o Hőegyenlet (Folyadék, z>s(t) \Rightarrow \eta>\lambda): $$ \frac{\partial T_l}{\partial t} = \kappa_l \frac{\partial^2 T_l}{\partial z^2} \implies -\frac{\eta}{2t} \frac{d\Theta_l}{d\eta} = \kappa_l \left(\frac{1}{4\kappa_l t} \frac{d^2\Theta_l}{d\eta^2} \right) \implies \frac{d^2\Theta_l}{d\eta^2} + 2\eta \frac{d\Theta_l}{d\eta} = 0 $$

o Hőegyenlet (Szilárd, 0 \le z < s(t) \Rightarrow 0 \le \eta < \lambda): $$ \frac{\partial T_s}{\partial t} = \kappa_s \frac{\partial^2 T_s}{\partial z^2} \implies -\frac{\eta}{2t} \frac{d\Theta_s}{d\eta} = \kappa_s \left(\frac{1}{4\kappa_l t} \frac{d^2\Theta_s}{d\eta^2} \right) \implies \frac{d^2\Theta_s}{d\eta^2} + \frac{2\eta}{\tilde{\kappa}} \frac{d\Theta_s}{d\eta} = 0 $$ ahol \tilde{\kappa} = \kappa_s/\kappa_l.

o Kémiai diffúziós egyenlet (Folyadék, \eta > \lambda): $$ \frac{\partial C_i}{\partial t} = D_i \frac{\partial^2 C_i}{\partial z^2} \implies -\frac{\eta}{2t} \frac{d\Gamma_i}{d\eta} = D_i \left(\frac{1}{4\kappa_l t} \frac{d^2\Gamma_i}{d\eta^2} \right) \implies \frac{d^2\Gamma_i}{d\eta^2} + 2\eta Le_i \frac{d\Gamma_i}{d\eta} = 0 $$ ahol Le_i = \kappa_l/D_i.

3. Határfelületi és peremfeltételek transzformációja: A határokon lévő feltételek a következőképpen transzformálódnak:

o Határfelületi hőmérséklet (\eta = \lambda): \Theta_s(\lambda) = \Theta_l(\lambda) = \tilde{T}_{liq}(\Gamma_1(\lambda), \Gamma_2(\lambda),...)

o Stefan-feltétel (\eta = \lambda): $$ \rho L \frac{ds}{dt} = k_s \frac{\partial T_s}{\partial z}|{s(t)} - k_l \frac{\partial T_l}{\partial z}|{s(t)} $$ Mivel ds/dt = \lambda\sqrt{\kappa_l/t}, ez a következővé válik: $$ \rho L \lambda \sqrt{\frac{\kappa_l}{t}} = \frac{T_0-T_{CMB}}{2\sqrt{\kappa_l t}} \left(k_s \frac{d\Theta_s}{d\eta}|{\lambda} - k_l \frac{d\Theta_l}{d\eta}|{\lambda} \right) $$ Dimenzió nélküli formában, a \text{Ste} = L/(c_p(T_0-T_{CMB})) és k=\rho c_p \kappa felhasználásával: $$ 2\lambda\text{Ste} = \tilde{k} \frac{d\Theta_s}{d\eta}|{\lambda} - \frac{d\Theta_l}{d\eta}|{\lambda} $$

o Kémiai megoszlási feltétel (\eta = \lambda): C_{i,l}(1-K_{D,i})\frac{ds}{dt} = -D_i \frac{\partial C_{i,l}}{\partial z}|_{s(t)} $$ (\tilde{C}{i,0} + \Delta\tilde{C}i \Gamma_i(\lambda))(1-K{D,i}) \lambda \sqrt{\frac{\kappa_l}{t}} = -D_i \frac{\Delta\tilde{C}i}{2\sqrt{\kappa_l t}} \frac{d\Gamma_i}{d\eta}|{\lambda} $$ $$ 2\lambda Le_i (1-K{D,i}) \Gamma_i(\lambda) = -\frac{d\Gamma_i}{d\eta}|{\lambda} \quad (\text{feltételezve, hogy } \tilde{C}{i,0} \text{ a referencia}) $$

o Peremfeltételek: \Theta_s(0) = \tilde{T}_{CMB} \lim_{\eta\to\infty} \Theta_l(\eta) = \tilde{T}_0 \lim_{\eta\to\infty} \Gamma_i(\eta) = \tilde{C}_{i,0}

Ez befejezi a transzformációt egy csatolt ODE-rendszerre, peremfeltételekkel 0, \lambda és \infty-ben.

B. Függelék: A Perturbációs Egyenletek Lépésről Lépésre Történő Megoldása

Ez a függelék vázolja a perturbációs sorfejtésből származó nulladrendű probléma megoldását.

1. Nulladrendű egyenletek (O(\text{Ste}^0)): A Stefan-feltétel nulladrendben a következővé válik: $$ \tilde{k} \frac{d\Theta_s^{(0)}}{d\eta}|{\lambda^{(0)}} - \frac{d\Theta_l^{(0)}}{d\eta}|{\lambda^{(0)}} = 0 $$ Ez azt jelenti, hogy a hőáram folytonos a határfelületen, ahogy az a látens hő hiányában várható.

2. Általános megoldások: A lineáris ODE-k általános megoldásai: \Theta_s^{(0)}(\eta) = A_s + B_s \cdot \text{erf}(\eta/\sqrt{\tilde{\kappa}}) \Theta_l^{(0)}(\eta) = A_l + B_l \cdot \text{erf}(\eta) \Gamma_i^{(0)}(\eta) = A_i + B_i \cdot \text{erf}(\eta\sqrt{Le_i})

3. Peremfeltételek alkalmazása: \Theta_s^{(0)}(0) = \tilde{T}_{CMB} \implies A_s = \tilde{T}_{CMB} \lim_{\eta\to\infty} \Theta_l^{(0)}(\eta) = \tilde{T}_0 \implies A_l + B_l = \tilde{T}_0 \lim_{\eta\to\infty} \Gamma_i^{(0)}(\eta) = \tilde{C}_{i,0} \implies A_i + B_i = \tilde{C}_{i,0}

4. Határfelületi feltételek alkalmazása \eta = \lambda^{(0)}-nál: Ez egy algebrai egyenletrendszert eredményez a fennmaradó állandókra (B_s, A_l, B_l, A_i, B_i) és az ismeretlen \lambda^{(0)}-ra.

o Hőmérséklet-folytonosságból: \Theta_s^{(0)}(\lambda^{(0)}) = \Theta_l^{(0)}(\lambda^{(0)})

o Hőáram-folytonosságból: \tilde{k} \frac{B_s}{\sqrt{\pi\tilde{\kappa}}} e^{-(\lambda^{(0)})^2/\tilde{\kappa}} = \frac{B_l}{\sqrt{\pi}} e^{-(\lambda^{(0)})^2}

o Kémiai megoszlásból: 2\lambda^{(0)}Le_i(1-K_{D,i})\Gamma_i^{(0)}(\lambda^{(0)}) = -B_i \frac{2\sqrt{Le_i}}{\sqrt{\pi}} e^{-(\lambda^{(0)})^2 Le_i}

o Likvidusz-relációból: \Theta_l^{(0)}(\lambda^{(0)}) = \tilde{T}_{liq}(\Gamma_1^{(0)}(\lambda^{(0)}),...)

Ez a rendszer még mindig nemlineáris az exponenciális és hibafüggvények, valamint a likvidusz-reláció miatt. A teljesen analitikus megoldáshoz további egyszerűsítés (pl. a likvidusz-függvény linearizálása és nagy Lewis-számok feltételezése a kémiai feltétel egyszerűsítésére) szükséges. Ennek a rendszernek a megoldása adja a \lambda^{(0)} értékét, amely a kristályosodási front sebességének vezető rendű közelítését biztosítja. A \lambda^{(1)}-re vonatkozó elsőrendű problémát ezután ezekkel a nulladrendű megoldásokkal mint forrástagokkal fogalmazzuk meg, és hasonló, bár algebrailag összetettebb módon oldjuk meg.

C. Függelék: Fogalmak és Szimbólumok Jegyzéke

Szimbólum Meghatározás Egység

t Idő s

z Térbeli koordináta (mélység) m

s(t) A kristályosodási front pozíciója m

T_s, T_l Hőmérséklet a szilárd, folyékony fázisban K

C_i Az i komponens koncentrációja \text{kg/m}^3 vagy tömeg%

\kappa_s, \kappa_l Termikus diffuzivitás a szilárd, folyékony fázisban \text{m}^2/\text{s}

D_i Az i komponens kémiai diffuzivitása \text{m}^2/\text{s}

k_s, k_l Hővezető-képesség a szilárd, folyékony fázisban \text{W/(m}\cdot\text{K)}

\rho Sűrűség \text{kg/m}^3

c_p Fajhőkapacitás \text{J/(kg}\cdot\text{K)}

L Fúzió látens hője J/kg

K_{D,i} Megoszlási hányados az i komponensre dimenzió nélküli

T_{liq} Likvidusz-hőmérséklet K

C_{sat} Egy illóanyag telítettségi koncentrációja tömeg%

\eta Hasonlósági változó dimenzió nélküli

\lambda Dimenzió nélküli frontsebesség-állandó dimenzió nélküli

Ste Stefan-szám dimenzió nélküli

Le_i Lewis-szám az i komponensre dimenzió nélküli

\Theta Dimenzió nélküli hőmérséklet dimenzió nélküli

\Gamma_i Dimenzió nélküli koncentráció dimenzió nélküli

fO_2 Oxigénfugacitás Pa

LLSVP Nagy, Alacsony Nyírásihullám-sebességű Tartomány

KREEP Kálium (Potassium), Ritkaföldfémek (Rare Earth Elements), Foszfor (Phosphorus)

CMB Mag-köpeny határ (Core-Mantle Boundary)

LMO Holdi Magmaóceán (Lunar Magma Ocean)

Idézett munkák

1. Goodness Gracious, Great Balls of Fire! - How Magma Oceans Shaped Planetary Evolution, hozzáférés dátuma: július 30, 2025, https://carnegiescience.edu/events/goodness-gracious-great-balls-fire-how-magma-oceans-shaped-planetary-evolution

2. Magma oceans as a critical stage in the tectonic development of rocky planets - PMC, hozzáférés dátuma: július 30, 2025, https://pmc.ncbi.nlm.nih.gov/articles/PMC6189560/

3. Magma oceans as a critical stage in the tectonic development of rocky planets - Journals, hozzáférés dátuma: július 30, 2025, https://royalsocietypublishing.org/doi/10.1098/rsta.2018.0109

4. planetology.ethz.ch, hozzáférés dátuma: július 30, 2025, https://planetology.ethz.ch/research/evolution-earth-magma-ocean.html#:~:text=Magma%20oceans%20are%20thought%20to,of%20core%2D%20and%20atmosphere%20formation.

5. Evolution of Earth's magma ocean - Experimental Planetology | ETH Zurich, hozzáférés dátuma: július 30, 2025, https://planetology.ethz.ch/research/evolution-earth-magma-ocean.html

6. The magma ocean stage in the formation of rocky-terrestrial planets | ACT of ESA, hozzáférés dátuma: július 30, 2025, https://www.esa.int/gsp/ACT/coffee/2020-01-16-%20Magma Oceans/

7. Petrology and Physics of Magma Ocean Crystallization - NASA Technical Reports Server (NTRS), hozzáférés dátuma: július 30, 2025, https://ntrs.nasa.gov/citations/20050157877

8. Magma Oceans in the Inner Solar System | Annual Reviews, hozzáférés dátuma: július 30, 2025, https://www.annualreviews.org/content/journals/10.1146/annurev-earth-042711-105503

9. Fractional crystallization (geology) - Wikipedia, hozzáférés dátuma: július 30, 2025, https://en.wikipedia.org/wiki/Fractional_crystallization_(geology

10. 4.4: Partial Melting and Crystallization - Geosciences Libre Texts, hozzáférés dátuma: július 30, 2025, https://geo.libretexts.org/Bookshelves/Geology/Book%3A_An_Introduction_to_Geology_(Johnson_Affolter_Inkenbrandt_and_Mosher)/04%3A_Igneous_Processes_and_Volcanoes/4.04%3A_Partial_Melting_and_Crystallization

11. Experimental Fractional Crystallization of the Lunar Magma Ocean, hozzáférés dátuma: július 30, 2025, https://ntrs.nasa.gov/api/citations/20120003028/downloads/20120003028.pdf

12. Geodynamics | Magma oceans - EGU Blogs, hozzáférés dátuma: július 30, 2025, https://blogs.egu.eu/divisions/gd/2018/06/22/magma-oceans/

13. Is lunar magma ocean (LMO) gone with the wind? - Oxford Academic, hozzáférés dátuma: július 30, 2025, https://academic.oup.com/nsr/article/3/1/12/2460231

14. [2301.07505] Redox evolution of the crystallizing terrestrial magma ocean and its influence on atmosphere outgassing - arXiv, hozzáférés dátuma: július 30, 2025, https://arxiv.org/abs/2301.07505

15. Magma Ocean, Water, and the Early Atmosphere of Venus - NASA Technical Reports Server, hozzáférés dátuma: július 30, 2025, https://ntrs.nasa.gov/api/citations/20230011915/downloads/MWaySpace_SciRev_Magma_Reprint.pdf

16. Thermal evolution of an early magma ocean in interaction with the atmosphere CORE, hozzáférés dátuma: július 30, 2025, https://core.ac.uk/download/pdf/52710935.pdf

17. An ocean of magma formed early in Earth's history and it may still influence our planet today, study finds - Live Science, hozzáférés dátuma: július 30, 2025, https://www.livescience.com/planet-earth/an-ocean-of-magma-formed-early-in-earths-history-and-it-may-still-influence-our-planet-today-study-finds

18. The three stages of magma ocean cooling - NASA Technical Reports Server (NTRS), hozzáférés dátuma: július 30, 2025, https://ntrs.nasa.gov/citations/19920019374

19. On the cooling of a deep terrestrial magma ocean, hozzáférés dátuma: július 30, 2025, https://www.ipgp.fr/~samuel/henrilPGP/Publications_files/Monteux_etal_2016.pdf

20. Dynamic crystallization in magmas | Request PDF - ResearchGate, hozzáférés dátuma: július 30, 2025, https://www.researchgate.net/publication/314244758_Dynamic_crystallization_in_magmas

21. Dynamic crystallization in magmas - SOEST Hawaii, hozzáférés dátuma: július 30, 2025, https://www.soest.hawaii.edu/GG/FACULTY/JHAMMER/documents/Mollo_EMU_2017.pdf

22. Analytical and numerical solutions to the three-phase Stefan problem with simultaneous occurrences of melting, solidification, boiling, and condensation phenomena - arXiv, hozzáférés dátuma: július 30, 2025, https://arxiv.org/html/2503.06360v2

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

24. Diffusion in Melts and Magmas - NASA Technical Reports Server, hozzáférés dátuma: július 30, 2025, https://ntrs.nasa.gov/api/citations/20230006048/downloads/Diffusion%20in%20Melts_rev_clean_final.pdf

25. Imperfect fractional crystallization of the lunar magma ocean and formation of the lunar mantle: A realistic chemical approach - NASA Technical Reports Server (NTRS), hozzáférés dátuma: július 30, 2025, https://ntrs.nasa.gov/citations/19920019367

26. Onset of solid-state mantle convection and mixing during magma.... hozzáférés dátuma: július 30, 2025, https://d-nb.info/1204478759/34

27. Magma Oceans and Primordial Mantle Differentiation - ResearchGate, hozzáférés dátuma: július 30, 2025, https://www.researchgate.net/publication/323799632_Magma_Oceans_and_Primordial_Mantle_Differentiation

28. Magma ocean fractional crystallization and cumulate overturn in terrestrial planets: Implications for Mars, hozzáférés dátuma: július 30, 2025, http://lunar.earth.northwestern.edu/courses/438/tanton2003.pdf

29. 2.2: Geological Applications of the Diffusion Equation - Geosciences Libre Texts, hozzáférés dátuma: július 30, 2025, https://geo.libretexts.org/Courses/University_of_California_Davis/GEL_056%3A_Introduction_to_Geophysics/Geophysics_is_everywhere_in_geology.../02%3A_Diffusion_and_Darcy's_Law/2.02%3A_Geological_Applications_of_the_Diffusion_Equation

30. Investigating Magma Ocean Solidification on Earth Through Laser-Heated Diamond Anvil Cell Experiments - PubMed Central, hozzáférés dátuma: július 30, 2025, https://pmc.ncbi.nlm.nih.gov/articles/PMC8244043/

31. Looking inside the Earth - The Australian Museum, hozzáférés dátuma: július 30, 2025, https://australian.museum/learn/minerals/shaping-earth/looking-inside-the-earth/

32. Reaction-Diffusion Equations: A Comprehensive Guide - Number Analytics, hozzáférés dátuma: július 30, 2025, https://www.numberanalytics.com/blog/ultimate-guide-reaction-diffusion-equations

33. Reaction-diffusion system - Wikipedia, hozzáférés dátuma: július 30, 2025, https://en.wikipedia.org/wiki/Reaction%E2%80%93diffusion_system

34. Volatiles in Magmas (Wallace), hozzáférés dátuma: július 30, 2025, https://www.geo.mtu.edu/EHaz/ConvergentPlates_Class/wallace/Wallace%20&%20Anderson%20%281999%29.pdf

35. Late-stage volatile saturation as a potential trigger for explosive volcanic eruptions | Request PDF - ResearchGate, hozzáférés dátuma: július 30, 2025, https://www.researchgate.net/publication/292673946_Late-stage_volatile_saturation_as_a_potential_trigger_for_explosive_volcanic_eruptions

36. PRE-ERUPTIVE WATER CONTENTS IN LUNAR SILICIC MAGMAS: CONSTRAINTS FROM VOLATILE SATURATION AT MAGMA STORAGE DEPTHS. C. Sun and J., hozzáférés dátuma: július 30, 2025, https://www.hou.usra.edu/meetings/lpsc2024/pdf/2744.pdf

37. Analytical solution of Stefan-type problems - ResearchGate, hozzáférés dátuma: július 30, 2025, https://www.researchgate.net/publication/377075238_Analytical_solution_of_Stefan-type_problems

38. Insights into magma ocean dynamics from the transport properties of basaltic melt - PMC, hozzáférés dátuma: július 30, 2025, https://pmc.ncbi.nlm.nih.gov/articles/PMC9731987/

39. existence of analytic solutions for the classical stefan problem - Department of Mathematics, hozzáférés dátuma: július 30, 2025, https://math.vanderbilt.edu/simoneg/preprints/Stefan-Annalen.pdf

40. Analytical Solutions for a Nonlinear Reaction-diffusion Equation - AIP Publishing, hozzáférés dátuma: július 30, 2025, https://pubs.aip.org/aip/acp/article-pdf/1558/1/1771/11527001/1771_1_online.pdf

41. On the one dimensional Stefan problem - DiVA portal, hozzáférés dátuma: július 30, 2025, https://www.diva-portal.org/smash/get/diva2:647481/FULLTEXT01.pdf

42. Similarity solutions - MEK4300 Lecture notes - Professor Mikael Mortensen, hozzáférés dátuma: július 30, 2025, https://mikaem.github.io/MEK4300/content/chapter3/similarity.html

43. Perturbation methods, Physics 2400, hozzáférés dátuma: július 30, 2025, https://www.phys.uconn.edu/~rozman/Courses/P2400_20F/downloads/perturbation-methods.pdf

44. Perturbation Methods GM01 Dr. Helen J. Wilson Autumn Term 2008 - UCL, hozzáférés dátuma: július 30, 2025, https://www.ucl.ac.uk/~ucahhwi/GM01/lecturenotes.pdf

45. Perturbation Methods, hozzáférés dátuma: július 30, 2025, https://www.iust.ac.ir/files/fnst/ssadeghzadeh_52bb7/perturbation.pdf

46. Reactive Melt Transport in the Mantle and Geochemical Signatures of Mantle-derived Magmas | Journal of Petrology | Oxford Academic, hozzáférés dátuma: július 30, 2025, https://academic.oup.com/petrology/article/39/5/1039/1445087

47. From Binary Mixing to Magma Chamber Simulator: Geochemical Modeling of Assimilation in Magmatic Systems, hozzáférés dátuma: július 30, 2025, https://magma.geol.ucsb.edu/papers/Heinonen-etal_Mixing-to-MCS_Ch07.pdf

48. How Did Magma Oceans Evolve on Early Earth and Mars? - Eos.org, hozzáférés dátuma: július 30, 2025, https://eos.org/research-spotlights/how-did-magma-oceans-evolve-on-early-earth-and-mars

49. Vertically Resolved Magma Ocean-Protoatmosphere Evolution: H2, H2O, CO2, CH4, CO, O2, and N2 as Primary Absorbers - PMC - PubMed Central, hozzáférés dátuma: július 30, 2025, https://pmc.ncbi.nlm.nih.gov/articles/PMC7988593/

50. WHY DO LUNAR MAGMA OCEAN COOLING MODELS STRUGGLE TO REPRODUCE A GRAIL-ERA CRUSTAL THICKNESS? | Request PDF - ResearchGate, hozzáférés dátuma: július 30, 2025, https://www.researchgate.net/publication/365033699_WHY_DO_LUNAR_MAGMA_OCEAN_COOLING_MODELS_STRUGGLE_TO_REPRODUCE_A_GRAIL-ERA_CRUSTAL_THICKNESS

51. Water in the Lunar Mantle: Results from Magma Ocean Modeling - ResearchGate, hozzáférés dátuma: július 30, 2025, https://www.researchgate.net/publication/241303726_Water_in_the_Lunar_Mantle_Results_from_Magma_Ocean_Modeling

52. Magma matters: modelling the complex generation and transport of magma beneath volcanoes | by Research Features | Medium, hozzáférés dátuma: július 30, 2025, https://medium.com/@Research_Features/magma-matters-modelling-the-complex-generation-and-transport-of-magma-beneath-volcanoes-c06a12e9c0dd

Works cited

1. evolution: from cosmogenesis to biogenesis - IAEA, https://inis.iaea.org/records/j52g6-nak94/files/22048985.pdf?download=1 2. The Characterization of Reaction-Convection-Diffusion Processes by Travelling Waves - https ://ris.utwen te.nl, https://ris.utwente.nl/ws/files/6544087/Gilding96characterization.pdf 3. Josef Stefan and His Contributions to Heat Transfer - ResearchGate, https://www.researchgate.net/publication/267575099_Josef_Stefan_and_His_Contributions_to_Heat_Transfer 4. Angol-magyar geológiai szótár – English-Hungarian Geological Dictionary - ELTE, https://kazmer.web.elte.hu/pubs/Kazmer_1995_English-Hungarian_Geological_Dictionary.pdf 5. ALEX STREKEISEN-Anorthosite-, https://www.alexstrekeisen.it/english/pluto/anorthosite.php 6. Large low-shear-velocity provinces - Wikipedia, https://en.wikipedia.org/wiki/Large_low-shear-velocity_provinces 7. The role of volatile cycles on Earth - Refubium - Freie Universität Berlin, https://refubium.fu-berlin.de/bitstream/handle/fub188/46613/Thesis_Vulpius.pdf?sequence=3&isAllowed=y 8. Free Boundary Problems, Theory and Applications, Herakleion, Crete, Greece, 8-14 June 1997 - DTIC, https://apps.dtic.mil/sti/tr/pdf/ADA327885.pdf 9. Cobalt and marine redox evolution - Archimer, https://archimer.ifremer.fr/doc/00174/28575/26986.pdf 10. Cobalt and marine redox evolution - The People of Earth & Planetary Sciences - Yale University, https://people.earth.yale.edu/sites/default/files/files/Planavsky/31_%20Swanner%20et%20al%20EPSL%202014.pdf 11. A területi geokémiai kutatás elméleti és gyakorlati módszerei - Magyar Elektronikus Könyvtár, http://mek.oszk.hu/19000/19095/19095.pdf 12. BÁNYÁSZATI és KOHÁSZATI LAPOK - OMBKE, https://ombke.hu/wp-content/uploads/2025/07/BKL_2025-2.pdf 13. ESA - Mapping for the Moon - European Space Agency, https://www.esa.int/ESA_Multimedia/Images/2025/07/Mapping_for_the_Moon 14. Földtani közlöny, https://upload.wikimedia.org/wikipedia/commons/c/ce/F%C3%B6ldtani_k%C3%B6zl%C3%B6ny_%28IA_foyldtanikoyzlo127magya%29.pdf 15. oktanol-víz megoszlási hányados - vilaglex - Kémia, http://www.vilaglex.hu/Kemia/Html/Oktanol.htm 16. Mafic glass compositions: a record of magma storage conditions, mixing and ascent - PMC, https://pmc.ncbi.nlm.nih.gov/articles/PMC6335480/ 17. Plagioclase - Wikipedia, https://en.wikipedia.org/wiki/Plagioclase 18. KREEP - Wikipedia, https://en.wikipedia.org/wiki/KREEP 19. UNDERSTANDING VARIATION IN PARTITION COEFFICIENT, Kd, VALUES - Environmental Protection Agency (EPA), https://www.epa.gov/sites/default/files/2015-05/documents/402-r-99-004b.pdf 20. Wide Range Variation of δ³⁴S Isotopes in the Çayırhan Oil Shales: Multiple Sulfur Sources and Environmental Influences | MAS Journal of Applied Sciences, https://masjaps.com/index.php/mas/article/view/643

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