2025-11-21T14:19:15.500377

Electrical and Thermal Conductivity of Earth's Iron-enriched Basal Magma Ocean

Dragulet, Stixrude

The Earth's earliest magnetic field may have originated in a basal magma ocean, a layer of silicate melt surround the core that could have persisted for billions of years. Recent studies show that the electrical conductivity of liquid with a bulk silicate Earth composition exceeds 10000 S/m at basal magma ocean conditions, potentially surprising the threshold for dynamo activity. Over most of its history however, the basal magma ocean is more enriched in iron than the bulk silicate Earth, due to iron's incompatibility in the mineral assemblages of the lower mantle. Using ab-initio molecular dynamics calculations, we examine how iron content affects the silicate dynamo hypothesis. We investigate how the electrical conductivity of silicate liquid changes with iron enrichment, at pressures and temperatures relevant for Earth's basal magma ocean. We also compute the electronic contribution to the thermal conductivity , to evaluate convective instability of basal magma oceans. Finally, we apply our results to model the thermal and magnetic evolution of Earth's basal magma ocean over time.

academic

पृथ्वी के लौह-समृद्ध आधार मैग्मा महासागर की विद्युत और तापीय चालकता

मूल जानकारी

पेपर ID: 2508.00791
शीर्षक: Electrical and Thermal Conductivity of Earth's Iron-enriched Basal Magma Ocean
लेखक: Francis Dragulet, Lars Stixrude (UCLA)
वर्गीकरण: astro-ph.EP cond-mat.mtrl-sci physics.geo-ph
प्रकाशन समय: 25 अगस्त 2025 (ड्राफ्ट संस्करण)
प्रस्तुत पत्रिका: Proceedings of the National Academy of Sciences
पेपर लिंक: https://arxiv.org/abs/2508.00791

सारांश

पृथ्वी के प्रारंभिक चुंबकीय क्षेत्र की उत्पत्ति आधार मैग्मा महासागर (BMO) से हुई हो सकती है—यह पृथ्वी के कोर के चारों ओर एक सिलिकेट पिघली हुई परत है, जो संभवतः अरबों वर्षों तक मौजूद रही है। नवीनतम शोध से पता चलता है कि आधार मैग्मा महासागर की स्थितियों में, पृथ्वी की समग्र सिलिकेट संरचना वाली तरल की विद्युत चालकता 10⁴ S/m से अधिक है, जो डायनेमो गतिविधि की सीमा तक पहुंच सकती है। हालांकि, निचले मेंटल खनिज संयोजन में लौह की असंगतता के कारण, आधार मैग्मा महासागर अपने इतिहास के अधिकांश समय में पृथ्वी की समग्र सिलिकेट संरचना की तुलना में अधिक लौह-समृद्ध रहा है। यह अध्ययन प्रथम-सिद्धांत आणविक गतिविज्ञान गणना का उपयोग करके, लौह सामग्री कैसे सिलिकेट डायनेमो परिकल्पना को प्रभावित करती है, इसकी जांच करता है।

अनुसंधान पृष्ठभूमि और प्रेरणा

वैज्ञानिक प्रश्न

पृथ्वी के प्रारंभिक चुंबकीय क्षेत्र की उत्पत्ति का रहस्य: पृथ्वी का चुंबकीय क्षेत्र कम से कम 3.5 अरब वर्षों से सक्रिय है, जो पृथ्वी की आवास क्षमता को बनाए रखने के लिए महत्वपूर्ण है, लेकिन वर्तमान में बाहरी कोर के तरल लौह संवहन द्वारा संचालित डायनेमो का इतिहास 1 अरब वर्ष से कम हो सकता है
आधार मैग्मा महासागर डायनेमो परिकल्पना: BMO प्रारंभिक चुंबकीय क्षेत्र के संभावित स्रोत के रूप में, चुंबकीय रेनॉल्ड्स संख्या Rm > 40 की शर्त को पूरा करने की आवश्यकता है
लौह-समृद्धि प्रभाव अज्ञात: मौजूदा चालकता गणना मुख्य रूप से सरलीकृत प्रणालियों पर आधारित है, BMO क्रिस्टलीकरण प्रक्रिया में लौह-समृद्धि प्रभाव पर विचार नहीं करती है

अनुसंधान का महत्व

पृथ्वी के प्रारंभिक चुंबकीय क्षेत्र द्वारा वायुमंडल की सुरक्षा और वायुमंडल के नुकसान को रोकने के तंत्र की व्याख्या करना
चट्टानी बहिर्ग्रह ग्रहों के चुंबकीय क्षेत्र विकास के लिए सैद्धांतिक आधार प्रदान करना
पृथ्वी की आंतरिक संरचना और तापीय विकास की समझ को गहरा करना

