فحص حسابي عالي الإنتاجية وتعلم آلي قابل للتفسير لالتقاط اليود من الهياكل العضوية المعدنية
High throughput computational screening and interpretable machine learning for iodine capture of metal-organic frameworks

المجلة: npj Computational Materials، المجلد: 11، العدد: 1
DOI: https://doi.org/10.1038/s41524-025-01617-2
تاريخ النشر: 2025-05-02
المؤلف: Haoyi Tan وآخرون
الموضوع الرئيسي: الأطر العضوية المعدنية: التركيب والتطبيقات

طرق

في هذه الدراسة، تم إجراء محاكاة مونت كارلو العامة (GCMC) باستخدام برنامج RASPA لتحليل سلوك امتصاص اليود (I₂) في الهياكل العضوية المعدنية (MOFs) تحت ظروف ذات صلة بالصناعة النووية، تحديدًا عند درجة حرارة 423 كلفن وضغط 1 بار. كانت بيئة المحاكاة تحاكي ظروف الرطوبة العالية التي تواجه أثناء مرحلة المعالجة اللاحقة للوقود النووي المستهلك، والتي تتكون من نظام غاز مختلط يحتوي على 300 جزء في المليون من I₂، و68.5% من N₂، و18.4% من O₂، و12.2% من H₂O، مما يحقق رطوبة نسبية تبلغ 100%. تم اعتبار MOFs كهياكل ثابتة صلبة مع شروط حدود دورية، وتم استخدام خلايا فائقة لضمان أن أبعاد النظام تتجاوز ضعف مسافة القطع البالغة 12 Å.

تم حساب انتقائية I₂ أثناء الامتصاص باستخدام الصيغة:

\[
\text{selectivity}_{I_2} = \frac{X_{I_2}/Y_{I_2}}{X_{\text{others}}/Y_{\text{others}}}
\]

حيث تمثل \(X_{I_2}\) و\(Y_{I_2}\) كميات الامتصاص وتركيز الغاز لI₂، بينما تتوافق \(X_{\text{others}}\) و\(Y_{\text{others}}\) مع تلك لمكونات الغاز الأخرى (N₂، O₂، وH₂O). شملت كل محاكاة GCMC مرحلة توازن من 50,000 دورة تلتها مرحلة إنتاج من 50,000 دورة، تضمنت حركات جزيئية متنوعة. تم نمذجة اليود ككيان كروي مع معلمات فان دير فالز المستمدة من لزوجته، بينما تم تمثيل الماء باستخدام نموذج TIP3P. تم نمذجة N₂ وO₂ كتمثيلات ثلاثية المواقع، وتم استخدام حقل القوة العالمي (UFF) لمعلمات لينارد-جونز لهياكل MOF. تم وصف التفاعلات بين الذرات باستخدام كل من دوال الطاقة المحتملة لينارد-جونز والكهربائية.

نتائج

يقدم قسم “النتائج” في ورقة البحث النتائج الرئيسية المستمدة من التجارب أو التحليلات التي تم إجراؤها. يوضح النتائج الناتجة عن اختبارات متنوعة، مع تسليط الضوء على الاتجاهات البيانية الهامة والتحليلات الإحصائية التي تدعم الفرضيات. غالبًا ما يتم توضيح النتائج من خلال الجداول والأشكال، التي توفر تمثيلات بصرية للبيانات، مما يسهل تفسير النتائج.

بالإضافة إلى ذلك، قد يتضمن القسم مقارنات بين المجموعات التجريبية، مما يظهر فعالية العلاجات أو التدخلات المختلفة. عادة ما يتم الإبلاغ عن الأهمية الإحصائية، غالبًا باستخدام قيم p أو فترات الثقة، للتحقق من النتائج. بشكل عام، يخدم هذا القسم لتأكيد أهداف البحث ويضع الأساس للنقاشات والاستنتاجات اللاحقة المستمدة من النتائج.

نقاش

في هذه الدراسة، بحث المؤلفون في أداء التقاط اليود لـ 1816 مادة من الهياكل العضوية المعدنية (MOF) تحت ظروف الهواء الرطب، مع التأكيد على أهمية إزالة اليود بكفاءة في إدارة النفايات النووية وتخفيف الحوادث. جمعت الدراسة بين الفحص الحسابي عالي الإنتاجية وتقنيات التعلم الآلي لاستكشاف العلاقة بين الخصائص الهيكلية (مثل قطر المسام المحدود، المساحة السطحية، والكثافة) وخصائص امتصاص اليود. أظهرت النتائج الرئيسية أن الامتصاص الأمثل لليود يحدث مع أقطار مسام محدودة تتراوح بين 4 Å و7.8 Å، وأن وجود هياكل حلقية سداسية وذرات النيتروجين في إطار MOF يعزز بشكل كبير من التقاط اليود. كما حددت الدراسة معامل هنري وحرارة الامتصاص كعوامل كيميائية حاسمة تؤثر على أداء امتصاص اليود.

