DOI: https://doi.org/10.1016/j.epsl.2025.119540
تاريخ النشر: 2025-07-22
المؤلف: Jac van Driel وآخرون
الموضوع الرئيسي: الجيولوجيا الفيزيائية عالية الضغط والمواد
نظرة عامة
تبحث هذه الدراسة في تكوين قلب المريخ باستخدام حسابات نظرية الكثافة الوظيفية المتقدمة (DFT)، مع التركيز على الخصائص المغناطيسية لسبائك الحديد لتحديد الكثافة وسرعات الموجات بدقة. يتم دمج النتائج في معادلة حالة جديدة لنموذج الانحدار Gaussian Process (GPR)، والتي تُستخدم لاستكشاف تكوينات مختلفة ضمن نظام Fe-S-O-C-H التي تتماشى مع الكثافة المرصودة وسرعات الموجات الضاغطة عبر عمق قلب المريخ. تأخذ الدراسة في الاعتبار نماذج داخلية مختلفة للمريخ، وخاصة تداعيات طبقة الانصهار العميقة في قاعدة الوشاح، مما يغير بشكل كبير التقديرات السابقة لحجم القلب وكثافته. تكشف التحليلات أن التغيرات في درجة الحرارة (التي تتراوح من 1900 كلفن إلى 2850 كلفن) لها تأثير ضئيل على النتائج التركيبية، حيث تتطلب جميع النماذج تركيزات عالية من العناصر الخفيفة (حوالي 50 مول%).
تؤكد الاستنتاجات على أهمية نمذجة التركيب الدقيق لقلب المريخ من خلال GPR-EOS المطور، الذي يدمج المبادئ الديناميكية الحرارية لتتوافق مع البيانات التجريبية حول سرعة الموجات الضاغطة والكثافة. تؤكد الدراسة على ضرورة وجود تركيزات عالية من العناصر الخفيفة (45-55 مول%) وتقترح أن التركيبات التي تناسب النماذج الجيوفيزيائية دون قيود كيميائية جيوكيميائية تتطلب عادة محتوى كبير من الكربون (20-30 مول%). ومع ذلك، عند تطبيق قيود من تجارب التقسيم وأقصى محتوى من الكبريت، يضيق الفضاء التركيبي بشكل كبير، مما يتطلب تركيزات عالية من الهيدروجين (>27 مول%) وأكسجين ضئيل. تشير الأبحاث إلى أنه بينما تعتبر درجة حرارة حدود القلب-الوشاح عاملاً حاسمًا، فإن النطاق المختبر لا يستبعد أي نماذج جيوفيزيائية مقترحة لقلب المريخ.
مقدمة
لقد عززت مهمة InSight بشكل ملحوظ فهمنا لخصائص قلب المريخ، بما في ذلك حجمه وكثافته وسرعات الموجات الزلزالية. أشارت النتائج الأولية إلى وجود قلب معدني سائل منخفض الكثافة بقطر يبلغ حوالي $1830 \pm 40$ كم وكثافة تتراوح بين $5.7$ و $6.3 \, \text{g/cm}^3$ (Stahler et al., 2021). ومع ذلك، قامت التحليلات اللاحقة بتحسين هذه التقديرات، مشيرة إلى كثافات قلب أعلى تتراوح بين $6.2$ إلى $6.3 \, \text{g/cm}^3$ وسرعة موجة P عند سطح القلب السائل تتراوح بين $4.9$ إلى $5 \, \text{km/s}$ (Irving et al., 2023; Le Maistre et al., 2023; Khan et al., 2022).
قدمت الدراسات الحديثة مفهوم طبقة سيليكات سائلة (LSL) في قاعدة وشاح المريخ، مقترحة سمكًا يبلغ حوالي $150$ إلى $170$ كم ومراجعة حجم القلب إلى $1,650 \pm 20$ كم أو $1,675 \pm 30$ كم، مع زيادة تقديرات الكثافة إلى $6.5 \, \text{g/cm}^3$ أو $6.65 \, \text{g/cm}^3$ (Samuel et al., 2023; Khan et al., 2023). يوفر دمج البيانات الجيوفيزيائية مع الرؤى التجريبية حول تقسيم العناصر – وخاصة تفاعلات الكبريت والأكسجين والكربون داخل سبائك الحديد تحت ظروف قلب المريخ – فهمًا أكثر شمولاً لتكوين القلب. يسمح هذا النهج بتقييم التركيبات المختلفة للقلب مقابل النماذج الجيوفيزيائية، مما يعزز فهمنا لخصائص قلب المريخ وتداعياته على تشكيله وتطوره (Tsuno et al., 2018; Gendre et al., 2022; Steenstra and van Westrenen, 2018).
الطرق
في هذه الدراسة، يهدف المؤلفون إلى تحديد نطاق التركيبات السائلة المتوافقة مع الكثافة المرصودة ($\rho$) ونسبة الحجم ($V_\phi$) لقلب المريخ، كما تم استنتاجه من بيانات InSight. تم هيكلة المنهجية في ثلاث خطوات رئيسية. في البداية، يتم استخدام محاكاة نظرية الكثافة الوظيفية (DFT) لحساب الضغط والطاقة الداخلية لـ 60 تركيبة متميزة عبر نطاق ضغط من 0-100 جيجا باسكال ودرجات حرارة من 2000-6000 كلفن، مما ينتج مجموعة بيانات تحتوي على 716 نقطة بيانات فريدة.
تتضمن الخطوة الثانية تحسين تمثيل الحالة المغناطيسية لنظام Fe-S-O-C-H لتقليل الفجوات بين الحسابات النظرية والنتائج التجريبية. يتم تحقيق ذلك من خلال تقييم الفروق في الطاقة الحرة لهيلمهولتز بين حالات مغناطيسية مختلفة، مما يؤدي إلى تحديد أدنى حالة طاقة حرة. ثم يتم تطوير تصحيح ضغط يعتمد على التركيب وتطبيقه على مجموعة بيانات المحاكاة غير المغناطيسية، والتي يمكن أيضًا استخدامها في دراسات أخرى. أخيرًا، يتم بناء معادلة حالة عملية Gaussian (GPR-EOS) لاستيفاء الكثافة والطاقة الداخلية عبر الفضاء متعدد الأبعاد للضغط ودرجة الحرارة والتركيب. يسمح هذا الإطار من GPR-EOS بالاستيفاء الذاتي المتسق للخصائص الديناميكية الحرارية دون الحاجة إلى شكل وظيفي محدد مسبقًا، مما يسهل اشتقاق الأديبات المعتمدة على التركيب التي تقيد خصائص قلب المريخ بشكل أكبر.
النتائج
تؤكد نتائج هذه الدراسة على الدور الحاسم للمغناطيسية في خصائص سبائك الحديد ذات الضغط المنخفض، خاصة في السياقات الكوكبية الصغيرة. ينخفض العزم المغناطيسي لذرات الحديد في هذه السبائك لوغاريتميًا من حوالي 2 مغناطيس بوهري لكل ذرة في الظروف العادية إلى حوالي 1 مغناطيس بوهري لكل ذرة عند 100 جيجا باسكال. من الجدير بالذكر أن FeO السائل يظهر أعلى العزوم المغناطيسية بين التركيبات المدروسة، متجاوزًا FeS وFeC وFeH. تشير النتائج إلى وجود ارتباط مباشر بين درجة الحرارة والعزم المغناطيسي، مما يبرز ضرورة وجود تصحيح ضغط مغناطيسي يعتمد على التركيب بدلاً من تعديل موحد، كما تم استخدامه سابقًا في دراسات أخرى.
تظهر معادلة حالة الدراسة (GPR-EOS) توافقًا قويًا مع البيانات التجريبية للحديد النقي ومختلف سبائك Fe-S وFe-C، خاصة عند الضغوط الأعلى ذات الصلة بقلب المريخ. بينما توجد بعض الفجوات بين مجموعات البيانات التجريبية، فإن الاتفاق العام مرضٍ، خاصة بالنسبة لسرعات الموجات الضاغطة والكثافات. كما تتناقض الأبحاث مع الدراسات السابقة من قبل Huang et al. (2023)، خاصة في معالجة تصحيحات الضغط المغناطيسي وتعديلات الطاقة الداخلية، والتي تعتبر حاسمة لتحديد الخصائص الديناميكية الحرارية بدقة واستنتاج تركيب قلب المريخ. تشير النتائج إلى أن العناصر الخفيفة، وخاصة الكربون، قد تؤثر بشكل كبير على كثافة وسرعة الموجات الضاغطة لسبائك FeX تحت ظروف المريخ، مع تداعيات لفهم تركيب القلب.
المناقشة
في هذا القسم، يوضح المؤلفون منهجيتهم الحاسوبية للتحقيق في تركيب قلب المريخ باستخدام حزمة المحاكاة الأولية فيينا (VASP) لكل من الحسابات غير المغناطيسية والمغناطيسية. استخدمت الحسابات غير المغناطيسية دالة تبادل-ترابط PBE وشملت مجموعة متنوعة من إمكانيات الموجات المعززة (PAW)، مع التركيز على تحقيق تقارب الطاقة ضمن 6.4 ميلي إلكترون فولت/ذرة. استخدم المؤلفون منظم حرارة Nose-Hoover لمحاكاة الديناميكا الجزيئية، مما يضمن التوازن المناسب والتحليل على مدى خطوات زمنية متغيرة. بالنسبة للحسابات المغناطيسية، تناولوا أهمية أخذ إنترتروبي المغناطيسية في الاعتبار، خاصة في سبائك الحديد السائلة عند الضغوط المنخفضة، مستخدمين التكامل الديناميكي الحراري لتحديد العزم المغناطيسي الصحيح وتطبيق تصحيحات على الضغط والطاقة الداخلية.
طور المؤلفون أيضًا معادلة حالة انحدار عملية Gaussian (GPR) لاستيفاء الضغط والطاقة الداخلية عبر كثافات ودرجات حرارة وتركيبات مختلفة. تم التحقق من صحة هذا النموذج ضد مجموعات بيانات غير مرئية، مما يظهر قدرته التنبؤية لمختلف تركيبات قلب المريخ. أشارت النتائج إلى أن التركيزات العالية من العناصر الخفيفة (45-55 مول%) ضرورية لمطابقة البيانات الزلزالية، مع كون الكربون مكونًا حاسمًا. كما تسلط الدراسة الضوء على تأثير القيود الجيوكيميائية على الفضاء التركيبي، مما يبرز الحاجة إلى تركيزات أعلى من الهيدروجين عندما يكون محتوى الكبريت محدودًا. بشكل عام، تسهم النتائج في فهم دقيق لخصائص قلب المريخ الحرارية والفيزيائية، مما يشير إلى أنه بينما تناسب تركيبات متعددة النماذج الجيوفيزيائية الحالية، فإن المزيد من التحقيقات مطلوبة لتحسين هذه التنبؤات.
DOI: https://doi.org/10.1016/j.epsl.2025.119540
Publication Date: 2025-07-22
Author(s): Jac van Driel et al.
Primary Topic: High-pressure geophysics and materials
Overview
This research investigates the composition of the Martian core using advanced density functional theory (DFT) calculations, focusing on the magnetic properties of iron alloys to accurately determine density and wave velocities. The findings are integrated into a novel Gaussian Process Regression (GPR) equation of state (EOS), which is employed to explore various compositions within the Fe-S-O-C-H system that align with the observed density and compressional wave speeds throughout the Martian core’s depth. The study considers different interior models of Mars, particularly the implications of a deep melt layer at the base of the mantle, which significantly alters previous estimates of core size and density. The analysis reveals that temperature variations (ranging from 1900 K to 2850 K) have a minimal impact on compositional outcomes, with all models necessitating high concentrations of light elements (approximately 50 mol%).
The conclusions emphasize the importance of accurately modeling the Martian core’s composition through the developed GPR-EOS, which incorporates thermodynamic principles to align with experimental data on compressional wave speed and density. The study confirms the necessity of high light element concentrations (45-55 mol%) and suggests that compositions fitting geophysical models without geochemical constraints typically require substantial carbon content (20-30 mol%). However, when applying constraints from partitioning experiments and maximum sulfur content, the compositional space narrows significantly, necessitating high hydrogen concentrations (>27 mol%) and minimal oxygen. The research indicates that while core-mantle boundary temperature is a crucial factor, the tested range does not exclude any proposed geophysical models for the Martian core.
Introduction
The InSight mission has markedly enhanced our comprehension of the Martian core’s characteristics, including its size, density, and seismic wave velocities. Initial findings indicated a low-density liquid metal core with a radius of approximately $1830 \pm 40$ km and a density between $5.7$ and $6.3 \, \text{g/cm}^3$ (Stahler et al., 2021). However, subsequent analyses refined these estimates, suggesting higher core densities of $6.2$ to $6.3 \, \text{g/cm}^3$ and a P-wave velocity at the liquid core’s surface ranging from $4.9$ to $5 \, \text{km/s}$ (Irving et al., 2023; Le Maistre et al., 2023; Khan et al., 2022).
Recent studies have introduced the concept of a liquid silicate layer (LSL) at the base of the Martian mantle, proposing a thickness of approximately $150$ to $170$ km and revising the core’s size to $1,650 \pm 20$ km or $1,675 \pm 30$ km, with density estimates increasing to $6.5 \, \text{g/cm}^3$ or $6.65 \, \text{g/cm}^3$ (Samuel et al., 2023; Khan et al., 2023). The integration of geophysical data with experimental insights into elemental partitioning—particularly the interactions of sulfur, oxygen, and carbon within iron alloys under Martian core conditions—provides a more comprehensive understanding of the core’s composition. This approach allows for the evaluation of various core compositions against geophysical models, thereby refining our understanding of the Martian core’s properties and the implications for its formation and evolution (Tsuno et al., 2018; Gendre et al., 2022; Steenstra and van Westrenen, 2018).
Methods
In this study, the authors aim to determine the range of liquid compositions consistent with the observed density ($\rho$) and volume fraction ($V_\phi$) of the Martian core, as inferred from InSight data. The methodology is structured into three primary steps. Initially, density functional theory (DFT) simulations are employed to calculate the pressure and internal energy for 60 distinct compositions across a pressure range of 0-100 GPa and temperatures from 2000-6000 K, yielding a dataset of 716 unique data points.
The second step involves refining the magnetic state representation of the Fe-S-O-C-H system to minimize discrepancies between theoretical calculations and experimental results. This is achieved by assessing the Helmholtz free energy differences between various magnetic states, leading to the identification of the lowest free energy state. A compositional-dependent pressure correction is then developed and applied to the dataset of non-magnetic simulations, which can also be utilized in other studies. Lastly, a Gaussian process equation of state (GPR-EOS) is constructed to interpolate density and internal energy across the multidimensional space of pressure, temperature, and composition. This GPR-EOS framework allows for the self-consistent interpolation of thermodynamic properties without requiring a predefined functional form, facilitating the derivation of compositionally dependent adiabats that further constrain the properties of the Martian core.
Results
The results of this study emphasize the critical role of magnetism in the properties of low-pressure iron alloys, particularly in small planetary contexts. The magnetic moment of iron atoms in these alloys decreases logarithmically from approximately 2 Bohr Magnetons per atom at ambient conditions to about 1 Bohr Magneton per atom at 100 GPa. Notably, liquid FeO exhibits the highest magnetic moments among the studied compositions, surpassing FeS, FeC, and FeH. The findings indicate a direct correlation between temperature and magnetic moment, highlighting the necessity for a compositionally dependent magnetic pressure correction rather than a uniform adjustment, as previously employed in other studies.
The study’s equation of state (GPR-EOS) demonstrates a strong fit with experimental data for pure iron and various Fe-S and Fe-C alloys, particularly at higher pressures relevant to the Martian core. While some discrepancies exist between experimental datasets, the overall agreement is satisfactory, especially for compressional wave speeds and densities. The research also contrasts with previous studies by Huang et al. (2023), particularly in the treatment of magnetic pressure corrections and internal energy adjustments, which are crucial for accurately determining thermodynamic properties and inferring the composition of the Martian core. The results suggest that light elements, particularly carbon, may significantly influence the density and compressional wave speed of FeX alloys under Martian conditions, with implications for understanding the core’s composition.
Discussion
In this section, the authors detail their computational methodology for investigating the composition of the Martian core using the Vienna Ab initio Simulation Package (VASP) for both non-magnetic and magnetic calculations. Non-magnetic calculations employed the PBE exchange-correlation function and included various projector augmented wave (PAW) pseudopotentials, with a focus on achieving energy convergence within 6.4 meV/atom. The authors utilized a Nose-Hoover thermostat for molecular dynamics simulations, ensuring proper equilibration and analysis over varying time steps. For magnetic calculations, they addressed the significance of accounting for magnetic entropy, particularly in liquid iron alloys at low pressures, employing thermodynamic integration to determine the correct magnetic moment and applying corrections to pressure and internal energy.
The authors further developed a Gaussian process regression (GPR) equation of state to interpolate pressure and internal energy across different densities, temperatures, and compositions. This model was validated against unseen datasets, demonstrating its predictive capability for various Martian core compositions. The results indicated that high concentrations of light elements (45-55 mol%) are necessary to match seismic data, with carbon being a critical component. The study also highlights the impact of geochemical constraints on the compositional space, emphasizing the need for higher hydrogen concentrations when sulfur content is limited. Overall, the findings contribute to a nuanced understanding of the Martian core’s thermal and physical properties, suggesting that while multiple compositions fit existing geophysical models, further investigation is required to refine these predictions.