मौजूदा विधि की सीमाएं

पिछली चालकता गणना MgO, MgSiO₃, SiO₂ जैसी सरलीकृत प्रणालियों पर केंद्रित थी
BMO क्रिस्टलीकरण प्रक्रिया में लौह-समृद्धि का व्यवस्थित अध्ययन अभाव
तापीय चालकता और विद्युत चालकता के युग्मन संबंध को पर्याप्त रूप से नहीं समझा गया है

मुख्य योगदान

पहली बार व्यवस्थित अध्ययन BMO विद्युत चालकता और तापीय चालकता पर लौह-समृद्धि के प्रभाव का
मात्रात्मक संबंध स्थापित करना: विद्युत चालकता और लौह सामग्री के बीच द्विघात संबंध, Mott-Ziman सिद्धांत के अनुरूप
महत्वपूर्ण तंत्र की खोज: लौह 3d इलेक्ट्रॉन और O 2p अवस्थाओं के संकरण से बनी विस्तृत ऊर्जा पट्टी विद्युत चालकता वृद्धि का मूल कारण है
विकास मॉडल का निर्माण: तापीय विकास मॉडल के साथ BMO डायनेमो के समय विकास की भविष्यवाणी करना
व्यवहार्यता का प्रमाण: लौह-समृद्धि सिलिकेट डायनेमो के जीवनकाल को 3.3 अरब वर्षों तक विस्तारित करती है

विधि विवरण

कार्य परिभाषा

विभिन्न लौह सामग्री (XFe = Fe/(Fe+Mg) = 0.12, 0.5, 1.0) वाली सिलिकेट तरल की विद्युत चालकता σel और तापीय चालकता kel की गणना BMO स्थितियों में (100-140 GPa, 4000-6000 K)।

गणना विधि आर्किटेक्चर

1. आणविक गतिविज्ञान सिमुलेशन

सैद्धांतिक आधार: घनत्व कार्यात्मक सिद्धांत (DFT) + PBEsol सन्निकटन + Hubbard U सुधार (U-J = 2.5 eV)
कार्यान्वयन मंच: VASP कोड, प्रक्षेपित संवर्धित समतल तरंग विधि
सिमुलेशन पैरामीटर:
- प्रणाली: 149 परमाणु, 6 तत्व (Mg, Fe, Si, O, Ca, Al)
- समय चरण: 1 fs, कुल समय: 10-15 ps
- ऊर्जा कटऑफ: 500 eV
- K-बिंदु नमूनाकरण: Gamma बिंदु
- समूह: NVT, Nosé-Hoover तापमान नियंत्रण

2. इलेक्ट्रॉन परिवहन गुण गणना

Chester-Thellung रूप में Kubo-Greenwood विधि का उपयोग:

विद्युत चालकता: $\sigma_{el}(\omega) = L_{11}(\omega)$

तापीय चालकता: $k_{el}(\omega) = \frac{1}{e^2T}\left(L_{22}(\omega) - \frac{L_{12}(\omega)L_{21}(\omega)}{L_{11}(\omega)}\right)$

जहां गतिविज्ञान गुणांक: $L_{\alpha\beta}(\omega) = \frac{2\pi e^2\hbar^2}{3m_e^2\omega\Omega}\sum_{i,j}(f_i-f_j)\delta(\epsilon_i-\epsilon_j-\hbar\omega)|\langle\psi_i|\nabla|\psi_j\rangle|^2(-1)^{\alpha+\beta}(\epsilon_i-\mu)^{\alpha-1}(\epsilon_j-\mu)^{\beta-1}$

3. स्पिन अवस्था उपचार

लौह के उच्च-स्पिन-निम्न-स्पिन संक्रमण पर विचार करते हुए, तापीय एकीकरण के माध्यम से मुक्त ऊर्जा अंतर की गणना: $\Delta F_{HS-LS} = \int_0^1 \langle\Delta U\rangle_\lambda d\lambda$

संतुलन विद्युत चालकता: $\frac{1}{\sigma_{el}^{eq}} = \frac{f_{eq}}{\sigma_{el}^{HS}} + \frac{1-f_{eq}}{\sigma_{el}^{LS}}$

तकनीकी नवाचार बिंदु

बहु-पैमाने युग्मन: परमाणु पैमाने DFT गणना और मैक्रोस्कोपिक तापीय विकास मॉडल का संयोजन
स्पिन अवस्था उपचार: उच्च दबाव में लौह के स्पिन संक्रमण प्रभाव को सटीकता से ध्यान में रखना
संरचना विकास: तत्व वितरण गुणांक के आधार पर BMO संरचना समय विकास
अरैखिक चरण आरेख: अरैखिक तरल चरण रेखा तापमान मॉडल का उपयोग, Fe-Mg प्रणाली को अधिक यथार्थवादी रूप से प्रतिबिंबित करना

