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
  • প্রকাশনার সময়: ২০২৫ সালের ১৭ আগস্ট (খসড়া সংস্করণ)
  • জমা দেওয়ার জার্নাল: Proceedings of the National Academy of Sciences
  • পত্রের লিঙ্ক: https://arxiv.org/abs/2508.00791

সারসংক্ষেপ

পৃথিবীর প্রাথমিক চৌম্বক ক্ষেত্র ভিত্তি ম্যাগমা মহাসাগর (BMO) থেকে উৎপন্ন হতে পারে—এটি পৃথিবীর মূলের চারপাশে একটি সিলিকেট গলিত স্তর যা সম্ভবত কয়েক বিলিয়ন বছর ধরে বিদ্যমান ছিল। সর্বশেষ গবেষণা দেখায় যে ভিত্তি ম্যাগমা মহাসাগরের অবস্থায়, পৃথিবীর সামগ্রিক সিলিকেট সংমিশ্রণ সহ তরল বৈদ্যুতিক পরিবাহিতা ১০⁴ S/m অতিক্রম করে এবং সম্ভবত ডায়নামো কার্যকলাপের প্রান্তিক মূল্যে পৌঁছায়। তবে, নিম্ন ম্যান্টেল খনিজ সমন্বয়ে লোহার অসামঞ্জস্যতার কারণে, ভিত্তি ম্যাগমা মহাসাগর তার ইতিহাসের বেশিরভাগ সময়ে পৃথিবীর সামগ্রিক সিলিকেটের চেয়ে লোহায় সমৃদ্ধ ছিল। এই গবেষণা প্রথম নীতি আণবিক গতিশীলতা গণনা ব্যবহার করে, লোহার সামগ্রী কীভাবে সিলিকেট ডায়নামো অনুমান প্রভাবিত করে তা পরীক্ষা করে।

গবেষণার পটভূমি এবং প্রেরণা

বৈজ্ঞানিক প্রশ্ন

  1. পৃথিবীর প্রাথমিক চৌম্বক ক্ষেত্রের উৎসের রহস্য: পৃথিবীর চৌম্বক ক্ষেত্র কমপক্ষে ৩.৫ বিলিয়ন বছর ধরে সক্রিয় রয়েছে এবং পৃথিবীর বাসযোগ্যতা বজায় রাখার জন্য অত্যন্ত গুরুত্বপূর্ণ, কিন্তু বর্তমানে বাহ্যিক মূলের তরল লোহার প্রবাহ দ্বারা চালিত ডায়নামো সম্ভবত ১ বিলিয়ন বছরেরও কম বয়সী
  2. ভিত্তি ম্যাগমা মহাসাগর ডায়নামো অনুমান: BMO প্রাথমিক চৌম্বক ক্ষেত্রের সম্ভাব্য উৎস হিসাবে, চৌম্বক রেনল্ডস সংখ্যা Rm > 40 শর্ত পূরণ করা প্রয়োজন
  3. লোহা সমৃদ্ধকরণ প্রভাব অজানা: বিদ্যমান পরিবাহিতা গণনা প্রধানত সরলীকৃত সিস্টেমের উপর ভিত্তি করে, BMO স্ফটিকীকরণ প্রক্রিয়ায় লোহা সমৃদ্ধকরণ প্রভাব বিবেচনা করে না

গবেষণার গুরুত্ব

  • পৃথিবীর প্রাথমিক চৌম্বক ক্ষেত্র বায়ুমণ্ডল রক্ষা করে এবং বায়ুমণ্ডল হ্রাস প্রতিরোধ করার প্রক্রিয়া ব্যাখ্যা করা
  • শিলাময় বহিঃগ্রহের চৌম্বক ক্ষেত্র বিবর্তনের জন্য তাত্ত্বিক ভিত্তি প্রদান করা
  • পৃথিবীর অভ্যন্তরীণ কাঠামো এবং তাপীয় বিবর্তন সম্পর্কে গভীর বোঝাপড়া

বিদ্যমান পদ্ধতির সীমাবদ্ধতা

  • পূর্ববর্তী পরিবাহিতা গণনা MgO, MgSiO₃, SiO₂ ইত্যাদি সরলীকৃত সিস্টেমে কেন্দ্রীভূত
  • BMO স্ফটিকীকরণ প্রক্রিয়ায় লোহা সমৃদ্ধকরণের সিস্টেমেটিক গবেষণার অভাব
  • তাপীয় পরিবাহিতা এবং বৈদ্যুতিক পরিবাহিতার সংযুক্তি সম্পর্ক যথাযথভাবে অন্বেষণ করা হয়নি

