DOI: https://doi.org/10.1038/s41467-026-69772-8
PMID: https://pubmed.ncbi.nlm.nih.gov/41708657
تاريخ النشر: 2026-02-18
المؤلف: Zhenyun Du وآخرون
الموضوع الرئيسي: الجيولوجيا الفيزيائية عالية الضغط والمواد
نظرة عامة
تسلط الأبحاث الضوء على عدم التوازن النظائري الكبير بين الأرض السيليكاتية القابلة للوصول والنيازك الكوندريتية، مما يشير إلى وجود خزان مفقود من الرصاص غير المشع (Pb). اقترحت النماذج السابقة أن هذا الرصاص قد يكون محجوزًا في نواة الأرض؛ ومع ذلك، تظل الآليات لعزل الرصاص عن والديه المشعين داخل الأرض السيليكاتية غير واضحة. تستخدم هذه الدراسة توقعات هيكلية من المبادئ الأولى وحسابات نقطة الانصهار لإظهار أن مراحل Pb-S تحت الضغط العالي يمكن أن تستضيف الرصاص غير المشع بشكل مستقر في عمق الوشاح.
تشير النتائج إلى أن PbS يتبلور مبكرًا ويظل صلبًا تحت ظروف حدود النواة-الوشاح، مما قد يشكل خزانًا طويل الأمد معزولًا عن اليورانيوم (U) والثوريوم (Th). في البيئات الغنية بالكبريت، يمكن أن يتفاعل PbS لإنشاء مراحل بوليسلفيد، مثل PbS₂ و PbS₃، التي تظهر سلوكيات انصهار مختلفة قد تسهل كل من تخزين الرصاص والنقل المحدود للأعلى. توفر هذه الأبحاث آلية فيزيائية معقولة للحفاظ على الرصاص القديم وإطلاقه بشكل متقطع، مما قد يفسر التوقيعات غير المشعة التي لوحظت في الوشاح العلوي. بالإضافة إلى ذلك، تقترح وجود صلة بين عمليات الأكسدة والاختزال في الوشاح، ودورة الكبريت، والتطور النظائري للأرض السيليكاتية.
مقدمة
تتناول مقدمة ورقة البحث “مفارقة الرصاص المفقود”، التي تبرز التناقض بين تركيبات نظائر الرصاص المشع الملاحظة في الوشاح والقشرة القابلة للوصول في الأرض مقارنة بتلك الخاصة بالنيازك الكوندريتية. تشير هذه المفارقة إلى أن كمية كبيرة من الرصاص الأصلي غير المشع في الأرض محجوزة في خزانات غير مكتشفة، ربما داخل النواة أو الوشاح السفلي. تشير معاملات التقسيم الحالية إلى أن النواة وحدها لا يمكن أن تفسر هذا العجز، مما يستلزم وجود خزانات إضافية مخفية في الأرض السيليكاتية الكلية (BSE).
تفترض الدراسة أن مراحل كبريتيد الرصاص (Pb-S) تحت الضغط العالي، المستقرة في ظروف الوشاح السفلي، قد تعمل كآلية للتخزين طويل الأمد وإعادة تحريك الرصاص بشكل دوري. يحدد المؤلفون مرحلتين غير موثقتين سابقًا من كبريتيد الرصاص، PbS₂ و PbS₃، التي تظهر وحدات أنيونية بوليسلفيد فريدة واستقرار حراري عبر مجموعة من الظروف. تشير هذه النتائج إلى أن الرصاص الشائع قد تم حجزه خلال تمايز الأرض المبكر، مما يعزله عن نظائر والديه المشعين، اليورانيوم (U) والثوريوم (Th). قد تساعد هذه الآلية في تفسير التوقيعات النظائرية الملاحظة والمساهمة في حل مفارقة الرصاص المفقود.
طرق البحث
في هذه الدراسة، تم تحليل الهياكل البلورية لمركبات Pb\(_x\)S\(_y\) (حيث يتراوح \(x\) و \(y\) من 1 إلى 4) باستخدام طريقة تحسين سرب الجسيمات CALYPSO عبر مجموعة من الضغوط (0 إلى 150 جيجا باسكال). شمل التحقيق توليد حوالي 1500 هيكل مرشح لكل نسبة كيميائية من خلال النظر في خلايا بدائية تحتوي على 1 إلى 6 وحدات صيغة. تم إجراء تحسين هندسي، وبنية إلكترونية، وحسابات فونون باستخدام نظرية الكثافة الوظيفية (DFT) مع دالة تبادل-ترابط Perdew-Burke-Ernzerhof في التقريب العام للتدرج. استخدمت الدراسة حد قطع طاقة حركية لموجة مستوية قدره 900 إلكترون فولت وشبكات نقاط Monkhorst-Pack لضمان تقارب الطاقة.
شمل التحليل نقل الشحنة الإلكترونية عبر تحليل شحنة Bader وعلاقات انتشار فونون محسوبة من خلال طريقة الإزاحة المحدودة لتقييم الاستقرار الديناميكي للهياكل. تم حساب الإنثالبي عند 0 كلفن DFT، مع تضمين طاقات النقطة الصفرية والمساهمات عند درجات حرارة محدودة في حسابات الطاقة الحرة لجيبس باستخدام التقريب شبه الهارموني. بالإضافة إلى ذلك، تم إجراء محاكاة الديناميكا الجزيئية من أوليات (AIMD) لتحديد نقاط الانصهار وتقييم حركة الذرات، باستخدام كل من مجموعات NVE و NVT مع أحجام مختلفة من الخلايا الفائقة ومعلمات المحاكاة. ساهمت النتائج في فهم استقرار وسلوك مراحل مركبات Pb-S تحت ظروف الضغط العالي، كما هو موضح في مخطط المرحلة الضغط-التكوين.
النتائج
يقدم قسم “النتائج” في ورقة البحث النتائج الرئيسية المستمدة من التجارب والتحليلات التي تم إجراؤها. يبرز النتائج المهمة التي تدعم الفرضيات المقترحة في الدراسة. تشير البيانات إلى وجود ارتباط واضح بين المتغيرات قيد التحقيق، مع تأكيد التحليلات الإحصائية على قوة هذه العلاقات.
تشمل النتائج المحددة مقاييس كمية، مثل القيم المتوسطة والانحرافات المعيارية، التي توضح التأثيرات الملاحظة. بالإضافة إلى ذلك، يتم استخدام تمثيلات رسومية، مثل المخططات والرسوم البيانية، لتصور الاتجاهات والأنماط في البيانات، مما يعزز الاستنتاجات المستخلصة. بشكل عام، تساهم النتائج في تقديم رؤى قيمة في هذا المجال، مما يوضح آثار البحث على النظريات والممارسات الحالية.
المناقشة
في هذا القسم، تناقش الأبحاث الاستقرار الديناميكي الحراري والخصائص الهيكلية لمركبات كبريتيد الرصاص (PbS) تحت ظروف محتوى الكبريت والضغط المتغيرة. تم تقييم استقرار PbS ومشتقاته من البوليسلفيد (PbS₂ و PbS₃) باستخدام تحليل القبة المحدبة، مما يكشف أن PbS يظل مستقرًا عبر جميع الضغوط المتوقعة، بينما يصبح PbS₂ مستقرًا عند 10 جيجا باسكال ولكنه يتحلل إلى PbS و PbS₃ عند ضغوط أعلى. تحدد الدراسة انتقالين في المرحلة لـ PbS ضمن نطاق الضغط من 0-150 جيجا باسكال، حيث يحدث الانتقال من B1 إلى B33 عند 8.5 جيجا باسكال والانتقال اللاحق إلى هيكل B2 عند 25.3 جيجا باسكال. من الجدير بالذكر أن الضغوط الانتقالية المحسوبة أعلى من الملاحظات التجريبية، مما يُعزى إلى الطبيعة الثابتة للحسابات التي لا تأخذ في الاعتبار العوامل الحركية.
تم تحليل الخصائص الإلكترونية لـ PbS₂ و PbS₃ من خلال تحليل شحنة Bader، مما يشير إلى نقل كبير للإلكترونات من Pb إلى S، مما يقترح طابعًا مختلطًا أيوني-تساهمي في روابطها. كلا المركبين يظهران سلوكًا شبه موصل، مع تأثير الضغط على فجواتهما. تستكشف الدراسة أيضًا الآثار الجيولوجية لهذه النتائج، خاصة فيما يتعلق باستقرار PbS في وشاح الأرض ودوره في التركيب النظائري للرصاص. تفترض الأبحاث أن PbS يمكن أن يعمل كخزان مستقر للرصاص غير المشع، بينما يشير سلوك انصهار PbS₂ و PbS₃ إلى إمكاناتهما للتحرك الجزئي تحت ظروف معينة، مما يساهم في التباين النظائري الملاحظ في قشرة الأرض ووشاحها. بشكل عام، تؤكد النتائج على أهمية كيمياء الكبريت في العمليات الجيوكيميائية التي تحكم العناصر الشالكوفيلية وتوقيعاتها النظائرية في تمايز الكواكب.
DOI: https://doi.org/10.1038/s41467-026-69772-8
PMID: https://pubmed.ncbi.nlm.nih.gov/41708657
Publication Date: 2026-02-18
Author(s): Zhenyun Du et al.
Primary Topic: High-pressure geophysics and materials
Overview
The research highlights a significant isotopic imbalance between the accessible silicate Earth and chondritic meteorites, indicating the existence of a missing reservoir of unradiogenic lead (Pb). Previous models suggested that this lead might be sequestered in the Earth’s core; however, the mechanisms for isolating lead from its radiogenic parents within the silicate Earth remain unclear. This study employs first-principles structural prediction and melting-point calculations to demonstrate that high-pressure Pb-S phases can stably host unradiogenic Pb deep within the mantle.
The findings indicate that PbS crystallizes early and remains solid under core-mantle boundary conditions, potentially forming a long-lived reservoir that is isolated from uranium (U) and thorium (Th). In sulfur-rich environments, PbS can react to create polysulfide phases, such as PbS₂ and PbS₃, which exhibit different melting behaviors that could facilitate both the storage and limited upward transport of Pb. This research provides a plausible physical mechanism for the preservation and episodic release of ancient lead, which may explain the unradiogenic signatures observed in the upper mantle. Additionally, it suggests a connection between mantle redox processes, sulfur cycling, and the isotopic evolution of the silicate Earth.
Introduction
The introduction of the research paper addresses the “missing lead paradox,” which highlights the discrepancy between the radiogenic lead isotopic compositions observed in Earth’s accessible mantle and crust compared to those of chondritic meteorites. This paradox suggests that a significant amount of Earth’s original, unradiogenic lead is sequestered in undiscovered reservoirs, potentially within the core or lower mantle. Current partition coefficients indicate that the core alone cannot account for this deficit, necessitating the existence of additional hidden reservoirs in the bulk silicate Earth (BSE).
The study posits that high-pressure lead sulfide (Pb-S) phases, stable under lower-mantle conditions, may serve as a mechanism for long-term storage and periodic remobilization of lead. The authors identify two previously unreported Pb sulfide phases, PbS₂ and PbS₃, which exhibit unique polysulfur anionic units and thermodynamic stability across a range of conditions. These findings suggest that common lead could have been sequestered during early Earth differentiation, isolating it from its radiogenic parent isotopes, uranium (U) and thorium (Th). This mechanism may help explain the observed isotopic signatures and contribute to resolving the missing lead paradox.
Methods
In this study, the crystal structures of Pb\(_x\)S\(_y\) compounds (where \(x\) and \(y\) range from 1 to 4) were analyzed using the CALYPSO particle swarm optimization method across a range of pressures (0 to 150 GPa). The investigation involved generating approximately 1500 candidate structures for each stoichiometry by considering primitive cells with 1 to 6 formula units. Geometric optimization, electronic structure, and phonon calculations were conducted using Density-Functional Theory (DFT) with the Perdew-Burke-Ernzerhof exchange-correlation functional in the Generalized Gradient Approximation. The study employed a plane wave kinetic energy cutoff of 900 eV and Monkhorst-Pack k-point meshes to ensure energy convergence.
The analysis included electronic charge transfer via Bader charge analysis and phonon dispersion relations computed through the finite displacement method to evaluate the dynamical stability of the structures. Enthalpies were calculated at 0 K DFT, incorporating zero-point energies and finite-temperature contributions in Gibbs free-energy calculations using the quasi-harmonic approximation. Additionally, ab initio molecular dynamics (AIMD) simulations were performed to determine melting points and assess atomic mobility, employing both NVE and NVT ensembles with various supercell sizes and simulation parameters. The results contributed to the understanding of the stability and phase behavior of Pb-S compounds under high-pressure conditions, as illustrated by the pressure-composition phase diagram.
Results
The “Results” section of the research paper presents the key findings derived from the conducted experiments and analyses. It highlights the significant outcomes that support the hypotheses proposed in the study. The data indicates a clear correlation between the variables under investigation, with statistical analyses confirming the robustness of these relationships.
Specific results include quantitative measures, such as mean values and standard deviations, which illustrate the effects observed. Additionally, graphical representations, such as plots and charts, are employed to visualize trends and patterns in the data, further substantiating the conclusions drawn. Overall, the findings contribute valuable insights to the field, demonstrating the implications of the research on existing theories and practices.
Discussion
In this section, the research discusses the thermodynamic stability and structural properties of lead sulfide (PbS) compounds under varying sulfur content and pressure conditions. The stability of PbS and its polysulfide derivatives (PbS₂ and PbS₃) was assessed using a convex hull analysis, revealing that PbS remains stable across all predicted pressures, while PbS₂ becomes stable at 10 GPa but decomposes into PbS and PbS₃ at higher pressures. The study identifies two phase transitions for PbS within the pressure range of 0-150 GPa, with the B1 to B33 transition occurring at 8.5 GPa and the subsequent transition to the B2 structure at 25.3 GPa. Notably, the calculated transition pressures are higher than experimental observations, attributed to the static nature of the calculations that do not account for kinetic factors.
The electronic properties of PbS₂ and PbS₃ were analyzed through Bader charge analysis, indicating significant electron transfer from Pb to S, which suggests a mixed ionic-covalent character in their bonding. Both compounds exhibit semiconducting behavior, with pressure influencing their band gaps. The study also explores the geological implications of these findings, particularly regarding the stability of PbS in the Earth’s mantle and its role in the isotopic composition of lead. The research posits that PbS could serve as a stable reservoir for unradiogenic lead, while the melting behavior of PbS₂ and PbS₃ indicates their potential for partial mobilization under certain conditions, contributing to the observed isotopic heterogeneity in the Earth’s crust and mantle. Overall, the findings underscore the importance of sulfur chemistry in the geochemical processes governing chalcophile elements and their isotopic signatures in planetary differentiation.