تم استخدام خوارزميات التعلم الآلي، وخاصة الغابة العشوائية وCatBoost، للتنبؤ بقدرات امتصاص اليود، حيث حقق الأخير دقة تنبؤ عالية (R² = 0.941) عند دمج الوصف الهيكلي والجزيئي والكيميائي. تم استخدام طريقة SHAP لترتيب أهمية الميزات، مما كشف أن قطر المسام المحدود كان الوصف الهيكلي الأكثر تأثيرًا. بالإضافة إلى ذلك، تم تقديم بصمات جزيئية لتوفير رؤى تفصيلية حول الميزات الهيكلية التي تؤثر على امتصاص اليود، حيث أظهرت بصمة MACCS أداءً تنبؤيًا متفوقًا. بشكل عام، تؤسس هذه العمل إطارًا شاملاً لتصميم مواد MOF عالية الأداء، مما يعزز فهم العوامل التي تحكم امتصاص اليود في البيئات الرطبة.

Journal: npj Computational Materials, Volume: 11, Issue: 1
DOI: https://doi.org/10.1038/s41524-025-01617-2
Publication Date: 2025-05-02
Author(s): Haoyi Tan et al.
Primary Topic: Metal-Organic Frameworks: Synthesis and Applications

Methods

In this study, Grand Canonical Monte Carlo (GCMC) simulations were conducted using RASPA software to analyze the adsorption behavior of iodine (I₂) in metal-organic frameworks (MOFs) under conditions relevant to the nuclear industry, specifically at a temperature of 423 K and a pressure of 1 bar. The simulation environment mimicked high humidity conditions encountered during the post-treatment phase of spent nuclear fuel, consisting of a mixed gas system with 300 ppm I₂, 68.5% N₂, 18.4% O₂, and 12.2% H₂O, achieving 100% relative humidity. The MOFs were treated as fixed rigid structures with periodic boundary conditions, and supercells were utilized to ensure system dimensions exceeded twice the cutoff distance of 12 Å.

The selectivity of I₂ during adsorption was calculated using the formula:

\[
\text{selectivity}_{I_2} = \frac{X_{I_2}/Y_{I_2}}{X_{\text{others}}/Y_{\text{others}}}
\]

where \(X_{I_2}\) and \(Y_{I_2}\) represent the uptake amounts and gas phase concentration of I₂, while \(X_{\text{others}}\) and \(Y_{\text{others}}\) correspond to those of other gas components (N₂, O₂, and H₂O). Each GCMC simulation included an equilibration phase of 50,000 cycles followed by a production phase of 50,000 cycles, involving various molecular movements. Iodine was modeled as a spherical entity with van der Waals parameters derived from its viscosity, while water was represented using the TIP3P model. N₂ and O₂ were modeled as three-site representations, and the Universal Force Field (UFF) was employed for the Lennard-Jones parameters of the MOF structures. Interatomic interactions were described using both Lennard-Jones and electrostatic potential energy functions.

Results

The “Results” section of the research paper presents the key findings derived from the conducted experiments or analyses. It details the outcomes of various tests, highlighting significant data trends and statistical analyses that support the hypotheses. The results are often illustrated through tables and figures, which provide visual representations of the data, making it easier to interpret the findings.

Additionally, the section may include comparisons between experimental groups, showcasing the effectiveness of different treatments or interventions. Statistical significance is typically reported, often using p-values or confidence intervals, to validate the results. Overall, this section serves to substantiate the research objectives and lays the groundwork for the subsequent discussion and conclusions drawn from the findings.

Discussion

In this study, the authors investigated the iodine capture performance of 1816 metal-organic framework (MOF) materials under humid air conditions, emphasizing the significance of efficient iodine removal in nuclear waste management and accident mitigation. The research combined high-throughput computational screening and machine learning techniques to explore the relationship between structural characteristics (such as pore limiting diameter, surface area, and density) and iodine adsorption properties. Key findings indicated that optimal iodine adsorption occurs with pore limiting diameters between 4 Å and 7.8 Å, and that the presence of six-membered ring structures and nitrogen atoms in the MOF framework significantly enhances iodine capture. The study also identified the Henry’s coefficient and heat of adsorption as critical chemical factors influencing iodine adsorption performance.

Machine learning algorithms, specifically Random Forest and CatBoost, were employed to predict iodine adsorption capabilities, with the latter achieving a high prediction accuracy (R² = 0.941) when incorporating structural, molecular, and chemical descriptors. The SHAP method was utilized to rank feature importance, revealing that the pore limiting diameter was the most influential structural descriptor. Additionally, molecular fingerprints were introduced to provide detailed insights into the structural features affecting iodine adsorption, with the MACCS fingerprint demonstrating superior predictive performance. Overall, this work establishes a comprehensive framework for the targeted design of high-performance MOF materials, enhancing the understanding of factors governing iodine adsorption in humid environments.