प्रायोगिक सेटअप

गणना पैरामीटर

दबाव सीमा: 100-140 GPa (BMO संबंधित दबाव)
तापमान सीमा: 4000-6000 K (BMO संबंधित तापमान)
लौह सामग्री: XFe = 0.12 (जैतून की चट्टान), 0.5, 1.0 (लौह-समृद्ध अंतिम सदस्य)
अभिसरण मानदंड: ऊर्जा 3 meV/atom, दबाव 0.2 GPa

मूल्यांकन संकेतक

विद्युत चालकता: σel (S/m), डायनेमो सीमा (~10⁴ S/m) के साथ तुलना
तापीय चालकता: kel (W/m/K), संवहन स्थिरता का मूल्यांकन
चुंबकीय रेनॉल्ड्स संख्या: Rm = μ₀vlσ, डायनेमो निर्णायक (Rm > 40)

तुलना बेंचमार्क

शुद्ध MgSiO₃ प्रणाली की आयनिक तापीय चालकता (~4-5 W/m/K)
पृथ्वी के कोर की तरल लौह विद्युत चालकता और तापीय चालकता
Wiedemann-Franz नियम द्वारा भविष्यवाणी मान

प्रायोगिक परिणाम

मुख्य परिणाम

1. लौह सामग्री के साथ विद्युत चालकता में परिवर्तन

XFe = 0.12: σel ~ 2×10⁴ S/m (6000K, 100GPa)
XFe = 1.0: σel > 10⁵ S/m, लगभग 5 गुना वृद्धि
तापमान निर्भरता: σel तापमान बढ़ने के साथ बढ़ता है
दबाव निर्भरता: BMO दबाव सीमा में प्रभाव अपेक्षाकृत कम है

2. तापीय चालकता विकास

XFe = 0.12: kel ~ 5 W/m/K
XFe = 1.0: kel ~ 30 W/m/K
अभी भी पृथ्वी के कोर से कम: तरल लौह बाहरी कोर से बहुत कम (>100 W/m/K)

3. इलेक्ट्रॉन अवस्था घनत्व विश्लेषण

लौह 3d इलेक्ट्रॉन और O 2p अवस्थाओं का संकरण विस्तृत ऊर्जा पट्टी बनाता है
फर्मी ऊर्जा पर अवस्था घनत्व g(EF) XFe के साथ रैखिक रूप से बढ़ता है
σel ∝ g(EF)², Mott-Ziman सिद्धांत के अनुरूप

चुंबकीय क्षेत्र विकास सिमुलेशन

समय विकास परिणाम

स्थिर संरचना (XFe=0.12): डायनेमो 1.4 Gyr तक चलता है
लौह-समृद्धि पर विचार: डायनेमो 3.3 Gyr तक चलता है (1.2 Ga पहले तक)
महत्वपूर्ण मोटाई: डायनेमो को बनाए रखने के लिए न्यूनतम BMO मोटाई 200 km से घटकर 100 km हो गई

डायनेमो अंतराल विश्लेषण

शासन आरेख के माध्यम से निर्धारित:

BMO हमेशा तापीय संवहन डायनेमो अंतराल में रहता है
उच्च तापीय चालकता द्वारा डायनेमो को दबाने वाले अंतराल तक नहीं पहुंचता है
पृथ्वी के कोर की तुलना में, BMO में अधिक अनुकूल σ-k पैरामीटर संयोजन है

विलोपन प्रयोग

स्पिन अवस्था प्रभाव: निम्न-स्पिन अवस्था विद्युत चालकता उच्च-स्पिन अवस्था का लगभग 2 गुना है
वेग स्केलिंग: CIA स्केलिंग MLT स्केलिंग की तुलना में 1.3 Gyr कम डायनेमो जीवनकाल देता है
आयनिक योगदान: इलेक्ट्रॉन चालकता कुल चालकता का 70-90% है

निष्कर्ष और चर्चा

मुख्य निष्कर्ष

लौह-समृद्धि विद्युत चालकता को महत्वपूर्ण रूप से बढ़ाती है: XFe 0.12 से 1.0 तक बढ़ने पर, σel लगभग एक दशक बढ़ता है
डायनेमो जीवनकाल विस्तारित होता है: लौह-समृद्धि पर विचार करने के बाद, सिलिकेट डायनेमो 3.3 अरब वर्षों तक बना रह सकता है
संवहन अभी भी बनाए रखा जा सकता है: तापीय चालकता में वृद्धि के बावजूद, BMO संवहन अस्थिरता बनाए रखता है
प्रारंभिक चुंबकीय क्षेत्र व्यवहार्य है: पृथ्वी के प्रारंभिक चुंबकीय क्षेत्र के लिए एक उचित भौतिक तंत्र प्रदान करता है