মূল অবদান

  1. প্রথমবারের মতো সিস্টেমেটিক গবেষণা BMO-তে লোহা সমৃদ্ধকরণের বৈদ্যুতিক এবং তাপীয় পরিবাহিতার প্রভাব
  2. পরিমাণগত সম্পর্ক প্রতিষ্ঠা: বৈদ্যুতিক পরিবাহিতা এবং লোহার সামগ্রীর মধ্যে দ্বিঘাত সম্পর্ক, Mott-Ziman তত্ত্বের সাথে সামঞ্জস্যপূর্ণ
  3. মূল প্রক্রিয়া আবিষ্কার: লোহার 3d ইলেকট্রন এবং O 2p অবস্থার সংকর যা প্রশস্ত শক্তি ব্যান্ড গঠন করে, এটি পরিবাহিতা বৃদ্ধির মূল কারণ
  4. বিবর্তন মডেল নির্মাণ: তাপীয় বিবর্তন মডেলের সাথে মিলিয়ে BMO ডায়নামোর সময়ের বিবর্তন পূর্বাভাস দেওয়া
  5. সম্ভাব্যতা প্রমাণ: লোহা সমৃদ্ধকরণ সিলিকেট ডায়নামোর আয়ু ৩.৩ বিলিয়ন বছর পর্যন্ত উল্লেখযোগ্যভাবে বৃদ্ধি করে

পদ্ধতির বিস্তারিত বিবরণ

কাজের সংজ্ঞা

BMO অবস্থায় (১০০-১৪০ GPa, ৪০০০-৬০০০ K) বিভিন্ন লোহার সামগ্রী (XFe = Fe/(Fe+Mg) = 0.12, 0.5, 1.0) সহ সিলিকেট তরলের বৈদ্যুতিক পরিবাহিতা σel এবং তাপীয় পরিবাহিতা kel গণনা করা।

গণনা পদ্ধতির স্থাপত্য

1. আণবিক গতিশীলতা সিমুলেশন

  • তাত্ত্বিক ভিত্তি: ঘনত্ব কার্যকরী তত্ত্ব (DFT) + PBEsol আনুমানিকতা + Hubbard U সংশোধন (U-J = 2.5 eV)
  • বাস্তবায়ন প্ল্যাটফর্ম: VASP কোড, প্রজেক্টেড অগমেন্টেড প্লেন ওয়েভ পদ্ধতি
  • সিমুলেশন প্যারামিটার:
    • সিস্টেম: ১৪৯ পরমাণু, ৬টি উপাদান (Mg, Fe, Si, O, Ca, Al)
    • সময় ধাপ: ১ fs, মোট সময়কাল: ১০-১৫ ps
    • শক্তি কাটঅফ: ৫০০ eV
    • K পয়েন্ট নমুনা: Gamma পয়েন্ট
    • সমষ্টি: NVT, Nosé-Hoover তাপমাত্রা নিয়ন্ত্রণ

2. ইলেকট্রন পরিবহন বৈশিষ্ট্য গণনা

Chester-Thellung ফর্মের Kubo-Greenwood পদ্ধতি ব্যবহার করা হয়েছে:

বৈদ্যুতিক পরিবাহিতা: σel(ω)=L11(ω)\sigma_{el}(\omega) = L_{11}(\omega)

তাপীয় পরিবাহিতা: kel(ω)=1e2T(L22(ω)−L12(ω)L21(ω)L11(ω))k_{el}(\omega) = \frac{1}{e^2T}\left(L_{22}(\omega) - \frac{L_{12}(\omega)L_{21}(\omega)}{L_{11}(\omega)}\right)