सीमाएं

संरचना सरलीकरण: क्रिस्टलीकरण के साथ SiO₂ सामग्री में कमी पर विचार नहीं किया गया है
संवहन पैटर्न: संरचना संवहन के योगदान को शामिल नहीं किया गया है
विकिरण ऊष्मा स्थानांतरण: विकिरण तापीय चालकता योगदान को नजरअंदाज किया गया है
त्रि-आयामी प्रभाव: एक-आयामी तापीय विकास मॉडल पर आधारित

भविष्य की दिशाएं

अति-पृथ्वी अनुप्रयोग: उच्च दबाव में विद्युत चालकता गुणों का अध्ययन करना
संरचना संवहन: लौह अवसादन द्वारा संचालित संवहन प्रभाव को मापना
बहु-तत्व प्रणाली: Al, Ca जैसे तत्वों के प्रभाव पर विचार करना
प्रायोगिक सत्यापन: उच्च तापमान उच्च दबाव प्रयोग सत्यापन

गहन मूल्यांकन

शक्तियां

उन्नत विधि: प्रथम-सिद्धांत गणना और मैक्रोस्कोपिक मॉडल का संयोजन, बहु-पैमाने अनुसंधान
स्पष्ट भौतिकी: इलेक्ट्रॉन संरचना के दृष्टिकोण से चालकता तंत्र की व्याख्या
महत्वपूर्ण परिणाम: प्रारंभिक पृथ्वी चुंबकीय क्षेत्र की उत्पत्ति के लिए मात्रात्मक बाधाएं प्रदान करता है
पूर्ण मॉडल: स्पिन संक्रमण, संरचना विकास आदि कई भौतिक प्रक्रियाओं पर विचार करता है

कमियां

गणना लागत: DFT+U विधि के पैरामीटर चयन परिणामों को प्रभावित कर सकते हैं
मॉडल धारणाएं: एक-आयामी तापीय विकास मॉडल बहुत सरलीकृत है
प्रायोगिक अभाव: उच्च तापमान उच्च दबाव प्रायोगिक डेटा सत्यापन की कमी
अनिश्चितता: कुछ पैरामीटर (जैसे वितरण गुणांक) में काफी अनिश्चितता है

प्रभाव

शैक्षणिक मूल्य: पृथ्वी भौतिकी और ग्रह विज्ञान के लिए महत्वपूर्ण सैद्धांतिक आधार प्रदान करता है
अनुप्रयोग संभावनाएं: बहिर्ग्रह ग्रह चुंबकीय क्षेत्र अनुसंधान तक विस्तारित किया जा सकता है
विधि प्रदर्शन: कम्प्यूटेशनल पृथ्वी भौतिकी अनुसंधान प्रतिमान प्रदर्शित करता है
नीति महत्व: पृथ्वी की आवास क्षमता विकास की समझ बढ़ाता है

लागू परिदृश्य

प्रारंभिक पृथ्वी चुंबकीय क्षेत्र विकास अनुसंधान
चट्टानी बहिर्ग्रह ग्रह चुंबकीय क्षेत्र भविष्यवाणी
उच्च तापमान उच्च दबाव में सिलिकेट गुण अनुसंधान
पृथ्वी की आंतरिक संरचना और गतिविज्ञान मॉडलिंग

संदर्भ

यह पेपर पृथ्वी भौतिकी, कम्प्यूटेशनल भौतिकी और ग्रह विज्ञान क्षेत्र के महत्वपूर्ण साहित्य का हवाला देता है, जिसमें शामिल हैं:

Labrosse et al. (2007): BMO तापीय विकास शास्त्रीय मॉडल
Christensen & Aubert (2006): डायनेमो चुंबकीय रेनॉल्ड्स संख्या सीमा
Pozzo et al. (2012): पृथ्वी के कोर तापीय चालकता गणना
McDonough & Sun (1995): पृथ्वी जैतून की चट्टान संरचना मॉडल

यह पेपर कठोर प्रथम-सिद्धांत गणना और तापीय विकास मॉडलिंग के माध्यम से, पृथ्वी के प्रारंभिक चुंबकीय क्षेत्र की उत्पत्ति इस महत्वपूर्ण वैज्ञानिक समस्या के लिए मात्रात्मक भौतिक बाधाएं प्रदान करता है, जो कम्प्यूटेशनल पृथ्वी भौतिकी क्षेत्र में एक महत्वपूर्ण योगदान है।