যেখানে গতিশীল সহগ: Lαβ(ω)=2πe2ℏ23me2ωΩ∑i,j(fi−fj)δ(ϵi−ϵj−ℏω)∣⟨ψi∣∇∣ψj⟩∣2(−1)α+β(ϵi−μ)α−1(ϵj−μ)β−1L_{\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. স্পিন অবস্থা পরিচালনা

লোহার উচ্চ স্পিন-নিম্ন স্পিন রূপান্তর বিবেচনা করে, তাপগতিবিদ্যা একীকরণের মাধ্যমে মুক্ত শক্তির পার্থক্য গণনা করা হয়েছে: ΔFHS−LS=∫01⟨ΔU⟩λdλ\Delta F_{HS-LS} = \int_0^1 \langle\Delta U\rangle_\lambda d\lambda

ভারসাম্য অবস্থার বৈদ্যুতিক পরিবাহিতা: 1σeleq=feqσelHS+1−feqσelLS\frac{1}{\sigma_{el}^{eq}} = \frac{f_{eq}}{\sigma_{el}^{HS}} + \frac{1-f_{eq}}{\sigma_{el}^{LS}}

প্রযুক্তিগত উদ্ভাবন পয়েন্ট

  1. বহু-স্কেল সংযোগ: পরমাণু স্কেল DFT গণনা এবং ম্যাক্রোস্কোপিক তাপীয় বিবর্তন মডেল একত্রিত করা
  2. স্পিন অবস্থা পরিচালনা: উচ্চ চাপে লোহার স্পিন রূপান্তর প্রভাব যথাযথভাবে বিবেচনা করা
  3. সংমিশ্রণ বিবর্তন: উপাদান বিতরণ সহগের উপর ভিত্তি করে BMO সংমিশ্রণের সময়ের বিবর্তন
  4. অরৈখিক পর্যায় চিত্র: অরৈখিক তরল লাইন তাপমাত্রা মডেল ব্যবহার করে, Fe-Mg সিস্টেম আরও বাস্তবসম্মতভাবে প্রতিফলিত করা

পরীক্ষামূলক সেটআপ

গণনা প্যারামিটার

  • চাপের পরিসীমা: ১০০-১৪০ GPa (BMO সম্পর্কিত চাপ)
  • তাপমাত্রার পরিসীমা: ৪০০০-৬০০০ K (BMO সম্পর্কিত তাপমাত্রা)
  • লোহার সামগ্রী: XFe = 0.12 (পেরিডোটাইট), 0.5, 1.0 (লোহা-সমৃদ্ধ শেষ সদস্য)
  • সংযোগ মানদণ্ড: শক্তি ৩ meV/atom, চাপ ০.২ GPa

মূল্যায়ন সূচক

  • বৈদ্যুতিক পরিবাহিতা: σel (S/m), ডায়নামো প্রান্তিক মূল্যের সাথে তুলনা (~10⁴ S/m)
  • তাপীয় পরিবাহিতা: kel (W/m/K), প্রবাহ স্থিতিশীলতা মূল্যায়ন
  • চৌম্বক রেনল্ডস সংখ্যা: Rm = μ₀vlσ, ডায়নামো মানদণ্ড (Rm > 40)

তুলনা মানদণ্ড

  • বিশুদ্ধ MgSiO₃ সিস্টেমের আয়নিক তাপীয় পরিবাহিতা (~৪-৫ W/m/K)
  • পৃথিবীর মূলের তরল লোহার বৈদ্যুতিক এবং তাপীয় পরিবাহিতা
  • Wiedemann-Franz আইন পূর্বাভাস মূল্য

পরীক্ষামূলক ফলাফল

প্রধান ফলাফল

1. লোহার সামগ্রীর সাথে বৈদ্যুতিক পরিবাহিতার পরিবর্তন

  • XFe = 0.12: σel ~ 2×10⁴ S/m (6000K, 100GPa)
  • XFe = 1.0: σel > 10⁵ S/m, প্রায় ৫ গুণ বৃদ্ধি
  • তাপমাত্রা নির্ভরতা: σ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 অবস্থার সংকর প্রশস্ত শক্তি ব্যান্ড গঠন করে
  • Fermi স্তরে অবস্থার ঘনত্ব g(EF) XFe এর সাথে রৈখিকভাবে বৃদ্ধি পায়
  • σel ∝ g(EF)², Mott-Ziman তত্ত্বের সাথে সামঞ্জস্যপূর্ণ

চৌম্বক ক্ষেত্র বিবর্তন সিমুলেশন

সময়ের বিবর্তন ফলাফল

  • ধ্রুবক সংমিশ্রণ (XFe=0.12): ডায়নামো ১.৪ Gyr স্থায়ী হয়
  • লোহা সমৃদ্ধকরণ বিবেচনা: ডায়নামো ৩.৩ Gyr স্থায়ী হয় (১.২ Ga আগে পর্যন্ত)
  • সমালোচনামূলক পুরুত্ব: ডায়নামো বজায় রাখার ন্যূনতম BMO পুরুত্ব ২০০ km থেকে ১০০ km এ হ্রাস পায়

ডায়নামো অন্তর্বর্তী বিশ্লেষণ

Regime diagram এর মাধ্যমে নির্ধারণ করা হয়েছে:

  • BMO সর্বদা তাপীয় প্রবাহ ডায়নামো অন্তর্বর্তীতে অবস্থিত
  • তাপীয় পরিবাহিতা অত্যধিক হয়ে ডায়নামো দমন করার অন্তর্বর্তীতে পৌঁছায় না
  • বাহ্যিক মূলের তুলনায়, BMO আরও অনুকূল σ-k প্যারামিটার সমন্বয় রয়েছে

বিলোপন পরীক্ষা

  • স্পিন অবস্থা প্রভাব: নিম্ন স্পিন অবস্থার বৈদ্যুতিক পরিবাহিতা উচ্চ স্পিন অবস্থার প্রায় ২ গুণ
  • বেগ স্কেলিং: CIA স্কেলিং MLT স্কেলিংয়ের চেয়ে ১.৩ Gyr কম ডায়নামো আয়ু প্রদান করে
  • আয়নিক অবদান: ইলেকট্রন পরিবাহিতা মোট পরিবাহিতার ৭০-৯০% গঠন করে

সম্পর্কিত কাজ

পরিবাহিতা গণনা

  • Holmström et al. (2018): MgO, MgSiO₃ সরল সিস্টেম
  • Scipioni et al. (2017): পৃথিবীর সামগ্রিক সিলিকেট সংমিশ্রণ
  • Soubiran & Militzer (2018): উচ্চ চাপে সিলিকেট তরল

BMO বিবর্তন মডেল

  • Labrosse et al. (2007): ক্লাসিক্যাল BMO তাপীয় বিবর্তন মডেল
  • Stixrude et al. (2020): সিলিকেট ডায়নামো তাত্ত্বিক কাঠামো
  • Ziegler & Stegman (2013): BMO ডায়নামো সংখ্যাগত সিমুলেশন

লোহার স্পিন রূপান্তর

  • Holmström & Stixrude (2016): তরল অক্সাইডে স্পিন রূপান্তর
  • Shim et al. (2023): সংঘর্ষ সংকোচন পরীক্ষা যাচাইকরণ

উপসংহার এবং আলোচনা

প্রধান উপসংহার

  1. লোহা সমৃদ্ধকরণ উল্লেখযোগ্যভাবে বৈদ্যুতিক পরিবাহিতা বৃদ্ধি করে: XFe ০.১২ থেকে ১.০ এ বৃদ্ধি পেলে, σel প্রায় এক দশক বৃদ্ধি পায়
  2. ডায়নামো আয়ু বৃদ্ধি: লোহা সমৃদ্ধকরণ বিবেচনা করার পরে, সিলিকেট ডায়নামো ৩.৩ বিলিয়ন বছর বজায় রাখতে পারে
  3. প্রবাহ এখনও বজায় রাখা যায়: তাপীয় পরিবাহিতা বৃদ্ধি সত্ত্বেও, BMO প্রবাহ অস্থিরতা বজায় রাখে
  4. প্রাথমিক চৌম্বক ক্ষেত্র সম্ভব: পৃথিবীর প্রাথমিক চৌম্বক ক্ষেত্রের জন্য যুক্তিসঙ্গত ভৌত প্রক্রিয়া প্রদান করে

সীমাবদ্ধতা

  1. সংমিশ্রণ সরলীকরণ: স্ফটিকীকরণের সাথে SiO₂ সামগ্রীর হ্রাস বিবেচনা করা হয়নি
  2. প্রবাহ প্যাটার্ন: সংমিশ্রণ প্রবাহের অবদান অন্তর্ভুক্ত করা হয়নি
  3. বিকিরণ তাপ স্থানান্তর: বিকিরণ তাপীয় পরিবাহিতা অবদান উপেক্ষা করা হয়েছে
  4. ত্রিমাত্রিক প্রভাব: এক-মাত্রিক তাপীয় বিবর্তন মডেলের উপর ভিত্তি করে

ভবিষ্যত দিকনির্দেশনা

  1. অতি-পৃথিবী প্রয়োগ: উচ্চতর চাপে বৈদ্যুতিক পরিবাহিতা বৈশিষ্ট্য গবেষণা করা
  2. সংমিশ্রণ প্রবাহ: লোহা অবসরণ দ্বারা চালিত প্রবাহ প্রভাব পরিমাণ করা
  3. বহু-উপাদান সিস্টেম: Al, Ca ইত্যাদি উপাদানের প্রভাব বিবেচনা করা
  4. পরীক্ষামূলক যাচাইকরণ: উচ্চ তাপমাত্রা উচ্চ চাপ পরীক্ষা পরিমাপ যাচাইকরণ

গভীর মূল্যায়ন

সুবিধা

  1. উন্নত পদ্ধতি: প্রথম নীতি গণনা এবং ম্যাক্রোস্কোপিক মডেল একত্রিত করে, বহু-স্কেল গবেষণা
  2. ভৌত স্পষ্টতা: ইলেকট্রন কাঠামো দৃষ্টিকোণ থেকে পরিবাহন প্রক্রিয়া ব্যাখ্যা করা
  3. ফলাফল গুরুত্ব: পৃথিবীর প্রাথমিক চৌম্বক ক্ষেত্রের উৎসের জন্য পরিমাণগত সীমাবদ্ধতা প্রদান করা
  4. মডেল সম্পূর্ণতা: স্পিন রূপান্তর, সংমিশ্রণ বিবর্তন ইত্যাদি একাধিক ভৌত প্রক্রিয়া বিবেচনা করা

অপূর্ণতা

  1. গণনা খরচ: DFT+U পদ্ধতির প্যারামিটার নির্বাচন ফলাফল প্রভাবিত করতে পারে
  2. মডেল অনুমান: এক-মাত্রিক তাপীয় বিবর্তন মডেল অত্যন্ত সরলীকৃত
  3. পরীক্ষামূলক অভাব: উচ্চ তাপমাত্রা উচ্চ চাপ পরীক্ষা ডেটা যাচাইকরণের অভাব
  4. অনিশ্চয়তা: কিছু প্যারামিটার (যেমন বিতরণ সহগ) উল্লেখযোগ্য অনিশ্চয়তা রয়েছে

প্রভাব

  1. একাডেমিক মূল্য: পৃথিবী পদার্থবিজ্ঞান এবং গ্রহ বিজ্ঞানের জন্য গুরুত্বপূর্ণ তাত্ত্বিক ভিত্তি প্রদান করা
  2. প্রয়োগ সম্ভাবনা: বহিঃগ্রহ চৌম্বক ক্ষেত্র গবেষণায় সম্প্রসারণ করা যায়
  3. পদ্ধতি প্রদর্শন: গণনামূলক পৃথিবী পদার্থবিজ্ঞানের গবেষণা প্যারাডাইম প্রদর্শন করা
  4. নীতি তাৎপর্য: পৃথিবীর বাসযোগ্যতা বিবর্তনের বোঝাপড়া বৃদ্ধি করা

প্রযোজ্য পরিস্থিতি

  • পৃথিবীর প্রাথমিক চৌম্বক ক্ষেত্র বিবর্তন গবেষণা
  • শিলাময় বহিঃগ্রহ চৌম্বক ক্ষেত্র পূর্বাভাস
  • উচ্চ তাপমাত্রা উচ্চ চাপে সিলিকেট বৈশিষ্ট্য গবেষণা
  • পৃথিবীর অভ্যন্তরীণ কাঠামো এবং গতিশীলতা মডেলিং

তথ্যসূত্র

এই পত্রটি পৃথিবী পদার্থবিজ্ঞান, গণনামূলক পদার্থবিজ্ঞান এবং গ্রহ বিজ্ঞান ক্ষেত্রের গুরুত্বপূর্ণ সাহিত্য উদ্ধৃত করে, যার মধ্যে রয়েছে:

  • Labrosse et al. (2007): BMO তাপীয় বিবর্তনের ক্লাসিক্যাল মডেল
  • Christensen & Aubert (2006): ডায়নামো চৌম্বক রেনল্ডস সংখ্যা প্রান্তিক মূল্য
  • Pozzo et al. (2012): পৃথিবীর মূল তাপীয় পরিবাহিতা গণনা
  • McDonough & Sun (1995): পৃথিবীর পেরিডোটাইট সংমিশ্রণ মডেল

এই পত্রটি কঠোর প্রথম নীতি গণনা এবং তাপীয় বিবর্তন মডেলিংয়ের মাধ্যমে, পৃথিবীর প্রাথমিক চৌম্বক ক্ষেত্রের উৎস এই প্রধান বৈজ্ঞানিক সমস্যার জন্য পরিমাণগত ভৌত সীমাবদ্ধতা প্রদান করে, এটি গণনামূলক পৃথিবী পদার্থবিজ্ঞান ক্ষেত্রের একটি গুরুত্বপূর্ণ অবদান।